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

Model2_1
Copy Line 4 into Line 5
 
#Line_4 Calibration 2+++++++++++Calibration 2+++++++++++Calibration 2+++++++++++
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 @@
-Lw/0002.5/000.0/000.0 @@
-Lm/0001.0/001.0/001.0 @@
-KK1_2_0.420/Knir4.500/mask_1 @@
-Z/MSL0.00/cZ1.000 @@
-B/tclNe/LBref200_100/Bmin0/cLM1.00 0.70 @@
-Extract/v/mBPL2 @@
@Extract/v/RawBDH/NIRband3/NIRmax5.0/FullBDH/mSOIL21/mBPL2 @@
-Calibrate/v/BdSNpzg/BDh_12_13 @@
-Model/mask_2 @@
-smooth/5/D/Smart+ @@
@Recode/201_InChannel_9_into_233_OutChannel_4

 

#Line_5 Model 2+++++++++++ Model 2+++++++++++ Model 2+++++++++++
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 @@
-Lw/0002.5/000.0/000.0 @@
-Lm/0001.0/001.0/001.0 @@
-KK1_2_0.420/Knir4.500/mask_1 @@
-Z/MSL0.00/cZ1.000 @@
-B/tclNe/LBref200_100/Bmin0/cLM1.00 0.70 @@
@Extract/v/mBPL2 @@
@Extract/v/RawBDH/NIRband3/NIRmax5.0/FullBDH/mSOIL21/mBPL2 @@
@Calibrate/v/BdSNpzg/BDh_12_13 @@
-Model/mask_2 @@
-smooth/5/D/Smart+ @@
@Recode/201_InChannel_9_into_233_OutChannel_4

 
  • Disable @Extract.... argument
  • Disable -Calibrate... argument
Run the command line script
shallow pixels are modeled by the Red solution
Note
  • -Lm/0001.0/001.0/255.0 :
    • very low thresholds for the green and red bands
    • very high threshold  for the NIR band
      • prevent any depth to be computed using the NIR solution
  • ?-M@000001/00002/00003 : shallow pixels are modeled by the Red solution
 


Unacceptable
Average bottom brightness results  are not acceptable
 

B average bottom brightness
-Lm/0001.0/001.0/255.0
  • As the Lm thresholds are very low, bottom brightness is seen to get darker as the depth increases
    • this is not acceptable in a "all sand" highly dynamic environment
  • Increasing from Lswred=14.2 to Lswred=15.2 
    • would be more representative of the waters in the main channel
    • would not solve any of our problems though
    • Try it NOW, see for yourself

wZ
  • Brown:......land pixels are coded at 248
  • Blue:.........deep water pixels
  • Cyan:........wave breaker pixels are coded at 233
  • Green:......all other shallow pixels
    • modeled using the Red solution

 
Bmin
We can impose a minimum average brightness B
  • Let's raise the threshold on average bottom brightness from Bmin=10 to Bmin=100
    • -Lm/0001.0/001.0/255.0
    • -M@000001/00002/00003
    • -B/tclNe/Bmin100/cLM1.00
  • Find your preferred value for Bmin


Quite neat!
  • On display now is B>=100 results 
    • as B<100 results have been ruled out




 
The image below provides the explanation

 


How threshold Bmin works
 
  • Left: the image is a false color composite of linearized data
    • ?with overlay of a red shapefile
  • Right: this plot of linearized Green vs Red pixels under the red shapefile shows 
    • a neat and curved Soil Line :
      • the water volume reflectance Lswgreen cannot be ignored :
      • its value must be set so that the Soil Line model model fits the Soil Line pixels
      • failing that, authors report that computed depth over dark bottoms are badly under-estimated
    • a fairly neat and linear Brightest Pixels Line
    • the scatter of foul water pixels plot below the red Bmin line
 
  • How Bmin works : 
  • pixels that plot below the red Bmin line shall 
    • either be mapped as deep water (see above)
    • or processed with the "one-band case" (see below) 

 



Créer un site
Créer un site