User's Guide

INPEDDR

Uses a header file created by FSTFMTIN to update a LAS DDR with valid projection corners.

Function:

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.

Parameters:

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).

Examples:

  1. LAS> inpeddr in=brazil

    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.

  2. LAS> inpeddr in=inimage sband=1

    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.

Description/Algorithm:

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) + 1

The 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.

Nonfatal Error Message:

  1. [inpeddr-calczone] The zone code has been calculated

    The zone code was not specified in the Fast Format header. An estimate zone code has been calculated and assigned.

  2. [inpeddr-calcangle] The oerinetaion angle has been calculated The orientation angle was not specified in the Fast Format header. An estimate orientation angle has been calculated and assigned.

Fatal Error Messages:

  1. [inpeddr-alloc] Error allocating dynamic memory

    An error was detected while trying to allocate dynamic memory. If the error persists, contact the system administrator.

  2. [inpeddr-fatal] Fatal error encountered

    A fatal error was encountered during processing. The error message that is displayed immediately preceding this message is the specific error encountered. Processing terminates.

  3. [inpeddr-getcor] Incorrect corners from scene

    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.

  4. [inpeddr-multi] Multiple images, sub-bands and windows not allowed

    Multiple images, sub-bands and windows may not be specified in INPEDDR (See User note 2). Re-run INPEDDR with the full image.

  5. [inpeddr-process] Invalid <xxxxx> in Fast Format Header

    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.

  6. [inpeddr-proj] Projection in Fast Format header is not UTM

    INPE tapes are expected to be in the UTM projection. Other projections are not supported. Verify that the image is from an INPE tape.

  7. [inpeddr-zone] Zone does not match corners

    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.

User Notes:

  1. INPEDDR utilizes images from the National Institute of Space Research located in Brazil. The actual name is Instituto Nacional de Pesquisas Espaciais (INPE).

  2. The original header must be retrieved from tape by FSTFMTIN. The full image (no windowing) from the Fast Format INPE tape must be ingested by FSTFMTIN, however, sub-banding is acceptable.

  3. INPEDDR, along with FSTFMTIN, handle the version A Fast Format header located on INPE tapes. FSTFMTIN will convert the header to the version B Fast Format header. Sections of the headers may be missing such as the USGS PROJECTION # (field 34), the USGS MAP ZONE (field 36), and the USGS PROJECTION PARAMETERS (field 38). INPEDDR will update the DDR with valid projection information.

  4. Images processed at INPE contain the projection coordinates of each corner of the data that was acquired by the satellite and not the corners of the entire image. Band SBAND is examined to locate the corners of the data acquired by the satellite.

    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.