首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 187 毫秒
1.
Remote sensing based biomass estimates in Arctic areas are usually produced using coarse spatial resolution satellite imagery, which is incapable of capturing the fragmented nature of tundra vegetation communities. We mapped aboveground biomass using field sampling and very high spatial resolution (VHSR) satellite images (QuickBird, WorldView-2 and WorldView-3) in four different Arctic tundra or peatland sites with low vegetation located in Russia, Canada, and Finland. We compared site-specific and cross-site empirical regressions. First, we classified species into plant functional types and estimated biomass using easy, non-destructive field measurements (cover, height). Second, we used the cover/height-based biomass as the response variable and used combinations of single bands and vegetation indices in predicting total biomass. We found that plant functional type biomass could be predicted reasonably well in most cases using cover and height as the explanatory variables (adjusted R2 0.21–0.92), and there was considerable variation in the model fit when the total biomass was predicted with satellite spectra (adjusted R2 0.33–0.75). There were dissimilarities between cross-site and site-specific regression estimates in satellite spectra based regressions suggesting that the same regression should be used only in areas with similar kinds of vegetation. We discuss the considerable variation in biomass and plant functional type composition within and between different Arctic landscapes and how well this variation can be reproduced using VHSR satellite images. Overall, the usage of VHSR images creates new possibilities but to utilize them to full potential requires similarly more detailed in-situ data related to biomass inventories and other ecosystem change studies and modelling.  相似文献   

2.
Abstract

The standing crop of herbaceous biomass produced during the 2-4?month summer rainy season by the annual grasses in the Sahel zone provides an indication of resource availability for livestock for the following 9-month dry season. Combined use of NOAA advanced very high resolution radiometer (AVHRR) local area coverage (LAC) satellite data and biomass data, obtained through vegetation sampling of 25-100 km2 areas, allowed the development of a method for biomass assessment in Niger. Vegetation sampling involved both visual estimates and clipped plots (double sampling). The relationship between time-integrated normalized difference vegetation index (NDVI) statistics derived from NOAA AVHRR LAC data (dependent variable) and total herbaceous biomass (independent variable) was obtained through regression analysis. An inverse prediction was used to estimate biomass from the satellite data. Biomass maps and statistics of the grasslands were produced for the end of each rainy season: 1986, 1987 and 1988. This information is being used for planning purposes by the pastoral resource managers of the Government of Niger.  相似文献   

3.
ABSTRACT

Satellite remote sensing has greatly facilitated the assessment of aboveground biomass in rangelands. Soil-adjusted vegetation indices have been developed to provide better predictions of aboveground biomass, especially for dryland regions. Semi-arid rangelands often complicate a remote sensing based assessment of aboveground biomass due to bright reflecting soils combined with sparse vegetation cover. We aim at evaluating whether soil-adjusted vegetation indices perform better than standard, i.e. unadjusted, vegetation indices in predicting dry aboveground biomass of a saline and semi-arid rangeland in NE-Iran. 672 biomass plots of 2 × 2 m were gathered and aggregated into 13 sites. Generalized Linear Regression Models (GLM) were compared for six different vegetation indices, three standard and three soil-adjusted vegetation indices. Vegetation indices were calculated from the MODIS MCD43A4 product. Model comparison was done using Akaike Information Criterion (AICc), Akaike weights and pseudo R2. Model fits for dry biomass showed that transformed NDVI and NDVI fitted best with R2 = 0.47 and R2 = 0.33, respectively. The optimized soil-adjusted vegetation index (OSAVI) behaved similar to NDVI but less precise. The soil-adjusted vegetation index (SAVI), the modified soil-adjusted vegetation index (MSAVI2) and the enhanced vegetation index (EVI) performed worse than a null model. Hence, soil-adjusted indices based on the soil-line concept performed worse than a simple square root transformation of the NDVI. However, more studies that compare MODIS based vegetation indices for rangeland biomass estimation are required to support our findings. We suggest applying a similar model comparison approach as performed in this study instead of relying on single vegetation indices in order to find optimal relationships with aboveground biomass estimation in rangelands.  相似文献   

4.
Pecan orchards are the largest agricultural water consumer in the lower part of the Mesilla Valley, NM, USA. Knowledge of fractional canopy (FC) cover allows better crop water use assessment and orchard management. FC can be estimated from vegetation indices (VIs), such as the normalized difference vegetation index (NDVI), the soil-adjusted vegetation index (SAVI), the simple ratio (SR), and the triangular vegetation index (TVI), using satellite imagery. The main objective of this research is to develop an approach to determine the FC from a simultaneous imagery campaign consisting of aerial imagery, orchard floor photographs, and satellite images. All the required data were collected based on satellite overpass times at three different times during the initial part of the growing season to enhance the quality of data and reduce errors. The data were processed using the software package Environment for Visualizing Images (ENVI® 4.6.1; ITT Research Systems Inc.). The orchard floor digital photographs were used as a ground truth data set that gave a good correlation to the aerial photography. The aerial images were then used to determine the relationship between the FC and the VIs using these ‘corrected FCs’. The results showed significant correlation between NDVI and FC (R 2 = 0.80; p < 0.0001). Likewise, the calculated SR not only showed good correlation to the FCs but also verified the calculated NDVI. The results indicated that the methodology of this research can be applied to other tree crops as an aid in estimating the FC.  相似文献   

5.
Much effort has been made in recent years to improve the spectral and spatial resolution of satellite sensors to develop improved vegetation indices reflecting surface conditions. In this study satellite vegetation indices from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Advanced Very High Resolution Radiometer (AVHRR) are evaluated against two years of in situ measurements of vegetation indices in Senegal. The in situ measurements are obtained using four masts equipped with self‐registrating multispectral radiometers designed for the same wavelengths as the satellite sensor channels. In situ measurements of the MODIS Normalized Difference Vegetation Index (NDVI) and AVHRR NDVI are equally sensitive to vegetation; however, the MODIS NDVI is consistently higher than the AVHRR NDVI. The MODIS Enhanced Vegetation Index (EVI) proved more sensitive to dense vegetation than both AVHRR NDVI and MODIS NDVI. EVI and NDVI based on the MODIS 16‐day constrained view angle maximum value composite (CV‐MVC) product captured the seasonal dynamics of the field observations satisfactorily but a standard 16‐day MVC product estimated from the daily MODIS surface reflectance data without view angle constraints yielded higher correlations between the satellite indices and field measurements (R 2 values ranging from 0.74 to 0.98). The standard MVC regressions furthermore approach a 1?:?1 line with in situ measured values compared to the CV‐MVC regressions. The 16‐day MVC AVHRR data did not satisfactorily reflect the variation in the in situ data. Seasonal variation in the in situ measurements is captured reasonably with R 2 values of 0.75 in 2001 and 0.64 in 2002, but the dynamic range of the AVHRR satellite data is very low—about a third to a half of the values from in situ measurements. Consequently the in situ vegetation indices were emulated much better by the MODIS indices than by the AVHRR NDVI.  相似文献   

6.
Woody biomass production is a critical indicator in evaluation of land use management and the dynamics of the global carbon cycle (sequestration/emission) in terrestrial ecosystems. The objective of the present study was to develop, through a case study in Sudan, an operational multiscale remote-sensing-based methodology for large-scale estimation of woody biomass in tropical savannahs. Woody biomass estimation models obtained by different authors from destructive field measurements in different tropical savannah ecosystems were expressed as functions of tree canopy cover (CC). The field-measured CC data were used for developing regression equations with atmospherically corrected and reflectance-based vegetation indices derived from Landsat ETM+ (Enhanced Thematic Mapper Plus) imagery. Among a set of vegetation indices, the normalized difference vegetation index (NDVI) provided the best correlation with CC (R2 = 0.91) and was hence selected for woodland woody biomass estimation. After validation of the CC-NDVI model and its applicability to Moderate Resolution Imaging Spectroradiometer (MODIS) data, time-series MODIS NDVI data (MOD13Q1) were used to partition the woody component from the herbaceous component for sparse woodlands, woodlands and forests defined by the Food and Agriculture Organization (FAO) of the United Nations Land Cover Map. Following the weighting of the estimation models based on the dominant woody species in each vegetation community, NDVI-based woody biomass models were applied according to their weighted ratios to the decomposed summer and autumn woody NDVI images in all vegetation communities in the whole of Sudan taking the year 2007, for example. The results were found to be in good agreement with those from other authors obtained by either field measurements or other remote sensing methods using MODIS and lidar data. It is concluded that the proposed approach is operational and can be applied for a reliable large-scale assessment of woody biomass at a ground resolution of 250 m in tropical savannah woodlands of any month or season.  相似文献   

7.
Vegetation indices (VIs) such as the Normalized Difference Vegetation Index (NDVI) are widely used for assessing vegetation cover and condition. One of the NDVI's significant disadvantages is its sensitivity to aerosols in the atmosphere, hence several atmospherically resistant VIs were formulated using the difference in the radiance between the blue and the red spectral bands. The state‐of‐the‐art atmospherically resistant VI, which is a standard Moderate Resolution Imaging Spectroradiometer (MODIS) product, together with the NDVI, is the Enhanced Vegetation Index (EVI). A different approach introduced the Aerosol‐free Vegetation Index (AFRI) that is based on the correlation between the shortwave infrared (SWIR) and the visible red bands. The AFRI main advantage is in penetrating an opaque atmosphere influenced by biomass burning smoke, without the need for explicit correction for the aerosol effect. The objective of this research was to compare the performance of these three VIs under smoke conditions. The AFRI was applied to the 2.1 µm SWIR channel of the MODIS sensor onboard the Earth Observing System (EOS) Terra and Aqua satellites in order to assess its functionality on these imaging platforms. The AFRI performance was compared with those of NDVI and EVI. All VIs were calculated on images with and without present smoke, using the surface‐reflectance MODIS product, for three case studies of fires in Arizona, California, and Zambia. The MODIS Fire Product was embedded on the images in order to identify the exact location of the active fires. Although good correlations were observed between all VIs in the absence of smoke (in the Arizona case R 2 = 0.86, 0.77, 0.88 for the NDVI–EVI, AFRI–EVI, and AFRI–NDVI, respectively) under smoke conditions a high correlation was maintained between the NDVI and the EVI, while low correlations were found for the AFRI–EVI and AFRI–NDVI (0.21 and 0.16, for the Arizona case, respectively). A time series of MODIS images recorded over Zambia during the summer of 2000 was tested and showed high NDVI fluctuations during the study period due to oscillations in aerosol optical thickness values despite application of aerosol corrections on the images. In contrast, the AFRI showed smoother variations and managed to better assess the vegetation condition. It is concluded that, beneath the biomass burning smoke, the AFRI is more effective than the EVI in observing the vegetation conditions.  相似文献   

8.
The fraction of photosynthetically active radiation (FPAR) absorbed by vegetation – a key parameter in crop biomass and yields as well as net primary productivity models – is critical to guiding crop management activities. However, accurate and reliable estimation of FPAR is often hindered by a paucity of good field-based spectral data, especially for corn crops. Here, we investigate the relationships between the FPAR of corn (Zea mays L.) canopies and vegetation indices (VIs) derived from concurrent in situ hyperspectral measurements in order to develop accurate FPAR estimates. FPAR is most strongly (positively) correlated to the green normalized difference vegetation index (GNDVI) and the scaled normalized difference vegetation index (NDVI*). Both GNDVI and NDVI* increase with FPAR, but GNDVI values stagnate as FPAR values increase beyond 0.75, as previously reported according to the saturation of VIs – such as NDVI – in high biomass areas, which is a major limitation of FPAR-VI models. However, NDVI* shows a declining trend when FPAR values are greater than 0.75. This peculiar VI–FPAR relationship is used to create a piecewise FPAR regression model – the regressor variable is GNDVI for FPAR values less than 0.75, and NDVI* for FPAR values greater than 0.75. Our analysis of model performance shows that the estimation accuracy is higher, by as much as 14%, compared with FPAR prediction models using a single VI. In conclusion, this study highlights the feasibility of utilizing VIs (GNDVI and NDVI*) derived from ground-based spectral data to estimate corn canopy FPAR, using an FPAR estimation model that overcomes limitations imposed by VI saturation at high FPAR values (i.e. in dense vegetation).  相似文献   

9.
Remotely sensed vegetation indices such as NDVI, computed using the red and near infrared bands have been used to estimate pasture biomass. These indices are of limited value since they saturate in dense vegetation. In this study, we evaluated the potential of narrow band vegetation indices for characterizing the biomass of Cenchrus ciliaris grass measured at high canopy density. Three indices were tested: Modified Normalized Difference Vegetation Index (MNDVI), Simple Ratio (SR) and Transformed Vegetation Index (TVI) involving all possible two band combinations between 350?nm and 2500?nm. In addition, we evaluated the potential of the red edge position in estimating biomass at full canopy cover. Results indicated that the standard NDVI involving a strong chlorophyll absorption band in the red region and a near infrared band performed poorly in estimating biomass (R 2=0.26). The MNDVIs involving a combination of narrow bands in the shorter wavelengths of the red edge (700–750?nm) and longer wavelengths of the red edge (750–780?nm), yielded higher correlations with biomass (mean R 2=0.77 for the highest 20 narrow band NDVIs). When the three vegetation indices were compared, SR yielded the highest correlation coefficients with biomass as compared to narrow band NDVI and TVI (average R 2=0.80, 0.77 and 0.77 for the first 20 ranked SR, NDVI and TVI respectively). The red edge position yielded comparable results to the narrow band vegetation indices involving the red edge bands. These results indicate that at high canopy density, pasture biomass may be more accurately estimated by vegetation indices based on wavelengths located in the red edge than the standard NDVI.  相似文献   

10.
Estimates of biomass production are important in a wildlife reserve such as Etosha National Park, Namibia, for assessment of fire risk and subsequent selection of sites for controlled burning. We present methodology for using locally acquired NOAA-AVHRR images to make estimates of biomass in near-to-realtime. To this end, techniques for rapid measurement of the biomass of herbaceous and woody vegetation were developed using a rising disc pasture meter and individual plant dimensions. A field sampling methodology is presented to make biomass estimates which were compatible with the scale of AVHRR spatial resolution and sufficiently close to the time of satellite overpasses to enable correlation with the NDVI from single images. Initial results show high correlations of biomass with NDVI for individual vegetation cover classes, which appear to be temporally stable. There seem to be different regression equations for the different savanna vegetation types although more field observations are needed to confirm this. The results were exploited to illustrate the potential application of this work for fire management. The combination of rapid field methods and real time image acquisition developed in this work provides a sound basis for biomass monitoring at local level.  相似文献   

11.
12.
ABSTRACT

Land-cover mapping in complex farming area is a difficult task because of the complex pattern of vegetation and rugged mountains with fast-flowing rivers, and it requires a method for accurate classification of complex land cover. Random Forest classification (RFC) has the advantages of high classification accuracy and the ability to measure variable importance in land-cover mapping. This study evaluates the addition of both normalized difference vegetation index (NDVI) time-series and the Grey Level Co-occurrence Matrix (GLCM) textural variables using the RFC for land-cover mapping in a complex farming region. On this basis, the best classification model is selected to extract the land-cover classification information in Central Shandong. To explore which input variables yield the best accuracy for land-cover classification in complex farming areas, we evaluate the importance of Random Forest variables. The results show that adding not only multi-temporal imagery and topographic variables but also GLCM textural variables and NDVI time-series variables achieved the highest overall accuracy of 89% and kappa coefficient (κ) of 0.81. The assessment of the importance of a Random Forest classifier indicates that the key input variables include the summer NDVI followed by the summer near-infrared band and the elevation, along with the GLCM-mean, GLCM-contrast.  相似文献   

13.
Various aspects of global environmental change affect plant photosynthesis, the primary carbon input in ecosystems. Thus, accurate methods of measuring plant photosynthesis are important. Remotely sensed spectral indices can monitor in detail the green biomass of ecosystems, which provides a measure of potential photosynthetic capacity. In evergreen vegetation types, however, such as Mediterranean forests, the amount of green biomass changes little during the growing season and, therefore, changes in green biomass are not responsible for changes in photosynthetic rates in those forests. This study examined the net photosynthetic rates and the diametric increment of stems in a Mediterranean forest dominated by Quercus ilex using three spectral indices (normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), and photochemical reflectance index (PRI)) derived from Moderate Resolution Imaging Spectroradiometer (MODIS) sensors. Average annual EVI accounted for 83% of the variability of the diametric increment of Q. ilex stems over a 10 year period. NDVI was marginally correlated with the diametric increment of stems. This study was the first to identify a significant correlation between net photosynthetic rates and radiation use efficiency at the leaf level using PRI derived from satellite data analysed at the ecosystem level. These results suggest that each spectral index provided different and complementary information about ecosystem carbon uptake in a Mediterranean Q. ilex forest.  相似文献   

14.
NOAA Advanced Very High Resolution Radiometer satellite data are applied to regional vegetation monitoring in East Africa. Normalized Difference Vegetation Index (NDVI) data for a one-year period from May 1983 are used to examine the phenology of a range of vegetation types. The integrated NDVI data for the same period are compared with an ecoclimatic zone map of the region and show marked similarities. Particular emphasis is placed on quantifying the phenology of the Acacia Commiphora bushlands. Considerable variation was found in the phenology of the bushlands as determined by the satellite NDVI, and is explained through the high spatial variability in the distribution of rainfall and the resulting green-up of the vegetation. The relationship between rainfall and NDVI is further examined for selected meteorological stations existing within the bushland. A preliminary estimate is made of the length of growing season using an NDVI thresholding technique  相似文献   

15.
A knowledge of the amount of pasture biomass available in farm paddocks is crucial for improving utilization and productivity in the Australian grazing industry. A method to quantitatively map the biomass of annual pastures under grazing has been developed using the Normalized Difference Vegetation Index (NDVI) derived from high-resolution satellite imagery. Relationships between field-measured pasture biomass and the NDVI were examined for different transects in paddocks under different grazing regimes across three geographically dispersed farm sites. A significant linear relationship (R 2 = 0.84) was observed when the NDVI was regressed against biomass. The slope of the relationship between the NDVI and biomass declined in a highly predictable (R 2 = 0.82) exponential form as the growing season progressed and this pattern was consistent across four separate seasons. This knowledge was used to formulate a reliable model to predict paddock average pasture biomass using the NDVI. The model estimates were validated against observed biomass in the range 500–4000 kilograms of dry matter per hectare (kg DM ha–1) with R 2 = 0.85 and a standard error of 315 (kg DM ha–1).  相似文献   

16.
Forest inventory data can be used along with remotely sensed data to estimate biomass and carbon stocks over large and inaccessible forested areas. In this study, the relationship between satellite-derived multispectral data and forest variables from intervened and non-intervened Nothofagus pumilio forest stands located in the Magellan region of Chile was examined, in order to quantify the over bark volume (OBV) and aboveground tree biomass (AGTB). Four vegetation parameters – the green normalised difference vegetation index (GNDVI), normalised difference vegetation index (NDVI), simple ratio (SR) and vegetation cover fraction (VCF) – were retrieved from an Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) image of the study area. The results indicate that only the VCF presents significant differences among intervened and non-intervened stands. The best OBV and AGTB models (R2 = 0.58) were found using the SR index and the VCF as predictors. This result could be transferred to estimate biomass and volume in other Nothofagus pumilio forests with similar conditions. Moreover, it can be used to assess temporal carbon changes.  相似文献   

17.
The productivity of semi-arid rangelands on the Arabian Peninsula is spatially and temporally highly variable, and increasing grazing pressure as well as the likely effects of climatic change further threatens vegetation resources. Using the Al Jabal al Akhdar mountains in northern Oman as an example, our objectives were to analyse the availability and spatial distribution of aboveground net primary production (ANPP) and the extent and causes of vegetation changes during the last decades with a remote sensing approach. A combination of destructive and non-destructive biomass measurements by life-form specific allometric equations was used to identify the ANPP of the ground vegetation (< 50 cm) and the leaf and twig biomass of phanerophytes. The ANPP differed significantly among the life forms and the different plant communities, and the biomass of the sparsely vegetated ground was more than 50 times lower (mean = 0.22 t DM ha− 1) than the biomass of phanerophytes (mean = 12.3 t DM ha− 1). Among the different vegetation indices calculated NDVI proved to be the best predictor for rangeland biomass.Temporal trend analysis of Landsat satellite images from 1986 to 2009 was conducted using a pixel-based least square regression with the annual maximum Normalized Differenced Vegetation Index (NDVImax) as a dependent variable. Additionally, linear relationships of NDVImax and annual rainfall along the time series were calculated. The extent of human-induced changes was analysed using the residual trends method. A strongly significant negative biomass trend detected for 83% of the study area reflected a decrease in annual rainfall but even without clear evidence of deforestation of trees and shrubs, human-induced vegetation degradation due to settlement activities were also important.  相似文献   

18.
Abstract

In order to obtain a model equation for the calculation of percentage plant cover by multi-spectral radiances remotely-sensed by satellites, a regression procedure is used to connect space remote-sensing data to ground plant cover measurement. A traditional linear regression model using the normalized difference vegetation index (NDVI) is examined by remote-sensing data of the SPOT satellite and ground measurement of LCTA project for a test site at Hohenfels. Germany. A relaxation vegetation index (RVI) is proposed in a non-linear regression modelling to replace the NDVI in linear regression modelling to get a better calculation of percentage plant cover. The definition of the RVI is

where X i is raw remote-sensing data in channel i. Using the RVI, the correlation coefficient between calculated and observed percentage plant cover for a test scene in 1989 reaches 0·9 while for the NDVI it is only 0·7; the coefficient of multiple determination R 2 reaches 0·8 for the RVI while it is only 0·5 for the NDVI. Numerical testing shows that the ability of using the RVI to predict percentage plant cover by space remote-sensing data for the same scene or the scene in other years is much stronger than the NDVI.  相似文献   

19.

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.  相似文献   

20.
Do flowers affect biomass estimate accuracy from NDVI and EVI?   总被引:1,自引:0,他引:1  
The Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI) are vegetation indices widely used in remote sensing of above-ground biomass. Because both indexes are based on spectral features of plant canopy, NDVI and EVI may suffer reduced accuracy in estimating above-ground biomass when flower signals are mixed in the plant canopy. This paper addresses how flowers influence the estimation of above-ground biomass using NDVI and EVI for an alpine meadow with mixed yellow flowers of Halerpestes tricuspis (Ranunculaceae). Field spectral measurements were used in combination with simulated reflectance spectra with precisely controlled flower coverage by applying a linear spectral mixture model. Using the reflectance spectrum for the in-situ canopy with H. tricuspis flowers, we found no significant correlation between above-ground biomass and EVI (p?=?0.17) or between above-ground biomass and NDVI (p?=?0.78). However, both NDVI and EVI showed very good prediction of above-ground biomass with low root mean square errors (RMSE?=?43 g m?2 for NDVI and RMSE?=?43 g m?2 for EVI, both p < 0.01) when all the flowers were removed from the canopies. Simulation analysis based on the in-situ measurements further indicated that high variation in flower coverage among different quadrats could produce more noise in the relationship between above-ground biomass and NDVI, or EVI, which results in an evident decline in the accuracy of above-ground biomass estimation. Therefore, the study suggests that attention should be paid both to the flower fraction and the heterogeneity of flower distribution in the above-ground biomass estimation via NDVI and EVI.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号