MeteoInfoLab脚本示例：CMORPH降水数据

chongzika · 发表于 2017-7-28 15:35:32

MeteoInfo 发表于 2017-4-17 23:42
1楼不是有脚本程序吗

王老师我想问下linux 下meteinfo怎么运行脚本，./milab.sh CMORPH.py 这样样子就可以了么？

MeteoInfo · 发表于 2017-7-28 16:29:51

chongzika 发表于 2017-7-28 15:35
王老师我想问下linux 下meteinfo怎么运行脚本，./milab.sh CMORPH.py 这样样子就可以了么？

是的

chongzika · 发表于 2017-9-5 17:56:31

王老师好，我想问一下我在服务器上运行，./milab.sh CMORPH.py后怎么自动设置保存输出图片呢？
每次运行完$ ./milab.sh CMORPH.py
1
Running Jython script...
mipylib is loaded...
$
什么都没看到？就结束了，不知道怎么保存图片

MeteoInfo · 发表于 2017-9-5 18:05:23

chongzika 发表于 2017-9-5 17:56
王老师好，我想问一下我在服务器上运行，./milab.sh CMORPH.py后怎么自动设置保存输出图片呢？
每次运行完 ...

savefig函数

chongzika · 发表于 2017-9-5 19:02:15

MeteoInfo 发表于 2017-9-5 18:05
savefig函数

直接使用savefig('cmo.png')，保存出来的图片分辨率特别低，只有200。设置savefig('cmo.png',dpi=600)保存的图片是空的？王老师这个像素要如何设置？

MeteoInfo · 发表于 2017-9-5 21:23:10

chongzika 发表于 2017-9-5 19:02
直接使用savefig('cmo.png')，保存出来的图片分辨率特别低，只有200。设置savefig('cmo.png',dpi=600)保 ...

可以存为eps或者pdf矢量格式图形

chongzika · 发表于 2017-9-5 22:10:39

MeteoInfo 发表于 2017-9-5 21:23
可以存为eps或者pdf矢量格式图形

试了下，linux下面运行有如下报错，window下面口令正常。
savefig('test.pdf')
  File "/data/c03n02/cliu/software/meteinfo/MeteoInfo/pylib/mipylib/plotlib/miplot.py", line 2292, in savefig
chartpanel.saveImage(fname, sleep)
      at org.freehep.graphicsio.pdf.ImageBytes.toZLIB(ImageBytes.java:80)
      at org.freehep.graphicsio.pdf.ImageBytes.<init>(ImageBytes.java:46)
      at org.freehep.graphicsio.pdf.PDFStream.image(PDFStream.java:627)
      at org.freehep.graphicsio.pdf.PDFImageDelayQueue.processAll(PDFImageDelayQueue.java:85)
      at org.freehep.graphicsio.pdf.PDFGraphics2D.processDelayed(PDFGraphics2D.java:470)
      at org.freehep.graphicsio.pdf.PDFGraphics2D.closePage(PDFGraphics2D.java:585)
      at org.freehep.graphicsio.pdf.PDFGraphics2D.writeTrailer(PDFGraphics2D.java:389)
      at org.freehep.graphicsio.AbstractVectorGraphicsIO.endExport(AbstractVectorGraphicsIO.java:260)
      at org.meteoinfo.chart.ChartPanel.saveImage(ChartPanel.java:1002)
      at org.meteoinfo.chart.ChartPanel.saveImage(ChartPanel.java:932)
      at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)
      at sun.reflect.NativeMethodAccessorImpl.invoke(NativeMethodAccessorImpl.java:57)
      at sun.reflect.DelegatingMethodAccessorImpl.invoke(DelegatingMethodAccessorImpl.java:43)
      at java.lang.reflect.Method.invoke(Method.java:606)
java.lang.ExceptionInInitializerError: java.lang.ExceptionInInitializerError

chongzika · 发表于 2017-9-6 02:41:24

本帖最后由 chongzika 于 2017-9-6 03:07 编辑

发现CMORPH 的数据已经做了更新调整，数据里只有一个变量CMORPH precipitation estimate

fn='C:/data/CMORPH_ds502_1/CMORPH_V0.x_RAW_0.25deg-3HLY_20170610'
nvar = 1
ntime = 8
nlat = 480
nlon = 1440
data = binread(fn, [nvar*ntime,nlat,nlon], 'float', byteorder='big_endian')
#comb = data[::2,::-1,:]
cpc = data[:,::-1,:]
#cpc[cpc==999.0] = nan
rain = cpc[1,:,:]
lon = arange1(0.125, nlon, 0.25)
lat = arange1(-59.875, nlat, 0.25)
#Plot
axesm()
#lworld = shaperead('C:/data/map/CHN_adm/CHN_adm2.shp')
lchina = shaperead('C:/data/map/CHN_adm/CHN_adm1.shp')
geoshow(lchina, edgecolor='k')
#geoshow(lworld, edgecolor='k')
levs = [1,2,4,8,10,20,60,80,120,160]
cols = makecolors(len(levs) + 1, cmap='BlAqGrYeOrReVi200')
cols[0] = 'w'
layer = imshowm(lon, lat, rain, levs, colors=cols)
colorbar(layer, label='mm')
title('Precipitation - CMORPH (2017-06-10 03:00 UTC)')
xlim(70, 140)
ylim(15, 60)

复制代码

把脚本稍做了修改，发现出来的图有些怪异，

又如此图，图中的点只有1.3对应的色标特别大

MeteoInfo · 发表于 2017-9-6 08:28:36

chongzika 发表于 2017-9-6 02:41
发现CMORPH 的数据已经做了更新调整，数据里只有一个变量CMORPH precipitation estimate

把 ...

你的数据从哪里下载的？数据说明是什么？

chongzika · 发表于 2017-9-6 14:19:31

本帖最后由 chongzika 于 2017-9-6 14:53 编辑

MeteoInfo 发表于 2017-9-6 08:28
你的数据从哪里下载的？数据说明是什么？

从ucar下载的数据，原来的数据为502.0，现在定义为502.1，更新了算法This dataset contains version 0.x of the NOAA CPC MORPHing technique (CMORPH) global precipitation analyses covering the period June 2014-present. Version 0.x comprises the original CMORPH precipitation analyses that CPC has been generating since CMORPH became operational in December 2002, and is generated using an improving algorithm with inputs of evolving versions. Older CMORPH version 0.x data (prior to June 2014) may be accessed from RDA dataset ds502.0.
网页为数据介绍及下载

CMORPH_V0.x_RAW_0.25deg-3HLY_20170613 (21.09 MB, 下载次数: 3)

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