I have a few questions about processing MSI and OLI, converting L1 to L2 with l2gen on a Linux machine. I have used SeaDAS/OCSSW from the GUI so far to get chlor-a from MODIS and OLI and have recently moved to using Linux - with command-line options. This is because I want to use MSI data.
- I have installed processors for MSI (2a and 2b) and OLI.
- I was able to convert OLI L1 to L2 with default parameters.
1) How can I find out what default processing options were used? I am want to know what option was used for aer_opt.
Also, is there a sample par file I can use? It should be a .txt file right?
2) On the GUI there is an option to choose 'Get Ancillary' -- is that the default?
3) in the GUI I use Rrs but also see rhos in the Product Selector:Radiance/Reflectance -- both are surface reflectances I believe - how are they
different? I tried look at the ATBD doc but am not clear!
4) Are the available products defined in some document (of course some are self-explanatory bu some are not to me)?
Lt, L_u, L_q,Lr, La, Zeu (_ee and _morel), zhl_morel, Zsd etc?
5) How do I get TOA reflectances?
Thanks so much for your help.
The default files are hierarchical, lowest common denominator is the $OCDATAROOT/common/msl12_defaults.par
They are ASCII text files
2) no, get ancillary is not run by default, but you should use it to get the appropriate ancillary data, otherwise climatologies are used
3) rhos is Rayliegh subtracted surface reflectance, it still includes aerosols.
4)The "L"s are TOA radiance components, Lt is the total, u and q are the polarized components of Rayliegh radiance, which is "r". "a" is aerosols. The "Z"s are various light penetration depths. "eu" is euphotic depth. "sd" is Secchi depth.
- Subject Matter Expert
- Posts: 225
- Joined: Thu Mar 05, 2009 10:25 am America/New_York
which will create the file A2008080195500.L1B_LAC.anc Now use the ancillary files in l2gen:
l2gen ifile=A2008080195500.L1B_LAC geofile=A2008080195500.GEO ofile=A2008080195500.L2 par=A2008080195500.L1B_LAC.anc par=<your_par_file>
Great! thanks very much for your reply. It worked fine.
Another question I have is about re-projection of l2gen output. I want to get the output file (which is .nc by default I believe) in WGS84/Geographic projection. In the GUI I just use the reproject option. How to do this in command line?
Thank you for the link to the 'gpt cookbook'! I have never used this before!
I am wondering what projection does the l2gen output file have. There are lat-lon values with each pixel.
For gpt.sh I am trying to follow Example 1 shown in the 'gpt cookbook'. I am not sure if there is anything I need to specify other than the resolution (pixelsize in .par file). Also, can I specify the output file to be a (.nc) file? is the (.dim) output file in NetCDF?
Thanks a lot for your help.
The next release of SeaDAS will include a more polished version of a python script that is in the current release (cleverly called l2mapgen.py) that used the l2bin/l3mapgen approach.
I had a recent exchange on this forum with another user with a similar question. It might help you as well :smile: