Spot Revisited: Accuracy Assessment, Dem Generation and Validation from Stereo Spot 5 Hrg Images

SPOT 5 HRG Level 1A and 1B stereo scenes covering Zonguldak testfield in north-west Turkey have been analysed. They comprise the left and right image components with base to height ratio of 0·54. The pixel size on the ground is 5 m. The bundle orientation was executed by the PCI Geomatica V9.1.4 software package and resulted in 3D geopositioning to sub-pixel accuracies in each axis provided that at least six control points were used in the computation. Root mean square error (rmse) values and vectors of residual errors for Levels 1A and 1B are similar, even for different control and check point configurations. Based on the scene orientation, Level 1A and 1B digital elevation models (DEMs) of the testfield have been determined by automatic matching and validated by the reference DEM digitised from the 1:25 000 scale topographic maps, interferometric DEMs from Shuttle Radar Topography Mission (SRTM) X- and C-band SAR data and the GPS profiles measured along the main roads in the testfield. Although the accuracies of reference data-sets are too similar to the generated SPOT DEMs, these are the only high quality reference materials available in this area. Sub-pixel height accuracy was indicated by the comparison with profile points. However, they are in favourable locations where matching is always successful, so such a result may give a biased measure of the accuracy of the corresponding DEMs.     

Assessment of Shuttle Radar Topography Mission Elevation Data Based on Topographic Maps in Turkey

Depending on scale, topographic maps depicting the shape of the land surfaces of the Earth are produced from different data sources. National topographic maps at a scale of 1:25 000 (25K maps) produced by General Command of Mapping are used as the base map set in Turkey. This map set, which consists of approximately 5500 sheets, covers the whole country and is produced using photogrammetric methods. Digital Elevation Models (DEMs) created from these maps are also available. Recently, another data source, Synthetic Aperture Radar (SAR) interferometric data, has become more important than those produced by conventional methods. The Shuttle Radar Topography Mission (SRTM) contains elevation data with 3 arc-second resolution and 16 m absolute height error (90 percent confidence level). These data are freely available via the Internet for approximately 80 percent of the Earth's land mass. In this study, SRTM DEM was compared with DEM derived from 25K topographic maps for different parts of Turkey. The study areas, each covering four neighboring 25K maps, and having an area of approximately 600 km², were chosen to represent various terrain characteristics. For the comparison, DEMs created from the 25K maps were obtained and organized as files for each map sheet in vector format, containing the digitized contour lines. From these data, DEMs in the resolution of 3 arc-second were created (25K-DEM), in the same structure as the SRTM DEM, allowing the 25K-DEMs and the SRTM DEM to be compared directly. The results show that the agreement of SRTM DEM to the 25K-DEM is within about 13 m, which is less than the SRTM's targeted error of 16 m. The spatial distribution of the height differences between SRTM-DEM and the 25K-DEM and correlation analysis show that the differences were mainly related to the topography of the test areas. In some areas, local height shifts were determined.     

Orthorectification of IKONOS and impact of different resolution DEM

IKONOS image has been wildly used in city planning, precision agriculture and emergence response. However, the accuracy of IKONOS Geo product is limited due to distortion caused by terrain relief. Orthorectification was performed to remove the distortion and the impact of different DEM on orthorectification were evaluated. 38 ground control points （GCPs） and 25 independent check points （ICPs） were collected. DEMs were generated from 1 ： 10 000 and 1 ： 50 000 topographic maps. Results show that RMS error at the check points is 1. 554 0 m using DEM generated from 1 ： 10 000 topographic map, which can meet the accuracy requirement of IKONOS Precision product （1.9 m RMSE）. While RMS error is 2. 572 4 m using DEM generated from 1 ： 50 000 topographic map.     

A Comparison Between Contour Elevation Data Sources for DEM Creation and Soil Carbon Prediction, Coshocton, Ohio

A raster and vector GIS was created for the North Appalachian Experimental Watershed (NAEW) from legacy (1960) 1:2,400‐scale contour maps. The intent of the study was to use terrain data for the spatial modeling of soil organic carbon. It was hypothesized that DEMs derived from these data would be more accurate and therefore more useful for terrain‐based soil modeling than those from USGS 1:24,000‐scale contour data. Central tasks for this study were to digitally capture the 1:2,400‐scale maps, convert digital contour data sources to raster DEMs at multiple resolutions, and derive terrain attributes. A flexible approach was adopted, using software outside of mainstream GIS sources where scientifically or practically advantageous. Elevation contours and streamlines were converted to raster DEMs using ANUDEM. DEMs ranging in resolution from 0.5–30 m were tested for accuracy against precision carrier‐phase GPS data. The residual standard deviation was 1.68 meters for the USGS DEM and 0.36 meters for the NAEW DEM. The optimal horizontal resolution for the NAEW DEM was 5 m and for the USGS 10 m. Five and 10 m resolution DEMs from both data sources were tested for carbon prediction. Multiple terrain parameters were derived as proxies for surficial processes. Soil samples (n = 184) were collected on four zero‐order watersheds (conventional tillage, no‐till, hay and pasture). Multiple least squares regressions (m.l.s.) were used to predict mass C (kg m^?2, 30 cm depth) from topographic information. Model residuals were not spatially autocorrelated. Statistically significant topographic parameters were attained most consistently from the 5 m NAEW DEM. However, topography was not a strong predictor of carbon for these watersheds, with r² ranging from 0.23 to 0.58.     

Evaluation of digital elevation models for delineation of hydrological response units in a Himalayan watershed

This study reports results from evaluation of the quality of digital elevation model (DEM) from four sources viz. topographic map (1:50,000), Shuttle Radar Topographic Mission (SRTM) (90 m), optical stereo pair from ASTER (15 m) and CARTOSAT (2.5 m) and their use in derivation of hydrological response units (HRUs) in Sitla Rao watershed (North India). The HRUs were derived using water storage capacity and slope to produce surface runoff zones. The DEMs were evaluated on elevation accuracy and representation of morphometric features. The DEM derived from optical stereo pairs (ASTER and CARTOSAT) provided higher vertical accuracies than the SRTM and topographic map-based DEM. The SRTM with a coarse resolution of 90 m provided vertical accuracy but better morphometry compared to topographic map. The HRU maps derived from the fine resolution DEM (ASTER and CARTOSAT) were more detailed but did not provide much advantage for hydrological studies at the scale of Sitla Rao watershed (5800 ha).     

Improved topographic mapping through high-resolution SAR interferometry with atmospheric effect removal

The application of SAR interferometry (InSAR) in topographic mapping is usually limited by geometric/temporal decorrelations and atmospheric effect, particularly in repeat-pass mode. In this paper, to improve the accuracy of topographic mapping with high-resolution InSAR, a new approach to estimate and remove atmospheric effect has been developed. Under the assumptions that there was no ground deformation within a short temporal period and insignificant ionosphere interference on high-frequency radar signals, e.g. X-bands, the approach was focused on the removal of two types of atmospheric effects, namely tropospheric stratification and turbulence. Using an available digital elevation model (DEM) of moderate spatial resolution, e.g. Shuttle Radar Topography Mission (SRTM) DEM, a differential interferogram was firstly produced from the high-resolution InSAR data pair. A linear regression model between phase signal and auxiliary elevation was established to estimate the stratified atmospheric effect from the differential interferogram. Afterwards, a combination of a low-pass and an adaptive filter was employed to separate the turbulent atmospheric effect. After the removal of both types of atmospheric effects in the high-resolution interferogram, the interferometric phase information incorporating local topographic details was obtained and further processed to produce a high-resolution DEM. The feasibility and effectiveness of this approach was validated by an experiment with a tandem-mode X-band COSMO-SkyMed InSAR data pair covering a mountainous area in Northwestern China. By using a standard Chinese national DEM of scale 1:50,000 as the reference, we evaluated the vertical accuracy of InSAR DEM with and without atmospheric effects correction, which shows that after atmospheric signal correction the root-mean-squared error (RMSE) has decreased from 13.6 m to 5.7 m. Overall, from this study a significant improvement to derive topographic maps with high accuracy has been achieved by using the proposed approach.     

THE IGN (FRANCE) TOPOGRAPHIC DATABASE

The French topographic database is a three dimensional database whose content is very near to that of the IGN 1:25 000 scale maps. It has a metric accuracy and is designed to have three aims (maps at scales of 1:25 000 and 1:50 000 and the topographic part of the cadastral map at 1:5000 scale, providing located structured data). The structure, including topology information, is described. The process of data capture is based on photogrammetry and field completion. Experiments concerned with revision are being made in order to define the methodology precisely.     

一种面向DEM多尺度表达的傅里叶能量谱模型

针对数字高程模型（digital elevation model,DEM）数据的多尺度表达问题,根据DEM格网数据在能量谱密度中“低频-高能-大尺度”的对应关系,在化简中关联地形语义特征,构建了DEM数据的多尺度表达模型。实验结果表明,该模型可以实时动态派生不同尺度下的DEM数据,通过等高线放样观察发现,该模型派生的DEM数据满足地形表达、空间认知和制图综合中的“保留主要地形特征、舍弃次要地形特征”的基本原则。与常用的DEM化简方法进行高程值统计以及坡形变化的定量对比分析,结果表明该方法在统计意义与结构意义上都具有较好的效果。     

The Laser Vegetation Imaging Sensor: a medium-altitude,digitisation-only,airborne laser altimeter for mapping vegetation and topography

The Laser Vegetation Imaging Sensor (LVIS) is an airborne, scanning laser altimeter, designed and developed at NASA's Goddard Space Flight Center (GSFC). LVIS operates at altitudes up to 10 km above ground, and is capable of producing a data swath up to 1000 m wide nominally with 25-m wide footprints. The entire time history of the outgoing and return pulses is digitised, allowing unambiguous determination of range and return pulse structure. Combined with aircraft position and attitude knowledge, this instrument produces topographic maps with dm accuracy and vertical height and structure measurements of vegetation. The laser transmitter is a diode-pumped Nd:YAG oscillator producing 1064 nm, 10 ns, 5 mJ pulses at repetition rates up to 500 Hz. LVIS has recently demonstrated its ability to determine topography (including sub-canopy) and vegetation height and structure on flight missions to various forested regions in the US and Central America. The LVIS system is the airborne simulator for the Vegetation Canopy Lidar (VCL) mission (a NASA Earth remote sensing satellite due for launch in year 2000), providing simulated data sets and a platform for instrument proof-of-concept studies. The topography maps and return waveforms produced by LVIS provide Earth scientists with a unique data set allowing studies of topography, hydrology, and vegetation with unmatched accuracy and coverage.     

利用TanDEM-X生成DEM的精度评定

目前,许多学者对TanDEM-X生成DEM开展了一些研究,其研究成果也显示了TanDEM-X生成高精度DEM的可行性。为了验证TanDEM-X/TerraSAR-X干涉生成的DEM能否满足测图要求,需要对其进行精度评价和分析。相对于C波段的ERS、ASAR和L波段的ALOS,X波段的高分辨率TerraSAR影像干涉条纹更密集,解缠更加困难。针对这一问题,本文设计了一种低分辨率SRTM辅助高分辨率的X波段的TerraSAR干涉相位解缠方案,提高了解缠的效率和精度。同时,本文提出了一种基于协方差函数的方法对TDX/TSX DEM进行精度分析和评价。该方法通过对各个距离上的协方差值进行拟合,消除了高程误差异常对InSAR DEM精度评价的影响,可以更加客观真实地反映DEM的精度。实例分析结果表明：采用协方差函数方法来评价DEM的精度是可行的,对于试验研究区域,TDX/TSX干涉生成的DEM总体精度为1.42 m,能够满足1：10 000测图要求,为我国空白地区的测图提供了有利条件。     

Study of dems derived from ERS-1/2 SAR and SRTM data

The DEM of the Bhuj earthquake affected area of 50 x 50 km was generated using the ERS-1/2 SAR tandem data (May 15—16,1996). Region growing algorithm coupled with path following approach was used for phase unwrapping. Phase to height conversion was done using D-GPS control points. Geocoding was done using GAMMA software. A sample data of DEM of Shuttle Radar Topography Mission (SRTM) of the Bhuj area is made available by DLR Germany. The intensity image, DEM and Error map are well registered. The spatial resolution of this DEM is about 25 m with height accuracy of a few meters. The DEM derived through ERS SAR data is prone to atmospheric affects as the required two images are acquired in different timings where as SRTM acquired the two images simultaneously. An RMS height error of 12.06 m is observed with reference to SRTM though some of the individual locations differ by as much as 35 m.     

Advanced analysis of differences between C and X bands using SRTM data for mountainous topography

By Interferometric Synthectic Aperture Radar (InSAR), during the Shuttle Radar Topography Mission (SRTM) height models have been generated, covering the earth surface from 56° south to 60.25° north. With the exception of small gaps in steep parts, dry sand deserts and water surfaces, the free available US C-band data cover the earth surface from 56° south to 60.25° north completely while the X-band data, distributed by the DLR (German Aerospace Center), cover it only partially. The C-band and X-band radar cannot penetrate the vegetation because of the short wavelength. Therefore, the height models are not Digital Elevation Models (DEM) representing bare Earth surface without any details, they are Digital Surface Models (DSM) representing the visible surface including vegetation and buildings. In the area of Zonguldak, Turkey, C-band and X-band DSMs are available and have been analysed in cooperation between Zonguldak Karaelmas University (ZKU) and Leibniz University of Hannover. The digitized contour lines from the 1:25,000 scale topographic maps and also a more precise height model derived directly from large scale photogrammetric mapping are used as reference height models. The terrain inclination influences the accuracy strongly, but also the directions of the inclination in relation to the radar view direction, the aspects, are important. Independent from the aspects, the analysed results do have root mean square differences against the reference data fitting very well to the Koppe formula SZ=a+b*tan α. The analyses are made separately for open and forest areas, with clear accuracy differences between both. Also, the analysis of X-band separately for three sub-areas is done and the positive effect of double observation to the accuracy has been clearly determined. The C-band data are only available with a spacing of 3 arcsec, corresponding to 92m × 70m, while the X-band data do have a spacing of 1 arcsec. This is important for the interpolation in the mountainous test area. The accuracy of the height points is approximately the same for the C- and the X-band data. But the C-band data which have three times larger spacing than Xband data, do not include the same morphological information. While C-band data contain very generalised contour lines X-band data have quite more details depending on 1 arcsec point spacing. The differential DEMs have been generated, separately, for displaying the differences between SRTM height models and reference DEMs of the test field.     

InSAR DEM精度与地形特征的关系分析

为研究InSARDEM与地形特征的关系,本文以从不同空间位置获取的两幅SAR影像作为实验数据,将InSARDEM与USGSDEM进行比较,分析了InSARDEM的精度,并研究其与坡度、坡向之间的关系。结果表明,本次实验InSARDEM与USGSDEM高程差异中误差为+19.11m,其精度与地形特征强烈相关,随着坡度的增加,InSARDEM精度降低,且前坡处高程精度高于后坡。     

基于SRTM DEM的InSAR高分辨率山区地表高程重建算法

山体的叠掩和阴影现象造成的信号去相关,一直是InSAR重建山区地表高程的瓶颈之一.为此,提出了一种新的基于粗分辨率SRTM DEM(约90m分辨率)辅助InSAR数据重建山区地表高程的方法.利用SRTM DEM模拟的干涉相位,对ERS-1/2干涉相位做去地形相位处理,得到残余相位.通过对解缠后的残余相位计算方差提取叠掩和阴影区域的噪声,并用平均相位近似恢复噪声区域的相位,然后将其转换为高程,并用SRTM DEM作高程补偿处理,从而实现地表高程重建.最后,定量比较了该方法与传统InSAR技术生成的DEM精度.实验表明,这种方法能有效提高传统InSAR技术生成地表高程的精度,这对提高星载雷达数据的使用效率具有重要意义.     

SPOT HRG影像无控制正射纠正实验与精度分析

本文根据影像正射纠正的基本原理和方法,提出了无控制点条件下的影像正射纠正理论估算公式,并据此公式对无控制条件点下的影像进行正射纠正理论精度估算,表明了在SRTMDEM支持下进行SPOT影像正射纠正具有较高的精度。并按此方法对杭州、重庆和大庆三个不同地形的区域制作DOM,然后利用我国1:1万地形图对其进行精度评定,实验结果证明正射纠正精度估算公式正确,同时利用SRTMDEM对SPOTHRG进行正射纠正得到的影像能够满足1:50000测图精度的要求。     

基于大比例尺地形图DEM快速建立方法的研究

利用现有大比例尺地形图快速生成DEM,通过对原数据格式的分析,利用程序采用基于高程点和等高线2种方式从dwg文件中批量提取出高程信息,高效、快速地获取大比例尺地形图中高精度的高程数据;使用VC^ 和ObjectARX编程技术进行DEM建模,在此过程中,通过采取优化数据结构等几种改进措施有效地减小了数据冗余并大大减少了程序运行的CPU时,快速制作出大区域高分辨率DEM。     

The influence correction of boreal forest vegetation on SRTM data

The digital elevation model based on SRTM is very convenient for a wide range of studies but requires correction due to the influence of forest vegetation. The present study was conducted to analyse the effect of boreal forests on altitudes, aspects and slopes calculated from the SRTM. A DEM based on topographic maps at 1:100 000 scale was used as a reference. The linear regression analysis showed low data correlation in forested areas. The presence of different types of forests and felling in the woods leads to a complex distribution of deviations from the SRTM. A simple correction method was proposed, using a forest mask, built according to Landsat, and forest heights indicated on the topographic maps. After correction, the correlation coefficient between the altitudes increased by 0.05–0.14, the share of matching aspects by 1–4% and the share of matching slopes by 2–8%.     

利用参数独立分解的星载SAR干涉测量检校方法

在星载干涉合成孔径雷达中,干涉参数的准确性对高程精度起着至关重要的作用。传统干涉测量检校方法往往将影响量级不同的干涉参数组合在一起解算,无法精确获得每项干涉参数的修正量。针对此问题,本文提出一种利用参数独立分解的干涉测量检校方法。首先,根据三维重建模型,确定与干涉SAR测高有关的参数;随后,在确保几何参数精度的前提下,对干涉参数的敏感度进行定量分析;最后,采用独立检校算法解算每一项干涉参数误差,完成干涉测量检校模型的建立。本文选择陕西渭南区域4对TerraSAR-X/TanDEM-X数据进行了试验分析。结果表明,对于该试验数据,采用本文提出的参数独立分解方法,干涉测量检校后干涉结果的高程精度优于2.54 m,平原地区获取DEM的绝对高程精度优于1.21 m,山地地区获取DEM的绝对高程精度优于3.11 m,验证了本文方法的有效性和正确性,为我国平原及山区1∶25 000比例尺的干涉SAR地形图测绘提供了技术基础。     

基于ICESat-2 ATLAS数据的DEM高程精度评价

为探究ASTER GDEMV3、SRTM1 DEM和AW3D30 DEM 3种开源DEM数据的高程精度,本文以高精度ICESat-2 ATLAS测高数据为参考数据,利用GIS统计分析、误差相关分析及数理统计对DEM的高程精度进行对比评价。结果表明：①AW3D30的质量最稳定;SRTM1 DEM在平原精度最高;在高原山地精度由高到低依次为AW3D30 DEM、ASTER GDEMV3、SRTM1 DEM。②DEM数据高程精度受地表覆盖影响较大,且与地形因素密切相关,在相同地表覆盖的两个研究区中DEM数据高程精度表现情况不一致,SRTM在平原地表覆盖下精度表现最好,平均误差为3.15 m,AW3D30 DEM在山地地表覆盖下精度表现最好,平均误差为7.61 m。③坡度对DEM数据的高程精度影响较大,在两个研究区3种DEM数据的高程误差均随坡度的增加而增加;坡向对DEM数据的高程精度影响较小,未发现明显的规律。     

北京市大比例尺地形图图幅号查询信息系统的研制

供应大比例尺地形图是城市测绘资料管理部门的基本职能。文中论述了在ARC／INFO基础上,以1：25000比例尺地形图上的点、线、面地物为索引,查询北京市1：500～1：10000各类比例尺地形图图幅号的算法。实验表明系统是实用的。