2
\ begingroup美元

我想reproject格丽来自NOAA预测数据到其他数据使用Python和pygrib项目包(和其他很多)。适当的为GFS格丽项目字符串是什么?

投影应该来自于一种文件格式文件,我以为GFS使用墨卡托,但pygribgrib.projparams ()返回{' a ': 6371229, ' b ': 6371229,“项目”:“共青团”}。“共青团”不是一个合法的项目CRS和其他圆柱投影(中央圆柱板Carree)没有或b参数。假设a和b是定义地球是一个球体半径为6371229米。

grib.projparams ()与一个HRRR格丽的回报{' a ': 6371229, ' b ': 6371229,“项目”:“lcc”,“lon_0”: 262.5,“lat_0”: 38.5,“lat_1”: 38.5,“lat_2”: 38.5}

这是完全合法的项目和reprojects数据准确。

\ endgroup美元
0

    1回答1

    3
    \ begingroup美元

    我拟设为共青团意味着可能等距圆柱投影,因为在老技术模块,“共青团”该投影的缩写吗

    但现实是,由于GFS全球谱模式与经度网格,它可能不需要CRS。也就是说,它是以一个球体,球面坐标CRS项目(经度、纬度)到一个形状(如锥)。

    \ endgroup美元
    5
    • \ begingroup美元 @BaroclinicCplusplus GFS光谱模型吗? \ endgroup美元
      - - - - - -gansub
      2020年6月25日,在1:26
    • 1
      \ begingroup美元 @gansub是的,但光谱模型描述了模型是解决,不一定CRS的模型。虽然更容易在全球纬度/经度网格,有一些区域光谱模型,如RSM或ALADIN使用傅里叶和/或勒让德变换。后者的模型,例如,使用一个墨卡托或兰伯特正形投影(umr-cnrm.fr gmapdoc / IMG / pdf / ykarpbasics46t1r1.pdf) \ endgroup美元
      - - - - - -BarocliniCplusplus
      2020年6月25日17:57
    • \ begingroup美元 @BaroncliniCplussplus所以我的项目应该是字符串“+项目= eqc + R = 6371229”eqc在哪里等距圆柱(板Carree)球的半径和R(给定的一个b底层的主要和次要的轴椭圆,相等)? \ endgroup美元
      - - - - - -兰斯Berc
      2020年6月26日,在18:37
    • \ begingroup美元 @LanceBerc如果这是你想要的,当然。但投影是没有实际意义的选择。CRS试图项目上一场球,就好像它是一个几何对象,说一个圆柱体。地图投影的事项如果模型的投影。但自从GFS没有投影,然后没有投影问题。这完全是你的选择。GFS的所有这一切,你的选择在地图上只会改变地图的样子,没有数据如何在地图上看起来。 \ endgroup美元
      - - - - - -BarocliniCplusplus
      2020年6月26日,在上19:15
    • 1
      \ begingroup美元 @BarocliniCplusplus我明白你的意思。但自从我使用pyproj.Transformer.from_crs()将经度/纬度x / y,它希望CRS来源;板Carree似乎是一个不错的选择。谢谢你! \ endgroup美元
      - - - - - -兰斯Berc
      2020年6月26日21:29

    你的答案

    通过点击“发布你的答案”,你同意我们服务条款,隐私政策饼干的政策

    不是你要找的答案?浏览其他问题标记问你自己的问题