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-boun...@postgis.refractions.net [mailto:postgis-users-boun...@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-boun...@postgis.refractions.net [mailto:postgis-users-boun...@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
