Recreating OCBG Quasi True Color VIIRS Images from Command Line
Posted: Thu Oct 19, 2017 11:30 am America/New_York
I'm currently unable to recreate the quasi true color image located on the Ocean Color Browse website. l1brsgen allows me to get pretty close, but the bowtie effect ruins the image. l1mapgen removes most of the bowtie effect, but the image isn't as clear as I would like and I can't seem to find the right projection. I'll provide a reproducible workflow and would really appreciate any help on this matter. I apologize if something similar had been posted before, I was unable to find anything.
*** Example ***
Goal: Recreate OCBG.png
Data found here:
https://oceancolor.gsfc.nasa.gov/cgi/browse.pl?sub=level1or2&per=8D&day=17447&prm=CHL&set=10&mon=17410&sen=v0&rad=0&frc=0&n=42.5390625&w=-70.6640625&id=1508264774.30949&cn=43.43&cs=41.69&cw=-71.57&ce=-69.72&dnm=D&file=V2017286183600.L2_SNPP_OC_CHLOR_A_BRS&ndx=4&wid=400&hgt=404
Command Line Syntax:
#Using l1brsgen and then R to convert l1brsgen_example.hdf to l1brsgen_example.png
l1brsgen ifile=V2017286183600.L1A_SNPP.nc geofile=V2017286183600.GEO-M_SNPP.nc ofile=l1brsgen_example.hdf subsamp=1 atmocor=1 resolution=1000 oformat=1 oformat_depth=24bit par=pars.txt
#Using l1mapgen
l1mapgen ifile=V2017286183600.L1A_SNPP.nc geofile=V2017286183600.GEO-M_SNPP.nc ofile=l1mapgen_example.png
***End Example Syntax***
Do I need to use l1brsgen in conjunction with l1mapgen? Does anyone have an example to follow for this?


*** Example ***
Goal: Recreate OCBG.png
Data found here:
https://oceancolor.gsfc.nasa.gov/cgi/browse.pl?sub=level1or2&per=8D&day=17447&prm=CHL&set=10&mon=17410&sen=v0&rad=0&frc=0&n=42.5390625&w=-70.6640625&id=1508264774.30949&cn=43.43&cs=41.69&cw=-71.57&ce=-69.72&dnm=D&file=V2017286183600.L2_SNPP_OC_CHLOR_A_BRS&ndx=4&wid=400&hgt=404
Command Line Syntax:
#Using l1brsgen and then R to convert l1brsgen_example.hdf to l1brsgen_example.png
l1brsgen ifile=V2017286183600.L1A_SNPP.nc geofile=V2017286183600.GEO-M_SNPP.nc ofile=l1brsgen_example.hdf subsamp=1 atmocor=1 resolution=1000 oformat=1 oformat_depth=24bit par=pars.txt
#Using l1mapgen
l1mapgen ifile=V2017286183600.L1A_SNPP.nc geofile=V2017286183600.GEO-M_SNPP.nc ofile=l1mapgen_example.png
***End Example Syntax***
Do I need to use l1brsgen in conjunction with l1mapgen? Does anyone have an example to follow for this?


