[postgis-users] transform() does not do work for NGO1948 zone III-> WGS84 longlat (epsg:27393->epsg:4326)
strk at refractions.net
strk at refractions.net
Thu Dec 22 04:22:13 PST 2005
mmm... I checked using transform_geometry() to exclude
the new proj4 caching code... It still gives the same
result:
strk=# select astext(transform('SRID=27393;POINT(2836.394208896 185526.991186176)',4326));
astext
--------------------------------------------
POINT(0.0464201635443672 59.6661916850467)
(1 row)
strk=# select astext(transform_geometry('SRID=27393;POINT(2836.394208896 185526.991186176)', get_proj4_from_srid(27393), get_proj4_from_srid(4326), 4326));
astext
--------------------------------------------
POINT(0.0464201635443672 59.6661916850467)
(1 row)
This is with Rel. 4.4.8, 3 May 2004, PostGIS 1.1.0.
PROJ4 text for 27393 is as follows:
+proj=tmerc +lat_0=58 +lon_0=-10.72291666666667 +k=1.000000 +x_0=0 +y_0=0 +a=6377492.018 +b=6356173.508712696 +towgs84=278.3,93,474.5,7.889,0.05,-6.61,6.21 +pm=oslo +units=m +no_defs
PROJ4 text for 4326 is as follows:
+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs
This is both in the db and in the spatial_ref_sys.sql script in 1.1.0.
Can you check corrispondent values in your old *working* setup ?
--strk;
On Thu, Dec 22, 2005 at 01:07:17PM +0100, Havard Tveite wrote:
> It is nice to provide warnings, but for epsg:27393 -> epsg:4326
> the following problems persist:
>
> 1) Longitude of origin (+lon_0=-10.72291666666667) is ignored
> (as demonstrated by strk's query).
> 2) No datum shift seems to be performed.
>
> Since the transformation runs smoothly using UMN Mapserver (
> with dynamically linked proj - so it should be 4.4.9), there
> is probably a problem somewhere.
>
>
> I have an old PostGIS running (0.8 over 7.4.1 - proj statically
> linked "old" version), and for that setup, the transformation
> seems to be OK (including datum shift).
>
> select AsText(
> transform('SRID=27393;POINT(2836.394208896 185526.991186176)',4326)
> );
> > POINT(10.7732463298862 59.665688397594)
>
>
> Håvard
More information about the postgis-users
mailing list