Hello community,
I'm trying to do a coordinate transformation from the Austrian EPSG:31256
to EPSG:4326 with the following code in a h2gis database:
*SELECT* ST_Transform(ST_GeomFromText('POINT(-38048.66 389405.66)', 31256),
4326) *FROM* dual;
This gives the following error message
"ST_Transform(conn48: url=jdbc:default:connection user=, POINT (-38048.66
389405.66), 4326): No transformation found from epsg:31256 to epsg:4326"
Exception calling user-defined function: "ST_Transform(conn48:
url=jdbc:default:connection user=, POINT (-38048.66 389405.66), 4326): No
transformation found from epsg:31256 to epsg:4326"; SQL statement:
SELECT ST_Transform(ST_GeomFromText('POINT(-38048.66 389405.66)', 31256),
4326) FROM dual [90105-197]
Other operations like ST_Within do work (if both geometries have the same
SRID).
Doing the same conversion in Java code using just the libraries imported by
the Maven artifact org.orbisgis:h2gis-functions:1.3.2 works perfectly fine:
import java.util.List;
import org.cts.CRSFactory;
import org.cts.IllegalCoordinateException;
import org.cts.crs.CRSException;
import org.cts.crs.CoordinateReferenceSystem;
import org.cts.crs.GeodeticCRS;
import org.cts.op.CoordinateOperation;
import org.cts.op.CoordinateOperationFactory;
import org.cts.registry.EPSGRegistry;
import org.junit.jupiter.api.Test;
import io.quarkus.test.junit.QuarkusTest;
// test code based on
//
https://github.com/orbisgis/cts/wiki/1.-Create-a-new-CoordinateReferenceSystem-from-a-reference-code
// and
https://geoinformatik.htw-dresden.de/anleitungen/gi-app/JavaKoordinatentransformationen.pdf
@QuarkusTest
public class CoordinatesTransformationTest {
@Test
void testConvertAustriaGkEastToGPS() throws CRSException,
IllegalCoordinateException {
var crsFactory = new CRSFactory();
crsFactory.getRegistryManager().addRegistry(new EPSGRegistry());
CoordinateReferenceSystem gps = crsFactory.getCRS("EPSG:4326");
CoordinateReferenceSystem austria = crsFactory.getCRS("EPSG:31256");
List<CoordinateOperation> ops = CoordinateOperationFactory
.createCoordinateOperations((GeodeticCRS) austria, (GeodeticCRS) gps);
double[] coord = new double[2];
coord[0] = -38048.66; // Longitude or Easting
coord[1] = 389405.66; // Latitude or Northing
for (var op : ops) {
var transformedCoord = op.transform(coord);
System.out.println("Longitude: " + transformedCoord[0] + " Latitude: " +
transformedCoord[1]);
}
}
}
It outputs
Longitude: 15.815805676616735 Latitude: 48.64153355841612
What am I doing wrong in the SQL version?
Any help is appreciated. Thank you very much!
best regards,
Christoph
--
You received this message because you are subscribed to the Google Groups "H2
Database" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to h2-database+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/h2-database/104321de-d151-44a1-8ba8-9849d375b147n%40googlegroups.com.
