Polymer running error for MSI

Post Reply
Fluquorr
Posts: 3
Joined: Tue Apr 26, 2022 5:01 pm
company / institution: City University of New York
Location: New York

Polymer running error for MSI

Post by Fluquorr »

Hi all,
I'm trying to use Polymer for MSI data. But it shows
>>> run_atm_corr(Level1_MSI('/Users/Downloads/S2A_MSIL1C_20190612T222541_N0207_R015_T03VVL_20190613T001333.SAFE/GRANULE/L1C_T03VVL_A020747_20190612T222536/'), Level2_NETCDF(outdir='/Users/Downloads/test', ext='.nc'))
Initialize MSI projection EPSG:32603
/Users/opt/anaconda3/envs/polymer/lib/python3.8/site-packages/pyproj/crs/crs.py:130: FutureWarning: '+init=<authority>:<code>' syntax is deprecated. '<authority>:<code>' is the preferred initialization method. When making the change, be mindful of axis order changes: https://pyproj4.github.io/pyproj/stable ... -in-proj-6
in_crs_string = _prepare_from_proj_string(in_crs_string)
Trying to download https://oceandata.sci.gsfc.nasa.gov/cgi ... MI_24h.hdf ...
2022-04-28 15:51:39 URL:https://oceandata.sci.gsfc.nasa.gov/ob/ ... MI_24h.hdf [201417/201417] -> "ANCILLARY/METEO/2019/164/N201916400_O3_AURAOMI_24h.hdf.tmp" [1]
Trying to download https://oceandata.sci.gsfc.nasa.gov/cgi ... 6h.hdf.bz2 ...
2022-04-28 15:51:39 URL:https://oceandata.sci.gsfc.nasa.gov/ob/ ... 6h.hdf.bz2 [854076/854076] -> "ANCILLARY/METEO/2019/163/N201916318_MET_NCEPR2_6h.hdf.bz2.tmp" [1]
Trying to download https://oceandata.sci.gsfc.nasa.gov/cgi ... 6h.hdf.bz2 ...
2022-04-28 15:51:40 URL:https://oceandata.sci.gsfc.nasa.gov/ob/ ... 6h.hdf.bz2 [849948/849948] -> "ANCILLARY/METEO/2019/164/N201916400_MET_NCEPR2_6h.hdf.bz2.tmp" [1]
Starting processing at 2022-04-28 15:51:40.912055
Initializing output file "/Users/Downloads/test/L1C_T03VVL_A020747_20190612T222536.nc"
Processing block: size (198, 1830), offset (0, 0)
Traceback (most recent call last):
File "/Users/Downloads/Satellite_data/polymer-v4.14/polymer/main.py", line 169, in get_no2
self.no2_tropo_data
AttributeError: 'InitCorr' object has no attribute 'no2_tropo_data'

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File "/Users/Downloads/Satellite_data/polymer-v4.14/polymer/main.py", line 517, in run_atm_corr
for block in block_iter:
File "/Users/Downloads/Satellite_data/polymer-v4.14/polymer/main.py", line 405, in process_block
c.gas_correction(block)
File "/Users/Downloads/CUNY/Satellite_data/polymer-v4.14/polymer/main.py", line 231, in gas_correction
no2_frac, no2_tropo, no2_strat = self.get_no2(block)
File "/Users/Downloads/Satellite_data/polymer-v4.14/polymer/main.py", line 171, in get_no2
self.read_no2_data(mon)
File "/Users/Downloads/Satellite_data/polymer-v4.14/polymer/main.py", line 126, in read_no2_data
hdf1 = SD(self.params.no2_climatology)
File "/Users/opt/anaconda3/envs/polymer/lib/python3.8/site-packages/pyhdf/SD.py", line 1429, in __init__
_checkErr('SD', id, "cannot open %s" % path)
File "/Users/opt/anaconda3/envs/polymer/lib/python3.8/site-packages/pyhdf/error.py", line 23, in _checkErr
raise HDF4Error(err)
pyhdf.error.HDF4Error: SD (7): Error opening file


Shows same question when I use ./polymer_cli.py. What might be the problem? Thanks.

Tong
User avatar
fsteinmetz
Site Admin
Posts: 306
Joined: Fri Sep 07, 2018 1:34 pm
company / institution: Hygeos
Location: Lille, France
Contact:

Re: Polymer running error for MSI

Post by fsteinmetz »

Dear Tong,

Can you confirm that the file auxdata/common/no2_climatology.hdf is present ?
Have you downloaded all auxiliary data with "make auxdata_all" ?
Cheers,
François
Post Reply