[gdal-dev] Shift of geoTiff
Paul Meems
bontepaarden at gmail.com
Mon Jun 7 10:15:03 EDT 2010
One of our user of MapWindow GIS noticed a shift of its geotiff by one pixel
width.
I took a look at it to see if the problem is with MapWindow or GDAL, because
MapWindow uses GDAL internally.
The provided metadata of the geotiff is:
GROUP = L1_METADATA_FILE
GROUP = METADATA_FILE_INFO
ORIGIN = "Image courtesy of the U.S. Geological Survey"
REQUEST_ID = "0101004096683_00045"
PRODUCT_CREATION_TIME = 2010-04-12T00:11:01Z
STATION_ID = "EDC"
LANDSAT7_XBAND = "2"
GROUND_STATION = "EDC"
LPS_PROCESSOR_NUMBER = 1
DATEHOUR_CONTACT_PERIOD = "0226604"
SUBINTERVAL_NUMBER = "08"
END_GROUP = METADATA_FILE_INFO
GROUP = PRODUCT_METADATA
PRODUCT_TYPE = "L1T"
ELEVATION_SOURCE = "GLS2000"
PROCESSING_SOFTWARE = "LPGS_11.0.0"
EPHEMERIS_TYPE = "DEFINITIVE"
SPACECRAFT_ID = "Landsat7"
SENSOR_ID = "ETM+"
SENSOR_MODE = "SAM"
ACQUISITION_DATE = 2002-09-23
SCENE_CENTER_SCAN_TIME = 00:09:22.3202747Z
WRS_PATH = 95
STARTING_ROW = 85
ENDING_ROW = 85
BAND_COMBINATION = "123456678"
PRODUCT_UL_CORNER_LAT = -35.0709086
PRODUCT_UL_CORNER_LON = 140.3495920
PRODUCT_UR_CORNER_LAT = -35.0547004
PRODUCT_UR_CORNER_LON = 143.0866731
PRODUCT_LL_CORNER_LAT = -37.0344333
PRODUCT_LL_CORNER_LON = 140.3332320
PRODUCT_LR_CORNER_LAT = -37.0170222
PRODUCT_LR_CORNER_LON = 143.1391164
PRODUCT_UL_CORNER_MAPX = 440700.000
PRODUCT_UL_CORNER_MAPY = -3881100.000
PRODUCT_UR_CORNER_MAPX = 690300.000
PRODUCT_UR_CORNER_MAPY = -3881100.000
PRODUCT_LL_CORNER_MAPX = 440700.000
PRODUCT_LL_CORNER_MAPY = -4098900.000
PRODUCT_LR_CORNER_MAPX = 690300.000
PRODUCT_LR_CORNER_MAPY = -4098900.000
PRODUCT_SAMPLES_PAN = 16641
PRODUCT_LINES_PAN = 14521
PRODUCT_SAMPLES_REF = 8321
PRODUCT_LINES_REF = 7261
PRODUCT_SAMPLES_THM = 8321
PRODUCT_LINES_THM = 7261
BAND1_FILE_NAME = "L71095085_08520020923_B10.TIF"
BAND2_FILE_NAME = "L71095085_08520020923_B20.TIF"
BAND3_FILE_NAME = "L71095085_08520020923_B30.TIF"
BAND4_FILE_NAME = "L71095085_08520020923_B40.TIF"
BAND5_FILE_NAME = "L71095085_08520020923_B50.TIF"
BAND61_FILE_NAME = "L71095085_08520020923_B61.TIF"
BAND62_FILE_NAME = "L72095085_08520020923_B62.TIF"
BAND7_FILE_NAME = "L72095085_08520020923_B70.TIF"
BAND8_FILE_NAME = "L72095085_08520020923_B80.TIF"
GCP_FILE_NAME = "L71095085_08520020923_GCP.txt"
METADATA_L1_FILE_NAME = "L71095085_08520020923_MTL.txt"
CPF_FILE_NAME = "L7CPF20020701_20020930_06"
END_GROUP = PRODUCT_METADATA
GROUP = MIN_MAX_RADIANCE
LMAX_BAND1 = 191.600
LMIN_BAND1 = -6.200
LMAX_BAND2 = 196.500
LMIN_BAND2 = -6.400
LMAX_BAND3 = 152.900
LMIN_BAND3 = -5.000
LMAX_BAND4 = 241.100
LMIN_BAND4 = -5.100
LMAX_BAND5 = 31.060
LMIN_BAND5 = -1.000
LMAX_BAND61 = 17.040
LMIN_BAND61 = 0.000
LMAX_BAND62 = 12.650
LMIN_BAND62 = 3.200
LMAX_BAND7 = 10.800
LMIN_BAND7 = -0.350
LMAX_BAND8 = 243.100
LMIN_BAND8 = -4.700
END_GROUP = MIN_MAX_RADIANCE
GROUP = MIN_MAX_PIXEL_VALUE
QCALMAX_BAND1 = 255.0
QCALMIN_BAND1 = 1.0
QCALMAX_BAND2 = 255.0
QCALMIN_BAND2 = 1.0
QCALMAX_BAND3 = 255.0
QCALMIN_BAND3 = 1.0
QCALMAX_BAND4 = 255.0
QCALMIN_BAND4 = 1.0
QCALMAX_BAND5 = 255.0
QCALMIN_BAND5 = 1.0
QCALMAX_BAND61 = 255.0
QCALMIN_BAND61 = 1.0
QCALMAX_BAND62 = 255.0
QCALMIN_BAND62 = 1.0
QCALMAX_BAND7 = 255.0
QCALMIN_BAND7 = 1.0
QCALMAX_BAND8 = 255.0
QCALMIN_BAND8 = 1.0
END_GROUP = MIN_MAX_PIXEL_VALUE
GROUP = PRODUCT_PARAMETERS
CORRECTION_METHOD_GAIN_BAND1 = "CPF"
CORRECTION_METHOD_GAIN_BAND2 = "CPF"
CORRECTION_METHOD_GAIN_BAND3 = "CPF"
CORRECTION_METHOD_GAIN_BAND4 = "CPF"
CORRECTION_METHOD_GAIN_BAND5 = "CPF"
CORRECTION_METHOD_GAIN_BAND61 = "CPF"
CORRECTION_METHOD_GAIN_BAND62 = "CPF"
CORRECTION_METHOD_GAIN_BAND7 = "CPF"
CORRECTION_METHOD_GAIN_BAND8 = "CPF"
CORRECTION_METHOD_BIAS = "IC"
BAND1_GAIN = "H"
BAND2_GAIN = "H"
BAND3_GAIN = "H"
BAND4_GAIN = "L"
BAND5_GAIN = "H"
BAND6_GAIN1 = "L"
BAND6_GAIN2 = "H"
BAND7_GAIN = "H"
BAND8_GAIN = "L"
BAND1_GAIN_CHANGE = "0"
BAND2_GAIN_CHANGE = "0"
BAND3_GAIN_CHANGE = "0"
BAND4_GAIN_CHANGE = "0"
BAND5_GAIN_CHANGE = "0"
BAND6_GAIN_CHANGE1 = "0"
BAND6_GAIN_CHANGE2 = "0"
BAND7_GAIN_CHANGE = "0"
BAND8_GAIN_CHANGE = "0"
BAND1_SL_GAIN_CHANGE = 0
BAND2_SL_GAIN_CHANGE = 0
BAND3_SL_GAIN_CHANGE = 0
BAND4_SL_GAIN_CHANGE = 0
BAND5_SL_GAIN_CHANGE = 0
BAND6_SL_GAIN_CHANGE1 = 0
BAND6_SL_GAIN_CHANGE2 = 0
BAND7_SL_GAIN_CHANGE = 0
BAND8_SL_GAIN_CHANGE = 0
SUN_AZIMUTH = 49.0012493
SUN_ELEVATION = 41.9440742
OUTPUT_FORMAT = "GEOTIFF"
END_GROUP = PRODUCT_PARAMETERS
GROUP = CORRECTIONS_APPLIED
STRIPING_BAND1 = "NONE"
STRIPING_BAND2 = "NONE"
STRIPING_BAND3 = "NONE"
STRIPING_BAND4 = "NONE"
STRIPING_BAND5 = "NONE"
STRIPING_BAND61 = "NONE"
STRIPING_BAND62 = "NONE"
STRIPING_BAND7 = "NONE"
STRIPING_BAND8 = "NONE"
BANDING = "N"
COHERENT_NOISE = "Y"
MEMORY_EFFECT = "N"
SCAN_CORRELATED_SHIFT = "N"
INOPERABLE_DETECTORS = "N"
DROPPED_LINES = "N"
END_GROUP = CORRECTIONS_APPLIED
GROUP = PROJECTION_PARAMETERS
REFERENCE_DATUM = "WGS84"
REFERENCE_ELLIPSOID = "WGS84"
GRID_CELL_SIZE_PAN = 15.000
GRID_CELL_SIZE_THM = 30.000
GRID_CELL_SIZE_REF = 30.000
ORIENTATION = "NUP"
RESAMPLING_OPTION = "CC"
MAP_PROJECTION = "UTM"
END_GROUP = PROJECTION_PARAMETERS
GROUP = UTM_PARAMETERS
ZONE_NUMBER = 54
END_GROUP = UTM_PARAMETERS
END_GROUP = L1_METADATA_FILE
END
When I perform a GDALINFO on the GeoTiff file, I get this report:
Driver: GTiff/GeoTIFF
Files: L71095085_08520020923_B10.TIF
Size is 8321, 7261
Coordinate System is:
PROJCS["WGS 84 / UTM zone 54N",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",141],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","32654"]]
Origin = (440699.999999999940000,-3881100.000000000000000)
Pixel Size = (30.000000000000000,-30.000000000000000)
Metadata:
AREA_OR_POINT=Point
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 440700.000,-3881100.000) (140d20'58.53"E, 35d 4'15.27"S)
Lower Left ( 440700.000,-4098930.000) (140d19'59.63"E, 37d 2'4.93"S)
Upper Right ( 690330.000,-3881100.000) (143d 5'13.21"E, 35d 3'16.90"S)
Lower Right ( 690330.000,-4098930.000) (143d 8'22.06"E, 37d 1'2.23"S)
Center ( 565515.000,-3990015.000) (141d43'38.60"E, 36d 3'8.96"S)
Band 1 Block=8321x1 Type=Byte, ColorInterp=Gray
As you can see the lower right Y should be -4098900, but GDALINFO reports
,-4098930
My question is:
Is this the correct behavior of GDAL or do we misinterpreted it?
Thanks,
Paul
--
Paul Meems
Release manager, configuration manager
and forum moderator of MapWindow GIS.
www.mapwindow.org
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20100607/d9836a1a/attachment-0001.html
More information about the gdal-dev
mailing list