Arcachon tutorial: Modeling 1
back to Arcachon tutorial     back to Calibration 1             next to Calibration 2

 



 
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
 
 


 
 
Model1_1
First modeling: copy
Line_1 Calibration 1          into          Line_2  Model 1

 
#Line_1 Calibration 1+++++++++++Calibration 1+++++++++++Calibration 1+++++++++++
nice -20  ./4SM.4.08 -Process/14_Jul_2011/06_JUL_2011  -OriginSPT @@
-DB/arcachon/27_1_0_0/3_1/500_630/0.010_22.610/1_1 @@
-Mis/Arcachon/France/Spot/XS/NA/UTM_1_008/0.020_0.020/1_AUG_1986 @@
-LS/0255.0/255.0/255.0 @@
-M/@000001/00002/00003 @@
-LsM/078.8/071.2/065.2_cLM=0.450 @@
-Lsw/030.4/014.2/008.0 @@
-Lw/0001.5/000.0/000.0 @@
-Lm/0001.0/001.0/001.0 @@
-KK1_2_0.3900/Knir4.500/mask_1 @@
-Z/MSL0.00/cZ1.000 @@
-B/tclNe/LBref200_100/Bmin0/cLM1.00 @@
-extract/v/mBPL2 @@
-Calibrate/v/BDSNpzg/BDh_12_13 @@
-Model/mask_2 @@
-Smooth/5/D/Smart+ @@

 

#Line_2 Model 1+++++++++++Model 1+++++++++++Model 1+++++++++++
nice -20  ./4SM.4.08 -Process/14_Jul_2011/06_JUL_2011  -OriginSPT @@
-DB/arcachon/27_1_0_0/3_1/500_630/0.010_22.610/1_1 @@
-Mis/Arcachon/France/Spot/XS/NA/UTM_1_008/0.020_0.020/1_AUG_1986 @@
-LS/0255.0/255.0/255.0 @@
-M/@000001/00002/00003 @@
-LsM/078.8/071.2/065.2_cLM=0.450 @@
-Lsw/030.4/014.2/008.0 @@
-Lw/0001.5/000.0/000.0 @@
-Lm/0001.0/001.0/255.0 @@
-KK1_2_0.3900/Knir4.500/mask_1 @@
-Z/MSL0.00/cZ1.000 @@
-B/tclNe/LBref200_100/Bmin0/cLM1.00 @@
-extract/v/mBPL2 @@
-calibrate/v/BDSNpzg/BDh_12_13 @@
-Model/mask_2 @@
-Smooth/5/D/Smart+ @@

 
Specify the command line for Modeling
  • Disable the -Calibrate.. argument??
    • -calibrate/v/BdSNpZg/BDh_12_23
  • Disable the NIR solution
    •  -Lm/0001.0/001.0/255.0 
  • Enable the -Model... argument?
    • -Model/mask_2
  • Enable the -Smooth argument
    • -Smooth/5/D/Smart+
  • Run the script
Run the script
  • Inspect the resulting image_B of the average bottom reflectance in channel_9?
    • see that the breaking waves are represented by festoons of saturated values 
    • we shall mask them into channel 4,
      • so that breaking wave pixels shall be excluded
      • fom the  Extraction, Calibration and Modeling processes
  • Do not bother about the results for now?
    • weshall deal with them in due time
  • At this stage, ?
    • only fairly shallow and bright bottoms are mapped at an acceptable depth
    • results over deeper bottoms would appear to be most questionable/unacceptable
Model 1: average bottom brightness B

 
 

 
Recode Magic
 

#Line_2 Model 1+++++++++++Model 1+++++++++++Model 1+++++++++++
nice -20  ./4SM.4.08 -Process/14_Jul_2011/06_JUL_2011  -OriginSPT @@
-DB/arcachon/27_1_0_0/3_1/500_630/0.010_22.610/1_1 @@
-Mis/Arcachon/France/Spot/XS/NA/UTM_1_008/0.020_0.020/1_AUG_1986 @@
-LS/0255.0/255.0/255.0 @@
-M/@000001/00002/00003 @@
-LsM/078.8/071.2/065.2_cLM=0.450 @@
-Lsw/030.4/014.2/008.0 @@
-Lw/0001.5/000.0/000.0 @@
-Lm/0001.0/001.0/255.0 @@
-KK1_2_0.3900/Knir4.500/mask_1 @@
-Z/MSL0.00/cZ1.000 @@
-B/tclNe/LBref200_100/Bmin0/cLM1.00 @@
-extract/v/mBPL2 @@
-calibrate/v/BDSNpzg/BDh_12_13 @@
-Model/mask_2 @@
-Smooth/5/D/Smart+ @@

-Recode/201_InChannel_9_into_233_OutChannel_4
  • Now we can mask out all pixels saturated at 201
  • Append a -Recode.. argument at the end of the current scritpt
    • -Recode/201_InChannel_9_into_233_OutChannel_4
    • this means: "recode all pixels mapped at 201 in channel 9            ?into 233 in channel 4"
  • Recode Magic!
    • ?all pixels coded at 201 in chanel_9
    • are recoded into 233 in chanel_4 (light blue)
  • From now on, wave breaker pixels shall be
    • ?excluded of the Extraction of calibration data
    • mapped at 233 in the Modeling process

Mask mSE in channel 4

 


 
 
Model1_2
 
Now let's check that the wave breakers are masked properly at 233 (light blue)




#Line_3 Model 1+++++++++++Model 1+++++++++++Model 1+++++++++++
nice -20  ./4SM.4.08 -Process/14_Jul_2011/06_JUL_2011  -OriginSPT @@
-DB/arcachon/27_1_0_0/3_1/500_630/0.010_22.610/1_1 @@
-Mis/Arcachon/France/Spot/XS/NA/UTM_1_008/0.020_0.020/1_AUG_1986 @@
-LS/0255.0/255.0/255.0 @@
-M/@000001/00002/00003 @@
-LsM/078.8/071.2/065.2_cLM=0.450 @@
-Lsw/030.4/014.2/008.0 @@
-Lw/0001.5/000.0/000.0 @@
-Lm/0001.0/001.0/255.0 @@
-KK1_2_0.3900/Knir4.500/mask_1 @@
-Z/MSL0.00/cZ1.000 @@
-B/tclNe/LBref200_100/Bmin0/cLM1.00 @@
-extract/v/mBPL2 @@
-calibrate/v/BDSNpzg/BDh_12_13 @@
-Model/mask_2 @@
-Smooth/5/D/Smart+ @@

@Recode/201_InChannel_9_into_233_OutChannel_4
 


Please disable the -Recode... argument : @Recode...
and run the script
All wave breaker pixels
are masked at 233 (light blue).


Now we want to
  • extract a "clean" calibration data textfile
  • perform a "clean" calibration
  • perform a "clean" modeling