From ecfe110922f3c542df7aaef9ac9e26cdb0b73e55 Mon Sep 17 00:00:00 2001
From: 13693261870 <252740454@qq.com>
Date: 星期六, 20 七月 2024 09:52:58 +0800
Subject: [PATCH] 1
---
SimuTools/Tools/Handle.cs | 77 +++++++++++++++++++++-----------------
1 files changed, 43 insertions(+), 34 deletions(-)
diff --git a/SimuTools/Tools/Handle.cs b/SimuTools/Tools/Handle.cs
index 7ada463..74796fc 100644
--- a/SimuTools/Tools/Handle.cs
+++ b/SimuTools/Tools/Handle.cs
@@ -36,9 +36,9 @@
if (!Directory.Exists(temp)) Directory.CreateDirectory(temp);
CopeTerrain(terrainFile, outPath, layer);
- CopeWater(waterPath, outPath, layer);
+ //CopeWater(waterPath, outPath, layer);
CopeFlow(flowPath, outPath, layer);
- //CopeLayerJson(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,7 +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.GetDescription(), tif, sizes[0], sizes[1]);
+ Resample(ds, tif, sizes[0], sizes[1]);
if (!File.Exists(tif))
{
continue;
@@ -105,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)
{
@@ -181,8 +208,8 @@
{
for (int y = 0; y < height; y++)
{
- int offset = x + y * height;
- if (float.IsNaN(buffer[offset]) || buffer[offset] == -9999)
+ int offset = x + y * width;
+ if (float.IsNaN(buffer[offset]) || buffer[offset] < -999)
{
//image.SetPixel(x, y, Color.Transparent);
continue;
@@ -388,20 +415,26 @@
Graphics graphic = Graphics.FromImage(image);
graphic.Clear(Color.Transparent); // 濉厖閫忔槑鑹�
+ // 鐢� R閫氶亾琛ㄧず锛屾祦鍚戜负褰掍竴鍖栫殑浜岀淮鍚戦噺锛坸,y锛夛紝G閫氶亾琛ㄧず涓� x *255 , B閫氶亾琛ㄧず涓� y * 255
for (int x = 0; x < width; x++)
{
for (int y = 0; y < height; y++)
{
- //int offset = x + y * height;
- //if (float.IsNaN(vxBuffer[offset]) || vxBuffer[offset] == -9999) continue;
+ int offset = x + y * width;
+ if (float.IsNaN(vxBuffer[offset]) || vxBuffer[offset] < -999 ||
+ float.IsNaN(vyBuffer[offset]) || vyBuffer[offset] < -999) continue;
+ float vx = vxBuffer[offset], vy = vyBuffer[offset];
//int val = Convert.ToInt32(vxBuffer[offset] * 100);
//int r = val / 65536;
//int g = (val - r * 65536) / 256;
//int b = val % 256;
+ int r = 0;
+ int g = 0;
+ int b = 0;
- //Color color = Color.FromArgb(127, r, g, b);
- //image.SetPixel(x, y, color);
+ Color color = Color.FromArgb(127, r, g, b);
+ image.SetPixel(x, y, color);
}
}
@@ -432,30 +465,6 @@
{
sw.Write(json);
}
- }
- #endregion
-
- #region 鏆傛椂涓嶇敤
- /// <summary>
- /// 閲嶉噰鏍� *
- /// </summary>
- private static void Resample(Dataset ds, string dest, int width, int height)
- {
- // 璁剧疆Warp鐨勯�夐」锛歨ttps://blog.csdn.net/qq_43210879/article/details/121350561
- string[] options = new string[] {
- "format=GTiff",
- "width=" + width,
- "height=" + height,
- "dstSRS=EPSG:4326"
- };
-
- OSGeo.GDAL.Driver driver = Gdal.GetDriverByName("GTiff");
- Dataset destDs = driver.Create(dest, width, height, ds.RasterCount, ds.GetRasterBand(1).DataType, null);
-
- GDALWarpAppOptions warpAppOptions = new GDALWarpAppOptions(options);
- Gdal.Warp(destDs, new Dataset[] { ds }, warpAppOptions, null, null);
-
- destDs.Dispose();
}
#endregion
}
--
Gitblit v1.9.3