我从http://ngdc.noaa.gov/mgg/global/global.html下载了一个文件,在那里我提取了一个自定义网格来解决我的海洋学问题。但是,此网格具有更高的分辨率(1分钟或1/60)。比我需要运行我的实验。因此,我想知道如果我想有一个10分钟的网格分辨率,我如何“调整”它使用for循环?我知道如何在MATLAB中实现它(代码粘贴在下面),但由于我想成为一个python用户,我不知道如何做到这一点。你知道吗
MATLAB代码:
ii=0;
for i=1:10:1681;
ii=ii+1;
jj=0;
for j=1:10:4561;
jj=jj+1;
n=j+((i-1)*4561);
long2(ii,jj)=lon(n);
latie2(ii,jj)=lat(n);
c2(ii,jj)=depth(n);
end
end
深度文件具有以下结构,其数组大小为1621x4561:
array([[ 200, 191, 192, ..., 1344, 1345, 1343],
[ 191, 178, 192, ..., 1330, 1332, 1337],
[ 178, 176, 177, ..., 1297, 1303, 1314],
...,
[-1806, -1853, -1897, ..., 295, 296, 295],
[-1803, -1845, -1887, ..., 294, 295, 295],
[-1806, -1844, -1881, ..., 294, 294, 295]], dtype=int16)
提前谢谢
你可以尝试用ncks每十个lon和lat超标签一次
相关问题 更多 >
编程相关推荐