[gdal-dev] Issue reprojecting geo tiff

[Even Rouault]

Le mercredi 14 août 2013 14:24:56, Reaves, Timothy a écrit :
> I'm using GDAL 1.10.0, released 2013/04/24 .  I am playing with TileMill,
> and want to add USGS topo maps. These topo maps are geo PDF.  I'm following
> along one of their documents, and the first thing I do is convert the PDF
> to TIFF.
> 
> gdal_translate -of GTiff -co "TILED=YES" -co "TFW=YES"
> PA_Tionesta_20130621_TM_geo.pdf PA_Tionesta_20130621_TM_geo.tiff
> ERROR 1: Couldn't find group for reference to set OFF
> Input file size is 13652, 17403
> 0...10...20...30...40...50...60...70...80...90...100 - done.
> 
> I do not like the ERROR; that should mean that the conversion failed;
> however, and output TIFF is generated.  So I use it.

You can probably ignore it, but I'm not 100% sure without investigation. This 
error seems to come from Poppler. Might be linked to code turning layers on. 
If this error has consequences, they will be layers not displayed in the PDF.

> 
> The next step is the re-translation:
> 
> gdalwarp -s_srs EPSG:4269 -t_srs EPSG:3785 -r bilinear
> PA_Tionesta_20130621_TM_geo.tiff PA_Tionesta_20130621_TM_geo_3785.tiff
> ERROR 1: Too many points (10201 out of 10201) failed to transform,
> unable to compute output bounds.
> 
> This ERROR does not create output.

Just drop the "-s_srs EPSG:4269". As the GeoPDF driver recognizes the input 
projection, it is completely useless to specify it. The error here is that the 
source data is in EPSG:26717 (UTM 17 / NAD 83), not EPSG:4269 (geodetic NAD 
83).


-- 
Geospatial professional services
http://even.rouault.free.fr/services.html