MeteoInfoLab脚本示例：CloudSAT HDF Swath数据

5xx

MeteoInfo 发表于 2016-1-4 11:30
要不你把hdf文件作为附件附上，我下载试试。

我还下了您例子中的数据，只能出反射率图，轨迹图出不来，是我下载版本有问题么？

MeteoInfo

5xx 发表于 2016-1-4 12:06
我还下了您例子中的数据，只能出反射率图，轨迹图出不来，是我下载版本有问题么？

不同时期的数据的格式还是有很大的区别的，示例脚本程序中用到的数据的信息：
File Name: D:/Temp/hdf/2010128055614_21420_CS_2B-GEOPROF_GRANULE_P_R04_E03.hdf
File type: Hierarchical Data Format, version 4 (HDF4-EOS)
Dimensions: 2
        nray = 37082;
        nbin = 125;
Global Attributes:
        : HDFEOSVersion = "HDFEOS_V2.5"
        : granule_number = 21420.0f
        : ID_SITE = "Cooperative Institute for Research in the Atmosphere"
        : ID_CENTER = "CloudSat Data Processing Center"
        : ID_CENTER_URL = "http://cloudsat.cira.colostate.edu"
        : ID_MACHINE = "GILA2"
        : ID_CREATED = "Sun May  9 15:11:59 2010"
        : algorithm_version = "5.3"
        : product_version = "011"
        : start_time = "20100508055615"
        : end_time = "20100508073508"
        : _History = "Direct read of HDF4 file through CDM library; HDF-EOS StructMetadata information was read"
        : HDF4_Version = "4.1.2 (NCSA HDF post Version 4.1 Release 2, NT DLL port November 1998)"
        : featureType = "PROFILE"
        : 2B-GEOPROF_Swath_Attributes_granule_number = 21420.0f
        : 2B-GEOPROF_Swath_Attributes_ID_SITE = "Cooperative Institute for Research in the Atmosphere"
        : 2B-GEOPROF_Swath_Attributes_ID_CENTER = "CloudSat Data Processing Center"
        : 2B-GEOPROF_Swath_Attributes_ID_CENTER_URL = "http://cloudsat.cira.colostate.edu"
        : 2B-GEOPROF_Swath_Attributes_ID_MACHINE = "GILA2"
        : 2B-GEOPROF_Swath_Attributes_ID_CREATED = "Sun May  9 15:11:59 2010"
        : 2B-GEOPROF_Swath_Attributes_algorithm_version = "5.3"
        : 2B-GEOPROF_Swath_Attributes_product_version = "011"
        : 2B-GEOPROF_Swath_Attributes_start_time = "20100508055615"
        : 2B-GEOPROF_Swath_Attributes_end_time = "20100508073508"
        : 2B-GEOPROF_Swath_Attributes_Profile_time\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Profile_time\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Profile_time\.long_name = "Seconds since the start of the granule."
        : 2B-GEOPROF_Swath_Attributes_Profile_time\.units = "seconds"
        : 2B-GEOPROF_Swath_Attributes_Profile_time\.valid_range = 0.0f
        : 2B-GEOPROF_Swath_Attributes_UTC_start\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_UTC_start\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_UTC_start\.long_name = "UTC seconds since 00:00 Z of the first profile"
        : 2B-GEOPROF_Swath_Attributes_UTC_start\.units = "seconds"
        : 2B-GEOPROF_Swath_Attributes_UTC_start\.valid_range = 0.0f
        : 2B-GEOPROF_Swath_Attributes_TAI_start\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_TAI_start\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_TAI_start\.long_name = "TAI time for the first profile."
        : 2B-GEOPROF_Swath_Attributes_TAI_start\.units = "seconds"
        : 2B-GEOPROF_Swath_Attributes_TAI_start\.valid_range = 0.0
        : 2B-GEOPROF_Swath_Attributes_Latitude\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Latitude\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Latitude\.long_name = "Spacecraft Latitude"
        : 2B-GEOPROF_Swath_Attributes_Latitude\.units = "degrees"
        : 2B-GEOPROF_Swath_Attributes_Latitude\.valid_range = -90.0f
        : 2B-GEOPROF_Swath_Attributes_Longitude\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Longitude\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Longitude\.long_name = "Spacecraft Longitude"
        : 2B-GEOPROF_Swath_Attributes_Longitude\.units = "degrees"
        : 2B-GEOPROF_Swath_Attributes_Longitude\.valid_range = -180.0f
        : 2B-GEOPROF_Swath_Attributes_Height\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Height\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Height\.long_name = "Height of range bin in Reflectivity/Cloud Mask above reference surface (~ mean sea level)."
        : 2B-GEOPROF_Swath_Attributes_Height\.units = "m"
        : 2B-GEOPROF_Swath_Attributes_Height\.valid_range = -5000S
        : 2B-GEOPROF_Swath_Attributes__FV_Height = -9999S
        : 2B-GEOPROF_Swath_Attributes_Height\.missing = -9999S
        : 2B-GEOPROF_Swath_Attributes_Height\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_Range_to_intercept\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Range_to_intercept\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Range_to_intercept\.long_name = "Range from the spacecraft to the CPR boresight intercept with the geoid"
        : 2B-GEOPROF_Swath_Attributes_Range_to_intercept\.units = "km"
        : 2B-GEOPROF_Swath_Attributes_Range_to_intercept\.valid_range = 600.0f
        : 2B-GEOPROF_Swath_Attributes__FV_DEM_elevation = -9999S
        : 2B-GEOPROF_Swath_Attributes_DEM_elevation\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_DEM_elevation\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_DEM_elevation\.long_name = "Digital Elevation Map"
        : 2B-GEOPROF_Swath_Attributes_DEM_elevation\.units = "meters"
        : 2B-GEOPROF_Swath_Attributes_DEM_elevation\.valid_range = -9999S
        : 2B-GEOPROF_Swath_Attributes_DEM_elevation\.missing = 9999S
        : 2B-GEOPROF_Swath_Attributes_DEM_elevation\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_Vertical_binsize\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Vertical_binsize\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Vertical_binsize\.units = "m"
        : 2B-GEOPROF_Swath_Attributes_Vertical_binsize\.missing = -9999.0f
        : 2B-GEOPROF_Swath_Attributes_Vertical_binsize\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_Pitch_offset\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Pitch_offset\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Pitch_offset\.long_name = "Nominal satellite pitch angle offset from nadir"
        : 2B-GEOPROF_Swath_Attributes_Pitch_offset\.units = "degrees"
        : 2B-GEOPROF_Swath_Attributes_Pitch_offset\.valid_range = -90.0f
        : 2B-GEOPROF_Swath_Attributes_Roll_offset\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Roll_offset\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Roll_offset\.long_name = "Nominal satellite roll angle offset from nadir"
        : 2B-GEOPROF_Swath_Attributes_Roll_offset\.units = "degrees"
        : 2B-GEOPROF_Swath_Attributes_Roll_offset\.valid_range = -90.0f
        : 2B-GEOPROF_Swath_Attributes__FV_Data_quality = 0B
        : 2B-GEOPROF_Swath_Attributes_Data_quality\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Data_quality\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Data_quality\.long_name = "Data Quality"
        : 2B-GEOPROF_Swath_Attributes_Data_quality\.units = "--"
        : 2B-GEOPROF_Swath_Attributes_Data_quality\.valid_range = 0B
        : 2B-GEOPROF_Swath_Attributes__FV_Data_status = 48B
        : 2B-GEOPROF_Swath_Attributes_Data_status\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Data_status\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Data_status\.long_name = "Data status flags"
        : 2B-GEOPROF_Swath_Attributes_Data_status\.units = "--"
        : 2B-GEOPROF_Swath_Attributes_Data_status\.valid_range = 0B
        : 2B-GEOPROF_Swath_Attributes__FV_Data_targetID = 0B
        : 2B-GEOPROF_Swath_Attributes_Data_targetID\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Data_targetID\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Data_targetID\.long_name = "CPR bus orientation (target ID)"
        : 2B-GEOPROF_Swath_Attributes_Data_targetID\.units = "--"
        : 2B-GEOPROF_Swath_Attributes_Data_targetID\.valid_range = 0B
        : 2B-GEOPROF_Swath_Attributes_SurfaceHeightBin\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_SurfaceHeightBin\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_SurfaceHeightBin\.long_name = "Location of Surface Bin as determined by 1B CPR algorithm.  The value here is shifted (as Height)."
        : 2B-GEOPROF_Swath_Attributes_SurfaceHeightBin\.valid_range = 1B
        : 2B-GEOPROF_Swath_Attributes_SurfaceHeightBin\.missing = -1B
        : 2B-GEOPROF_Swath_Attributes_SurfaceHeightBin\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_SurfaceHeightBin_fraction\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_SurfaceHeightBin_fraction\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_SurfaceHeightBin_fraction\.long_name = "SurfaceHeightBin_fraction"
        : 2B-GEOPROF_Swath_Attributes_SurfaceHeightBin_fraction\.units = "none"
        : 2B-GEOPROF_Swath_Attributes__FV_SurfaceHeightBin_fraction = 0.0f
        : 2B-GEOPROF_Swath_Attributes_SurfaceHeightBin_fraction\.missing = 0.0f
        : 2B-GEOPROF_Swath_Attributes_CPR_Cloud_mask\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_CPR_Cloud_mask\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_CPR_Cloud_mask\.long_name = "CPR Cloud Mask"
        : 2B-GEOPROF_Swath_Attributes_CPR_Cloud_mask\.valid_range = 0B
        : 2B-GEOPROF_Swath_Attributes__FV_CPR_Cloud_mask = -99B
        : 2B-GEOPROF_Swath_Attributes_CPR_Cloud_mask\.missing = -9B
        : 2B-GEOPROF_Swath_Attributes_CPR_Cloud_mask\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_Gaseous_Attenuation\.factor = 100.0f
        : 2B-GEOPROF_Swath_Attributes_Gaseous_Attenuation\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Gaseous_Attenuation\.long_name = "Gaseous_Attenuation"
        : 2B-GEOPROF_Swath_Attributes_Gaseous_Attenuation\.units = "dBZe"
        : 2B-GEOPROF_Swath_Attributes_Gaseous_Attenuation\.valid_range = 0S
        : 2B-GEOPROF_Swath_Attributes__FV_Gaseous_Attenuation = 15360S
        : 2B-GEOPROF_Swath_Attributes_Gaseous_Attenuation\.missing = -9999S
        : 2B-GEOPROF_Swath_Attributes_Gaseous_Attenuation\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_Radar_Reflectivity\.factor = 100.0f
        : 2B-GEOPROF_Swath_Attributes_Radar_Reflectivity\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Radar_Reflectivity\.long_name = "Radar Reflectivity Factor"
        : 2B-GEOPROF_Swath_Attributes_Radar_Reflectivity\.units = "dBZe"
        : 2B-GEOPROF_Swath_Attributes_Radar_Reflectivity\.valid_range = -4000S
        : 2B-GEOPROF_Swath_Attributes__FV_Radar_Reflectivity = 15360S
        : 2B-GEOPROF_Swath_Attributes_Radar_Reflectivity\.missing = -8888S
        : 2B-GEOPROF_Swath_Attributes_Radar_Reflectivity\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes__FV_Sigma-Zero = -9999S
        : 2B-GEOPROF_Swath_Attributes_Sigma-Zero\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Sigma-Zero\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Sigma-Zero\.long_name = "Sigma-Zero"
        : 2B-GEOPROF_Swath_Attributes_Sigma-Zero\.units = "dB*100"
        : 2B-GEOPROF_Swath_Attributes_Sigma-Zero\.valid_range = -1000S
        : 2B-GEOPROF_Swath_Attributes_Sigma-Zero\.missing = -9999S
        : 2B-GEOPROF_Swath_Attributes_Sigma-Zero\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_MODIS_cloud_flag\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_MODIS_cloud_flag\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_MODIS_cloud_flag\.long_name = "MOD35_bit_2and3_cloud_flag"
        : 2B-GEOPROF_Swath_Attributes_MODIS_cloud_flag\.units = "None"
        : 2B-GEOPROF_Swath_Attributes_MODIS_cloud_flag\.valid_range = 0B
        : 2B-GEOPROF_Swath_Attributes__FV_MODIS_cloud_flag = 99B
        : 2B-GEOPROF_Swath_Attributes_MODIS_cloud_flag\.missing = 99B
        : 2B-GEOPROF_Swath_Attributes_MODIS_cloud_flag\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_MODIS_Cloud_Fraction\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_MODIS_Cloud_Fraction\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_MODIS_Cloud_Fraction\.long_name = "MODIS 250m Cloud Fraction"
        : 2B-GEOPROF_Swath_Attributes_MODIS_Cloud_Fraction\.valid_range = 0B
        : 2B-GEOPROF_Swath_Attributes__FV_MODIS_Cloud_Fraction = -99B
        : 2B-GEOPROF_Swath_Attributes_MODIS_Cloud_Fraction\.missing = -99B
        : 2B-GEOPROF_Swath_Attributes_MODIS_Cloud_Fraction\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_MODIS_scene_char\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_MODIS_scene_char\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_MODIS_scene_char\.long_name = "MODIS scene characterizations"
        : 2B-GEOPROF_Swath_Attributes_MODIS_scene_char\.valid_range = 0B
        : 2B-GEOPROF_Swath_Attributes__FV_MODIS_scene_char = -99B
        : 2B-GEOPROF_Swath_Attributes_MODIS_scene_char\.missing = -9B
        : 2B-GEOPROF_Swath_Attributes_MODIS_scene_char\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_MODIS_scene_var\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_MODIS_scene_var\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_MODIS_scene_var\.long_name = "MODIS scene variability"
        : 2B-GEOPROF_Swath_Attributes_MODIS_scene_var\.valid_range = 0B
        : 2B-GEOPROF_Swath_Attributes__FV_MODIS_scene_var = -99B
        : 2B-GEOPROF_Swath_Attributes_MODIS_scene_var\.missing = -9B
        : 2B-GEOPROF_Swath_Attributes_MODIS_scene_var\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_CPR_Echo_Top\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_CPR_Echo_Top\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_CPR_Echo_Top\.long_name = "CPR Echo Top Characterizations"
        : 2B-GEOPROF_Swath_Attributes_CPR_Echo_Top\.valid_range = 0B
        : 2B-GEOPROF_Swath_Attributes__FV_CPR_Echo_Top = -99B
        : 2B-GEOPROF_Swath_Attributes_CPR_Echo_Top\.missing = -9B
        : 2B-GEOPROF_Swath_Attributes_CPR_Echo_Top\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseFloor\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseFloor\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseFloor\.long_name = "Noise Floor calculated in Cloud_Masking Routine"
        : 2B-GEOPROF_Swath_Attributes__FV_sem_NoiseFloor = 0.0f
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseFloor\.missing = 0.0f
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseFloor\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseFloorVar\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseFloorVar\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseFloorVar\.long_name = "Noise Variance in measured return (i.e. Echo) power"
        : 2B-GEOPROF_Swath_Attributes__FV_sem_NoiseFloorVar = 0.0f
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseFloorVar\.missing = 0.0f
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseFloorVar\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseGate\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseGate\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseGate\.long_name = "Noise Floor and Noise Variance estimate  based on window of a fixed size centered at this range bin."
        : 2B-GEOPROF_Swath_Attributes__FV_sem_NoiseGate = 0B
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseGate\.missing = 0B
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseGate\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes__FV_Navigation_land_sea_flag = 2B
        : 2B-GEOPROF_Swath_Attributes_Navigation_land_sea_flag\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Navigation_land_sea_flag\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Navigation_land_sea_flag\.long_name = "Land Sea Flag"
        : 2B-GEOPROF_Swath_Attributes_Navigation_land_sea_flag\.units = "--"
        : 2B-GEOPROF_Swath_Attributes_Navigation_land_sea_flag\.valid_range = 1B
        : 2B-GEOPROF_Swath_Attributes_Clutter_reduction_flag\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Clutter_reduction_flag\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Clutter_reduction_flag\.long_name = "Clutter_reduction_flag"
Variations: 31
        char StructMetadata.0(null);
        float Profile_time(nray);
        float UTC_start);
        double TAI_start);
        float Latitude(nray);
                Latitude: _CoordinateAxisType = "Lat"
                Latitude: units = "degrees_north"
        float Longitude(nray);
                Longitude: _CoordinateAxisType = "Lon"
                Longitude: units = "degrees_east"
        short Height(nray,nbin);
                Height: _FillValue = -9999S
                Height: factor = 1.0f
                Height: offset = 0.0f
                Height: long_name = "Height of range bin in Reflectivity/Cloud Mask above reference surface (~ mean sea level)."
                Height: units = "m"
                Height: valid_range = -5000S, 30000S
                Height: missing = -9999S
                Height: missop = "=="
        float Range_to_intercept(nray);
        short DEM_elevation(nray);
        float Vertical_binsize);
        float Pitch_offset);
        float Roll_offset);
        byte Data_quality(nray);
                Data_quality: _Unsigned = "true"
        byte Data_status(nray);
                Data_status: _Unsigned = "true"
        byte Data_targetID(nray);
                Data_targetID: _Unsigned = "true"
        byte SurfaceHeightBin(nray);
        float SurfaceHeightBin_fraction(nray);
        byte CPR_Cloud_mask(nray,nbin);
                CPR_Cloud_mask: _FillValue = -99B
                CPR_Cloud_mask: factor = 1.0f
                CPR_Cloud_mask: offset = 0.0f
                CPR_Cloud_mask: long_name = "CPR Cloud Mask"
                CPR_Cloud_mask: valid_range = 0B, 40B
                CPR_Cloud_mask: missing = -9B
                CPR_Cloud_mask: missop = "=="
        short Gaseous_Attenuation(nray,nbin);
                Gaseous_Attenuation: _FillValue = 15360S
                Gaseous_Attenuation: factor = 100.0f
                Gaseous_Attenuation: offset = 0.0f
                Gaseous_Attenuation: long_name = "Gaseous_Attenuation"
                Gaseous_Attenuation: units = "dBZe"
                Gaseous_Attenuation: valid_range = 0S, 1000S
                Gaseous_Attenuation: missing = -9999S
                Gaseous_Attenuation: missop = "=="
        short Radar_Reflectivity(nray,nbin);
                Radar_Reflectivity: _FillValue = 15360S
                Radar_Reflectivity: factor = 100.0f
                Radar_Reflectivity: offset = 0.0f
                Radar_Reflectivity: long_name = "Radar Reflectivity Factor"
                Radar_Reflectivity: units = "dBZe"
                Radar_Reflectivity: valid_range = -4000S, 5000S
                Radar_Reflectivity: missing = -8888S
                Radar_Reflectivity: missop = "=="
        short Sigma-Zero(nray);
        byte MODIS_cloud_flag(nray);
        byte MODIS_Cloud_Fraction(nray);
        byte MODIS_scene_char(nray);
        byte MODIS_scene_var(nray);
        byte CPR_Echo_Top(nray);
        float sem_NoiseFloor(nray);
        float sem_NoiseFloorVar(nray);
        byte sem_NoiseGate(nray);
        byte Navigation_land_sea_flag(nray);
                Navigation_land_sea_flag: _Unsigned = "true"
        byte Clutter_reduction_flag(nray);


你提供的数据的信息：
File Name: D:/Temp/hdf/2015216230303_49308_CS_2B-GEOPROF_GRANULE_P_R04_E06.hdf
File type: Hierarchical Data Format, version 4 (HDF4-EOS)
Dimensions: 2
        nray = 20679;
        nbin = 125;
Global Attributes:
        : HDFEOSVersion = "HDFEOS_V2.5"
        : _History = "Direct read of HDF4 file through CDM library; HDF-EOS StructMetadata information was read"
        : HDF4_Version = "4.2.1 (NCSA HDF Version 4.2 Release 1, February 17, 2005)"
        : featureType = "PROFILE"
        : 2B-GEOPROF_Swath_Attributes_granule_number = 49308.0f
        : 2B-GEOPROF_Swath_Attributes_ID_SITE = "Cooperative Institute for Research in the Atmosphere"
        : 2B-GEOPROF_Swath_Attributes_ID_CENTER = "CloudSat Data Processing Center"
        : 2B-GEOPROF_Swath_Attributes_ID_CENTER_URL = "http://cloudsat.cira.colostate.edu"
        : 2B-GEOPROF_Swath_Attributes_ID_MACHINE = "VS15"
        : 2B-GEOPROF_Swath_Attributes_ID_CREATED = "Fri Sep 11 13:07:27 2015"
        : 2B-GEOPROF_Swath_Attributes_algorithm_name = "2B-GEOPROF"
        : 2B-GEOPROF_Swath_Attributes_algorithm_version = "5.3"
        : 2B-GEOPROF_Swath_Attributes_product_version = "011"
        : 2B-GEOPROF_Swath_Attributes_start_time = "20150804233503"
        : 2B-GEOPROF_Swath_Attributes_end_time = "20150805003011"
        : 2B-GEOPROF_Swath_Attributes_Profile_time\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Profile_time\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Profile_time\.long_name = "Seconds since the start of the granule."
        : 2B-GEOPROF_Swath_Attributes_Profile_time\.units = "seconds"
        : 2B-GEOPROF_Swath_Attributes_Profile_time\.valid_range = 0.0f
        : 2B-GEOPROF_Swath_Attributes_UTC_start\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_UTC_start\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_UTC_start\.long_name = "UTC seconds since 00:00 Z of the first profile"
        : 2B-GEOPROF_Swath_Attributes_UTC_start\.units = "seconds"
        : 2B-GEOPROF_Swath_Attributes_UTC_start\.valid_range = 0.0f
        : 2B-GEOPROF_Swath_Attributes_TAI_start\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_TAI_start\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_TAI_start\.long_name = "TAI time for the first profile."
        : 2B-GEOPROF_Swath_Attributes_TAI_start\.units = "seconds"
        : 2B-GEOPROF_Swath_Attributes_TAI_start\.valid_range = 0.0
        : 2B-GEOPROF_Swath_Attributes_Latitude\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Latitude\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Latitude\.long_name = "Spacecraft Latitude"
        : 2B-GEOPROF_Swath_Attributes_Latitude\.units = "degrees"
        : 2B-GEOPROF_Swath_Attributes_Latitude\.valid_range = -90.0f
        : 2B-GEOPROF_Swath_Attributes_Longitude\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Longitude\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Longitude\.long_name = "Spacecraft Longitude"
        : 2B-GEOPROF_Swath_Attributes_Longitude\.units = "degrees"
        : 2B-GEOPROF_Swath_Attributes_Longitude\.valid_range = -180.0f
        : 2B-GEOPROF_Swath_Attributes_Height\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Height\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Height\.long_name = "Height of range bin in Reflectivity/Cloud Mask above reference surface (~ mean sea level)."
        : 2B-GEOPROF_Swath_Attributes_Height\.units = "m"
        : 2B-GEOPROF_Swath_Attributes_Height\.valid_range = -5000S
        : 2B-GEOPROF_Swath_Attributes__FV_Height = -9999S
        : 2B-GEOPROF_Swath_Attributes_Height\.missing = -9999S
        : 2B-GEOPROF_Swath_Attributes_Height\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_Range_to_intercept\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Range_to_intercept\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Range_to_intercept\.long_name = "Range from the spacecraft to the CPR boresight intercept with the geoid"
        : 2B-GEOPROF_Swath_Attributes_Range_to_intercept\.units = "km"
        : 2B-GEOPROF_Swath_Attributes_Range_to_intercept\.valid_range = 600.0f
        : 2B-GEOPROF_Swath_Attributes__FV_DEM_elevation = -9999S
        : 2B-GEOPROF_Swath_Attributes_DEM_elevation\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_DEM_elevation\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_DEM_elevation\.long_name = "Digital Elevation Map"
        : 2B-GEOPROF_Swath_Attributes_DEM_elevation\.units = "meters"
        : 2B-GEOPROF_Swath_Attributes_DEM_elevation\.valid_range = -9999S
        : 2B-GEOPROF_Swath_Attributes_DEM_elevation\.missing = 9999S
        : 2B-GEOPROF_Swath_Attributes_DEM_elevation\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_Vertical_binsize\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Vertical_binsize\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Vertical_binsize\.units = "m"
        : 2B-GEOPROF_Swath_Attributes_Vertical_binsize\.missing = -9999.0f
        : 2B-GEOPROF_Swath_Attributes_Vertical_binsize\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_Pitch_offset\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Pitch_offset\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Pitch_offset\.long_name = "Nominal satellite pitch angle offset from nadir"
        : 2B-GEOPROF_Swath_Attributes_Pitch_offset\.units = "degrees"
        : 2B-GEOPROF_Swath_Attributes_Pitch_offset\.valid_range = -90.0f
        : 2B-GEOPROF_Swath_Attributes_Roll_offset\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Roll_offset\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Roll_offset\.long_name = "Nominal satellite roll angle offset from nadir"
        : 2B-GEOPROF_Swath_Attributes_Roll_offset\.units = "degrees"
        : 2B-GEOPROF_Swath_Attributes_Roll_offset\.valid_range = -90.0f
        : 2B-GEOPROF_Swath_Attributes__FV_Data_quality = 0B
        : 2B-GEOPROF_Swath_Attributes_Data_quality\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Data_quality\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Data_quality\.long_name = "Data Quality"
        : 2B-GEOPROF_Swath_Attributes_Data_quality\.units = "--"
        : 2B-GEOPROF_Swath_Attributes_Data_quality\.valid_range = 0B
        : 2B-GEOPROF_Swath_Attributes__FV_Data_status = 48B
        : 2B-GEOPROF_Swath_Attributes_Data_status\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Data_status\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Data_status\.long_name = "Data status flags"
        : 2B-GEOPROF_Swath_Attributes_Data_status\.units = "--"
        : 2B-GEOPROF_Swath_Attributes_Data_status\.valid_range = 0B
        : 2B-GEOPROF_Swath_Attributes__FV_Data_targetID = 0B
        : 2B-GEOPROF_Swath_Attributes_Data_targetID\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Data_targetID\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Data_targetID\.long_name = "CPR bus orientation (target ID)"
        : 2B-GEOPROF_Swath_Attributes_Data_targetID\.units = "--"
        : 2B-GEOPROF_Swath_Attributes_Data_targetID\.valid_range = 0B
        : 2B-GEOPROF_Swath_Attributes_SurfaceHeightBin\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_SurfaceHeightBin\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_SurfaceHeightBin\.long_name = "Location of Surface Bin as determined by 1B CPR algorithm.  The value here is shifted (as Height)."
        : 2B-GEOPROF_Swath_Attributes_SurfaceHeightBin\.valid_range = 1B
        : 2B-GEOPROF_Swath_Attributes_SurfaceHeightBin\.missing = -1B
        : 2B-GEOPROF_Swath_Attributes_SurfaceHeightBin\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_SurfaceHeightBin_fraction\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_SurfaceHeightBin_fraction\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_SurfaceHeightBin_fraction\.long_name = "SurfaceHeightBin_fraction"
        : 2B-GEOPROF_Swath_Attributes_SurfaceHeightBin_fraction\.units = "none"
        : 2B-GEOPROF_Swath_Attributes__FV_SurfaceHeightBin_fraction = 0.0f
        : 2B-GEOPROF_Swath_Attributes_SurfaceHeightBin_fraction\.missing = 0.0f
        : 2B-GEOPROF_Swath_Attributes_CPR_Cloud_mask\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_CPR_Cloud_mask\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_CPR_Cloud_mask\.long_name = "CPR Cloud Mask"
        : 2B-GEOPROF_Swath_Attributes_CPR_Cloud_mask\.valid_range = 0B
        : 2B-GEOPROF_Swath_Attributes__FV_CPR_Cloud_mask = -99B
        : 2B-GEOPROF_Swath_Attributes_CPR_Cloud_mask\.missing = -9B
        : 2B-GEOPROF_Swath_Attributes_CPR_Cloud_mask\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_Gaseous_Attenuation\.factor = 100.0f
        : 2B-GEOPROF_Swath_Attributes_Gaseous_Attenuation\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Gaseous_Attenuation\.long_name = "Gaseous_Attenuation"
        : 2B-GEOPROF_Swath_Attributes_Gaseous_Attenuation\.units = "dBZe"
        : 2B-GEOPROF_Swath_Attributes_Gaseous_Attenuation\.valid_range = 0S
        : 2B-GEOPROF_Swath_Attributes__FV_Gaseous_Attenuation = 15360S
        : 2B-GEOPROF_Swath_Attributes_Gaseous_Attenuation\.missing = -9999S
        : 2B-GEOPROF_Swath_Attributes_Gaseous_Attenuation\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_Radar_Reflectivity\.factor = 100.0f
        : 2B-GEOPROF_Swath_Attributes_Radar_Reflectivity\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Radar_Reflectivity\.long_name = "Radar Reflectivity Factor"
        : 2B-GEOPROF_Swath_Attributes_Radar_Reflectivity\.units = "dBZe"
        : 2B-GEOPROF_Swath_Attributes_Radar_Reflectivity\.valid_range = -4000S
        : 2B-GEOPROF_Swath_Attributes__FV_Radar_Reflectivity = 15360S
        : 2B-GEOPROF_Swath_Attributes_Radar_Reflectivity\.missing = -8888S
        : 2B-GEOPROF_Swath_Attributes_Radar_Reflectivity\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes__FV_Sigma-Zero = -9999S
        : 2B-GEOPROF_Swath_Attributes_Sigma-Zero\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Sigma-Zero\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Sigma-Zero\.long_name = "Sigma-Zero"
        : 2B-GEOPROF_Swath_Attributes_Sigma-Zero\.units = "dB*100"
        : 2B-GEOPROF_Swath_Attributes_Sigma-Zero\.valid_range = -1000S
        : 2B-GEOPROF_Swath_Attributes_Sigma-Zero\.missing = -9999S
        : 2B-GEOPROF_Swath_Attributes_Sigma-Zero\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_MODIS_cloud_flag\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_MODIS_cloud_flag\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_MODIS_cloud_flag\.long_name = "MOD35_bit_2and3_cloud_flag"
        : 2B-GEOPROF_Swath_Attributes_MODIS_cloud_flag\.units = "None"
        : 2B-GEOPROF_Swath_Attributes_MODIS_cloud_flag\.valid_range = 0B
        : 2B-GEOPROF_Swath_Attributes__FV_MODIS_cloud_flag = 99B
        : 2B-GEOPROF_Swath_Attributes_MODIS_cloud_flag\.missing = 99B
        : 2B-GEOPROF_Swath_Attributes_MODIS_cloud_flag\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_MODIS_Cloud_Fraction\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_MODIS_Cloud_Fraction\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_MODIS_Cloud_Fraction\.long_name = "MODIS 250m Cloud Fraction"
        : 2B-GEOPROF_Swath_Attributes_MODIS_Cloud_Fraction\.valid_range = 0B
        : 2B-GEOPROF_Swath_Attributes__FV_MODIS_Cloud_Fraction = -99B
        : 2B-GEOPROF_Swath_Attributes_MODIS_Cloud_Fraction\.missing = -99B
        : 2B-GEOPROF_Swath_Attributes_MODIS_Cloud_Fraction\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_MODIS_scene_char\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_MODIS_scene_char\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_MODIS_scene_char\.long_name = "MODIS scene characterizations"
        : 2B-GEOPROF_Swath_Attributes_MODIS_scene_char\.valid_range = 0B
        : 2B-GEOPROF_Swath_Attributes__FV_MODIS_scene_char = -99B
        : 2B-GEOPROF_Swath_Attributes_MODIS_scene_char\.missing = -9B
        : 2B-GEOPROF_Swath_Attributes_MODIS_scene_char\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_MODIS_scene_var\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_MODIS_scene_var\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_MODIS_scene_var\.long_name = "MODIS scene variability"
        : 2B-GEOPROF_Swath_Attributes_MODIS_scene_var\.valid_range = 0B
        : 2B-GEOPROF_Swath_Attributes__FV_MODIS_scene_var = -99B
        : 2B-GEOPROF_Swath_Attributes_MODIS_scene_var\.missing = -9B
        : 2B-GEOPROF_Swath_Attributes_MODIS_scene_var\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_CPR_Echo_Top\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_CPR_Echo_Top\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_CPR_Echo_Top\.long_name = "CPR Echo Top Characterizations"
        : 2B-GEOPROF_Swath_Attributes_CPR_Echo_Top\.valid_range = 0B
        : 2B-GEOPROF_Swath_Attributes__FV_CPR_Echo_Top = -99B
        : 2B-GEOPROF_Swath_Attributes_CPR_Echo_Top\.missing = -9B
        : 2B-GEOPROF_Swath_Attributes_CPR_Echo_Top\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseFloor\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseFloor\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseFloor\.long_name = "Noise Floor calculated in Cloud_Masking Routine"
        : 2B-GEOPROF_Swath_Attributes__FV_sem_NoiseFloor = 0.0f
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseFloor\.missing = 0.0f
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseFloor\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseFloorVar\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseFloorVar\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseFloorVar\.long_name = "Noise Variance in measured return (i.e. Echo) power"
        : 2B-GEOPROF_Swath_Attributes__FV_sem_NoiseFloorVar = 0.0f
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseFloorVar\.missing = 0.0f
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseFloorVar\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseGate\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseGate\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseGate\.long_name = "Noise Floor and Noise Variance estimate  based on window of a fixed size centered at this range bin."
        : 2B-GEOPROF_Swath_Attributes__FV_sem_NoiseGate = 0B
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseGate\.missing = 0B
        : 2B-GEOPROF_Swath_Attributes_sem_NoiseGate\.missop = "=="
        : 2B-GEOPROF_Swath_Attributes__FV_Navigation_land_sea_flag = 2B
        : 2B-GEOPROF_Swath_Attributes_Navigation_land_sea_flag\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Navigation_land_sea_flag\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Navigation_land_sea_flag\.long_name = "Land Sea Flag"
        : 2B-GEOPROF_Swath_Attributes_Navigation_land_sea_flag\.units = "--"
        : 2B-GEOPROF_Swath_Attributes_Navigation_land_sea_flag\.valid_range = 1B
        : 2B-GEOPROF_Swath_Attributes_Clutter_reduction_flag\.factor = 1.0f
        : 2B-GEOPROF_Swath_Attributes_Clutter_reduction_flag\.offset = 0.0f
        : 2B-GEOPROF_Swath_Attributes_Clutter_reduction_flag\.long_name = "Clutter_reduction_flag"
Variations: 31
        char StructMetadata.0(null);
        float Profile_time(nray);
        float UTC_start);
        double TAI_start);
        float Latitude(nray);
                Latitude: _CoordinateAxisType = "Lat"
                Latitude: units = "degrees_north"
        float Longitude(nray);
                Longitude: _CoordinateAxisType = "Lon"
                Longitude: units = "degrees_east"
        short Height(nray,nbin);
                Height: _FillValue = -9999S
        float Range_to_intercept(nray);
        short DEM_elevation(nray);
        float Vertical_binsize);
        float Pitch_offset);
        float Roll_offset);
        byte Data_quality(nray);
                Data_quality: _Unsigned = "true"
        byte Data_status(nray);
                Data_status: _Unsigned = "true"
        byte Data_targetID(nray);
                Data_targetID: _Unsigned = "true"
        byte SurfaceHeightBin(nray);
        float SurfaceHeightBin_fraction(nray);
        byte CPR_Cloud_mask(nray,nbin);
                CPR_Cloud_mask: _FillValue = -99B
        short Gaseous_Attenuation(nray,nbin);
                Gaseous_Attenuation: _FillValue = 15360S
        short Radar_Reflectivity(nray,nbin);
                Radar_Reflectivity: _FillValue = 15360S
        short Sigma-Zero(nray);
        byte MODIS_cloud_flag(nray);
        byte MODIS_Cloud_Fraction(nray);
        byte MODIS_scene_char(nray);
        byte MODIS_scene_var(nray);
        byte CPR_Echo_Top(nray);
        float sem_NoiseFloor(nray);
        float sem_NoiseFloorVar(nray);
        byte sem_NoiseGate(nray);
        byte Navigation_land_sea_flag(nray);
                Navigation_land_sea_flag: _Unsigned = "true"
        byte Clutter_reduction_flag(nray);

MeteoInfo

5xx 发表于 2016-1-4 12:06
我还下了您例子中的数据，只能出反射率图，轨迹图出不来，是我下载版本有问题么？

针对你的数据修改脚本如下：

1. # Add file
   
2. f = addfile('D:/Temp/hdf/2015216230303_49308_CS_2B-GEOPROF_GRANULE_P_R04_E06.hdf')
   
3. 
   
4. # Read data
   
5. vname = 'Radar_Reflectivity'
   
6. v_data = f[vname]
   
7. data = v_data[:,:]
   
8. v_height = f['Height']
   
9. height = v_height[0,:]
   
10. time = f['Profile_time'][:]
    
11. lon = f['Longitude'][:]
    
12. lat = f['Latitude'][:]
    
13. 
    
14. # Read attributes
    
15. long_name = f.attrvalue('2B-GEOPROF_Swath_Attributes_Radar_Reflectivity\.long_name')[0]
    
16. scale_factor = f.attrvalue('2B-GEOPROF_Swath_Attributes_Radar_Reflectivity\.factor')[0]
    
17. valid_min = f.attrvalue('2B-GEOPROF_Swath_Attributes_Radar_Reflectivity\.valid_range')[0]
    
18. valid_max = 5000
    
19. fill_value = v_data.attrvalue('_FillValue')[0]
    
20. units = f.attrvalue('2B-GEOPROF_Swath_Attributes_Radar_Reflectivity\.units')[0]
    
21. units_h = f.attrvalue('2B-GEOPROF_Swath_Attributes_Height\.units')[0]
    
22. 
    
23. # Apply scale factor
    
24. valid_max = valid_max / scale_factor
    
25. valid_min = valid_min / scale_factor
    
26. fill_value = fill_value / scale_factor
    
27. data = data / scale_factor
    
28. data[data>valid_max] = nan
    
29. data[data<valid_min] = nan
    
30. #data[data==fill_value] = nan
    
31. data = transpose(data)
    
32. data = data[::-1,:]
    
33. 
    
34. # Make a split window plot
    
35. subplot(2, 1, 1)
    
36. 
    
37. # Contour the data
    
38. levs = arange(-38, 50, 2)
    
39. layer = imshow(time, height[::-1], data, levs)
    
40. colorbar(layer)
    
41. title('Radar Reflectivity Factor')
    
42. xlabel('Seconds since the start of the granule. (seconds)')
    
43. ylabel('Height (m)')
    
44. 
    
45. # The 2nd plot is the trajectory
    
46. subplot(2, 1, 2)
    
47. axesm()
    
48. lworld = shaperead('D:/Temp/map/country1.shp')
    
49. geoshow(lworld, edgecolor='k')
    
50. plotm(lon, lat, '-b', linewidth=4)
    
51. #scatterm(lon, lat, lon, size=4, edge=False, facecolor='b')
    
52. scatterm(lon[0], lat[0], size=6, facecolor='r')
    
53. xlim(-180, 180)
    
54. ylim(-90, 90)
    
55. title('Trajectory of Flight Path (starting point in red)')

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

5xx

MeteoInfo 发表于 2016-1-4 13:03
针对你的数据修改脚本如下：

谢谢老师指导，请问下就是这个属性名称写错了么，我把您新加的属性fill-value消去之后也能出图。

MeteoInfo

5xx 发表于 2016-1-4 14:10
谢谢老师指导，请问下就是这个属性名称写错了么，我把您新加的属性fill-value消去之后也能出图。

应该没有写错，不过fill_value实际上没有用到，注释掉也可以。

5xx

MeteoInfo 发表于 2016-1-4 14:39
应该没有写错，不过fill_value实际上没有用到，注释掉也可以。

哦哦，谢谢老师，请问下脚本里面data = data[::-1,:]和 height[::-1]分别代表什么啊~

MeteoInfo

5xx 发表于 2016-1-4 18:47
哦哦，谢谢老师，请问下脚本里面data = data[::-1,:]和 height[::-1]分别代表什么啊~

::-1 指的是该维反向。

yuyuxuan

您好  请问一下  怎么才可以把出的图中 后面那些 蓝色的花的换成白色的  可以换吗

a930013

謝謝大大分享!!