如何在python中用format.nc打开/获取Modis数据(二级)的变量?

2024-05-23 19:06:20 发布

您现在位置:Python中文网/ 问答频道 /正文

我在python中使用.nc格式打开/获取modis数据(级别2)的变量时遇到问题。通常,{cd4>{cd4>

from netCDF4 import dataset

f = r"F:\WORK\L3_modis\input\A2016006.L3m_DAY_SST_sst_4km.nc"    
fh = Dataset(f, mode='r')    
lons = fh.variables['lon'][:]
lats = fh.variables['lat'][:]  
sst = fh.variables['sst'][:]

但是这个脚本对二级modis数据不起作用。脚本如下:

^{pr2}$

问题是它找不到二级modis数据的变量(chlo_asst,等等…)。在

据我所知,当我尝试使用hdf格式的modislevel2时,它是有效的。氯霉素——一种称为“氯-a”的变量。脚本如下:

import pyhdf.SD import SD
f=SD('D:\\ORS\cilacap_pacitan_2013079\A2013079064500.L2_LAC_OC.x.hdf')
print f.datasets().keys()
['l2_flags', 'elat', 'clon', 'vcal_offset', 'msec', 'vcal_gain', 'chlor_a', 'year', 'elon', 'wavelength', 'tilt', 'cntl_pt_rows', 'latitude', 'slon', 'slat', 'Tau_r', 'day', 'F0', 'k_oz', 'longitude', 'cntl_pt_cols', 'clat', 'csol_z']
data=array(F.select('chlor_a').get()) #trying to get chlor_a variable
print data

Tags: 数据import脚本格式sdvariablesmodisprint