洪涝淹没分析输出淹没范围图、深度图及面积体积等信息【转】

接上一篇博客《洪涝有源淹没算法及结果分析》。可以输出洪涝淹没范围图、深度图以及淹没面积等信息，下一步输出历时图。代码很简单，下面直接给出源代码，欢迎大家留言交流。

[csharp] view plain​ copy

1. 1. 
      
2. /// <summary>
3. /// 输出洪涝淹没范围图
4. /// </summary>
5. public void OutPutFloodRegion()
6. {
7. 
   
8. //创建洪涝淹没范围影像
9. string m_FloodRegionPath = System.IO.Path.GetDirectoryName(System.Windows.Forms.Application.ExecutablePath) + "\\FloodSimulation\\FloodedRegion.tif";
10. if (System.IO.File.Exists(m_FloodRegionPath))
11. {
12. System.IO.File.Delete(m_FloodRegionPath);
13. }
14. 
    
15. //在GDAL中创建影像,先需要明确待创建影像的格式,并获取到该影像格式的驱动
16. driver = Gdal.GetDriverByName("GTiff");
17. //调用Creat函数创建影像
18. m_FloodSimulationRegionDataSet = driver.Create(m_FloodRegionPath, m_XSize, m_YSize, 1, DataType.GDT_Float32, null);
19. //设置影像属性
20. m_FloodSimulationRegionDataSet.SetGeoTransform(m_adfGeoTransform); //影像转换参数
21. m_FloodSimulationRegionDataSet.SetProjection(m_DEMDataSet.GetProjectionRef()); //投影
22. 
    
23. //输出影像
24. m_FloodSimulationRegionDataSet.GetRasterBand(1).WriteRaster(0, 0, m_XSize, m_YSize, m_FloodRegionBuffer, m_XSize, m_YSize, 0, 0);
25. m_FloodSimulationRegionDataSet.GetRasterBand(1).FlushCache();
26. m_FloodSimulationRegionDataSet.FlushCache();
27. 
    
28. }
29. 
    
30. /// <summary>
31. /// 输出洪涝淹没深度图
32. /// </summary>
33. public void OutPutFloodDepth()
34. {
35. 
    
36. for (int i = 0; i < m_XSize; i++)
37. {
38. for (int j = 0; j < m_YSize; j++)
39. {
40. Point m_point = new Point();
41. m_point.X = i;
42. m_point.Y = j;
43. if (m_FloodRegionBuffer[getIndex(m_point)] == 1)
44. {
45. int evaluation = m_DEMdataBuffer[getIndex(m_point)]; //淹没格网高程值
46. int value = pFloodLevel - evaluation;
47. m_FloodDepthBuffer[getIndex(m_point)] = value;
48. }
49. }
50. 
    
51. }
52. //输出洪涝淹没深度图
53. string m_FloodDepthPath = System.IO.Path.GetDirectoryName(System.Windows.Forms.Application.ExecutablePath) + "\\FloodSimulation\\FloodedDepth.tif";
54. if (System.IO.File.Exists(m_FloodDepthPath))
55. {
56. System.IO.File.Delete(m_FloodDepthPath);
57. }
58. 
    
59. //在GDAL中创建影像,先需要明确待创建影像的格式,并获取到该影像格式的驱动
60. driver = Gdal.GetDriverByName("GTiff");
61. 
    
62. //调用Creat函数创建影像
63. m_FloodSimulationDepthDateSet = driver.Create(m_FloodDepthPath, m_XSize, m_YSize, 1, DataType.GDT_Float32, null);
64. //设置影像属性
65. m_FloodSimulationDepthDateSet.SetGeoTransform(m_adfGeoTransform); //影像转换参数
66. m_FloodSimulationDepthDateSet.SetProjection(m_DEMDataSet.GetProjectionRef()); //投影
67. 
    
68. //输出影像
69. m_FloodSimulationDepthDateSet.GetRasterBand(1).WriteRaster(0, 0, m_XSize, m_YSize, m_FloodDepthBuffer, m_XSize, m_YSize, 0, 0);
70. m_FloodSimulationDepthDateSet.GetRasterBand(1).FlushCache();
71. m_FloodSimulationDepthDateSet.FlushCache();
72. 
    
73. }
74. /// <summary>
75. /// 输出洪涝淹没面积及水量
76. /// </summary>
77. public void OutPutFloodInfo()
78. {
79. int count = 0;
80. for (int i = 0; i < m_XSize; i++)
81. {
82. for (int j = 0; j < m_YSize; j++)
83. {
84. Point m_point = new Point();
85. m_point.X = i;
86. m_point.Y = j;
87. if (m_FloodRegionBuffer[getIndex(m_point)] == 1)
88. {
89. count++;
90. }
91. }
92. 
    
93. }
94. //统计面积
95. double S = count * 90; //平方米
96. if (S > 10000)
97. {
98. S = S / 10000; //公顷
99. }
100. MessageBox.Show("淹没面积：" + S.ToString() + "公顷");
101. }
     1. 
        

结果图：范围图与深度图