Friday 23 August 2019

gdal - Reprojection problem with gdalwarp


I have a GeoTiff file with WGS84 datum projection and I'm trying to reproject it using gdalwarp. Using this command


gdalwarp -t_srs '+proj=tpers +azi=19 +tilt=30 +lon_0=-109.5 +lat_0=23 +h=1274199.4 +wktext' atnight.tif atnighttilted.tif

returns this


ERROR 1: Too many points (10000 out of 10000) failed to transform,
unable to compute output bounds.

I'm using GDAL version 1.10.0. This is the contents of the .tfw file:



0.0066667901
0.0000000000
0.0000000000
-0.0066669136
-32.8672753199
72.8560317049

I've just started using gdalwarp. Thanks in advance, Tom


The output of atnight.tif using gdalinfo is:


Driver: GTiff/GeoTIFF

Files: atnight.tif
atnight.tfw
Size is 9123, 4957
Coordinate System is:
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"]]
Origin = (-32.870608714976186,72.859365161672642)
Pixel Size = (0.006666790125743,-0.006666913589395)
Metadata:
AREA_OR_POINT=Area
TIFFTAG_DATETIME=2012:11:26 15:13:57
TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
TIFFTAG_SOFTWARE=Adobe Photoshop CS6 (Macintosh)
TIFFTAG_XRESOLUTION=72

TIFFTAG_YRESOLUTION=72
Image Structure Metadata:
COMPRESSION=LZW
INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left ( -32.8706087, 72.8593652) ( 32d52'14.19"W, 72d51'33.71"N)
Lower Left ( -32.8706087, 39.8114745) ( 32d52'14.19"W, 39d48'41.31"N)
Upper Right ( 27.9505176, 72.8593652) ( 27d57' 1.86"E, 72d51'33.71"N)
Lower Right ( 27.9505176, 39.8114745) ( 27d57' 1.86"E, 39d48'41.31"N)
Center ( -2.4600456, 56.3354198) ( 2d27'36.16"W, 56d20' 7.51"N)

Band 1 Block=9123x115 Type=Byte, ColorInterp=Red
Band 2 Block=9123x115 Type=Byte, ColorInterp=Green
Band 3 Block=9123x115 Type=Byte, ColorInterp=Blue


No comments:

Post a Comment

arcpy - Changing output name when exporting data driven pages to JPG?

Is there a way to save the output JPG, changing the output file name to the page name, instead of page number? I mean changing the script fo...