[gdal-dev] Re: gdal to define projection, EPSG:3412

[Frank Warmerdam]

Hermann Peifer wrote:
> John Callahan wrote:
>> Thanks for responding Hermann. Yes, I tried using -s_srs EPSG:3412 as 
>> well as the full PROJ4 string with the same result. When specifying 
>> -s_srs, does gdalwarp ignore (override) the project information in the 
>> hdr file?
>>
> 
> Hmm. AFAIK: gdal_translate -a_srs does override the file's projection 
> information. I have a random.tif in LAEA projection (EPSG:3035), so I 
> tried the following with gdal 1.6.2:
> 
> gdal_translate random.tif out1.tif -a_srs EPSG:3412  # this results in [1]
> 
> gdal_translate random.tif out2.tif -a_srs "+proj=stere +lat_0=-90 
> +lat_ts=-70 +lon_0=0 +k=1 +x_0=0 +y_0=0 +a=6378273 +b=6356889.449 
> +units=m +no_defs"  # this gives [2]
> 
> To me as an innocent GDAL user, this behaviour looks similar to 
> http://trac.osgeo.org/gdal/ticket/3016. I might be wrong, though, in 
> particular as you wrote that using the PROJ.4 string didn't help (and 
> your source file isn't a GeoTIFF in the first place).

Hermann,

I have confirmed this problem and filed it as:

   http://trac.osgeo.org/gdal/ticket/3220

It is presumably a problem with translation of a specific projection method,
possibly in libgeotiff itself.

Best regards,
-- 
---------------------------------------+--------------------------------------
I set the clouds in motion - turn up   | Frank Warmerdam, warmerdam at pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush    | Geospatial Programmer for Rent