添加“使用GDAL获取EPSG编码”

13693261870

2023-11-29 a3d8b2a5a449f2f633973786573b5306f29a517e

 src/main/java/com/lf/server/service/data/RasterService.java

@@ -1,11 +1,18 @@
package com.lf.server.service.data;

import com.lf.server.entity.all.StaticData;
import com.lf.server.helper.StringHelper;
import com.lf.server.helper.WebHelper;
import org.apache.commons.logging.Log;
import org.apache.commons.logging.LogFactory;
import org.gdal.gdal.Dataset;
import org.gdal.gdal.gdal;
import org.gdal.gdalconst.gdalconst;
import org.gdal.ogr.Geometry;
import org.gdal.ogr.ogr;
import org.gdal.osr.SpatialReference;
import org.osgeo.proj4j.*;
import org.springframework.beans.factory.annotation.Value;
import org.springframework.stereotype.Service;

import java.io.File;
@@ -17,7 +24,28 @@
 */
@Service
public class RasterService {
    private SpatialReference sr4326;

    private SpatialReference sr4490;

    @Value("${sys.gdal_path}")
    private String gdalPath;

    private final static Log log = LogFactory.getLog(RasterService.class);

    /**
     * 初始化空间参考
     */
    public void initSr() {
        if (null == sr4326) {
            sr4326 = new SpatialReference();
            sr4326.ImportFromEPSG(StaticData.I4326);
        }
        if (null == sr4490) {
            sr4490 = new SpatialReference();
            sr4490.ImportFromEPSG(StaticData.I4490);
        }
    }

    /**
     * 获取栅格数据的EPSG编码
@@ -55,4 +83,77 @@
            }
        }
    }

    /**
     * 使用GDAL获取EPSG编码
     */
    public Integer getEpsgByGdal(String file) {
        try {
            File f = new File(file);
            if (!f.exists() || f.isDirectory()) {
                return null;
            }

            String cmd = String.format("%s\\gdalsrsinfo.exe \"%s\" -o epsg", gdalPath, file);
            String rs = WebHelper.exec2(cmd);
            if (StringHelper.isEmpty(rs)) {
                return null;
            }

            String[] strs = rs.split("\n");
            for (String str : strs) {
                if (str.contains("EPSG:")) {
                    String epsg = str.replace("EPSG:", "");

                    return Integer.parseInt(epsg);
                }
            }

            return 0;
        } catch (Exception ex) {
            log.error(ex.getMessage(), ex);
            return null;
        }
    }

    /**
     * 坐标转换
     */
    public Object csTransform(double x, double y, int epsg) {
        this.initSr();

        SpatialReference srs = new SpatialReference();
        srs.ImportFromEPSG(epsg);

        Geometry point = new Geometry(ogr.wkbPoint);
        point.AddPoint(x, y, 0);
        point.AssignSpatialReference(srs);

        if (srs.GetName().contains(StaticData.CGCS2000)) {
            point.TransformTo(sr4490);
        } else {
            point.TransformTo(sr4326);
        }
        point.SwapXY();

        return new double[]{point.GetX(), point.GetY()};
    }

    /**
     * 坐标转换-使用Proj4j库
     */
    public String transformByProj4j(double x, double y, int epsg) {
        CRSFactory crsFactory = new CRSFactory();
        CoordinateReferenceSystem fromCrs = crsFactory.createFromName("EPSG:" + epsg);
        CoordinateReferenceSystem toCrs = crsFactory.createFromName("EPSG:4326");

        CoordinateTransformFactory ctf = new CoordinateTransformFactory();
        CoordinateTransform transform = ctf.createTransform(fromCrs, toCrs);

        ProjCoordinate fromCoord = new ProjCoordinate(x, y);
        ProjCoordinate toCoord = new ProjCoordinate();
        transform.transform(fromCoord, toCoord);

        return String.format("%f,%f", toCoord.x, toCoord.y);
    }
}