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()