Huabing Huang

Modelling spatial-temporal change of Poyang Lake using multitemporal Landsat imagery

Poyang Lake is a seasonal lake, exchanging water with the lower branch of the Yangtze River. During the spring and summer flooding season it inundates a large area while in the winter it shrinks considerably, creating a large tract of marshland for wild migratory birds. A better knowledge of the water coverage duration and the beginning and ending dates for the vast range of marshlands surrounding the lake is important for the measurement, modelling and management of marshland ecosystems. In addition, the abundance of a special type of snail (Oncomelania hupensis), the intermediate host of parasite schistosome (Schistosoma japonicum) in this region, is also heavily dependent on the water coverage information. However, there is no accurate digital elevation model (DEM) for the lake bottom and the inundated marshland, nor is there sufficient water level information over this area. In this study, we assess the feasibility of the use of multitemporal Landsat images for mapping the spatial‐temporal change of Poyang Lake water body and the temporal process of water inundation of marshlands. Eight cloud‐free Landsat Thematic Mapper images taken during a period of one year were used in this study. We used the normalized difference water index (NDWI) and the modified normalized difference water index (MNDWI) methods to map water bodies. We then examined the annual spatial‐temporal change of the Poyang Lake water body. Finally we attempted to obtain the duration of water inundation of marshlands based on the temporal sequence of water extent determined from the Landsat images. The results showed that although the images can be used to capture the snapshots of water coverage in this area, they are insufficient to provide accurate estimation of the spatial‐temporal process of water inundation over the marshlands through linear interpolation.

Spatio-Temporal Distribution of Malaria in Yunnan Province, China

The spatio-temporal distribution pattern of malaria in Yunnan Province, China was studied using a geographic information system technique. Both descriptive and temporal scan statistics revealed seasonal fluctuation in malaria incidences in Yunnan Province with only one peak during 1995-2000, and two apparent peaks from 2001 to 2005. Spatial autocorrelation analysis indicated that malaria incidence was not randomly distributed in the province. Further analysis using spatial scan statistics discovered that the high risk areas were mainly clustered at the bordering areas with Myanmar and Laos, and in Yuanjiang River Basin. There were obvious associations between Plasmodium vivax and Plasmodoium falciparum malaria incidences and climatic factors with a clear 1-month lagged effect, especially in cluster areas. All these could provide information on where and when malaria prevention and control measures would be applied. These findings imply that countermeasures should target high risk areas at suitable times, when climatic factors facilitate the transmission of malaria.

Earth science applications of ICESat/GLAS: a review

The Ice, Cloud, and Land Elevation Satellite (ICESat) completed 19 successful campaigns for Earth observation missions following its launch in 2003. The Geoscience Laser Altimeter System (GLAS) on board ICESat provided data of high quality with unprecedented accuracy over the globe. The three laser sensors of GLAS acquired a large volume of data between 2003 and 2009. These data were used widely to detect changes in Greenland and Antarctic ice sheets and to determine forest heights, sea-ice freeboard heights and the distribution of cloud and aerosols. Here, we provide a review of these applications, describe the methodology involved in GLAS data processing and summarize some of the challenges to make better use of GLAS data. Other applications, including ice-sheet slope extraction, distinguishing between water, bare land, urban building and high forest, urban building height extraction, changes in glaciers and ice caps and water levels in lakes are discussed more briefly.

Automated Methods for Measuring DBH and Tree Heights with a Commercial Scanning Lidar

Accurate forest structural parameters are crucial to forest inventory, and modeling of the carbon cycle and wildlife habitat. Lidar (Light Detection and Ranging) is particularly suitable to the measurement of forest structural parameters. In this paper, we describe a pilot study to extract forest structural parameters, such as tree height, diameter at breast height (DBH), and position of individual tree using a terrestrial lidar (LMS-Z360i; Riegel, Inc.). The lidar was operated to acquire both vertical and horizontal scanning in the field in order to obtain a point cloud of the whole scene. An Iterative Closet Point (ICP) algorithm was introduced to obtain the transformation matrix of each range image and to mosaic multiple range images together. Based on the mosaiced data set, a variable scale and threshold filtering method was used to separate ground from the vegetation. Meanwhile, a Digital Elevation Model (DEM) and a Canopy Height Model (CHM) were generated from the classified point cloud. A stem detection algorithm was used to extract the location of individual trees. A slice above 1.3 m from the ground was extracted and rasterized. A circle fitting algorithm combined with the Hough transform was used to retrieve the DBH based on the rasterized grid. Tree heights were calculated using the height difference between the minimum and maximum Z values within the position of each individual tree with a 1 m buffer. All of the 26 trees were detected correctly, tree height and DBH were determined with a precision of 0.76 m and 3.4 cm, respectively, comparing with those visually measured in the lidar data. Our methods and results confirm that terrestrial lidar can provide nondestructive, high-resolution, and automatic determination of parameters required in forest inventory.

Reduction of atmospheric and topographic effect on Landsat TM data for forest classification

The incident radiance in forested areas with rugged terrain varies greatly with the changes in solar elevation and azimuth, slope and aspect of the terrain, and the relative position of trees. The geotropic nature must be considered in the course of topographic correction. The Sun‐Canopy‐Sensor (SCS) model is introduced to substitute the cosine correction in a physical model. We used an atmospheric simulation code, MODTRAN, and a digital elevation model (DEM) to calculate the path radiance, downwards diffuse radiance and two‐way transmittance of direct and diffuse light at different altitudes. Based on the atmospheric parameters derived above and the Lambertian assumption, surface reflectance in a forested area was retrieved from Landsat Thematic Mapper (TM) imagery using a revised physical model. Meanwhile, a smoothed DEM was used to assess the effect of noise on the DEM and misregistration between the DEM and the satellite imagery. Correlation analysis, spectral comparison between sunlit and shaded slopes and a support vector machine (SVM) classification were performed to assess the effect of the revised radiometric correction algorithm. Results indicate that the revised physical model with smoothed DEM is more adequate for forested terrain and more consistent spectra for similar vegetation under different illuminations can be obtained. Finally, higher classification accuracy of forested land can be achieved with the revised correction algorithm compared with the SCS correction and the original physical correction model.

Improving measurement of forest structural parameters by co-registering of high resolution aerial imagery and low density LiDAR data.

Forest structural parameters, such as tree height and crown width, are indispensable for evaluating forest biomass or forest volume. LiDAR is a revolutionary technology for measurement of forest structural parameters, however, the accuracy of crown width extraction is not satisfactory when using a low density LiDAR, especially in high canopy cover forest. We used high resolution aerial imagery with a low density LiDAR system to overcome this shortcoming. A morphological filtering was used to generate a DEM (Digital Elevation Model) and a CHM (Canopy Height Model) from LiDAR data. The LiDAR camera image is matched to the aerial image with an automated keypoints search algorithm. As a result, a high registration accuracy of 0.5 pixels was obtained. A local maximum filter, watershed segmentation, and object-oriented image segmentation are used to obtain tree height and crown width. Results indicate that the camera data collected by the integrated LiDAR system plays an important role in registration with aerial imagery. The synthesis with aerial imagery increases the accuracy of forest structural parameter extraction when compared to only using the low density LiDAR data.

ICESat GLAS Data for Urban Environment Monitoring

Although the Geoscience Laser Altimeter System (GLAS) onboard the NASA Ice, Cloud and Land Elevation Satellite was not designed for urban applications, its 3-D measurement capability over the globe makes it a nice feature for consideration in monitoring urban heights. However, this has not been previously done. In this paper, we report some preliminary assessment of the GLAS data for building height and density estimation in a suburb of Beijing, China. Building heights can be directly calculated from a GLAS data product (GLA14). Because GLA14 limits height levels to six in each ground footprint, we developed a new method to remove this restriction by processing the raw GLAS data. The maximum heights measured in the field at selected GLAS footprints were used to validate the GLAS measurement results. By assuming a constant incident energy and surface reflectance within a GLAS footprint, the building density can be estimated from GLA14 or from our newly processed GLAS data. The building density determined from high-resolution images in Google Earth was used to validate the GLAS estimation results. The results indicate that the newly developed method can produce more accurate building height estimation within each GLAS footprint (R2 = 0.937, rmse = 6.4 m, and n = 26) than the GLA14 data product (R2 = 0.808, rmse = 11.5 m, and n = 26). However, satisfactory estimation results on building density cannot be obtained from the GLAS data with the methods investigated in this paper. Forest cover could be a challenge to building height and density estimation from the GLAS data. It should be addressed in future research.

A robust multisource image automatic registration system based on the SIFT descriptor

Image registration is an essential step in many remote-sensing (RS) applications. This article presents a study of a multisource image automatic registration system (MIARS) based on the scale-invariant feature transform (SIFT), which has been demonstrated to be the most robust local invariant feature descriptor for automatically registering various RS images. The SIFT descriptor has two shortcomings: it is unsuitable for extremely large images and has an irregular distribution of feature points. Therefore, three steps are proposed for the MIARS: image division, histogram equalization and the elimination of false point matches by a subregion least squares iteration. Image division makes it possible to use the SIFT descriptor to extract control points from an extremely large RS image. Histogram equalization in prematching improves the contrast sensitivity of RS images. The subregion least squares iteration refines the registration accuracy. Images from multisensor systems, including Quickbird, IRS-P6, Landsat/TM, HJ-CCD, HJ-IRS, light detection and ranging (LiDAR) intensity images and aerial data, were selected to test the reliability of the MIARS. The results indicated that better registration accuracy was achieved, which will be very helpful in the future development of a registration model.

Mapping blue-ice areas in Antarctica using ETM+ and MODIS data

Abstract Blue-ice areas (BIAs) and their geographical distribution in Antarctica were mapped using Landsat-7 ETM+ images with 15 m spatial resolution obtained during the 1999–2003 austral summers and covering the area north of 82.5° S, and a snow grain-size image of the MODIS-based Mosaic of Antarctica (MOA) dataset with 125 m grid spacing acquired during the 2003/04 austral summer from 82.5°S to the South Pole. A map of BIAs was created with algorithms of thresholds based on band ratio and reflectance for ETM+ data and thresholds based on snow grain size for the MOA dataset. The underlying principle is that blue ice can be separated from snow or rock by their spectral discrepancies and by different grain sizes of snow and ice. We estimate the total area of BIAs in Antarctica during the data acquisition period is 234 549 km2, or 1.67% of the area of the continent. Blue ice is scattered widely over the continent but is generally located in coastal or mountainous regions. The BIA dataset presented in this study is the first map covering the entire Antarctic continent sourced solely from ETM+ and MODIS data. This dataset can potentially benefit other studies in glaciology, meteorology, climatology and paleoclimate, meteorite collection and airstrip site selection.

Quantifying Multi-Decadal Change of Planted Forest Cover Using Airborne LiDAR and Landsat Imagery

Continuous monitoring of forest cover condition is key to understanding the carbon dynamics of forest ecosystems. This paper addresses how to integrate single-year airborne LiDAR and time-series Landsat imagery to derive forest cover change information. LiDAR data were used to extract forest cover at the sub-pixel level of Landsat for a single year, and the Landtrendr algorithm was applied to Landsat spectral data to explore the temporal information of forest cover change. Four different approaches were employed to model the relationship between forest cover and Landsat spectral data. The result shows incorporating the historic information using the temporal trajectory fitting process could infuse the model with better prediction power. Random forest modeling performs the best for quantitative forest cover estimation. Temporal trajectory fitting with random forest model shows the best agreement with validation data (R2 = 0.82 and RMSE = 5.19%). We applied our approach to Youyu county in Shanxi province of China, as part of the Three North Shelter Forest Program, to map multi-decadal forest cover dynamics. With the availability of global time-series Landsat imagery and affordable airborne LiDAR data, the approach we developed has the potential to derive large-scale forest cover dynamics.