| 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
- 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
|
| -DB/... -DB/tmnov/31_3s_0s_0s/4_1/1265_1560/208.350_2398.770 - syntax_DB
- We want to create a subset database named tmnov.pix
- The structure of the database for a U8 4 bands image is
- 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
|
| -Mis/... -Mis/Bahamas/Caicos/Landsat/TM/FUGRO-NPA/UTM_19_008/0.030_0.030/22_NOV_1990 - syntax_Mis
- Name is Bahamas
- Area is Caicos
- Sensor is Landsat
- Mode is TM
- Archive is FUGRO-NPA
- Datum:
- Projection is UTM
- UTMzone is 19
- Ellipsoid is 8
- ==> .../UTM_19_008/...
- Pixel size in kilometers is 0.03 km
- The scene was acquired on November 22nd 1990
|
| -Import/... -Import/data*TM_009_045_1990-11-22/dbnc_6_0_0_0/R4341_L4462/Origin_129.510_2448.060_0_0/chIn1,4/chOut1,4/import - 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
- this original image has 4341 rows and 4462 lines
- 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 ./import.sh
|
|
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
|
