[GRASSLIST:10246] exporting to postgis - reference system problem

ivan marchesini marchesini at unipg.it
Fri Feb 10 08:30:33 EST 2006


Dear users,
I have a problem exporting maps to postgis db.
I have a map with GaussBoaga Coordinate System (an Italian Projected
System)...
I would want to export this map to a postgres/postgis database.
I have tried with 

v.out.ogr input=rupi_rec type=line 'dsn=PG:host=localhost dbname=pieve
user=brunori' olayer=rupi_def3 layer=1 format=PostgreSQL

in the geometry_columns table I can see that the srid is 32767 that
contain this srtext:
_____________________________________________________________________
PROJCS["TransverseMercator",GEOGCS["international",DATUM["Monte_Mario",
SPHEROID["international",6378388,297],TOWGS84[-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68]],
PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",15],PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",2520000],PARAMETER["false_northing",0],UNIT["meter",1]] 
_____________________________________________________________________


that seems correct .. 
I have tried visualize the postgis layer (with srtext=32767) using qgis
but it doesn't recognise the reference system.



but I know that one correct value for the srid is surely the 3004 for my
system. This is the corresponding srtext as defined in the
geometry_columns table:
_________________________________________________________________
PROJCS["Monte Mario / Italy zone
2",GEOGCS["MonteMario",DATUM["Monte_Mario",
SPHEROID["International
1924",6378388,297,AUTHORITY["EPSG","7022"]],AUTHORITY["EPSG","6265"]],
PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,
AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4265"]],PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",15],PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",2520000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AUTHORITY["EPSG","3004"]]  
_________________________________________________________________


and this is the corresponding proj4text:

__________________________________________________________________
+proj=tmerc +lat_0=0 +lon_0=15 +k=0.999600 +x_0=2520000 +y_0=0
+ellps=intl +units=m +no_defs
__________________________________________________________________




if I export from grass gis to shp file using the -e option I obtain the
next srtext in the .prj file:
 _________________________________________________________________
PROJCS["TransverseMercator",GEOGCS["international",DATUM["D_Monte_Mario",
SPHEROID["international",6378388,297]],
PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",15],PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",2520000],PARAMETER["false_northing",0],UNIT["Meter",1]]
__________________________________________________________________



then, importing from Spit Plugin in Qgis and defining the srid=3004 all
is working well!!!


the problem is surely my fault but I can't understand where I'm wrong!!!
thank you very much!!!!
Ivan





-- 
Ivan Marchesini
Department of Civil and Environmental Engineering
University of Perugia
Via G. Duranti 93/a 
06125
Perugia (Italy)
e-mail: marchesini at unipg.it
        ivan.marchesini at gmail.com
tel: +39(0)755853760
fax: +39(0)755853756







More information about the grass-user mailing list