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

冷不防 2022-08-14 01:45 231阅读 0赞

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

/// <summary>
            /// 输出洪涝淹没范围图
            /// </summary>
            public void OutPutFloodRegion()
            {
    
                //创建洪涝淹没范围影像
                string m_FloodRegionPath = System.IO.Path.GetDirectoryName(System.Windows.Forms.Application.ExecutablePath) + "\\FloodSimulation\\FloodedRegion.tif";
                if (System.IO.File.Exists(m_FloodRegionPath))
                {
                    System.IO.File.Delete(m_FloodRegionPath);
                }
    
                //在GDAL中创建影像,先需要明确待创建影像的格式,并获取到该影像格式的驱动
                driver = Gdal.GetDriverByName("GTiff");
                //调用Creat函数创建影像
                m_FloodSimulationRegionDataSet = driver.Create(m_FloodRegionPath, m_XSize, m_YSize, 1, DataType.GDT_Float32, null);
                //设置影像属性
                m_FloodSimulationRegionDataSet.SetGeoTransform(m_adfGeoTransform); //影像转换参数
                m_FloodSimulationRegionDataSet.SetProjection(m_DEMDataSet.GetProjectionRef()); //投影
    
                //输出影像
                m_FloodSimulationRegionDataSet.GetRasterBand(1).WriteRaster(0, 0, m_XSize, m_YSize, m_FloodRegionBuffer, m_XSize, m_YSize, 0, 0);
                m_FloodSimulationRegionDataSet.GetRasterBand(1).FlushCache();
                m_FloodSimulationRegionDataSet.FlushCache();
    
            }
    
            /// <summary>
            /// 输出洪涝淹没深度图
            /// </summary>
            public void OutPutFloodDepth()
            {
    
                for (int i = 0; i < m_XSize; i++)
                {
                    for (int j = 0; j < m_YSize; j++)
                    {
                        Point m_point = new Point();
                        m_point.X = i;
                        m_point.Y = j;
                        if (m_FloodRegionBuffer[getIndex(m_point)] == 1)
                        {
                            int evaluation = m_DEMdataBuffer[getIndex(m_point)]; //淹没格网高程值
                            int value = pFloodLevel - evaluation;
                            m_FloodDepthBuffer[getIndex(m_point)] = value;
                        }
                    }
    
                }
                //输出洪涝淹没深度图
                string m_FloodDepthPath = System.IO.Path.GetDirectoryName(System.Windows.Forms.Application.ExecutablePath) + "\\FloodSimulation\\FloodedDepth.tif";
                if (System.IO.File.Exists(m_FloodDepthPath))
                {
                    System.IO.File.Delete(m_FloodDepthPath);
                }
    
                //在GDAL中创建影像,先需要明确待创建影像的格式,并获取到该影像格式的驱动
                driver = Gdal.GetDriverByName("GTiff");
              
                //调用Creat函数创建影像
                m_FloodSimulationDepthDateSet = driver.Create(m_FloodDepthPath, m_XSize, m_YSize, 1, DataType.GDT_Float32, null);
                //设置影像属性
                m_FloodSimulationDepthDateSet.SetGeoTransform(m_adfGeoTransform); //影像转换参数
                m_FloodSimulationDepthDateSet.SetProjection(m_DEMDataSet.GetProjectionRef()); //投影
       
                //输出影像
                m_FloodSimulationDepthDateSet.GetRasterBand(1).WriteRaster(0, 0, m_XSize, m_YSize, m_FloodDepthBuffer, m_XSize, m_YSize, 0, 0);
                m_FloodSimulationDepthDateSet.GetRasterBand(1).FlushCache();
                m_FloodSimulationDepthDateSet.FlushCache();
    
            }
            /// <summary>
            /// 输出洪涝淹没面积及水量
            /// </summary>
            public void OutPutFloodInfo()
            {
                int count = 0;
                for (int i = 0; i < m_XSize; i++)
                {
                    for (int j = 0; j < m_YSize; j++)
                    {
                        Point m_point = new Point();
                        m_point.X = i;
                        m_point.Y = j;
                        if (m_FloodRegionBuffer[getIndex(m_point)] == 1)
                        {
                            count++;
                        }
                    }
    
                }
                //统计面积
                double S = count * 90; //平方米
                if (S > 10000)
                {
                    S = S / 10000;  //公顷
                }
                MessageBox.Show("淹没面积：" + S.ToString() + "公顷");
            }

结果图：范围图与深度图

![Center][]

![Center 1][]

[Link 1]: http://blog.csdn.net/giser_whu/article/details/41288761
[Center]: /images/20220810/ad0de20cfc4443afa2ba2d1a28849b24.png
[Center 1]: /images/20220810/99deccd47b32476a968f8bf2f9cdf9d8.png