Search

Glacier movement is a crucial factor for assessing cryospheric climate change. Traditional methods of field surveys for studying glacier movement and velocity are often not possible owing to inaccessibility and harsh terrains. Furthermore, as it is not feasible to physically monitor and survey many glaciers around the globe, these traditional methods are limited in their global coverage. Remote sensing is an ideal tool to study such phenomena on a global scale. Optical remote sensing employs techniques such as feature tracking and pixel tracking, whereas, microwave remote sensing uses intensity tracking, speckle tracking, Interferometric SAR and Differential InSAR (DInSAR). This study focuses on estimation of glacier velocity and its seasonal variations using the image-matching technique for optical images for the Fisher glacier, a tributary glacier of the Amery ice shelf in Antarctica. The tool used in this study is the COSI-Corr module embedded in ENVI which provides the velocity in both azimuth and range resolution. The principle of estimating velocity using this tool is pixel tracking wherein similar pixels on two images are tracked where one is the master image and the other is a slave. This technique correlates the master and slave images over a time period and generates three outputs: displacements in the East-West and North-South directions and signal-to-noise ratio (SNR) image. Landsat 8 image pairs were used for cross correlation over a time span of four years spanning 2013–2017. With a resolution of 15 m, the panchromatic band (Band 8) was the ideal choice for pixel tracking as the resolution of other bands is coarser. The initial window size for correlation was 64 while the final window size was 16. The resolution of the displacement images produced is dependent on the value assigned for the step size, which was set to 8. The resultant velocity was derived using the result of the two displacement images. The SNR image was used to remove all the pixels from the velocity output having the value of SNR less than 0.9, in order to reduce the effect of noise. The annual velocity of the Fisher glacier was estimated to be around 600 to 650 myr−1 near the tongue where it merges with the Amery Ice Shelf, which was reduced to 150 myr−1 as it recedes. The resultant velocity images have a resolution of 120 m. The seasonal variation in velocity for the year 2013–2014 is 1.8 myr−1, while in the succeeding year 2014-2015 it subdued to 1.7 myr−1. The seasonal variation for the year 2015–2016 was estimated to be 7.9 myr−1. The seasonal variation for 2016–2017 was 17.4 myr−1.

This work uses the canopy height model (CHM) based workflow for individual tree crown delineation and 3D feature extraction approach (Overwatch Geospatial's proprietary algorithm) for building feature delineation from high-density light detection and ranging (LiDAR) point cloud data in an urban environment and evaluates its accuracy by using very high-resolution panchromatic (PAN) (spatial) and 8-band (multispectral) WorldView-2 (WV-2) imagery. LiDAR point cloud data over San Francisco, California, USA, recorded in June 2010, was used to detect tree and building features by classifying point elevation values. The workflow employed includes resampling of LiDAR point cloud to generate a raster surface or digital terrain model (DTM), generation of a hill-shade image and an intensity image, extraction of digital surface model, generation of bare earth digital elevation model (DEM) and extraction of tree and building features. First, the optical WV-2 data and the LiDAR intensity image were co-registered using ground control points (GCPs). The WV-2 rational polynomial coefficients model (RPC) was executed in ERDAS Leica Photogrammetry Suite (LPS) using supplementary *.RPB file. In the second stage, ortho-rectification was carried out using ERDAS LPS by incorporating well-distributed GCPs. The root mean square error (RMSE) for the WV-2 was estimated to be 0.25 m by using more than 10 well-distributed GCPs. In the second stage, we generated the bare earth DEM from LiDAR point cloud data. In most of the cases, bare earth DEM does not represent true ground elevation. Hence, the model was edited to get the most accurate DEM/ DTM possible and normalized the LiDAR point cloud data based on DTM in order to reduce the effect of undulating terrain. We normalized the vegetation point cloud values by subtracting the ground points (DEM) from the LiDAR point cloud. A normalized digital surface model (nDSM) or CHM was calculated from the LiDAR data by subtracting the DEM from the DSM. The CHM or the normalized DSM represents the absolute height of all aboveground urban features relative to the ground. After normalization, the elevation value of a point indicates the height from the ground to the point. The above-ground points were used for tree feature and building footprint extraction. In individual tree extraction, first and last return point clouds were used along with the bare earth and building footprint models discussed above. In this study, scene dependent extraction criteria were employed to improve the 3D feature extraction process. LiDAR-based refining/ filtering techniques used for bare earth layer extraction were crucial for improving the subsequent 3D features (tree and building) feature extraction. The PAN-sharpened WV-2 image (with 0.5 m spatial resolution) was used to assess the accuracy of LiDAR-based 3D feature extraction. Our analysis provided an accuracy of 98 % for tree feature extraction and 96 % for building feature extraction from LiDAR data. This study could extract total of 15143 tree features using CHM method, out of which total of 14841 were visually interpreted on PAN-sharpened WV-2 image data. The extracted tree features included both shadowed (total 13830) and non-shadowed (total 1011). We note that CHM method could overestimate total of 302 tree features, which were not observed on the WV-2 image. One of the potential sources for tree feature overestimation was observed in case of those tree features which were adjacent to buildings. In case of building feature extraction, the algorithm could extract total of 6117 building features which were interpreted on WV-2 image, even capturing buildings under the trees (total 605) and buildings under shadow (total 112). Overestimation of tree and building features was observed to be limiting factor in 3D feature extraction process. This is due to the incorrect filtering of point cloud in these areas. One of the potential sources of overestimation was the man-made structures, including skyscrapers and bridges, which were confounded and extracted as buildings. This can be attributed to low point density at building edges and on flat roofs or occlusions due to which LiDAR cannot give as much precise planimetric accuracy as photogrammetric techniques (in segmentation) and lack of optimum use of textural information as well as contextual information (especially at walls which are away from roof) in automatic extraction algorithm. In addition, there were no separate classes for bridges or the features lying inside the water and multiple water height levels were also not considered. Based on these inferences, we conclude that the LiDAR-based 3D feature extraction supplemented by high resolution satellite data is a potential application which can be used for understanding and characterization of urban setup.

Available digital elevation models (DEMs) of Antarctic region generated by using radar altimetry and the Antarctic digital database (ADD) indicate elevation variations of up to hundreds of meters, which necessitates the generation of local DEM and its validation by using ground reference. An enhanced digital elevation model (eDEM) of the Schirmacher oasis region, east Antarctica, is generated synergistically by using Cartosat-1 stereo pair-derived photogrammetric DEM (CartoDEM)-based point elevation dataset and multitemporal radarsat Antarctic mapping project version 2 (RAMPv2) DEM-based point elevation dataset. In this study, we analyzed suite of interpolation techniques for constructing a DEM from RAMPv2 and CartoDEM-based point elevation datasets, in order to determine the level of confidence with which the interpolation techniques can generate a better interpolated continuous surface, and eventually improves the elevation accuracy of DEM from synergistically fused RAMPv2 and CartoDEM point elevation datasets. RAMPv2 points and CartoDEM points were used as primary data for various interpolation techniques such as ordinary kriging (OK), simple kriging (SK), universal kriging (UK), disjunctive kriging (DK) techniques, inverse distance weighted (IDW), global polynomial (GP) with power 1 and 2, local polynomial (LP) and radial basis functions (RBF). Cokriging of 2 variables with second dataset was used for ordinary cokriging (OCoK), simple cokriging (SCoK), universal cokriging (UCoK) and disjunctive cokriging (DCoK). The IDW, GP, LP, RBF, and kriging methods were applied to one variable, while Cokriging experiments were employed on two variables. The experiment of dataset and its combination produced two types of point elevation map categorized as (1) one variable (RAMPv2 Point maps and CartoDEM Point maps) and (2) two variables (RAMPv2 Point maps + CartoDEM Point maps). Interpolated surfaces were evaluated with the help of differential global positioning system (DGPS) points collected from study area during the Indian Scientific Expedition to Antarctic (ISEA). Accuracy assessment of the RAMPv2 DEM, CartoDEM, and combined eDEM (RAMPv2 + CartoDEM) by using DGPS as ground reference data shows that eDEM achieves much better accuracy (average elevation error 8.44 m) than that of existing DEM constructed by using only CartoDEM (13.57 m) or RAMPv2 (41.44 m) alone. The newly constructed eDEM achieves a vertical accuracy of about 7 times better than RAMPv2 DEM and 1.5 times better than CartoDEM. After using accurate DGPS data for accuracy assessment, the approximation to the actual surface of the eDEM extracted here is much more accurate with least mean root mean square error (RMSE) of 9.22 m than that constructed by using only CartoDEM (RMSE = 14.15 m) point elevation data and RAMPv2 (RMSE = 69.48 m) point elevation data. Our results indicate that, the overall trend of accuracy for the interpolation methods for generating continuous elevation surface from CartoDEM + RAMPv2 point elevation data, based on RMSE, is as follows: GP1 > IDW > GP2 > OK > LP2 > DK > LP1 > RBF > SK > UK. In case of cokriging interpolation methods, OCoK yields more accurate eDEM with the least RMSE of 8.16 m, which can be utilized to generate a highly accurate DEM of the research area.. Based on this work, it is inferred that GP2 and OCok interpolation methods and synergistic use of RAMPv2 and CartoDEM-based point elevation datasets lead to a highly accurate DEM of the study region. This research experiment demonstrates the stability (w.r.t multi-temporal datasets), performance (w.r.t best interpolation technique) and consistency (w.r.t all the experimented interpolation techniques) of synergistically fused eDEM. On the basis of average elevation difference and RMSE mentioned in present research, the newly constructed eDEM may serve as a benchmark for future elevation models such as from the ICESAT-II mission to spatially monitor ice sheet elevation.

This study presents an operational case of advancements in urban land cover classification and change detection by using very high resolution spatial and multispectral information from 4-band QuickBird (QB) and 8-band WorldView-2 (WV-2) image sequence. Our study accentuates quantitative, pixel based, image difference approach for operational change detection using very high resolution pansharpened QB and WV-2 images captured over San Francisco city, California, USA (37° 44" 30N', 122° 31" 30' W and 37° 41" 30° N ,122° 20" 30' W). In addition to standard QB image, we compiled three multiband images from eight pansharpened WV-2 bands: (1) multiband image from four traditional spectral bands, i.e., Blue, Green, Red and near-infrared 1 (NIR1) (henceforth referred as "QB equivalent WV-2"), (2) multiband image from four new spectral bands, i.e., Coastal, Yellow, Red Edge and NIR2 (henceforth referred as "new band WV-2"), and (3) multiband image consisting of four traditional and four new bands (henceforth referred as "standard WV-2"). All the four multiband images were classified using support vector machine (SVM) classifier into four most abundant land cover classes, viz, hard surface, vegetation, water and shadow. The assessment of classification accuracy was performed using random selection of 356 test points. Land cover classifications on "standard QB" image (kappa coeffiecient, κ = 0.93), "QB equivalent WV-2" image (κ = 0.97), and "new band WV-2" image (κ = 0.97) yielded overall accuracies of 96.31 %, 98.03 % and 98.31 %, respectively, while "standard WV-2" image (κ = 0.99) yielded an improved overall accuracy of 99.18 %. It is concluded that the addition of four new spectral bands to the existing four traditional bands improved the discrimination of land cover targets, due to increase in the spectral characteristics of WV-2 satellite. Consequently, to test the validity of improvement in classification process for implementation in operational change detection application, comparative assessment of transition of various landcover classes in three WV-2 images with respect to "standard QB" image was carried out using image difference method. As far as waterbody class is concerned, there was no significant transition observed in all the three WorldView-2 Images, whereas, hard surface class showed lowest transition in "standard WV-2" image and highest in case of "new band WV-2". The most significant transition was occurred in vegetation class in all of the three images, showing positive change (increase) in standard WV-2 image (0.31 Sq. Km) and negative change (decrease) in other two images (−0.12 Sq. Km for "QB equivalent WV-2" image and −31.15 Sq. Km in "new band WV-2" image) with considerable amount. Similar case was observed with the shadow class, but the difference is, transition from shadow to other classes was negative in all the three WV-2 images which can be attributed to the fact that, "standard QB" image had more shadow area (based on acquisition time and sun position) than WV-2, that means all the band combinations of WV-2 succeeded in extracting the features hidden below the shadow in "standard QB" image. These trends indicate that the overall bandwise transition in landcover classes in case of "standard WV-2" is more precise than other two images. We note that "QB equivalent WV-2" image had narrower band widths than those of "standard QB" image but the observed vegetation change is not prominent as in case of other two images, but at the same time, transition in hard surface and waterbody was discerned more efficiently than "new band WV-2" image. The addition of new bands in WV-2 enabled more effective vegetation analysis, so the vegetation transition results shown by "new band WV-2" image were at par with the "standard WV-2" image, showing the importance of these newly added bands in the WV-2 imagery, with comparatively lower transitions in other classes. In a nutshell, it can be claimed that incorporation of new bands along with even narrower Red, Green, Blue and Near Infrared-1 bands in WV-2 image holds remarkable importance which leads to enhancement in the potential of WV-2 imagery in change detection and other feature extraction studies.