4SM Rangiroa OLIP

rangiroa_SPOT
rangiroa_ETM
rangiroa_OLI
rangiroa_CASI

Tikehau bathymetry
Rangiroa bathymetry

3454*1590, 30 m pixel size, UTM zone 6
scene LC80530702013233LGN00, downloaded from USGS
August 21rst 2013
Work done september 2013

This image is a mess.
Waters inside the lagon become distinctly less clear:in 2013, 4SM can't cope
home
The data

Deglinting
Optical calibration
Retrieved depth
Command Line

The data
Seven bands are used: Purple, Blue, Green, PAN, Red, NIR and SWIR1

FCC	TCC
GREEN band	Panchromatic band, sum-aggregated to 30 m,

Deglinting

Glint regression based on SWIR1 band	Glint regression based on NIR band
Deglinted TCC	Note: PAN not usable that the PAN band exhibits very poor correlation this means that the PAN band may not be deglinted to satisfaction, and therefore shall not be usable for shallow water work unlike at Caicos

Optical calibration:
Water type OIB of Jerlov
Spectral reflectance of the brightest bottoms (LsM):
0.207 0.252 0.294 0.314 0.329 0.463 0.343
on the scale 0-1, as calibrated using the _MTL.txt textfile
This is distinctly darker than at Caicos Bank

Calibration diagram
for bands Blue, Green, Red and NIR

Calibration diagram
for bands Purple, Green, Red and NIR

Spectral operational 2K

This calibration yields the following operational 2K:

2K/0.113/0.106/0.242/0.393/0.807/4.490/15.700
2K_blue =0.106 m^-1
2K_green=0.242 m^-1
2K_red =0.807 m^-1
2K_nir =0.450 m^-1 a customary value for NIR

Quite possibly, seatruth should reveal that these values are realistic, i.e. that retrieved depths shall be correctly estimated.

Final_Z = CoefZ*Z - Tide_to_chart_datum

Operational wavelengths

These diagrams demonstrate that a very good fit of the BPL mode has been achieved among all visible bands
This was obtained by setting
- WLpurple at 450 nm
- WLgreen at 577 nm
- WLred at 655 nm

This is much more straightforward and acceptable than what was adopted at Caicos

Water column correction: retrieved depth along Profile Red
Performance of the Green band reaches ~30 m over bright bottoms
in this OIB +0.6 water type of Jerlov

Retrieved depth shows Profile Red Retrieved depth shows Profile Green	Retrieved depth along Profile Red also shows water column corrected Blue band
Plot of Z3 vs Z5 along Profile Green Z5 uses Purple, Blue, Green vs Red Z3 uses Purple, Blue vs Green Deglinting and Smart-Smoothing applied This plot is used to fine-tune several critical model parameters, like WL: wavelengths Lsw: deep water radiances LsM: brightest shallow bottom Lm: threshold value until a satisfactory fit is obtained	Plot of Z3 and Z5 along Profile Green Z5 uses Purple, Blue, Green vs Red Z3 uses Purple, Blue vs Green This plot also shows the watercolumn-corrected reflectances (scaled 0-250) in the Purple, Blue and Green bands Deglinting and Smart-Smoothing applied

Command Line

nice -20 4SM -Process    -Origin/USGS
-DB/rangiroa_20130821/40_3S_7S_0s/7_44/3454_1590/360.825_-1647.015/1_1
-Mis/Tuamotu/Rangiroa/Landsat/OLIP/TIF/UTM_6_008/0.030_-0.030/21_AUG_2013
@WLM/0450.0/510.0/590.0/680.0/670.0/880.0/1650
@WLm/0430.0/450.0/530.0/500.0/640.0/850.0/1570
-WL/00449.8/480.0/550.0/590.0/655.0/865.0/1610
-MakePIX
-import/data*LC80530702013233LGN00/dbnc_0_0s_1S_0s/R3454_L1590/Origin_360.810_-1647.000_0_0/chIn1/chOut44
-SCL/02243/02397/03003/02892/02939/03207/02742    @-Prepare/Scale
-LS0030000/30000/30000/30000/30000/30000/30000
-LsM/200.0/200.0/200.0/200.0/200.0/200.0/200.0_new
-M/@000001/@0002/00003/@0004/00005/00006/00007
-Lm/0001.0/001.0/001.0/001.0/001.0/003.0/255.0
-KK2_3_0.690/Knir15.70/Bm400_BG520_GR590_RN700/mask_3
-Z/MSL0.00d/cSL1.5/N_2/nDR_53_zDTM/mask_3            -LL10
-extract/v/NIRband7/NIRmax5.0/mBPL2/mSOIL21/mask_3
-AutoCAL/print_2_3/mapBPL/MakeNewMask/getGlint/GlintM58/offset5/slope1.10/Clouds @-DNshp

nice -20 4SM -Process -Origin/USGS
-DB/rangiroa_20130821/40_3S_7S_0s/7_44/3454_1590/360.795000_-1646.955/1_1
-Mis/Tuamotu/Rangiroa/Landsat/OLIP/TIF/UTM_6_008/0.030_-0.030/21_AUG_2013
-import/data*LC80530702013233LGN00
-LS/030000/30000/30000/30000/30000/30000/30000
@WLM/0450.0/510.0/590.0/680.0/670.0/880.0/1650
@WLm/0430.0/450.0/530.0/500.0/640.0/850.0/1570
-WL/00450.0/480.0/577.0/590.0/655.0/865.0/1610
-M/@000001/@0002/00003/@0004/00005/@0006/00007
-SCL/02243/02397/03003/02892/02939/03207/02742_new
-LsM/200.0/200.0/190.0/200.0/200.0/200.0/200.0_new
-Lsw/110.0/090.6/057.4/056.8/051.5/043.1/047.5_new
-dLsw000.0/000.0/000.0/00-.2/00-.2/000.0/000.0_new
-Lw/0008.0/004.5/000.8/000.0/000.0/000.0/000.0_new
-Lm/0001.0/001.0/000.1/001.0/001.0/001.0/001.0_new
-Veg/114.2/088.6/070.0/147.3/075.6/186.8/165.5_new
-KK2_3_0.44/Knir15.70/Bm400_BG520_GR590_RN700/mask_4    OIB+0.6
-Z/MSL0.00d/cSL1.5/N_2/nDR_53_zDTM/no_mask
-B/tclNe/cNorm2.50/LBref200_100/Bot_0_250_254_253/cLM1 .5
-deglint/vRbaD/GlintM68
-extract/v/rawBDH/FullBDH/NIRband7/NIRmax5/mapBPL/mSOIL21/mBPL3/mask_3 Eshp
-calibrate/v/BdSNpzg/BC_7_5_3_2/BDh_12_13_14
-Model/v/runWV2/mask_3        @-Tp/2/4/text -LL1
-Smooth/5/D/Smart+
@ProfileAB/LVW/profile_red/chAB_1_2_3_4_5_6_7
@ProfileAB/ZrvwB/profile_red/dTM/chAB_2 1_2_3_5
-ProfileAB/RvB/profile_green/dTM/3_5/0_201/0_50_50/*