MeteoInfoLab脚本示例：站点数据绘制等值线

gjj123

MeteoInfo 发表于 2016-9-18 10:44
主要是插值后格点数据坐标的问题。脚本中
layer = contourfm(prg, levs, colors=cols)
改为

好的，确实是之前没有指定经纬度坐标，谢谢老师

wuwei2163

王老师，您好！今天学习站点插值，以前基本上都是Fortran+grads作图，相对来说MeteoInfoLab是比较简便的，折腾了一个多小时，终于实现了站点数据插值，并且与其他数据填色叠加，这是王老师帮助的结果。

但是这其中有几个小问题，在help(griddata)以及教程的帮助下，图是画出来了，但看到您上面脚本中插值函数是这个：

另外又看到http://bbs.06climate.com/forum.php?mod=viewthread&tid=44833&extra=page%3D1贴上关于读取2015年1月的温度数据（0.1度为单位），插值为格点数据并绘制等值线上的插值：


而我是按照这个插值的，这两个有何不同吗？请王老师指导~

MeteoInfo

wuwei2163 发表于 2016-10-22 18:05
王老师，您好！今天学习站点插值，以前基本上都是Fortran+grads作图，相对来说MeteoInfoLab是比较简便的， ...

Python是面向对象的语言，脚本中pr是一个PyStationData对象，有griddata函数将站点插值为格点。minum模块中也有griddata函数（minum模块在软件启动时已经对其所有的函数全部加载，因此可以直接调用），这个适用性更广，建议用这个函数。

wuwei2163

MeteoInfo 发表于 2016-10-22 20:49
Python是面向对象的语言，脚本中pr是一个PyStationData对象，有griddata函数将站点插值为格点。minum模块 ...

王老师，有个小问题请教，设置等值线间隔和线条颜色如下：levs = [0,1 , 5, 10, 25, 50, 100]
layer= contourm(x, y, pr,levs,colors='black')，但 等值线太细了，在help(contourm)说明中也没找到如何在设置等值线粗细啊，比如加粗线条？标注是用clabel()函数，怎么在所有等值线上全部标注？clabel()函数似乎也没有说明

MeteoInfo

wuwei2163 发表于 2016-10-22 22:13
王老师，有个小问题请教，设置等值线间隔和线条颜色如下：levs = [0,1 , 5, 10, 25, 50, 100]
layer=  ...

clabel()函数里加avoidcoll=False参数来标注所有等值线。
等值线粗细等的设置我会尽快发布一个改进版本来解决。

wuwei2163

MeteoInfo 发表于 2016-10-23 11:10
clabel()函数里加avoidcoll=False参数来标注所有等值线。
等值线粗细等的设置我会尽快发布一个改进版本 ...

好的，谢谢王老师，否则等值线太细了，要是要叠加上色斑图，就不容易分辨了~

wuwei2163

MeteoInfo 发表于 2016-10-23 11:10
clabel()函数里加avoidcoll=False参数来标注所有等值线。
等值线粗细等的设置我会尽快发布一个改进版本 ...

王老师，您好！现在无论是发文章还是科研，对图表的要求越来越高，除了上次跟您提过关于设置等值线粗细的问题（坐等王老师新版本）外，还有关于图片保存提高图的分辨率的问题，看了帮助文档savefig()可以存成 'png', 'bmp', 'jpg'等格式，图片质量不是很高，有些美中不足。王老师您看有什么好的办法吗？

MeteoInfo

wuwei2163 发表于 2016-10-25 16:33
王老师，您好！现在无论是发文章还是科研，对图表的要求越来越高，除了上次跟您提过关于设置等值线粗细的 ...

已经上传了一个改进版本（1.3.6R4），给contourm函数添加size参数可以设置等值线粗细，比如：
layer = contourm(x, y, pr, levs, colors='k', size=2)

可以保存为eps或者pdf格式，这两种都是矢量格式，质量最高。

wuwei2163

MeteoInfo 发表于 2016-10-25 23:08
已经上传了一个改进版本（1.3.6R4），给contourm函数添加size参数可以设置等值线粗细，比如：
layer = c ...

好的，谢谢王老师~

wuwei2163

王老师，您好：
      之前一直通过站点数据作图并mask本省数据，没出现问题；今天上午想mask自定义地图，出现问题，先上图吧：

如图，我把代码中的23行换成24行，要么仅仅出现一个空白的自定义地图，要么出现下面的问题：
王老师，这是什么情况啊？自定义的地图shp也没问题啊？