I am trying to regrid the MODIS_L2 10km Aerosol data from Swath to regular grid.
As suggested on hdfeos forum (http://hdfeos.org/zoo/LAADS_MOD_py.php), I am using pyresample (a geo-resample package in python) to do the same. Since the daily MODIS_L2 swath data comes in huge numbers (particularly in case of a decade of analysis), so post-processing during data order is another option to subset the SDS and ask for the Mosaic from multiple swath to single region of interest swath file.
But in doing so in post-processing the geolocation field (Geo-2D latitude & longitude as in Level 2 hdf files) get removed which adds an another level of challenge to determine the latitude and longitude of the hdf files. While the same post-processed .hdf files in panoply shows virtual lat lon in Data_Fields as 1D.
Is there any robust option to retain the same geolocation field as regular Level2 hdf files and ask for doing post-processing? or any other way to sub-set and mosaic the multiple swath files to single?
Thanks in anticipation!
- User Services
- Posts: 266
- Joined: Mon Sep 30, 2019 8:33 am America/New_York
- Has thanked: 1 time
LAADS User Services
To receive news from LAADS DAAC direct to your inbox, email email@example.com with “subscribe” in the subject line.
I am following up on this query, as I am undertaking a similar task to the OP. I want to analyse several years of L2 10km aerosol data (from MYDATML2/MODATML2 files) for a specified (large) region, and am interested by the LAADS mosaic and reproject post-processing options. However, there is a limit of 30 files per order, and this does not cover a single day even when restricted to files that only cover my region of interest (via the geographic search utility when ordering).
I am therefore wondering if there is a workaround to facilitate post-processing on a larger number of files; or whether the only LAADS solution is to submit a huge number of orders (which would be on the order of thousands for several years of data).
Note, I am also exploring the option of bulk-downloading swath files and processing them myself following python examples on the hdfeos forum (http://hdfeos.org/zoo/LAADS_MOD_py.php), but am encountering some issues with that too (specifically with mosaicing overwriting earlier swath data with later - an issue that I have raised on the forum there).
Thanks in advance for any help or guidance.
It's quite easy using python.
I have composed a python script for multiple days with global coverage.
My script is on github. Hope it would help you. The link to my github repository is https://github.com/kmmrao/Mosaicing-MODIS-Level-2-Global-AOD-and-Visualizing-/blob/main/MODIS%20L2%2010KM%20AOD%20Regridding%20Single%20File%20and%20plotting.ipynb