首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 453 毫秒
1.
Spatial differences in drought proneness and intensity of drought caused by differences in cropping patterns and crop growing environments within a district indicate the need for agricultural drought assessment at disaggregated level. The objective of this study is to use moderate resolution satellite images for detailed assessment of the agricultural drought situation at different administrative units (blocks) within a district. Monthly time composite NDVI images derived from moderate resolution AWiFS (60 m) and WiFS (180 m) images from Indian Remote Sensing satellites were analysed along with ground data on rainfall and crop sown areas for the kharif seasons (June – November) of 2002 (drought year), 2004 (early season drought) and 2005 (good monsoon year). The impact of the 2002 meteorological drought on crop area in different blocks of the district was assessed. The amplitude of crop condition variability in a severe drought year (2002) and a good year (2005) was used to map the degree of vulnerability of different blocks in the district to agricultural drought. The impact of early season deficit rainfall in 2004 on the agricultural situation and subsequent recovery of the agricultural situation was clearly shown. Agricultural drought assessment at disaggregated level using moderate resolution images is useful for prioritizing the problem areas within a district to undertake, in season drought management plans, such as alternate cropping strategies, as well as for end of the season drought relief management actions. The availability of ground data on rainfall, cropping pattern, crop calendar, irrigation, soil type etc., is very crucial in order to interpret the seasonal NDVI patterns at disaggregated level for drought assessment. The SWIR band of AWiFS sensor is a potential data source for assessing surface drought at the beginning of the season.  相似文献   

2.
Developing a robust drought monitoring tool is vital to mitigate the adverse impacts of drought. A drought monitoring system that integrates multiple agrometeorological variables into a single drought indicator is lacking in areas such as Ethiopia, which is extremely susceptible to this natural hazard. The overarching goal of this study is to develop a combined drought indicator (CDI-E) to monitor the spatial and temporal extents of historic agricultural drought events in Ethiopia. The CDI-E was developed by combining four satellite-based agrometeorological input parameters – the Standardized Precipitation Index (SPI), Land Surface Temperature (LST) anomaly, Standardized Normalized Difference Vegetation Index (stdNDVI) and Soil Moisture (SM) anomaly – for the period from 2001 to 2015. The method used to combine these indices is based on a quantitative approach that assigns a weight to each input parameter using Principal Component Analysis (PCA). The CDI-E results were evaluated using satellite-based gridded rainfall (3-month SPI) and crop yield data for 36 intra-country crop growing zones for a 15-year period (2001 to 2015). The evaluation was carried out for the main rainfall season, Kiremt (June-September), and the short rainfall season, Belg (February-May). The results showed that moderate to severe droughts were detected by the CDI-E across the food insecure regions reported by FEWS NET during Kiremt and Belg rainfall seasons. Relatively higher correlation coefficient values (r > 0.65) were obtained when CDI-E was compared with the 3-month SPI across the majority of Ethiopia. The spatial correlation analyses of CDI-E and cereal crop yields showed relatively good correlations (r > 0.5) in some of the crop growing zones in the northern, eastern and southwestern parts of the country. The CDI-E generally mapped the spatial and temporal patterns of historic drought and non-drought years and hence the CDI-E could potentially be used to develop an agricultural drought monitoring and early warning system in Ethiopia. Moreover, decision makers and donors may potentially use CDI-E to more accurately monitor crop yields across the food-insecure regions in Ethiopia.  相似文献   

3.
The untimely onset and uneven distribution of south-west monsoon rainfall lead to agricultural drought causing reduction in food-grain production with high vulnerability over semi-arid tract (SAT) of India. A combined deficit index (CDI) has been developed from tri-monthly sum of deficit in antecedent rainfall and deficit in monthly vegetation vigor with a lag period of one month between the two. The formulation of CDI used a core biophysical (e.g., NDVI) and a hydro-meteorological (e.g., rainfall) variables derived using observation from Indian geostationary satellites. The CDI was tested and evaluated in two drought years (2009 and 2012) within a span of five years (2009–2013) over SAT. The index was found to have good correlation (0.49–0.68) with standardized precipitation index (SPI) computed from rain-gauge measurements but showed lower correlation with anomaly in monthly land surface temperature (LST). Significant correlations were found between CDI and reduction in agricultural carbon productivity (0.67–0.83), evapotranspiration (0.64–0.73), agricultural grain yield (0.70–0.85). Inconsistent correlation between CDI and ET reduction was noticed in 2012 in contrast to consistent correlation between CDI and reduction in carbon productivity both in 2009 and 2012. The comparison of CDI-based drought-affected area with those from existing operational approach showed 75% overlapping regions though class-to-class matching was only 40–45%. The results demonstrated that CDI is a potential indicator for assessment of late-season regional agricultural drought based on lag-response between water supply and crop vigor.  相似文献   

4.
Evapotranspiration (ET) is a vital process in land surface atmosphere research. In this study, Surface Energy Balance Algorithm for Land (SEBAL) for the assessment of ET (for 23 December 2010, 8 January 2011, 24 January 2011, 9 February 2011, 25 February 2011, 29 March 2011 and 14 April 2011) from LANDSAT7-ETM+ and validation with Lysimeter data set is illustrated. It is based on the evaporative fraction concept, and it has been applied to LANDSAT7-ETM + (30 m resolution) data acquired over the Indian Agricultural Research Institute’s agricultural farm land. The ET from SEBAL was compared with Lysimeter ET using four statistical tests (root-mean-square error (RMSE), relative root-mean-square error (R-RMSE), mean absolute error (MAE), and normalized root-mean square error (NRMSE)), and each test showed a good correlation between the predicted and observed ET values. Results from this study revealed that the RMSE of crop-growing period was 0.51 mm d?1 for ETSEBAL, i.e. ETSEBAL having good accuracy with respect to observed ETLysimeter. Results were also validated using R-RMSE test, which also proved that ETSEBAL data are having good accuracy with respect to observed ETLysimeter as R-RMSE of crop-growing period is 0.19 mm d?1. MAE (0.19), NRMSE (0.21) and r2 (0.91) tests indicated that model prediction is significant, and model can be effectively used for the estimation of ET from SEBAL as input of remote sensing data sets. Finally, the SEBAL has been useful for remote agricultural land where ground-based data (Lysimeter data) are not available for daily ET (ET24 h) estimation. The temporal study of the ET24 h values analysed has revealed that the highest ET24 h values are owing to the higher development (high greenness) of crops, whereas the lower values are related to the lower development (low greenness) or null crop.  相似文献   

5.
Numerous efforts have been made to develop various indices using remote sensing data such as normalized difference vegetation index (NDVI), vegetation condition index (VCI) and temperature condition index (TCI) for mapping and monitoring of drought and assessment of vegetation health and productivity. NDVI, soil moisture, surface temperature and rainfall are valuable sources of information for the estimation and prediction of crop conditions. In the present paper, we have considered NDVI, soil moisture, surface temperature and rainfall data of Iowa state, US, for 19 years for crop yield assessment and prediction using piecewise linear regression method with breakpoint. Crop production environment consists of inherent sources of heterogeneity and their non-linear behavior. A non-linear Quasi-Newton multi-variate optimization method is utilized, which reasonably minimizes inconsistency and errors in yield prediction.  相似文献   

6.
遥感技术具备实时快速、时空连续、广覆盖尺度等独特优势,在全球气候恶化大背景下,利用遥感干旱监测方法相比于传统地面监测手段,能够提供实时、准确、稳定的旱情信息,辅助科学决策。目前常用遥感旱情监测方法大多依赖全域性数学模型建模,假定了旱情模式的空间平稳特性,因而难以准确反映旱情模式的局部差异特征。本文提出利用地理加权回归模型GWR (Geographically Weighted Regression),考虑旱情模式的空间非平稳特性,综合多种遥感地面旱情监测指数,以实现传统全域旱情监测模型的局部优化。以美国大陆为研究区,监测2002年—2011年共10年的旱情状态。研究表明,GWR模型能够提供空间变化的局部最佳估计模型参数,监测结果更加吻合标准美国旱情监测USDM (U.S Drought Monitor)验证数据,且与地面实测值的最高相关系数R达到0.8552,均方根误差RMSE达到0.972,显著优于其他遥感旱情监测模型。GWR模型具备空间非平稳探测优势,实现了旱情模式的局部精细探测,能够显著提升遥感旱情监测精度,具备较好的应用前景。  相似文献   

7.
This study investigated rice cropping practices and rice growing areas in the Vietnamese Mekong Delta using MODIS 250 × 250 m normalized difference vegetation index (NDVI) data acquired during the 2002 and 2007 rice cropping seasons. Data processing was conducted in five main steps: (1) constructing time-series MODIS NDVI data; (2) noise filtering of the time-series MODIS NDVI data using empirical mode decomposition (EMD); (3) extracting and evaluating phenological rice training patterns from the smooth time profiles of NDVI; (4) classifying rice cropping systems using support vector machines (SVMs); and (5) conducting an error analysis using ground reference data and government rice statistics. The results indicated that EMD was an efficient filter for noise removal in the time-series MODIS NDVI data. The filtered temporal NDVI profile characterized the distinct behaviors of the rice cropping systems. The estimated sowing and harvesting dates were compared with the field-survey data and indicated root mean square error (RMSE) values of 7.5 and 8.2 days, respectively. The comparison results between the 2002 classification map and the ground reference data indicated that the overall accuracy for the 2002 data was 92.9% with a Kappa coefficient of 0.89, while in 2007 these values were 93.8% and 0.90, respectively. At the district level, there was good agreement between the MODIS-based estimated areas and government rice statistics for 2002 and 2007 (R 2 ≥ 0.85). An investigation of changes in cropping practices from 2002 to 2007 showed that 12.9% of the area used for double-cropped irrigated rice in 2002 had been converted to triple-cropped irrigated rice by 2007, whereas 27.4% of the area used for triple-cropped irrigated rice in 2002 had been converted to double-cropped irrigated rice by 2007.  相似文献   

8.
Rice crop occupies an important aspect of food security and also contributes to global warming via GHGs emission. Characterizing rice crop using spatial technologies holds the key for addressing issues of global warming and food security as different rice ecosystems respond differently to the changed climatic conditions. Remote sensing has become an important tool for assessing seasonal vegetation dynamics at regional and global scale. Bangladesh is one of the major rice growing countries in South Asia. In present study we have used remote sensing data along with GIS and ancillary map inputs in combination to derive seasonal rice maps, rice phenology and rice cultural types of Bangladesh. The SPOT VGT S10 NDVI data spanning Aus, Aman and Boro crop season (1st May 2008 to 30th April 2009) were used, first for generating the non-agriculture mask through ISODATA clustering and then to generate seasonal rice maps during second classification. The spectral rice profiles were modelled and phenological parameters were derived. NDVI growth profiles were modelled and crop calendar was derived. To segregate the rice cultural types of Bangladesh into IPCC rice categories, we used elevation, irrigated area, interpolated rainfall maps and flood map through logical modelling in GIS. The results indicated that the remote sensing derived rice area was 9.99 million ha as against the reported area of 11.28 million ha. The wet and dry seasons accounted for 64% and 36 % of the rice area, respectively. The flood prone, drought prone and deep water categories account for 7.5%, 5.56% and 2.03%, respectively. The novelty of current findings lies in the spatial outcome in form of seasonal and rice cultural type maps of Bangladesh which are helpful for variety of applications.  相似文献   

9.
Pre-harvest crop production forecast has been successfully provided by remote sensing technique. However, the probability to get cloud-free optical remote sensing data during kharif season is poor. Microwave data having the capability to penetrate cloud is used in the absence of cloud free optical remote sensing data. Yield models in broad band frequency range are in development stage. Meteorological yield models are developed and predicted yield is combined with area estimated by remote sensing data to provide rice production forecast. This paper describes the methodology adopted for improving the predictability of rice yield before harvest of the crop in Bihar province by taking into consideration meteorological parameters during its growth cycle upto October. Models developed using fortnightly meteorological data have been found to give reasonably fair indications of expected yield of rice in advance of harvest. The yield predictions have been made based on meteorological data and effective rainfall based on water requirement calculations representing a group of districts under similar agro-climatic zones, which could be further improved by incorporating meteorological data of individual districts within each group.  相似文献   

10.
Large scale adoption of input intensive rice–wheat cropping system in the centrally located Jalandhar district of Indian Punjab has led to over-exploitation of ground water resources, intensive use of chemical fertilizers and deterioration of soil health. To overcome these shortfalls, in the present study, agricultural area diversification plan has been generated from agricultural area and crop rotation maps derived from remote sensing data (IRS P6-AWiFS and RADARSAT ScanSAR) along with few agro-physical parameters in GIS environment. Cropping system indices (area diversity, multiple cropping and cultivated land utilization) were also worked out from remote sensing data .Analysis of remote sensing data (2004–05) revealed that rice and wheat individually remained the dominant crops, occupy 57.8% and 64.9% of total agricultural area (TAA), respectively. Therefore, in the diversified plan, it is suggested that at least 39% of the current 40% TAA under rice–wheat rotation should be replaced by other low water requiring, high value and soil enriching crops, particularly in coarse textured alluvial plain having good quality ground water zones with low annual rainfall(<700 mm). This will reduce water requirement to the tune of 15,660 cm depth while stabilizing the production and profitability by crop area diversification without further degradation of natural resources.  相似文献   

11.
Evapotranspiration is a key parameter for water stress assessment as it is directly related to the moisture status of the soil-vegetation system and describes the moisture transfer from the surface to the atmosphere. With the launch of the Meteosat Second Generation geostationary satellites and the setup of the Satellite Application Facilities, it became possible to operationally produce evapotranspiration data with high spatial and temporal evolution over the entire continents of Europe and Africa. In the frame of this study we present an evaluation of the potential of the evapotranspiration (ET) product from the EUMETSAT Satellite Application Facility on Land Surface Analysis (LSA-SAF) for drought assessment and monitoring in Europe.To assess the potential of this product, the LSA-SAF ET was used as input for the ratio of ET to reference evapotranspiration (ET0), the latter estimated from the ECMWF interim reanalysis. In the analysis two case studies were considered corresponding to the drought episodes of spring/summer 2007 and 2011. For these case studies, the ratio ET/ET0 was compared with meteorological drought indices (SPI, SPEI and Sc-PDSI for 2007 and SPI for 2011) as well as with the anomalies of the fraction of absorbed photosynthetic active radiation (fAPAR) derived from remote sensing data. The meteorological and remote sensing indicators were taken from the European Drought Observatory (EDO) and the CARPATCLIM climatological atlas.Results show the potential of ET/ET0 to characterize soil moisture variability, and to give additional information to fAPAR and to precipitation distribution for drought assessment. The main limitations of the proposed ratio for drought characterization are discussed, including options to overcome them. These options include the use of filters to discriminate areas with a low percentage vegetation cover or areas that are not in their growing period and the use of evapotranspiration without water restriction (ETwwr), obtained as output of the LSA-SAF model instead of ET0. The ET/ETwwr ratio was tested by comparing its accumulated values per growing period with the winter wheat yield values per country published by Eurostat. The results point to the potential of using the remote sensing based LSA-SAF evapotranspiration and the ET/ETwwr ratio for vegetation monitoring at large scale, especially in areas where data is generally lacking.  相似文献   

12.
In semiarid regions the occurrence of alternating long drought and heavy rainfall periods directly impacts water availability, affecting human water supply, agriculture development and the provision of ecosystem services. Because of that, research on the water input and output fluxes at the basin scale is of paramount importance. In this sense, rainfall-evapotranspiration (ET) dynamics play a critical role in water, soil and vegetation interactions, in hydrometeorological modelling and in the energy fluxes dynamics of semiarid regions. Therefore, the objective of this study was to quantify daily ET during a wet year and a dry year in a watershed located in the Brazilian Semiarid, by using remote sensing data and formulations based on the Simplified Surface Energy Balance Index (S-SEBI) and the Simplified Surface Energy Balance (SSEB) algorithms. Land surface temperature, albedo and NDVI data from MODIS sensor and solar radiation data from weather stations located in the basin were used. Rainfall analysis indicated 2009 and 2012 as being representatives of anomalously wet and dry years respectively, which were selected for the quantification of ET. The proposed algorithm was adjusted and verified with data from a flux tower equipped with eddy covariance system. Daily remote sensing ET estimates showed good agreement with observed values (RMSE = 0.79 mm.d−1) and the annual ET relative error was of 7.7% (35.4 mm.year−1). Results showed that the native vegetation can delay its dormant state for five months during wet years. During the wet year, ET differences between land cover classes were less noticeable due to soil saturation and the urgency of vegetated surfaces to meet their physiological needs. In dry year, however, differences were more evident, with bare soil presenting lower ET rates and vegetation classes showing higher ET values.  相似文献   

13.
Human diets strongly rely on wheat, maize, rice and soybean; research on the potential crop productivity of these four main crops could provide the basis for increasing global crop yields. The evaluation model of realistic potential crop productivity based on remote sensing and agro-ecological zones was proposed in this study to provide reliable reference data for world food security. The statistical data on these four main crops yields were obtained from the FAO. The model was used to investigate the potential production of four staple crops in the world. The distributions of the realistic potential productivity of four staple crops (winter wheat, maize, rice and soybean) were produced. In the main producing countries of the four staple crops, statistical analysis was conducted on the realistic potential productivity (RPP) of the four staple crops, the highest productivity (HP) during the period 1983–2011 and the gap between RPP and HP.  相似文献   

14.
The paper presents a detailed understanding of nitrogenous fertilizer use in Indian agriculture and estimation of seasonal nitrogen loosses from rice crop in Indo-Gangetic plain region, the ‘food bowl’ of the Indian sub-continent. An integrated methodology was developed for quantification of different forms of nitrogen losses from rice crop using remote sensing derived inputs, field data of fertilizer application, collateral data of soil and rainfall and nitrogen loss coefficients derived from published nitrogen dynamics studies. The spatial patterns of nitrogen losses in autumn or ‘kharif’ and spring or ‘rabi’ season rice at 1 × 1 km grid were generated using image processing and GIS. The nitrogen losses through leaching in form of urea-N, ammonium-N (NH4-N) and nitrate-N (NO3-N) are dominant over ammonia volatilization loss. The study results indicate that nitrogen loss through leaching in kharif and rabi rice is of the order of 34.9% and 39.8% of the applied nitrogenous fertilizer in the Indo-Gangetic plain region. This study provides a significant insight to the role of nitrogenous fertilizer as a major non-point source pollutant from agriculture.  相似文献   

15.
The study aims to simulate the peri-urban growth dynamics in a growing region of India using Weights of Evidence (WOE) based cellular automata model. The growth process was expressed as a function of four causative variables corresponding to which seven data layers were generated in a Geographic Information Systems environment. The model was calibrated for the period 2000–2005 using Kappa indices and fuzzy set theory based two way comparison method. The Kappa value was 0.7, while the value of Klocation and Khisto were 0.81 and 0.93, respectively. The fuzzy similarity values increased for small to large neighbourhood sizes which showed that the model was able to simulate the contiguous and dense growth. However, for dispersed and isolated growth the model showed less accuracy. The model was validated for period 2005–2010 and revealed a Kappa value of 0.88, while value of Klocation and Khisto were 0.91 and 0.96, respectively.  相似文献   

16.
Monitoring crop conditions and forecasting crop yields are both important for assessing crop production and for determining appropriate agricultural management practices; however, remote sensing is limited by the resolution, timing, and coverage of satellite images, and crop modeling is limited in its application at regional scales. To resolve these issues, the Gramineae (GRAMI)-rice model, which utilizes remote sensing data, was used in an effort to combine the complementary techniques of remote sensing and crop modeling. The model was then investigated for its capability to monitor canopy growth and estimate the grain yield of rice (Oryza sativa), at both the field and the regional scales, by using remote sensing images with high spatial resolution. The field scale investigation was performed using unmanned aerial vehicle (UAV) images, and the regional-scale investigation was performed using RapidEye satellite images. Simulated grain yields at the field scale were not significantly different (= 0.45, p = 0.27, and p = 0.52) from the corresponding measured grain yields according to paired t-tests (α = 0.05). The model’s projections of grain yield at the regional scale represented the spatial grain yield variation of the corresponding field conditions to within ±1 standard deviation. Therefore, based on mapping the growth and grain yield of rice at both field and regional scales of interest within coverages of a UAV or the RapidEye satellite, our results demonstrate the applicability of the GRAMI-rice model to the monitoring and prediction of rice growth and grain yield at different spatial scales. In addition, the GRAMI-rice model is capable of reproducing seasonal variations in rice growth and grain yield at different spatial scales.  相似文献   

17.
The goal of this research was to conduct an initial investigation into whether a time-series NDVI reference curve library for crops over a growing season for one year could be used to map crops for a different year. Time-series NDVI libraries of curves for 2001 and 2005 were investigated to ascertain whether or not the 2001 dataset could be used to map crops for 2005. The 2005 16-day composite MODIS 250 m NDVI data were used to extract NDVI values from 1,615 field sites representing alfalfa, corn, sorghum, soybeans, and winter wheat. A k-means cluster analysis of NDVI values from the field sites was performed to identify validation sites with time-series NDVI spectral profiles characteristic of the major crop types grown in Kansas. After completing the field site refinement process, there were 1,254 field sites retained for further analysis, referred to as "final" field sites. The methods employed to evaluate whether the MODIS-based NDVI profiles for major crops in Kansas are stable from year-to-year involved both graphical and statistical analyses. First, the time-series NDVI values for 2005 from the final field sites were aggregated by crop type and the crop NDVI profiles were then visually assessed and compared to the profiles of 2001 to ascertain if each crop's unique phenological pattern was consistent between the two years. Second, separability within each crop class in the time-series NDVI data between 2001 and 2005 was investigated numerically using the Jeffries-Matusita (JM) distance statistic. The results seem to suggest that time-series NDVI response curves for crops over a growing period for one year of valid ground reference data may be useful for mapping crops for a different year when minor temporal shifts in the NDVI values (resulting from inter-annual climate variations or changes in agricultural management practices) are taken into account.  相似文献   

18.
The rice land is linked to the climate change due to its methane emission potential. The systems of growing rice and associated soil and crop management practices that have evolved are varied and complex. However, from the methane emission point of view, water regime is a crucial parameter. According to IPCC guidelines the rice ecosystem need to be categorized into four strata for methane emission study. The remote sensing based stratification map previously developed was used for in-situ weekly/monthly measurements of methane emission from the representative ecosystems, samples were collected and analysed using gas chromatography following the IPCC standards for three consecutive years; 2003, 2004 and 2005. This paper highlights the results of methane emission measurement and pattern from rice lands of India based on in-situ measurements. The CH4 emission pattern of irrigated crop in dry season showed a steady increase in the beginning which peaks during flowering stage, decreasing gradually thereafter. The results were consistent for different varieties and across the years. The emission pattern of irrigated wet season crop showed two peaks. The emission pattern also showed the influence of crop variety as well as year (of observation). The mean emission coefficient derived from all categories and all samples (n = 471) weighted for the Indian rice crop was 74.05 + 43.28 kg/ha.  相似文献   

19.
Abstract

This study aims to develop an approach to characterize cropland drought conditions in El Salvador, Central America. The data were processed for 2016–2017 through three main steps: (1) reconstructing MODIS land-surface temperature (LST), (2) Landsat-MODIS data fusion and (3) drought delineation using the temperature vegetation dryness index (TVDI). The results of LST reconstruction using the random forests (RF) indicated the median RMSE value of 0.5?°C. The fusion results achieved from the STARFM compared with the reference Landsat data revealed close agreement with the correlation coefficient (r) values higher than 0.84. The TVDI results verified with that from the reference Landsat data indicated r values of 0.85 and 0.75 for 2016 and 2017, respectively. The larger very dry area was observed for the 2016 primera season due to prolonged droughts. Approximately 11.5% and 10.7% of croplands were, respectively, associated with very dry moisture condition in the 2016 and 2017 primera seasons.  相似文献   

20.
The Upper Tons River Basin of North India has been selected for prioritisation of sub-watersheds (SW) based on morphometric parameters with respect to groundwater derived from topographic sheets and CARTOSAT data. There are 10 SW have been delineated in the region, high stream frequency (Fs) values of SW (1–5) and SW-9 indicated the occurrence of steep slopes, less permeable rocks, greater runoff, less infiltration possibility. Further, these regions have been predicted as poor groundwater potentialities. SW-2 has been identified as poorest groundwater potential zone, whereas SW-4 and SW (6–8) regions possess good permeable bed rocks. The Drainage density (Dd) map demonstrated that the middle south-west region possesses higher Dd whereas northeastern regions contain lower Dd. Further, the areal parameters indicate elongated shape of the basin, hilly region has moderate to steeper ground slope. The outcomes of work have potential to manage groundwater and to ameliorate the flash flood and droughts.  相似文献   

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

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