SRTM‐assisted block adjustment for stereo pushbroom imagery

KOMPSAT-2와 같은 고해상 위성영상은 대상영역의 3차원 위치결정을 위하여 RPC(Rational Polynomial Coefficient)가 포함된 자료를 제공한다. 그러나 RPC로 계산된 영상기하는 일정량의 편이(systematic errors)를 지니고 있는 상태이며, 이를 보정하기 위해서는 수 개 이상의 지상기준점(ground control point)이 필요하다. 이에 본 논문에서는 지상기준점 없이 입체영상(stereo pair)과 SRTM(Shuttle Radar Topography Mission) DEM(Digital Elevation Model) 사이의 대응점(tie point)만을 이용하여 자동으로 영상 기하를 보정하는 효과적인 방법을 제안하였다. 이러한 방법은 4가지 단계를 포함 한다: 1) 대응점 추출, 2) 대응점에 대한 지상좌표 결정, 3) SRTM DEM을 이용한 지상좌표의 보정, 4) RPC 보정 모델의 파라미터 결정. 우리는 KOMPSAT-2 입체영상을 이용하여 제안된 방법의 성과를 입증하였다. 검사점(check point)을 통해 계산된 RMSE(Root Mean Square Error)는 X와 Y, Z방향으로 각각 약 3.55 m, 9.70 m, 3.58 m를 나타냈다. 이는 SRTM DEM을 이용하여 RPC가 지닌 편이를 X, Y 및 Z 모든방향에 대하여 10 m이내의 정확도로 자동보정할 수 있다는 것을 의미한다. A high resolution satellite imagery such as KOMPSAT-2 includes a material containing rational polynomial coefficient (RPC) for three-dimensional geopositioning. However, image geometries which are calculated from the RPC must have inevitable systematic errors. Thus, it is necessary to correct systematic errors of the RPC using several ground control points (GCPs). In this paper, we propose an efficient method for automatic correction of image geometries using tie points of a stereo pair and the Shuttle Radar Topography Mission (SRTM) Digital Elevation Model (DEM) without GCPs. This method includes four steps: 1) tie points extraction, 2) determination of the ground coordinates of the tie points, 3) refinement of the ground coordinates using SRTM DEM, and 4) RPC adjustment model parameter estimation. We validates the performance of the proposed method using KOMPSAT-2 stereo pair. The root mean square errors (RMSE) achieved from check points (CPs) were about 3.55 m, 9.70 m and 3.58 m in X, Y;and Z directions. This means that we can automatically correct the systematic error of RPC using SRTM DEM.

Uncrewed aircraft systems (UASs) and structure-from-motion/multi-view stereo (SfM/MVS) photogrammetry are efficient methods for mapping terrain at local geographic scales. Traditionally, indirect georeferencing using ground control points (GCPs) is used to georeference the UAS image locations before further processing in SfM software. However, this is a tedious practice and unsuitable for surveying remote or inaccessible areas. Direct georeferencing is a plausible alternative that requires no GCPs. It relies on global navigation satellite system (GNSS) technology to georeference the UAS image locations. This research combined field experiments and simulation to investigate GNSS-based post-processed kinematic (PPK) as a means to eliminate or reduce reliance on GCPs for shoreline mapping and charting. The study also conducted a brief comparison of real-time network (RTN) and precise point positioning (PPP) performances for the same purpose. Ancillary experiments evaluated the effects of PPK base station distance and GNSS sample rate on the accuracy of derived 3D point clouds and digital elevation models (DEMs). Vertical root mean square errors (RMSEz), scaled to the 95% confidence interval using an assumption of normally-distributed errors, were desired to be within 0.5 m to satisfy National Oceanic and Atmospheric Administration (NOAA) requirements for nautical charting. Simulations used a Monte Carlo approach and empirical tests to examine the influence of GNSS performance on the quality of derived 3D point clouds. RTN and PPK results consistently yielded RMSEz values within 10 cm, thus satisfying NOAA requirements for nautical charting. PPP did not meet the accuracy requirements but showed promising results that prompt further investigation. PPK experiments using higher GNSS sample rates did not always provide the best accuracies. GNSS performance and model accuracies were enhanced when using base stations located within 30 km of the survey site. Results without using GCPs observed a direct relationship between point cloud accuracy and GNSS performance, with R2 values reaching up to 0.97.

The past decade has witnessed the rapid development of the SRTM (Shuttle Radar Topography Mission) DEM (digital elevation model) in engineering applications and scientific research. The near-global SRTM DEM was generated based on radar interference theory. The latest version of the SRTM DEM with a resolution of 1 arc-second has been widely used in various applications. However, many studies have shown the poor elevation accuracy of the SRTM DEM in forested areas. Recent developments in the field of spaceborne lidar have provided an additional chance to correct the elevation error of the SRTM DEM in forested areas. We developed an easy-to-use method to correct the elevation error of the SRTM DEM based on the spatial interpolation method using the recent Ice, Cloud and land Elevation Satellite-2 data. First, an ICESat-2 terrain control point selection criterion was proposed to reject some erroneous ICESat-2 terrains caused by many factors. Second, we derived the elevation correction surface based on the interpolation method using the refined ICESat-2 terrain. Finally, a corrected SRTM DEM of forested areas was generated through the obtained elevation correction surface. The proposed method was tested in the typical forested area located in Massachusetts, USA. The results show that the RMSE of the selected terrain control points in vegetation areas and non-vegetation areas are 1.03 and 0.68 m, respectively. The corrected SRTM DEM have an RMSE of 4.2 m which is significantly less than that of the original SRTM DEM with an RMSE of 9.8 m, which demonstrates the proposed method is feasible to correct the elevation error in forested areas. It can be concluded that the proposed method obviously decreases the elevation error of the original SRTM DEM.

Tianwen-1 is the first Mars probe launched by China and the first mission in the world to successfully complete the three steps of exploration (orbiting, landing, and roving) at the one time. Based on the unverifiable descent images which cover the full range of the landing area, trajectory recovery and fine terrain reconstruction are important parts of the planetary exploration process. In this paper, a novel trajectory recovery and terrain reconstruction (TR-TR) algorithm employing descent images is proposed for the dual-restrained conditions: restraints of the flat terrain resulting in an unstable solution of the descent trajectory and of the parabolic descent trajectory causing low accuracy of terrain reconstruction, respectively. A landing simulation experiment on a landing field with Mars-like landform was carried out to test the robustness and feasibility of the algorithm. The experiment result showed that the horizontal error of the recovered trajectory didn’t exceed 0.397 m, and the elevation error of the reconstructed terrain was no more than 0.462 m. The algorithm successfully recovered the descent trajectory and generated high-resolution terrain products using in-orbit data of Tianwen-1, which provided effective support for the mission planning of the Zhurong rover. The analysis of the results indicated that the descent trajectory has parabolic properties. In addition, the reconstructed terrain contains abundant information and the vertical root mean square error (RMSE) of ground control points is smaller than 1.612 m. Terrain accuracy obtained by in-orbit data is lower than that obtained by field experiment. The work in this paper has made important contributions to the surveying and mapping of Tianwen-1 and has great application value.

Satellite-derived Digital Elevation Model (DEM) selection, preparation and correction for hydrodynamic modelling in large, low-gradient and data-sparse catchments

The Shuttle Radar Topography Mission (SRTM) provided scientists with digital elevation data on a nearly global scale and with highly consistent accuracy. This paper compares elevation values of the C-band SRTM 30-m digital elevation model (DEM) with pointwise elevations from the Ice, Cloud, and land-Elevation Satellite (ICESat) laser altimetry for Otter Tail County, Minnesota. The accuracy of SRTM DEM is measured as a function of land covers and geomorphologic characteristics. The typical mean vertical difference between the SRTM DEM and ICESat elevations in this paper was determined in each classified land-use type and is approximately 1.5 m over bare ground, with the SRTM measuring lower elevations. Significant changes in the SRTM DEM uncertainties have been identified over different surface types classified from Landsat-7 imagery, e.g., bare ground, urban, and forested areas. Based on this result, the difference of the SRTM 30-m DEM and ICESat elevations has been removed from the DEM and made available for improved hydrological applications

Along with the appearance of high resolution satellite images, image correction using Rational Polynomial Coefficients (RPCs) has become common. Location accuracy of Korea Multipurpose Satellite (KOMPSAT) standard images is still not adequate, so, in order to leverage the KOMPSAT images for applications such as mapping and change detection, it is necessary to orthorectify the images. In this study, using updated RPCs, we performed orthorectification of KOMPSAT-2, KOMPSAT-3, and KOMPSAT-3A images using various data. Through this study, we discovered that the orthorectification result using precise Ground Control Points (GCPs) and Digital Elevation Model (DEM) is the best, but it was found that the correction results through image matching are also excellent. In particular, it was confirmed that orthoimages with a planimetric accuracy around 3 m (Root Mean Square Error (RMSE)) can be generated by using well-known matching algorithms with open data such as OpenStreetMap (OSM) and Shuttle Radar Topography Mission (SRTM) DEM, which can be acquired by anyone. Although the accuracy was low in some mountainous terrain, it was confirmed that it could be used for generating KOMPSAT orthoimages using open data. This paper describes the results for orthorectifying high resolution KOMPSAT optical images using various reference data.

Assimilation of Clear-sky FY-4A AGRI radiances within the WRFDA system for the prediction of a landfalling Typhoon Hagupit (2020)

Accurate measurements of glacier surface topography and their changes play an essential role in various glaciological studies related to glacier dynamics and mass balance. The focus of this study is on mapping glacier digital elevation model (DEM) and elevation changes in the western Qilian Mountains, northern Tibetan Plateau, by synergistically using the TanDEM-X (TDX) bistatic Interferometric Synthetic Aperture Radar (InSAR) data in 2013 and Shuttle Radar Topography Mission (SRTM) DEM in 2000. The first high-resolution and high-precision glacier DEM is derived in this region by a TDX InSAR procedure with a nonlocal (NL) filter. Validated against the Ice, Cloud, and land Elevation Satellite (ICESat) height references, the absolute height error of the TanDEM-X DEM derived with the NL filter and the Goldstein filter with the parameters investigated is, respectively, 1.493 ± 0.747 and 1.857 ± 1.709 m. Further, four combinations of differential phase method (DiffPha) and DEM differencing method (DiffDem) with Goldstein filter and NL filter are applied to estimate glacier elevation changes between 2000 and 2013. The synergistic use of the DiffPha method and the NL filter is superior to other three combinations in terms of uncertainty and noise reduction. Generally, a clear surface thinning can be found in most glacier tongue regions, the maximum value of elevation lowering up to approximately −40 m, whereas a slight thickening is detected in accumulation areas, which are in agreement with the height difference results between GPS measurements and SRTM DEM over Laohugou Glacier No.12. This study demonstrates the potential of the TanDEM-X bistatic InSAR in mapping surface topography and elevation changes of valley glaciers in the Tibetan Plateau.

Aiming to the difficult area of surveying and mapping, such as Antarctic, it is difficult to guarantee favourable observation condition for traditional mapping methods. Large-scale block adjustment (BA) without ground control point (GCP) is carried out by using domestic high-resolution satellite remote sensing imagery. The key technologies of BA model construction based on virtual control point (VCP), gross error detection and elimination, and robust and fast method of large-scale BA are studied in this paper. On this basis, the experimental analysis and validation of 13pairs (39 images) of ZY-3 satellite imagery in Antarctic are carried out. The results show that the sparse matrix technology can effectively reduce the memory requirement. The combined matrix block and GPU parallel technology can solve the problem of large-scale BA computational efficiency. In addition, after BA, the maximum residual is 3.920 pixel, the root mean square error (RMSE) is 0.169 pixel in the X (flight) direction, the maximum residual error is 5.933 pixel, and the RMSE is 0.191 pixel in the Y (scan) direction. The proposed method has certain accuracy and stability in large-scale BA without GCP using high-resolution satellite remote sensing imagery in Antarctica. The relative positioning accuracy can reach sub-pixel level, which can meet the requirements of cartographic splicing.

Digital elevation models (DEMs) are the most obvious data sources in landslide susceptibility assessment. Many landslide casual factors are often generated from DEMs. Most studies on landslide susceptibility assessments rely on freely available DEMs. However, very little is known about the performance of different DEMs with varying spatial resolutions on the accurate assessment of landslide susceptibility. This study compared the performance of four different DEMs including 30 m Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Global Digital Elevation Model (GDEM), 30–90 m Shuttle Radar Topographic Mission (SRTM), 12.5 m Advanced Land Observation Satellite (ALOS) Phased Array Type L band Synthetic Aperture Radar (PALSAR), and 25 m Survey of Bangladesh (SOB) DEM in landslide susceptibility assessment in the Rangamati district in Bangladesh. This study used three different landslide susceptibility assessment techniques: modified frequency ratio (bivariate model), logistic regression (multivariate model), and random forest (machine-learning model). This study explored two scenarios of landslide susceptibility assessment: using only DEM-derived causal factors and using both DEM-derived factors as well as other common factors. The success and prediction rate curves indicate that the SRTM DEM provides the highest accuracies for the bivariate model in both scenarios. Results also reveal that the ALOS PALSAR DEM shows the best performance in landslide susceptibility mapping using the logistics regression and the random forest models. A relatively finer resolution DEM, the SOB DEM, shows the lowest accuracies compared to other DEMs for all models and scenarios. It can also be noted that the performance of all DEMs except the SOB DEM is close (72%–84%) considering the success and prediction accuracies. Therefore, anyone of the three global DEMs: ASTER, SRTM, and ALOS PALSAR can be used for landslide susceptibility mapping in the study area.

The correct representation of the topography of terrain is an important requirement to generate photogrammetric products such as orthoimages and maps from high-resolution (HR) or very high-resolution (VHR) satellite datasets. The refining of the digital elevation model (DEM) for the generation of an orthoimage is a vital step with a direct effect on the final accuracy achieved in the orthoimages. The refined DEM has potential applications in various domains of earth sciences such as geomorphological analysis, flood inundation mapping, hydrological analysis, large-scale mapping in an urban environment, etc., impacting the resulting output accuracy. Manual editing is done in the presented study for the automatically generated DEM from IKONOS data consequent to the satellite triangulation with a root mean square error (RMSE) of 0.46, using the rational function model (RFM) and an optimal number of ground control points (GCPs). The RFM includes the rational polynomial coefficients (RPCs) to build the relation between image space and ground space. The automatically generated DEM initially represents the digital surface model (DSM), which is used to generate a digital terrain model (DTM) in this study for improving orthoimages for an area of approximately 100 km2. DSM frequently has errors due to mass points in hanging (floating) or digging, which need correction while generating DTM. The DTM assists in the removal of the geometric effects (errors) of ground relief present in the DEM (i.e., DSM here) while generating the orthoimages and thus improves the quality of orthoimages, especially in areas such as Dehradun that have highly undulating terrain with a large number of natural drainages. The difference image of reference, i.e., edited IKONOS DEM (now representing DTM) and automatically generated IKONOS DEM, i.e., DSM, has a mean difference of 1.421 m. The difference DEM (dDEM) for the reference IKONOS DEM and generated Cartosat-1 DEM at a 10 m posting interval (referred to as Carto10 DEM) results in a mean difference of 8.74 m.

The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) aboard the Terra satellite was designed to generate along‐track stereo images. The data are available at low cost, providing a feasible opportunity for generating digital elevation models (DEMs) in areas where little or no elevation data are yet available. This study evaluates the accuracy of DEMs extracted from ASTER data covering mountainous terrain. For an assessment of the achieved accuracies in the Andean study site, comparisons were made to similar topographical conditions in Switzerland, where reference data were available. All raw DEMs were filtered and interpolated by the post‐processing tools included with PCI Geomatica, the software package used. After carefully checking the DEM quality, further post‐processing was undertaken to eliminate obvious artefacts such as peaks and sinks. Accuracy was tested by comparing the DEMs in the Swiss Alps to three reference models. The achieved results of the generated DEMs are promising, considering the extreme terrain. Given accurate and well‐distributed ground control points (GCPs), it is possible to generate DEMs with a root mean square (RMS) error between 15 m and 20 m in hilly terrain and about 30 m in mountainous terrain. The DEMs are very accurate in nearly flat regions and on smooth slopes with southern expositions: errors are generally within ±10 m in those cases. Larger errors do appear in forested, snow covered or shady areas and at steep cliffs and deep valleys with extreme errors of a few hundred metres. The evaluation showed that the quality of the DEMs is sufficient for enabling atmospheric, topographic and geometric correction to various satellite datasets and for deriving additional products.

The launch of the very high resolution (VHR) sensor satellites has paved the way for further exploitation of the capabilities of stereo satellite imagery for many applications and most importantly in the extraction of topographical details. Advances in digital technologies in the photogrammetric technique provide more accurate and up to date topographic representation, while traditional methods of topographic mapping are time-consuming and require large technical support. This paper explores the methods of production and an analysis of the capabilities of the high-resolution stereo imaging system for the extraction of features (topographic map) and generation of digital elevation model (DEM). This study uses high-resolution GeoEye-1 stereo pair images. For qualitative assessment, topographic features were extracted by digitizing the georeferenced stereo pair images and overlapped to form a DEM and orthophoto within the Summit Evolution software. Georeferencing was carried out using the existing Ground Control Points (GCPs) and Check Points (CPs) established using Real-Time Kinematic Global Navigation Satellite System (RTK-GNSS) technique. For quantitative assessment, a total of 7 existing GCPs were considered suitable for use and 5 CPs were established as a check for accuracy of the output. The accuracy of the CPs was evaluated by computing Root Mean Square Error (RMSE). The average horizontal RMSE for the left image and right image were 1.17150098m and 0.819334672m respectively and the average vertical RMSE was 1.605898565m. Additionally, the standard deviation for elevation data from generated DEM using photogrammetric technique and existing elevation data of RTK-GNSS technique were computed, which resulted in 15.91286667 and 16.59529228 respectively. The accuracy was evaluated by comparing the results of the data extracted using stereo satellite images and those extracted from ground measurement techniques (Total station and RTK-GNSS). The comparative analysis showed that stereo pair images provide a viable alternative to other ground mapping techniques providing results within the acceptable range. Based on error estimation and analysis, it concludes that, if strict photogrammetric processing model and ground control points are employed, high-resolution satellite imagery can be used for the accurate generation and update of the topographic maps, DEM and digital orthophotos, which can be further be used as a reference for comparative analysis applications.

Recently, there has been increasing attention to the digital elevation model (DEM) because of its ability to learn the Earth’s surface’s topography. Freely accessible DEMs such as the Cartosat-1, shuttle radar topography mission (SRTM), light detection and ranging, and advanced spaceborne thermal emission and reflection radiometer DEM, contain large vertical errors. The errors are aggravated over multifaceted geography and cannot rectify microgeographic deviations in the moderately flat landscape. As high-accuracy DEMs have limited availability, dated low-accuracy DEMs are still used in various models, specifically in the data-sparse areas. However, it is necessary to enhance the quality of these DEMs before their use in the geomorphometric analysis. We aim to investigate the effect of noise reduction filters on DEM’s accuracy and quality. The noise reduction filters such as weighted average filter, median filter, sharpen filter, Lee sigma filter, and local sigma filter are used for noise reduction in the stereo images. The quality and accuracy of the generated DEMs are further improved by selecting an optimum number of tie points in the image matching process. The effects are assessed by correlating the surface profiles for the final DEM obtained and SRTM DEM as a reference with a resolution of 1 arc sec.