Optimization of Optical Image Geometric Modeling, Application to Topography Extraction and Topographic Change Measurements Using PlanetScope and SkySat Imagery

The volume of data generated by earth observation satellites has increased tremendously over the last few decades and will increase further in the coming decade thanks in particular to the launch of nanosatellites constellations. These data should open new avenues for Earth surface monitoring due to highly improved spectral, spatial and temporal resolution. Many applications depend, however, on the accuracy of the image geometric model. The geometry of optical images, whether acquired from pushbroom or frame systems, is now commonly represented using a Rational Function Model (RFM). While the formalism has become standard, the procedures used to generate these models and their accuracies are diverse. As a result, the RFM models delivered with commercial data are commonly not accurate enough for 3-D extraction, subpixel registration or ground deformation measurements. In this study, we present a methodology for RFM optimization and demonstrate its potential for 3D reconstruction using tri-stereo and multi-date Cubesat images provided by SkySat and PlanetScope, respectively. We use SkySat data over the Morenci Mine, Arizona, which is the largest copper mine in the United States. The re-projection error after the RFM refinement is 0.42 pix without using ground control points (GCPs). Comparison of our Digital Elevation Model (DEM with ~3 m GSD) with a reference DEM obtained from an airborne LiDAR survey (with ~1 m GSD) over stable areas yields a standard deviation of the elevation differences of ~3.9 m. The comparison of the two DEMs allows detecting and measuring the topographic changes due to the mine activity (excavation and stockpiles). We assess the potential of PlanetScope data, using multi-date DOVE-C images from the Shisper glacier, located in the Karakoram (Pakistan), which is known for its recent surge. We extracted DEMs in 2017 and 2019 before and after the surge. The re-projection error after the RFM refinement is 0.38 pix without using GCPs. The accuracy of our DEMs (with ~9 m GSD) is evaluated through comparison with the SRTM DEM (GSD ~30 m) and with a DEM (GSD ~2 m) calculated from Geoeye-1 (GE-1) and World-View-2 (WV-2) stereo images. The standard deviation of the elevation differences in stable areas between the PlanetScope DEM and SRTM is ~12 m, and ~7 m with the GE-1&WV-2 DEM. The mass transfer due to the surge is clearly revealed from a comparison of the 2017 and 2019 DEMs. The study demonstrates that, with the proposed scheme for RFM optimization, times series of DEM extracted from SkySat and PlanetScope images can be used to measure topographic changes due to mining activities or ice flow, and could also be used to monitor geomorphic processes such as landslides, or coastal erosion for example

Optimization of Optical Image Geometric Modeling, Application to Topography Extraction and Topographic Change Measurements Using PlanetScope and SkySat Imagery

The volume of data generated by earth observation satellites has increased tremendously over the last few decades and will increase further in the coming decade thanks in particular to the launch of nanosatellites constellations. These data should open new avenues for Earth surface monitoring due to highly improved spectral, spatial and temporal resolution. Many applications depend, however, on the accuracy of the image geometric model. The geometry of optical images, whether acquired from pushbroom or frame systems, is now commonly represented using a Rational Function Model (RFM). While the formalism has become standard, the procedures used to generate these models and their accuracies are diverse. As a result, the RFM models delivered with commercial data are commonly not accurate enough for 3-D extraction, subpixel registration or ground deformation measurements. In this study, we present a methodology for RFM optimization and demonstrate its potential for 3D reconstruction using tri-stereo and multi-date Cubesat images provided by SkySat and PlanetScope, respectively. We use SkySat data over the Morenci Mine, Arizona, which is the largest copper mine in the United States. The re-projection error after the RFM refinement is 0.42 pix without using ground control points (GCPs). Comparison of our Digital Elevation Model (DEM with ~3 m GSD) with a reference DEM obtained from an airborne LiDAR survey (with ~1 m GSD) over stable areas yields a standard deviation of the elevation differences of ~3.9 m. The comparison of the two DEMs allows detecting and measuring the topographic changes due to the mine activity (excavation and stockpiles). We assess the potential of PlanetScope data, using multi-date DOVE-C images from the Shisper glacier, located in the Karakoram (Pakistan), which is known for its recent surge. We extracted DEMs in 2017 and 2019 before and after the surge. The re-projection error after the RFM refinement is 0.38 pix without using GCPs. The accuracy of our DEMs (with ~9 m GSD) is evaluated through comparison with the SRTM DEM (GSD ~30 m) and with a DEM (GSD ~2 m) calculated from Geoeye-1 (GE-1) and World-View-2 (WV-2) stereo images. The standard deviation of the elevation differences in stable areas between the PlanetScope DEM and SRTM is ~12 m, and ~7 m with the GE-1&WV-2 DEM. The mass transfer due to the surge is clearly revealed from a comparison of the 2017 and 2019 DEMs. The study demonstrates that, with the proposed scheme for RFM optimization, times series of DEM extracted from SkySat and PlanetScope images can be used to measure topographic changes due to mining activities or ice flow, and could also be used to monitor geomorphic processes such as landslides, or coastal erosion for example

Accuracy comparison of Pléiades satellite ortho-images using GPS device based GCPs against TerraSAR-X-based GCPs

Conducting single frame orthorectification on satellite images to create an ortho-image requires four basic components, namely an image, a geometric sensor model, elevation data (for example a digital elevation model (DEM)) and ground control points (GCPs). For this study, orthorectification entailed the use of a single scene Pléiades primary panchromatic image, applying the Pléiades rigorous geometric model, utilising a high-quality 2 m DEM and using GCPs that were acquired from two different collection methods. The application of these different GCPs to the execution of orthorectification encompassed the aim of this paper, which was to investigate and compare the positional accuracies of ortho-images under two scenarios. Firstly, GCPs were manually collected through fieldwork utilising a Trimble GeoExplorer 6000 series handheld GPS device and secondly, by utilising TerraSAR-X based GCPs that were acquired from Airbus Defence and Space. The objective of this study was to determine the geolocation accuracy of a high-resolution satellite ortho-image when different types of ground control are used. This required the execution of two orthorectification tests where only the type of GCPs differed. The results of these tests were interesting since it highlighted the difference in positional accuracy when utilising various sources of ground control to perform orthorectification on satellite imagery. The comparison results showed that utilising the manual GCPs produced a better positional accurate ortho-image as opposed to using the TerraSAR-X based GCPs. Nonetheless, the TerraSAR-X based GCPs still produced a sub 2 m accurate ortho-image, which is more than sufficient for the production of most geospatial products.Keywords: orthorectification, digital elevation model (DEM), ground control point (GCP), high-resolution satellite imagery, TerraSAR-X based GCPs, WorldDEM™, Airbus Defence and Spac

Outlet Glacier Dynamics in East Greenland and East Antarctica

Ice mass from the interior of Greenland and Antarctica is transported to the ocean by numerous large, fast-flowing outlet glaciers. Changes in the flow configuration of these outlet glaciers modulate ice sheet mass balance and sea level. Several recent studies have highlighted rapid increases in glacier speed in both Greenland and Antarctica, implying that the near-term contribution to sea level from ice sheets is under-estimated by current models. Here, the mass balance and force budget of several large outlet glaciers in East Greenland and East Antarctica are investigated using remote-sensing and field-based measurements. Recent estimates show that Greenland’s contribution to sea level more than doubled in the past decade, and that the majority of this additional mass loss is due to changes in the dynamics of a few large outlet glaciers. Our measurements indicate that up to ~10% of global sea level rise over the period 2001 – 2006 was contributed by just two glaciers, Helheim and Kangerdlugssuaq, in Southeast Greenland. We also find a latitudinal pattern of glacier behavior in East Greenland, where large and rapid changes are taking place south of 70°N while glaciers north of 70°N are stable. The East Antarctic Ice Sheet is Earth’s largest source of freshwater and has the potential to raise sea level by 57 m. The dynamics of outlet glaciers draining the ice sheet through the Transantarctic Mountains are largely unknown, but the glaciers are often assumed to be stable. In this study we investigate the dynamics of four large East Antarctic outlet glaciers. Together, these glaciers drain ~1,500,000 km2, or 12% by area of the entire Antarctic Ice Sheet. Mass balance calculations show modest imbalances for some glaciers, and a large imbalance for Byrd Glacier. Observations indicate a possible recent increase in flow speed, but this is insufficient to explain the large imbalance. We argue that catchment–wide estimates of accumulation rate contain large errors. This research provides new insights into the dynamic character of ice sheet outlet glaciers. In addition to quantifying recent changes, it also provides baseline data against which future behavior can be assessed

A mixed spaceborne sensor approach for surface modelling of an urban scene

Three-dimensional (3D) surface models are vital for sustainable urban management studies, and there is a nearly unlimited range of possible applications. Along-or across-track pairs from the same set of sensor imagery may not always be available or economical for a certain study area. Therefore, a photogrammetric approach is proposed in which a digital surface model (DSM) is extracted from a stereo pair of satellite images, acquired by different sensors. The results demonstrate that a mixed-sensor approach may offer a sound alternative to the more established along-track pairs. However, one should consider several criteria when selecting a suitable stereo pair. Two cloud-free acquisitions are selected from the IKONOS and QuickBird image archives, characterized by sufficient overlap and optimal stereo constellation in terms of complementarity of the azimuth and elevation angles. A densely built-up area in Istanbul, Turkey, covering 151 km(2) and with elevations ranging between sea level and approximately 160 m is presented as the test site. In addition to the general complexity of modelling the surface and elevation of an urban environment, multi-sensor image fusion has other particular difficulties. As the images are acquired from a different orbital pass, at a different date or instant and by a different sensor system, radiometric and geometric dissimilarities can occur, which may hamper the image-matching process. Strategies are presented for radiometric and geometric normalization of the multi-temporal and multi-sensor imagery and to deal with the differences in sensor characteristics. The accuracy of the generated surface model is assessed in comparison with 3D reference points, 3D rooftop vector data and surface models extracted from an along-track IKONOS stereo pair and an IKONOS triplet. When compared with a set of 35 reference GPS check points, the produced mixed-sensor model yields accuracies of 1.22, 1.53 and 2.96 m for the X, Y and Z coordinates, respectively, expressed in terms of root mean square errors (RMSEs). The results show that it is feasible to extract the DSM of a highly urbanized area from a mixed-sensor pair, with accuracies comparable with those observed from the DSM extracted from an along-track pair. Hence, the flexibility of reconstructing valuable elevation models is greatly increased by considering the mixed-sensor approach

Assessment of the CORONA series of satellite imagery for landscape archaeology: a case study from the Orontes valley, Syria

In 1995, a large database of satellite imagery with worldwide coverage taken from 1960 until 1972 was declassified. The main advantages of this imagery known as CORONA that made it attractive for archaeology were its moderate cost and its historical value. The main disadvantages were its unknown quality, format, geometry and the limited base of known applications. This thesis has sought to explore the properties and potential of CORONA imagery and thus enhance its value for applications in landscape archaeology. In order to ground these investigations in a real dataset, the properties and characteristics of CORONA imagery were explored through the case study of a landscape archaeology project working in the Orontes Valley, Syria. Present-day high-resolution IKONOS imagery was integrated within the study and assessed alongside CORONA imagery. The combination of these two image datasets was shown to provide a powerful set of tools for investigating past archaeological landscape in the Middle East. The imagery was assessed qualitatively through photointerpretation for its ability to detect archaeological remains, and quantitatively through the extraction of height information after the creation of stereomodels. The imagery was also assessed spectrally through fieldwork and spectroradiometric analysis, and for its Multiple View Angle (MVA) capability through visual and statistical analysis. Landscape archaeology requires a variety of data to be gathered from a large area, in an effective and inexpensive way. This study demonstrates an effective methodology for the deployment of CORONA and IKONOS imagery and raises a number of technical points of which the archaeological researcher community need to be aware of. Simultaneously, it identified certain limitations of the data and suggested solutions for the more effective exploitation of the strengths of CORONA imagery

Research Issues in Image Registration for Remote Sensing

Image registration is an important element in data processing for remote sensing with many applications and a wide range of solutions. Despite considerable investigation the field has not settled on a definitive solution for most applications and a number of questions remain open. This article looks at selected research issues by surveying the experience of operational satellite teams, application-specific requirements for Earth science, and our experiments in the evaluation of image registration algorithms with emphasis on the comparison of algorithms for subpixel accuracy. We conclude that remote sensing applications put particular demands on image registration algorithms to take into account domain-specific knowledge of geometric transformations and image content