修改栅格读取服务

13693261870

2023-08-25 dfded458c5ee662a2dc4a51a6d0640e1a4f4a4bb

修改栅格读取服务

 src/main/java/com/moon/server/entity/all/StaticData.java

@@ -76,34 +76,16 @@

    public static final int I104903 = 104903;

    /**
     * 字符1
     */
    public final static String S1 = "1";

    /**
     * 等号
     */
    public final static String EQ = "=";

    /**
     * 字符点
     */
    public final static String POINT = ".";

    /**
     * 逗号
     */
    public final static String COMMA = ",";

    /**
     * 波浪号
     */
    public final static String TILDE = "~";

    /**
     * 单引号
     */
    public final static String SINGLE_QUOTES = "'";

    public final static String BBOREHOLE = "bborehole";
@@ -114,9 +96,6 @@

    public final static String TEXT_XML = "text/xml";

    /**
     * 正斜杠
     */
    public final static String SLASH = "/";

    public final static String IN = "in";
@@ -157,53 +136,26 @@

    public final static String SYS_META = "sysmeta";

    /**
     * 版本号
     */
    public final static String VERSION = "1.0.0";

    /**
     * 令牌键
     */
    public final static String TOKEN_KEY = "token";

    /**
     * Cookie中令牌键
     */
    public final static String TOKEN_COOKIE_KEY = "token";

    /**
     * 文本编码方式
     */
    public final static String TEXT_ENCODER = "UTF-8";

    /**
     * 总质检
     */
    public final static String CHECK_MAIN = "checkMain";

    /**
     * Object对象
     */
    public final static String OBJECT = "java.lang.Object";

    /**
     * Cookie中druid键
     */
    public final static String DRUID_COOKIE_KEY = "JSESSIONID";

    public final static String YES = "YES";

    public final static String NO = "NO";

    /**
     * 线路
     */
    public final static String ROUTE = "线路";

    /**
     * 穿跨越
     */
    public final static String CROSSING = "穿跨越";

    public final static String LINESTRING = "LINESTRING";
@@ -234,6 +186,12 @@
            "        AUTHORITY[\"EPSG\",\"9122\"]],\r\n" +
            "    AUTHORITY[\"ESRI\",\"104903\"]]";

    public final static String CGCS2000 = "CGCS2000";

    public final static String MOON200 = "GCS_Moon_2000";

    public final static List<String> EPSGS = new ArrayList<>(Arrays.asList("EPSG:4326", "EPSG:4490", "ESRI:104903"));

    /**
     * 16进制
     */

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

@@ -1,10 +1,15 @@
package com.moon.server.service.data;

import com.moon.server.entity.all.StaticData;
import com.moon.server.entity.data.MetaEntity;
import com.moon.server.helper.StringHelper;
import org.apache.commons.logging.Log;
import org.apache.commons.logging.LogFactory;
import org.gdal.gdal.ColorTable;
import org.gdal.gdal.Dataset;
import org.gdal.gdal.gdal;
import org.gdal.ogr.Geometry;
import org.gdal.ogr.ogr;
import org.gdal.osr.SpatialReference;
import org.springframework.stereotype.Service;

@@ -17,12 +22,33 @@
 */
@Service
public class RasterService {
    private SpatialReference sr4326;

    private SpatialReference sr4490;

    private SpatialReference sr104903;

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

    /**
     * 初始化空间引用
     */
    private void initSpatialReference() {
        if (null == sr4326) {
            sr4326 = new SpatialReference();
            sr4326.ImportFromEPSG(StaticData.I4326);

            sr4490 = new SpatialReference();
            sr4490.ImportFromEPSG(StaticData.I4490);

            sr104903 = new SpatialReference(StaticData.MOON_2000_WKT);
        }
    }

    /**
     * 读取栅格信息
     */
    public void ReadRasterInfo(MetaEntity me, String file) {
    public void readRasterInfo(MetaEntity me, String file) {
        Dataset ds = null;
        try {
            File f = new File(file);
@@ -35,7 +61,35 @@
                return;
            }

            SpatialReference sr = ds.GetSpatialRef();
            // EPSG编码
            String epsg = getEpsg(ds);
            // 波段数
            String bands = "" + ds.getRasterCount();
            //vd.Meta.band_type = Enum.GetName(typeof(DataType), ds.GetRasterBand(1).DataType); // 数据类型
            ColorTable colorTable = ds.GetRasterBand(1).GetRasterColorTable();
            // 数据颜色表
            String ct = null == colorTable ? null : colorTable.toString();
            // 高程基准
            String hDatum = null;

            double[] transform = new double[6];
            ds.GetGeoTransform(transform);
            // 分辨率
            String resolution = String.format("%f,%f", transform[1], transform[5]);

            if (!StaticData.EPSGS.contains(epsg)) {
                return;
            }

            Geometry minPoint = getMinPoint(ds);
            Geometry maxPoint = getMaxPoint(ds);
            double xmin = minPoint.GetX(0);
            double ymin = minPoint.GetY(0);
            double xmax = maxPoint.GetX(0);
            double ymax = maxPoint.GetY(0);

            // 四至范围
            String geom = String.format("ST_GeomFromText('POLYGON ((%f %f,%f %f,%f %f,%f %f,%f %f))')", xmin, ymax, xmax, ymax, xmax, ymin, xmin, ymin, xmin, ymax);
        } catch (Exception ex) {
            log.error(ex.getMessage(), ex);
        } finally {
@@ -44,4 +98,86 @@
            }
        }
    }

    /**
     * 获取EPSG编码
     */
    private String getEpsg(Dataset ds) {
        SpatialReference sr = ds.GetSpatialRef();
        if (sr != null) {
            // 坐标系统
            String coorSys = sr.GetName();
            if (StaticData.MOON200.equals(coorSys)) {
                return "ESRI:" + StaticData.I104903;
            } else {
                String code = sr.GetAuthorityCode(null);
                return StringHelper.isEmpty(code) ? null : "EPSG:" + code;
            }
        }

        return null;
    }

    /**
     * 获取Dataset的最小点
     */
    private Geometry getMinPoint(Dataset ds) {
        double[] transform = new double[6];
        ds.GetGeoTransform(transform);

        String epsg = ds.GetSpatialRef().GetAuthorityCode(null);
        double xMin = transform[0];
        double yMin = transform[3] - ds.getRasterYSize() * transform[1];

        Geometry point = new Geometry(ogr.wkbPoint);
        point.AddPoint(xMin, yMin, 0);

        return transform(ds, point, epsg);
    }

    /**
     * 获取Dataset的最大点
     */
    private Geometry getMaxPoint(Dataset ds) {
        double[] transform = new double[6];
        ds.GetGeoTransform(transform);

        String epsg = ds.GetSpatialRef().GetAuthorityCode(null);
        double xMax = transform[0] + (ds.getRasterXSize() * transform[1]);
        double yMax = transform[3];

        Geometry point = new Geometry(ogr.wkbPoint);
        point.AddPoint(xMax, yMax, 0);

        return transform(ds, point, epsg);
    }

    /**
     * 坐标转换
     */
    private Geometry transform(Dataset ds, Geometry point, String epsg) {
        this.initSpatialReference();

        if (StringHelper.isEmpty(epsg)) {
            point.AssignSpatialReference(sr104903);
            return point;
        }
        if (String.valueOf(StaticData.I4326).equals(epsg)) {
            point.AssignSpatialReference(sr4326);
            return point;
        }
        if (String.valueOf(StaticData.I4490).equals(epsg)) {
            point.AssignSpatialReference(sr4490);
            return point;
        }

        point.AssignSpatialReference(ds.GetSpatialRef());
        if (ds.GetSpatialRef().GetName().contains(StaticData.CGCS2000)) {
            point.TransformTo(sr4490);
        } else {
            point.TransformTo(sr4326);
        }

        return point;
    }
}