利用GIS技术计算洪水淹没区范围一直是灾害评估研究中的一个热点问题。通过给出两种情形下基于种子蔓延算法的淹没区计算方法 ,即有源淹没和无源淹没。淹没区计算精度主要取决于空间数据精度的优劣 ;种子蔓延算法及探测分辨率决定了整个模型的效率。文中最后给出了该模型在“水利综合管理信息系统”中得到验证和实现的实例。
In this paper, a model of flood area display and damage evaluation is discussed, which adopts seed spread algorithm to calculate and visualize flood area based on a DEM data set and under the condition of a given water level of flood. The raster image containing flood area is overlaid with raster data converted from thematic vector maps such as land use map and related to other socio economic data, and then quantitative evaluation of flood damage is achieved. There are two different cases, the so called “non source flood” and “source flood”, for determining flood area based on DEM and under the given water level condition. Non source flood corresponds to the case with well distributed and large area rainfall where all the low lying land may have a flood disaster, and in this case all the points with elevations below the given water level should be included in the flood area. The “source flood” is the case of flood (e.g. from bank burst) flushing through the surrounding lands that needs to consider “circulating” condition. Because flood may be obstructed by ring structures or high lands, thus it can only cover the place where it flows and reaches. One of the computer processing methods suitable for source flood is the seed spread algorithm. The principle of this method is to select a representative pixel as a seed and to examine its contiguous pixels outward along 4 or 8 directions. The pixels, contiguous to it and meeting the given conditions, become seeds, and then contiguous pixels to the new seeds will be examined in the same way. This process repeats until an aggregation of the pixels which meet the given conditions on a certain plane are obtained. Adopting seed spread algorithm for source flood area calculation is just to achieve an aggregation of pixels, which satisfies the connectivity requirement and the given conditions of water level. The initial position of seed pixel should usually be selected on the characteristic points, e.g. adjacent to a reservoir dam or stream boundary. The precision of flood area calculation is mainly determined by spatial resolution of the DEM, the area on the ground is represented by each pixel; whereas the accuracy of flood area calculation, to a great extent, is dependent on the precision of the DEM’s pixel value or elevation. Once the area below the flood water level has been calculated, the flood damage evaluation can be calculated by overlaying the flood range with thematic maps and by relating to other socio economic data. All thematic maps were converted to raster data because algorithm of the overlay analysis implemented in raster format is relatively simple and computing is fast. In order to improve accuracy of the flood damage evaluation it is necessary to accurately register and georefer all raster layers before overlay analysis. The calculation results of flood damage evaluation are generally represented by tabulation, which is a kind of classified statistics of flood damage data, e.g. a table representing total area of every land use type in different administrative divisions. The results of flood area calculation and damage evaluation can also be visually and dynamically provided in the forms of graph, image, 3D or virtual reality. The model has been successfully applied in the “Zhejiang Provincial Water Conservancy Management System” developed for a national project.
[1] Zhou Chenghu.F lood H azard A ssessm en t Info rm ation System[M].B eijing:Ch ina Science and T echno logy P ress,1993.58-85.(In Ch inese)[周成虎.洪水灾害评估信息系统研究[M].北京:中国科学技术出版社,1993.58~85.]
[2] L iu Xue,W ang X ingku i,W ang Guangqian.Study of G IS-based m odeling of sp atial p rocesses smi u lation[J].J ou rna l of Im ag e and G rap h ics,1999,4(6):476-480.(In Ch inese)[刘学,王兴奎,王光谦.基于G IS的空间过程模拟建模方法研究[J].中国图像图形学报,1999,4(6):476-480.]
[3] M arinellil,M ichel R,B eaudo in A et a.lF lood m app ing u sing ER S tandem coherence mi age:a case study in Sou thern F rance[J].ESA S P,1997(414):531-536.
[4] W ah lstrom E,L oague K,Kyriak id is P C.H yd ro log ic respon se:Kaho’o law e,H aw aii[J].J ou rna l of E nv ironm en ta l Q ua lity,1999,28(2):481-492.
[5] H allam C A,L anfear K J,Salisbu ry J M et a.lGeog rap h ic Info rm ation System s(G IS)and Im age P rocessing fo r Sto rm-w ater M anagem en t M odeling U sing XPSWMM in T allahassee,F lo rida[A].In:P roceed ing s AW RA A nnual Sym po sium’98[C].G IS and W ater R esou rces,1998.305-313.
[6] Ke Zhengy.iD ig ital T errain M odel[M].B eijing:Ch ina Science and T echno logy P ress,1993.51-67.(In Ch inese)[柯正谊.数字地面模型[M].北京:中国科学技术出版社,1993.51-67.]