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'].variablesa=b['sensor_band_parameters'].variables
f=b['processing_control'].variablesfor key in d:print(key)chlor_a=d['chlor_a'][:].data#[:].data 大小 2030, 1354
chl_ocx=d['chl_ocx']
flags=d['l2_flags'][:].datayear=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'][:].datalon=e['longitude'][:].data   #记录了航线的经度
lat=e['latitude'][:].data    #记录了航线的纬度plt.plot(slat)
plt.plot(clat)
plt.plot(elat)
plt.show()

本文来自互联网用户投稿，文章观点仅代表作者本人，不代表本站立场，不承担相关法律责任。如若转载，请注明出处。 如若内容造成侵权/违法违规/事实不符，请点击【内容举报】进行投诉反馈！