MeteoInfoLab脚本示例：Sea ice netCDF数据

MeteoInfo

NOAA/NSIDC Climate Data Record of Passive Microwave Sea Ice Concentration数据是netCDF格式，极地投影，这里给出一个数据读取和绘图的示例脚本：

1. f = addfile('D:/Temp/nc/seaice_conc_daily_sh_f17_20080830_v02r00.nc')
   
2. data = f['seaice_conc_cdr'][0,:,:]
   
3. #Plot
   
4. axesm(projinfo=f.proj, gridline=True)
   
5. lworld = shaperead('D:/Temp/map/country1.shp')
   
6. geoshow(lworld)
   
7. layer = imshowm(data, 20, cmap='WhBlGrYeRe', proj=f.proj)
   
8. colorbar(layer)
   
9. t = f.gettime(0)
   
10. title('Sea ice concentration (' + t.strftime('%Y-%m-%d') + ')')

f = addfile('D:/Temp/nc/seaice_conc_daily_sh_f17_20080830_v02r00.nc')<br /> data = f['seaice_conc_cdr'][0,:,:]<br /> #Plot<br /> axesm(projinfo=f.proj, gridline=True)<br /> lworld = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(lworld)<br /> layer = imshowm(data, 20, cmap='WhBlGrYeRe', proj=f.proj)<br /> colorbar(layer)<br /> t = f.gettime(0)<br /> title('Sea ice concentration (' + t.strftime('%Y-%m-%d') + ')')

MeteoInfo

real-slimshady 发表于 2017-5-4 16:24
老师，我用这个程序可以画出日平均数据的图，但对于月平均数据会报错，同时想请问如何求文件中数据的平均 ...

参考此脚本：

1. fn = 'C:/Temp/seaice_conc_monthly_sh_f17_200801_v02r00.nc'
   
2. f = addfile(fn)
   
3. data = f['seaice_conc_monthly_cdr'][0,:,:]
   
4. data[data>1] = nan
   
5. #Plot
   
6. axesm(projinfo=f.proj, gridline=True)
   
7. lworld = shaperead('D:/Temp/map/country1.shp')
   
8. geoshow(lworld, edgecolor='k', facecolor='lightgray')
   
9. layer = imshowm(data, 20, cmap='WhBlGrYeRe', proj=f.proj)
   
10. colorbar(layer)

fn = 'C:/Temp/seaice_conc_monthly_sh_f17_200801_v02r00.nc'<br /> f = addfile(fn)<br /> data = f['seaice_conc_monthly_cdr'][0,:,:]<br /> data[data>1] = nan<br /> #Plot<br /> axesm(projinfo=f.proj, gridline=True)<br /> lworld = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(lworld, edgecolor='k', facecolor='lightgray')<br /> layer = imshowm(data, 20, cmap='WhBlGrYeRe', proj=f.proj)<br /> colorbar(layer)

微微一笑

,我按照这个脚本怎么运行不了呢，是版本太低了吗？

MeteoInfo

微微一笑 发表于 2016-4-1 10:36
,我按照这个脚本怎么运行不了呢，是版本太低了吗？

在MeteoInfoLab里运行，参考此贴：MeteoInfoLab脚本汇总贴
http://bbs.06climate.com/forum.p ... 151&fromuid=106
(出处: 气象家园)

微微一笑

谢谢老师！我在MeteoInfoLab里运行了，和您一样的脚本，就修改了路径，可是图出来不太对呢？

MeteoInfo

微微一笑 发表于 2016-4-1 15:08
谢谢老师！我在MeteoInfoLab里运行了，和您一样的脚本，就修改了路径，可是图出来不太对呢？

数据Y轴需要反向：
data = data[::-1,:]

real-slimshady

MeteoInfo 发表于 2016-4-1 16:00
数据Y轴需要反向：
data = data[::-1,:]

老师，数据中大于2.5的部分该如何处理？colorbar该如何修改成？

MeteoInfo

real-slimshady 发表于 2016-4-1 16:47
老师，数据中大于2.5的部分该如何处理？colorbar该如何修改成？

没看明白你的问题

real-slimshady

MeteoInfo 发表于 2016-4-1 17:07
没看明白你的问题

数据在0-1间，对于大于2.5的缺测数据（包括陆地数据），您是如何设为0以下的，此外如何把color bar调整成您图中的那样

MeteoInfo

real-slimshady 发表于 2016-4-1 17:16
数据在0-1间，对于大于2.5的缺测数据（包括陆地数据），您是如何设为0以下的，此外如何把color bar调整成 ...

原始数据就是这样呀，没有特殊处理。图例用cmap参数控制。

ynkmyzy

楼主
http://www.meteothinker.com