首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 796 毫秒
1.
In modeling forest stand growth and yield,crown width,a measure for stand density,is among the parameters that allows for estimating stand timber volumes.However,accurately measuring tree crown size in the field,in particu-lar for mature trees,is challenging.This study demonstrated a novel method of applying machine learning algorithms to aerial imagery acquired by an unmanned aerial vehi-cle (UAV) to identify tree crowns and their widths in two loblolly pine plantations in eastern Texas,USA.An ortho mosaic image derived from UAV-captured aerial photos was acquired for each plantation (a young stand before canopy closure,a mature stand with a closed canopy).For each site,the images were split into two subsets:one for training and one for validation purposes.Three widely used object detection methods in deep learning,the Faster region-based convolutional neural network (Faster R-CNN),You Only Look Once version 3 (YOLOv3),and single shot detection(SSD),were applied to the training data,respectively.Each was used to train the model for performing crown recogni-tion and crown extraction.Each model output was evaluated using an independent test data set.All three models were successful in detecting tree crowns with an accuracy greater than 93%,except the Faster R-CNN model that failed on the mature site.On the young site,the SSD model performed the best for crown extraction with a coefficient of determination(R2) of 0.92,followed by Faster R-CNN (0.88) and YOLOv3(0.62).As to the mature site,the SSD model achieved a R2 as high as 0.94,follow by YOLOv3 (0.69).These deep leaning algorithms,in particular the SSD model,proved to be successfully in identifying tree crowns and estimat-ing crown widths with satisfactory accuracy.For the pur-pose of forest inventory on loblolly pine plantations,using UAV-captured imagery paired with the SSD object deten-tion application is a cost-effective alternative to traditional ground measurement.  相似文献   

2.
Crown shyness or canopy disengagement, the phenomenon wherein gaps around trees develop from swaying, whipping and shading, has been identified in the literature since the 1920s. Recent results by researchers at the University of Alberta have clearly described many of the processes involved for lodgepole pine [e.g. Rudnicki, M., Silins, U., Lieffers, V.J., Josi, G., 2001. Measure of simultaneous tree sways and estimation of crown interactions among a group of trees. Trees 15, 83–90; Rudnicki, M., Lieffers, V.J., Silins, U., 2003. Stand structure governs the crown collisions of lodgepole pine, Canadian Journal of Forestry Research 33, 1238–1244; Rudnicki, M., Silins, U., Lieffers, V.J., 2004. Crown cover is correlated with relative density, tree slenderness, and tree height in lodgepole pine. Forest Science 50, 356–363; Fish, H., Lieffers, V.J., Silins, U., Hall, R.J., 2006. Crown shyness in lodgepole pine stands of varying stand height, density, and site index in the upper foothills of Alberta. Canadian Journal of Forestry Research 9, 2104–2111]. However, explicit models of crown shyness are sparse in the literature. This paper describes the development of empirical models of crown shyness in lodgepole pine for British Columbia (BC). We measured crown area and neighbour locations on 60 trees growing in 13 stands in central BC. We estimated potential crown area (AV) using stem maps and Voronoi polygons constrained by estimates of maximum crown width, and then related observed crown area (AC) to AV and additional individual tree variables. One of the nine prediction equations was coded into a spatially explicit tree growth model modified to evaluate the effects of crown shyness at the stand level. Crown shyness models validated well against two independent sources and when linked with a light model tRAYci [Brunner, A., 1998. A light model for spatially explicit forest stand models. Forest Ecology and Management 107, 19–46], increased the below-canopy light by 0.07–0.11.  相似文献   

3.
A crown profile model was developed for radiata pine (Pinus radiata D. Don) in Galicia (northwestern Spain). Data from 443 trees located in 56 permanent plots, established by the Unidade de Xestión Forestal Sostible (UXFS) of the University of Santiago de Compostela in plantations of this species in the region, were used. The crowns of the trees were measured by a visual method based on similar triangles. Both simple geometric shapes and mathematical equations were used to describe the crown profiles. As crown profile models usually require variables that are expensive to measure, equations to estimate the maximum crown radius and the height to the maximum crown radius were also developed, using other easily measured tree and stand variables. Several models were fitted using a system of equations approach and accounting for an autocorrelated, heteroscedastic error structure. The selected crown profile model consists of a system of two allometric equations for the crown below (primarily shade needles) and above (primarily sun needles) maximum crown radius. The model explained 88% of the variability in crown radius with a mean error of 0.24 m.  相似文献   

4.
Cork oak (Quercus suber) forests are acknowledged for their biodiversity and economic (mainly cork production) values. Wildfires are one of the main threats contributing to cork oak decline in the Mediterranean Basin, and one major question that managers face after fire in cork oak stands is whether the burned trees should be coppiced or not. This decision can be based on the degree of expected crown regeneration assessed immediately after fire. In this study we carried out a post-fire assessment of the degree of crown recovery in 858 trees being exploited for cork production in southern Portugal, 1.5 years after a wildfire. Using logistic regression, we modelled good or poor crown recovery probability as a function of tree and stand variables. The main variables influencing the likelihood of good or poor crown regeneration were bark thickness, charring height, aspect and tree diameter. We also developed management models, including simpler but easier to measure variables, which had a lower predictive power but can be used to help managers to identify, immediately after fire, trees that will likely show good crown regeneration, and trees that will likely die or show poor regeneration (and thus, potential candidates for trunk coppicing).  相似文献   

5.
Forest floor characteristics influence nutrient cycling and energy flow properties of forest ecosystems, and determine quality of habitat for many forest plants and animals. Differential crown recession and crown development among stands of differing density suggest that an opportunity may exist to control the input of fine woody litter into the system by manipulating stand density. The objective was to measure the rate of branch mortality among stands of differing density and to estimate the range in total per hectare necromass inputs. Although litter traps are reliable for estimating per hectare rates of litterfall, branch mortality dating on sectioned stems uniquely allows assessment of several other litterfall components: (1) individual tree contributions to total litterfall; (2) the amount of branch material released by mortality, regardless of whether the branches are shed to the forest floor; (3) the distribution of basal diameters characterizing the litterfall from a given tree and stand. Twenty-four trees were felled and sectioned on permanent plots that were part of a silvicultural study of stand density regimes in Douglas-fir (Pseudotsuga menziesii (Mirbel) Franco.). Whorl branches were dissected out of bole sections to determine the dates of mortality, and a branch biomass equation was applied to estimate potential rate of litterfall. Periodic annual rates were expressed in four ways: (1) number of branches per tree; (2) mass of branches per tree; (3) mass of branches per unit of crown projection area; (4) mass of branches per hectare. For the growth periods investigated, larger trees and trees growing on denser plots tended to release a greater necromass through branch mortality. Average branch basal diameter generally decreased with increasing stand density. Annual branch mortality ranged from 33 to 430 g m−2 crown projection area for individual trees, and from 236 to 1035 kg ha−1 for individual plots. These rates approached the low end of the range of previously published fine litterfall rates for Douglas-fir. Rates on these plots were relatively low owing to the temporary delay in crown recession imposed by artificial thinning. A conceptual model of branch litter dynamics is presented to depict consistencies with crown development among stands managed under different density regimes.  相似文献   

6.
Crown dimensions are important for the quantification of tree interactions in some growth models. This study investigates the potential for structural indices and other spatial measures to improve the prediction of crown radius and crown length for birch (Betula spp.) and Sitka spruce (Picea sitchensis (Bong.) Carr.) in forests in Wales. Crown dimensions were measured for 125 birch and 154 spruce in six fully stem-mapped research plots. These data were used to test the performance of a crown radius model and a crown length model which estimated crown dimensions on the basis of allometric relationships with stem dimensions. Spatial data from the six plots were used to calculate the structural indices mean directional index, diameter correlation index, species mingling, dbh and height dominance, and dbh differentiation, as well as the Hegyi competition index, and basal area of neighbours and larger neighbours, for each crown measurement sample tree, using various numbers of nearest neighbours. Two non-spatial indices, BAL and BALMOD, were also calculated for all sample trees for comparison. These spatial and non-spatial variables were then incorporated into modified crown dimension models. Model performances, in terms of efficiency and relative bias, were compared to determine whether the inclusion of spatial or non-spatial variables resulted in any improvements over models using tree dimensions alone. Crown length and radius were found to be correlated with most of the spatial measures studied. Models incorporating spatial variables gave improvements in performance over allometric models for every data set, and performed more consistently than models containing non-spatial variables. The greatest improvements were achieved for suppressed birch in unthinned forests which had irregularly shaped and strongly displaced crowns. The spatial variable contributing to the most efficient model for each data set varied widely. This points to the complexity of tree spatial interactions and indicates that there is a great deal of scope for investigating other structural indices and crown dimension model forms.  相似文献   

7.
Trees are able to respond to their local biotic and abiotic environment with morphological adjustments which improve resource acquisition and, thus, growth. In forests, light is broadly recognised as one of the major factors determining growth, and morphological responses comprise changes in crown architecture and stem stature. On sloping terrain, the interplay of phototropism and gravitropism may further affect morphological growth characteristics. However, different tree species are expected to show species-specific responses. In this study, we analysed three growth characteristics of tree individuals belonging to four species of two functional groups (evergreen: Schima superba, Castanopsis eyrei, deciduous: Quercus serrata var. brevipetiolata, Castanea henryi) in a species-rich Chinese subtropical forest. Crown projection area, relative crown displacement and stem inclination were related to biotic (local species richness, functional richness, competition, stand age) and abiotic (slope aspect and inclination, soil depth) variables in the local neighbourhood of the tree individuals. We hypothesised that (i) there are species-specific differences in the morphological response of crown architecture and stem stature and (ii) that crown size and asymmetry as well as stem inclination are influenced by both, biotic and abiotic factors. In contrast to our expectations we were unable to reveal any species-specific differences in any of the three growth characteristics. The results of mixed effects models showed that crown area was mainly affected by the target tree's dbh and biotic variables related to neighbours (competition, functional diversity), whereas stem inclination was mainly influenced by slope. Relative crown displacement was influenced by both, biotic and abiotic variables. We conclude that growth responses resulting in crown displacement and stem inclination seem to be an important mechanism to ameliorate foraging for light in our study area, but that these responses appear to be species-independent. The interplay of stem inclination and crown displacement allows for a plastic response of tree individuals in biotically and abiotically heterogeneous environments. Our results indicate that forest management in this region should focus on functionally diverse stands which are promoting crown area positively resulting in increased growth rates of individual trees.  相似文献   

8.
参数化林木个体及林分场景可视化模拟技术   总被引:1,自引:0,他引:1       下载免费PDF全文
以湖南攸县黄丰桥国有林场为试验区,以1块林分边界规整的典型杉木人工同龄纯林为研究对象,测定林木胸径、树高、冠幅、冠高与活枝下高等测树因子及林木位置信息。研究冠形曲线函数,利用Direct 3D技术,实现冠形控制下的参数化林木模拟。研究林木模型格式转换方法,为林分场景构建提供模型数据,结合MOGRE技术,研究林分场景模拟技术方法,实现林分场景可视化模拟。结果表明:参数化林木个体建模算法简单适用,测树因子信息可直接用于林木个体可视化模拟,所建模型可体现林木个体差异,形态逼真;利用MOGRE技术可快速、有效模拟林分场景,所建场景真实感较强。此方法所需数据易于获取,适用性强,可对林木个体及林分场景进行逼真模拟,可为研究林木竞争关系、林分生长、林分结构与林分经营提供新思路。  相似文献   

9.
【目的】无人机机载激光雷达能够准确地测定单木、林分乃至大尺度森林结构参数(树高和树冠因子)。为应用无人机激光雷达技术准确估测森林蓄积量、生物量和碳储量提供计量依据和技术支撑。【方法】以150株实测马尾松生物量样本数据为研究对象,采用非线性回归估计方法和度量误差联立方程组方法,分析立木材积和地上生物量与树高、树冠因子的相关性,并在此基础上研究建立基于树高和树冠因子的立木材积与地上生物量相容模型。【结果】单株材积和地上生物量与树高因子的相关性最为紧密,其次才是树冠因子;基于树高和冠幅因子的二元材积和地上生物量模型预估精度较高,达到92%以上,再考虑冠长因子的三元模型预估精度改进不大;基于树高和冠幅因子的二元立木材积与地上生物量相容模型估计效果更好,相对于一元相容模型系统而言,二元相容模型拟合效果有较大幅度提高,预估精度达到92%以上。【结论】采用度量误差联立方程组方法可以有效解决基于树高和树冠因子的立木材积与地上生物量相容问题,并且预估精度达到92%以上,所建二元立木材积与地上生物量相容模型可为应用激光雷达技术反演森林蓄积量和生物量提供计量依据。  相似文献   

10.
Korean larch(Larix olgensis) is one of the main tree species for afforestation and timber production in northeast China. However, its timber quality and growth ability are largely infl uenced by crown size, structure and shape. The majority of crown models are static models based on tree size and stand characteristics from temporary sample plots, but crown dynamic models has seldom been constructed. Therefore, this study aimed to develop height to crown base( HCB) and crown length( CL) dynamic models using the branch mortality technique for a Korean larch plantation. The nonlinear mixed-effects model with random effects, variance functions and correlation structures, was used to build HCB and CL dynamic models. The data were obtained from 95 sample trees of 19 plots in Meng JiaGang forest farm in Northeast China. The results showed that HCB progressively increases as tree age, tree height growth(HT growth) and diameter at breast height growth( DBH growth).The CL was increased with tree age in 20 years ago, and subsequently stabilized. HT growth, DBH growth stand basal area( BAS) and crown competition factor( CCF) signifi cantly infl uenced HCB and CL. The HCB was positively correlated with BAS, HT growth and D BH growth, but negatively correlated with CCF. The CL was positively correlated with BAS and CCF, but negatively correlated with D BH growth. Model fi tting and validation confi rmed that the mixed-effects model considering the stand and tree level random effects was accurate and reliable for predicting the HCB and CL dynamics. However, the models involving adding variance functions and time series correlation structure could not completely remove heterogeneity and autocorrelation, and the fi tting precision of the models was reduced. Therefore, from the point of view of application, we should take care to avoid setting up over-complex models. The HCB and CL dynamic models in our study may also be incorporated into stand growth and yield model systems in China.  相似文献   

11.
The crown profile, volume, and surface area obtained from 25 trees in a stand of 29-year-old Japanese cypress (Chamaecyparis obtusa Endl.) using a Spiegel relascope were compared with those obtained by direct measurement after tree felling. In 17 of the 25 trees, the crown profile measured with the relascope was identical to the crown profile obtained by direct measurement. For the 17 trees, the mean relative errors of crown volume and surface area tended to be unbiased. On the other hand, the crown profiles obtained by the two methods were different in five of the 25 trees, while it could not been judged whether the crown profiles obtained by the two methods were identical or different in three of the 25 trees. However, for each group of five or three trees, the mean relative errors of crown volume and surface area were also unbiased. Moreover, for both crown volume and surface area, no significant differences were detected in variance of error among the groups of 17, five, and three trees. Therefore, the discrepancy between crown profiles obtained by the two methods for the five and three trees might be negligible in terms of the accuracy and precision of crown volume and surface area measurements. The standard errors for both crown volume and surface area were larger in each of the five and three trees than in the 17 trees. However, even for the five and three trees, the crown volume and surface area obtained using the relascope were more precise than those based on the assumption that crown form was a cone or paraboloid. We concluded that using a Spiegel relascope would enable us to obtain an accurate and precise crown profile, volume, and surface area of Japanese cypress trees for practical use, especially when the main objective of measurement is determining crown volume or surface area.  相似文献   

12.
Static models of individual tree crown attributes such as height to crown base and maximum branch diameter profile have been developed for several commercially important species. Dynamic models of individual branch growth and mortality have received less attention, but have generally been developed retrospectively by dissecting felled trees; however, this approach is limited by the lack of historic stand data and the difficulty in determining the exact timing of branch death. This study monitored the development of individual branches on 103 stems located on a variety of silvicultural trials in the Pacific Northwest, USA. The results indicated that branch growth and mortality were significantly influenced by precommercial thinning (PCT), commercial thinning, fertilization, vegetation management, and a foliar disease known as Swiss needle cast [caused by Phaeocryptopus gaeumannii (T. Rohde) Petr.]. Models developed across these datasets accounted for treatment effects through variables such as tree basal area growth and the size of the crown. Insertion of the branch growth and mortality equations into an individual-tree modeling framework, significantly improved short-term predictions of crown recession on an independent series of silvicultural trials, which increased mean accuracy of diameter growth prediction (reduction in mean bias). However, the static height to crown base equation resulted in a lower mean square error for the tree diameter and height growth predictions. Overall, individual branches were found to be highly responsive to changes in stand conditions imposed by silvicultural treatments, and therefore represent an important mechanism explaining tree and stand growth responses.  相似文献   

13.
毛竹混交林林分生长状况和结构特征研究   总被引:3,自引:0,他引:3  
林桂明 《世界竹藤通讯》2012,10(2):18-22,26
对福建省建阳交溪毛竹阔叶树混交林、毛竹杉木混交林、毛竹杉木阔叶树混交林及毛竹纯林等6种不同混交模式林分生长状况进行比较,结果表明:混交林毛竹个体生物量、毛竹胸径和竹高较大,不同模式林分间的胸径、竹高、冠幅等生长因子均达显著差异;竹阔混交林较有利于竹鞭孕芽发笋;各模式竹林结构综合评判为一般经营竹林;竹林竹鞭长、鞭径、鞭节生长以竹阔比例为5:5林分竹鞭较好。  相似文献   

14.
Decreasing growth rate (Feeley et al., 2007) and large die-back due to drought (Phillips et al., 2009) suggest that tropical forests are suffering recent climate changes. Forest vulnerability to external factors (e.g. air pollution, acid rain) is widely studied in northern countries, while only a few attempts have investigated crown integrity in the Tropics. The method needs to be generic enough to account for the large number of species and crown shapes encountered in tropical forests. In the present study, we developed and tested a novel field method that estimates crown fragmentation (main branch mortality (MB) and secondary branch mortality (SB)), liana infestation (LI) and crown position (CP) in the canopy. The relationship between crown fragmentation and annual growth rate (agr) was investigated through multiple regression. Six out of eight canopy tree species showed significant growth decline with increasing crown fragmentation. Higher probability of death was also found in trees with severe crown fragmentation. The capacity of such crown assessment to depict tree vitality in a forest stand is discussed along with potential applications in both forest science and management.  相似文献   

15.
长白落叶松人工林树冠形状的模拟   总被引:17,自引:3,他引:17  
李凤日 《林业科学》2004,40(5):16-24
以长白山地区 2 6a生长白落叶松人工林为研究对象 ,采用枝解析的方法 ,测定了 2 5株林木 (直径 10 5~2 4 9cm)的树冠变量 ,并建立了预测树冠外侧形状的冠形模型。基于枝条着枝深度 (DINC)和林木变量所建立的树冠形状模型包括 :基径、枝长、着枝角度和弦长等预估模型。对于大小相同树木的主要枝条来讲 ,这些树冠变量是随着DINC的增加而增大 ;而林木的胸径 (DBH)和树高 (HT)变量很好地反映了不同大小树木的冠形变化。冠形预测模型预测效果良好 ,充分体现了树冠结构的变化趋势 :树冠形状在树冠的中上部呈抛物线体 ,而在树冠的下部则为近圆柱体。文中所建模型 ,可以合理地描述长白落叶松人工林的树冠形状及其变化规律。  相似文献   

16.
Competition for canopy space is a fundamental structuring feature of forest ecosystems and remains an enduring focus of research attention. We used a spatial neighborhood approach to quantify the influence of local competition on the size of individual tree crowns in north-central British Columbia, where forests are dominated by subalpine fir (Abies lasiocarpa), lodgepole pine (Pinus contorta) and interior spruce (Picea glauca × engelmanii). Using maximum likelihood methods, we quantified crown radius and length as functions of tree size and competition, estimated by the species identity and spatial arrangement of neighboring trees. Tree crown size depended on tree bole size in all species. Given low levels of competition, pine displayed the widest, shortest tree crowns compared to the relatively long and narrow crowns found in spruce and fir. Sensitivity to crowding by neighbors declined with increasing tree height in all but the pine crown radius model. Five of the six selected best models included separate competition coefficients for each neighboring tree species, evidence that species generally differ in their competitive effects on neighboring tree crowns. The selected crown radius model for lodgepole pine, a shade-intolerant species, treated all neighbors as equivalent competitors. In all species, competition from neighbors exerted an important influence on crown size. Per-capita effects of competition across different sizes and species of neighbors and target trees varied, but subalpine fir generally displayed the strongest competitive effects on neighbors. Results from this study provide evidence that species differ both in their response to competition and in their competitive influence on neighbors, factors that may contribute to maintaining coexistence.  相似文献   

17.
Since individual tree leaf area is an important measure for productivity as well as for site occupancy, it is of high interest in many studies about forest growth. The exact determination of leaf area is nearly impossible. Thus, a common way to get information about leaf area is to use substitutes. These substitutes are often variables which are collected in a destructive way which is not feasible for long term studies. Therefore, this study aimed at testing the applicability of using substitutes for leaf area which could be collected in a non-destructive way, namely crown surface area and crown projection area. In 8 stands of Norway spruce (Picea abies L. Karst.), divided into three age classes and two thinning treatments, a total of 156 trees were felled in order to test the relationship between leaf area and crown surface area and crown projection area, respectively. Individual tree leaf area of the felled sample trees was estimated by 3P-branch sampling with an accuracy of ±10%. Crown projection area and crown surface area were compared with other, more commonly used, but destructive predictors of leaf area, namely sapwood area at different heights on the bole. Our investigations confirmed findings of several studies that sapwood area is the most precise measure for leaf area because of the high correlation between sapwood area and the leaf area. But behind sapwood area at crown base and sapwood area at three tenth of the tree height the predictive ability of crown surface area was ranked third and even better than that of sapwood area at breast height (R2 = 0.656 compared with 0.600). Within the stands leaf area is proportional to crown surface area. Using the pooled data of all stands a mixed model approach showed that additionally to crown surface area dominant height and diameter at breast height (dbh) improved the leaf area estimates. Thus, taking dominant height and dbh into account, crown surface area can be recommended for estimating the leaf area of individual trees. The resulting model was in line with many other findings on the leaf area and leaf mass relationships with crown size. From the additional influence of dominant height and dbh in the leaf area model we conclude that the used crown model could be improved by estimating the position of the maximum crown width and the crown width at the base of the crown depending on these two variables.  相似文献   

18.
天然林分生长模型在小班数据更新中的应用   总被引:20,自引:0,他引:20  
当小班调查得到各组成树种的平均高,平均直径以及小班的年龄,每公顷株数,郁闭度,散生木蓄积等因子后,应用天然林分生长模型可以实现无人为干预小班数据的全面更新和连续预测。这些模型应用于小班数据更新时的特点是分别组成树种,以地位级指数评价立地质量,以预估间隔期代替年龄,更适用于天然混交异龄林小班数据的更新。通过吉林省汪清林业局108块复位样地的验证,表明天然林分生长模型的预估结果优于目前生产单位所用的生  相似文献   

19.
There is currently great interest in improving the applicability of forest gap models to changing environmental conditions, in order to facilitate the assessment of possible impacts of climatic change on forest ecosystems. Moreover, for the development of mitigation strategies, it is necessary to include forest management options in the models. Both the simulation of transient effects of climatic change and of forest management regimens require a realistic representation of stand structure in gap models, since tree species respond to variations in stand density in characteristic ways, depending on their ecological strategies.In this study, we compared the effect of five different height growth functions that are sensitive to stand density on simulated stand structure of the FORSKA forest gap model. We used long term observation data from a beech thinning trial at Fabrikschleichach, Bavaria, to test the alternative functions. First, we compared simulation results of the original FORSKA model with measured stand development from 1870 to 1990. Whereas simulated stand level variables (e.g. biomass, mean diameter and height) showed good correspondence with observations, individual tree dimensions and simulated stand structure were quite unrealistic. After calibrating parameters of the height growth functions with data from a lightly thinned plot at Fabrikschleichach, we ran the model with data from a heavily thinned plot for validation. All five functions considerably improved the simulation of height/diameter relationships and stand structure. However, there were distinct differences between functions. The best correspondence with measurements was shown by a function which uses the relative radiation intensity in the centre of a tree crown as an indicator of the competition status of the tree. This function is rather simple and needs only two growth parameters, which can be derived for different functional types of species, according to their shade tolerance.With the new, flexible height growth function it should be possible to extend the applicability of gap models to more realistic simulation experiments including forest management and natural disturbance. To our knowledge, this was the first attempt to employ long term forest observation data for the calibration and validation of a forest gap model. The results suggest that such data could be very useful in model testing and improvement.  相似文献   

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

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