Print all MODIS C6 MYD03 SDS names and attributes in python 2.7

Outputs of a code (see bottom of the page) written in python 2.7 to print all MODIS C6 MYD03 SDS names and attributes in python 2.7. All SDS infos (nd: number of dimensions, t: type, na: number of attributes):

Latitude

units degrees
_FillValue 999.0
valid_range [90.0, 90.0]

Longitude

units degrees
_FillValue 999.0
valid_range [180.0, 180.0]

Scan Offset

units km IFOV
_FillValue 128
scale_factor 0.006
valid_range [127, 127]

Track Offset

units km IFOV
_FillValue 128
scale_factor 0.006
valid_range [127, 127]

Height Offset

units km
_FillValue 128
scale_factor 0.006
valid_range [127, 127]

Height

units meters
_FillValue 32767
valid_range [400, 10000]

SensorZenith

units degrees
_FillValue 32767
scale_factor 0.01
valid_range [0, 18000]

SensorAzimuth

units degrees
_FillValue 32767
scale_factor 0.01
valid_range [18000, 18000]

Range

units meters
_FillValue 0
scale_factor 25.0
valid_range [27000, 65535]

SolarZenith

units degrees
_FillValue 32767
scale_factor 0.01
valid_range [0, 18000]

SolarAzimuth

units degrees
_FillValue 32767
scale_factor 0.01
valid_range [18000, 18000]

Land/SeaMask

_FillValue 221
valid_range [0, 7]

WaterPresent

_FillValue 255
valid_range [0, 8]

gflags

_FillValue 255

Scan number

_FillValue 0

EV frames

valid_range [0, 1400]

SD frames

SV frames

EV start time

units seconds
_FillValue 2000000000.0

SD start time

units seconds
_FillValue 2000000000.0

SV start time

units seconds
_FillValue 2000000000.0

EV center time

units seconds
_FillValue 2000000000.0

Mirror side

_FillValue 65535
valid_range [0, 1]

SD Sun zenith

units radians
_FillValue 9.96920996839e+36

SD Sun azimuth

units radians
_FillValue 9.96920996839e+36

Moon Vector

_FillValue 9.96920996839e+36

L1 scan quality

_FillValue 1

Geo scan quality

_FillValue 127

orb_pos

units meters
_FillValue 9.96920996839e+36
valid_range [7200000.0, 7200000.0]

orb_vel

units meters per second
_FillValue 9.96920996839e+36
valid_range [7600.0, 7600.0]

T_inst2ECR

_FillValue 9.96920996839e+36
valid_range [1.0, 1.0]

attitude_angles

units radians
_FillValue 9.96920996839e+36
pitch_element 1
roll_element 0
yaw_element 2

sun_ref

_FillValue 9.96920996839e+36
valid_range [1.0, 1.0]

num_impulse

valid_range [0, 25]

impulse_enc

units encoder pulses
_FillValue 9.96920996839e+36
valid_range [0.0, 16384.0]

impulse_time

units seconds
_FillValue 9.96920996839e+36
valid_range [0.0, 1.5]

Scan Type

thermal_correction

_FillValue 999.0
pitch_element 1
yaw_element 2
valid_range [180.0, 180.0]
roll_element 0
units degrees

attitude_quality

_FillValue 4096
interpolated_point 16384
overall 1
data_summary 2
yellow_limit_low 8
repaired_data_point 1024
quality_flag_problem 2048
short_gap_precedes 256
red_limit_low 4
long_gap_follows 64
yellow_limit_high 16
valid_range [0, 20139]
short_gap_follows 128
red_limit_high 32
long_gap_precedes 512

ephemeris_quality

_FillValue 4096
interpolated_point 16384
overall 1
data_summary 2
yellow_limit_low 8
repaired_data_point 1024
quality_flag_problem 2048
short_gap_precedes 256
red_limit_low 4
long_gap_follows 64
yellow_limit_high 16
valid_range [0, 18091]
short_gap_follows 128
red_limit_high 32
long_gap_precedes 512

Focal_length

units meters

band_position

units frame sample time

detector_space

units meters

detector_offsets

units meters
scandim 0
trackdim 1

T_offset

units frame sample time

num_samples

Code source:

#!/usr/bin/env python

from pyhdf.SD import SD, SDC

import pprint

file_name = 'MYD03.A2008190.1420.006.2012069052937.hdf'

file = SD(file_name, SDC.READ)

file_info = file.info()
print file_info  # number of sds and metadata

#----------------------------------------------------------------------------------------#
# print sds names

datasets_dic = file.datasets()

#for idx,sds in enumerate(datasets_dic.keys()):
#   print idx,sds

sds_dic = {}
for key, value in datasets_dic.iteritems():
    #print key, value, value[3]
    sds_dic[value[3]] = key

pprint.pprint( sds_dic )

#----------------------------------------------------------------------------------------#
# print SDS info

for i in range(file_info[0]):
    sds_obj = file.select(sds_dic[i])
    sds_info = sds_obj.info()
    print i, ' | ', sds_dic[i], \
             ' | ', sds_info[1] , \
             ' | ', sds_info[2] , \
             ' | ', sds_info[3] , \
             ' | ', sds_info[4]

#----------------------------------------------------------------------------------------#
# print SDS Attributes

for i in range(file_info[0]):
    sds_obj = file.select(sds_dic[i])
    sds_attributes = sds_obj.attributes()
    attributes = ''
    for key, value in sds_attributes.iteritems():
        attributes = attributes + key + ' ' + str(value) + '\n'
    print '#### ',  sds_dic[i]
    print
    print attributes.replace('-', '')
    print