4SM SPOT arcachon model

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
  - saturation is coded at 201
  - see legend for image B
  - code 201 in light blue color
- 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

Please note that the -Recode... argument takes precedence over other actions

Syntax_Recode

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