MODISA L1A to L1B questions

Use this Forum to find information on, or ask a question about, NASA Earth Science data.
Post Reply
alaroy
Posts: 78
Joined: Tue Nov 21, 2006 3:33 pm America/New_York
Answers: 0

MODISA L1A to L1B questions

by alaroy » Tue Jun 26, 2018 2:34 pm America/New_York

By way of background I'm looking to get atmospheric corrected and geolocated radiance for the 645nm, 859nm and 1240nm channels on MODISA.

I was able to convert an L1A file (A2018175065500.L1A_LAC to be precise) to L1B using the getanc.py, modis_GEO.py and modis_L1B.py scripts.

The results seem to be in HDF. Is there some option I can use to get netCDF L1B output instead?

Also three output files are created - ...L1B_HKM, ...L1B_LAC and ...L1B_QKM. What are the differences between them?

-Andrew L.

Tags:

OB.DAAC - SeanBailey
User Services
User Services
Posts: 1464
Joined: Wed Sep 18, 2019 6:15 pm America/New_York
Answers: 1
Been thanked: 5 times

MODISA L1A to L1B questions

by OB.DAAC - SeanBailey » Tue Jun 26, 2018 3:32 pm America/New_York

Andrew,

The MODIS L1B code outputs HDF as that is what it outputs - no, there is no option to have it output netCDF...although you *can* run the resulting L1B HDF file through l1bgen_generic and generate a generic L1B file in netCDF format if you desire....but it's not necessary if you're going to use l2gen to process for atmospheric correction.

The three files are because there are three resolutions 1km, 500m and 250m - not all bands are at all resolutions, e.g. only 645 and 859 are at 250m.

Sean

alaroy
Posts: 78
Joined: Tue Nov 21, 2006 3:33 pm America/New_York
Answers: 0

MODISA L1A to L1B questions

by alaroy » Tue Jun 26, 2018 4:21 pm America/New_York

Thank you Sean, that helps! :)

I'm getting an error when trying to use l1bgen_generic:

>l1bgen_generic ifile=A2018175065500.L1B_LAC ofile=A2018175065500.L1B_LAC.nc


Input file A2018175065500.L1B_LAC is MODIS Aqua Level-1B HDF-EOS product.
...
Reading Thuillier_F0.dat
-W- /home/seadas/ocssw/src/l2gen/l1_hmodis_hdf.c:399: Cannot open file:
Error reading ; please specify a valid geolocation file.

There is a A2018175065500.GEO in the same directory as A2018175065500.L1B_LAC. Is there a command line parameter I'm missing?

-Andrew L.

OB.DAAC - SeanBailey
User Services
User Services
Posts: 1464
Joined: Wed Sep 18, 2019 6:15 pm America/New_York
Answers: 1
Been thanked: 5 times

MODISA L1A to L1B questions

by OB.DAAC - SeanBailey » Wed Jun 27, 2018 8:39 am America/New_York

geofile=A2018175065500.GEO

alaroy
Posts: 78
Joined: Tue Nov 21, 2006 3:33 pm America/New_York
Answers: 0

MODISA L1A to L1B questions

by alaroy » Thu Jun 28, 2018 8:55 am America/New_York

That worked great, thank you Sean! :)

This is no longer an L1B question but since I already have this thread going I'll ask it anyway...

I want to select Rrs_645, Rrs_859, Rrs_1240 an l2_flags as my custom L2 products. I see the l2prod flag in the l2gen command line help but don't know what to set it to in order to get those. Can you please give me a pointer? Thank you.

OB SeaDAS - knowles
Subject Matter Expert
Subject Matter Expert
Posts: 269
Joined: Mon Apr 07, 2008 4:40 pm America/New_York
Answers: 0

MODISA L1A to L1B questions

by OB SeaDAS - knowles » Thu Jun 28, 2018 9:31 am America/New_York

Andrew,

Here's a sample l2gen command:
l2gen ifile=A2018175065500.L1B_LAC geofile=A2018175065500.GEO ofile=A2018175065500.L2_LAC_RRS l2prod="Rrs_645,Rrs_859,Rrs_1240" resolution=500

The product "l2_flags" will automatically get created so you don't need to specify that.

Note the native instrument resolution: 645 and 859 are 250m; 1240 is 500m.  Running l2gen at 500 resolution will preserve 1240 and downsample 645 and 859.  Likewise, running l2gen at 250 resolution will upsample (interpolate) 1240 and preserve 645 and 859.  No specification for resolution will downsample everything to 1km resolution (default).

Danny

Post Reply