投影方式转换

ymy_00000 · 发表于 2021-8-18 14:35:36

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

您需要登录才可以下载或查看，没有帐号？立即注册

x

请问实现Lambert投影转换成等经纬度投影的函数是什么？

ymy_00000 · 发表于 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()

cc123 · 发表于 2021-8-23 21:32:49

{:5_213:}{:5_213:}{:5_213:}{:5_213:}

wsc123856 · 发表于 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的实现方式吗

		自动登录	找回密码
密码			立即注册

投影方式转换

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

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