修复时间轴问题

dcb

2025-05-22 f984496e9bb3612ce38040a5dd0be548b181e971

 src/main/java/com/se/nsl/service/TestService.java

@@ -22,6 +22,10 @@
import java.io.File;
import java.io.IOException;
import java.text.SimpleDateFormat;
import java.time.Instant;
import java.time.LocalDateTime;
import java.time.ZoneId;
import java.time.format.DateTimeFormatter;
import java.util.*;
import java.util.List;

@@ -99,10 +103,14 @@
    public void setTerrainInfo(Dataset ds, LayerDto layer) {
        Geometry minPoint = GdalHelper.getMinPoint(ds);
        Geometry maxPoint = GdalHelper.getMaxPoint(ds);
        double minx = ComHelper.getMinVal(minPoint.GetX(0), 10000000);
        double miny = ComHelper.getMinVal(minPoint.GetY(0), 10000000);
        double maxx = ComHelper.getMaxVal(maxPoint.GetX(0) + MAX_X_OFFSET, 10000000);
        double maxy = ComHelper.getMaxVal(maxPoint.GetY(0), 10000000);
//        double minx = ComHelper.getMinVal(minPoint.GetX(0), 10000000);
//        double miny = ComHelper.getMinVal(minPoint.GetY(0), 10000000);
//        double maxx = ComHelper.getMaxVal(maxPoint.GetX(0) + MAX_X_OFFSET, 10000000);
//        double maxy = ComHelper.getMaxVal(maxPoint.GetY(0), 10000000);
        double minx = Double.MAX_VALUE;
        double miny = Double.MAX_VALUE;
        double maxx = Double.MIN_VALUE;
        double maxy = Double.MIN_VALUE;
        //layer.setExtension(new ExtensionDto(minx, miny, maxx, maxy, Double.MAX_VALUE, Double.MIN_VALUE));

        Band band = ds.GetRasterBand(1);
@@ -193,6 +201,21 @@
    }

    public void setWaterData(LayerDto layer, List<String> files) {
        DateTimeFormatter formatter = DateTimeFormatter.ofPattern("yyyyMMddHHmmss");
        for (String file : files) {
            String fileName = ComHelper.getNameWithExt(file);
            // 解析为 LocalDateTime（默认基于系统时区，可能需指定时区）
            LocalDateTime dateTime = LocalDateTime.parse(fileName, formatter);
            // 转换为 UTC 时间戳（推荐，避免时区歧义）
            Instant utcInstant = dateTime.atZone(ZoneId.of("UTC")).toInstant();
            long timestamp = utcInstant.toEpochMilli();
            layer.getWaters().getData().add(timestamp);
        }
        layer.getDuration().setStart(layer.getWaters().getData().get(0));
        layer.getDuration().setEnd(layer.getWaters().getData().get(layer.getWaters().getData().size() - 1));
    }

    /*public void setWaterData(LayerDto layer, List<String> files) {
        Calendar calendar = Calendar.getInstance();
        calendar.setTime(new Date());
        calendar.set(Calendar.MILLISECOND, 0);
@@ -217,15 +240,11 @@
        }
        layer.getDuration().setStart(layer.getWaters().getData().get(0));
        layer.getDuration().setEnd(layer.getWaters().getData().get(layer.getWaters().getData().size() - 1));
    }
    }*/

    public void setWaterHeight(LayerDto layer, List<String> files) {
        int c = files.size(), step = files.size() / 10;
        if (step < 1) step = 1;
        double minx = Double.MAX_VALUE;
        double miny = Double.MAX_VALUE;
        double maxx = Double.MIN_VALUE;
        double maxy = Double.MIN_VALUE;
        for (int i = 0; i < c; i += step) {
            Dataset ds = null;
            try {
@@ -235,16 +254,6 @@
                double[] mm = new double[2];
                ds.GetRasterBand(1).ComputeRasterMinMax(mm, 0);
                layer.getWaters().setHeight(mm[0], mm[1]);

                double[] bbox = readTifBbox(ds);
                double minLon = bbox[0];
                double maxLon = bbox[1];
                double minLat = bbox[2];
                double maxLat = bbox[3];
                if (minx > minLon) minx = minLon;
                if (miny > minLat) miny = minLat;
                if (maxx < maxLon) maxx = maxLon;
                if (maxy < maxLat) maxy = maxLat;
            } finally {
                if (null != ds) ds.delete();
            }
@@ -254,29 +263,6 @@
        layer.getExtension().setMaxHeight(ComHelper.getMaxVal(layer.getExtension().getMaxHeight(), 1000000));
        layer.getExtension().setMinHeight(ComHelper.getMaxVal(layer.getExtension().getMinHeight(), 1000000));
        layer.getExtension().setDiffer();
    }

    private double[] readTifBbox(Dataset ds) {
        // 1. 获取图像尺寸
        int width = ds.getRasterXSize();
        int height = ds.getRasterYSize();

        // 2. 获取GeoTransform参数
        double[] geoTransform = new double[6];
        ds.GetGeoTransform(geoTransform);

        // 3. 解析GeoTransform参数（经纬度坐标）
        double originLon = geoTransform[0];  // 左上角经度
        double originLat = geoTransform[3];  // 左上角纬度
        double pixelWidth = geoTransform[1];  // 经度方向分辨率（度/像素）
        double pixelHeight = geoTransform[5]; // 纬度方向分辨率（度/像素，通常为负）

        // 4. 计算四至范围（经纬度）
        double minLon = Math.min(originLon, originLon + width * pixelWidth);
        double maxLon = Math.max(originLon, originLon + width * pixelWidth);
        double minLat = Math.min(originLat, originLat + height * pixelHeight);
        double maxLat = Math.max(originLat, originLat + height * pixelHeight);
        return new double[] {minLon, maxLon, minLat, maxLat};
    }

    public void processWaters(ResultDto dto, List<String> files, LayerDto layer) {
@@ -444,7 +430,7 @@
        for (int[] sizes : layer.getTerrain().getSize()) {
            String name = ComHelper.getNameWithExt(ds.GetDescription()) + "_" + sizes[0] + "_" + sizes[1];
            String tif = dto.getTemp() + File.separator + name + ".tif";
            ComHelper.Resample(ds, tif, sizes[0], sizes[1], layer);
//            ComHelper.Resample(ds, tif, sizes[0], sizes[1], layer);
            if (!new File(tif).exists()) continue;

            String png = flowPath + File.separator + sizes[0] + "_" + sizes[1] + ".png";