编辑-地球科学江南体育网页版堆栈交换江南电子竞技平台

你不是登录。你的编辑将被放置在一个队列,直到同行评议。

我们欢迎编辑,使读者更容易理解和更有价值。因为社区成员审查编辑,请尽量使文章更好比你怎么发现,例如,通过修正语法或添加额外的资源和超链接。

牧师

标题

身体

我想选择一个地区的地中海ETOPO1数据。#深度测量法从netCDF4导入数据集数据=数据集(“ETOPO1_Bed_g_gdal.grd”、“r”)打印(data.variables.keys ()) # lon_range =数据获取数据。变量[' x_range '] [:] lat_range =数据。变量[' y_range '] [:] topo_range =数据。变量[' z_range '][:]间距=数据。变量=[“间距”][:]维度数据。变量(“维度”][:]z =数据。变量[' z '] [:] lon_num维度[0]= lat_num =维度[1]#准备数组lon_input = np.zeros (lon_num);lat_input = np.zeros (lat_num)我的范围(lon_num): lon_input[我]= lon_range [0] + i *间距[0]我的范围(lat_num): lat_input[我]= lat_range [0] + i *间距[1]#创建二维数组经度,纬度= np。meshgrid (lon_input lat_input) # 2 d数组从一维数组转换为z值威尼斯平底渔船= np。重塑(z (lat_num lon_num) #选择地区lon_idx_min = np.abs (lon_input - 12.6) .argmin () lon_idx_max = np.abs (lon_input - 13.4) .argmin () lat_idx_min = np.abs (lat_input - 43.8) .argmin () lat_idx_max = np.abs (lat_input - 44.2) .argmin() #威尼斯平底渔船获取数据(lon_idx_min: lon_idx_max lat_idx_min: lat_idx_max]。数据返回一个空数组:数组> >([],形状= (0,24),dtype = float64)的一个主要问题:* * 1 -我做错了什么?* *当我检查经度纬度值我切片(lat_input [lat_idx_max]等)他们是有意义的。然而,lon_idx_max lon_idx_min高于len(威尼斯平底渔船)(11556 vs 10801和11604年)。我想这是问题,但我不能找到它的原因。一个额外的问题:2 -为什么如果我翻转经度和纬度给我一些数据吗? Note that the data do not make sense (-4000ish meters in the Adriatic Sea...)

标签

编辑总结

取消

\ begingroup美元嗨,x_range是从哪里来的,因为当我试图重建你的例子python不知道它对我来说 \ endgroup美元
- - - - - -维斯
2022年6月22日,星期日
\ begingroup美元当我打印数控文件的变量,它是沿着y_range和z_range \ endgroup美元
- - - - - -shamalaia
2022年6月22日23点

添加一个评论 |

正确的小错误或错误
澄清的意思,而不更改它
添加相关资源或链接
总是尊重作者的意图
不使用编辑与作者回复