laparguera-20130926-3

Modeling results

LIDAR depth, rotated LIDAR depth, rotated see depth legend	Retrieved depth, PAN solution image Z see depth legend
Average bottom brightness image B VERY DARK: one can see that the green-purple-black hues signal severe water turbidity and therefore severe depth underestimation off Mayaguez and at La Parguera Marine Reserve BRIGHT: pink hues are rendered where bottom is very bright depth is over-estimated, that is where waters becomes clearer than assumed	GREEN solution vs PAN solution GREEN solution vs PAN solution for Profile Purple Waters are far from homogeneous !!!!!!!!!!!
Blue-pink: isolated very bright pixels/patches, possibly of clean coral sands Brown: topograhic highs, probably coral/seaweed covered hard bottom	Retrieved depth is too deep over the ridge
TCC BOA, water column corrected reddish areas signal waters contaminated in the yellow-red range by slightly turbid coastal waters

LIDAR/4SM comparison

GOOD: results are obviously good in the surrounding of Pole Ojea.
BAD: retrieved depths are badly underestimated at La Parguera Marine Reserve and off the city of Mayaguez.
Access to the LIDAR image has helped me a lot to understand/decipher this situation.

CoefK

HICO at La Parguera Marine Reserve, Puerto Rico:
K_blue/K_green=0.60 ==> Water type OIB+0.4 ==> CoefK=0.5

Command Line may 2016

nice -20 ./4SM   -Process/19_DEC_2014 -Origin/USGS           @-Prepare/Scale
-DB/prhico_20130926/75_3S_28S_0s/28_79/500_900/0.000_900.000/0_0
-Mis/PuertoRico/Mayaguez/HICO/HICO/OSU/bOA/UTM_+00_008/0.1000_0.1000/26_SEP_2013
-LS/10000/10000/10000/10000/10000/10000/10000/10000/10000/10000/10000/10000/10000/10000/10000/10000/10000/10000/10000/10000/10000/10000/10000/10000/10000/10000/10000/10000
@WLM/412.7/424.1/435.6/447.0/458.5/469.9/481.4/492.9/504.3/515.8/527.2/538.7/550.1/561.6/573.0/584.5/596.0/607.4/618.9/630.3/641.8/653.2/664.7/676.1/687.6/699.1/710.5/848.0
@WLm/401.2/412.7/424.1/435.6/447.1/458.5/470.0/481.4/492.9/504.3/515.8/527.2/538.7/550.2/561.6/573.1/584.5/596.0/607.4/618.9/630.4/641.8/653.3/664.7/676.2/687.6/699.1/836.6
@WL/0406.9/418.4/429.9/441.3/452.8/464.2/475.7/487.1/498.6/510.0/521.5/533.0/544.4/555.9/567.3/578.8/590.2/601.7/613.2/624.6/636.1/647.5/659.0/670.4/681.9/693.3/704.8/842.3
-SCL/00088/00108/00111/00132/00155/00151/00163/00160/00158/00147/00157/00172/00180/00194/00190/00187/00180/00190/00192/00188/00188/00176/00177/00186/00174/00163/00187/00194_new
-LsM/701.0/668.3/613.4/614.9/607.5/618.8/613.3/623.2/628.8/630.8/627.6/622.5/623.9/614.3/613.9/611.6/641.2/609.4/609.3/610.0/610.3/612.2/611.8/611.6/612.1/612.5/612.7/602.8_cLM=3.000
-Lsw/159.5/136.0/120.8/106.8/091.3/091.9/080.7/075.3/066.3/056.2/042.3/030.2/023.7/018.8/017.9/017.1/017.0/015.9/015.4/015.7/015.3/015.6/015.4/014.7/015.7/016.1/019.8/007.3_Lsw
-Lsw/159.5/136.0/120.8/106.8/091.3/091.9/080.7/075.3/066.3/056.2/042.3/030.2/023.7/018.8/017.9/017.1/017.0/015.9/015.4/015.7/015.3/015.6/015.4/014.7/015.7/016.1/019.8/007.3_Lsw
-La/0060.1/045.6/042.8/039.5/026.0/034.1/024.2/030.0/031.0/030.6/026.8/021.0/017.1/013.4/012.7/012.2/012.9/012.7/012.8/013.6/013.7/015.2/015.0/014.7/15.65/016.1/011.9/007.3_La
-dLs/000.0/000.0/000.0/000.0/000.0/000.0/000.0/000.0/000.0/000.0/000.0/000.0/000.0/000.0/000.0/000.0/000.0/000.0/000.0/000.0/000.0/000.0/000.0/000.0/000.0/000.0/000.0/000.0_dLsw
-dLs/-37.3/-30.2/-27.2/-23.8/-18.3/-17.0/-13.3/-09.5/-04.9/-00.8/000.9/001.6/001.7/001.8/002.0/001.4/000.6/-00.3/-00.5/-00.4/-00.4/-00.3/-00.4/000.1/000.3/000.0/000.5/000.0_dLsw
-Veg/097.0/090.1/085.8/081.2/074.3/078.7/024.2/074.6/075.7/081.0/081.0/078.9/050.6/073.3/071.6/067.9/064.1/059.0/056.4/054.9/053.1/052.5/049.1/045.8/048.0/058.1/38.0/152.9_Veg
-M/@000001/@0002/@0003/@0004/@0005/@0006/@0007/@0008/@0009/@0010/@0011/@0012/00013/@0014/@0015/@0016/@0017/@0018/@0019/@0020/@0021/@0022/@0023/@0024/@0025/@0026/00027/00028-M
-Lm/0255.0/255.0/255.0/001.0/001.0/001.0/001.0/001.0/001.0/001.0/001.0/001.0/000.1/001.0/001.0/001.0/001.0/001.0/001.0/001.0/001.0/001.0/001.0/001.0/001.0/255.0/000.1/255.0_Lm
-KK7_13_0.70/Knir4.500/CoefK_0/Bm400_BG500_GR620_RN700/mask_3
-Z/MSL0.00d/cSL1.18/N_9/nDR_78_1000_zDTM/mask_3
-B/TclNe/cNorm1.00/Bmin0/cLM 0.28 3 2.5 0.7
-deglint/vrbaD/Q/F1/L26/FN/mDEGLINT0/GlintM50 255
-extract/v/rawBDH/FullBDH/NIRband28/NIRmax5.0/mSOIL21/mBPL2 /Mask_6
-calibrate/vaI/BdSMpzg/plot_BPL/BC_27_22_13_7/BDh_12_13_14
-Model/v/runWV2/DTM/makeClouds /Mask_6   -LL10
-Smooth/1/D/Csmooth13/smart+
-Tpsv/7/27/Text/Bmin150/RLD4/MSE_2/PanCode
@-E1/500/1/900 @-E1/500/1/900
@-DN300/100/3/comment @-DDshp @-DNshp @ -Eshp
@Classify/typerp_1.05/Classifier_40_2_scef_LastChType15/mode_5_16_25
@ProfileAB/LzrvWb/profile_blue/dTM/chAB_2_4_7_13_26 _27 7_13_26_27
@ProfileAB/lZrvWB/profile_white/DTM/chAB_3_7_13 _16_20 22_26
@ProfileAB/lZrvWB/profile_red/dTM/13_27/0_201/0_50_5/*
@ProfileAB/lzRvWC/profile_purple/dTM/13_27/0_201/0_30_2/*
-RegressZZ/v/p/DTM_2/CXY_1.0_0_0/CZR_1.000_0.000/ZR/WZ_13_27/B_0_201/ZC_0_30_50/ALL
@Recode/234_InChannel_29_into_248_OutChannel_29/248
@Recode/240_InChannel_34_into_222_OutChannel_34 /mask_5 /land
@Recode/LT0_InChannel_76_into_0_OutChannel_76 /mask_3 /land
@Recode/2_InChannel_29_into_1_OutChannel_29/mask_1 /mSOIL
@Recode/LT28_InChannel_29_into_29_OutChannel_29/dLsw
@Recode/GT2_InChannel_29_into_5_OutChannel_29/mask_5
@Recode/GT21_InChannel_26_into_234_OutChannel_29/shallow
@Recode/LT20_InChannel_29_into_20_OutChannel_29/dLsw
@Recode/244_InChannel_29_into_240_OutChannel_29 /deep

Seatruth on depths
retrieved by Hernandez et al. from a WV2 image

This is area located in the surrounding of Pole Ojea.
Hernandez et al. had to go through atmospheric correction,
and used the LIDAR data to calibrate their very special "second order polynomial equation" method.
We want to see such seatruth regression over the yellow study site.

4SM work clearly obtained good results in this area,
using just the image and Jerlov's data.
I wish I could work on this WV2 image!