NOAA-AVHRR data processing for the mapping of vegetation cover

The NOAA-AVHRR images have been widely used for global studies due to their low cost, suitable wavebands and high temporal resolution. Data from the AVHRR sensor (Bands 1 and 2) transformed to the Normalized Difference Vegetation Index (NDVI) are the most common product used in global land cover studies. The purpose of this Letter is to present the vegetation, soil, and shade fraction images derived from AVHRR, in addition to NDVI, to monitor land cover. Six AVHRR images from the period of 21 to 26 June 1993 were composed and used to obtain the above mentioned products over Sa o Paulo State, in the south-east of Brazil. Vegetation fraction component values were strongly correlated with NDVI values ( r 0.95; n 60). Also, the fraction image presented a good agreement with the available global vegetation map of Sao Paulo State derived from Landsat TM images.     

Crop type identification by integration of high-spatial resolution multispectral data with features extracted from coarse-resolution time-series vegetation index data

Crop type identification is the basis of crop acreage estimation and plays a key role in crop production prediction and food security analysis. However, the accuracy of crop type identification using remote-sensing data needs to be improved to support operational agriculture-monitoring tasks. In this paper, a new method integrating high-spatial resolution multispectral data with features extracted from coarse-resolution time-series vegetation index data is proposed to improve crop type identification accuracy in Hungary. Four crop growth features, including peak value, date of peak occurrence, average rate of green-up, and average rate for the senescence period were extracted from time-series Moderate Resolution Imaging Spectroradiometer (MODIS) normalized difference vegetation index (NDVI) profiles and spatially enhanced to 30 m resolution using resolution merge tools based on a multiplicative method to match the spatial resolution of Landsat Thematic Mapper (TM) data. A maximum likelihood classifier (MLC) was used to classify the TM and merged images. Independent validation results indicated that the average overall classification accuracy was improved from 92.38% using TM to 94.67% using the merged images. Based on the classification results using the proposed method, acreages of two major summer crops were estimated and compared to statistical data provided by the United States Department of Agriculture (USDA). The proposed method was able to achieve highly satisfactory crop type identification results.     

A technique for extrapolating and validating forest cover across large regions Calibrating AVHRR data with TM data

Abstract

An approach to extending high-resolution forest cover information across large regions is presented and validated. Landsat Thematic Mapper (TM) data were classified into forest and nonforest for a portion of Jackson County, Illinois. The classified TM image was then used to determine the relationship between forest cover and the spectral signature of Advanced Very High Resolution Radiometer (AVHRR) pixels covering the same location. Regression analysis was used to develop an empirical relationship between AVHRR spectral signatures and forest cover. The regression equation developed from data from the single county calibration area in southern Illinois was then applied to the entire AVHRR scene, which covered all or parts of ten states, to produce a regional map of forest cover. This map was used to derive estimates of forest cover, within a geographical information system (GIS), for each of the 428 counties located within the boundaries of the original AVHRR scene. The validity of the overall regional map was tested by comparing the AVHRR/TM-derived estimates of county forest cover with independent estimates of county forest cover developed by the U.S. Forest Service (USFS). The overall correlation coefficient of the AVHRR/TM and USFS county forest cover estimates was r=0-89 (n=428 counties). Not surpris0ingly, some individual states and the areas nearer to the southern Illinois calibration centre had higher correlation coefficients. Absolute estimates of forest cover percentages were also significantly well predicted. With the future inclusion of multiple calibration centres representing a number of physiographic regions, the method shows promise for predicting continental and global estimates of forest cover.     

Subpixel measurement of tropical forest cover using AVHRR data

Abstract

Tropical forest assessment using data from the Advanced Very High Resolution Radiometer (AVHRR) may lead to inaccurate estimates of forest cover in regions of small subpixel forest or non-forest patches and in regions where the pattern of clearance is particularly convoluted. Test sites typifying these two patterns were chosen in Ghana and Rondonia, respectively. To capture the subpixel proportions of forest cover, a linear mixture model was applied to two AVHRR test images over the test sites. The model produced image outputs in which pixel intensities indicated the proporton of forest cover per km². For comparison, supervised maximum likelihood classifications were also performed. The outputs were assessed against classified Landsat TM scenes, converted to proportions maps and coregistered to the AVHRR images. An empirical method was applied for determining the critical forest cover per km² needed for an AVHRR pixel to be classified as forest. The critical values exceeded 50 per cent, indicating a tendency for AVHRR classification to underestimate forest cover. This was confirmed by comparing estimates of total forest cover obtained from the AVHRR and TM classifications. In the case of Ghana, a more accurate estimate of forest cover was obtained from the AVHRR mixture model than from the classification. Both mixture model outputs were found to be well correlated with those from Landsat TM. Further work should test the robustness of the approach adopted here when applied to much larger areas.     

Trends in vegetation NDVI from 1 km AVHRR data over Canada for the period 1985–2006

Long‐term changes in the Normalized Difference Vegetation Index (NDVI) have been evaluated in several studies but results have not been conclusive due to differences in data processing as well as the length and time of the analysed period. In this research a newly developed 1 km Advanced Very High Resolution Radiometer (AVHRR) satellite data record for the period 1985–2006 was used to rigorously evaluate NDVI trends over Canada. Furthermore, climate and land cover change as potential causes of observed trends were evaluated in eight sample regions. The AVHRR record was generated using improved geolocation, cloud screening, correction for sun‐sensor viewing geometry, atmospheric correction, and compositing. Results from both AVHRR and Landsat revealed an increasing NDVI trend over northern regions where comparison was possible. Overall, 22% of the vegetated area in Canada was found to have a positive NDVI trend based on the Mann–Kendal test at the 95% confidence level. Of these, 40% were in northern ecozones. The mean absolute difference of NDVI measurements between AVHRR and Landsat data was <7%. When compared with results from other studies, similar trends were found over northern areas, while in southern regions the results were less consistent. Local assessment of potential causes of trends in each region revealed a stronger influence of climate in the north compared to the south. Southern regions with strong positive trends appeared to be most influenced by land cover change.     

Agreement assessment of NOAA/AVHRR NDVI with Landsat TM NDVI for mapping burned forested areas

Meteorological satellites are appropriate for operational applications related to early warning, monitoring and damage assessment of forest fires. Environmental or resources satellites, with better spatial resolution than meteorological satellites, enable the delineation of the affected areas with a higher degree of accuracy. In this study, the agreement of two datasets, coming from National Oceanic and Atmospheric Administration/Advanced Very High Resolution Radiometer (NOAA/AVHRR) and Landsat TM, for the assessment of the burned area, was investigated. The study area comprises a forested area, burned during the forest fire of 21-24 July 1995 in Penteli, Attiki, Greece. Based on a colour composite image of Landsat TM a reference map of the burned area was produced. The scatterplot of the multitemporal Normalized Difference Vegetation Index (NDVI) images, from both Landsat TM and NOAA/AVHRR sensors, was used to detect the spectral changes due to the removal of vegetation. The extracted burned area was compared to the digitized reference map. The synthesis of the maps was carried out using overlay techniques in a Geographic Information System (GIS). It is illustrated that the NOAA/AVHRR NDVI accuracy is comparable to that from Landsat TM data. As a result NOAA/AVHRR data can, operationally, be used for mapping the extent of the burned areas.     

Intercalibration of vegetation indices from different sensor systems

Spectroradiometric measurements were made over a range of crop canopy densities, soil backgrounds and foliage colour. The reflected spectral radiances were convoluted with the spectral response functions of a range of satellite instruments to simulate their responses. When Normalised Difference Vegetation Indices (NDVI) from the different instruments were compared, they varied by a few percent, but the values were strongly linearly related, allowing vegetation indices from one instrument to be intercalibrated against another. A table of conversion coefficents is presented for AVHRR, ATSR-2, Landsat MSS, TM and ETM+, SPOT-2 and SPOT-4 HRV, IRS, IKONOS, SEAWIFS, MISR, MODIS, POLDER, Quickbird and MERIS (see Appendix A for glossary of acronyms). The same set of coefficients was found to apply, within the margin of error of the analysis, for the Soil Adjusted Vegetation Index SAVI. The relationships for SPOT vs. TM and for ATSR-2 vs. AVHRR were directly validated by comparison of atmospherically corrected image data. The results indicate that vegetation indices can be interconverted to a precision of 1-2%. This result offers improved opportunities for monitoring crops through the growing season and the prospects of better continuity of long-term monitoring of vegetation responses to environmental change.     

Landscape-scale characterization of cropland in China using Vegetation and Landsat TM images

In this landscape-scale study we explored the potential for multitemporal 10-day composite data from the Vegetation sensor to characterize land cover types, in combination with Landsat TM image and agricultural census data. The study area (175 km by 165 km) is located in eastern Jiangsu Province, China. The Normalized Difference Vegetation Index (NDVI ) and the Normalized Difference Water Index (NDWI ) were calculated for seven 10-day composite (VGT-S10) data from 11 March to 20 May 1999. Multi-temporal NDVI and NDWI were visually examined and used for unsupervised classification. The resultant VGT classification map at 1 km resolution was compared to the TM classification map derived from unsupervised classification of a Landsat 5 TM image acquired on 26 April 1996 at 30 m resolution to quantify percent fraction of cropland within a 1 km VGT pixel; resulting in a mean of 60% for pixels classified as cropland, and 47% for pixels classified as cropland/natural vegetation mosaic. The estimates of cropland area from VGT data and TM image were also aggregated to county-level, using an administrative county map, and then compared to the 1995 county-level agricultural census data. This landscape-scale analysis incorporated image classification (e.g. coarse-resolution VGT data, fineresolution TM data), statistical census data (e.g. county-level agricultural census data) and a geographical information system (e.g. an administrative county map), and demonstrated the potential of multi-temporal VGT data for mapping of croplands across various spatial scales from landscape to region. This analysis also illustrated some of the limitations of per-pixel classification at the 1 km resolution for a heterogeneous landscape.     

Decomposition of vegetation cover into woody and herbaceous components using AVHRR NDVI time series

A method is developed to separate Normalised Difference Vegetation Index (NDVI) time series data into contributions from woody (perennial) and herbaceous (annual) vegetation, and thereby to infer their separate leaf area indices and cover fractions. The method is formally consistent with fundamental linearity requirements for such a decomposition, and is capable of rejecting contaminated NDVI data. In this study, estimates of annual averaged woody cover and monthly averaged herbaceous cover over Australia are determined using Pathfinder AVHRR Land series (PAL) Global Area Coverage (GAC) Advanced Very High Resolution Radiometer (AVHRR) NDVI data from 1981 to 1994, together with ground-based measurements of leaf area index (LAI) and foliage projective cover (FPC).     

Multi-temporal analysis of deforestation in Rondônia state in Brazil using Landsat MSS,TM, ETM+ and NOAA AVHRR imagery and its relationship to changes in the local hydrological environment

Deforestation in Rondônia state in the south-western part of the Brazilian Legal Amazon was analysed using Landsat Multi-Spectral Scanner (MSS), Thematic Mapper (TM), Enhanced Thematic Mapper Plus (ETM+), National Oceanic & Atmospheric Administration (NOAA) Advanced Very High Resolution Radiometer (AVHRR) and hydrological data. The Landsat sensor data coverage was supplemented with Pathfinder AVHRR Land (PAL) Normalized Difference Vegetation Index (NDVI) datasets. The results from the Landsat-based analysis show that more than 30% of the natural vegetation in the study area was subject to deforestation between 1973 and 1999, a finding reinforced by analysis of the PAL NDVI data. In addition, it was established that trends in the PAL NDVI datasets were coincident with the pattern of deforestation. Apart from imagery analysis, time variations in the hydrological data between 1982 and 1988 were used to estimate the evapotranspiration. A decreasing trend was identified in the rate of evapotranspiration, suggesting that deforestation has a significant impact on the local hydrological cycle.     

Multisensor comparison of NDVI for a semi‐arid environment in Spain

The joint use of multiresolution sensors from different satellites offers many opportunities to describe vegetation and its dynamics. This paper introduces the concept of a virtual constellation (defined as an ensemble of all Earth Observation satellites in orbit that satisfy common requirements) for agricultural applications and contributes to providing the necessary inter‐sensor calibration methodology for spectral reflectances and NDVI. For this purpose, we performed an observational study, comparing reflectances and the Normalized Difference Vegetation Index (NDVI), from near‐synchronous image pairs of Landsat 7 Enhanced Thematic Mapper Plus (ETM+), as the reference sensor and Landsat 5 Thematic Mapper (TM), IRS 1C/D LISS‐III (LISS), Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), QuickBird, and NOAA Advanced Very High‐resolution Radiometer (AVHRR). Linear relationships were found for the intercalibration of reflectances and NDVI from one sensor to another, for all sensors, provided that some spatial aggregation was performed. The main source of data dispersion in our linear cross‐sensor translation equations is the geolocation uncertainty inherent in the process of geometric correction. Consequently, spatial aggregation always needs to be performed if (different or the same) sensors are to be used to derive time‐series of biogeophysical parameters over heterogeneous areas. The homogenous zone approach developed here is recommended as an excellent tool for deriving robust new cross‐sensor relationships, provided that the selected homogeneous crops cover the full NDVI range. The linear cross‐sensor relationships derived from one image pair were shown to be valid for the whole season and for all areas with similar vegetation and climate.     

Analysis of time-series MODIS 250 m vegetation index data for crop classification in the U.S. Central Great Plains

The global environmental change research community requires improved and up-to-date land use/land cover (LULC) datasets at regional to global scales to support a variety of science and policy applications. Considerable strides have been made to improve large-area LULC datasets, but little emphasis has been placed on thematically detailed crop mapping, despite the considerable influence of management activities in the cropland sector on various environmental processes and the economy. Time-series MODIS 250 m Vegetation Index (VI) datasets hold considerable promise for large-area crop mapping in an agriculturally intensive region such as the U.S. Central Great Plains, given their global coverage, intermediate spatial resolution, high temporal resolution (16-day composite period), and cost-free status. However, the specific spectral-temporal information contained in these data has yet to be thoroughly explored and their applicability for large-area crop-related LULC classification is relatively unknown. The objective of this research was to investigate the general applicability of the time-series MODIS 250 m Enhanced Vegetation Index (EVI) and Normalized Difference Vegetation Index (NDVI) datasets for crop-related LULC classification in this region. A combination of graphical and statistical analyses were performed on a 12-month time-series of MODIS EVI and NDVI data from more than 2000 cropped field sites across the U.S. state of Kansas. Both MODIS VI datasets were found to have sufficient spatial, spectral, and temporal resolutions to detect unique multi-temporal signatures for each of the region's major crop types (alfalfa, corn, sorghum, soybeans, and winter wheat) and management practices (double crop, fallow, and irrigation). Each crop's multi-temporal VI signature was consistent with its general phenological characteristics and most crop classes were spectrally separable at some point during the growing season. Regional intra-class VI signature variations were found for some crops across Kansas that reflected the state's climate and planting time differences. The multi-temporal EVI and NDVI data tracked similar seasonal responses for all crops and were highly correlated across the growing season. However, differences between EVI and NDVI responses were most pronounced during the senescence phase of the growing season.     

Evaluating different NDVI composite techniques using NOAA-14 AVHRR data

Normalized Difference Vegetation Index (NDVI) data derived from Advanced Very High Resolution Radiometer (AVHRR) data are influenced by cloud contamination, which is common in individual AVHRR scenes. Maximum value compositing (MVC) of NDVI data has been employed to minimize cloud contamination. Two types of weekly NDVI composites were built for crop seasons in summer: one from all available AVHRR data (named the traditional NDVI composite) and the other from solely cloud-free AVHRR data (named the conditional NDVI composite). The MVC method was applied to both composites. The main objective of this study was to compare the two types of NDVI composites using Texas data. The NDVI seasonal profiles produced from the conditional NDVI composites agreed with the field measured leaf area index (LAI) data, reaching maximum values at similar times. However, the traditional NDVI composites showed irregular patterns, primarily due to cloud contamination. These study results suggest that cloud detection for individual AVHRR scenes should be strongly recommended before producing weekly NDVI composites. Appropriate AVHRR data pre-processing is important for composite products to be used for short-term vegetation condition and biomass studies, where the traditional NDVI composite data do not eliminate cloud-contaminated pixels. In addition, this study showed that atmosphere composition affected near-infrared reflectance more than visible reflectance. The near-infrared reflectance was increasingly adjusted through atmospheric correction.     

Land‐cover classification in the Brazilian Amazon with the integration of Landsat ETM+ and Radarsat data

Land‐cover classification with remotely sensed data in moist tropical regions is a challenge due to the complex biophysical conditions. This paper explores techniques to improve land‐cover classification accuracy through a comparative analysis of different combinations of spectral signatures and textures from Landsat Enhanced Thematic Mapper Plus (ETM+) and Radarsat data. A wavelet‐merging technique was used to integrate Landsat ETM+ multispectral and panchromatic data or Radarsat data. Grey‐level co‐occurrence matrix (GLCM) textures based on Landsat ETM+ panchromatic or Radarsat data and different sizes of moving windows were examined. A maximum‐likelihood classifier was used to implement image classification for different combinations. This research indicates the important role of textures in improving land‐cover classification accuracies in Amazonian environments. The incorporation of data fusion and textures increases classification accuracy by approximately 5.8–6.9% compared to Landsat ETM+ data, but data fusion of Landsat ETM+ multispectral and panchromatic data or Radarsat data cannot effectively improve land‐cover classification accuracies.     

Subpixel forest cover in central Africa from multisensor, multitemporal data

Seven Landsat Multispectral Scanner (MSS) scenes in central Africa were coregistered with 8 km resolution data from the 1987 AVHRR Pathfinder Land data set. Percent forest cover in each 8 km grid cell was derived from the classified MSS scenes. Linear relationships between percent forest cover and 30 multitemporal metrics derived from all AVHRR optical and thermal channels were determined. Correlations were strongest for the mean annual normalized difference vegetation index (NDVI) and mean annual brightness temperature (AVHRR Channel 3) and weakest for those metrics, besides NDVI, based on near-infrared reflectances (AVHRR Channel 2). The relationships were used to estimate percent forest cover in various locations in the study area using multiple linear regression and regression trees. Overall, the multiple linear regression provided more accurate results. Predicted percent forest cover estimates were within 20% of the “actual” percent forest cover (derived from the MSS data) for approximately 90% of the grid cells. The RMS error for the prediction was 12% forest cover. RMS errors above 18% forest cover were obtained when using AVHRR data from a single month to derive predictive relationships. The results demonstrate that multitemporal data reflecting vegetation phenology can be used to estimate subpixel forest cover at coarse spatial resolutions.     

基于ISODATA法的冬小麦产量分级监测预报

以江苏省姜堰市为例,进行了基于TM卫星遥感技术和小麦估产模型的冬小麦产量监测研究。在利用GPS实地采样调查和建立解译标志的基础上,通过影像校正、采用优化的ISODATA分类方法,结合人机交互式判读解译等操作,将样点的作物信息数据贯穿到整个校验分类过程中,信息解译精度在90%以上。利用分类提取的冬小麦数据,反演叶面积指数、生物量信息等,结合冬小麦估产模型,计算单点产量信息,经过线性转换,对整个区域的冬小麦产量进行监测预报,并制作了冬小麦产量分级专题图。     

Drought assessment and monitoring through remote sensing and GIS in western tracts of Tamil Nadu,India

Drought is an insidious hazard of nature and is considered to be the most complex but least understood of all natural hazards. Large historical datasets are required to study drought and these involve complex interrelationships between climatological and meteorological data. Rainfall is an important meteorological parameter; the amount and distribution influence the type of vegetation in a region. To analyse the changes in vegetation cover due to variation in rainfall and identify the land-use areas facing drought risk, rainfall data from 1981 to 2003 were categorized into excess, normal, deficit and drought years. The Advanced Very High Resolution Radiometer (AVHRR) sensor's composite dataset was used for analysing the temporal and interannual behaviour of surface vegetation. The various land-use classes – crop land (annual, perennial crops), scrub land, barren land, forest land, degraded pasture and grassland – were identified using satellite data for excess, normal, deficit and drought years. Normalized Difference Vegetation Indices (NDVIs) were derived from satellite data for each land-use class and the highest NDVI mean values were 0.515, 0.436 and 0.385 for the tapioca crop in excess, normal and deficit years, respectively, whereas in the drought year, the groundnut crop (0.267) showed the maximum. Grassland recorded the lowest value of NDVI in all years except for the excess year. Annual crops, such as groundnut (0.398), pulses (0.313), sorghum (0.120), tapioca (0.436) and horse gram (0.259), registered comparatively higher NDVI values than the perennial crops for the normal year. The Vegetation Condition Index (VCI) was used to estimate vegetation health and monitor drought. Among land-use classes, the maximum VCI value of 92.1% was observed in onions for the excess year, whereas groundnut witnessed the maximum values of 78.2, 64.5 and 55.2% for normal, deficit and drought years, respectively. Based on the VCI classification, all land-use classes fall into the optimal or normal vegetation category in excess and normal years, whereas in drought years most of the land-use classes fall into the drought category except for sorghum, groundnut, pulses and grasses. These crops (sorghum 39.7%, groundnut 55.2%, pulses 38.5% and grassland 38.6%) registered maximum VCI values, revealing that they were sustained under drought conditions. It is suggested that the existing crop pattern be modified in drought periods by selecting the suitable crops of sorghum, groundnut and pulses and avoiding the cultivation of onion, rice and tapioca.     

Use of NDVI/AVHRR time-series profiles for soybean crop monitoring in Brazil

In Brazil there is a need for less subjective, more efficient and less expensive methodologies for crop yield forecast. Owing to the continental dimensions of the country, orbital images have been used to estimate the productive potential of crops. In this study, NDVI (Normalized Difference Vegetation Index) time-series, derived from AVHRR/NOAA (Advanced Very High Resolution Radiometer/National Oceanic and Atmospheric Administration) imagery were used for the soybean crop monitoring in a large production region in Brazil in the 2002/2003 and 2003/2004 cropping seasons. NDVI temporal profiles describing the biomass condition of crops throughout the phenological stages were generated in 18 municipalities. Quantitative parameters were measured from the temporal profiles, based on the full time or partial phenological cycle. Linear regressions between the quantitative parameters and the municipal average yields in both seasons have shown that the most significant correlations occurred when the full time period was considered. When considering periods prior to harvest, the correlations showed a tendency to decline. The NDVI monitoring during these two cropping seasons, which presented different weather conditions, could explain a major part of the soybean yield variability at the municipal level. Results showed the potential of the NDVI time-series analysis in generating parameters to be employed by agrometeorological–spectral models for soybean yield estimations. The automatic system for temporal profiles generation developed in this study sped up the analysis and can be used for further studies at a regional scale.     

Scene selection and the use of NASA's global orthorectified Landsat dataset for land cover and land use change monitoring

This study examines the utility of NASA's circa 1990 and circa 2000 global orthorectified Landsat dataset for land cover and land use change mapping and monitoring across Africa. This is achieved by comparing the temporal and spatial variation of NDVI, measured independently by the NOAA‐AVHRR at the time of Landsat scene acquisition, against the seasonal mean for each Landsat scene extent. Decadal sequences of drift‐corrected NOAA‐AVHRR imagery were used to calculate NDVI means and standard deviations for the periods covered by the scenes composing the c.1990 and c.2000 Landsat datasets. The specific NOAA‐AVHRR NDVI values at the acquisition date of each individual Landsat scene were also calculated and the differences, both from the mean and scaled by standard deviation, were mapped for the Landsat scene footprints in the c.1990 and c.2000 datasets. The resulting maps show the temporal position of each Landsat scene within the seasonal NDVI cycle, and provide a valuable guide to assist in quantifying uncertainty and interpreting land cover and land use changes inferred from these Landsat data.     

基于AVHRR和TM数据的时间序列较高分辨率NDVI数据集重构方法

由于技术条件的限制,一个传感器很难同时具有高空间分辨率和高时间分辨率。然而,在高分辨率尺度上监测地表景观季节性变化的能力是全球的迫切需要,融合周期短、覆盖范围广与分辨率高、周期长的遥感数据是一种较好的方法。基于AVHRR时间分辨率高和TM空间分辨率高及其数据积累时间长的特点,选择若尔盖高原为研究区域,在改进ESTARFM方法的基础上,对TM NDVI和AVHRR NDVI进行融合,构建高时空分辨率的NDVI数据集。研究结果表明:该方法能有机结合AVHRR NDVI的时间变化信息与TM NDVI的空间差异信息,有效实现高时空分辨率NDVI数据集的重构,3景预测高分辨率NDVI与MODIS NDVI产品相关系数分别达到了0.89、0.91和0.85。该方法能够在时间上保留高时间分辨率数据的时间变化信息,同时在空间上反映高空间分辨率数据的空间差异信息,从而为有效构建相对高分辨率时间序列NDVI数据集提供了可能的方法。