[Gdal-dev] sterea projection

[Jan Hartmann]

Frank Warmerdam wrote:
> 
> Yes, it should.  Could you perhaps try processing these
> through "testepsg"?  If you give it a PROJ.4 coordinate
> system as input it should report various representations,
> including what it looks like after going to WKT and then
> back to PROJ.4.
> 

You are right for "testepsg":

- testepsg "+proj=sterea ..."

and

- testepsg "+proj=sterea +towgs84=..."

gives the correct output for the generated proj-strings. The towgs 
parameter is present in the second case. "testepsg -t" with a test 
coordinate gives correct different results in both cases.

However,

- gdalwarp -s_srs "+proj=sterea..." -t_srs "+proj=longlat..."

and

- gdalwarp -s_srs "+proj=sterea +towgs84=..." -t_srs "+proj=longlat..."

results in two identical output files, both without the towgs84 transform.

Jan