首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 109 毫秒
1.
Canada is dominated by forested ecosystems which are subject to various inventory and management practices, with more northern boreal forests subject to neither. Our objectives were to measure the capacity of temporal trajectory metrics for estimating selected forest attributes in a northern Canadian boreal forest context using Landsat imagery and investigate the importance of different types of temporal trajectory metrics. Results indicated that Wetness was the best Tasseled Cap (TC) component for aboveground biomass estimation (R2 = 50%, RMSE% = 56%), and the combination of simple and complex metrics from all TC components produced the highest R2 (62%) and lowest RMSE% (49%). Using a similar combination of variables, other forest attributes were estimated equally reliably with lower RMSE% values. The most important temporal trajectory metrics were simple and described TC component values at each point of change in the temporal trajectory, however the most important variables overall were environmental variables.  相似文献   

2.
Airborne LiDAR techniques can provide accurate measurements of tree height, from which estimates of stem volume and forest woody biomass can be obtained. These techniques, however, are still expensive to apply repeatedly over large areas. The current paper presents a methodology which first transforms mean stand heights obtained from LiDAR over small strips into relevant stem volume estimates. These are then extended over an entire forest by applying two estimation methods (k-NN and locally calibrated regression) to Landsat ETM+ images. The methodology is tested over a coastal area covered by pine forest in the Regional Park of San Rossore (Central Italy). The results are evaluated by comparison with the ground stem volumes of a recent forest inventory, taking into consideration the effect of stand size. In general, the accuracies of two estimation methods are dependent on the size of the forest stands and are satisfactory only when considering stands larger than 5-10 ha. The outputs of the parametric regression procedure are slightly more stable than those of k-NN and more faithfully reproduce the spatial patterns of the ground data.  相似文献   

3.
Coffee is a commodity of international trade significance, and its value chain can benefit from age-specific thematic maps. This study aimed to assess the potential of Landsat 8 OLI to develop these maps. Using field-collected samples with the random forest classifier, splitting coffee into three age classes (Scheme A) was compared with running the classification with one compound coffee class (Scheme B). Higher overall classification accuracy was obtained in Scheme B (90.3% for OLI and 86.8% for ETM+) than in Scheme A (86.2% for OLI and 81.0% for ETM+). The NIR band of OLI was the most important band in intra-class discrimination of coffee. Landsat 8 OLI mapped area closely matched farm records (R2?=?0.88) compared to that of Landsat 7 ETM+ (R2?=?0.78). It was concluded that Landsat 8 OLI data can be used to produce age-specific thematic maps in coffee production areas although disaggregating coffee classes reduces overall accuracy.  相似文献   

4.
Locally computed statistics of image texture and a case-based reasoning (CBR) system were evaluated for mapping of forest attributes. Cluster analysis was preferred to regression models, as a pre-selection method of features. The best stand-based accuracy using satellite sensor images was 74.64 m−3 ha−1 (36%) RMSE for stand volume, 1.98 m−3 ha−1 a−1 (49%) for annual increase in stand volume, where κ = 0.23 for stand growth classes and κ = 0.41 for dominant tree species in stands. The top pixel-based accuracy using orthophotos was 76.54 m−3 ha−1 (41%) RMSE for stand volume, 1.87 m−3 ha−1 a−1 (44%) for annual increase in stand volume, where κ = 0.24 for stand growth classes and κ = 0.38 for dominant tree species in stands. Mean saturation in 30 m radius was the most useful feature when orthophotos were used, and standard deviation of Landsat ETM 6.2 values in 80 m radius was the best when satellite sensor images were used. The most valuable feature components (radii, channels and local statistics) for orthophotos were: 30 m kernel radius, lightness and the mean of pixel values; for satellite sensor images: 80 m kernel radius, near-infrared channel (ETM 4) and the mean of pixel values. Locally computed statistics.  相似文献   

5.
Tree species composition of forest stand is an important indicator of forest inventory attributes for assessing ecosystem health, understanding successional processes, and digitally displaying forest biodiversity. In this study, we acquired high spatial resolution multispectral and RGB imagery over a subtropical natural forest in southwest China using a fixed-wing UAV system. Digital aerial photogrammetric (DAP) technique was used to generate multi-spectral and RGB derived point clouds, upon which individual tree crown (ITC) delineation algorithms and a machine learning classifier were used to identify dominant tree species. To do so, the structure-from-motion method was used to generate RGB imagery-based DAP point clouds. Then, three ITC delineation algorithms (i.e., point cloud segmentation (PCS), image-based multiresolution segmentation (IMRS), and advanced multiresolution segmentation (AMRS)) were used and assessed for ITC detection. Finally, tree-level metrics (i.e., multispectral, texture and point cloud metrics) were used as metrics in the random forest classifier used to classify eight dominant tree species. Results indicated that the accuracy of the AMRS ITC segmentation was highest (F1-score = 82.5 %), followed by the segmentation using PCS (F1-score = 79.6 %), the IMRS exhibited the lowest accuracy (F1-score = 78.6 %); forest types classification (coniferous and deciduous) had a higher accuracy than the classification of all eight tree species, and the combination of spectral, texture and structural metrics had the highest classification accuracy (overall accuracy = 80.20 %). In the classification of both eight tree species and two forest types, the classification accuracies were lowest when only using spectral metrics, indicated that the texture metrics and point cloud structural metrics had a positive impact on the classification (the overall accuracy and kappa accuracy increased by 1.49–4.46 % and 2.86–6.84 %, respectively).  相似文献   

6.
  相似文献   

7.
Present study attempts to utilise the spatial information from IRS-IC LISS-III satellite data, through, tie use of Semivartogram apprqatfu for analyzing the overstorey. diversity In forest stands, Semivariance has been calculated for all the pixels in I5x15 pixel window for each of the forest stand using digital numbers from near infrared channel Variation in semivariance has been correlated to stand density and species composition. Analysis of averaged semivariance for overall diversity of forest stands suggested, decreasing semi variance in the order, mixed dry (feciduous forests, mixed scrub forests followed by dry deciduous forest Semivariance has been found to be highly Correlated to tree density (R2 = 0.96) suggesting semivariograms as one of the measures for studies on tree density, canapy cover and diversity patterns.  相似文献   

8.
机载激光雷达及高光谱的森林乔木物种多样性遥感监测   总被引：1，自引：0，他引：1  
利用机载LiDAR和高光谱数据并结合37个地面调查样本数据,基于结构差异与光谱变异理论,通过相关分析法分别筛选了3个最优林冠结构参数和6个最优光谱指数,在单木尺度上利用自适应C均值模糊聚类算法,在神农架国家自然保护区开展森林乔木物种多样性监测,实现了森林乔木物种多样性的区域成图。研究结果表明,(1)基于结合形态学冠层控制的分水岭算法可以获得较高精度的单木分割结果(R~2=0.88,RMSE=13.17,P0.001);(2)基于LiDAR数据提取的9个结构参数中,95%百分位高度、冠层盖度和植被穿透率为最优结构参数,与Shannon-Wiener指数的相关性达到R~2=0.39—0.42(P0.01);(3)基于机载高光谱数据筛选的16个常用的植被指数中,CRI、OSAVI、Narrow band NDVI、SR、Vogelmann index1、PRI与Shannon-Wiener指数的相关性最高(R~2=0.37—0.45,P0.01);(4)在研究区,利用以30 m×30 m为窗口的自适应模糊C均值聚类算法可预测的最大森林乔木物种数为20,物种丰富度的预测精度为R~2=0.69,RMSE=3.11,Shannon-Wiener指数的预测精度为R~2=0.70,RMSE=0.32。该研究在亚热带森林开展乔木物种多样性监测,是在区域尺度上进行物种多样性成图的重要实践,可有效补充森林生物多样性本底数据的调查手段,有助于实现生物多样性的长期动态监测及科学分析森林物种多样性的现状和变化趋势。  相似文献   

9.
Assessment of the susceptibility of forests to mountain pine beetle (Dendroctonus ponderosae Hopkins) infestation is based upon an understanding of the characteristics that predispose the stands to attack. These assessments are typically derived from conventional forest inventory data; however, this information often represents only managed forest areas. It does not cover areas such as forest parks or conservation regions and is often not regularly updated resulting in an inability to assess forest susceptibility. To address these shortcomings, we demonstrate how a geometric optical model (GOM) can be applied to Landsat-5 Thematic Mapper (TM) imagery (30 m spatial resolution) to estimate stand-level susceptibility to mountain pine beetle attack. Spectral mixture analysis was used to determine the proportion of sunlit canopy and background, and shadow of each Landsat pixel enabling per pixel estimates of attributes required for model inversion. Stand structural attributes were then derived from inversion of the geometric optical model and used as basis for susceptibility mapping. Mean stand density estimated by the geometric optical model was 2753 (standard deviation ± 308) stems per hectare and mean horizontal crown radius was 2.09 (standard deviation ± 0.11) metres. When compared to equivalent forest inventory attributes, model predictions of stems per hectare and crown radius were shown to be reasonably estimated using a Kruskal–Wallis ANOVA (p < 0.001). These predictions were then used to create a large area map that provided an assessment of the forest area susceptible to mountain pine beetle damage.  相似文献   

10.
Abstract

A linear regression‐based model for mapping forest age using Landsat Thematic Mapper data is evaluated in the lodgepole pine forests of Yellowstone National Park. Regression models predicting age (R2=0.62) and a logarithmic transformation of age (R2 = 0.90) used a combination of visible, near‐infrared, and middle‐infrared TM bands. Forest age maps produced using the regression method match broad‐scale patterns of forest age within the Yellowstone Central Plateau study area. Per‐pixel estimates of forest age, however, may depart substantially from actual forest age, particularly for older stands, and the maps are most appropriate for depicting regional patterns of forest age.  相似文献   

11.
Understanding forest biomass dynamics is crucial for carbon and environmental monitoring, especially in the context of climate change. In this study, we propose a robust approach for monitoring aboveground forest biomass (AGB) dynamics by combining Landsat time-series with single-date inventory data. We developed a Random Forest (RF) based kNN model to produce annual maps of AGB from 1988 to 2017 over 7.2 million ha of forests in Victoria, Australia. The model was internally evaluated using a bootstrapping technique. Predictions of AGB and its change were then independently evaluated using multi-temporal Lidar data (2008 and 2016). To understand how natural and anthropogenic processes impact forest AGB, we analysed trends in relation to the history of disturbance and recovery. Specifically, change metrics (e.g., AGB loss and gain, Years to Recovery - Y2R) were calculated at the pixel level to characterise the patterns of AGB change resulting from forest dynamics. The imputation model achieved a RMSE value of 132.9 Mg ha−1 (RMSE% = 46.3%) and R2 value of 0.56. Independent assessments of prediction maps in 2008 and 2016 using Lidar-based AGB data achieved relatively high accuracies, with a RMSE of 108.6 Mg ha−1 and 135.9 Mg ha−1 for 2008 and 2016, respectively. Annual validations of AGB maps using un-changed, homogenous Lidar plots suggest that our model is transferable through time (RMSE ranging from 109.65 Mg ha−1 to 112.27 Mg ha−1 and RMSE% ranging from 25.38% to 25.99%). In addition, changes in AGB values associated with forest disturbance and recovery (decrease and increase, respectively) were captured by predicted maps. AGB change metrics indicate that AGB loss and Y2R varied across bioregions and were highly dependent on levels of disturbance severity (i.e., a greater loss and longer recovery time were associated with a higher severity disturbance). On average, high severity fire burnt from 200 Mg ha−1 to 550 Mg ha−1 of AGB and required up to 15 years to recover while clear-fell logging caused a reduction in 250 Mg ha−1 to 600 Mg ha−1 of AGB and required nearly 20 years to recover. In addition, AGB within un-disturbed forests showed statistically significant but monotonic trends, suggesting a mild gradual drop over time across most bioregions. Our methods are designed to support forest managers and researchers in developing forest monitoring systems, especially in developing regions, where only a single date forestry inventory exists.  相似文献   

12.
Recently there have been reports of forest regrowth occurring in different regions across the world. There is also a growing recognition of the potential beneficial impact that secondary forests may have on the global environment: providing crucial ecosystem services such as soil conservation, stabilization of hydrological cycles, carbon sequestration, and support for forest dependent communities. Consequently, there is a growing awareness of the need to recognize that landscapes are complex shifting mosaics wherein forest clearing and reforestation take place. In this study, the rates of reforestation, deforestation, forest regrowth and degradation were measured using multi-temporal Landsat images of Danjiangkou, China. Landsat data from 1990, 1999 and 2007 were (1) classified as dense forest, open forest and non-forest areas and (2) compared between years to identify forest cutting, regeneration and degradation. The results showed that there was a net gain of 29,315 ha of forest area (including dense and open forest) from 1990 to 2007, showing a clear trend of reforestation in the study area. Forest modification (degradation and regrowth) and change categories (deforestation and reforestation) occurred simultaneously during the observation time period. Socioeconomic data from public statistics and environmental attributes allowed the assessment of the socioeconomic factors and the environmental conditions that caused these changes using non-metric multidimensional scaling (NMDS). The research showed that the socioeconomic factors due to different policies were major driving forces of forest transition, whereas environmental attributes of the underlying landscape constrained forest cover changes. These findings have led to a better understanding of forest transition at a local scale in our study region. Comprehensive knowledge of these relationships may be useful to reconstruct past forest transitions and predict future changes, and may help to enhance sustainable management practices aimed at preserving essential ecological functions.  相似文献   

13.
曹林  徐婷  申鑫  佘光辉 《遥感学报》2016,20(4):665-678
以亚热带天然次生林为研究对象,借助一个条带的少量LiDAR点云数据和覆盖整个研究区的免费Landsat OLI多光谱数据,并结合地面实测数据,探索森林生物量低成本高精度制图方法。首先,提取了OLI和LiDAR特征变量,并与地上和地下生物量进行相关分析以筛选变量;然后,借助LiDAR数据覆盖区的样地和条带LiDAR数据构建"LiDAR生物量模型";再从LiDAR反演生物量的结果中进行采样,结合OLI特征变量构建"LiDAR-OLI模型";最后,与单独使用OLI多光谱数据建立的"OLI估算模型"结果进行比较,分析精度并验证新方法的效果。结果表明,"LiDAR-OLI模型"对地上和地下生物量的模型拟合效果较好且均优于"OLI模型",且其交叉验证的精度也较高并优于"OLI模型",从而证明了新方法的可靠性及有效性。本研究为主、被动遥感技术在中小尺度上协同反演森林参数提供了实验基础,也为基于全覆盖免费OLI多光谱数据及条带LiDAR数据的低成本森林生物量制图探索了技术路线。  相似文献   

14.
In remote sensing–based forest aboveground biomass (AGB) estimation research, data saturation in Landsat and radar data is well known, but how to reduce this problem for improving AGB estimation has not been fully examined. Different vegetation types have their own species composition and stand structure, thus they have different data saturation values in Landsat or radar data. Optical and radar data also have different characteristics in representing forest stand structures, thus effective use of their features may improve AGB estimation. This research examines the effects of Landsat Thematic Mapper (TM) and ALOS PALSAR L-band data and their integrations in forest AGB estimation of Zhejiang Province, China, and the roles of textural images from both datasets. The linear regression models of AGB were conducted by using (1) Landsat TM alone, (2) ALOS PALSAR data alone, (3) their combination as extra bands, and (4) their data fusion, based on non-stratification and stratification of vegetation types, respectively. The results show that (1) overall, Landsat TM data perform better than PALSAR data, but the latter can produce more accurate estimates for bamboo and shrub, and for forests with AGB values less than 60 Mg/ha; (2) the combination of TM and PALSAR data as extra bands can greatly improve AGB estimation performance, but their fusion using the modified high-pass filter resolution-merging technique cannot; (3) textures are indeed valuable in AGB estimation, especially for forests with complex stand structures such as mixed forests and pine forests with understories of broadleaf species; (4) stratification of vegetation types can improve AGB estimation performance; and (5) the results from the linear regression models are characterized by overestimation and underestimation for the smaller and larger AGB values, respectively, and thus, selecting non-linear models or non-parametric algorithms may be needed in future research.  相似文献   

15.
Remote sensing technologies are an ideal platform to examine the extent and impact of fire on the landscape. In this study we assess that capacity of the RapidEye constellation and Landsat (Thematic Mapper and Operational Land Imager to map fine-scale burn attributes for a small, low severity prescribed fire in a dry Western Canadian forest. Estimates of burn severity from field data were collated into a simple burn index and correlated with a selected suite of common spectral vegetation indices. Burn severity classes were then derived to map fire impacts and estimate consumed woody surface fuels (diameter ≥2.6 cm). All correlations between the simple burn index and vegetation indices produced significant results (p < 0.01), but varied substantially in their overall accuracy. Although the Landsat Soil Adjusted Vegetation Index provided the best regression fit (R2 = 0.56), results suggested that RapidEye provided much more spatially detailed estimates of tree damage (Soil Adjusted Vegetation Index, R2 = 0.51). Consumption estimates of woody surface fuels ranged from 3.38 ± 1.03 Mg ha−1 to 11.73 ± 1.84 Mg ha−1, across four derived severity classes with uncertainties likely a result of changing foliage moisture between the before and after fire images. While not containing spectral information in the short wave infrared, the spatial variability provided by the RapidEye imagery has potential for mapping and monitoring fine scale forest attributes, as well as the potential to resolve fire damage at the individual tree level.  相似文献   

16.

Background

A simulation model that relies on satellite observations of vegetation cover from the Landsat 7 sensor and from the Moderate Resolution Imaging Spectroradiometer (MODIS) was used to estimate net primary productivity (NPP) of forest stands at the Bartlett Experiment Forest (BEF) in the White Mountains of New Hampshire.

Results

Net primary production (NPP) predicted from the NASA-CASA model using 30-meter resolution Landsat inputs showed variations related to both vegetation cover type and elevational effects on mean air temperatures. Overall, the highest predicted NPP from the NASA-CASA model was for deciduous forest cover at low to mid-elevation locations over the landscape. Comparison of the model-predicted annual NPP to the plot-estimated values showed a significant correlation of R2 = 0.5. Stepwise addition of 30-meter resolution elevation data values explained no more than 20% of the residual variation in measured NPP patterns at BEF. Both the Landsat 7 and the 250-meter resolution MODIS derived mean annual NPP predictions for the BEF plot locations were within ± 2.5% of the mean of plot estimates for annual NPP.

Conclusion

Although MODIS imagery cannot capture the spatial details of NPP across the network of closely spaced plot locations as well as Landsat, the MODIS satellite data as inputs to the NASA-CASA model does accurately predict the average annual productivity of a site like the BEF.  相似文献   

17.
Field surveys are often a primary source of aboveground biomass (AGB) data, but plot-based estimates of parameters related to AGB are often not sufficiently precise, particularly not in tropical countries. Remotely sensed data may complement field data and thus help to increase the precision of estimates and circumvent some of the problems with missing sample observations in inaccessible areas. Here, we report the results of a study conducted in a 15,867 km² area in the dry miombo woodlands of Tanzania, to quantify the contribution of existing canopy height and biomass maps to improving the precision of canopy height and AGB estimates locally. A local and a global height map and three global biomass maps, and a probability sample of 513 inventory plots were subject to analysis. Model-assisted sampling estimators were used to estimate mean height and AGB across the study area using the original maps and then with the maps calibrated with local inventory plots. Large systematic map errors – positive or negative – were found for all the maps, with systematic errors as great as 60–70 %. The maps contributed nothing or even negatively to the precision of mean height and mean AGB estimates. However, after being calibrated locally, the maps contributed substantially to increasing the precision of both mean height and mean AGB estimates, with relative efficiencies (variance of the field-based estimates relative to the variance of the map-assisted estimates) of 1.3–2.7 for the overall estimates. The study, although focused on a relatively small area of dry tropical forests, illustrates the potential strengths and weaknesses of existing global forest height and biomass maps based on remotely sensed data and universal prediction models. Our results suggest that the use of regional or local inventory data for calibration can substantially increase the precision of map-based estimates and their applications in assessing forest carbon stocks for emission reduction programs and policy and financial decisions.  相似文献   

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.
Spatial predictions of forest variables are required for supporting modern national and sub-national forest planning strategies, especially in the framework of a climate change scenario. Nowadays methods for constructing wall-to-wall maps and calculating small-area estimates of forest parameters are becoming essential components of most advanced National Forest Inventory (NFI) programs. Such methods are based on the assumption of a relationship between the forest variables and predictor variables that are available for the entire forest area. Many commonly used predictors are based on data obtained from active or passive remote sensing technologies. Italy has almost 40% of its land area covered by forests. Because of the great diversity of Italian forests with respect to composition, structure and management and underlying climatic, morphological and soil conditions, a relevant question is whether methods successfully used in less complex temperate and boreal forests may be applied successfully at country level in Italy.For a study area of more than 48,657 km2 in central Italy of which 43% is covered by forest, the study presents the results of a test regarding wall-to-wall, spatially explicit estimation of forest growing stock volume (GSV) based on field measurement of 1350 plots during the last Italian NFI. For the same area, we used potential predictor variables that are available across the whole of Italy: cloud-free mosaics of multispectral optical satellite imagery (Landsat 5 TM), microwave sensor data (JAXA PALSAR), a canopy height model (CHM) from satellite LiDAR, and auxiliary variables from climate, temperature and precipitation maps, soil maps, and a digital terrain model.Two non-parametric (random forests and k-NN) and two parametric (multiple linear regression and geographically weighted regression) prediction methods were tested to produce wall-to-wall map of growing stock volume at 23-m resolution. Pixel level predictions were used to produce small-area, province-level model-assisted estimates. The performances of all the methods were compared in terms of percent root mean-square error using a leave-one-out procedure and an independent dataset was used for validation. Results were comparable to those available for other ecological regions using similar predictors, but random forests produced the most accurate results with a pixel level R2 = 0.69 and RMSE% = 37.2% against the independent validation dataset. Model-assisted estimates were more precise than the original design-based estimates provided by the NFI.  相似文献   

20.
Land use and land cover (LULC) change detection associated with oil and gas activities plays an important role in effective sustainable management practices, compliance monitoring, and reclamation assessment. In this study, a mapping methodology is presented for quantifying pre- and post-disturbance LULC types with annual Landsat Best-Available-Pixel multispectral data from 2005 to 2013. Annual LULC and land disturbance maps were produced for one of the major conventional oil and gas production areas in West-Central Alberta with an accuracy of 78% and 87%, respectively. The highest rate of vegetation loss (178 km2/year) was observed in coniferous forest compared to broadleaf forest, mixed forest, and native vegetation. Integration of ancillary oil and gas geospatial data with annual land disturbances indicated that less than 20% of the total land disturbances were attributable to oil and gas activities. In 2013, approximately 44% of oil and gas disturbances from 2005 to 2013 showed evidence of vegetation recovery. In the future, geospatial data related to wildfire, logging activities, insect defoliation, and other natural and anthropogenic factors can be integrated to quantify other causes of land disturbances.  相似文献   

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

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