设置地形信息

13693261870

2024-07-18 9f044bbbf3fa5b67fd1d120140f15f0b3729f3e8

设置地形信息

 SimuTools/Domain/Extension.cs

@@ -45,5 +45,11 @@
            this.minHeight = minHeight;
            this.maxHeight = maxHeight;
        }

        public void SetHeight(double minHeight, double maxHeight)
        {
            this.minHeight = minHeight;
            this.maxHeight = maxHeight;
        }
    }
}

 SimuTools/Tools/GdalHelper.cs

@@ -86,7 +86,7 @@
        /**
         * 获取Dataset的最小点
         */
        private static Geometry GetMinPoint(Dataset ds)
        public static Geometry GetMinPoint(Dataset ds)
        {
            double[] transform = new double[6];
            ds.GetGeoTransform(transform);
@@ -103,7 +103,7 @@
        /**
         * 获取Dataset的最大点
         */
        private static Geometry GetMaxPoint(Dataset ds)
        public static Geometry GetMaxPoint(Dataset ds)
        {
            /*
             * transform[0] 左上角x坐标

 SimuTools/Tools/Handle.cs

@@ -1,4 +1,6 @@
using Newtonsoft.Json;
using OSGeo.GDAL;
using OSGeo.OGR;
using SimuTools.Domain;
using System;
using System.Collections.Generic;
@@ -18,10 +20,8 @@
        /// </summary>
        public static void Run(string terrainFile, string waterPath, string flowPath, string outPath)
        {
            Layer layer = new Layer();
            layer.terrain = new Terrain();
            Domain.Layer layer = new Domain.Layer();
            layer.duration = new Duration();
            layer.extension = new Extension();
            layer.waters = new Water();

            CopeTerrain(terrainFile, outPath, layer);
@@ -33,15 +33,41 @@
        /// <summary>
        /// 处理地形
        /// </summary>
        private static void CopeTerrain(string terrainFile, string outPath, Layer layer)
        private static void CopeTerrain(string terrainFile, string outPath, Domain.Layer layer)
        {
            Dataset ds = null;
            try
            {
                ds = Gdal.Open(terrainFile, Access.GA_ReadOnly);
                if (null == ds || 0 == ds.RasterCount || null == ds.GetSpatialRef()) return;

                SetTerrainInfo(ds, layer);
            }
            finally
            {
                if (null != ds) ds.Dispose();
            }
        }

        /// <summary>
        /// 设置地形信息
        /// </summary>
        private static void SetTerrainInfo(Dataset ds, Domain.Layer layer)
        {
            Geometry minPoint = GdalHelper.GetMinPoint(ds);
            Geometry maxPoint = GdalHelper.GetMaxPoint(ds);
            layer.extension = new Extension(minPoint.GetX(0), minPoint.GetY(0), maxPoint.GetX(0), maxPoint.GetY(0));

            Band band = ds.GetRasterBand(1);
            double[] mm = new double[2];
            band.ComputeRasterMinMax(mm, 0);
            layer.extension.SetHeight(mm[0], mm[1]);
        }

        /// <summary>
        /// 处理水面
        /// </summary>
        private static void CopeWater(string waterPath, string outPath, Layer layer)
        private static void CopeWater(string waterPath, string outPath, Domain.Layer layer)
        {

        }
@@ -49,7 +75,7 @@
        /// <summary>
        /// 处理流速流向
        /// </summary>
        private static void CopeFlow(string flowPath, string outPath, Layer layer)
        private static void CopeFlow(string flowPath, string outPath, Domain.Layer layer)
        {

        }
@@ -57,7 +83,7 @@
        /// <summary>
        /// 处理元数据
        /// </summary>
        private static void CopeLayerJson(string outPath, Layer layer)
        private static void CopeLayerJson(string outPath, Domain.Layer layer)
        {
            if (null == layer) return;