tmnov tutorial : 
Prepare your image
 
 
 
Following is a copy/paste of the AutoCal.txt text file
Copy it into your script,
replace @@ by two antislashes,
then run the Calibration

#Line_1 Model_1=============Model_1=============Model_1=============
#Following is the simplified AutoCal CommandLine=====mBPL=2 ======on Sat Nov 25 05:01:20 2017
nice -20  ./4SM.8.02 -Process/25_Nov_2017  -Origin/FUGRO-NPA @@
-DB/tmnov/31_3s_0s_0s/4_1/1265_1560/208.350_2398.770 @@
-Mis/Bahamas/Caicos/Landsat/TM/FUGRO-NPA/bOA/UTM_+19_008/0.0300_0.0300/22_NOV_1990 @@
@WLM/518.0/610.0/693.0/905.0 @@
@WLm/452.0/529.0/624.0/776.0 @@
@WL/0485.0/569.5/658.5/840.5 @@
-Lsw/060.5/017.4/011.1/006.9_Lsw_AutoCal @@
@dLsw/0-1.3/00-0.5/00-0.3/00-0.1_dLsw @@
-dLsw/000.0/000.0/000.0/000.0 @@
-Lw/0010.0/001.5/000.0/000.0_Lw_AutoCal @@
-M/@000001/00002/00003/00004_M @@
-Lm/0001.0/001.0/001.0/255.0_Lm @@
-LsM/203.6/112.8/136.1/111.3_LsM_AutoCal @@
@CP/192.58/99.02/87.90/10.43_0.75m @@
-KK1_2_0.5141/Knir4.5/CoefK_1.0/mask_3 @@
-Z/NA0.00d/N_1/nDR_34_1000_zDTM/mask_3 @@
-B/tclNe/cNorm2.50/Bmin0/cLM1.00 @@
-deglint/vRbaD/GlintM25 255 @@
@AutoCAL/R2_0.95/Land_-5.0_1.18_4s/GlintM25 @@
-extract/v/rawBDH/FullBDH/NIRband4/NIRmax255.0/mBPL2 /mask_3 @@
-Calibrate/voAI/BDSNpZG/BDh_12_13_34 @@
-Model/v/rt1.00/dTM/ZM30/mask_3 @@
-smooth/5/D/Smart+ -LL10 @@
@-DN300/100/3/comment @-DDshp @-DNshp @-Eshp @-E/1/1265/1/1560 @@
 
exit
Note:
arguments in italic are ignored, may be removed
Please run this script as it is, straight from the AutoCalibration:
  • The AutoCalibration has done a pretty good job
    • not much need to be changed/improved
    • let's show off a bit
  • Next we shall proceed to modeling straight away:
    • we shall just have to disable @Calibrate/....
    • then run the script
    • come on: do it now!


 
-DN926/340
Syntax_DN
  • Please enable a last argument -DN926/340
    • and run the script, see what happens
  • This reads the raw data (DNs) at Row_926 and Line_340
    • deglint and smooth as appropriate
    • then prints the results on screen

M_002 1  2   926   340 204.0 112.0 127.0  88.0

This shows/proves that the calibration data for BPL are real georeferenced pixels.
Handy to check on any particular pixel in the image.




 
-DD926/340
  • Please change for -DD926/340,
    • and run the script, see what happens
  • This models the pixel located at Row_926 and Line_340
    • and displays a detailed diagram of the water column correction process
    • see Syntax_D
On a scale 0-200:  
  • Lw: the Cyan profile denotes
    • the BOA water volume reflectance
  • L: the Blue profile denotes:
    • the BOA spectral radiance at Row_926 Line_340
  • LB: the Red profile denotes
    • the BOA water column corrected reflectance
 
  • Retrieved depth is 0.13 m
  • LBblue~=LBgreen~=LBred~=200
    • this is the result of the optimization of LsM in the AutoCal process



 
-DD300/500
This models the pixel at Row_300/Line_500 and displays a diagram of the water column correction process
The choice of -Lm... threshold values is very important           see Syntax_Lm
-DN300/500
8U Ls   300  500 /118.0/041.0/012.0/006.0
  • Please change for -DN300/500
    • and run the script (-smooth/... is disabled)
    • 8U Ls   300  500 /118.0/041.0/012.0/006.0
    • Lsred=12.0
    • Lsred-Lswred=12.0-11.1=0.9

===================================================
===================================================
 
  • Now change for -DD300/500
    • and run the script (-smooth/... is disabled)
    • Lr=(Lsred-Lswred)*CNred=1.4
    • Lm*CNred=1.6
    • Lr<Lm ==>Z_GreenSolution=5.03 m
  • Now, enable -Smooth/....
    • and run the script
    • Lr<Lm ==>Z_RedSolution=6.61 m
  • Now, change Lmred  from 1 to 3: 
    • -Lm/0001.0/001.0/003.0/255.0_Lm
    • and run the script
    • Lr<Lm ==>Z_GreenSolution=5.09 m
  • Now enable -Model/v/Rt1.00/....           Syntax_Model
    • and run the script
    • Lr<Lm ==>Z_GreenSolution=5.74 m
 -DD300/500
Lso   118.0  41.0  12.0   6.0
 
-Lm/0001.0/001.0/001.0/255.0 
Red solution: Z=6.61 m
  • The Red solution should be prefered,
  • or the PAN solution when available.
=============================================
-Lm/0001.0/001.0/003.0/255.0 
Green Solution: Z=5.09<m
===================================================
  • The last trick  -Model/v/Rt1.00/....
    • is an attempt to bridge that gap
 
 
 

 



Créer un site
Créer un site