Uses a header file created by FSTFMTIN to update a LAS DDR with valid projection corners.
Updates the DDR for a LAS image ingested from an Instituto Nacional de Pesquisas Espaciais (INPE) Fast Format tape. Fast Format tapes created at INPE are missing the required projection information. This function will examine the ;FHDR file ingested by FSTFMTIN and determine the projection information necessary to update the DDR.
- IN
- Input image. The name of the LAS image ingested by FSTFMTIN. The image must not have been sub-windowed, however, it may have been sub-banded. The full image from FSTFMTIN is required to create a valid DDR from the associated ;FHDR file. Multiple images, windows, and sub-bands are not allowed in INPEDDR.
- SBAND(2)
- Source band. The physical band to use when determining the line/sample corners of the image. Any band may be used, however, Landsat band 1 should be avoided to eliminate erroneous results in calculating the file corners (See User note 4).
Updates BRAZIL;DDR with information retrieved from
BRAZIL;FHDR. By default, band 2 will be used to find
image corners, therefore, IN must contain at least two
bands.
Updates INIMAGE;DDR with information retrieved from INIMAGE;FHDR. The physical band 1 will be used to find image corners. To avoid complications when calculating the file corners, band 1 of INIMAGE should not be the Landsat band 1.
The projection information and the projection corners in geographic space (latitude and longitude) in the Universal Transverse Mercator (UTM) projection space, pixel size, and orientataion angle of the scene are read from the header. The projection is checked to ensure it is UTM. If the UTM zone code is not included in the header, a zone is estimated by approximating the center longitude and calculated as follows:
zone code = ((center longitude + 180) / 6) + 1The estimated UTM zone code is tested by converting the corner latitudes and longitudes found in ;FHDR to the UTM projection space and comparing these values to the UTM projection corners in the header. If the values do not differ by more than 5.0 meters, the estimated zone is assumed to be the actual zone. Otherwise, the zones one larger and one smaller than the original zone are examined. The projection corners of the new zones are converted and compared again. If the corner projection values of either of the new zones differ by less than 5.0 meters, that zone is assumed to be the actual zone. If neither of the new zones is correct, the program will exit with a fatal error.
When the zone code is included in the header, the corner coordinates will be converted as before and tested for accuracy. This is done as a quality assurance measure for the INPE tapes, If the converted corners do not match the projection corners in the header, the program will exit with a fatal error.
The projection corners just found represent the image data acquired by the satellite and not the fill or file data. The projection corners of the file must be determined. If the orientation angle is not included in the .FHDR then it is calculated. The slope and y intercept for the line between the upper left and lower right corners are found. Then the slope and y intercept for the line beween the lower left and upper right corners are found. Using the slopes and intercepts the center projection information is calculated.
The projection information and file corners are entered into the image DDR. The spheroid is the Australian National which is not supported for the UTM projection and therefore the projection will be set to Transverse Mercator (TM) and the appropriate TM projection information will be specified. The UTM zone code will be filled in for future reference.
The zone code was not specified in the Fast Format header. An estimate zone code has been calculated and assigned.
An error was detected while trying to allocate dynamic memory. If the error persists, contact the system administrator.
A fatal error was encountered during processing. The error message that is displayed immediately preceding this message is the specific error encountered. Processing terminates.
One of the following happened: the upper left x coordinate = the lower right x coordinate, the lower left x coordinate = the upper right x coordinate, or the slope between the upper left and lower right scene corners equals the slope between the lower left and upper right scene corners. Any of these problems would have caused a divide by zero error in calculating the center coordinates of the scene. The contents of the header will have to be examined to resolve this problem. Contact the system administrator.
Multiple images, sub-bands and windows may not be specified in INPEDDR (See User note 2). Re-run INPEDDR with the full image.
The <xxxxx> field in the Fast Format header was blank or contained an illegal value. This information is necessary to update the LAS DDR. Contact the system administrator after confirming that the image was ingested by FSTFMTIN from an INPE tape.
INPE tapes are expected to be in the UTM projection. Other projections are not supported. Verify that the image is from an INPE tape.
The zone code in the Fast Format Header does not correspond to the corner values in the Fast Format Header. The contents of the header will have to be examined to resolve this problem. Contact the system administrator.
Landsat band one of the INPE images has intermittent lines on the right edge of the image data. This may produce erroneous results when determining the corners of the image data. Any other bands will work effectively.