[prev in list] [next in list] [prev in thread] [next in thread] 

List:       postgis-users
Subject:    Re: [postgis-users] Coordinate transformation with custom recordsofspatial_sys_ref
From:       DUPREZ_Cédric <Cedric.DUPREZ () ign ! fr>
Date:       2012-06-29 14:31:47
Message-ID: 5EB13D1B2A28494FBA69E90063F4ED84741C75 () populus3 ! ifn ! fr
[Download RAW message or body]

Thank you for the answer.

I will test with AddGeometryColumm.

Unfortunately, the coordinate transformation between 27582 and 2154 SRID gives wrong \
results, since it does not consider the correction grid specified by the French IGN. \
The custom SRIDs does.

Regards,

Cedric

-----Message d'origine-----
De : postgis-users-bounces@postgis.refractions.net \
[mailto:postgis-users-bounces@postgis.refractions.net] De la part de Francois Hugues \
Envoyé : vendredi 29 juin 2012 16:21 À : PostGIS Users Discussion
Objet : Re: [postgis-users] Coordinate transformation with custom \
recordsofspatial_sys_ref


Hello,

I think you use wrong srid : for lambert93 you should use 2154 and for lambert 2 \
carto it is 27582. 

Anyway, if you want to use the SRID 320002120 and 310024140, your problem is related \
to the geometry_column view and constraints checking on geometry columns. You will \
have to use AddGeometryColumm with use_typmod option set to false.

http://postgis.refractions.net/docs/AddGeometryColumn.html said : "Changed: 2.0.0 If \
you require the old behavior of constraints use the default use_typmod, but set it to \
false."

Hugues.

-----Message d'origine-----
De : postgis-users-bounces@postgis.refractions.net \
[mailto:postgis-users-bounces@postgis.refractions.net] De la part de DUPREZ Cédric \
Envoyé : vendredi 29 juin 2012 15:46 À : postgis-users@postgis.refractions.net
Objet : [postgis-users] Coordinate transformation with custom records \
ofspatial_sys_ref


Hi all list members,

I am facing a problem with coordinate transformation, using ST_Transform() function.

First of all, my PostGIS version is 2.0.1 : POSTGIS="2.0.1 r9979" \
GEOS="3.3.5-CAPI-1.7.5" PROJ="Rel. 4.8.0, 6 March 2012" GDAL="GDAL 1.9.1, released \
2012/05/15" LIBXML="2.7.6" LIBJSON="UNKNOWN" TOPOLOGY RASTER

Here are the steps to reproduce the problem:
1/ Create a table containing points
CREATE TABLE pts (
	code INT,
	xl2 INT,
	yl2 INT,
	geom GEOMETRY,
	CONSTRAINT pts_pkey PRIMARY KEY (code)
)
WITHOUT oids;

INSERT INTO pts (code, xl2, yl2)
VALUES (1, 556866, 1902345)
, (2, 689174, 2448195)
, (3, 729095, 2528577)
, (4, 335755, 2059702)
, (5, 623207, 2210468)
, (6, 850337, 1871486)
, (7, 915275, 1889534)
, (8, 639551, 1778327)
, (9, 584365, 2273523)
, (10, 649599, 2228833);

2/ Import the custom records of spatial_sys_ref specified by the French IGN. The file \
is provided here: http://lambert93.ign.fr/fileadmin/files/IGNF-spatial_ref_sys.sql

3/ Populate the geom column:
UPDATE pts
SET geom = ST_SetSRID(ST_MakePoint(xl2, yl2), 320002120);
I get the following notices: NOTICE:  SRID value 320002120 > SRID_MAXIMUM converted \
to 999442

SELECT DISTINCT ST_SRID(geom)
FROM pts;
Gives: 999442

But SELECT * FROM geometry_columns;
Gives: 0

4/ Update geometry SRID:
SELECT UpdateGeometrySRID('pts', 'geom', 320002120);

SELECT * FROM geometry_columns;
Gives: 999442

5/ Try the transformation from 320002120 SRID (Lambert II) to 310024140 SRID (Lambert \
93): SELECT code, ST_Transform(geom, 310024140)
FROM pts;
I get the following error message: Invalid reserved SRID (999442).

With PostGIS 1.5.3, this worked perfectly, without any error message and gave the \
good coordinate transformation.

What should I change with PostGIS 2.0.1 in order to get the same result again?

Thanks in advance for your help.

Regards,

Cedric Duprez
_______________________________________________
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users
_______________________________________________
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users
_______________________________________________
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


[prev in list] [next in list] [prev in thread] [next in thread]