Optical calibration, bathymetry, water column correction and bottom typing of shallow marine areas, using passive remote sensing imageries
water column correction since 1994
4SM : an operational suite of tools 
written in C under Ubuntu Linux (60,000 lines)
  • no need for field data
  • no need for atmospheric correction
  • works with DNs
  • "same day" service
  • from raw spectral data to final deliverable products 
  • through one single commandline


 
1 - NO NEED for field data, nor for atmospheric correction
2 - this is demonstrated in this website, using a variety of hyper/multi spectral data
 
Requirements are
1 - homogeneous water body and atmosphere
2 - some coverage of optically deep water
3 - some coverage of dry land
 
Problems are
1 - the precision on estimated depth is found wanting, because the noise-equivalent change in radiance  of accessible data is too high for shallow water column correction work 
2 - radiance data should be preprocessed by the provider at level 1 in order to improve S/N ratio
3 - exponential decay: the deeper/darker the bottom, the poorer the performances
 
So
I keep digging
until suitable data
become available
 
 
The simplified RTE http://www.watercolumncorrection.com/sld016.php
JERLOV wrote http://www.watercolumncorrection.com/sld036.php
KIRK       wrote http://www.watercolumncorrection.com/sld037.php
A family of curves http://www.watercolumncorrection.com/4sm-illustrations-calibration-process.php#porquerolles
The SL and BPL http://www.watercolumncorrection.com/4sm-illustrations-calibration-process.php
Calibration diagram http://www.watercolumncorrection.com/4sm-illustrations-calibration-process.php#yes
Deglinting http://www.watercolumncorrection.com/4sm-illustrations-the-deglinting-process.php#lsi
Modeling http://www.watercolumncorrection.com/waimanalo-pollution.php#profiles
Seatruthing http://www.watercolumncorrection.com/waimanalo-pollution.php#regression_west
RMSE http://www.watercolumncorrection.com/waimanalo-pollution.php#Lmyellow=4
Bottom typing http://www.watercolumncorrection.com/waimanalo-bottom-typing.php#July_25th
The 4SM command line  http://www.watercolumncorrection.com/waimanalo-pollution.php#commandline













 

slide 1

  4SM: a ratio method 
for water column correction of shallow spectral images

 
Following work by Lyzenga, Philpot, Maritorena and Jerlov
 Simplified TOA RTE : Ls – Lsw =          (LsB – Lsw) / exp (2K * Z)
 Inverted    BOA RTE : LB           = Lw + (Ls    – Lw  ) *exp (2K * Z)
where Lsw=La+Lw, LsB=LB+La, Ls=L+La
 
No need for atmospheric correction: works with Dns
From raw data to final deliverables “on the fly” in one single commandline
Aims at “Same day” service 





 
slide 2
Deglinting in 4SM
Hedley's method  
  • Get the linear regressions of Lsi vs LsNIR
    • for all pairs of bands i and NIR over an optically deep area
  • The amount of glint at the current pixel is LsNIR-LswNIR
  • Use the regressions to estimate the amount of glint in other bands,
 
 
For WV2 images
 






 
slide 3

WV2 at Waimanalo Bay, Oahu, Hawaii

LIDAR seatruth
Zcomputed-ZLidar
Red tones for over-estimated depths
Blue tones for under-estimated depths
Dark_blue tones delineate artifacts caused by ”pollution”                   
  
Profile_black                                                                                          Enhanced TCC
of water column corrected image

                     Deep purple tones signal "pollution"






 
slide 4

WV2: Profile_black at Waimanalo Bay, Oahu, Hawaii
         





 
slide 5

 
 
 






 
slide 6
WV2 at waimanalo Bay 
Reducing the error on estimated depth

LIDAR DTM is in plain meters: I need a DTM in centimeters

Away from suspected pollution
from row_1 to row_1000

from 0 to 10 m along profile_black

96.2% of all estimated depths  
are within +-2 m of DTM

 
RMSE=0.84-0.50=0.34 m

average ZComputed-ZLidar = 0.61-0.50=0.11 m
 

This image is a "worst case"
fully glinted image
Sensor is looking 37° sideways
right into the sun


Thanks you, Ron!
 
 September 1rst 2012