我怎么纬度/经度坐标转换为NOAA GFS网格单元的位置吗?

Question

看这个GFShttp://nomads.ncep.noaa.gov: 9090 /国防部gfs_0p50 / gfs20170707 / gfs_0p50_12z.ascii ? tmax2m [0:2] [124:128] [279:283]

范围定义的第二和第三组括号指定网格坐标,不过我不确定如何将NOAA之间使用和常见的地理坐标系统。

NOAA当量是一对坐标:(36.25,-78.75)?

Answer 1

* *所有链接的时间需要调整到当前日期如果阅读这个在以后的时代里,否则会死* *

从前我办理这些文件,但不是在好几年……和谷歌不是最好的帮助首先找到的文档.ascii正是……但由于凯文Goebbert的代码在这里,我能够发现问题的描述页面的模型位于.ascii和输入(通过删除,因此在这里):

所以现在列表:

经度:0.00000000000°E 359.50000000000°E(720分,平均0.5°)》纬度:-90.00000000000°N 90.00000000000°N(361分,平均0.5°)》高度:1000.00000000000 - 1.00000000000(47分,平均> 21.717)时间:12 z07jul2017 12 z17jul2017(81分,平均0.125天》)

只有一半的p50表明网格空间学位该文件,所以你不能在这个moel实际使用你的观点。

然而,使用0.25解决而不是(描述):

首先是时间第一是时间,第二是经度点开始在本初子午线和东,第三个是纬度开始在南极和北。(GFS时间步骤总是一个泥潭……但是基于这完整列表数据文件和信息页面之前,它确实出现这里提供的数据是每3小时前10天。)

选择你的第二个输入值,这将是类似

ifpositive:(经度)* 4 - 1 ifnegative:经度(180 +)* 4 + 719

和纬度:

360 +(纬度)* 4

所以如果我没有犯错误,会这一点(现纬度现505年,经度1124)。

Answer 2

中定义的坐标全球预测系统(GFS -0.5º学位)数据定义有些不同寻常。GFS,经度和纬度包含值分别为720人和361人。因此,将它们转换为一个共同的地理坐标系统遵循以下程序:

假设你正在寻找协调双指数(36.5,-78.5)GFS数据:

your_lat_start = x your_lat_end = x your_lon_start = x your_lon_end = x your_dataset = half_degree_gfs_dataset lat_index_start = np.where (np.array (your_dataset[“纬度”][:]. data) = = your_lat) [0] [0] lat_index_end = np.where (np.array (your_dataset[“纬度”][:]. data) = = your_lat_end) [0] [0] lon_index_start = np.where (np.array (your_dataset[“朗”][:]. data) = = your_lon_start) [0] lon_index_end = np.where (np.array (your_dataset[“朗”][:]. data) = = your_lon_end)[0][0]如果lat_index_start > lat_index_end: lat_index_start, lat_index_end = lat_index_end lat_index_start + 1 lon_index_start, lon_index_end = lon_index_start lon_index_end + 1数据集[var_name] [: lat_index_start: lat_index_end, lon_index_start: lon_index_end] . data