[Gdal-dev] sterea projection

[Frank Warmerdam]

On 11/10/05, Jan Hartmann <j.l.h.hartmann at uva.nl> wrote:
> 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,

OK, that is weird.  Could you submit a bugzilla entry with
exact commands, and an input file that demonstrates the problem?

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