_____________________________________________________________________
IN
INTER
Sub to interpolate maps according to inter- and intra-map arshifts
Each map to be interpolated has a specifimap shift
There are common offsets for scans within a map
NX LENGTH OF MAP IN X
NY LENGTH OF MAP IN Y
OFFSET INTRA-MAP OFFSETS ie relative shifts from row to row
MAPS NOS. OF MAPS TO BE INTERPOLATED
NEND TOTAL NO. OF MAPS TO BE INTERPOLATED ( < NMAP)
______________________________________________________________________
TYPE OPTIE 24-OCT-1981 11:51:33.75
>>>IN
OPTION IN
INTER, VERSION 1
TYPE FILENAME (<CR> FOR DEFAULT)
FILI=VANMAN
There is one map in the file
Reading from the file
1 rows of 2040 pixels in R4 format
TYPE UNIT NO FOR OFFSETS (1 FOR KEYBOARD)
OFF=1
TYPE OFFSETS FOR INDIVIDUAL SPECTRA
OFFSETS=25.2 ! there should be as many as rows shifted
TYPE MAPNUMBERS TO BE INTERPOLATED
MAPI=1
TYPE OUTPUT MAPNUMBERS (def = MAPI)
MAPO=2
ALL SET ?
OK=Y
File extended to map 2
MAP NUMBER 1
TYPE SPECIFIC SHIFT FOR MAP
MAPSHIFT=22.2 ! if there was more than one row each
! row (J) would have a shift of*
! mapshift + offset(j)
_____________________________________________________________________