MeteoInfoLab脚本示例：CloudSAT HDF Swath数据

MeteoInfo · 发表于 2015-10-14 23:14:05

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由 MeteoInfo 于 2017-5-11 22:39 编辑

读取CloudSAT HDF Swath数据，绘图分上下两部分，上面是时间和高度维的Radar Reflectivity Factor二维图，下面是卫星轨迹图。

示例程序：

# Add file
f = addfile('D:/Temp/hdf/2010128055614_21420_CS_2B-GEOPROF_GRANULE_P_R04_E03.hdf')
# Read data
vname = 'Radar_Reflectivity'
v_data = f[vname]
data = v_data[:,:]
v_height = f['Height']
height = v_height[0,:]
time = f['Profile_time'][:]
lon = f['Longitude'][:]
lat = f['Latitude'][:]
# Read attributes
long_name = v_data.attrvalue('long_name')[0]
scale_factor = v_data.attrvalue('factor')[0]
valid_min = v_data.attrvalue('valid_range')[0]
valid_max = v_data.attrvalue('valid_range')[1]
units = v_data.attrvalue('units')[0]
units_h = v_height.attrvalue('units')[0]
# Apply scale factor
valid_max = valid_max / scale_factor
valid_min = valid_min / scale_factor
data = data / scale_factor
data[data>valid_max] = nan
data[data<valid_min] = nan
data = transpose(data)
data = data[::-1,:]
# Make a split window plot
subplot(2, 1, 1)
# Contour the data
levs = arange(-38, 50, 2)
layer = imshow(time, height[::-1], data, levs)
colorbar(layer)
title('Radar Reflectivity Factor')
xlabel('Seconds since the start of the granule. (seconds)')
ylabel('Height (m)')
# The 2nd plot is the trajectory
subplot(2, 1, 2)
axesm(newaxes=False)
lworld = shaperead('D:/Temp/map/country1.shp')
geoshow(lworld, edgecolor='k')
plotm(lon, lat, '-b', linewidth=4)
#scatterm(lon, lat, lon, size=4, edge=False, facecolor='b')
scatterm(lon[0], lat[0], size=6, facecolor='r')
xlim(-180, 180)
ylim(-90, 90)
title('Trajectory of Flight Path (starting point in red)')

# Add file 
f = addfile('D:/Temp/hdf/2010128055614_21420_CS_2B-GEOPROF_GRANULE_P_R04_E03.hdf') 
 
# Read data 
vname = 'Radar_Reflectivity' 
v_data = f[vname] 
data = v_data[:,:] 
v_height = f['Height'] 
height = v_height[0,:] 
time = f['Profile_time'][:] 
lon = f['Longitude'][:] 
lat = f['Latitude'][:] 
 
# Read attributes 
long_name = v_data.attrvalue('long_name')[0] 
scale_factor = v_data.attrvalue('factor')[0] 
valid_min = v_data.attrvalue('valid_range')[0] 
valid_max = v_data.attrvalue('valid_range')[1] 
units = v_data.attrvalue('units')[0] 
units_h = v_height.attrvalue('units')[0] 
 
# Apply scale factor 
valid_max = valid_max / scale_factor 
valid_min = valid_min / scale_factor 
data = data / scale_factor 
data[data>valid_max] = nan 
data[data<valid_min] = nan 
data = transpose(data) 
data = data[::-1,:] 
 
# Make a split window plot 
subplot(2, 1, 1) 
 
# Contour the data 
levs = arange(-38, 50, 2) 
layer = imshow(time, height[::-1], data, levs) 
colorbar(layer) 
title('Radar Reflectivity Factor') 
xlabel('Seconds since the start of the granule. (seconds)') 
ylabel('Height (m)') 
 
# The 2nd plot is the trajectory 
subplot(2, 1, 2) 
axesm(newaxes=False) 
lworld = shaperead('D:/Temp/map/country1.shp') 
geoshow(lworld, edgecolor='k') 
plotm(lon, lat, '-b', linewidth=4) 
#scatterm(lon, lat, lon, size=4, edge=False, facecolor='b') 
scatterm(lon[0], lat[0], size=6, facecolor='r') 
xlim(-180, 180) 
ylim(-90, 90) 
title('Trajectory of Flight Path (starting point in red)')

freekiller · 发表于 2015-10-20 22:54:19

赞王老师，和之前matlab的几乎一样，功能越来越强大了。把第一幅图的横坐标设置为经纬就更好了。

大花园 · 发表于 2015-10-20 23:47:59

{:eb502:}

5xx · 发表于 2016-1-4 09:41:40

file:///C:/Users/cma5x/Desktop/捕获.PNG 请问下用这个脚本老提示这个是什么原因

MeteoInfo · 发表于 2016-1-4 09:58:34

5xx 发表于 2016-1-4 09:41
请问下用这个脚本老提示这个是什么原因

你修改过这个脚本吗？怎么修改的？

5xx · 发表于 2016-1-4 10:28:25

MeteoInfo 发表于 2016-1-4 09:58
你修改过这个脚本吗？怎么修改的？

我下载的也是这个反射率数据，他里面数据存放的格式不都应该一样的么？

MeteoInfo · 发表于 2016-1-4 10:34:22

5xx 发表于 2016-1-4 10:28
我下载的也是这个反射率数据，他里面数据存放的格式不都应该一样的么？

你运行脚本后在Console里输入：
>>> f

然后回车可以看到文件的详细信息。看看变量名还有属性之类的是不是和脚本不同？

597442388 · 发表于 2016-1-4 10:36:32

。。。。。。。。楼主

5xx · 发表于 2016-1-4 11:26:26

MeteoInfo 发表于 2016-1-4 10:34
你运行脚本后在Console里输入：
>>> f

从读取属性那里开始报错，但是我属性名都是相同的

MeteoInfo · 发表于 2016-1-4 11:30:03

5xx 发表于 2016-1-4 11:26
从读取属性那里开始报错，但是我属性名都是相同的

要不你把hdf文件作为附件附上，我下载试试。

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