请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

新浪微博登陆

只需一步, 快速开始

QQ登录

只需一步,快速开始

搜索
查看: 1031|回复: 3

关于WRF中lambert投影...

[复制链接] |关注本帖

新浪微博达人勋

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
楼主首先非常惭愧,一个根正苗红的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,又该怎么对应起来?


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

已有2人关注本帖

问题少女香草拿铁
密码修改失败请联系qq:937062711

新浪微博达人勋

 楼主| 发表于 2018-1-12 19:13:23 | 显示全部楼层 |取消关注该作者的回复
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里对这个投影进行如上定义?
密码修改失败请联系qq:937062711

新浪微博达人勋

发表于 2018-1-13 08:01:41 | 显示全部楼层 |取消关注该作者的回复
不懂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的字眼。
密码修改失败请联系qq:937062711

新浪微博达人勋

 楼主| 发表于 2018-1-13 12:14:14 | 显示全部楼层 |取消关注该作者的回复
packard 发表于 2018-1-13 08:01
不懂ArcGIS,不懂WGS84,只能说说你的第四个问题:
可以看看 wrf_ll_to_ij 和 wrf_map_resources 这两个nc ...

多谢回复,还要好好修炼,是应该多看看wrf里面的源码...
密码修改失败请联系qq:937062711
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表