tmnov tutorial : Prepare your image
return to tmnov tutorial         next         next to Model 1 
Prepare your image for the AutoCalibration process
tmnov tutorial: the AutoCAL process
now takes care of almost everything for you:
  • Prepare shapefiles
    • nothing to do here, unless you want to experiment with any particular shapefile
  • Prepare your script
  • Run the AutoCAL script
    • Create the database tmnov.pix
    • Compute glint regressions
    • Write the special effects mask in channel 5
    • Extract the calibration data
    • Map the BPL pixels
    • Perform an automatic calibration and displays calibration diagram
    • Deglint the image (does not apply here)
  • Output of the AutoCAL process
    • Display a calibration diagram
    • Write and display the 4SM command line for the next step
Please use OpenEV to
  • Export 
    • data/TM_009_045_1990-11-22.tif
      • into PCIDSK format:
    • data/TM_009_045_1990-11-22.pix 
    • this is a six U8 spectral bands TM image
    • we only shall import/use the following spectral bands
      • 1 Blue   TM1
      • 2 Green TM2
      • 3 Red     TM3
      • 4 NIR    TM4
  • Display : 
    • use OpenEV 
      • to display a FCC  and a TCC of TM_009_045_1990-11-22.pix
      • with overlay of seatruth depth points seatruth/depth_points_reproject.shp
    • this is a six U8 bands image  ==> /dbnc_6_0_0_0/      in -Import... argument

Now, customize the AutoCAL command line
  • In a real life scenario, you would modify an existing command line for a Landsat TM image,
    • rather than create a new one from scratch
  • Raw data is stored in ~/4SM_tutorials/4sm_tutorial_tmnov/data/TM_009_045_1990-11-22.pix
  • Customize the parameters in -DB..., -Mis..., and -Import... arguments
  • syntax_DB
  • We want to create a subset database named tmnov.pix
    • ==> .../tmnov/...                                                       
  • The structure of the database for a U8 4 bands image is
    • ==> .../31_3s_0s_0s/...                                               
  • Bands: there shall be 4 raw data channels: Blue=TM1, Green=TM2, Red=TM3 and NIR=TM4
    • and the FirstRawChannel shall be channel_1
    • ==> .../4_1/...          
  • At this stage, the coordinates of the UL corner of the working database are not known
    • they shall be elaborated by 4SM from reading the import.shp shapefile
      • .../1265_1560/...             nbRow and nbLin of the working database tmnov.pix
      • .../208.350_2398.770/... West and North are UTM coordinates of the UL corner of the UL pixel of the database??
  • 4SM shall not proceed until the -DB... argument is correct                    Syntax_Import
    • either for a full image if no ./import.shp shapefile is available
    • or for just a subset of it, specified by shapefile ./import.shp
  • Georeferencing :
    • in Remote Sensing, the origin of an image is the UL corner of the UL pixel
  • syntax_Mis
  • Name is Bahamas
    • ==> .../Bahamas/...                                                  
  • Area is Caicos
    • ==> .../Caicos/...                                                     
  • Sensor is Landsat
    • ==> .../Landsat/...                                                    
  • Mode is TM
    • ==> .../TM/...                                                            
  • Archive is FUGRO-NPA
    • ==> .../FUGRO-NPA/...                                          
  • Datum:
    • Projection is UTM
    • UTMzone is 19
    • Ellipsoid is 8
    • ==> .../UTM_19_008/...
  • Pixel size in kilometers is 0.03 km
    • ==> .../0.030_0.030/...                                
  • The scene was acquired on November 22nd 1990
    • ==>.../22_NOV_1990/...   
  • Syntax_Import
  • database_In: the original image is data*TM_009_045_1990-11-22.pix
    • remember that data*TM... stands for data/TM...
  • this image has 6 U8 channels   
    • .../dbnc_6_0_0_0/...
  • this original image has 4341 rows and 4462 lines 
    • .../R4341_L4462/...
  • this image's UTM origin in kilometers is 
    • .../Origin_129.510_2448.060_0_0/...
      • Origin_West =  129.510 km
      • Origin_North=2448.060 km
      • False_Easting=False_Northing=0
  • shall import from channels 1 to 4 
    • .../chIn1,4/... of original image data*TM_009_045_1990-11-22.pix
  • shall................write to channels 1 to 4 
    • .../chOut1,4/...  of working database tmnov.pix
  • shall use shapefile ./import.shp
    • .../import
    • draw a rectangular polygon to specify your area of interest
    • and save it in your work directory as ./

Your updated AutoCAL command line arguments
  • -DB/tmnov/31_3s_0s_0s/4_1/1265_1559/208.350_2398.770
    • the georeferencing in the -DB... argument shall be computed and formated by the AutoCAL process 
    • 4SM won't proceed until you have updated the -DB... argument as required
  • -Mis/Bahamas/Caicos/Landsat/TM/FUGRO-NPA/UTM_19_008/0.030_0.030/22_NOV_1990
  • -Import/v/dTM/data*TM_009_045_1990-11-22/dbnc_6_0s_0s_0s/R4341_L4462/Origin_129.510_2448.060_0_0/chIn1,4/chOut1,4/import

shapefile import.shp is shown as a white rectangle
seatruth depth points are shown
Now you need to draw several shapefiles inside your area of interest import.shp
and save them in the work directory,
ready for use in the AutoCAL process:
  • baddata.shp
  • Lsw.shp
  • dLsw.shp
  • LsM.shp
  • glint.shp
  • vegetation.shp


Now, please run the AutoCalibration process

#Line_0 - create the database tmnov.pix and run the AutoCalibrator
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.030_0.030/22_NOV_1990 @@
-MakePIX/Import @@
-Import/v/dTM/data*TM_009_045_1990-11-22/dbnc_6_0S_0s_0s/R4341_L4462/Origin_129.510_2448.060_0_0/chIn1,4/chOut1,4/import @@

Créer un site
Créer un site