PDF(2839 KB)
基于Copernicus DEM 30的大尺度河道数字高程模型重构方法
李玉建, 赵明成, 李琳, 戴文鸿, 安鹏
长江科学院院报 ›› 2026, Vol. 43 ›› Issue (2) : 201-210.
PDF(2839 KB)
PDF(2839 KB)
基于Copernicus DEM 30的大尺度河道数字高程模型重构方法
Method of Reconstructing Digital Elevation Model for Large-scale River Based on Copernicus DEM 30
为解决在大尺度河道数值模拟中河道遥感影像信息不全和实测地形高程资料部分缺失导致河道数字高程模型难以构建的问题,应用塔里木河阿拉尔—新渠满河段2011年相关水文资料与实测地形高程资料,通过ArcGIS软件结合Google Earth历史影像与Copernicus DEM 30数据绘制河道临水线与外缘线,基于二次插值、平均差值法与局部加权回归算法,结合断面高程数据,完善河道地形高程。通过Mesh Generator组件,重构该河段的河道数字高程模型,并检验数据在MIKE 21水动力-泥沙模块数值模拟中的可行性。结果表明:通过模拟得到的流速、流量-水位关系及河道冲淤变化与实测资料对比,各项指标的误差均符合相关技术规程的允许偏差要求;平均差值法可以弥补主槽与河漫滩高程补衔接不自然的问题;局部加权回归算法能有效平滑河道主槽沿程断面高程数据。研究成果旨在丰富河道地形高程数据的重构方法,为解决大尺度河道DEM难以构建问题提供一种新思路。
[Objective] To address the challenges of incomplete remote sensing imagery and insufficient measured terrain data in large-scale river numerical simulations, this study proposes a novel method for river channel terrain reconstruction that integrates Copernicus DEM 30 data with limited measured data, aiming to solve the problem of constructing large-scale river channel DEMs in data-scarce areas and to verify its applicability in MIKE 21 hydrodynamic-sediment numerical simulations. [Methods] This study took the Alar-Xinquman river section of the Tarim River as the study area and utilized relevant hydrological data and measured terrain elevation data from 2011 as the basic dataset. The specific reconstruction process was as follows: (1) Using ArcGIS software, combined with Google Earth historical imagery and Copernicus DEM 30, the inner and outer bank lines of the river channel were delineated by scaling overlapping imagery proportionally to determine the main channel boundary.(2) The river DEM was corrected stepwise. First, the quadratic interpolation method was applied to densify the elevation data of the main channel cross-sections. Second, the mean difference method was used to calculate the deviations in elevation between measured points and the Copernicus DEM at corresponding locations, enabling an overall vertical correction of the Copernicus DEM data. Next, a locally weighted regression (LOESS) algorithm was introduced to correct the longitudinal cross-section elevations of the main channel, smoothing the riverbed terrain and generating the main channel DEM. Finally, the main channel and floodplain DEM data were merged to construct the complete river channel DEM.(3) The reconstructed DEM was imported into the MIKE 21 hydrodynamic-sediment module. Measured data such as water level, flow, and sediment concentration were selected as boundary conditions for numerical simulation. The accuracy and reliability of the reconstructed DEM were evaluated by comparing the errors between the simulated and measured values. [Results] Comparison between numerical simulation results and measured data revealed that the results for flow velocity, flow-stage relationships, and river channel erosion-deposition tests all met the allowable deviation requirements specified in relevant technical standards. This indicated that the river channel DEM constructed in this study was reasonably suitable for two-dimensional hydrodynamic-sediment numerical simulations. [Conclusion] (1) The river channel terrain reconstruction method proposed in this study, integrating Copernicus DEM 30, quadratic interpolation, mean difference method, and locally weighted regression (LOESS) algorithm, can effectively address the lack of measured underwater terrain data for large-scale rivers. It demonstrates superior continuity and accuracy compared to traditional single-interpolation methods.(2) The hydrodynamic model established based on the reconstructed DEM achieves simulation accuracy within the allowable deviation limits specified by relevant technical standards. This method is characterized by low cost, high efficiency, and strong applicability, providing a practical new approach and technical support for the numerical simulation of large-scale rivers with scarce data, such as the Tarim River.
大尺度河道 / 数字高程模型(DEM) / Copernicus DEM 30 / 塔里木河干流
large-scale river / digital elevation model (DEM) / Copernicus DEM 30 / main stream of Tarim River
| [1] |
苏程佳, 陈晓宏. 一种基于Google Earth影像的河道地形数据提取方法[J]. 水资源保护, 2019, 35(2): 80-84.
(
|
| [2] |
胡朝阳, 王新强, 叶丽清, 等. 利用MIKE21前后处理模块实现河床冲淤可视化分析[J]. 水电能源科学, 2015, 33(11): 143-146.
(
|
| [3] |
张文婷, 李祉璇, 张行南, 等. 基于DEM的河道断面构造改进方法及洪水演进精度评估[J]. 南水北调与水利科技(中英文), 2022, 20(3): 563-572.
(
|
| [4] |
|
| [5] |
|
| [6] |
孙中强, 谷硕, 王厚杰, 等. 一种基于断面高程数据的河道冲淤计算方法:曲线正交网格地形法[J]. 中国海洋大学学报(自然科学版), 2018, 48(9):90-99.
(
|
| [7] |
赵科锋, 曹慧群, 罗平安. 河道二维数值模拟地形数据前处理精细化研究[J]. 长江科学院院报, 2021, 38(3):16-19.
数值模型是分析河流水动力水质演变规律的重要方法。通过对MIKE21模型边界等间距取节点,并进行编码,实现实测河道断面地形沿河道弯曲方向的精细化插值,为河道数值建模提供优化后的地形数据。为了说明此方法的优越性,选取荆门城区某河流部分河段作为研究对象,利用Python完成程序的编写,对实测河道断面数据进行精细化插值运算,并对比地形精细化处理前后与实际地形的差别。结论认为,经过对6组断面插值运算,河道地形得以细化,减少了模型运算错误率,使得模型更趋近于真实情况。
(
|
| [8] |
|
| [9] |
孔乔, 韩露, 刘兴坡, 等. 一种采用SRTMDEM重构河床数字地形的方法[J]. 地球信息科学学报, 2021, 23(3):385-394.
数字高程模型（DEM）包含的信息常作为重要的水文水动力研究基础数据,但是由于公共源DEM数据精度不能完整表达河床地形,所以无法应用于河流泛洪分析等研究工作。因此,本研究开展了基于DEM数据构建数字河床的工作,首先对提取的纵向河网高程数据引入了强局部加权回归算法进行平滑处理以消除畸点;然后以河面要素文件为掩膜,采用反距离加权方法进行横向的空间插值得到3D数字化河面;最后结合河面高程与河深数据完成河床的整体构建。本研究以SRTM DEM为数据地形基础,以位于我国吉林省永吉县境内的温德河为研究实例,对其下游局部区域内的河段重构了矩形、梯形和V型 3种河床断面数字地形。为评价重构数据的合理性,分别对所得到的数字河面进行等值线分析;对河床构建前后DEM水文分析提取的河网进行河道偏移量计算;对构建的地形数据进行河流水动力模拟应用分析,结果表明：① 应用反距离加权插值可以很好地实现河网高程向河面高程的横向延展;② 基于河床重构DEM数据的水文分析误差得到了很好地消除;③ 本研究的重构数据在河流分析中具有较好地应用性和可靠性。
(
|
| [10] |
Multidimensional hydrodynamic modelling becomes tricky when lacking the bathymetric data representing the continuous underwater riverbed surface. Light detection and ranging (LiDAR)-based and radar-based digital elevation models (DEMs) are often used to build the high-accuracy floodplain topography, while in most cases the submerged riverbed could not be detected because both radar and LiDAR operate at wavelengths that cannot penetrate the water. Data from other sources is therefore required to establish the riverbed topography. The inundated river channel is often surveyed with an echo sounder to obtain discrete cross-section data. In this context, an improved algorithm based on the classic flow-oriented coordinates transformation is proposed to generate the riverbed topography using surveyed cross-sections. The dimensionless channel width (DCW) processing method is developed within the algorithm to largely increase the prediction accuracy, especially for the meandering reaches. The generated riverbed topography can be merged with the floodplain DEM to create an integrated DEM for 2D and 3D hydrodynamic simulations. Two case studies are carried out: a benchmark test in the Baxter River, United States, with carefully surveyed channel–floodplain topographic data to validate the algorithm, and a 3D hydrodynamic modelling-based application in Three Gorges Reservoir (TGR) area, China. Results from the benchmark case demonstrate very good consistency between the created topography and the surveyed data with root mean square error (RMSE) = 0.17 m and the interpolation accuracy was increased by 55% compared to the traditional method without DCW processing. 3D hydrodynamic modelling results match the observed field data well, indicating that the generated DEM of the TGR area was good enough not only to predict water depths along the tributary, but also to allow the hydrodynamic model to capture the typical features of the complex density currents caused by both the topography of the tributary estuary and the operation rules of TGR.
|
| [11] |
|
| [12] |
王延贵, 胡春宏, 周文浩, 等. 塔里木河干流河道的河型及其成因[J]. 泥沙研究, 2002, 27(6):19-25.
(
|
| [13] |
袁寄望, 宗全利, 冯博. 塔里木河河槽形态调整特点及影响因素[J]. 水利水电科技进展, 2019, 39(3): 24-32.
(
|
| [14] |
王延贵, 胡春宏, 周文浩, 等. 塔里木河干流的河床演变特点[J]. 水利学报, 2003, 34(12):27-33.
(
|
| [15] |
Flood hazard and risk analysis in developing countries is a difficult task due to the absence or scarce availability of flow data and digital elevation models (DEMs) with the necessary quality. Up to eight DEMs (ALOS Palsar, Aster GDEM, Bare Earth DEM, SRTM DEM, Merit DEM, TanDEM-X DEM, NASA DEM, and Copernicus DEM) of different data acquisition, spatial resolution, and data processing were used to reconstruct the January 2015 flood event. The systematic flow rate record from the Mocuba city gauge station as well as international aid organisms and field data were used to define both the return period peak flows in years for different flood frequencies (Tyear) and the January 2015 flooding event peak flow. Both visual and statistical analysis of flow depth values at control point locations give us a measure of the different hydraulic modelling performance. The results related to the Copernicus DEM, both in visual and statistical approach, show a clear improvement over the results of the other free global DEMs. Under the assumption that Copernicus DEM provides the best results, a flood hazard analysis was carried out, its results being in agreement with previous data of the effects of the January 2015 flooding event in the Mocuba District. All these results highlight the step forward that Copernicus DEM represents for flood hazard analysis in developing countries, along with the use of so-called “citizen science” in the form of flooding evidence field data acquisition.
|
| [16] |
Quantitative geomorphic research depends on accurate topographic data often collected via remote sensing. Lidar, and photogrammetric methods like structure-from-motion, provide the highest quality data for generating digital elevation models (DEMs). Unfortunately, these data are restricted to relatively small areas, and may be expensive or time-consuming to collect. Global and near-global DEMs with 1 arcsec (∼30 m) ground sampling from spaceborne radar and optical sensors offer an alternative gridded, continuous surface at the cost of resolution and accuracy. Accuracy is typically defined with respect to external datasets, often, but not always, in the form of point or profile measurements from sources like differential Global Navigation Satellite System (GNSS), spaceborne lidar (e.g., ICESat), and other geodetic measurements. Vertical point or profile accuracy metrics can miss the pixel-to-pixel variability (sometimes called DEM noise) that is unrelated to true topographic signal, but rather sensor-, orbital-, and/or processing-related artifacts. This is most concerning in selecting a DEM for geomorphic analysis, as this variability can affect derivatives of elevation (e.g., slope and curvature) and impact flow routing. We use (near) global DEMs at 1 arcsec resolution (SRTM, ASTER, ALOS, TanDEM-X, and the recently released Copernicus) and develop new internal accuracy metrics to assess inter-pixel variability without reference data. Our study area is in the arid, steep Central Andes, and is nearly vegetation-free, creating ideal conditions for remote sensing of the bare-earth surface. We use a novel hillshade-filtering approach to detrend long-wavelength topographic signals and accentuate short-wavelength variability. Fourier transformations of the spatial signal to the frequency domain allows us to quantify: 1) artifacts in the un-projected 1 arcsec DEMs at wavelengths greater than the Nyquist (twice the nominal resolution, so > 2 arcsec); and 2) the relative variance of adjacent pixels in DEMs resampled to 30-m resolution (UTM projected). We translate results into their impact on hillslope and channel slope calculations, and we highlight the quality of the five DEMs. We find that the Copernicus DEM, which is based on a carefully edited commercial version of the TanDEM-X, provides the highest quality landscape representation, and should become the preferred DEM for topographic analysis in areas without sufficient coverage of higher-quality local DEMs.
|
| [17] |
Many biochemical processes and dynamics are strongly controlled by terrain topography, making digital elevation models (DEM) a fundamental dataset for a range of applications. This study investigates the quality of four pan-Arctic DEMs (Arctic DEM, ASTER DEM, ALOS DEM and Copernicus DEM) within the Kalix River watershed in northern Sweden, with the aim of informing users about the quality when comparing these DEMs. The quality assessment focuses on both the vertical accuracy of the DEMs and their abilities to model two fundamental elevation derivatives, including topographic wetness index (TWI) and landform classification. Our results show that the vertical accuracy is relatively high for Arctic DEM, ALOS and Copernicus and in our study area was slightly better than those reported in official validation results. Vertical errors are mainly caused by tree cover characteristics and terrain slope. On the other hand, the high vertical accuracy does not translate directly into high quality elevation derivatives, such as TWI and landform classes, as shown by the large errors in TWI and landform classification for all four candidate DEMs. Copernicus produced elevation derivatives with results most similar to those from the reference DEM, but the errors are still relatively high, with large underestimation of TWI in land cover classes with a high likelihood of being wet. Overall, the Copernicus DEM produced the most accurate elevation derivatives, followed by slightly lower accuracies from Arctic DEM and ALOS, and the least accurate being ASTER.
|
| [18] |
|
| [19] |
The changes in catchments can be analyzed through the generation of DEM, which is important as input data in hydrologic modeling. This study aims to analyze the effect of anthropogenic activities on hydrological studies based on DEM comparison and GIUH hydrographs. The four DEM datasets (SRTM, ALOS, Copernicus, Sentinel-1) were compared to the topographic map of Makkah City and GPS data in order to assess the quality of the DEM elevation. The GIS Arc Hydro toolbox was used to extract morphometric and Horton–Strahler ratio characteristics to generate a GIUH hydrograph of the catchments of Wadi Nouman and Wadi Ibrahim inside Makkah City. Based on the DEM comparison, Copernicus and SRTM have the highest accuracy, with R2 = 0.9788 and 0.9765, and the lowest RMSE, 3.89 m and 4.23 m, respectively. ALOS and Sentinel-1 have the lowest R2, 0.9687 and 0.9028, and the highest RMSE, 4.27 m and 6.31 m, respectively. GIUH Copernicus DEM on Wadi Nouman has a higher qp and lower tp (0.21 1/h and 2.66 h) than SRTM (0.20 1/h and 2.75 h), respectively. On Wadi Ibrahim, the SRTM has a greater qp and lower tp than Copernicus due to the wadi having two shapes. Based on the anthropogenic effect, the stream network in the mountain area is quite similar for SRTM and Copernicus due to the dominant influence of the mountainous relief and relatively inconsequential influence of anthropogenic activities and DEM noise. In the urban area, the variation of the stream network is high due to differing DEM noise and significant anthropogenic activities such as urban redevelopment. The Copernicus DEM has the best performance of the others, with high accuracy, less RMSE, and stream flow direction following the recent condition.
|
| [20] |
|
| [21] |
|
| [22] |
The interferometric synthetic aperture radar (InSAR) data set, acquired by the TanDEM-X (TerraSAR-X add-on for Digital Elevation Measurement) mission (TDM), represents a unique data source to derive geo-information products at a global scale. The complete Earth’s landmasses have been surveyed at least twice during the mission bistatic operation, which started at the end of 2010. Examples of the delivered global products are the TanDEM-X digital elevation model (DEM) (at a final independent posting of 12 m × 12 m) or the TanDEM-X global Forest/Non-Forest (FNF) map. The need for a reliable water product from TanDEM-X data was dictated by the limited accuracy and difficulty of use of the TDX Water Indication Mask (WAM), delivered as by-product of the global DEM, which jeopardizes its use for scientific applications, as well. Similarly as it has been done for the generation of the FNF map; in this work, we utilize the global data set of TanDEM-X quicklook images at 50 m × 50 m resolution, acquired between 2011 and 2016, to derive a new global water body layer (WBL), covering a range from −60∘ to +90∘ latitudes. The bistatic interferometric coherence is used as the primary input feature for performing water detection. We classify water surfaces in single TanDEM-X images, by considering the system’s geometric configuration and exploiting a watershed-based segmentation algorithm. Subsequently, single overlapping acquisitions are mosaicked together in a two-step logically weighting process to derive the global TDM WBL product, which comprises a binary averaged water/non-water layer as well as a permanent/temporary water indication layer. The accuracy of the new TDM WBL has been assessed over Europe, through a comparison with the Copernicus water and wetness layer, provided by the European Space Agency (ESA), at a 20 m × 20 m resolution. The F-score ranges from 83%, when considering all geocells (of 1∘ latitudes × 1∘ longitudes) over Europe, up to 93%, when considering only the geocells with a water content higher than 1%. At global scale, the quality of the product has been evaluated, by intercomparison, with other existing global water maps, resulting in an overall agreement that often exceeds 85% (F-score) when the content in the geocell is higher than 1%. The global TDM WBL presented in this study will be made available to the scientific community for free download and usage.
|
| [23] |
|
| [24] |
|
| [25] |
|
| [26] |
刘佳嘉, 周祖昊, 贾仰文, 等. 河网提取算法中河道断裂现象原因剖析及算法改进[J]. 中国水利水电科学研究院学报, 2013, 11(3): 195-200.
(
|
| [27] |
杨发相. 塔里木河流域地貌过程引起的环境变化[C]// 中国地理学会地貌与第四纪专业委员会.地貌·环境·发展:2004丹霞山会议文集. 北京: 中国环境科学出版社, 2004:62-66.
(
|
| [28] |
刘伟. 塔里木河上游段河道形态初步分析与对策[J]. 水利建设与管理, 2007, 27(4):70-72.
(
|
| [29] |
周庆庆. 水力插板透水丁坝缓流促淤效果及主要影响因素试验研究[D]. 乌鲁木齐: 新疆农业大学, 2015.
(
|
| [30] |
黄志涛. 某江段冲淤演变及防护方案数值模拟研究[D]. 哈尔滨: 哈尔滨工程大学, 2016.
(
|
| [31] |
中华人民共和国水利部. 中华人民共和国水文年鉴: 第10卷内陆河湖水文资料[R]. 北京: 中华人民共和国水利部, 2011.
(Ministry of Water Resources of the People’s Republic of China. Hydrological Yearbook of the People’s Republic of China: Volume 10 Hydrological Data of Inland Rivers and Lakes[R]. Beijing: Ministry of Water Resources of the People’s Republic of China, 2011.(in Chinese))
|
| [32] |
中工武大设计集团有限公司. 新疆维吾尔自治区塔里木河干流沙雅段防洪治理项目可行性研究报告[R]. 武汉: 中工武大设计集团有限公司, 2022.
(China Industry Wuda Design Group Co., Ltd. Feasibility Study Report on Flood Control Project of Shaya Section of the Main Stream of Tarim River in Xinjiang Uygur Autonomous Region[R]. Wuhan: China Industry Wuda Design Group Co., Ltd., 2022.(in Chinese))
|
| [33] |
郭庆超, 王新平, 周军, 等. 塔里木河干流近期治理后水沙运动规律研究[M]. 北京: 中国水利水电出版社, 2016.
(
|
| [34] |
JTS/T 231-4—2018,内河航道与港口水流泥沙模拟技术规程[S]. 北京: 人民交通出版社, 2018.
(JTS/T 231-4—2018, Technical Regulation of Modelling for Flow and Sediment in Inland Waterway and Harbour[S]. Beijing: China Communications Press, 2018.(in Chinese))
|
/
| 〈 |
|
〉 |
