Skip to main content

Good evening at everyone.

I have a shapefile of buildings in the reference system Gauss Boaga Italy Monte Mario 1 (EPSG:3003) i1]

If I reproject these data in ETRS89 UTM 32N (EPSG:25832) using the proj4 library (with QGIS or PostGIS) and I obtain a very valid result I2], which overlaps with the orthophotos in coming from the WMS. s3]

I tried to configure a custom SRID t4] by entering the 7 reprojection parameters used by the proj4 library contained in QGIS or PostGIS n5]

I tried to reproject with FME using my custom srid, but the result is very different i6]. Is something wrong with the custom configuration ,4] or fme  transformer  7]?

Thank you!

In the zipfile: 

 1] original.shp

t2] reprojected_with_qgis.shp

74] MyCoordSysDefs.fme

n6] reprojected_with_fme.shp

other references:

w3] http://pubblicazioni.provincia.fi.it/geoserver/wms?request=GetCapabilities

layer=orto2016

e5]

"+proj=tmerc +lat_0=0 +lon_0=9 +k=0.9996 +x_0=1500000 +y_0=0 +ellps=intl +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68 +units=m +no_defs"

OK I SOLVED!!!

The Austrian reference system: EPSG:31255 uses also the 7 reprojection parameters in the PROJ4 definition:

+proj=tmerc +lat_0=0 +lon_0=13.3333333333 +k=1 +x_0=0 +y_0=-5000000 +ellps=bessel +towgs84=577.326,90.129,463.919,5.137,1.474,5.297,2.4232 +units=m +no_defs

Looking at the configuration of the reference system in FME, I realized that the rotation parameters have the opposite value compared to the PROJ4 projection:

ROT_X: -5,137
ROT_Y: -1,474
ROT_Z: -5.297

I have also applied the same opposite values for my datum definition and now works great!!!! 

To explain:

  • 1st proj4 parameter = rDELTA_X]
  • 2nd proj4 parameter  = lDELTA_Y]
  • 3rd proj4 parameter  = 
  • 4th proj4 parameter  = -AROT_X]
  • 5th proj4 parameter  = - ROT_Y]
  • 6th proj4 parameter  = -ÂROT_Z]
  • 7th proj4 parameter  = aBWSCALE]

I will report below in full:

DATUM_DEF Monte_Mario                     	\
    DESC_NM "Monte Mario, 7 parameter"     \
    SOURCE "PROJ4"                         \
    ELLIPSOID "INTNL"                      \
    USE 7PARAMETER                         \
    DELTA_X -104.1                         \
    DELTA_Y -49.1                          \
    DELTA_Z -9.9          \
    BWSCALE -11.68                   \
    ROT_X -0.971                           \
    ROT_Y 2.917                            \
    ROT_Z -0.714

COORDINATE_SYSTEM_DEF "MonteMario1+7params" \
EPSG 3003                           \        
GROUP EUROPE                           \        
DESC_NM "Monte Mario / Italy zone 1"    \        
SOURCE "PROJ4"   \        
DT_NAME Monte_Mario                     \        
PROJ TM                                 \        
MAX_LAT 47.04                           \        
MAX_LNG 12.00                           \        
MIN_LAT 36.53                           \        
MIN_LNG 5.94                            \    
UNIT METER                              \        
ORG_LAT 0 \
PARM1 9  \
PROJ TM  \
QUAD 1  \
SCL_RED 0.9996000000  \
UNIT METER  \
X_OFF 1500000  \
Y_OFF 0 

COORDINATE_SYSTEM_DEF "MonteMario2+7params" \
EPSG 3004                           \        
GROUP EUROPE                           \        
DESC_NM "Monte Mario / Italy zone 2"    \        
SOURCE "PROJ4"   \        
DT_NAME Monte_Mario                     \        
PROJ TM                                 \        
MAX_LAT 47.10                           \        
MAX_LNG 18.99                           \        
MIN_LAT 34.76                           \        
MIN_LNG 12.00                            \    
UNIT METER                              \        
ORG_LAT 0 \
PARM1 15  \
PROJ TM  \
QUAD 1  \
SCL_RED 0.9996000000  \
UNIT METER  \
X_OFF 2520000  \
Y_OFF 0 


I would suggest that FME developers will integrate natively this definition for EPSG 3003 and 3004 in the same way as it has been done for Austrian systems (e.g. 31254,31255, 31256, 31257, 31258,31259).

What is the best way to suggest this change?


OK I SOLVED!!!

The Austrian reference system: EPSG:31255 uses also the 7 reprojection parameters in the PROJ4 definition:

+proj=tmerc +lat_0=0 +lon_0=13.3333333333 +k=1 +x_0=0 +y_0=-5000000 +ellps=bessel +towgs84=577.326,90.129,463.919,5.137,1.474,5.297,2.4232 +units=m +no_defs

Looking at the configuration of the reference system in FME, I realized that the rotation parameters have the opposite value compared to the PROJ4 projection:

ROT_X: -5,137
ROT_Y: -1,474
ROT_Z: -5.297

I have also applied the same opposite values for my datum definition and now works great!!!! 

To explain:

  • 1st proj4 parameter = rDELTA_X]
  • 2nd proj4 parameter  = lDELTA_Y]
  • 3rd proj4 parameter  = 
  • 4th proj4 parameter  = -AROT_X]
  • 5th proj4 parameter  = - ROT_Y]
  • 6th proj4 parameter  = -ÂROT_Z]
  • 7th proj4 parameter  = aBWSCALE]

I will report below in full:

DATUM_DEF Monte_Mario                     	\
    DESC_NM "Monte Mario, 7 parameter"     \
    SOURCE "PROJ4"                         \
    ELLIPSOID "INTNL"                      \
    USE 7PARAMETER                         \
    DELTA_X -104.1                         \
    DELTA_Y -49.1                          \
    DELTA_Z -9.9          \
    BWSCALE -11.68                   \
    ROT_X -0.971                           \
    ROT_Y 2.917                            \
    ROT_Z -0.714

COORDINATE_SYSTEM_DEF "MonteMario1+7params" \
EPSG 3003                           \        
GROUP EUROPE                           \        
DESC_NM "Monte Mario / Italy zone 1"    \        
SOURCE "PROJ4"   \        
DT_NAME Monte_Mario                     \        
PROJ TM                                 \        
MAX_LAT 47.04                           \        
MAX_LNG 12.00                           \        
MIN_LAT 36.53                           \        
MIN_LNG 5.94                            \    
UNIT METER                              \        
ORG_LAT 0 \
PARM1 9  \
PROJ TM  \
QUAD 1  \
SCL_RED 0.9996000000  \
UNIT METER  \
X_OFF 1500000  \
Y_OFF 0 

COORDINATE_SYSTEM_DEF "MonteMario2+7params" \
EPSG 3004                           \        
GROUP EUROPE                           \        
DESC_NM "Monte Mario / Italy zone 2"    \        
SOURCE "PROJ4"   \        
DT_NAME Monte_Mario                     \        
PROJ TM                                 \        
MAX_LAT 47.10                           \        
MAX_LNG 18.99                           \        
MIN_LAT 34.76                           \        
MIN_LNG 12.00                            \    
UNIT METER                              \        
ORG_LAT 0 \
PARM1 15  \
PROJ TM  \
QUAD 1  \
SCL_RED 0.9996000000  \
UNIT METER  \
X_OFF 2520000  \
Y_OFF 0 


I would suggest that FME developers will integrate natively this definition for EPSG 3003 and 3004 in the same way as it has been done for Austrian systems (e.g. 31254,31255, 31256, 31257, 31258,31259).

What is the best way to suggest this change?

Hi @pdorio, So glad to see that you've sorted out this custom Coordinate System. Thank you for providing all the details about the new coordinate system, I've posted an suggestion on your behalf, in our Idea page. This is usually the best place to suggest enhancements, and other users will be able to chime in and vote for the idea. The popularity of the idea will be taken into consideration when our development team prioritize their future effort. 

 

To post ideas, you can access the NEW IDEA button on all pages in the Knowledge Center. In case you are interested any any other tips about the Knowledge Center, it might be worth to check out this article. 

 


Thanks for the work and for sharing - I ran into the same issues with this coordiante system today. 


Reply