爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4424|回复: 3

投影方式转换

[复制链接]

新浪微博达人勋

发表于 2021-8-18 14:35:36 | 显示全部楼层 |阅读模式

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

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

x
请问实现Lambert投影转换成等经纬度投影的函数是什么?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2021-8-19 09:31:42 | 显示全部楼层

问题解决啦!现将示例脚本附上,供有需要的人参考呀~

本帖最后由 ymy_00000 于 2021-8-19 09:33 编辑


fn = 'D:/Temp/grib/Z_NAFP_C_BCSH_20190301000000_P_surface-warms-f00.BIN'
f = addfile(fn)
data = f['Temperature_surface'][:]
yn, xn = data.shape

f.proj = projinfo(proj='lcc', lat_0=35, lon_0=105, lat_1=30, lat_2=60)
lon0 = 78.30664
lat0 = 7.29969
x0, y0 = geolib.project(lon0, lat0, toproj=f.proj)
xx = arange1(x0, xn, 9000)
yy = arange1(y0, yn, 9000)

ds_,lon,lat = geolib.reproject(data, xx, yy, fromproj=f.proj,
    toproj=projinfo())

axesm()
geoshow('cn_province', edgecolor='gray')
geoshow('country')
imshow(lon, lat, ds_, 28)
colorbar()
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2021-8-23 21:32:49 | 显示全部楼层
{:5_213:}{:5_213:}{:5_213:}{:5_213:}
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-12-9 10:55:52 | 显示全部楼层
ymy_00000 发表于 2021-8-19 09:31
fn = 'D:/Temp/grib/Z_NAFP_C_BCSH_20190301000000_P_surface-warms-f00.BIN'
f = addfile(fn)
data =  ...

您好,这段代码有java的实现方式吗
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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