首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 609 毫秒
1.
Land use and land cover change are of prime concern due to their impacts on CO2 emissions, climate change and ecological services. New global land cover products at 300 m resolution from the European Space Agency (ESA) Climate Change Initiative Land Cover (CCI LC) project for epochs centered around 2000, 2005 and 2010 were analyzed to investigate forest area change and land cover transitions. Plant functional types (PFTs) fractions were derived from these land cover products according to a conversion table. The gross global forest loss between 2000 and 2010 is 172,171 km2, accounting for 0.6% of the global forest area in year 2000. The forest changes are mainly distributed in tropical areas such as Brazil and Indonesia. Forest gains were only observed between 2005 and 2010 with a global area of 9844 km2, mostly from crops in Southeast Asia and South America. The predominant PFT transition is deforestation from forest to crop, accounting for four-fifths of the total increase of cropland area between 2000 and 2010. The transitions from forest to bare soil, shrub, and grass also contributed strongly to the total areal change in PFTs. Different PFT transition matrices and composition patterns were found in different regions. The highest fractions of forest to bare soil transitions were found in the United States and Canada, reflecting forest management practices. Most of the degradation from grassland and shrubland to bare soil occurred in boreal regions. The areal percentage of forest loss and land cover transitions generally decreased from 2000–2005 to 2005–2010. Different data sources and uncertainty in the conversion factors (converting from original LC classes to PFTs) contribute to the discrepancy in the values of change in absolute forest area.  相似文献   

2.
Tongyu County in Northeast China is highly prone to land degradation due to its fragile physical settings characterized by a flat topography, a semi-arid climate, and a shallow groundwater table. This study aims to determine the causes of land degradation through detecting the long-term trend of land cover changes. Degraded lands were mapped from satellite images recorded in 1992 and 2002. These land cover maps revealed that the area subject to land degradation in the form of soil salinization, waterlogging and desertification increased from 2400 to 4214 km2, in sharp contrast to most severely degraded land that decreased by 122.5 km2. Newly degraded land stems from productive farmland (263 km2), harvested farmland (551 km2), and grassland (468 km2). Therefore, the worsened degradation situation is attributed to excessive reclamation of grassland for farming, over cultivation, overgrazing, and deforestation. Mechanical, biological, ecological and engineering means should be adopted to rehabilitate the degraded land.  相似文献   

3.
Istanbul is the largest city in Turkey with an area of around 5750 km2 and a population of around 10.8 M (2000). In 1980, the population was only around 4.7 M and so has more than doubled in only 2 decades. In 2000, around 65% of the population were living on the European side of the city with its large industrial/commercial and trade centres. The population is increasing as a result of both births exceeding deaths and mass immigration. Consequently, planned and unplanned housing are increasing while green areas are decreasing in area. Monitoring urban growth will enable the Municipality of Istanbul to better manage this complex urban area.  相似文献   

4.
A tree survey and an analysis of high resolution satellite data were performed to characterise the woody vegetation within a 10 × 10 km2 area around a site located close to the town of Dahra in the semi-arid northern part of Senegal. The surveyed parameters were tree species, height, tree crown radius, and diameter at breast height (DBH), for which allometric models were determined. An object-based classification method was used to determine tree crown cover (TCC) from Quickbird data. The average TCC from the tree survey and the respective TCC from remote sensing were both about 3.0%. For areas beyond the surveyed areas TCC varied between 3.0% and 4.5%. Furthermore, an empirical correction factor for tree clumping was obtained, which considerably improved the estimated number of trees and the estimated average tree crown area and radius. An allometric model linking TCC to tree stem crosssectional area (CSA) was developed, which allows to estimate tree biomass from remote sensing. The allometric models for the three main tree species found performed well and had r2-values of about 0.7–0.8.  相似文献   

5.
This study attempts to identify and forecast future land cover (LC) by using the Land Transformation Model (LTM), which considers pixel changes in the past and makes predictions using influential spatial features. LTM applies the Artificial Neural Networks algorithm) in conducting the analysis. In line with these objectives, two satellite images (Spot 5 acquired in 2004 and 2010) were classified using the Maximum Likelihood method for the change detection analysis. Consequently, LC maps from 2004 to 2010 with six classes (forest, agriculture, oil palm cultivations, open area, urban, and water bodies) were generated from the test area. A prediction was made on the actual soil erosion and the soil erosion rate using the Universal Soil Loss Equation (USLE) combined with remote sensing and GIS in the Semenyih watershed for 2004 and 2010 and projected to 2016. Actual and potential soil erosion maps from 2004 to 2010 and projected to 2016 were eventually generated. The results of the LC change detections indicated that three major changes were predicted from 2004 to 2016 (a period of 12 years): (1) forest cover and open area significantly decreased at rates of almost 30 and 8 km2, respectively; (2) cultivated land and oil palm have shown an increment in sizes at rates of 25.02 and 5.77 km2, respectively; and, (3) settlement and Urbanization has intensified also by almost 5 km2. Soil erosion risk analysis results also showed that the Semenyih basin exhibited an average annual soil erosion between 143.35 ton ha?1 year?1 in 2004 and 151 in 2010, followed by the expected 162.24 ton ha?1 year?1. These results indicated that Semenyih is prone to water erosion by 2016. The wide range of erosion classes were estimated at a very low level (0–1 t/ha/year) and mainly located on steep lands and forest areas. This study has shown that using both LTM and USLE in combination with remote sensing and GIS is a suitable method for forecasting LC and accurately measuring the amount of soil losses in the future.  相似文献   

6.
The studies on forest cover change can reveal the status of forests and facilitate for its conservation planning. Idukki is the largest district in the state of Kerala having a total geographical area of 5019 km2. The objectives of the present study are to map forest cover in Idukki district using multi-temporal remote sensing data (1975, 1990, 2001 and 2012) and topographical maps (1925), to analyze the trends in deforestation and land use changes. Overall statistics for the period of 1925 indicate that about 4675.7 km2 (93.2 %) of the landscape was under forest. The forest cover in 2012 was estimated as 2613.4 km2 (52.1 %). Recently, due to the implementation of policies and protection efforts, the rate of deforestation was greatly reduced. The commencement of hydroelectric projects during 1925–1990 responsible for an increase of area under water bodies by inundating other land uses. The long term analysis shows agricultural area been decreasing and commercial plantations been increasing in the district. There has been a significant increase in the area of plantations from 1236.2 km2 (1975) to 1317.3 km2 (2012).  相似文献   

7.
Sagebrush ecosystems of the western US provide important habitat for several ungulate and vertebrate species. As a consequence of energy development, these ecosystems in Wyoming have been subjected to a variety of anthropogenic disturbances. Land managers require methodology that will allow them to consistently catalog sagebrush ecosystems and evaluate potential impact of proposed anthropogenic activities. This study addresses the utility of remotely sensed and ancillary geospatial data to estimate sagebrush cover using ordinal logistic regression. We demonstrate statistically significant prediction of ordinal sagebrush cover categories using spectral (χ2 = 113; p < 0.0001) and transformed indices (χ2 = 117; p < 0.0001). Both Landsat spectral bands (c-value = 0.88) and transformed indices (c-value = 0.89) can distinguish sites with closed, moderate and open cover sagebrush cover categories from no cover. The techniques described in this study can be used for estimating categories of sagebrush cover in arid ecosystems.  相似文献   

8.
Sikkim is a small, mountainous, Indian state (7,096 km2) located in the eastern Himalayan region. Though a global biodiversity hotspot, it has been relatively less studied. A detailed forest type, density and change dynamics study was undertaken, using SATELLITE remote sensing data and intensive field verification. The landscape was found to be dominated by alpine and nival ecosystems, with a large portion above the tree line, considerable snow cover, and a sizeable area under forest cover (72%, 5,094 km2). A total of 18 landscape components including 14 vegetation classes were delineated, with the major ones being oak forest, alpine meadow, alpine scrub, conifer forest and alder-cardamom agro-forestry. Of the 3,154 km2 of forests below the tree line, 40% were found to be dense (>40% tree canopy cover). A sizeable portion of the non dense forests below the tree line was contributed by the degradation of oak forests, which was confirmed by change detection analysis. However on a positive front over the past decade, ban on grazing and felling of trees in forests has been implemented. In order to expand the extent of dense forests, further efforts are needed for the restoration of oak forests such as fire protection, providing alternatives to firewood use, promotion of alder-cardamom agro-forestry in the private lands and protection of the small-sized, fragmented forest patches in the subtropical belt.  相似文献   

9.
Soil erosion rates in alpine regions are related to high spatial variability complicating assessment of risk and damages. A crucial parameter triggering soil erosion that can be derived from satellite imagery is fractional vegetation cover (FVC). The objective of this study is to assess the applicability of normalized differenced vegetation index (NDVI), linear spectral unmixing (LSU) and mixture tuned matched filtering (MTMF) in estimating abundance of vegetation cover in alpine terrain. To account for the small scale heterogeneity of the alpine landscape we used high resolved multispectral QuickBird imagery (pixel resolution = 2.4 m) of a site in the Urseren Valley, Central Swiss Alps (67 km2). A supervised land-cover classification was applied (total accuracy 93.3%) prior to the analysis in order to stratify the image. The regression between ground truth FVC assessment and NDVI as well as MTMF-derived vegetation abundance was significant (r2 = 0.64, r2 = 0.71, respectively). Best results were achieved for LSU (r2 = 0.85). For both spectral unmixing approaches failed to estimate bare soil abundance (r2 = 0.39 for LSU, r2 = 0.28 for MTMF) due to the high spectral variability of bare soil at the study site and the low spectral resolution of the QuickBird imagery. The LSU-derived FVC map successfully identified erosion features (e.g. landslides) and areas prone to soil erosion. FVC represents an important but often neglected parameter for soil erosion risk assessment in alpine grasslands.  相似文献   

10.
The Western Ghats constitute one of the three biodiversity hot spots in India, which is under constant threat from various quarters. Among the several anthropogenic causes, fire is one of the important anthropogenic factor, which plays a pivotal role in vegetation succession and ecosystem processes. It is very important to understand the ecological changes due to fire and other anthropogenic factors for conservation and management of biodiversity. Because of its synoptic, multi-spectral and multi-temporal nature remote sensing data can be a good source for forest fire monitoring. In the present study, an effort has been made to monitor the burnt areas using March 2000 and 2004 IRS LISS — III data. The study revealed that an area of 2.15 km2 and 4.46 km2 was affected by fire in 2000 and 2004 respectively. Repeated drought, followed by mass flowering and dying of bamboo accelerated the spread of fire from ground to canopy in areas with high bamboo density.  相似文献   

11.
Defoliation is a key parameter of forest health and is associated with reduced productivity and tree mortality. Assessing the health of forests requires regular observations over large areas. Satellite remote sensing provides a cost-effective alternative to traditional ground-based assessment of forest health, but assessing defoliation can be difficult due to mixed pixels where vegetation cover is low or fragmented. In this study we apply a novel spectral unmixing technique, referred to as weighted Multiple Endmember Spectral Mixture Analysis (wMESMA), to Landsat 5-TM and EO-1 Hyperion data acquired over a Eucalyptus globulus (Labill.) plantation in southern Australia. This technique combines an iterative mixture analysis cycle allowing endmembers to vary on a per pixel basis (MESMA) and a weighting algorithm that prioritizes wavebands based on their robustness against endmember variability. Spectral mixture analysis provides an estimate of the physically interpretable canopy cover, which is not necessarily correlated with defoliation in mixed-aged plantations due to natural variation in canopy cover as stands age. There is considerable variability in the degree of defoliation as well as in stand age among sites and in this study we found that results were significantly improved by the inclusion of an age correction algorithm for both the multi-spectral (R2no age correction = 0.55 vs R2age correction = 0.73 for Landsat) and hyperspectral (R2no age correction = 0.12 vs R2age correction = 0.50 for Hyperion) image data. The improved accuracy obtained from Landsat compared to the Hyperion data illustrates the potential of applying SMA techniques for analysis of multi-spectral datasets such as MODIS and SPOT-VEGETATION.  相似文献   

12.
The Moderate Resolution Imaging Spectrometer (MODIS) Normalized Difference Vegetation Index (NDVI) 16-day composite data product (MOD12Q) was used to develop annual cropland and crop-specific map products (corn, soybeans, and wheat) for the Laurentian Great Lakes Basin (GLB). The crop area distributions and changes in crop rotations were characterized by comparing annual crop map products for 2005, 2006, and 2007. The total acreages for corn and soybeans were relatively balanced for calendar years 2005 (31,462 km2 and 31,283 km2, respectively) and 2006 (30,766 km2 and 30,972 km2, respectively). Conversely, corn acreage increased approximately 21% from 2006 to 2007, while soybean and wheat acreage decreased approximately 9% and 21%, respectively. Two-year crop rotational change analyses were conducted for the 2005–2006 and 2006–2007 time periods. The large increase in corn acreages for 2007 introduced crop rotation changes across the GLB. Compared to 2005–2006, crop rotation patterns for 2006–2007 resulted in increased corn–corn, soybean–corn, and wheat–corn rotations. The increased corn acreages could have potential negative impacts on nutrient loadings, pesticide exposures, and sediment-mediated habitat degradation. Increased in US corn acreages in 2007 were related to new biofuel mandates, while Canadian increases were attributed to higher world-wide corn prices. Additional study is needed to determine the potential impacts of increases in corn-based ethanol agricultural production on watershed ecosystems and receiving waters.  相似文献   

13.
The focus of soil erosion research in the Alps has been in two categories: (i) on-site measurements, which are rather small scale point measurements on selected plots often constrained to irrigation experiments or (ii) off-site quantification of sediment delivery at the outlet of the catchment. Results of both categories pointed towards the importance of an intact vegetation cover to prevent soil loss. With the recent availability of high-resolution satellites such as IKONOS and QuickBird options for detecting and monitoring vegetation parameters in heterogeneous terrain have increased. The aim of this study is to evaluate the usefulness of QuickBird derived vegetation parameters in soil erosion models for alpine sites by comparison to Cesium-137 (Cs-137) derived soil erosion estimates. The study site (67 km2) is located in the Central Swiss Alps (Urseren Valley) and is characterised by scarce forest cover and strong anthropogenic influences due to grassland farming for centuries. A fractional vegetation cover (FVC) map for grassland and detailed land-cover maps are available from linear spectral unmixing and supervised classification of QuickBird imagery. The maps were introduced to the Pan-European Soil Erosion Risk Assessment (PESERA) model as well as to the Universal Soil Loss Equation (USLE). Regarding the latter model, the FVC was indirectly incorporated by adapting the C factor. Both models show an increase in absolute soil erosion values when FVC is considered. In contrast to USLE and the Cs-137 soil erosion rates, PESERA estimates are low. For the USLE model also the spatial patterns improved and showed “hotspots” of high erosion of up to 16 t ha−1 a−1. In conclusion field measurements of Cs-137 confirmed the improvement of soil erosion estimates using the satellite-derived vegetation data.  相似文献   

14.
We estimated urbanization rates (2001–2006) in the Gulf of Mexico region using the National Land Cover Database (NLCD) 2001 and 2006 impervious surface products. An improved method was used to update the NLCD impervious surface product in 2006 and associated land cover transition between 2001 and 2006. Our estimation reveals that impervious surface increased 416 km2 with a growth rate of 5.8% between 2001 and 2006. Approximately 1110.1 km2 of non-urban lands were converted into urban land, resulting in a 3.2% increase in the region. Hay/pasture, woody wetland, and evergreen forest represented the three most common land cover classes that transitioned to urban. Among these land cover transitions, more than 50% of the urbanization occurred within 50 km of the coast. Our analysis shows that the close-to-coast land cover transition trend, especially within 10 km off the coast, potentially imposes substantial long-term impacts on regional landscape and ecological conditions.  相似文献   

15.
Estimation of forest structural parameters by field-based data collection methods is both expensive and time consuming. Satellite remote sensing is a low-cost alternative in modeling and mapping structural parameters in large forest areas. The current study investigates the potential of using WordView-2 multispectral satellite imagery for predicting forest structural parameters in a dryland plantation forest in Israel. The relationships between image texture features and the several structural parameters such as Number of Trees (NT), Basal Area (BA), Stem Volume (SV), Clark-Evans Index (CEI), Diameter Differentiation Index (DDI), Contagion Index (CI), Gini Coefficient (GC), and Standard Deviation of Diameters at Breast Heights (SDDBH) were examined using correlation analyses. These variables were obtained from 30 m × 30 m square-shaped plots. The Standard Deviation of Gray Levels (SDGL) as a first order texture feature and the second order texture variables based on Gray Level Co-occurrence Matrix (GLCM) were calculated for the pixels that corresponds to field plots. The results of the correlation analysis indicate that the forest structural parameters are significantly correlated with the image texture features. The highest correlation coefficients were calculated for the relationships between the SDDBH and the contrast of red band (r = 0.75, p < 0.01), the BA and the entropy of blue band (r = 0.73, p < 0.01), and the GC and the contrast of blue band (r = 0.71, p < 0.01). Each forest structural parameter was modeled as a function of texture measures derived from the satellite image using stepwise multi linear regression analyses. The determination coefficient (R2) and root mean square error (RMSE) values of the best fitting models, respectively, are 0.38 and 109.56 ha−1 for the NT; 0.54 and 1.79 m2 ha−1 for the BA; 0.42 and 27.18 m3 ha−1 for the SV; 0.23 and 0.16 for the CEI; 0.32 and 0.05 for the DDI; 0.25 and 0.06 for the CI; 0.50 and 0.05 for the GC; and 0.67 and 0.70 for the SDDBH. The leave-one-out cross-validation technique was applied for validation of the best-fitted models (R2 > 0.50). In conclusion, cross-validated statistics confirmed that the structural parameters including the BA, SDDBH, and GC can be predicted and mapped with a reasonable accuracy using the texture features extracted from the spectral bands of WorldView-2 image.  相似文献   

16.
This paper evaluates the potential of a terrestrial laser scanner (TLS) to characterize forest canopy fuel characteristics at plot level. Several canopy properties, namely canopy height, canopy cover, canopy base height and fuel strata gap were estimated. Different approaches were tested to avoid the effect of canopy shadowing on canopy height estimation caused by deployment of the TLS below the canopy. Estimation of canopy height using a grid approach provided a coefficient of determination of R2 = 0.81 and an RMSE of 2.47 m. A similar RMSE was obtained using the 99th percentile of the height distribution of the highest points, representing the 1% of the data, although the coefficient of determination was lower (R2 = 0.70). Canopy cover (CC) was estimated as a function of the occupied cells of a grid superimposed upon the TLS point clouds. It was found that CC estimates were dependent on the cell size selected, with 3 cm being the optimum resolution for this study. The effect of the zenith view angle on CC estimates was also analyzed. A simple method was developed to estimate canopy base height from the vegetation vertical profiles derived from an occupied/non-occupied voxels approach. Canopy base height was estimated with an RMSE of 3.09 m and an R2 = 0.86. Terrestrial laser scanning also provides a unique opportunity to estimate the fuel strata gap (FSG), which has not been previously derived from remotely sensed data. The FSG was also derived from the vegetation vertical profile with an RMSE of 1.53 m and an R2 = 0.87.  相似文献   

17.
Countries like Iran, which are geographically situated in a rather arid and warm regions, will suffer more from global warming than countries located in humid and semi-humid regions. In such environments, analyzing the variations of mountain glaciers can reveal several aspects of climate change patterns more efficiently in comparison to the other geo-indicators. The present study reports some evidence of changes for Alamkouh glacier between 1955 and 2010 based on several mediums to high-resolution satellite images. Considering that most part of the Alamkouh glacier is covered by debris and delineating its actual area is not possible, planimetric change analysis was restricted to the clean-ice regions. The object-oriented classification approach was used to estimate the clean ice areas. This technique takes into account the shapes of the features along with their spectral patterns. Results revealed that clean ice regions of Alamkouh glacier shrank since 1955 with an overall area reduction of about 59 %. Although the general observed trend is a clean ice area reduction, some advancement was detected over the period from 2000 to 2010. During 1992–2000, the maximum reduction in the clean ice area was observed (0.08 km2.a?1). However, clean ice area of the case study has partially increased about 0.028 km2.a?1 from 2000 to 2010. Supra-glacial lake change analysis illustrated that at the surface of the glacier, lakes have been enlarged remarkably in the past 55 years (about 4.75 times greater). In addition, clean ice area and the surface area of supra-glacial lakes oscillated in compliance with each other. The findings revealed that the maximum expansion of supra-glacial lake occurred during 1992–2000, which demonstrate the glacier maximum reduction during this period. This shrinkage in the Alamkouh glacier caused an extensive glacial lake outburst flood in Jun 2011. The results of this study agree with documented changes in other mountain glaciers located in arid and semi-arid environments and they also confirm the application of mountain glaciers in climate variations monitoring over such regions.  相似文献   

18.
Estimating forest structural attributes using multispectral remote sensing is challenging because of the saturation of multispectral indices at high canopy cover. The objective of this study was to assess the utility of hyperspectral data in estimating and mapping forest structural parameters including mean diameter-at-breast height (DBH), mean tree height and tree density of a closed canopy beech forest (Fagus sylvatica L.). Airborne HyMap images and data on forest structural attributes were collected from the Majella National Park, Italy in July 2004. The predictive performances of vegetation indices (VI) derived from all possible two-band combinations (VI(i,j) = (Ri − Rj)/(Ri + Rj), where Ri and Rj = reflectance in any two bands) were evaluated using calibration (n = 33) and test (n = 20) data sets. The potential of partial least squares (PLS) regression, a multivariate technique involving several bands was also assessed. New VIs based on the contrast between reflectance in the red-edge shoulder (756–820 nm) and the water absorption feature centred at 1200 nm (1172–1320 nm) were found to show higher correlations with the forest structural parameters than standard VIs derived from NIR and visible reflectance (i.e. the normalised difference vegetation index, NDVI). PLS regression showed a slight improvement in estimating the beech forest structural attributes (prediction errors of 27.6%, 32.6% and 46.4% for mean DBH, height and tree density, respectively) compared to VIs using linear regression models (prediction errors of 27.8%, 35.8% and 48.3% for mean DBH, height and tree density, respectively). Mean DBH was the best predicted variable among the stand parameters (calibration R2 = 0.62 for an exponential model fit and standard error of prediction = 5.12 cm, i.e. 25% of the mean). The predicted map of mean DBH revealed high heterogeneity in the beech forest structure in the study area. The spatial variability of mean DBH occurs at less than 450 m. The DBH map could be useful to forest management in many ways, e.g. thinning of coppice to promote diameter growth, to assess the effects of management on forest structure or to detect changes in the forest structure caused by anthropogenic and natural factors.  相似文献   

19.
A suitable index is proposed to evaluate the natural short–medium-term recovery capability of vegetation in burnt areas. The study area covers 2450 km2 in western Tuscany (Province of Pisa, Italy). This region is characterized by a typical Mediterranean climate and is subject to fire damage during the dry summer season. Damage is mitigated where a natural rapid regrowth of vegetation prevents soil erosion, supporting the return to a natural pre-fire state.  相似文献   

20.
Floodplain wetlands in the China side of the Amur River Basin (CARB) undergone consistent decreases because of both natural and anthropogenic drivers. Monitoring floodplain wetlands dynamics and conversions over long-time periods is thus fundamental to sustainable management and protection. Due to complexity and heterogeneity of floodplain environments, however, it is difficult to map wetlands accurately over a large area as the CARB. To address this issue, we developed a novel and robust classification approach integrating image compositing algorithm, objected-based image analysis, and hierarchical random forest classification, named COHRF, to delineate floodplain wetlands and surrounding land covers. Based on the COHRF classification approach, 4622 Landsat images were applied to produce a 30-m resolution dataset characterizing dynamics and conversions of floodplain wetlands in the CARB during 1990–2018. Results show that (1) all floodplain land cover maps in 1990, 2000, 2010, and 2018 had high mapping accuracies (ranging from 90 %±0.001–97%±0.005), suggesting that COHRF is a robust classification approach; (2) CARB experienced an approximately 25 % net loss of floodplain wetlands with an area declined from 8867 km2 to 6630 km2 during 1990–2018; (3) the lost floodplain wetlands were mostly converted into croplands, while, there were 111 km2 and 256 km2 of wetlands rehabilitated from croplands during periods of 2000–2010 and 2010–2018, respectively. To our knowledge, this study is the first attempt that focus on delineating floodplain wetlands at a large-scale and produce the first 30-m spatial resolution dataset demonstrating long-term dynamics of floodplain wetlands in the CARB. The COHRF classification approach could be used to classify other ecosystems readily and robustly. The resultant dataset will contribute to sustainable use and conservation of wetlands in the Amur River Basin and provide essential information for related researches.  相似文献   

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

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