MODIS L2级产品处理

 

MODIS L2数据说明

https://oceancolor.gsfc.nasa.gov/docs/format/l2nc/

详细数据例子

https://oceancolor.gsfc.nasa.gov/docs/format/l2oc_modis/

import netCDF4
import matplotlib.pyplot as plt
file='H:\\MODIS L2 data\\A2003001061000.L2_LAC_OC.nc'
# file='D:\\A2011009003500.L2_LAC_IOP.nc'
data=netCDF4.Dataset(file)
b=data.groups
# for key in b:
#     print(key)


# numbers of band? 24或者26应该怎么处理
c=b['scan_line_attributes'].variables
d=b['geophysical_data'].variables
e=b['navigation_data'].variables

a=b['sensor_band_parameters'].variables
f=b['processing_control'].variables

for key in d:
    print(key)

chlor_a=d['chlor_a'][:].data#[:].data 大小 2030, 1354
chl_ocx=d['chl_ocx']
flags=d['l2_flags'][:].data

year=c['year'][:].data
day=c['day'][:].data
slon=c['slon'][:].data  
clon=c['clon'][:].data
elon=c['elon'][:].data
slat=c['slat'][:].data
clat=c['clat'][:].data
elat=c['elat'][:].data
csol_z=c['csol_z'][:].data

lon=e['longitude'][:].data   #记录了航线的经度
lat=e['latitude'][:].data    #记录了航线的纬度

plt.plot(slat)
plt.plot(clat)
plt.plot(elat)
plt.show()