From 4e7ae7e903a07389607cf47ad8f2bbd5e8ccd8d6 Mon Sep 17 00:00:00 2001
From: 13693261870 <252740454@qq.com>
Date: 星期三, 30 八月 2023 10:18:50 +0800
Subject: [PATCH] 元数据查询接口添加传感器类型、元数据类型、开始日期、结束日期参数
---
src/main/java/com/moon/server/service/data/RasterService.java | 127 +++++++++++++++++++++++++++---------------
1 files changed, 82 insertions(+), 45 deletions(-)
diff --git a/src/main/java/com/moon/server/service/data/RasterService.java b/src/main/java/com/moon/server/service/data/RasterService.java
index be905d5..21f89d2 100644
--- a/src/main/java/com/moon/server/service/data/RasterService.java
+++ b/src/main/java/com/moon/server/service/data/RasterService.java
@@ -1,13 +1,14 @@
package com.moon.server.service.data;
import com.moon.server.entity.all.StaticData;
-import com.moon.server.entity.data.MetaEntity;
+import com.moon.server.entity.data.MetaFileEntity;
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.gdalconst.gdalconst;
import org.gdal.ogr.Geometry;
import org.gdal.ogr.ogr;
import org.gdal.osr.SpatialReference;
@@ -48,7 +49,7 @@
/**
* 璇诲彇鏍呮牸淇℃伅
*/
- public void readRasterInfo(MetaEntity me, String file) {
+ public void readRasterInfo(MetaFileEntity mf, String file) {
Dataset ds = null;
try {
File f = new File(file);
@@ -56,31 +57,44 @@
return;
}
- ds = gdal.Open(file, 0);
- if (null == ds || ds.getRasterCount() == 0) {
+ ds = gdal.Open(file, gdalconst.GA_ReadOnly);
+ if (null == ds || ds.getRasterCount() == 0 || null == ds.GetSpatialRef()) {
return;
}
- // EPSG缂栫爜
- String epsg = getEpsg(ds);
+ SpatialReference sr = ds.GetSpatialRef();
+ // 鍧愭爣绯荤粺
+ mf.setCoorSys(sr.GetName());
+ if (StaticData.MOON200.equals(mf.getCoorSys())) {
+ // EPSG缂栫爜
+ mf.setEpsg("EPSG:" + StaticData.I104903);
+ } else {
+ // EPSG缂栫爜锛歅ROJCS銆丟EOGCS銆丟EOGCS|UNIT 鎴� NULL
+ String code = sr.GetAuthorityCode(null);
+ mf.setEpsg(StringHelper.isEmpty(code) ? null : "EPSG:" + code);
+ }
+
+ // 琛屽垪鏁�
+ mf.setGridsize(String.format("%d,%d", ds.getRasterXSize(), ds.getRasterYSize()));
// 娉㈡鏁�
- String bands = "" + ds.getRasterCount();
- //vd.Meta.band_type = Enum.GetName(typeof(DataType), ds.GetRasterBand(1).DataType); // 鏁版嵁绫诲瀷
- ColorTable colorTable = ds.GetRasterBand(1).GetRasterColorTable();
+ mf.setBands(String.valueOf(ds.getRasterCount()));
+ // 鏁版嵁绫诲瀷
+ int dataType = ds.GetRasterBand(1).GetRasterDataType();
+ mf.setBandType(getDataType(dataType));
// 鏁版嵁棰滆壊琛�
- String ct = null == colorTable ? null : colorTable.toString();
+ ColorTable colorTable = ds.GetRasterBand(1).GetRasterColorTable();
+ mf.setCt(null == colorTable ? null : colorTable.toString());
// 楂樼▼鍩哄噯
- String hDatum = null;
+ mf.sethDatum(null);
- double[] transform = new double[6];
- ds.GetGeoTransform(transform);
// 鍒嗚鲸鐜�
- String resolution = String.format("%f,%f", transform[1], transform[5]);
+ double[] tr = new double[6];
+ ds.GetGeoTransform(tr);
+ mf.setResolution(String.format("%f,%f", tr[1], tr[5]));
- if (!StaticData.EPSGS.contains(epsg)) {
+ if (tr[StaticData.I0] == 0.0 && tr[StaticData.I1] == 1.0 && tr[StaticData.I2] == 0.0 && tr[StaticData.I3] == 0.0 && tr[StaticData.I4] == 0.0 && tr[StaticData.I5] == 1.0) {
return;
}
-
Geometry minPoint = getMinPoint(ds);
Geometry maxPoint = getMaxPoint(ds);
double xmin = minPoint.GetX(0);
@@ -90,6 +104,7 @@
// 鍥涜嚦鑼冨洿
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);
+ mf.setGeom(geom);
} catch (Exception ex) {
log.error(ex.getMessage(), ex);
} finally {
@@ -100,19 +115,41 @@
}
/**
- * 鑾峰彇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;
- }
+ private String getDataType(int dataType) {
+ if (dataType == gdalconst.GDT_Byte) {
+ return "GDT_Byte";
+ }
+ if (dataType == gdalconst.GDT_UInt16) {
+ return "GDT_UInt16";
+ }
+ if (dataType == gdalconst.GDT_Int16) {
+ return "GDT_Int16";
+ }
+ if (dataType == gdalconst.GDT_UInt32) {
+ return "GDT_UInt32";
+ }
+ if (dataType == gdalconst.GDT_Int32) {
+ return "GDT_Int32";
+ }
+ if (dataType == gdalconst.GDT_Float32) {
+ return "GDT_Float32";
+ }
+ if (dataType == gdalconst.GDT_Float64) {
+ return "GDT_Float64";
+ }
+ if (dataType == gdalconst.GDT_CInt16) {
+ return "GDT_CInt16";
+ }
+ if (dataType == gdalconst.GDT_CInt32) {
+ return "GDT_CInt32";
+ }
+ if (dataType == gdalconst.GDT_CFloat32) {
+ return "GDT_CFloat32";
+ }
+ if (dataType == gdalconst.GDT_CFloat64) {
+ return "GDT_CFloat64";
}
return null;
@@ -125,58 +162,58 @@
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);
+ return transform(ds, point);
}
/**
* 鑾峰彇Dataset鐨勬渶澶х偣
*/
private Geometry getMaxPoint(Dataset ds) {
+ /**
+ * transform[0] 宸︿笂瑙抶鍧愭爣
+ * transform[1] 涓滆タ鏂瑰悜鍒嗚鲸鐜�
+ * transform[2] 鏃嬭浆瑙掑害, 0琛ㄧず鍥惧儚 "鍖楁柟鏈濅笂"
+ *
+ * transform[3] 宸︿笂瑙抷鍧愭爣
+ * transform[4] 鏃嬭浆瑙掑害, 0琛ㄧず鍥惧儚 "鍖楁柟鏈濅笂"
+ * transform[5] 鍗楀寳鏂瑰悜鍒嗚鲸鐜�
+ */
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);
+ return transform(ds, point);
}
/**
* 鍧愭爣杞崲
*/
- private Geometry transform(Dataset ds, Geometry point, String epsg) {
+ private Geometry transform(Dataset ds, Geometry point) {
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);
+ point.AssignSpatialReference(ds.GetSpatialRef());
+ if (ds.GetSpatialRef().IsGeographic() > 0) {
return point;
}
- point.AssignSpatialReference(ds.GetSpatialRef());
- if (ds.GetSpatialRef().GetName().contains(StaticData.CGCS2000)) {
+ String srsName = ds.GetSpatialRef().GetName();
+ if (srsName.contains(StaticData.CGCS2000)) {
point.TransformTo(sr4490);
} else {
point.TransformTo(sr4326);
}
+ point.SwapXY();
return point;
}
--
Gitblit v1.9.3