首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 265 毫秒
1.
The Normalized Difference Vegetation Index (NDVI) is an important vegetation index, widely applied in research on global environmental and climatic change. However, noise induced by cloud contamination and atmospheric variability impedes the analysis and application of NDVI data. In this work, a simplified data assimilation method is proposed to reconstruct high-quality time-series MODIS NDVI data. We extracted 16-Day L3 Global 1 km SIN Grid NDVI data sets for western China from MODIS vegetation index (VI) products (MOD13A2) for the period 2003–2006. NDVI data in the first three years (2003–2005) were used to generate the background field of NDVI based on a simple three-point smoothing technique, which captures annual features of vegetation change. NDVI data for 2006 were used to test our method. For every time step, the quality assurance (QA) flags of the MODIS VI products were adopted to empirically determine the weight between the background field and NDVI observations. Ultimately, more reliable NDVI data can be produced. The results indicate that the newly developed method is robust and effective in reconstructing high-quality MODIS NDVI time-series.  相似文献   

2.
Vegetation typically elicits dynamics at the seasonal and annual level. Time-series of normalized difference vegetation index (NDVI) datasets, such as the pathfinder AVHRR land (PAL) NDVI dataset, have proven to be appropriate for the detection of long-term vegetation cover changes. It has been applied in modelling experiments for terrestrial ecosystems at the global, continental, and regional scales. But some PAL NDVI time series remain significant residual effects and noise levels. A simple method, the mean-value iteration filter (MVI), has been developed to reduce the noise and to enable the reconstruction of high quality NDVI time-series. A comparison between the newly developed method and other existing methods (the modified BISE algorithm and a fast Fourier transform algorithm) indicates that the newly developed method is an effective tool for reconstructing high-quality time series of PAL NDVI time series.  相似文献   

3.
This research explores the sensitivity of vegetation in China to El-Niño/Southern Oscillation (ENSO) events from 1982 to 2006. The ENSO events are defined by the Multivariate ENSO Index (MEI), and variation in vegetation cover is captured by the Advanced Very High Resolution Radiometer (AVHRR) Normalized Difference Vegetation Index (NDVI). Pearson’s χ2 test was used to identify the areas where the variation in vegetation was sensitive to El Niño and La Niña events. The difference in the sensitivity of various ecosystems was investigated using the Moderate Resolution Imaging Spectroradiometer (MODIS) land cover product in 2000. Composite NDVI graphs during El Niño, La Niña and non-ENSO years were also produced to investigate the ENSO relationship with the six vegetation ecosystems during El Niño, La Niña and normal phases. The results show that most of the ENSO-sensitive land in China is only affected by one of the two phases of ENSO events, and the area of El Niño-sensitive vegetation is much larger than that of La Niña-sensitive vegetation. North China and the Hengduan Mountains are the two cores of the El Niño-sensitive areas, while the La Niña-sensitive areas are mainly distributed in the central, northwest and northeast regions of China. The sensitivity of vegetation varies across ecosystems: grassland and shrubland had the largest share of El Niño-sensitive areas, and sparse vegetation and savanna were the most sensitive to La Niña events. Overall, the impacts of El Niño events on vegetation in China had regular seasonal variation, while the impacts of La Niña events had regular zonal distribution.  相似文献   

4.
Remote sensing applications have greatly enhanced ability to monitor and manage in the areas of forestry. Accurate measurements of regional and global scale vegetation dynamics (phenology) are required to improve models and understanding of inter-annual variability in terrestrial ecosystem carbon exchange and climate–biosphere interactions. Study of vegetation phenology is required for understanding of variability in ecosystem. In this paper, monitoring of vegetation dynamics using time series of satellite data is presented. Vegetation variability (vegetation rate) in different topoclimatic areas is investigated. Original software using IDL interactive language for processing of satellite long-term data series was developed. To investigate growth dynamics vegetation rate inferred from remote sensing was used. All estimations based on annual time series of Moderate Resolution Imaging Spectroradiometer (MODIS) imagery. Vegetation rate for Enhanced Vegetation Index (EVI) and Normalized Difference Vegetation Index (NDVI) was calculated using MODIS data. The time series covers spring seasons of each of 9 years, from 2000 to 2008. Comparison of EVI and NDVI derived growth rates has shown that NDVI derived rates reveal spatial structure better. Using long-term data of vegetation rates variance was estimated that helps to reveal areas with anomalous growth rate. Such estimation shows sensitivity degree of different areas to different topoclimatic conditions. Woods of heights depend on spatial topoclimatic variability unlike woods of lowlands. Principal components analysis shows vegetation with different rate conditions. Also it reveals vegetation of same type in areas with different conditions. It was demonstrated that using of methods for estimating the dynamic state of vegetation based on remote sensing data enables successful monitoring of vegetation phenology.  相似文献   

5.
The visible and near infrared channels, Ch1 and CH2 respectively, on the Advanced Very High Resolution Radiometer (AVHRR) provide daily information for monitoring changes in vegetation and crops. Data from these channels are used to create a normalized vegetation index (NVI) that is sensitive to changes in green leaf biomass and is represented mathematically by:
NVI = CH2 ? CH1CH2 + CH1
Operational products generated at NOAA include full-scale 1-km resolution images of the NVI covering areas viewed in a single swath of the polar-orbiting NOAA satellite. Global scale NVI images are also produced by compositing over a seven-day period, saving the maximum NVI created daily for each local array (resolution of 15 km at the equator to 30 km at the poles). Such seven-day mapping reduces the effect of cloud contamination. The global vegetation indices are used by foreign and U.S. government agencies for operational and experimental purposes such as assessment of crop conditions, monitoring potential desert locust breeding grounds, forest fire danger models, and monitoring range lands for forage availability. Examples include changes in the NVI in the Lake Chad vicinity, 1981–1982 and 1984; western United States NVI; and seasonal variations of the NVI in the Sahel using the global operational data base, 1982–1983.  相似文献   

6.
The Dongting Lake Basin is an important hydrological regulation and flood storage area in the Yangtze River Basin, which plays an important role in maintaining regional ecological security. The watershed vegetation and its carbon sequestration capacity have changed dramatically due to climate change and human activities in the last two decades. In this paper, the monthly and annual vegetation net primary productivity (NPP) of the Dongting Lake basin during 2000 to 2020 was firstly estimated using the improved Carnegie-Ames-Stanford Approach (CASA) model. Then the multi-year NPP change trend and its significance were analyzed based on Theil-Sen median and Mann-Kendall method. Subsequently, the Hurst index was used to simulate the vegetation NPP persistence in the study area. Finally, the driving mechanisms of vegetation NPP changes in the study area were explored using partial correlation coefficients and residual analysis. The results demonstrated that: 1) The annual average NPP in the basin showed a fluctuating and increasing trend from 273.54 to 718.30 g C/m2·yr1 during 2000 to 2020, and except for autumn, all seasons (spring, summer, winter) and annual NPP series showed an upward trend. The spatial distribution of vegetation NPP is characterized by an asymmetrical horseshoe shape in general, with high values in the west, south and east parts, and low values in the Lake area; 2) During the study years, about 84.38 % of the basin area showed a significant and extremely significant increase of NPP; 3) The future trend of vegetation NPP in the basin shows that the area of decrease (22.79 %) is more than the area of increase (11.35 %). The Chang-Zhu-Tan urban agglomeration will generally show a continuous and extremely significant reduction trend, while parts of Dongting Lake will show a continuous and extremely significant increase trend; 4) The correlation between solar radiation and NPP is stronger than the other two meteorological factors (precipitation and temperature). Temperature has a significant inhibitory effect on NPP, solar radiation has a significant promoting effect on NPP, and the effect of precipitation on NPP was relatively complicated. The relative importance of meteorological factors on vegetation NPP was ranked as follows: solar radiation > precipitation > temperature. The impacts of both climate change and human activities on NPP changes showed great spatial variability, and the positive contributions (89.81 % for climate change and 82.87 % for human activities) were both greater than the negative contributions.  相似文献   

7.
Utilizing freely available MODIS NDVI and Natural color imageries of 250 m spatial resolution produced by NASA, an experiment was made to map land-cover and its change with an emphasis on vegetation cover in southeastern Sri Lanka, which plays a vital role for control of green house gas. For the change detection purpose, 1987 land cover map made by present authors from Landsat MSS image and extensive ground truth survey data was used as the base map. The result of the experiment shows that MODIS data are useful to make a land cover map of 250 m spatial resolution for tropical areas with high cloud coverage like Sri Lanka. It was found that the forest cover decrease amounted as large as 21% in 19 years time span in southeastern Sri Lanka, the prominent forest region of the country. On the other hand homestead/vegetation and mixed vegetation/scrub dominant categories increased by 13.7% and 7.1%, respectively. These changes are considered due to a large clearance of forest areas for agriculture and building houses to accommodate increasing inhabitants.  相似文献   

8.
Drought is an important natural disaster that causes devastating impacts on the ecosystem, livestock, environment, and society. So far, various remote-sensing methods have been developed to estimate drought conditions, each of which has advantages and restrictions. This study aims to monitor the real-time drought indices at the field scales via the integration of various earth observations. Our proposed method consists of two steps. In the first step, the relationships between long-term standardized precipitation indices (SPI) derived from PERSIANN-CDR rainfall data and two drought-dependent parameters derived from MODIS products, including normalized NDVI and soil-air temperature gradient, are obtained at the spatial resolution of PERSIANN-CDR grid (approximately 25 km). As the next step, the corresponding relationships are applied to estimate the drought index maps at the spatial resolution of MODIS products (1 km). Numerous analyses are carried out to evaluate the proposed method. The results revealed that, from various drought indices, including SPIs of different timescales (1, 3, 6, and 12-months), SPI-3 and SPI-6 are more appropriate to the proposed method in terms of correlation with temperature and vegetation parameters. The findings also demonstrate the competency of the proposed method in estimating SPI indices with average RMSE 0.67 and the average correlation coefficient of 0.74.  相似文献   

9.
This report summarizes a major paper reviewing Canadian work using Landsat imagery for studying changes in lakes and in coastal environments. The nature of environmental change is discussed. For lakes and coastal environments, it is suggested that change is either seasonal, long term, short term or constant; examples of each are given. There is also an important distinction between natural and man-induced change. Outlined in the paper are studies of the filling of the LG2 Reservoir, water level and vegetation changes in the Peace-Athabasca Delta, possible vegetation changes due to the construction of Roberts Bank Port near Vancouver and measurement of surface suspended sediment concentration in the Bay of Fundy using chromaticity analysis.  相似文献   

10.
In recent years, land surface temperature (LST) has become critical in environmental studies and earth science. Remote sensing technology enables spatiotemporal monitoring of this parameter on large scales. This parameter can be estimated by satellite images with at least one thermal band. Sentinel-3 SLSTR data provide LST products with a spatial resolution of 1 km. In this research, direct and indirect validation procedures were employed to evaluate the Sentinel-3 SLSTR LST products over the study area in different seasons from 2018 to 2019. The validation method was based on the absolute (direct) evaluation of this product with field data and comparison (indirect) evaluation with the MODIS LST product and the estimated LST using the non-linear split-window (NSW) algorithm. Also, two emissivity estimation methods, (1) NDVI thresholding method (NDVI-THM) and (2) classification-based emissivity method (CBEM), were used to estimate the LST using the NSW method according to the two thermal bands of Sentinel-3 images. Then, the accuracy of these methods in estimating LST was evaluated using field data and temporal changes of vegetation, which the NDVI-THM method generated better results. For indirect evaluation between the Sentinel-3 LST product, MODIS LST product, and LST estimated using NSW, four filters based on spatial and temporal separates between pairs of pixels and pixel quality were used to ensure the accuracy and consistency of the compared pairs of a pixel. In general, the accuracy results of the LST products of MODIS and Sentinel-3, and LST estimated using NSW showed a similar trend for LST changes during the seasons. With respect to the two absolute and comparative validations for the Sentinel-3 LST products, summer with the highest values of bias (?1.24 K), standard deviation (StDv = 2.66 K), and RMSE (2.43 K), and winter with the lowest ones (bias of 0.14 K, StDv of 1.13 K, and RMSE of 1.12 K) provided the worst and best results for the seasons in the period of 2018–2019, respectively. According to both absolute and comparative evaluation results, the Sentinel-3 SLSTR LST products provided reliable results for all seasons on a large temporal and spatial scale over our studied area.  相似文献   

11.
There are 1.2 × 106 ha of heath vegetation on the island of Newfoundland. Some of this heath is of climatic or edaphic origin, but most of it has originated because of increasing frequency of fire and cutting following settlement of the area by European man since discovery in 1497. Modern fire suppression has been successful in reducing the size and frequency of fires in the last 40–50 years, however, natural forest regeneration is still not apparent in the anthropogenic heath.The ecotone between the forest-heath appears to have stabilized following an initial influx of trees up to 10 m from the forest margin. Analysis of the age structure of the invading trees demonstrates that after a period of 10 years new tree establishment is reduced to less than 5 percent of the total establishment. This decreased establishment over time is correlated with a coincidental increase in dwarf shrub heath vegetation that creates undesirable germination conditions for tree seed entering the heath. Even with ideal seedbed conditions for germination, the short distance of seed dispersal (10–20 m) and the excessively slow growth rates of trees, limit forest encroachment to a rate of 10 cm per year. At this rate more than a thousand years would be required for successional encroachement of even 2–3 ha of heathland. The successional status of the heath vegetation is considered to be permanent even though traditional concepts of secondary succession and chronosequence would predict the recovery of a climax forest vegetation.  相似文献   

12.
Land surface temperature (LST) is an important factor in global change studies, heat balance and as control for climate change. A comparative study of LST over parts of the Singhbhum Shear Zone in India was undertaken using various emissivity and temperature retrieval algorithms applied on visible and near infrared (VNIR), and thermal infrared (TIR) bands of high resolution Landsat-7 ETM+ imagery. LST results obtained from satellite data of October 26, 2001 and November 2, 2001 through various algorithms were validated with ground measurements collected during satellite overpass. In addition, LST products of MODIS and ASTER were compared with Landsat-7 ETM+ and ground truth data to explore the possibility of using multi-sensor approach in LST monitoring. An image-based dark object subtraction (DOS3) algorithm, which is yet to be tested for LST retrieval, was applied on VNIR bands to obtain atmospheric corrected surface reflectance images. Normalized difference vegetation index (NDVI) was estimated from VNIR reflectance image. Various surface emissivity retrieval algorithms based on NDVI and vegetation proportion were applied to ascertain emissivities of the various land cover categories in the study area in the spectral range of 10.4–12.5 μm. A minimum emissivity value of about 0.95 was observed over the reflective rock body with a maximum of about 0.99 over dense forest. A strong correlation was established between Landsat ETM+ reflectance band 3 and emissivity. Single channel based algorithms were adopted for surface radiance and brightness temperature. Finally, emissivity correction was applied on ‘brightness temperature’ to obtain LST. Estimated LST values obtained from various algorithms were compared with field ground measurements for different land cover categories. LST values obtained after using Valor’s emissivity and single channel equations were best correlated with ground truth temperature. Minimum LST is observed over dense forest as about 26 °C and maximum LST is observed over rock body of about 38 °C. The estimated LST showed that rock bodies, bare soils and built-up areas exhibit higher surface temperatures, while water bodies, agricultural croplands and dense vegetations have lower surface temperatures during the daytime. The accuracy of the estimated LST was within ±2 °C. LST comparison of ASTER and MODIS with Landsat has a maximum difference of 2 °C. Strong correlation was found between LST and spectral radiance of band 6 of Landsat-7 ETM+. Result corroborates the fact that surface temperatures over land use/land cover types are greatly influenced by the amount of vegetation present.  相似文献   

13.
This paper proposes a particular approach to assess information about soil degradation, based on a methodology to calculate soil color from NOAA/AVHRR data. As erosive processes change physical and chemical properties of the soil, altering, consequently, the superficial color, monitoring the change in color over time can help to identify and analyze those processes. A relationship among the soil color (described in the Munsell Color System, i.e., in terms of Hue, Value and Chroma), vegetation indices, surface temperature and emissivity has been established, which is based on the theoretical model. The methodology has three main phases: determination of the regression models among soil color and vegetation indices, emissivity and surface temperature; generation of digital soil color models; and statistical evaluation of the estimated color. The tests showed that the methodology is efficient in determining soil color using the various vegetation indices (i.e., Normalized vegetation index NDVI, Modified soil adjusted vegetation index MSAVI). One vegetation index, i.e., Purified adjusted vegetation index (PAVI) is proposed to subsidies the effect of vegetation over the soil. Best results were obtained for the Hue color component. To further test the methodology, the estimated digital color models were compared with the characteristic color of soil classes in the test area. The results of this application confirmed the methodology’s capacity to determine the soil color from NOAA/AVHRR data. This type of study is quite helpful to know the erosion of soil as well as some abrupt change in soil due to natural hazards by space borne or air-borne sensors.  相似文献   

14.
Updated information of rubber plantations is essential for assessing socioeconomic and environmental impacts, especially in the emerging region of northern tropics. Here, a phenological method was modified to detect rubber plantations using Landsat Operational Land Imager (OLI) imagery in Phongsaly Province of northern Laos, where it begun a rubber boom in the mid-2000s due to geo-economic cooperation. It highlighted the landscape and pixel differences of deciduous rubber plantations in the tri-temporal phases (i.e., pre-defoliation, defoliation, and foliation) during the dry season due to phenological changes. Six commonly used vegetation indices (VIs), including the Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI), Land Surface Water Index (LSWI), Atmospherically Resistant Vegetation Index (ARVI), Normalized Burn Ratio (NBR and NBR2) derived from OLI imagery during 2013–2016 were compared to determine the most suitable VI for discriminating the phenological differences of rubber plantations from natural forests. Then, the Differences of Normalized Burn Ratio (DNBR) was applied to generate the 30 m map of rubber plantations in 2016, by combining two masks of Landsat-derived forest and suitable elevation for rubber trees cultivation. The resultant map of rubber plantations had a classification accuracy of 93.7% and the Kappa coefficient of 0.848. Our study demonstrated the usefulness of the Landsat-derived tri-temporal phenological DNBR approach in an emerging region of northern Laos, despite requiring more scenes compared with single- and double temporal window methods.  相似文献   

15.
The Nile River Basin (NRB) is facing extreme demand for its water resources due to an alarming increase in population and the changing climate. The NRB is not compatible with ground-based in-situ observations owing to its large basin area size and limited hydrological data access from basin countries. Thus, it lends itself to remotely sensed approaches with high spatial resolution and extended temporal coverage. The Gravity Recovery and Climate Experiment (GRACE) avails a unique opportunity to investigate the changes in key components of terrestrial water storage (TWS). GRACE TWS solutions have specific tuning parameters and processing strategies that result in regionally specific variations and error patterns. We explored the TWS time series spatiotemporal changes, trends, uncertainties, and signal-to-noise ratio among different GRACE TWS data. We had also investigated the key terrestrial water storage components such as surface water, soil moisture, and groundwater storage changes. The results show that GRACE spherical harmonic solutions' uncertainty is higher than the mass concentration (mascon) over the NRB, and the Center for Space Research-mascons had the best performance. The evapotranspiration correlation (R2 = 0.85) has the highest correlation with GRACE’s TWS, whereas the normalized difference vegetation index (R2 = 0.82) has the second highest correlation. Notably, significant long-term (2003–2017) negative groundwater and soil moisture trends demonstrate a potential depletion of the NRB. Despite an increase in precipitation and the TWS time series, the rate of decline increased rapidly after 2008, thereby indicating the possibility of human-induced change (e.g. for irrigation purposes). Therefore, the results of this study provide a guide for future studies related to hydro-climatic change over the NRB and similar basins.  相似文献   

16.
The global developments of the stratospheric events (~20–50 km) are briefly described using balloon and satellite data. Winds data from L.F. drift (52°N, 15°E, Europe) for heights of 90–100 km, and from M.F. radar (52°N, 107°W, Canada) for heights of 60–110 km are then compared with the stratospheric morphology.Data for 1982/3 and 1983/4 show that the planetary wave activity and warmings produced strong westward and southward perturbations in the radar winds. Satellite data from 0.1, 0.01 hPa are consistent with these winds; and also show smaller scale structures in the mesosphere than the stratosphere. The semi-diurnal tide responded strongly to the atmospheric disturbances in Europe and Canada: for the latter vertical wavelength changes occurred for heights of 70–100 km. However the correlation between these tidal fluctuations was not high indicating that the tidal adjustments were continental rather than hemispheric.  相似文献   

17.
Information on rice growing areas is important for policymakers to devise agricultural plans. This research explores the monitoring of rice cropping intensity in the upper Mekong Delta, Vietnam (from 2001 to 2007) using time-series MODIS NDVI 250-m data. Data processing includes three steps: (1) noise is filtered from the time-series NDVI data using empirical mode decomposition (EMD); (2) endmembers are extracted from the filtered time-series data and trained in a linear mixture model (LMM) for classification of rice cropping systems; and (3) classification results are verified by comparing them with the ground-truth and statistical data. The results indicate that EMD is a good filter for noise removal from the time-series data. The classification results confirm the validity of LMM, giving an overall accuracy of 90.1% and a Kappa coefficient of 0.7. The lowest producer and user accuracies were associated with single crop rain-fed rice class due to the mixed pixel problems. A strong yearly correlation at the district level was revealed in the MODIS-derived areas (R2 ? 0.9). Investigation of interannual changes in rice cropping intensity from 2001 to 2007 showed a remarkable conversion from double to triple crop irrigated rice from 2001 to 2003, especially in the Thoai Son and Phu Tan districts. A big conversion from triple crop rice back to double crop rice cultivation was also observed in Phu Tan from 2005 to 2006. These changes were verified by visual interpretation of Landsat images and examination of NDVI profiles.  相似文献   

18.
太阳黑子的分维数及其可预报性   总被引:1,自引:0,他引:1  
本文利用1932年至1982年的逐日国际太阳黑子相对数的9日平均时间序列来分析太阳黑子演变的动力学特性,发现太阳黑子是至少需要七个独立变量来描述的复杂动力系统,其关联维数为6.3±0.1,二阶Renyi熵和最大Lyapunov指数分别为(0.37±0.02)bit/yr和(0.34±0.01)bit/yr,它们表征的可预报时间尺度,即初始误差增长一倍所需的时间分别为(2.7±0.2)年和(2.9±0.l)年。   相似文献   

19.
Recent studies of the vegetation fluorescence show that it can be successfully used as an intrinsic indicator of plant photosynthetic activity. With respect to the vegetation spectral reflectance, the chlorophyll (Chl) fluorescence is more specific as an observable of basic biophysical processes in the plant cells. Laser induced fluorescence is widely used in near range remote sensing, but it is not suitable for the global monitoring of vegetation. Decades of active fluorometry studies have collected useful information of leaf reaction to natural and anthropogenic stress. Still the passive fluorescence, the one that could be registered from satellite orbit has still to prove its advantage over widely used reflectance signature. The weakness of the signal and the lack of experience with passive fluorescence measurements require extensive technical, theoretical and experimental studies. New imaging fluorometres are to be designed for measuring steady state fluorescence in controlled and natural conditions.

In order to compare reflectance and steady state fluorescence sensitivity to stress impact, a set of experiments have been conducted under controlled illumination conditions in a bio-chamber, designed by the author’s team. The equipment allows plant vitality to be monitored both by passive fluorescence and spectral reflectance imaging. Different types of stress factors (heat and drought stress, acid impact) were investigated to demonstrate equipments ability in monitoring changes of fluorescence signal. Selected fluorescence images of foliage illustrate an early detection of plant dysfunction and the temporal and spatial spreading of the stress impact. Analysis shows that fluorescence imaging of green plants can be developed as a highly effective early warning remote sensing method, which could have application for an ecosystems’ monitoring along with high-spectral reflectance imagery.  相似文献   


20.
The purpose of this paper is to present the results of research to identify, by satellite imagery, parameters of the environment affecting health on Earth. Thus, we suggest expanding the application of space technology to preventive medicine, as a new field in the peaceful uses of outer space. The scope of the study includes all parts of the environment, natural and man-made, and all kinds of protection of life: human, animal and vegetation health. The general objective is to consider and classify those factors, detectable from space, that affect or are relevant to health and may be found in the air, water, sea, soil, land, vegetation, as well as those linked to climate, industry, energy production, development works, irrigation systems, and human settlements. The special objective is the classification of environmental factors detectable from space, that are linked to communicable or chronic endemic diseases or health problems. The method of identifying the factors affecting health was the parallel study of environmental epidemiological and biological parameters. The role of environmental factors common to both human and animal populations is discussed. Conclusive findings are formulated and possible applications, both scientific and practical, in other sectors are also discussed.  相似文献   

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

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