如何用JTS对Geometry空间数据进行坐标系投影转换

蔚落 2021-09-22 02:34 1273阅读 0赞

如何利用jts对空间数据进行坐标系投影转换，以EPSG:4326转EPSG:3857为例，其他转换调换参数即可。

import org.geotools.geometry.jts.JTS;
    import org.geotools.referencing.CRS;
    import org.locationtech.jts.geom.Envelope;
    import org.locationtech.jts.geom.Geometry;
    import org.opengis.referencing.FactoryException;
    import org.opengis.referencing.crs.CoordinateReferenceSystem;
    import org.opengis.referencing.operation.MathTransform;
    
    public class Reproject {
        private CoordinateReferenceSystem sourceCRS;
        private CoordinateReferenceSystem targetCRS;
    
        public Reproject(String srcCrsCode, String tarCrsCode, boolean longFirst) {
            try{
                sourceCRS = CRS.decode(srcCrsCode, longFirst);
                targetCRS = CRS.decode(tarCrsCode, longFirst);
            }catch (FactoryException e){
                System.out.println(e);
            }
        }
    
        public Geometry reproject(Geometry geometry){
            Geometry reprojectGeometry = null;
            try {
                MathTransform transform = CRS.findMathTransform(sourceCRS, targetCRS);
                reprojectGeometry = JTS.transform(geometry, transform);
            }catch (Exception e){
                System.out.println(e);
            }
            return reprojectGeometry;
        }
    
        public Envelope reproject(Envelope envelope){
            Envelope reprojectEnvelope = null;
            try {
                MathTransform transform = CRS.findMathTransform(sourceCRS, targetCRS);
                reprojectEnvelope = JTS.transform(envelope, transform);
            }catch (Exception e){
                System.out.println(e);
            }
            return reprojectEnvelope;
        }
    
        public static void main(String[] args) {
            String srcCrsCode = "EPSG:4326";
            String tarCrsCode = "EPSG:3857";
            Reproject reproject = new Reproject(srcCrsCode,tarCrsCode,true);
            Geometry geometry = reproject.reproject(geo);
        }
        
    }