It is a reference to the file : /usr/share/proj/ntf_r93.gsb The permissions are the following : -rw-r--r--
Cedric -----Message d'origine----- De : postgis-users-boun...@postgis.refractions.net [mailto:postgis-users-boun...@postgis.refractions.net] De la part de Stephen Woodbridge Envoyé : lundi 22 octobre 2012 17:27 À : postgis-users@postgis.refractions.net Objet : Re: [postgis-users] Coordinatetransformationwithcustomrecordsofspatial_sys_ref Hi Cédric, +nadgrids=ntf_r93.gsb is a reference to the file: /usr/share/proj/ntf_r93.gsb or /usr/local/share/proj/ntf_r93.gsb Have you verified that you have this file and the permissions are something like 644? -Steve W On 10/22/2012 11:16 AM, DUPREZ Cédric wrote: > Hi all, > > After several tries, I am still unable to get this feature to work within > Postgis 2.0... > And with PostgreSQL 9.2, it is no more possible to downgrade to Postgis 1.5... > I am completely locked. > > I try, once again, to give the steps to reproduce my bug: > > 1/ insert spatial_ref_sys entries which used to work in Postgis 1.5 : > insert into spatial_ref_sys values > (320002120,'IGNF',320002120,'PROJCS["Lambert II étendu",GEOGCS["Nouvelle > Triangulation Française Paris grades",DATUM["NTF",SPHEROID["Clarke 1880 > IGN",6378249.2000,293.4660210000000,AUTHORITY["IGNF","ELG010"]],TOWGS84[-168.0000,-60.0000,320.0000,0,0,0,0],AUTHORITY["IGNF","REG002"]],PRIMEM["Paris",2.337229167,AUTHORITY["IGNF","LGO02"]],UNIT["grad",0.01570796326794897],AXIS["Longitude",EAST],AXIS["Latitude",NORTH],AUTHORITY["IGNF","NTFP"]],PROJECTION["Lambert_Conformal_Conic_1SP",AUTHORITY["IGNF","PRC0120"]],PARAMETER["semi_major",6378249.2000],PARAMETER["semi_minor",6356515.0000],PARAMETER["latitude_of_origin",46.800000000],PARAMETER["central_meridian",0.000000000],PARAMETER["scale_factor",0.99987742],PARAMETER["false_easting",600000.000],PARAMETER["false_northing",2200000.000],UNIT["metre",1],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["IGNF","LAMBE"]]','+init=IGNF:LAMBE > +title=Lambert II etendu +proj=lcc +nadgrids=ntf_r93.gsb,null +wkte xt +towgs84=-168.0000,-60.0000,320.0000 +a=6378249.2000 +rf=293.4660210000000 +pm=2.337229167 +lat_0=46.800000000 +lon_0=0.000000000 +k_0=0.99987742 +lat_1=46.800000000 +x_0=600000.000 +y_0=2200000.000 +units=m +no_defs'); > insert into spatial_ref_sys values > (310024140,'IGNF',310024140,'PROJCS["Lambert 93",GEOGCS["Réseau géodésique > français 1993",DATUM["Réseau géodésique français 1993",SPHEROID["IAG GRS > 1980",6378137.0000,298.2572221010000,AUTHORITY["IGNF","ELG037"]],TOWGS84[0.0000,0.0000,0.0000,0,0,0,0],AUTHORITY["IGNF","REG024"]],PRIMEM["Greenwich",0.000000000,AUTHORITY["IGNF","LGO01"]],UNIT["degree",0.01745329251994330],AXIS["Longitude",EAST],AXIS["Latitude",NORTH],AUTHORITY["IGNF","RGF93G"]],PROJECTION["Lambert_Conformal_Conic_2SP",AUTHORITY["IGNF","PRC0140"]],PARAMETER["semi_major",6378137.0000],PARAMETER["semi_minor",6356752.3141],PARAMETER["latitude_of_origin",46.500000000],PARAMETER["central_meridian",3.000000000],PARAMETER["standard_parallel_1",44.000000000],PARAMETER["standard_parallel_2",49.000000000],PARAMETER["false_easting",700000.000],PARAMETER["false_northing",6600000.000],UNIT["metre",1],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["IGNF","LAMB93"]]','+init=IGNF:LAMB93'); > > 2/ This raises an when adding a geometry column : "ERROR: AddGeometryColumn() > - SRID must be <= 998999" > > 3/ I reaffect a neuw SRID to those 2 entries : > insert into spatial_ref_sys values (910002,'IGNF',320002120,'PROJCS["Lambert > II étendu",GEOGCS["Nouvelle Triangulation Française Paris > grades",DATUM["NTF",SPHEROID["Clarke 1880 > IGN",6378249.2000,293.4660210000000,AUTHORITY["IGNF","ELG010"]],TOWGS84[-168.0000,-60.0000,320.0000,0,0,0,0],AUTHORITY["IGNF","REG002"]],PRIMEM["Paris",2.337229167,AUTHORITY["IGNF","LGO02"]],UNIT["grad",0.01570796326794897],AXIS["Longitude",EAST],AXIS["Latitude",NORTH],AUTHORITY["IGNF","NTFP"]],PROJECTION["Lambert_Conformal_Conic_1SP",AUTHORITY["IGNF","PRC0120"]],PARAMETER["semi_major",6378249.2000],PARAMETER["semi_minor",6356515.0000],PARAMETER["latitude_of_origin",46.800000000],PARAMETER["central_meridian",0.000000000],PARAMETER["scale_factor",0.99987742],PARAMETER["false_easting",600000.000],PARAMETER["false_northing",2200000.000],UNIT["metre",1],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["IGNF","LAMBE"]]','+init=IGNF:LAMBE > +title=Lambert II etendu +proj=lcc +nadgrids=ntf_r93.gsb,null +wktext +towgs84=-168.0000,-60.0000,320.0000 +a=6378249.2000 +rf=293.4660210000000 +pm=2.337229167 +lat_0=46.800000000 +lon_0=0.000000000 +k_0=0.99987742 +lat_1=46.800000000 +x_0=600000.000 +y_0=2200000.000 +units=m +no_defs'); > insert into spatial_ref_sys values (910001,'IGNF',310024140,'PROJCS["Lambert > 93",GEOGCS["Réseau géodésique français 1993",DATUM["Réseau géodésique > français 1993",SPHEROID["IAG GRS > 1980",6378137.0000,298.2572221010000,AUTHORITY["IGNF","ELG037"]],TOWGS84[0.0000,0.0000,0.0000,0,0,0,0],AUTHORITY["IGNF","REG024"]],PRIMEM["Greenwich",0.000000000,AUTHORITY["IGNF","LGO01"]],UNIT["degree",0.01745329251994330],AXIS["Longitude",EAST],AXIS["Latitude",NORTH],AUTHORITY["IGNF","RGF93G"]],PROJECTION["Lambert_Conformal_Conic_2SP",AUTHORITY["IGNF","PRC0140"]],PARAMETER["semi_major",6378137.0000],PARAMETER["semi_minor",6356752.3141],PARAMETER["latitude_of_origin",46.500000000],PARAMETER["central_meridian",3.000000000],PARAMETER["standard_parallel_1",44.000000000],PARAMETER["standard_parallel_2",49.000000000],PARAMETER["false_easting",700000.000],PARAMETER["false_northing",6600000.000],UNIT["metre",1],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["IGNF","LAMB93"]]','+init=IGNF:LAMB93'); > > 4/ Performing a ST_Transform fails : > select xl, yl, ST_X(geom)::INT AS xl93, ST_Y(geom)::INT AS yl93 > from ( > select *, ST_Transform(ST_SetSRID(ST_MakePoint(xl, yl), 910143), 910142) AS > geom from (values (556866, 1902345) > , (689174, 2448195) > , (729095, 2528577) > , (335755, 2059702) > , (623207, 2210468) > , (850337, 1871486) > , (915275, 1889534) > , (639551, 1778327) > , (584365, 2273523) > , (649599, 2228833)) as t(xl, yl) > )g order by xl; > Gives : " PostGIS was unable to transform the point because either no grid > shift files were found, or the point does not lie within the range for which > the grid shift is defined. Refer to the ST_Transform() section of the PostGIS > manual for details on how to configure PostGIS to alter this behaviour." > > I repeat that this conversion works good under proj4 : cs2cs -I > +init=IGNF:LAMBE +to +init=IGNF:LAMB93 > > I have tried to specify the proj4text of IGNF:LAMBE with its nadgrids, but I > get the same result. > > Any help would be appreciated. > > Thanks in advance, > > Cedric > > > > -----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é : jeudi 19 juillet 2012 09:22 > À : PostGIS Users Discussion > Objet : Re: [postgis-users] Coordinate > transformationwithcustomrecordsofspatial_sys_ref > > Thank you for your help. > I have modified the SRIDs of each entry in the IGNF file, by incrementing > them from 910001 to 910320. > > But now, when I try the conversion with ST_Transform, I get the following > error message: > > "transform: couldn't project point (556866 1.90234e+06 0): failed to load > datum shift file (-38)" > With the following notice: "PostGIS was unable to transform the point because > either no grid shift files were found, or the point does not lie within the > range for which the grid shift is defined. Refer to the ST_Transform() > section of the PostGIS manual for details on how to configure PostGIS to > alter this behaviour." > > I tried the conversion directly in Proj4 with cs2cs and it works correctly. > This should mean that Proj.4 is build with the grid-shift files, as it is > said in the ST_Transform documentation. > Here is the command I execute in Proj4 : cs2cs -I +init=IGNF:LAMBE +to > +init=IGNF:LAMB93 > Then, giving Lambert II coordinates transform them correctly. > But in Postgis, it does not work. > > What am I doing wrong? > > Regards, > > Cedric > > > > -----Message d'origine----- > De : postgis-users-boun...@postgis.refractions.net > [mailto:postgis-users-boun...@postgis.refractions.net] De la part de > robert_cl...@doh.state.fl.us > Envoyé : mardi 17 juillet 2012 14:49 > À : postgis-users@postgis.refractions.net > Objet : Re: [postgis-users] Coordinate transformation > withcustomrecordsofspatial_sys_ref > > One easy way to generate the syntax required to add a custom projection is to > select something similar at spatialreference.org and click "PostGIS > spatial_ref_sys INSERT statement," then copy and paste into a text or query > editor to revise. Blinding flash of the obvious, perhaps, but hope it's > helpful to someone. > > -Rob > > -----Original Message----- > From: postgis-users-boun...@postgis.refractions.net > [mailto:postgis-users-boun...@postgis.refractions.net] On Behalf Of Sandro > Santilli > Sent: Tuesday, July 17, 2012 5:49 AM > To: DUPREZ Cédric > Cc: PostGIS Users Discussion > Subject: Re: [postgis-users] Coordinate transformation with > customrecordsofspatial_sys_ref > > On Tue, Jul 17, 2012 at 11:34:46AM +0200, DUPREZ Cédric wrote: > >> SELECT AddGeometryColumn('public', 'pts', 'geom', 320002120, 'POINT', >> 2, false); >> >> This gives me the following error message: "ERROR: AddGeometryColumn() >> - SRID must be <= 998999 >> >> How can I use the correct 320002120 SRID specified by French IGN? > > Is 320002120 in your spatial_ref_sys already ? > You'll need to assign it a different SRID, the spatial_ref_sys has provision > for "auth_srid" so you can put 320002120 in there and use something valid for > PostGIS' "srid" field. > > Unfortunately we never got around writing a page about SRID ranges in PostGIS > manual (contributions welcome!) but the idea is that you put your "custom" > entries somewhere from 910000 and below 999000. > > --strk; > > ,------o-. > | __/ | Delivering high quality PostGIS 2.1 > | / 2.1 | http://strk.keybit.net - http://vizzuality.com > `-o------' > > _______________________________________________ > 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 > _______________________________________________ > 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
