members-wv2-princesscays-multi-pan-3

August 2014
Optical calibration of the PAN band

Now that the PAN solution using WV2 imagery is established, it still has to be illustrated and made into a practical operational process
The key element is the precise description of the variation of the operational diffuse attenuation coefficient for the WV2 Panchromatic band as water depth increases in this image
- this is a dimensionless relationship, which is based on optical path in the Blue and PAN bands
I now have designed in 4SM code an automated way to do that
- it requires a carefull preparation of the special effect mask
- using this WV2 image of Princess Cays, this in the end yielded a fairly reasonable result, as illustrated below
  - this plot mainly shows the curved "exponential decrease" of the linearized X_pan=log(Ls_pan-Lsw_pan) and X_blue=log(ls_blue-Lsw_blue), as water depth increases over the brightest bottom type
  - this plot also shows all isobath lines at 1 m interval, from 0 to 30 m
- this result still needs to be refined using a few more WV2 images over scenes where clean bright sands are more prevailing: it should be a very smooth curve, without any "road humps"

The Blue vs PAN observed optical model

As of september 2014,
the Blue vs PAN optical model
is now computed, using

WV02 response curve for the Panchro band,
Jerlov's data for the observed water type

The Blue vs PAN computed optical model
september 2014

August_2014
Results

WZ=all solutions	WZ=PAN solution only
image WZ Black....solution by the RedEdge band Red.......solution by the Red band Yellow...solution by the Yellow band Grey......solution by the PAN band	image WZ Grey.....solution by the PAN band
normalized TCC water column corrected using all solutions some areas exhibit weird results, likely caused by locally turbid waters: these areas yield underestimated depth results (and therefore darker water column corrected radiances)	normalized TCC water column corrected using only the PAN solution some areas exhibit weird results, likely caused by locally turbid waters this is a remarkable feature of the PAN soliution in 4SM: areas of questionable results are highlighted in dark purple in this normalized water column corrected display the PAN solution yields more homogeneous water column corrected results this improves the potential for bottom typing
image Z250 using all solutions	image Z250 using only the PAN solution

September_2014
Benefits of the PAN solution

Using the PAN solution is feasible and profitable: only uses Coastal, Blue, Green and PAN bands for retrieving depth
Based on the optical calibration for the MULTI bands, the optical calibration of the PAN band is derived, using the specific PAN response curve and Jerlov's data.
- this is a now generic and automatic process, hard coded in 4SM, which requires minimal manual control
- it may be used for any brand of imageries, only requires the specific PAN response curve
It yields very detailed results in the shallow range 0-5 m
It yields reasonably smooth results over the deeper range, thanks to the following:
- the level of system noise is quite low in this image, apart from vertical striping
- the computed optical model for the PAN band is a continuous smooth curve
No more of these annoying gaps that show up when using the MULTI solutions: removal of the water column effect is much improved

princesscayswvmp_201409_TCC_watercolumncorrected_normalized_PAN

BOA normalized TCC water column corrected (blue, green and red)
using only the PAN solution
The footprint of UKHO's MBES coverage is shown

princesscayswvmp_201409_TCC_watercolumncorrected_normalized_PAN_MBES

BOA normalized TCC water column corrected (blue, green and red)
over the MBES area

Retrieved depth

Note that

The result is remarkably smooth, thanks to the Blue vs PAN computed optical model for PANchro band: no more "road humps"
Water column correction is complete, as the retrieved depth image does not exhibit any trace of the bottom type variations

I wish I had access to UKHO'S MBES dataset, property of Princess Cruise Lines

Download from my GoogleDrive
~~princesscayswvmp_sept2014_results.zip (123 MB)~~

princesscayswvmp_sept2014_WCC.tif....................BOA watercolumn corrected radiances, from 1 to 9: channel4 is PAN
princesscayswvmp_sept2014_Zcm.tif......................Retrieved depth in centimeters
princesscayswvmp_sept2014_BottomType.tif..........Bottom type classification (step 2 final)

Command Line september 2014

nice -20 ./4SM -Process -Origin/DigitalGlobe
-DB/princesscayswvmp/50_3S_9S_0s/9_54/6001_5701/373.655_2732.957/1_1
-Mis/Bahamas/PrincessCays/WV02/MultiPan/PIX/bOA/UTM_+18_008/0.002_0.002/11_NOV_2011
-LS/005000/05000/05000/05000/05000/05000/05000/05000/05000
-dK000.000/0.000/0.015/0.000/0.000/0.150/0.040/0.030/0.000
@WL/0427.0/477.5/546.0/570.0/607.5/679.0/724.0/831.0/908.0
-LsM/200.2/192.0/187.5/215.4/202.5/210.4/210.3/180.7/175.6
-Lsw114.01/74.15/38.23/30.46/25.14/17.89/13.39/07.13/06.08_glint
dLsw/3.20/02.20/01.50/00.85/01.50/00.25/01.000/0000/00000
-Lw/0013.0/008.3/000.5/001.5/000.5/000.0/000.0/000.0/000.0
-SCL/00095/00107/00167/02500/00171/00116/00166/00164/00149
-Lm/0000.5/000.5/001.5/000.5/002.0/002.0/010.0/001.0/255.0
-Veg/109.4/072.8/070.2/072.0/037.5/027.8/093.1/139.3/131.6
-M/@000001/@0002/00003/00004/@0005/@0006/@0007/@0008/00009
-KK2_3_0.5500/Knir5.50/mask_3
-Z/MSL0.00d/cSL1.50/N_3/PAN_1_9/NDR_53_1000_zDTM/mask_3
-B/tclNe/cNorm5.00/cLM1
-deglint/VRbaD/q/F1/L9/FN/mDEGLINT1/GlintM28
-extract/v/rawBDH/FullBDH/NIRband8/NIRmax5/mapBPL/mSOIL21/mBPL2
-Calibrate/Va/BdSMpzg/plot_BPL/BC_6_5_4_2/bDh_12_23_24
-Model/v/runWV2/makeClouds/mask_3 -Tpsv/2/4/Text/MSE_3/Bmin100/panCode1.08
-Smooth/5/Csmooth4/D/Smart+    @-Eshp    -LL10    @-DNshp-DDshp
@-DD1686/1743    @-E1/6001/1/200 @-E1/6001/900/5701    @-DN4456/493
-Classify/Typer26_1.05/Classifier_3.5_0_scef_LastChType3/mode_5_16_25
@ProfileAB/LzrvWb/profile_white/dTM/chAB_6_7 _2_3_4
@-ProfileAB/lZrvWB/profile_purple/DTM/chAB_1_2_3 _5_6
@ProfileAB/lzRvWb/profile_purple/DTM/3_4/0_201/0_50_50/*