MeteoInfoLab脚本示例：CMORPH降水数据

MeteoInfo

chongzika 发表于 2017-9-6 14:19
从ucar下载的数据，原来的数据为502.0，现在定义为502.1，更新了算法This dataset contains version 0.x  ...

这个数据是little_endian，参考下面的脚本：

1. fn = 'D:/Temp/binary/CMORPH_V0.x_RAW_0.25deg-3HLY_20170610'
   
2. ntime = 8
   
3. nlat = 480
   
4. nlon = 1440
   
5. data = binread(fn, [ntime,nlat,nlon], 'float', byteorder='little_endian')
   
6. rain = data[2,:,:]
   
7. lon = arange1(0.125, nlon, 0.25)
   
8. lat = arange1(-59.875, nlat, 0.25)
   
9. #Plot
   
10. axesm()
    
11. geoshow('country', edgecolor='k')
    
12. geoshow('cn_province', edgecolor='k')
    
13. levs = [0.1,0.5,1,2,3,4,5,6,8,10,20,40]
    
14. cols = makecolors(len(levs) + 1, cmap='BlAqGrYeOrReVi200')
    
15. cols[0] = 'w'
    
16. layer = imshowm(lon, lat, rain, levs, colors=cols)
    
17. colorbar(layer, label='mm/3hr')
    
18. title('Precipitation - CMORPH (2017-06-10 06:00)')
    
19. #xlim(70, 140)
    
20. #ylim(15, 55)

fn = 'D:/Temp/binary/CMORPH_V0.x_RAW_0.25deg-3HLY_20170610'<br /> ntime = 8<br /> nlat = 480<br /> nlon = 1440<br /> data = binread(fn, [ntime,nlat,nlon], 'float', byteorder='little_endian')<br /> rain = data[2,:,:]<br /> lon = arange1(0.125, nlon, 0.25)<br /> lat = arange1(-59.875, nlat, 0.25)<br /> #Plot<br /> axesm()<br /> geoshow('country', edgecolor='k')<br /> geoshow('cn_province', edgecolor='k')<br /> levs = [0.1,0.5,1,2,3,4,5,6,8,10,20,40]<br /> cols = makecolors(len(levs) + 1, cmap='BlAqGrYeOrReVi200')<br /> cols[0] = 'w'<br /> layer = imshowm(lon, lat, rain, levs, colors=cols)<br /> colorbar(layer, label='mm/3hr')<br /> title('Precipitation - CMORPH (2017-06-10 06:00)')<br /> #xlim(70, 140)<br /> #ylim(15, 55)

chongzika

MeteoInfo 发表于 2017-9-6 15:02
这个数据是little_endian，参考下面的脚本：

可以了，多谢王老师指导。数据格式不同了么？两者什么差别呢？另外发现，win10下面的Meteoinfo不能正确识别默认shp文件目录File not exists: C:/Users/lch/software/MeteoInfo/map/contoury.shp

MeteoInfo

chongzika 发表于 2017-9-6 17:01
可以了，多谢王老师指导。数据格式不同了么？两者什么差别呢？另外发现，win10下面的Meteoinfo不能正确识 ...

你拼写错误吧

jerrychu

王老师，问一个弱弱的问题，这个软件可以显示，怎么把数据导出来？

可以再写一个脚本文件是循环读取导出数据的，供我们参考吗？谢谢

MeteoInfo

jerrychu 发表于 2017-11-21 12:24
王老师，问一个弱弱的问题，这个软件可以显示，怎么把数据导出来？

可以再写一个脚本文件是循环读取导出 ...

导出数据可以用asciiwrite, binwrite或者ncwrite函数，参考这里：http://www.meteothinker.com/docs/meteoinfolab/dataset/midata.html

fresh_fish

谢谢分享！！

不说话就很美

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

e327741

liujie2462 发表于 2017-6-15 16:15
王老师，您好。请问该数据的第一个点，经度0.125，纬度-59.875，对应的是matlab读取后的第一个数据吗？期待 ...

请问，你是怎么用matlab打开的这个数据？这个数据到底是什么格式的文件？

有盐呀喜欢安静

老师您好，冒昧打扰了，看到您2016年12月31日发布的关于《中国自动站与CMORPH降水产品融合的逐时降水量网格数据集（1.0版）》的帖子，所以想来问问您中国自动气象站分布图可以在哪里获取？就是《中国自动站与CMORPH降水产品融合的逐时降水量网格数据集（1.0版）》里面用到的全国自动气象站点数据。打扰了，祝好

大雪纷飞

王老师，请问pscf分析时如何使用脚本将输出文件转化为.dat格式呢？