From 1f22a86346c18ccb29a4a712c4acb05c013b33d4 Mon Sep 17 00:00:00 2001
From: 13693261870 <252740454@qq.com>
Date: 星期一, 22 七月 2024 10:50:57 +0800
Subject: [PATCH] 修改高程和水面的计算方法
---
SimuTools/Tools/Handle.cs | 80 +++++++++++++++++----------------------
1 files changed, 35 insertions(+), 45 deletions(-)
diff --git a/SimuTools/Tools/Handle.cs b/SimuTools/Tools/Handle.cs
index 8b45ace..6882a03 100644
--- a/SimuTools/Tools/Handle.cs
+++ b/SimuTools/Tools/Handle.cs
@@ -36,8 +36,8 @@
if (!Directory.Exists(temp)) Directory.CreateDirectory(temp);
CopeTerrain(terrainFile, outPath, layer);
- //CopeWater(waterPath, outPath, layer);
- CopeFlow(flowPath, outPath, layer);
+ CopeWater(waterPath, outPath, layer);
+ //CopeFlow(flowPath, outPath, layer);
CopeLayerJson(outPath, layer);
if (Directory.Exists(temp)) Directory.Delete(temp, true);
@@ -56,7 +56,7 @@
if (null == ds || 0 == ds.RasterCount || null == ds.GetSpatialRef()) return;
SetTerrainInfo(ds, layer);
- //CreateTerrainPng(ds, layer, outPath);
+ CreateTerrainPng(ds, layer, outPath);
}
finally
{
@@ -92,8 +92,8 @@
foreach (int[] sizes in layer.terrain.size)
{
string tif = Path.Combine(tempPath, DateTime.Now.Ticks.ToString() + ".tif");
- Resample(ds.GetDescription(), tif, sizes[0], sizes[1]);
- //Resample(ds, tif, sizes[0], sizes[1]);
+ //Resample(ds.GetDescription(), tif, sizes[0], sizes[1]);
+ Resample(ds, tif, sizes[0], sizes[1]);
if (!File.Exists(tif))
{
continue;
@@ -106,6 +106,32 @@
/// <summary>
/// 閲嶉噰鏍�
+ /// </summary>
+ private static void Resample(Dataset ds, string dest, int width, int height)
+ {
+ // https://blog.51cto.com/u_16099346/6691820
+ //string args = string.Format("-t_srs {0} -ts {1} {2} -r {3} -of {4}", "EPSG:4326", width, height, "bilinear", "GTiff");
+ //string[] options = Gdal.ParseCommandLine(args);
+ string[] options = new string[]
+ {
+ "-t_srs",
+ "EPSG:4326",
+ "-ts",
+ width.ToString(),
+ height.ToString(),
+ "-r",
+ "bilinear",
+ "-of",
+ "GTiff"
+ };
+ GDALWarpAppOptions warpAppOptions = new GDALWarpAppOptions(options);
+
+ Dataset destDs = Gdal.Warp(dest, new Dataset[] { ds }, warpAppOptions, null, null);
+ destDs.Dispose();
+ }
+
+ /// <summary>
+ /// 閲嶉噰鏍� *
/// </summary>
private static void Resample(string source, string dest, int width, int height)
{
@@ -174,22 +200,16 @@
Graphics graphic = Graphics.FromImage(image);
graphic.Clear(Color.Transparent); // 濉厖閫忔槑鑹�
- //double perHeight = (layer.extension.maxHeight - layer.extension.minHeight) * 100 / 65536;
- //int val = Convert.ToInt32(buffer[x * y]* 100); //Convert.ToInt32((buffer[x * y] - layer.extension.minHeight) * 100);
- //int val = Convert.ToInt32((buffer[x * y] - layer.extension.minHeight) * 100 / perHeight);
-
+ double perHeight = (layer.extension.maxHeight - layer.extension.minHeight) * 100 / 65536;
for (int x = 0; x < width; x++)
{
for (int y = 0; y < height; y++)
{
int offset = x + y * width;
- if (float.IsNaN(buffer[offset]) || buffer[offset] < -999)
- {
- //image.SetPixel(x, y, Color.Transparent);
- continue;
- }
+ if (float.IsNaN(buffer[offset]) || buffer[offset] < -999) continue;
- int val = Convert.ToInt32(buffer[offset] * 100);
+ //int val = Convert.ToInt32(buffer[offset] * 100);
+ int val = Convert.ToInt32((buffer[offset] - layer.extension.minHeight) * 100 / perHeight);
int r = val / 65536;
int g = (val - r * 65536) / 256;
int b = val % 256;
@@ -198,7 +218,6 @@
image.SetPixel(x, y, color);
}
}
-
image.Save(png, System.Drawing.Imaging.ImageFormat.Png);
}
finally
@@ -439,35 +458,6 @@
{
sw.Write(json);
}
- }
- #endregion
-
- #region 鏆傛椂涓嶇敤
- /// <summary>
- /// 閲嶉噰鏍� *
- /// </summary>
- private static void Resample(Dataset ds, string dest, int width, int height)
- {
- // https://blog.51cto.com/u_16099346/6691820
- //string args = string.Format("-t_srs {0} -ts {1} {2} -r {3} -of {4}", "EPSG:4326", width, height, "bilinear", "GTiff");
- //string[] options = Gdal.ParseCommandLine(args);
-
- string[] options = new string[]
- {
- "-t_srs",
- "EPSG:4326",
- "-ts",
- width.ToString(),
- height.ToString(),
- "-r",
- "bilinear",
- "-of",
- "GTiff"
- };
- GDALWarpAppOptions warpAppOptions = new GDALWarpAppOptions(options);
-
- Dataset destDs = Gdal.Warp(dest, new Dataset[] { ds }, warpAppOptions, null, null);
- destDs.Dispose();
}
#endregion
}
--
Gitblit v1.9.3