基于CryoSat-2的东南极PANDA断面考察沿线DEM制作及精度分析

基于卫星测高数据,国际上先后发布了ICESat DEM、Bamber DEM等全南极DEM.相对于传统高度计,大轨道倾角、长重访周期的设计使CryoSat-2具有更大的数据覆盖范围以及更加密集的轨道覆盖;同时,SARIn模式的启用也提高了CryoSat-2对于南极边缘区域的监测能力.针对利用CryoSat-2数据提取DEM问题,对CroySat-2的轨迹覆盖特征、数据滤波方法、最优内插参数选取以及DEM精度验证等方面进行了探讨,并利用克里金插值法生成了东南极PANDA断面考察沿线1 km分辨率的DEM.结果表明：通过与ICESat数据对比,发现PANDA断面考察沿线DEM整体高程精度约为（1.57±3.30）m;但是,局部区域高程精度分析表明,DEM精度具有不均一性,随着坡度的增加DEM高程精度逐渐下降,高程稳定性也同时下降.     

Monitoring surface elevation changes in Jharia coalfield,India using synthetic aperture radar interferometry

This paper investigates surface elevation changes that occurred during 1996–2004 in the Jharia coalfield through the digital elevation model (DEM) generated using synthetic aperture radar interferometry (InSAR) using ERS-1/2 (European Remote Sensing Satellite) tandem and RADARSAT-1 data. The comparison of elevation values derived from the InSAR DEM and topographic height data shows a bias of 23.08 m with root-mean-square error of ±2.31 m (5.8 %). The accuracy of the DEM was investigated by comparing the elevation profiles with the digitized elevation contour data at four different locations. The profile comparison shows a mean bias of 22.68 m. Local topography shows changes in elevation up to ±40.00 m due to mining activities on the 8-year time period. The results of InSAR-derived heights and topographic heights were comparable and well-matched except at a few locations where topographic data were unavailable. DEM generated using InSAR due to its high spatial details is ideal for the detection and estimation of surface elevation changes in mining areas.     

Geostatistical Mapping with Continuous Moving Neighborhood

An issue that often arises in such GIS applications as digital elevation modeling (DEM) is how to create a continuous surface using a limited number of point observations. In hydrological applications, such as estimating drainage areas, direction of water flow is easier to detect from a smooth DEM than from a grid created using standard interpolation programs. Another reason for continuous mapping is esthetic; like a picture, a map should be visually appealing, and for some GIS users this is more important than map accuracy. There are many methods for local smoothing. Spline algorithms are usually used to create a continuous map, because they minimize curvature of the surface. Geostatistical models are commonly used approaches to spatial prediction and mapping in many scientific disciplines, but classical kriging models produce noncontinuous surfaces when local neighborhood is used. This motivated us to develop a continuous version of kriging. We propose a modification of kriging that produces continuous prediction and prediction standard error surfaces. The idea is to modify kriging systems so that data outside a specified distance from the prediction location have zero weights. We discuss simple kriging and conditional geostatistical simulation, models that essentially use information about mean value or trend surface. We also discuss how to modify ordinary and universal kriging models to produce continuous predictions, and limitations using the proposed models.     

Glacier volume estimation from ice-thickness data,applied to the Muz Taw glacier,Sawir Mountains,China

In this work, ground penetrating radar (GPR) data, global position system (GPS) data, Landsat remote sensing images and digital elevation model (DEM) SRTM were used to study the Muz Taw glacier in the Sawir Mountains. By a variety of interpolation methods, spatial pattern of thickness of the Muz Taw glacier was obtained. The results show that there are obvious differences between the available GPR transverse profiles. The lower transverse profile has typical “U”-type topographic feature, the maximum glacier thickness is up to 116.4 m; glacier thickness distribution in the higher profile is relatively flat, with an average about 70–90 m. The longitudinal profile forms a ladder-like distribution, which in the middle of the ice bedrock topography has obvious depression. The average longitudinal ice thickness is about 80.89 m; the maximum ice thickness is up to 122.61 m. In addition, the preliminary estimation of the average thickness and total terminus ice volume were approximately 60.5 m and 0.195 km³, respectively.     

Modelling the geoid and sea-surface topography in coastal areas

This paper looks at the relation between the time-averaged level of the sea surface and a gravimertic geoid, as determined in coastal areas. Measurements in local regions can now be accurate enough to demonstrate that the geoid and mean sea level are not even parallel to each other, let alone identical. The accuracy and pattern structure of surface gravity data in some shelf seas is comparable with those on land, so that a marine geoid can be derived from surface data without using satellite altimetry. The geodetic objective is then to combine the two to determine sea surface topography. In principle, gravimetric studies provide the absolute datum so that local oceanographic models on the shelf can be combined with sea surface topography models related to the global ocean circulation. In contrast, sea surface topography information near deep ocean coasts must come from external sources and satellite altimetry used to give the gravity data needed to offset the less good coverage by ship-borne gravimetry.Marine Bouguer anomalies enable two specific problems of gravity anomaly patterns near the continent ocean transition to be overcome. The necessary extension of Stokes' condensation reduction is developed and illustrated along a north-south profile from the Mediterranean across the Cote d'Azur. The effect on gravity of deep ocean water introduces a geoid correction in the form of a dipolar ridge whose amplitude at the shore is about 11 cm. In addition to geostrophic currents, a semi-quantitative model for the thermohaline effects on sea surface topography is discussed in relation to sea level differences between the Atlantic and Mediterranean.In considering appropriate algorithms for local geoid computation, Kirby's Iterative Fourier Combination routine for combining altimetry and surface gravity is extended to account for global sea surface topography. The impact of very fast spherical harmonic analysis algorithms is discussed and a simple physical model is given which explains the short coherence lengths found for the global gravity field. This necessary assumption for any local geoid computation was hitherto purely empirical.Finally, the use of land data such as tide gauges, ellipsoidal heights from GPS, and orthometric heights from first order levelling are reviewed as ways of corroborating geodetic estimates of sea surface topography and its relation to levelling datums. Successful examples are given from southern England.     

基于遥感专题解译的白洋淀湿地底部DEM构建

基于多时相Landsat遥感影像,利用适合白洋淀湿地信息提取的综合提取方法提取湿地范围。将湿地范围视为水体淹没区,淹没区的边界线视为对应地表水位高程的等高线,并根据遥感影像、高程控制点等对淀区内受人为影响的局部区域进行校正,插值生成白洋淀湿地底部的数字高程模型(DEM)。选用遥感影像和地面高程控制点对构建的数字高程模型进行验证,遥感影像验证精度在80%以上,地面控制点验证误差在±0.5 m以内的点达到80%。这种由一系列遥感影像提取等高线生成数字高程模型的方法可以弥补白洋淀湿地内数字高程信息不足的缺陷,对于提高白洋淀湿地的水均衡计算和构建湖泊与地下水耦合模型的精度有着重要意义。     

基于ICESat数据的南极冰盖DEM插值方法比较及精度分析

南极数字高程模型(DEM)是从事南极地学和环境变化研究的基础. 内插是建立数字高程模型的重要技术点, 插值方法有多种, 根据不同的适用情况, 不同的插值方法各有优劣. 利用克里格、 距离反权、 三角网剖分、 最小曲率以及移动平均5种插值方法分别建立南极冰盖小范围区域的DEM, 通过抽取部分观测数据作为验证值对各插值方法进行了比较. 结果表明: 克里格插值方法的可靠性最好, 稳定性最高. 然后, 利用克里格插值方法, 基于ICESat测高卫星的GLA12数据建立了南极冰盖的DEM. 由于南极大陆实测数据有限, 缺乏对DEM的检核. 为了分析所建DEM的可靠性, 利用中国南极内陆冰盖考察所采集的GPS实测数据, 对所建立的DEM进行了验证分析. 结果显示, DEM在坡度较缓的南极内陆冰盖区域精度较高, 符合度在3 m以内; 距离卫星轨道越近的区域精度越高, 可达到1 m 以内. 在坡度较大, 高程变化较为显著的区域如沿海地区, 精度较低, 差距最大的点超过40 m.     

A modified HASM algorithm and its application in DEM construction

In many spatial interpolation fields, high accuracy surface modeling (HASM) has yielded better accuracy than classical interpolation methods. The Gaussian equation is the core of the HASM algorithm; The current version of the HASM method builds the Gaussian equation in Cartesian coordinates and, computes the two partial derivatives of the surface in the horizontal and vertical directions for each grid. In this paper, a modified HASM method is proposed that integrates flow paths to improve the original HASM methodology. The modified HASM approach involves two steps. The first step generates an initial DEM, which is used to compute the flow path. Then, the second step is conducted based on scatter points and the flow direction. The output from this step is better than the initial DEM. First, we used a theoretical mathematical surface to validate the correctness of the modified model. Then, we chose a small study area where the topography is affected by hydrological erosion for analysis. The test results showed that the modified HASM method constructed a DEM with low MAE and RMSE values compared to those of traditional methods, and it more accurately characterized topographic features. Finally, a relatively gently sloping area was selected to validate that the applicability of the new method in other areas.     

基于ICESat/GLAS数据的中国典型区域SRTM与ASTER GDEM高程精度评价

航天飞机雷达地形测绘(shuttle radar topography mission, SRTM)和先进星载热发射和反射辐射成像仪全球数字高程模型(advanced spaceborne thermal emission and reflection radiometer global digital elevation model, ASTER GDEM)提供了全球覆盖面积最广的数字高程模型(digital elevation model, DEM)数据, 但其高程精度还未得到充分验证, 传统地面测量方法很难适用于验证大面积范围的DEM精度.以冰、云和陆地高程卫星/地学激光测高系统(ICESat/GLAS)高程数据为参考, 综合利用地理信息系统(geographic information system, GIS)空间分析、三维可视化与统计分析方法, 对中国典型低海拔沿海平原地区和高海拔山地的两种DEM数据高程精度进行了对比分析.结果表明, 高程值小于20m的低海拔地区, SRTM高程精度达到2.39m, ASTER GDEM的精度达到4.83m, 均远远高于这两种数据的标称精度; 而在西南山地, 这两种DEM的精度大约为20m, 与标称精度相当.最后, 建立了ICESat/GLAS与SRTM和ASTER GDEM的一元线性回归模型, 该模型具有较高的拟合度和显著线性关系, 可用于改善这两种DEM的高程精度.      

四种南极数字高程模型的精度比较与分析

南极数字高程模型（DEM）是南极冰盖研究的基础数据,目前国际通用的全南极DEM数据主要有JLB97 DEM、RAMPv2 DEM、ICESat DEM以及Bamber 1km DEM. 利用DEM对DEM验证的方式对四种DEM的精度进行比较分析. 结果表明：Bamber 1km DEM和ICESat DEM之间的高程差异最小,平均高程差小于1.8 m,二者均有较高的可靠性. RAMPv2 DEM与Bamber 1km DEM的高程差大于1.9 m,在81.5° S以南和坡度较大的区域,高程差异更为明显,高程可靠性较低. JLB97 DEM与上述三种DEM的偏差超过10 m,高程可靠性最低.     

东昆仑山煤矿冰川雷达测厚及冰储量估算

基于2015年5月东昆仑山煤矿冰川冰厚测量资料,结合2015年Landsat8 OLI影像,利用Ordinary Kriging插值方法对冰川非测厚区进行插值计算,绘制了该冰川厚度等值线并对该冰川冰储量进行了估算。2015年煤矿冰川最大厚度为87 m,位于海拔4 952 m主流线附近,平均厚度为25.77 m,冰储量为0.0242 km³。利用煤矿冰川冰面地形图与冰厚度分布图,获得该冰川冰床地形图。结果显示,冰川上宽下窄沿山谷分布,冰床地形复杂,在冰厚较大区域,地形呈近“V”字分布,显示了冰斗冰川的形貌特征。     

Accuracy Analysis of Digital Elevation Model Relating to Spatial Resolution and Terrain Slope by Bilinear Interpolation

This study investigates the accuracy analysis of the digital elevation model (DEM) with respect to the following two major factors that strongly affect the interpolated accuracy: (1) spatial resolution of a DEM and (2) terrain slope. Unlike existing studies based mainly on a simulation approach, this research first provides an analytical approach in order to build the relationship between the interpolated DEM accuracy and its influencing factors. The bi-linear interpolation model was adopted to produce this analytic model formalized as inequalities. Then, our analytic models were verified and further rectified by means of experimental studies in order to derive a practical formula for estimating the DEM accuracy together with an optimization model for calculating the required resolution when a prescribed upper bound to the DEM accuracy is given. Moreover, this analytic approach can cope with either a grid-based DEM or a randomly scattered scenario whose efficacies have been validated by the experiments using both synthetic and realistic data sets. In particular, these findings first establish the rules for directly correlating the horizontal resolution of DEM data with vertical accuracy.     

基于合成孔径雷达技术及DEM的公格尔山冰川动力特征分析

冰川流速是冰川动力状况的重要标志,利用合成孔径雷达技术能快速获得大范围冰川的表面流速.利用日本高级陆地观测卫星(ALOS)相控阵型合成孔径雷达L波段(PALSAR)及欧洲太空局的ENVISAT/ASAR数据的特征匹配方法获得帕米尔高原公格尔山区冰川表面流速,并结合合成孔径雷达干涉相干与不同时期数字高程模型对公格尔山区典型冰川动力进行分析,获得研究区不同时间基线冰川表面相干性、表面流速以及基于不同时相DEM的典型冰川表面高程变化信息.结果表明:30 a来克拉牙依拉克冰川表面高程下降了(15±12.1)m,表碛区域近期运动速度变化不大;其木干冰川平均表面高程几乎无变化,但2007-2011年夏季表面流速明显减缓,靠近末端附近部分区域可能已经演化为非活动区;姜满加尔冰川位于西风带的迎风坡,积累区面积大,冰川流速较快,无表碛覆盖,但表面高程仍下降了(8.8±12.7)m.编号为5Y663D0009的冰川冰舌表碛覆盖区可能已经演化为非活动区,30 a来表面高程下降(8.6±12.0)m.综合分析表明,冰川规模特别是积累区面积的大小及所处位置、地形对冰川演化具有重要影响.     

Evaluation of intermediate TanDEM-X digital elevation data products over Tasmania using other digital elevation models and accurate heights from the Australian National Gravity Database

The successful operation of the TanDEM-X satellite mission is the start of a new era of globally consistent and accurate digital elevation data for planet Earth. In this work available 12 m-resolution intermediate TanDEM-X products (DEM: digital elevation model; HEM: height error map; COV: coverage map; WAM: water indication mask) are evaluated over Tasmania. Elevations from the TanDEM-X intermediate digital elevation model (IDEM) are compared with (a) other global DEMs (30 m-resolution SRTM1 USGS v3 and 30 m-resolution Advanced Spaceborne Thermal Emission Reflectometer (ASTER GDEM2), (b) the local 25 m-resolution DEM made available by Tasmanian environmental authority (DPIPWE), and (c) over 15 000 accurate ground-control-points (GCPs) from the Australian National Gravity Database (ANGD). The comparison with ASTER and SRTM over the area of Tasmania involves over 500 million valid TanDEM-X IDEM elevations. The root-mean-square (RMS) of 8.8 m indicates a reasonable to good agreement of TanDEM-X IDEM and SRTM, while ASTER shows almost twice the disagreement in terms of RMS (～16.5 m). Both, ASTER and SRTM show a (mean) offset of –1.9 m and –2.3 m w.r.t. TanDEM-X IDEM, respectively. By comparisons with GCPs, we find that SRTM and ASTER overestimate the terrain height. The comparison with the AGND GCPs also allows an estimate of the absolute accuracy of the IDEM, which is found to be superior to that of SRTM or ASTER. The RMS error of 6.6 m shows that the IDEM is close to the officially denoted 4 m absolute vertical accuracy considering that the GCPs are not error free. The height error map information layer is found to a suitable first indicator of the (local) accuracy of the IDEM in a relative sense. However, we find that the HEM tends to underestimate observed differences to the GCPs. Terrain-type analyses reveal that the TanDEM-X IDEM is a very consistent elevation database over Tasmania. In conclusion, our study demonstrates that the new TanDEM-X elevation data sets provide improved high-resolution terrain information over Tasmania and beyond.     

Continental/Cordilleran ice interactions: a dominant cause of westward super-elevation of the last glacial maximum continental ice limit in southwestern Alberta, Canada

In southwestern Alberta, Canada, a westward-rising last-glacial-maximum continental ice limit has been identified. This limit is defined by the upper elevation of Canadian Shield erratics deposited by last-glacial-maximum continental ice along the flanks of prominent ridges and buttes within the region. The interpolation between ice-limit data points has produced two distinct slope profiles: 2.9 m/km to the east, and 4.2 m/km to the west of Mokowan Butte. Three hypotheses are proposed to explain this westward rise of the last-glacial-maximum continental ice limit: (1) regional tectonic uplift, (2) glacio-isostatic uplift, and (3) continental ice-flow convergence due to topographic obstacles and interaction with montane ice. Inferred long-term rates of tectonic uplift and glacio-isostatic modelling show that these two mechanisms account for less than 25% of the observed absolute elevation increase of the limit between the Del Bonita uplands and Cloudy Ridge in southwestern Alberta. The remaining rise in elevation of the continental ice-sheet margin in this region is thought to result from continental ice-flow convergence due to the combined effects of the regional topography and interaction with montane glaciers to the west. The steeper rise in the former continental ice surface west of Mokowan Butte can be explained by the topographic obstruction and interaction with montane glaciers in the area of the Rocky Mountain front.     

念青唐古拉山扎当冰川冰储量估算及冰下地形特征分析

冰川体积估算对水资源以及冰川变化研究具有重要的意义. 但是实测的冰川厚度数据十分稀少,限制了冰川体积的估算. 2011年5月对念青唐古拉山北坡扎当冰川进行了雷达测厚工作,获取了该冰川的厚度分布状况. 基于该冰川的厚度数据,测量点的GPS数据,1970年的地形图和2010年Landsat TM影像,在ArcGIS技术的支持下,采用简单Kriging插值方法对冰川非测厚区域的厚度进行了插值计算,绘制出了冰川厚度等值线图并估算了冰川的冰储量. 结果表明：冰川最大厚度出现于海拔约5 748 m靠近主流线的位置,最大冰厚度为108 m,冰川平均厚度为38.1 m,2010年冰川面积为1.73 km²,扎当冰川的冰储量为0.066 km³. 将扎当冰川表面DEM与冰川厚度分布图相结合,绘制出了该冰川的冰床地形图. 结果显示,在冰川厚度大的区域,冰床地形呈现近V字形分布,这与其相对平缓的冰面地形形成明显对比;同时,在冰表地形较陡区域,冰川厚度不大,冰床地形呈现U形分布.     

Determining the optimum cell size of digital elevation model for hydrologic application

Scale is one of the most important but unsolved issues in various scientific disciplines that deal with spatial data. The arbitrary choice of grid cell size for contour interpolated digital elevation models (DEM) is one of the major sources of uncertainty in the hydrologic modelling process. In this paper, an attempt was made to identify methods for determining an optimum cell size for a contour interpolated DEM in prior to hydrologic modelling. Twenty-meter interval contour lines were used to generate DEMs of five different resolutions, viz., 30, 45, 60, 75, and 90 m using TOPOGRID algorithm. The obtained DEMs were explored for their intrinsic quality using four different methods, i.e., sink analysis, fractal dimension of derived stream network, entropy measurement and semivariogram modelling. These methods were applied to determine the level artifacts (interpolation error) in DEM surface as well as derived stream network, spatial information content and spatial variability respectively. The results indicated that a 90 m cell size is sufficient to capture the terrain variability for subsequent hydrologic modelling in the study area. The significance of this research work is that it provides methods which DEM users can apply to select an appropriate DEM cell size in prior to detailed hydrologic modelling.     

Incorporation of ice sheet models into an Earth system model: Focus on methodology of coupling

Elaboration of a modern Earth system model (ESM) requires incorporation of ice sheet dynamics. Coupling of an ice sheet model (ICM) to an AOGCM is complicated by essential differences in spatial and temporal scales of cryospheric, atmospheric and oceanic components. To overcome this difficulty, we apply two different approaches for the incorporation of ice sheets into an ESM. Coupling of the Antarctic ice sheet model (AISM) to the AOGCM is accomplished via using procedures of resampling, interpolation and assigning to the AISM grid points annually averaged meanings of air surface temperature and precipitation fields generated by the AOGCM. Surface melting, which takes place mainly on the margins of the Antarctic peninsula and on ice shelves fringing the continent, is currently ignored. AISM returns anomalies of surface topography back to the AOGCM. To couple the Greenland ice sheet model (GrISM) to the AOGCM, we use a simple buffer energy- and water-balance model (EWBM-G) to account for orographically-driven precipitation and other sub-grid AOGCM-generated quantities. The output of the EWBM-G consists of surface mass balance and air surface temperature to force the GrISM, and freshwater run-off to force thermohaline circulation in the oceanic block of the AOGCM. Because of a rather complex coupling procedure of GrIS compared to AIS, the paper mostly focuses on Greenland.