读取CloudSAT HDF Swath数据,绘图分上下两部分,上面是时间和高度维的Radar Reflectivity Factor二维图,下面是卫星轨迹图。
示例程序:# Add filef = addfile('D:/Temp/hdf/2010128055614_21420_CS_2B-GEOPROF_GRANULE_P_R04_E03.hdf')# Read datavname = '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 attributeslong_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 factorvalid_max = valid_max / scale_factorvalid_min = valid_min / scale_factordata = data / scale_factordata[data>valid_max] = nandata[data