modis_atteph.py is downloading the wrong ATT/EPH pair

[user_gis]

If it's possible, let me know about the completion of work with the new API.
Thanks...

Regards,
Dmitry

[OB SeaDAS - xuanyang02]

Sure, we will notify you once the new API and script work as expected.

Meanwhile, you can still use the V2021.2 or earlier OCSSW

[user_gis]

Thanks,

I don't know your plans to update the API, but it would be nice if the new API could check MODIS data coverage with att/eph data before generating a GEO file using a given threshold, rather than just perform a data coverage check on their relevance to use.

What I mean:

Code: Select all

Aqua

Input file: A2022149111601.L1A_LAC
Sensor    : aqua
Start time: 2022149111601
End time  : 2022149112500

OBPG session started
Downloading A2022149111601.anc.server (    0.00 MBs)
[==================================================]
...Done
  Found: /opt/rki/ocssw/var/anc/2022/149/PM1ATTNR_NRT.A2022149.1110.061
  Found: /opt/rki/ocssw/var/anc/2022/149/PM1ATTNR_NRT.A2022149.1115.061
  Found: /opt/rki/ocssw/var/anc/2022/149/PM1ATTNR_NRT.A2022149.1120.061
  Found: /opt/rki/ocssw/var/anc/2022/149/PM1ATTNR_NRT.A2022149.1125.061
  Found: /opt/rki/ocssw/var/anc/2022/148/PM1EPHND.P2022148.1200.003
All required attitude and ephemeris files successfully determined and downloaded.
att1=/opt/rki/ocssw/var/anc/2022/149/PM1ATTNR_NRT.A2022149.1110.061
att2=/opt/rki/ocssw/var/anc/2022/149/PM1ATTNR_NRT.A2022149.1115.061
att3=/opt/rki/ocssw/var/anc/2022/149/PM1ATTNR_NRT.A2022149.1120.061
att4=/opt/rki/ocssw/var/anc/2022/149/PM1ATTNR_NRT.A2022149.1125.061
eph1=/opt/rki/ocssw/var/anc/2022/148/PM1EPHND.P2022148.1200.003

- All optimal ancillary data files were determined and downloaded. -

Code: Select all

echo $?

1

Then I create a GEO file and do modis_geocheck.

Code: Select all

modis_GEO --enable-dem A2022149111601.L1A_LAC

Code: Select all

modis_geocheck --verbose --threshold 100 A2022149111601.GEO
Percentage of pixels with missing geolocation:

0.00

Code: Select all

Validation test passed for geolocation file A2022149111601.GEO

It turns out that output status 1 has nothing to do with the actual ATT/EPH data coverage.

Thanks, again.

Regards,
Dmitry

[OB SeaDAS - xuanyang02]

True, modis_atteph.py does not have knowledge of failed pixel geolocation. It looks at the L1 file to determine the scene start and stop time which it then uses to select and download the correct ancillary data.