Optical calibration, bathymetry, water column correction and bottom typing of shallow marine areas, using passive remote sensing imageries
WorldView 2 image Cockburn Sound, Western Australia
June 2nd 2012, 11210*14500, 2 m ground resolution, 22.42 km * 29 km
"(c) DigitalGlobe, Inc. All Rights Reserved"
"Includes copyrighted material of DigitalGlobe, Inc., All Rights Reserved"
image courtesy of Gregory Miecznik, DigitalGlobe, under a Data Exchange agreement

 
the "PAN solution"
uses the PANCHRO band,
4x4 sum-aggregated and coregistered
bands 1, 2, 3, Panchro, 4, 5, 6, 7 and 8       
become      
bands 1, 2, 3, 4, 5, 6, 7, 8 and 9
home
 
 
 


Modeling: exclusively by the "PAN solution"
that is Purple, Blue and Green against PAN (other bands not used at all)

Retrieved depth in centimeters
  • I have adjusted Lswpan manually so as to obtain ~20m+, as shown in the bathymetric map of Cockburn Sound for the Central Basin.
  • This was the only use of existing field data.

 

Water column corrected reflectances
RGB color composite
for bands Green, Blue and Purple

image B
This scene exhibits a lot of bright bottoms,
unlike the marmion scene

image NDR
NormalizedDifferenceRatio=(LB1-LB3)/(LB1+LB3)
  • negative values are mapped in shades of green
  • positive values are mapped in shades of purple
 
Pure fantasy:
Perth River and ground water outflow
  • I think polluted Perth River and aquifer waters outflow causes foul waters to linger in the very dark area north of Success Bank.
  • Shallow retrieved depth and extremely dark bottoms  results there are pure fantasy, as the  radiances are at or below adopted deep water radiance
  • This area ought to be masked/thresholded out of the results alltogether, as depths there are reported to be similar to Cockburn Sound.





Profiles

Profile RED from West to East
  • Profile RED B north of Carnac Island, right through Owen Anchorage at depths 10-15 m
    • the sudden drop at km 4 is related to the sharp contact front between clean oceanic waters and turbid coastal waters: most probably an artifact: this image really is shit, but these results still might useful(?).
  • Profile RED A crosses Cockburn Sound.
    • It reaches in excess of 20 m of depth where expected at ~km 10

Profile YELLOW from West to East
  • Profile YELLOW B shows depth ~15 m at the narrow dredged boat channel (waypoint B1)
  • Profile YELLOW B shows depth ~15 m at the narrow dredged boat channel (waypoint A1)
   
   





Command Line
nice -19  4SM -Process    -Origin/DigitalGlobe                               
@DB/cockburnwvmp/50_3s_9S_0s/9_54/11210_14500/362.567_6456.625/0_0/1_1                   
@DB/cockburnwvmp_Zcm/00_1s_0S_0s/1_1/11210_14500/362.567_6456.625/0_0/1_1               
-DB/cockburnwvmp_LB/04_0s_0S_0s/1_1/11210_14500/362.567_6456.625/0_0/1_1               
-Mis/WesternAustralia/CockburnMP/WV02/MultiPan/PIX/UTM_50_008/0.002_-0.002/02_JUN_2012           
-MakePIX                                               
@Import/data*perth_cockburn/dbnc_0_0s_9S_0s/R11210_L14500/Origin_362.566_6456.626_500_10000/chIn1,9/chOut54,62   
@WLm/401.0/447.0/511.0/PAN.0/588.0/629.0/704.0/772.0/772.0                       
@-WL/427.0/477.5/546.0/570.0/595.8/641.0/724.0/831.0/908.0                       
@WLM/453.0/508.0/581.0/PAN.0/627.0/689.0/744.0/890.0/954.0                       
@cWL/0.500/0.500/0.570/0.500/0.500/0.700/0.500/0.500/0.500                       
-LS/005000/05000/05000/05000/05000/05000/05000/05000/05000                       
-M/@000001/@0002/@0003/00004/00005/00006/00007/00008/00009                       
-SCL/00055/00074/00114/00080/00104/00096/00104/00118/00078     @-Prepare/Scale           
-LsM/196.3/220.8/197.0/305.4/180.1/193.8/191.7/207.9/197.3                       
-Lsw/117.8/081.8/043.1/046.3/024.0/018.7/012.2/006.8/006.0                       
-dLsw000.0/000.0/000.0/000.2/000.0/000.0/000.0/000.0/000.0                       
-Lw/0008.4/004.2/002.5/004.5/000.5/000.0/000.0/000.0/000.0                       
@extract/v/NIRband4/NIRmax5.0/mBPL2/mSOIL21                                
@deglint/VRbaD/DegTol1.00/F1/L8/N1/FN/mDEGLINT1/GlintM255                       
-AutoCAL/print_2_3/mapBPL/MakeNewMask/offset5/slope1.7/GetGlint/GlintM16.0/clouds    /offset5/slope1.5               
@Recode/0_InChannel_57_into_10_OutChannel_57                                
@Recode/2_InChannel_10_into_248_OutChannel_10/land                            
@Recode/2_InChannel_10_into_1_OutChannel_10/mask_1                            
@-DN7038/781/2bright-DNshp-DN8874/14172-DDshp 

exit

#Following is the simplified AutoCAL CommandLine=====mBPL=2 ======on Tue Nov  5 09:11:47 2013
#/home/ymorel/4SM/compile.sh
nice -20  4SM   -Process  -Origin/DigitalGlobe                       
-DB/cockburnwvmp/50_3s_9S_0s/9_54/11210_14500/362.567_6456.625/1_1           
-Mis/WesternAustralia/CockburnMP/WV02/MultiPan/PIX/UTM_50_008/0.002_-0.002/2_JUN_2012   
-LS/0255.0/255.0/255.0/255.0/255.0/255.0/255.0/255.0/255.0           
@WLM/453.0/508.0/581.0/575.0/627.0/689.0/744.0/890.0/954.0           
@WLm/401.0/447.0/511.0/565.0/588.0/629.0/704.0/772.0/862.0           
@WL/427.0/477.5/574.0/570.0/607.5/669.0/724.0/831.0/908.0           
-cWL/0.500/0.500/0.500/0.500/0.200/0.500/0.900/0.500/0.500           
@CP/0000.0/000.0/000.0/000.0/000.0/000.0/000.0/000.0/000.0_1.33m       
-SCL/00055/00074/00114/00080/00104/00096/00104/00118/00078           
-LsM/201.3/200.8/197.0/295.4/180.1/193.8/191.7/207.9/197.3           
-Lsw/117.8/081.8/043.1/046.3/024.0/018.7/012.2/006.8/006.0           
-dLsw000.0/000.0/001.0/001.0/001.0/000.0/000.0/000.0/000.0           
-Lw/0008.4/004.2/002.5/004.5/000.5/000.0/000.0/000.0/000.0           
-M/@000001/@0002/00003/00004/@0005/@0006/@0007/@0008/00009           
-Lm/0001.0/001.0/000.5/000.2/002.0/001.0/001.0/255.0/255.0           
-KK2_3_0.750/Knir6.200/Bm400_BG510_GR600_RN700/mask_3               
-Z/MSL0.00d/cSL1.5/N_4/NDR_53_1000_zDTM/mask_3                   
-B/TclNe/cNorm5.00/cLM1.00                           
-deglint/vRbaD/q/GlintM0 16.0                           
-extract/v/rawBDH/fullBDH/NIRband9/NIRmax5.0/mBPL3/mask_3           
-calibrate/v/BDSMpzg/plot_BPL/BC_9_8_7_1/BDh_12_13_14               
-calibrate/v/BdSNpzg/plot_BPL/BC_6_5_4_2/BDh_12_23_24               
-Model/v/runWV2/mask_3    -TPsv/2/4/Text/100    -LL10               
-Smooth/5/D/Smart+                               
@ProfileAB/LvW/profile_red/chAB_1_2_3_4_5_6_7_8_9               
@-ProfileAB/ZrvWb/profile_yellow/chAB_1_2_3_4                   
@-E1/11210/12700/14500  @-E1/11210/1/14500                       
@-DN300/100/3/comment @-DDshp @-DNshp  @-DN8874/14172 @-Eshp           
@-Recode/2_InChannel_10_into_3_OutChannel_10/mask_3               
 
 
 



Download
cockburnwvmp_results_2013nov9th.zip (412 MB)
from UbuntuOne at http://ubuntuone.com/5ECQUu46CDBelIz40aCHm1

 
It contains the following:
  • cockburnwvmp_Zcm.tif....S16...one channel......Retrieved depth in centimeters
  • cockburnwvmp_LB.tif......U8....four channels....TOA Water column corrected reflectance, 0-250, for bands Purple, Blue, Green and PAN
This is a twin stripes dataset
I used 4SM to mosaick and balance the following tiles

12JUN02025337-M3XS_R05C1-052870327010_01_P001.TIF and 12JUN02025403-M3XS_R05C1-052870327010_01_P001.TIF
12JUN02025337-M3XS_R06C1-052870327010_01_P001.TIF and 12JUN02025403-M3XS_R06C1-052870327010_01_P001.TIF
12JUN02025337-M3XS_R07C1-052870327010_01_P001.TIF and 12JUN02025403-M3XS_R07C1-052870327010_01_P001.TIF




Georeferencing

of tile 12JUN02025337-M3XS_R05C1-052870327010_01_P001.TIF
BEGIN_GROUP = TILE_5
    filename = "12JUN02025337-M3XS_R05C1-052870327010_01_P001.TIF";
    ULColOffset = 0;
    ULRowOffset = 25000;
    URColOffset = 9167;
    URRowOffset = 25000;
    LRColOffset = 9167;
    LRRowOffset = 31249;
    LLColOffset = 0;
    LLRowOffset = 31249;
    ULLon =    115.54477698;
    ULLat =    -32.01816857;
    URLon =    115.73886499;
    URLat =    -32.02024724;
    LRLon =    115.73731623;
    LRLat =    -32.13297307;
    LLLon =    115.54299006;
    LLLat =    -32.13088533;
    ULX = 362566.99999812;
    ULY = 6456625.00000058;
    URX = 380900.99999930;
    URY = 6456625.00000055;
    LRX = 380900.99999929;
    LRY = 6444127.00000054;
    LLX = 362566.99999810;
    LLY = 6444127.00000057;
END_GROUP = TILE_5



 



Créer un site
Créer un site