MeteoInfoLab脚本示例：Hamawari-8 netCDF data

MeteoInfo

示例数据：ftp://ftp.bom.gov.au/anon/sample/catalogue/Satellite/IDE00220.201507140300.nc

该数据的分辨率很高（22000*22000），如果全部读入一个数组需要很大的内存，可以在读入时设置step（比如4，以间隔4读取数据），或者限定更小的空间范围。

脚本程序：

1. #Add data file
   
2. fn = 'D:/Temp/nc/IDE00220.201507140300.nc'
   
3. f = addfile(fn)
   
4. #Get data variable
   
5. v = f['channel_0003_brf']
   
6. data = v[0,::4,::4]
   
7. data = data[::-1,:]
   
8. #Plot
   
9. ax, proj = axesm(proj='geos', lon_0=104.7, h=35785863, gridlabel=True, gridline=True, frameon=False)
   
10. mlayer = shaperead('D:/Temp/map/country1.shp')
    
11. geoshow(mlayer)
    
12. levs = arange(0, 1, 0.1)
    
13. layer = imshowm(data, levs, proj=proj)
    
14. colorbar(layer)

#Add data file<br /> fn = 'D:/Temp/nc/IDE00220.201507140300.nc'<br /> f = addfile(fn)<br /> #Get data variable<br /> v = f['channel_0003_brf']<br /> data = v[0,::4,::4]<br /> data = data[::-1,:]<br /> #Plot<br /> ax, proj = axesm(proj='geos', lon_0=104.7, h=35785863, gridlabel=True, gridline=True, frameon=False)<br /> mlayer = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(mlayer)<br /> levs = arange(0, 1, 0.1)<br /> layer = imshowm(data, levs, proj=proj)<br /> colorbar(layer)

federa123

不明帮顶～～～～

ldepn

那个板凳坐坐

xf小飞飞侠fx

这应该是可见光通道吧，为什么我用HSD数据读出来画图默认的色标数值范围在0-2000呢？请教老师，应该如何修改？@MeteoInfo

xf小飞飞侠fx

看了您的《处理FY2E圆盘标称投影数据》帖子之后，想请教一下老师，葵花卫星的标称投影如何转换为等经纬度投影或者其他投影方式呢？为什么我处理出来的很怪？下图是提取了中国区域并转了投影的，边缘都有弯曲变形。

dymost

xf小飞飞侠fx 发表于 2017-5-16 14:42
看了您的《处理FY2E圆盘标称投影数据》帖子之后，想请教一下老师，葵花卫星的标称投影如何转换为等经纬度投 ...

您好，我也正在做葵花8卫星圆盘数据转为 经纬度坐标，能麻烦提供一份圆盘数据坐标文件吗？38298579@qq.com

xf小飞飞侠fx

dymost 发表于 2017-7-12 17:21
您好，我也正在做葵花8卫星圆盘数据转为 经纬度坐标，能麻烦提供一份圆盘数据坐标文件吗？

我这里也只有一份坛友提供的4KM分辨率的经纬度对照表，其他分辨率也没有，你需要的话我可以发给你。