关于WRF中lambert投影...

ShuyingChenUCAS

楼主首先非常惭愧，一个根正苗红的GIS生，上手wrf几个月了吧，逛了逛家园，发现自己对投影、对坐标的理解远不及众大气科学的各位...决心起码要搞清自己现在在用的lambert投影，看了一些资料后还是有些恍惚，不甘心，自己来开一个帖子希望能和大家讨论一下~


先说ArcGIS，是我们GISer的看家本领，每个学GIS的人对这个软件应该都很熟悉，但不幸的是这里面有关投影这部分东西对我而言一直是灰色地带，学艺不精...
说回WRF，最近几个月一直在做WRF的模拟，刚上手就一直想着赶紧跑起来，现在有些成果了，开始要分析，就有些傻眼。我的模拟区在中亚，看wrf手册推荐使用lambert投影，所以我用了... domain我是用domain_wizard画的，最后生成namelist.wps对lamber投影的描述大概是下面这几行代码：
   map_proj =  'lambert',     投影名称，兰勃特投影，正轴等角割圆锥（第一、二标准纬线相等：正轴等角切圆锥）
   ref_lat   = 42,    中心点纬度
   ref_lon   = 80,   中心点经度
   truelat1  = 42,   第一标准纬线
   truelat2  = 42,   第二标准纬线
   stand_lon = 80,  中央经线
   ref_x = 50.0,   中心点的格网x坐标
   ref_y = 25.0,   中心点的格网y坐标

那么，
1.我想把这个投影在ArcGIS中定义出来应该怎么做呢？这个投影的datum是WGS-84地理坐标系吗？
2.lambert投影坐标的wrfout文件（加.nc以后）不能直接在ArcGIS中打开，lat-long坐标、mercator坐标都可以，有办法解决吗？
3.我想把上面domian_wizard画出来的格网（domain）在ArcGIS里用fishnet同样画出来的话，fishnet的参数应该怎么设置就能和wrf中的格网完全一致？（finsnet在画格网的时候可以选择定义起始点或者区域范围，现在这个投影下我只知道一个中心点的经纬度坐标、格子长宽（m）、上下左右的格子数，既无法准确推出起始点的经纬度值也无法准确推出区域范围的经纬度值...）
另外，
4.ncl有没有可以对地理坐标投影的函数？我的wrf输出是lambert，当我想把WGS-84的shp文件叠到wrfout上时，空间位置没有对应起来。
5.在验证结果是要用到一些比如说TRMM降水产品是WGS-84坐标系，我的WRFout(降水结果)又是lambert，又该怎么对应起来？


一堆问题，无解，欢迎吐槽、闲聊~求大神

ShuyingChenUCAS

Projection: Lambert_Conformal_Conic
False_Easting: 0.0
False_Northing: 0.0
Central_Meridian: 80.0
Standard_Parallel_1: 42.0
Standard_Parallel_2: 42.0
Scale_Factor: 1.0
Latitude_Of_Origin: 42.0
Linear Unit: Meter (1.0)

Geographic Coordinate System: GCS_WGS_1984
Angular Unit: Degree (0.0174532925199433)
Prime Meridian: Greenwich (0.0)
Datum: D_WGS_1984
  Spheroid: WGS_1984
    Semimajor Axis: 6378137.0
    Semiminor Axis: 6356752.314245179
    Inverse Flattening: 298.257223563
是不是应该在ArcGIS里对这个投影进行如上定义？

packard

不懂ArcGIS，不懂WGS84，只能说说你的第四个问题：
可以看看 wrf_ll_to_ij 和 wrf_map_resources 这两个ncl 函数，以及ncl 官网里面名字类似的函数。
别的我都不懂。

WRF的相关代码在 share/module_llxy.F
我扫了一眼，如果没有特别指明是PS_WGS84的话，只是一般的 Lambert Conformal 的话，那么地球长短轴都是6370000m。不知道我理解那个code对不对，你也可以看看，我的WRF版本跟你可能也有差异。

   REAL, PARAMETER :: EARTH_RADIUS_M = 6370000.
      proj%re_m     = EARTH_RADIUS_M
            proj%rebydx = proj%re_m / dx

   SUBROUTINE ijll_lc( i, j, proj, lat, lon)  里面没有WGS84字眼，函数看起来简单一些
      r = SQRT(r2)/proj%rebydx

   SUBROUTINE llij_ps_wgs84(lat,lon,proj,i,j)  这个函数看起来复杂一些，也有WGS84的字眼。

ShuyingChenUCAS

packard 发表于 2018-1-13 08:01
不懂ArcGIS，不懂WGS84，只能说说你的第四个问题：
可以看看 wrf_ll_to_ij 和 wrf_map_resources 这两个nc ...

多谢回复，还要好好修炼，是应该多看看wrf里面的源码...

一介白衣

我想请问一下 第一、二标准纬线和中心经线是怎么确定下来的？

香草拿铁

一介白衣 发表于 2019-9-23 15:38
我想请问一下 第一、二标准纬线和中心经线是怎么确定下来的？

可以以你的研究中心作为中心，设置中心经线，标准纬线

Gumho

香草拿铁 发表于 2019-9-25 14:48
可以以你的研究中心作为中心，设置中心经线，标准纬线

谢谢您，懂了。{:loveliness:}{:loveliness:}

xyy123

请问你怎么解决的啊

Rainzmh

Gumho 发表于 2019-9-25 16:31
谢谢您，懂了。

什么意思，我是小白，可以请教一下具体是什么意思吗