!----------------------------------------------------------------------------------! 
!----------------------------------------------------------------------------------!
!                                                                                  !
!        Automated SPOT Ortho-Mosaic WorkFlow [SPOT_AUTOIMG.eas]                   !
!                                                                                  !
!  This script is an automatic workflow used to process several Satellite          !
!  images into one mosaicked image and was demonstrated during a recent PCI        !
!  Webinar. It will generate orthorectified images using an existing geocoded      !
!  image with the AUTOIMG program and then mosaic them together into one           !
!  color balanced output image.                                                    !
!                                                                                  !
!  Inputs for this model are SPOT satellite images with orbital data, a            !
!  digital elevation model and an existing geocoded image for the area             !
!  [Note: the Geocoded image does not need to be from the same sensor, and         !
!  in this demo the blue band of a multispectral image was used]. All data         !
!  came from the Geomatica demo data and the script should also work with          !
!  other similair data sets.                                                       !
!                                                                                  !
!  This script is the EASI equivalent of SPOT_AUTOIMG.mod                          !
!                                                                                  !
!----------------------------------------------------------------------------------! 
!----------------------------------------------------------------------------------! 

!--------------------------------------------------------------------------------- 
! Variables
!--------------------------------------------------------------------------------- 
  local string imagefile, mosaicfile , in_files, out_files
  local string geocoded, dem, refimage, temp_dir, proj
  local mstring dirlist
  local integer i, px, py  
  local $Z
  
!--------------------------------------------------------------------------------- 
! Define location of files and directories necessary for this script
! This is set up this way so that the details can be easily changed to 
! use other images with
!--------------------------------------------------------------------------------- 

in_files = "C:\Program Files\Geomatica_V100\user\input_images"
out_files = "C:\Program Files\Geomatica_V100\user\output_images"
geocoded = "C:\Program Files\Geomatica_V100\Demo\l7_ms.pix"
dem = "C:\Program Files\Geomatica_V100\user\SPOTDEM.PIX"
mosaicfile = "C:\Program Files\Geomatica_V100\user\output_images\SPOTmosaic.pix"
refimage = "C:\Program Files\Geomatica_V100\user\input_images\SPOTLEFT.PIX"
temp_dir = "C:\Program Files\Geomatica_V100\user"
proj = "UTM 11 E000"
px = 10
py = 10

!--------------------------------------------------------------------------------- 
! Clear the EASI window & get the contents of the directory
!--------------------------------------------------------------------------------- 
   
  PRINT @(1 ,1,CLREOS)

  dirlist = getdirectory(in_files)
  let $Z = "\

!--------------------------------------------------------------------------------- 
! Set up the parameters and execute the commands, then continue again until
! all of the input files in the directory have been processed
!--------------------------------------------------------------------------------- 

for i = 1 to f$len(dirlist)
   imagefile = in_files + $Z + dirlist[i]
	
   FILE = imagefile
   MMSEG =
   ORBIT = 2
   DBIC =
   DBIB =
   DBGC =
   GEOFILE = geocoded
   DBCC =
   BACKVAL =
   FILEDEM =
   DBEC =
   NUMGCPS = 40
   MATSCORE =
   SRCHRAD =
   CLOUDTHR =
   REPORT =
   VERVAL =

   run AUTOIMG

   FILE = imagefile
   DBGC = 3
   MMSEG =
   ORBIT = 2
   MODELTYP = "SRIT"
   ORDER =
   REJECT = 4, 1, 1
   REPORT =

   run GCPREFN

   MFILE  = imagefile
   DBGC =
   ORBIT = 2
   MMSEG =
   MAPUNITS = 
   REPORT = 

   run SATMODEL

   MFILE  = imagefile
   DBICLIST = 
   MMSEG  =
   DBIW  =
   FILO  = 
   OUTDIR = out_files
   DBOCLIST =
   FTYPE  =
   FOPTIONS =
   UPLEFT  =
   LORIGHT =
   MAPUNITS = proj
   PXSZOUT = px, py
   FILEDEM = dem
   DBEC = 1
   BACKELEV = -150
   ELEVREF =
   ELEVUNIT =
   ELFACTOR =
   SAMPLING = 4
   RESAMPLE ="NEAR"
   FLSZ =
   SIGMSQ =
   FLTFIL =
   NLOOK =
   IMAGEFMT =
   EDGECLIP =
   MEMSIZE = 512
   REPORT =
 	
   run ORTHO
 
endfor

!--------------------------------------------------------------------------------- 
! Set up the parameters and execute the Auto Mosaic command using the 
! processed ortho images that were produced with the commands above
!--------------------------------------------------------------------------------- 

   MFILE = out_files
   DBICLIST =
   FILO    = mosaicfile
   TILES   =
   TFILE   =
   FTYPE   =
   FOPTIONS =
   DBLUTI  =
   DBLUTO  =
   FILVIN  =
   DBIV    =
   FILVOUT =
   DBOV    =
   CLRMOS  =
   GENERATE = "YES"
   PREVIEW =
   REFIMG1 = refimage
   HSPOT   =
   BALANCE = "ENTIRE"
   CUTMTHD = "MINDIFF"	
   BLEND   =
   BACKVAL =
   TEMPDIR = temp_dir

run AUTOMOS

  PRINT @(1 ,1,CLREOS)
  print "-------------------------------------------------------------------"

  print ""
  print " Your output files are stored in the following directory:"
  print "	", out_files	
  print ""
  print ""
  print @reverse,"                  SPOT_AUTOIMG.eas Script Finished              ",@alloff
  print ""
  print "-------------------------------------------------------------------"

return

! end of script

!---------------------------------------------------------------------------------! 
!---------------------------------------------------------------------------------!