In [17]:
import netCDF4
import numpy as np

fd = netCDF4.Dataset('geodetic_an.nc', 'r')
for name in fd.variables:
    print(name)


latitude_an
longitude_an
elevation_an
latitude_orphan_an
longitude_orphan_an
elevation_orphan_an

In [6]:
for d in fd.dimensions:
    print(d)


rows
columns
orphan_pixels

In [9]:
latitude = fd.variables['latitude_an']
print(latitude)


<type 'netCDF4.Variable'>
int32 latitude_an(rows, columns)
    standard_name: latitude
    long_name: Latitude of detector FOV centre on the earth's surface
    units: degrees_north
    scale_factor: 1e-06
    add_offset: 0.0
    valid_min: -91
    valid_max: 91
    _FillValue: -999
unlimited dimensions: 
current shape = (540, 3136)
filling on

In [11]:
print(latitude.ncattrs())


[u'standard_name', u'long_name', u'units', u'scale_factor', u'add_offset', u'valid_min', u'valid_max', u'_FillValue']

In [13]:
print(latitude.units)


degrees_north

In [18]:
latitude[:]


Out[18]:
masked_array(data =
 [[-- -- -- ..., 57.515093 57.513109 57.511126]
 [-- -- -- ..., 57.511027999999996 57.509045 57.507061]
 [-- -- -- ..., 57.506963999999996 57.504979999999996 57.502997]
 ..., 
 [-- -- -- ..., 55.317406999999996 55.315528 55.313649999999996]
 [-- -- -- ..., 55.313289999999995 55.311412 55.309532999999995]
 [-- -- -- ..., 55.309172999999994 55.307294999999996 55.305416]],
             mask =
 [[ True  True  True ..., False False False]
 [ True  True  True ..., False False False]
 [ True  True  True ..., False False False]
 ..., 
 [ True  True  True ..., False False False]
 [ True  True  True ..., False False False]
 [ True  True  True ..., False False False]],
       fill_value = -999)