Creation of DTM with ASTER Data and Statistical Verification of the Accuracy of the Model (Western Peloponnese,Greece)

Abstract

In this paper we present the use of ASTER data for the creation of a Digital Terrain Model (DMT) of high accuracy. Using a stereo pair of ASTER satellite images with 15m resolution we created two DMTs: one with a 30m pixel size and another one with a 15m pixel size. Then we made a statistical verification of the two DTMs accuracy. We created another DTM with 30m pixel size from digitized contours of 1:50000 scale topographic maps. We first made an optical comparison of the two DTMs with 30m pixel size. Then we subtracted the two DTMs and we presented their difference. Finally, we verified the DTMs accuracy using 68 points of a well‐known elevation. All the results demonstrated that DTMs derived from ASTER data have very good accuracy.     

SPIRIT. SPOT 5 stereoscopic survey of Polar Ice: Reference Images and Topographies during the fourth International Polar Year (2007–2009)

Monitoring the evolution of polar glaciers, ice caps and ice streams is of utmost importance because they constitute a good indicator of global climate change and contribute significantly to ongoing sea level rise. Accurate topographic surveys are particularly relevant as they reflect the geometric evolution of ice masses. Unfortunately, the precision and/or spatial coverage of current satellite missions (radar altimetry, ICESat) or field surveys are generally insufficient. Improving our knowledge of the topography of Polar Regions is the goal of the SPIRIT (SPOT 5 stereoscopic survey of Polar Ice: Reference Images and Topographies) international polar year (IPY) project. SPIRIT will allow (1) the acquisition of a large archive of SPOT 5 stereoscopic images covering most polar ice masses and, (2) the delivery of digital terrain models (DTM) to the scientific community.Here, we present the architecture of this project and the coverage achieved over northern and southern polar areas during the first year of IPY (July 2007 to April 2008). We also provide the first accuracy assessments of the SPIRIT DTMs. Over Jakobshavn Isbrae (West Greenland), SPIRIT elevations are within ±6 m of ICESat elevations for 90% of the data. Some comparisons with ICESat profiles over Devon ice cap (Canada), St Elias Mountains (Alaska) and west Svalbard confirm the good overall quality of the SPIRIT DTMs although large errors are observed in the flat accumulation area of Devon ice cap. We then demonstrate the potential of SPIRIT DTMs for mapping glacier elevation changes. The comparison of summer-2007 SPIRIT DTMs with October-2003 ICESat profiles shows that the thinning of Jakobshavn Isbrae (by 30–40 m in 4 years) is restricted to the fast glacier trunk. The thinning of the coastal part of the ice stream (by over 100 m) and the retreat of its calving front (by up to 10 km) are clearly depicted by comparing the SPIRIT DTM to an ASTER April-2003 DTM.     

Averaging and Formulation Impact on GB-SAR Topographic Mapping

A novel processor for the generation of a digital terrain model (DTM) using ground-based synthetic aperture radar (GB-SAR) imagery is presented. The performance of the processor has been assessed using the archive of C-band GB-SAR images available at the Joint Research Centre from the Sion Valley (CH) test site, from which a high-resolution ground-truth DTM is available. In order to reduce the standard deviation of the differences between the two DTMs, several techniques have been applied: interferogram averaging, phase referencing, and differential path calibration. In addition, different formulations have been reviewed, and its impact has been quantified in the DTM accuracy obtained. Results show that the standard deviation of the differences between the ground-truth DTM and that of the C-band GB-SAR can be reduced from the 5 m that was previously reported down to 3 m with the new processor.      

Improving Cartosat-1 DEM accuracy using synthetic stereo pair and triplet

Cartosat–1 is the first Indian Remote Sensing Satellite capable of providing along-track stereo images. Cartosat–1 provides forward stereo images with look angles +26° and −5° with respect to nadir for generating Digital Elevation Models (DEMs), Orthoimages and value added products for various applications. A pitch bias of −21° to the satellite resulted in giving reverse tilt mode stereo pair with look angles of +5° and −26° with respect to nadir. This paper compares DEMs generated using forward, reverse and other possible synthetic stereo pairs for two different types of topographies. Stereo triplet was used to generate DEM for Himalayan mountain topography to overcome the problem of occlusions.For flat to undulating topography it was shown that using Cartosat-1 synthetic stereo pair with look angles of −26° and +26° will produce improved version of DEM. Planimetric and height accuracy (Root Mean Square Error (RMSE)) of less than 2.5 m and 2.95 m respectively were obtained and qualitative analysis shows finer details in comparison with other DEMs. For rugged terrain and steep slopes of Himalayan mountain topography simple stereo pairs may not provide reliable accuracies in DEMs due to occlusions and shadows. Stereo triplet from Cartosat-1 was used to generate DEM for mountainous topography. This DEM shows better reconstruction of elevation model even at occluded region when compared with simple stereo pair based DEM. Planimetric and height accuracy (RMSE) of nearly 3 m were obtained and qualitative analysis shows reduction of outliers at occluded region.     

About the International GIS Workshop on Antarctic King George Island

Digital Elevation Models (DEM) of a hilly–valley region are prepared using stereo images of Cartosat-1 and Shuttle Radar Topography Mission (SRTM) images. The procedure of ortho-image generation from Cartosat-1 stereo images and the estimation of ground features from ortho-image are elaborated in the paper. Comparison of DEMs prepared from both images is discussed in terms of the quality of ground features detection, hydrological applications and geometrical calculations. It is found that DEM prepared from Cartosat-1 images are more accurate in the valley region and hence it is better suited for hydrological applications. On the contrary, for hilly region, SRTM images produce better DEM. However, if ground control points and Rational Polynomial Coefficients can be obtained in the hilly region, more accurate DEM can be prepared using Cartosat-1 stereo images.     

Geological study of sedimentary basins using SPOT data

This paper illustrates the benefits of DTMs created from SPOT images for the exploration of sedimentary basins. We chose an example located in the Ebro sedimentary basin in Spain, characterized by good outcropping conditions and slight deformations. The data used consist of a pair of SPOT panchromatic images and a SPOT XS image. The work consists of making up a 3D database, followed by interpretation of stereo pairs computed from orthoimages and the DTM. This interpretation is made on a stereoscopic desk and results in a digital file containing in the form of vectors all the observed faults and lithology. These vectors are then used to make calculations on the geometry of the objects they represent: we demonstrate that we can accurately measure layer directions and dips, sediment thicknesses and fault throws. Synthesis of perspective views made at the same time provide a good understanding of the structures and help to test their geometric consistency. Thus knowledge about relief given by DTMs helps to interpret remote detection images in 3D space, and particularly to accurately quantify the results of this interpretation.     

Cartosat-1 height product and ICESat/GLAS data for digital elevation surface generation

Space born systems like Geoscience Laser Altimeter System (GLAS) onboard collect data for ice, cloud and Land. Elevation satellite (ICESat) collects an unparalleled data set as waveform over terrestrial targets, helps in evaluating the global elevation data. In this study we compared the Digital Elevation Surface (DES) generated by Cartosat-1 point data and DES generated by merging the Cartosat-1 data with ICESat data. Outputs in the form of interpolated surfaces were evaluated with the help of differential global positioning system (DGPS) points collected from study area. The study showed the results that the DES generated from Cartosat — 1 data had less elevation accuracy when compared with the DGPS data. While merging Cartosat-1 point height data with ICESat/GLAS data resulted in better accuracy. On the practical side for processing the interpolation, based on the research the ICESat /GLAS with Cartosat-1 height data can produce better DES compared to the Cartosat-1 stereo data. The DES was generated using geostatistical interpolation methods in which the global polynomial method proved to be the better for generating the surface compare to other interpolation techniques studied in this work. For co-kriging method, the accuracy decreases compare to the kriging interpolation, due to the complexity of parameters that were used for interpolation. On the theory side, based on this research the statement of which interpolation technique is better than the other cannot be mentioned easily, because these are based on the data type, parameters and also on method of interpolation. So research experiment should be more intensely and with more focused.     

应用遥感影像制作矿区大比例尺地形图

本文在专业的遥感和地理信息系统软件平台支持下,利用Aster立体像对和IKONOS高分辨率遥感影像分别提取矿区地形信息和地物信息,结合矿区的实际特点,参照地形图制作规范对矿区大比例尺地形图的制作进行了分析和探讨,结果表明利用遥感影像进行矿区大比例尺地形图制作不但周期短、速度快、时效性强,而且精度完全能够满足生产需要,是一种实用可行的方法。     

Timing of erosion and satellite data: A multi-resolution approach to soil erosion risk mapping

Erosion reduces soil productivity and causes negative downstream impacts. Erosion processes occur on areas with erodible soils and sloping terrain when high-intensity rainfall coincides with limited vegetation cover. Timing of erosion events has implications on the selection of satellite imagery, used to describe spatial patterns of protective vegetation cover. This study proposes a method for erosion risk mapping with multi-temporal and multi-resolution satellite data. The specific objectives of the study are: (1) to determine when during the year erosion risk is highest using coarse-resolution data, and (2) to assess the optimal timing of available medium-resolution images to spatially represent vegetation cover during the high erosion risk period. Analyses were performed for a 100-km² pasture area in the Brazilian Cerrados. The first objective was studied by qualitatively comparing three-hourly TRMM rainfall estimates with MODIS NDVI time series for one full year (August 2002–August 2003). November and December were identified as the months with highest erosion risk. The second objective was examined with a time series of six available ASTER images acquired in the same year. Persistent cloud cover limited image acquisition during high erosion risk periods. For each ASTER image the NDVI was calculated and classified into five equally sized classes. Low NDVI was related to high erosion risk and vice versa. A DEM was used to set approximately flat zones to very low erosion risk. The six resulting risk maps were compared with erosion features, visually interpreted from a fine-resolution QuickBird image. Results from the October ASTER image gave highest accuracy (84%), showing that erosion risk mapping in the Brazilian Cerrados can best be performed with images acquired shortly before the first erosion events. The presented approach that uses coarse-resolution temporal data for determining erosion periods and medium-resolution data for effective erosion risk mapping is fast and straightforward. It shows good potential for successful application in other areas with high spatial and temporal variability of vegetation cover.     

Characterization of ASTER GDEM elevation data over vegetated area compared with lidar data

Current researches based on areal or spaceborne stereo images with very high resolutions (<1 m) have demonstrated that it is possible to derive vegetation height from stereo images. The second version of the Advanced Spaceborne Thermal Emission and Reflection Radiometer Global Digital Elevation Model (ASTER GDEM) is the state-of-the-art global elevation data-set developed by stereo images. However, the resolution of ASTER stereo images (15 m) is much coarser than areal stereo images, and the ASTER GDEM is compiled products from stereo images acquired over 10 years. The forest disturbances as well as forest growth are inevitable in 10 years time span. In this study, the features of ASTER GDEM over vegetated areas under both flat and mountainous conditions were investigated by comparisons with lidar data. The factors possibly affecting the extraction of vegetation canopy height considered include (1) co-registration of DEMs; (2) spatial resolution of digital elevation models (DEMs); (3) spatial vegetation structure; and (4) terrain slope. The results show that the accurate coregistration between ASTER GDEM and national elevation dataset (NED) is necessary over mountainous areas. The correlation between ASTER GDEM minus NED and vegetation canopy height is improved from 0.328 to 0.43 by degrading resolutions from 1 arc-second to 5 arc-second and further improved to 0.6 if only homogenous vegetated areas were considered.     

DEM约束的卫星影像定位法

作为"云控制"摄影测量理论和方法的发展,研究了DEM约束的立体卫星影像区域网平差方法。与DEM仅作为高程控制信息使用,或者是通过DEM表面匹配实现绝对定向的间接定位方法不同,DEM作为平高控制信息被直接引入至基于RFM模型的卫星影像区域网平差之中。本文方法将连接点地面高程与DEM格网内插高程之差作为虚拟观测值构建约束方程,不仅利用了DEM高程信息,并且利用了其地形曲面包含的平面信息,以"云控制"方式在区域网平差过程中有效消除卫星影像RPC参数中包含的整体偏移及区域网内部的扭曲变形,实现了无地面控制点条件下卫星影像平面及高程绝对定位精度的大幅提升。使用覆盖山东全境的330景天绘一号立体卫星影像进行试验,分别以AW3D30、ASTER GDEM和SRTM GL3共3种开源DEM作为控制信息,并使用100个外业实测控制点进行精度评测。试验表明,以DEM作为控制可显著提高区域网平差的平面与高程精度,卫星影像绝对定位精度与DEM自身精度有关。当使用AW3D30作为控制时,可以取得与使用100个外业控制点平差同等精度,平面中误差为5.0 m(约1像素),高程中误差为2.9 m。试验结果证明了DEM替代外业控制点作为平差控制信息的有效性与可行性。     

不同软件在ASTER数据中提取DEM的精度对比

利用3个不同的软件对四川省龙门山中段ASTER 15 m分辨率的立体像对进行了DEM提取,并对其精度进行了初步评价。分别使用立体测量法和干涉测量法提取DEM,并通过检验点法和剖面线法对比分析。结果表明,利用ERDAS的干涉测量法提取出的DEM效果较好,高程精度可达30 m,对后续数据深挖掘和高层次地形分析具有应用价值。     

Applying GNSS and DTM Technologies to Monitor the Ice Balance of the Horcones Inferior Glacier, Aconcagua Region, Argentina

This work provides results and analysis of the behavior in a reconstituted debris covered glacier, Horcones Inferior (HIG). This glacier is located at 32° 41′ S and 69° 57′ W in the Aconcagua region, Mendoza, Argentina. The HIG has experienced surges events in 1984 and 2003. GNSS techniques and Digital Terrain Models (DTMs) derived from ASTER optical imagery, were used to detect movements and altimetric changes of the HIG in the 2001–2012 period. Based on a GNSS semi - continuous station, N and E mean velocities of 0.9 cm/d and 1.1 cm/d and in U direction were 1.1 cm/d and 1.5 cm/d were obtained respectively in 2009–2010 period. Additionally, GNSS profiles showed a velocity range between 0.5 cm/d and 2.8 cm/d. Finally, through a DTM comparison the altimetric change before and after the last surge in 2003 was obtained. The applied techniques allowed accurate and reliable change detection of the HIG.     

Toward accountable land use mapping: Using geocomputation to improve classification accuracy and reveal uncertainty

The classification of satellite imagery into land use/cover maps is a major challenge in the field of remote sensing. This research aimed at improving the classification accuracy while also revealing uncertain areas by employing a geocomputational approach. We computed numerous land use maps by considering both image texture and band ratio information in the classification procedure. For each land use class, those classifications with the highest class-accuracy were selected and combined into class-probability maps. By selecting the land use class with highest probability for each pixel, we created a hard classification. We stored the corresponding class probabilities in a separate map, indicating the spatial uncertainty in the hard classification. By combining the uncertainty map and the hard classification we created a probability-based land use map, containing spatial estimates of the uncertainty. The technique was tested for both ASTER and Landsat 5 satellite imagery of Gorizia, Italy, and resulted in a 34% and 31% increase, respectively, in the kappa coefficient of classification accuracy. We believe that geocomputational classification methods can be used generally to improve land use and land cover classification from imagery, and to help incorporate classification uncertainty into the resultant map themes.     

TERRA: Terrain Extraction from elevation Rasters through Repetitive Anisotropic filtering

Over the past decades, several filters have been developed to derive a Digital Terrain Model (DTM) from a Digital Surface Model (DSM), by means of filtering out aboveground objects such as vegetation. In this filtering process, however, one of the major challenges remains to precisely distinguish sharp terrain features, e.g. ridges, agricultural terraces or other anthropogenic geomorphology such as open-pit mines, riverbanks or road ramps. Hence, loss of elevation data around terrain edges (and consequent smoothing) is very common with existing algorithms. In terraced landscapes, the preservation of precise geomorphology is of key importance in digital terrain analyses, such as hydrologic and erosion modelling, or automatic feature recognition and inventorying. In this work, we propose a new filtering method called TERRA (Terrain Extraction from elevation Rasters through Repetitive Anisotropic filtering). The novelty of the algorithm lies within its usage of terrain aspect to guide the anisotropic filtering direction, therefore maximising the preservation of terrain edges. We derived six DTMs from DSMs using UAV Structure from Motion (SfM) photogrammetry, laser altimetry and satellite sources (grid resolutions ranging from 0.1–1.0 m). The results indicated a close agreement of DTMs filtered using the TERRA algorithm and reference DTMs, while terrace risers were well preserved even under thick canopies of vines and trees. Compared to existing filtering approaches, TERRA performed well in minimising Type I errors (false ground removal), while Type II errors occurred locally where vegetation was covering the terrace edges. Given the promising filtering performance, and supported by the minimal requirements of parameterisation and computation, the TERRA algorithm could be a useful tool in DTM preparation for digital terrain analysis of agricultural terraces and similar hillslopes characterised by a complex mosaic of sharp terrain and non-terrain features.     

Lithological mapping and fuzzy set theory: Automated extraction of lithological boundary from ASTER imagery by template matching and spatial accuracy assessment

Lithological boundaries provide information useful for activities such as mineral and hydrocarbon exploration, water resource surveys, and natural hazard evaluation. Automated detection of lithological boundaries reduces bias inherent in expert interpretation of boundaries and thus improves the reliability of lithological mapping. The Rotation Variant Template Matching (RTM) algorithm was applied to ASTER imagery to detect pre-defined lithological boundaries. Templates incorporating the mineral combinations gypsum–calcite and calcite–illite were designed to detect boundaries between evaporites, marly limestone, and sandstone. The RTM algorithm successfully detected lithological boundaries by rotating the templates over the ASTER imagery. The accuracy of the detected boundaries was spatially assessed using fuzzy set theory. Boundaries from a published geological map and boundaries interpreted from a stereo pair of aerial photos by five experts were used as references for assessing the accuracy. A confidence region unifying spatial errors was defined for the geological map and stereo-pair interpretation to provide boundary zones from these references. The correspondence between detected boundaries and the boundary zones of the aerial photo was better than between detected boundaries and boundary zones of the geological map.     

Comparing the terrain reversal effect in satellite images and in shaded relief maps: an examination of the effects of color and texture on 3D shape perception from shading

Terrain reversal effect (TRE) causes reversed 3D shape perception in satellite images and shaded relief maps (SRMs), and introduces difficulties in identifying landforms such as valleys and ridges. With this paper, in a controlled laboratory experiment, we compare how well 27 participants could identify valleys and ridges over 33 locations using SRMs, color satellite images and grayscale satellite images. The main depth cue is shadow both in vertical-view images and SRMs. However, the presence of texture and color in images also affect 3D shape perception. All our participants experience the illusion strongly: with the SRMs, it is very severe (2% accuracy), with grayscale images low but considerably better than SRMs (17.6% accuracy), and slightly worse with color imagery (15.3% accuracy). These differences between SRMs and imagery suggest that the participants who are able to bypass the illusion consciously or subconsciously interpret the photographic information. We support this observation further with a cue-strength analysis. Furthermore, we provide exploratory analyses of the effects of expertise, global convexity bias, and bistable perception. Our original empirical observations serve towards a better understanding of this visual illusion, and contribute towards nuanced and appropriate solutions to correcting for TRE differently for satellite images and SRMs.     

CBERS-02B卫星遥感影像的区域网平差

针对中巴资源一号卫星(CBERS-02B)卫星遥感影像姿态角误差较大的特点,提出了利用区域网平差方法提高其对地目标定位精度的策略和具体计算过程。首先对参与平差的每景影像选取4个地面控制点进行影像姿态角常差检校,然后采用与地形无关方案解求各自的RPC参数,最后选取带仿射变换项的有理函数模型(RFM)进行多重覆盖影像的区域网平差。对两个地区的0级CBERS-02B单条CCD立体影像对的区域网平差试验表明,对地目标定位的平面和高程精度均达到了±3个像元的水平,且高程精度明显优于平面精度。相对于常规的卫星遥感影像区域网平差方法,平面和高程精度均有明显提升,几乎达到国外同等高分辨率卫星遥感影像的几何定位精度。这说明中国卫星遥感影像亦具有较好的几何定位潜力,在区域网平差之前进行系统误差预改正是必要和可行的。     

高分辨率卫星影像几何精度真实性检验方法

以日本ALOS卫星为例,分析了线阵推扫式卫星影像成像特征,从理论上阐述利用严密成像几何模型建立卫星影像像点与物点之间的成像几何关系,并利用光束法区域网平差技术对影像几何精度进行真实性检验的方法。在实验部分,基于中国嵩山遥感几何定标场的参考数据,对ALOS卫星影像的几何精度进行了验证,获取了它的有控和无控的定位精度,并用解算的外方位元素建立立体模型,检测不同地形生成DEM精度。     

资源三号测绘卫星影像平面和立体区域网平差比较

针对弱交会条件下卫星遥感影像区域网平差无法正确求解的问题,本文提出了利用数字高程模型（DEM）作为高程约束的平面区域网平差方法提高其对地目标定位精度的策略。首先,选取带仿射变换项的有理函数模型（RFM）作为卫星影像平面区域网平差的数学模型。其次,在平差过程中更新连接点的地面坐标时仅求解地面点的平面坐标,高程值利用DEM进行内插获得。最后,在布设少量控制点的情况下通过平面区域网平差求解所有参与平差的卫星影像定向参数和连接点的地面平面坐标。利用两个地区的资源三号正视影像的平面区域网平差以及前正后三视影像的立体区域网平差的对比试验表明,对于资源三号卫星影像在1:50000DEM的支持下,平面平差可以达到和立体平差相当平面精度。对于近似垂直正视的资源三号影像,全球1km格网的DEM和90m格网的SRTM可以取代1:50000DEM作为高程控制,平面精度几乎没有损失。最终,试验结果证明了平面区域网平差方法的有效性和可行性。