MeteoinfoLab脚本示例：中国地形数据提取及3维地形图

OomaulenOo

MeteoInfo 发表于 2020-5-29 21:24
ax.text()

好的 非常感谢

油油切克闹

学习了，赞！

HANSEN

请问一下 直接用NCL读取您这个地形的nc文件 画地形 为什么报错说没有经纬度坐标信息呢？报错如下
Variable: height
Type: short
Total Size: 4830422 bytes
            2415211 values
Number of Dimensions: 2
Dimensions and sizes:   [lat | 1201] x [lon | 2011]
Coordinates:
            lat: [15..55]
            lon: [70..137]
Number Of Attributes: 3
  long_name :   Longitude
  units :       degrees_east
  name :        z
(0)     check_for_y_lat_coord: Warning: Data either does not contain
(0)     a valid latitude coordinate array or doesn't contain one at all.
(0)     A valid latitude coordinate array should have a 'units'
(0)     attribute equal to one of the following values:
(0)         'degrees_north' 'degrees-north' 'degree_north' 'degrees north' 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north'
(0)     check_for_lon_coord: Warning: Data either does not contain
(0)     a valid longitude coordinate array or doesn't contain one at all.
(0)     A valid longitude coordinate array should have a 'units'
(0)     attribute equal to one of the following values:
(0)         'degrees_east' 'degrees-east' 'degree_east' 'degrees east' 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'

kasal2010

wow 太给力了 学习了！

Tiangy

老师您好，如何将地形设置为灰色呢，这个切面图的地形颜色，谢谢

yanmuyan

很给力！谢谢

凉白开

谢谢分享

小浪

这个精细程度，不像一度的分辨率的啊

Finnick

MeteoInfo 发表于 2020-4-30 10:44
MeteoInfoLab软件Figure区域有专门的Rotate工具按钮，选中后可以手动旋转三维图形。或者用三维图形坐标系 ...

王老师 我用一楼一模一样的脚本跑出来的盆地图横纵坐标标注是有问题的，想问一下王老师出现这个问题是什么原因啊？
#read data
fn = r'D:\test/dixing-china.nc'
f = addfile(fn)
elev = f['z']['26:34','102:108']/1000.0
lchina = shaperead('D:/Temp/map/bou2_4p.shp')

#Plot
figure(figsize=[1400,850],newfig=False)
ax = axes3d(tickfontsize=20,position=[0.1, 0.13, 0.75, 0.85],bbox=False)
levs = arange(0,6.1,0.1)
cols = makecolors(len(levs) + 1, cmap='temp_diff_18lev',reverse=False)
ls = ax.plot_surface(elev, levs, colors=cols, edge=False)
colorbar(ls,ticks=arange(0,6.1,1),orientation='horizontal',label='Altitude(km)',labelloc='right',fontsize=24,shrink=0.8,aspect=40,yshift=35)
ylabel('Latitude (degrees)',fontsize=20)
xlabel('Longitude (degrees)',fontsize=20)
zlabel('Altitude (km)',fontsize=20)
zlim(0, 6)
xlim(102, 108)
ylim(26,34)
picpath = 'E:/topographicmap.png'
plt.savefig(picpath,dpi=1200)
plt.clf()

MeteoInfo

Finnick 发表于 2022-2-13 16:55
王老师 我用一楼一模一样的脚本跑出来的盆地图横纵坐标标注是有问题的，想问一下王老师出现这个问题是什 ...

你用的什么版本的MeteoInfo ？dpi调低些试试。