One more thing, although the data should be NAD83, QGIS is interpreting them as an unknown CRS.
I do get a nice black-and-white raster which disappears when I reproject to WGS84. Metadata does survive the process but won't go through the GeoTIF to INF conversion.
Below is the USGS TIF file's data found in QGIS layer properties; maybe someone else can make more sense out of this than I can. Additionally, I switched from a relatively flat terrain in Missouri to Lake Tahoe, California. That way I can more easily see the difference between original mesh and my high-definition scenery.
Thanks,
Chris
< Table not included in quote >
Hi Chris:
You are merely seeing the internal Metadata for the GIS data file in the format that it is distributed by
ex: USGS;
change it to EPSG:4326.
Blank workspace task sessions of most GIS applications open a data file in the GIS projection format (CRS / SRS) specified by its Metadata.
Normally, all subsequent GIS data files opened into that work-space configuration display in the same CRS as the
first such data file opened.
If you follow Dick's steps specified above
you will re-project the source data to EPSG:4326 when it is written out to disk by QGIS.
FYI: These (2) source files were
exported by QGIS as Dick specified above; note they are now properly formatted for use with SDK Resample
FILENAME=[path]\USGS_13_n38w090_EPSG-4326.tif
DESCRIPTION=USGS_13_n38w090_EPSG-4326.tif
UPPER LEFT X=-90.0005092597
UPPER LEFT Y=38.0005092597
LOWER RIGHT X=-88.9994907413
LOWER RIGHT Y=36.9994907413
WEST LONGITUDE=90° 00' 01.8333" W
NORTH LATITUDE=38° 00' 01.8333" N
EAST LONGITUDE=88° 59' 58.1667" W
SOUTH LATITUDE=36° 59' 58.1667" N
UL CORNER LONGITUDE=90° 00' 01.8333" W
UL CORNER LATITUDE=38° 00' 01.8333" N
UR CORNER LONGITUDE=88° 59' 58.1667" W
UR CORNER LATITUDE=38° 00' 01.8333" N
LR CORNER LONGITUDE=88° 59' 58.1667" W
LR CORNER LATITUDE=36° 59' 58.1667" N
LL CORNER LONGITUDE=90° 00' 01.8333" W
LL CORNER LATITUDE=36° 59' 58.1667" N
PROJ_DESC=Geographic (Latitude/Longitude) / WGS84 / arc degrees
PROJ_DATUM=WGS84
PROJ_UNITS=arc degrees
EPSG_CODE=EPSG:4326
COVERED AREA=9837.7 sq km
GDAL_NO_DATA_VALUE=-999999
NUM COLUMNS=10812
NUM ROWS=10812
NUM BANDS=1
PIXEL WIDTH=0.0000926 arc degrees
PIXEL HEIGHT=0.0000926 arc degrees
MIN ELEVATION=-2.5 meters
MAX ELEVATION=313.873 meters
ELEVATION UNITS=meters
BIT DEPTH=32
SAMPLE TYPE=32-bit Floating Point
GEOG_CITATION=WGS 84
PHOTOMETRIC=Greyscale (Min is Black)
BIT_DEPTH=32
SAMPLE_FORMAT=Floating Point
ROWS_PER_STRIP=1
COMPRESSION=None
PIXEL_SCALE=( 0.00009259259258, 0.00009259259258, 1 )
TIEPOINTS=( 0.00, 0.00, 0.00 ) --> ( -90.0005555560, 38.0005555560, 0.0000000000 )
MODEL_TYPE=Geographic lat-long system
RASTER_TYPE=Pixel is Area
VERT_DATUM=None Specified
FILENAME=[path]\USGS_13_n38w091_EPSG-4326.tif
DESCRIPTION=USGS_13_n38w091_EPSG-4326.tif
UPPER LEFT X=-91.0005092597
UPPER LEFT Y=38.0005092597
LOWER RIGHT X=-89.9994907363
LOWER RIGHT Y=36.9994907363
WEST LONGITUDE=91° 00' 01.8333" W
NORTH LATITUDE=38° 00' 01.8333" N
EAST LONGITUDE=89° 59' 58.1667" W
SOUTH LATITUDE=36° 59' 58.1667" N
UL CORNER LONGITUDE=91° 00' 01.8333" W
UL CORNER LATITUDE=38° 00' 01.8333" N
UR CORNER LONGITUDE=89° 59' 58.1667" W
UR CORNER LATITUDE=38° 00' 01.8333" N
LR CORNER LONGITUDE=89° 59' 58.1667" W
LR CORNER LATITUDE=36° 59' 58.1667" N
LL CORNER LONGITUDE=91° 00' 01.8333" W
LL CORNER LATITUDE=36° 59' 58.1667" N
PROJ_DESC=Geographic (Latitude/Longitude) / WGS84 / arc degrees
PROJ_DATUM=WGS84
PROJ_UNITS=arc degrees
EPSG_CODE=EPSG:4326
COVERED AREA=9837.7 sq km
GDAL_NO_DATA_VALUE=-999999
NUM COLUMNS=10812
NUM ROWS=10812
NUM BANDS=1
PIXEL WIDTH=0.0000926 arc degrees
PIXEL HEIGHT=0.0000926 arc degrees
MIN ELEVATION=92.463 meters
MAX ELEVATION=539.919 meters
ELEVATION UNITS=meters
BIT DEPTH=32
SAMPLE TYPE=32-bit Floating Point
GEOG_CITATION=WGS 84
PHOTOMETRIC=Greyscale (Min is Black)
BIT_DEPTH=32
SAMPLE_FORMAT=Floating Point
ROWS_PER_STRIP=1
COMPRESSION=None
PIXEL_SCALE=( 0.00009259259304, 0.00009259259304, 1 )
TIEPOINTS=( 0.00, 0.00, 0.00 ) --> ( -91.0005555560, 38.0005555560, 0.0000000000 )
MODEL_TYPE=Geographic lat-long system
RASTER_TYPE=Pixel is Area
VERT_DATUM=None Specified
Color Hill Shader View in Global Mapper:
Above cited 1/3 Arc Second USGS DEMs draped with ESRI ArcGIS World Imagery in Global Mapper - 3D Viewer:
0T3
Hope this helps.
GaryGB