Hi all,
I am using geotools to convert coordinates for seismic data. I've never
had any issue before but this week I had to convert data using EPSG
22192 and got wrong results (my data is supposed to be onshore Argentina
and ended up close to the antartic).
Geotools does the following conversion (EPSG 22192 -> EPSG 4326)
x = 2572506.0 -> Lat = -57.53722563
y = 5757839.0 -> Lon = -38.73009392
The correct conversion (according to epsg.io
http://epsg.io/transform#s_srs=22192&t_srs=4326&x=2572506.0000000&y=5757839.0000000 should
be)
Transform coordinates - GPS online convertor
<http://epsg.io/transform#s_srs=22192&t_srs=4326&x=2572506.0000000&y=5757839.0000000>
Transform your coordinates online easily with epsg.io Transform
coordinates Online convertor for lat & long coordinates, geodetic datums
and projected systems
epsg.io
Lat = -38.328903
Lon = -68.1718113
After digging through the source code I found that WKT are defined (at
least) under
https://raw.githubusercontent.com/geotools/geotools/master/modules/plugin/epsg-wkt/src/main/resources/org/geotools/referencing/epsg/wkt/epsg.properties
I then compared it to the definition on https://epsg.io/22192 and found
different values.
Definition from Geotools:
PROJCS["Campo Inchauspe / Argentina 2",
GEOGCS["Campo Inchauspe",
DATUM["Campo Inchauspe",
SPHEROID["International 1924", 6378388.0, 297.0,
AUTHORITY["EPSG","7022"]],
TOWGS84[-154.5,150.7,100.4,0.0,0.0, 0.0,0.0],
AUTHORITY["EPSG","6221"]],
PRIMEM["Greenwich",0.0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.017453292519943295],
AXIS["Geodetic longitude",EAST],
AXIS["Geodetic latitude",NORTH],
AUTHORITY["EPSG","4221"]],
PROJECTION["Transverse_Mercator",
AUTHORITY["EPSG","9807"]],
PARAMETER["central_meridian", -69.0],
PARAMETER["latitude_of_origin",-90.0],
PARAMETER["scale_factor",1.0],
PARAMETER["false_easting",2500000.0],
PARAMETER["false_northing",0.0],
UNIT["m",1.0],
AXIS["Easting",EAST],
AXIS["Northing",NORTH],
AUTHORITY["EPSG","22192"]]
Definition from epsg.io:
PROJCS["Campo Inchauspe / Argentina 2",
GEOGCS["Campo Inchauspe",
DATUM["Campo_Inchauspe",
SPHEROID["International 1924",6378388,297,
AUTHORITY["EPSG","7022"]],
TOWGS84[-148,136,90,0,0,0,0],
AUTHORITY["EPSG","6221"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4221"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",-90],
PARAMETER["central_meridian",-69],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",2500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","22192"]]
I am not an expert in coordinate system. Would this difference be the
cause of my issue? If so, how can I fix this (scala code that depends on
geospark, a package that uses geotools)?
Is there anything in my code below (scala) that is wrong otherwise?
val sourceCRS = CRS.decode(s"EPSG:22192")
val targetCRS = CRS.decode(s"EPSG:4326")
val transform = CRS.findMathTransform(sourceCRS, targetCRS, true)
val cdp_x =2572506.0
val cdp_y = 5757839.0
val ptSrc = new DirectPosition2D(sourceCRS, cdp_x, cdp_y)
val ptDst = new DirectPosition2D()
transform.transform(ptSrc, ptDst)
println(ptDst.getCoordinate()(0))
println(ptDst.getCoordinate()(1))
Campo Inchauspe / Argentina 2 - EPSG:22192 <https://epsg.io/22192>
EPSG:22192 Projected coordinate system for Argentina - between 70°30'W
and 67°30'W, onshore. Large and medium scale topographic mapping and
engineering survey.
epsg.io
Thanks,
Julien Bluteau
_______________________________________________
GeoTools-GT2-Users mailing list
GeoTools-GT2-Users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/geotools-gt2-users
