Robust multi-scale image matching for deriving ice surface velocity field from sequential satellite images

2011 ◽  
Vol 33 (6) ◽  
pp. 1799-1822 ◽  
Author(s):  
Hongxing Liu ◽  
Lei Wang ◽  
Sheng-Jung Tang ◽  
Kenneth C. Jezek
2020 ◽  
Vol 66 (260) ◽  
pp. 905-915
Author(s):  
Bas Altena ◽  
Andreas Kääb

AbstractVelocities within an icefall are typically the fastest within a glacier system and experience complex flow. The combination of convergent and fast flow, and steep slope generate a quickly changing and intensely fractured surface. This complicates velocity extraction from repeat satellite images, especially when common pattern matching procedures are used. In this study, we exploit the high temporal revisit of medium-resolution satellite images using a novel image matching technique, ensemble matching, making it possible to generate a high-resolution (30 m) velocity field from high-repeat image sequences despite challenging image conditions. We demonstrate this technique for the first time in the glaciology domain using repeat Sentinel-2 optical data over the famous Khumbu icefall, situated on the southern slopes of Mount Everest. Estimates of velocity go just over 1 m d−1, which is slower than summer velocities from noisy single pair image matching. This icefall is frequently crossed by high-altitude mountaineers who use a route confined by fixed ropes and ladders set out every season. The mountain climbers typically record their trajectory on their personal satellite navigation device. We use such volunteered geographic information to verify our velocity estimates, confirming our underestimation with ensemble matching. Besides unprecedented remotely sensed surface velocities over the icefall, we also note that the generated velocity field can aid with the planning of a safe passage through this icefall.


2010 ◽  
Vol 4 (2) ◽  
pp. 161-178 ◽  
Author(s):  
T. Haug ◽  
A. Kääb ◽  
P. Skvarca

Abstract. We investigate the velocity field of the Larsen C ice shelf, Antarctic Peninsula, over the periods 2002–2006 and 2006–2009 based on repeat optical satellite data. The velocity field of the entire ice shelf is measured using repeat low resolution MODIS data (250 m spatial resolution). The measurements are validated for two ice shelf sections against repeat medium resolution Landsat 7 ETM + pan data (15 m spatial resolution). Horizontal surface velocities are obtained through image matching using both orientation correlation operated in the frequency domain and normalized crosscorrelation operated in the spatial domain, and the two methods compared. The uncertainty in the displacement measurements turns out to be about one fourth of the pixel size for the MODIS derived data, and about one pixel for the Landsat derived data. The difference between MODIS and Landsat based speeds is −15.4 m a−1 and 13.0 m a−1, respectively, for the first period for the two different validation sections on the ice shelf, and −26.7 m a−1 and 27.9 m a−1 for the second period for the same sections. This leads us to conclude that repeat MODIS images are well suited to measure ice shelf velocity fields and monitor their changes over time. Orientation correlation seems better suited for this purpose because it produces fewer mismatches, is able to match images with regular noise and data voids, and is faster. Since it can match images with regular data voids it is possible to match Landsat 7 ETM+ images even after the 2003 failure of the Scan Line Corrector (SLC off) that leaves significant image stripes with no data. Image matching based on the original 12-bit radiometric resolution MODIS data produced slightly better results than using the 8-bit version of the same images. Streamline interpolation from the obtained surface velocity field on Larsen C indicates ice travel times of up to 450 to 550 years between the inland boundary and the ice shelf edge. In a second step of the study we test our method successfully on 10 other ice shelves around Antarctica demonstrating that the approach presented could in fact be used for large scale monitoring of ice shelf dynamics.


2019 ◽  
Vol 60 (78) ◽  
pp. 1-7
Author(s):  
Hao Ke ◽  
Yuande Yang ◽  
Fei Li ◽  
Zemin Wang ◽  
Bo Sun ◽  
...  

ABSTRACTIn January 2016, static GPS measurements were carried out in a 30 × 30 km2 area centered around Kunlun station at Dome Argus (Dome A), East Antarctica, to acquire high-precision 3-D geodetic coordinates at 49 sites. By comparing the coordinates with previous GPS measurements in 2008 and 2013 at the same sites, we constructed a detailed and long-term record of the ice-surface velocity field, 2008–2016, around Dome A. During this time span, the estimated ice-surface velocity ranges from 0.8 ± 0.3 to 28.7 ± 1.6 cm a−1, with a mean of 10.4 ± 0.3 cm a−1. From 2013 to 2016, the surface elevation of most Dome A areas exhibits a rising trend, and the maximum increase of snow surface elevation is 84.8 cm. The mean snow surface elevation change rate at Dome A is estimated to be 6.6 ± 0.7 cm a−1. The difference of 1.0 cm a−1 between the snow surface change rate derived from GPS and pole-height change rate from surface mass balance is suspected to be a result of a combination of firn densification and basal melt under Dome A.


2010 ◽  
Vol 4 (1) ◽  
pp. 31-75 ◽  
Author(s):  
T. Haug ◽  
A. Kääb ◽  
P. Skvarca

Abstract. We investigate the velocity field of the Larsen C ice shelf, Antarctic Peninsula, over the periods 2002–2006 and 2006–2009 based on repeat optical satellite data. The velocity field of the entire ice shelf is measured using repeat low resolution MODIS data (250 m spatial resolution). The measurements are validated for two ice shelf sections against repeat medium resolution Landsat 7 ETM+ pan data (15 m spatial resolution). Horizontal surface velocities are obtained through image matching in both frequency and spatial domain, and the two methods compared. The uncertainty in the displacement measurements turns out to be less than 70 m for the MODIS derived data, and less than 15 m for the Landsat derived ones. The difference between MODIS and Landsat based speeds is −15.4 m a−1 and 13.0 m a−1, respectively, for the first period for the two different validation sections on the ice shelf, and −26.7 m a−1 and 27.9 m a−1 for the second period for the same sections. This leads us to conclude that repeat MODIS images are well suited to measure ice shelf velocity fields and monitor their changes over time. The frequency domain image correlation method seems better suited for this purpose because it is faster, produces fewer mismatches, and is able to match images with regular noise and data voids. The latter makes it possible to match Landsat 7 ETM+ images even after the 2003 failure of the Scan Line Corrector (SLC off) that leaves significant image sections with no data. Image matching based on the original 12-bit radiometric resolution MODIS data produced slightly better results than using the 8-bit version of the same images. Streamline interpolation from the obtained surface velocity field on Larsen C indicates ice travel times of up to 450 to 550 a between the inland boundary and the ice shelf edge. In a second step of the study we test our method successfully on 10 other ice shelves around Antarctica demonstrating that the approach presented could in fact be used for large scale monitoring of ice shelf dynamics.


2021 ◽  
Vol 13 (12) ◽  
pp. 2425
Author(s):  
Yiheng Cai ◽  
Dan Liu ◽  
Jin Xie ◽  
Jingxian Yang ◽  
Xiangbin Cui ◽  
...  

Analyzing the surface and bedrock locations in radar imagery enables the computation of ice sheet thickness, which is important for the study of ice sheets, their volume and how they may contribute to global climate change. However, the traditional handcrafted methods cannot quickly provide quantitative, objective and reliable extraction of information from radargrams. Most traditional handcrafted methods, designed to detect ice-surface and ice-bed layers from ice sheet radargrams, require complex human involvement and are difficult to apply to large datasets, while deep learning methods can obtain better results in a generalized way. In this study, an end-to-end multi-scale attention network (MsANet) is proposed to realize the estimation and reconstruction of layers in sequences of ice sheet radar tomographic images. First, we use an improved 3D convolutional network, C3D-M, whose first full connection layer is replaced by a convolution unit to better maintain the spatial relativity of ice layer features, as the backbone. Then, an adjustable multi-scale module uses different scale filters to learn scale information to enhance the feature extraction capabilities of the network. Finally, an attention module extended to 3D space removes a redundant bottleneck unit to better fuse and refine ice layer features. Radar sequential images collected by the Center of Remote Sensing of Ice Sheets in 2014 are used as training and testing data. Compared with state-of-the-art deep learning methods, the MsANet shows a 10% reduction (2.14 pixels) on the measurement of average mean absolute column-wise error for detecting the ice-surface and ice-bottom layers, runs faster and uses approximately 12 million fewer parameters.


2013 ◽  
Vol 54 (63) ◽  
pp. 221-228 ◽  
Author(s):  
James Turrin ◽  
Richard R. Forster ◽  
Chris Larsen ◽  
Jeanne Sauber

AbstractBering Glacier, Alaska, USA, has a ∼20 year surge cycle, with its most recent surge reaching the terminus in 2011. To study this most recent activity a time series of ice velocity maps was produced by applying optical feature-tracking methods to Landsat-7 ETM+ imagery spanning 2001-11. The velocity maps show a yearly increase in ice surface velocity associated with the down-glacier movement of a surge front. In 2008/09 the maximum ice surface velocity was 1.5 ±0.017 km a-1 in the mid-ablation zone, which decreased to 1.2 ±0.015 km a-1 in 2009/10 in the lower ablation zone, and then increased to nearly 4.4 ± 0.03 km a-1 in summer 2011 when the surge front reached the glacier terminus. The surge front propagated down-glacier as a kinematic wave at an average rate of 4.4 ±2.0 km a-1 between September 2002 and April 2009, then accelerated to 13.9 ± 2.0 km a-1 as it entered the piedmont lobe between April 2009 and September 2010. The wave seems to have initiated near the confluence of Bering Glacier and Bagley Ice Valley as early as 2001, and the surge was triggered in 2008 further down-glacier in the mid-ablation zone after the wave passed an ice reservoir area.


2014 ◽  
Vol 2 (2) ◽  
pp. 911-933 ◽  
Author(s):  
N. F. Glasser ◽  
S. J. A. Jennings ◽  
M. J. Hambrey ◽  
B. Hubbard

Abstract. Continent-wide mapping of longitudinal ice-surface structures on the Antarctic Ice Sheet reveals that they originate in the interior of the ice sheet and are arranged in arborescent networks fed by multiple tributaries. Longitudinal ice-surface structures can be traced continuously down-ice for distances of up to 1200 km. They are co-located with fast-flowing glaciers and ice streams that are dominated by basal sliding rates above tens of m yr-1 and are strongly guided by subglacial topography. Longitudinal ice-surface structures dominate regions of converging flow, where ice flow is subject to non-coaxial strain and simple shear. Associating these structures with the AIS' surface velocity field reveals (i) ice residence times of ~ 2500 to 18 500 years, and (ii) undeformed flow-line sets for all major flow units analysed except the Kamb Ice Stream and the Institute and Möller Ice Stream areas. Although it is unclear how long it takes for these features to form and decay, we infer that the major ice-flow and ice-velocity configuration of the ice sheet may have remained largely unchanged for several thousand years, and possibly even since the end of the last glacial cycle. This conclusion has implications for our understanding of the long-term landscape evolution of Antarctica, including large-scale patterns of glacial erosion and deposition.


2019 ◽  
Author(s):  
Wentao Zhu ◽  
Yufang Huang ◽  
Mani A Vannan ◽  
Shizhen Liu ◽  
Daguang Xu ◽  
...  

AbstractEchocardiography has become routinely used in the diagnosis of cardiomyopathy and abnormal cardiac blood flow. However, manually measuring myocardial motion and cardiac blood flow from echocar-diogram is time-consuming and error-prone. Computer algorithms that can automatically track and quantify myocardial motion and cardiac blood flow are highly sought after, but have not been very successful due to noise and high variability of echocardiography. In this work, we propose a neural multi-scale self-supervised registration (NMSR) method for automated myocardial and cardiac blood flow dense tracking. NMSR incorporates two novel components: 1) utilizing a deep neural net to parameterize the velocity field between two image frames, and 2) optimizing the parameters of the neural net in a sequential multi-scale fashion to account for large variations within the velocity field. Experiments demonstrate that NMSR yields significantly better registration accuracy than the state-of-the-art methods, such as advanced normalization tools (ANTs) and Voxel Morph, for both myocardial and cardiac blood flow dense tracking. Our approach promises to provide a fully automated method for fast and accurate analyses of echocardiograms.


2011 ◽  
Vol 5 (1) ◽  
pp. 299-313 ◽  
Author(s):  
G. E. Flowers ◽  
N. Roux ◽  
S. Pimentel ◽  
C. G. Schoof

Abstract. Glacier surges are a well-known example of an internal dynamic oscillation whose occurrence is not a direct response to the external climate forcing, but whose character (i.e. period, amplitude, mechanism) may depend on the glacier's environmental or climate setting. We examine the dynamics of a small (∼5 km2) valley glacier in Yukon, Canada, where two previous surges have been photographically documented and an unusually slow surge is currently underway. To characterize the dynamics of the present surge, and to speculate on the future of this glacier, we employ a higher-order flowband model of ice dynamics with a regularized Coulomb-friction sliding law in both diagnostic and prognostic simulations. Diagnostic (force balance) calculations capture the measured ice-surface velocity profile only when non-zero basal water pressures are prescribed over the central region of the glacier, coincident with where evidence of the surge has been identified. This leads to sliding accounting for 50–100% of the total surface motion in this region. Prognostic simulations, where the glacier geometry evolves in response to a prescribed surface mass balance, reveal a significant role played by a bedrock ridge beneath the current equilibrium line of the glacier. Ice thickening occurs above the ridge in our simulations, until the net mass balance reaches sufficiently negative values. We suggest that the bedrock ridge may contribute to the propensity for surges in this glacier by promoting the development of the reservoir area during quiescence, and may permit surges to occur under more negative balance conditions than would otherwise be possible. Collectively, these results corroborate our interpretation of the current glacier flow regime as indicative of a slow surge that has been ongoing for some time, and support a relationship between surge incidence or character and the net mass balance. Our results also highlight the importance of glacier bed topography in controlling ice dynamics, as observed in many other glacier systems.


Sign in / Sign up

Export Citation Format

Share Document