首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 608 毫秒
1.
Optimizing nitrogen (N) fertilization in crop production by in-season measurements of crop N status may improve fertilizer N use efficiency. Hyperspectral measurements may be used to assess crop N status by estimating leaf chlorophyll content. This study evaluated the ability of the PROSAIL canopy-level reflectance model to predict leaf chlorophyll content. Trials were conducted with two potato cultivars under different N fertility rates (0–300 kg N ha−1). Canopy reflectance, leaf area index (LAI) and leaf chlorophyll and N contents were measured. The PROSAIL model was able to predict leaf chlorophyll content with reasonable accuracy later in the growing season. The low estimation accuracy earlier in the growing season could be due to model sensitivity to non-homogenous canopy architecture and soil background interference before full canopy closure. Canopy chlorophyll content (leaf chlorophyll content × LAI) was predicted less accurately than leaf chlrophyll content due to the low estimation accuracy of LAI for values higher than 4.5.  相似文献   

2.
The retrieval of canopy biophysical variables is known to be affected by confounding factors such as plant type and background reflectance. The effects of soil type and plant architecture on the retrieval of vegetation leaf area index (LAI) from hyperspectral data were assessed in this study. In situ measurements of LAI were related to reflectances in the red and near-infrared and also to five widely used spectral vegetation indices (VIs). The study confirmed that the spectral contrast between leaves and soil background determines the strength of the LAI–reflectance relationship. It was shown that within a given vegetation species, the optimum spectral regions for LAI estimation were similar across the investigated VIs, indicating that the various VIs are basically summarizing the same spectral information for a given vegetation species. Cross-validated results revealed that, narrow-band PVI was less influenced by soil background effects (0.15 ≤ RMSEcv ≤ 0.56). The results suggest that, when using remote sensing VIs for LAI estimation, not only is the choice of VI of importance but also prior knowledge of plant architecture and soil background. Hence, some kind of landscape stratification is required before using hyperspectral imagery for large-scale mapping of vegetation biophysical variables.  相似文献   

3.
The current development of satellite technology particularly in the sensors like POLDER and MISR, has emphasized more on directional reflectance measurements (i.e. spectral reflectance of the target measured from different view zenith and azimuth angles) of the earth surface features mainly the vegetation for retrieval of biophysical parameters at regional scale using radiative transfer models. This approach being physical process based and uses directional reflectance measurement has been found to better and more reliable compared to the conventional statistical approach used till date and takes care of anisotropic nature (i.e. reflectance from the target is different if measured from different view angles) of the target. Keeping this in view a field experiment was conducted in mustard crop to evaluate the radiative transfer model for biophysical parameter retrieval through its inversion with the objectives set as (i) to relate canopy biophysical parameters and geometry to its bidirectional reflectance, (ii) to evaluate a canopy reflectance model to best represent the radiative transfer within the canopy for its inversion and (iii) to retrieve crop biophysical parameters through inversion of the model. Two varieties of the mustard crop (Brassica juncea L) were grown with two nitrogen treatments. The bidirectional reflectance data obtained at 5 nm interval for a range of 400–1100 nm were integrated to IRS LISS–II sensor’s four band values using Newton Cotes Integration technique. Biophysical parameters like leaf area index, leaf chlorophyll content, leaf length, plant height and average leaf inclination angle, biomass etc were estimated synchronizing with the bi-directional reflectance measurements. Radiative transfer model PROSAIL model was validated and its inversion was done to retrieve LAI and ALA. Look Up Table (LUT) of Bidirectional reflectance distribution function (BRDF) was prepared simulating through PROSAIL model varying only LAI (0.2 interval from 1.2 to 5.4 ) and ALA (5° interval from 40° to 55°) parameters and inversion was done using a merit function and numerical optimization technique given by Press et al. (1986). The derived LAI and ALA values from inversion were well matched with observed one with RMSE 0.521 and 5.57, respectively.  相似文献   

4.
高光谱反演水稻叶面积指数的主成分分析法   总被引:1,自引:0,他引:1  
为了通过水稻冠层反射光谱来提取水稻叶面积指数信息,尝试利用辐射传输模型PROSPECT+SAIL来模拟水稻冠层反射光谱, 比较了各植被指数中叶面积指数(LAI)和叶绿素浓度的相关性。在观察光谱曲线后发现,红边位置光谱可以较好地区分LAI和叶绿素 浓度二者引起光谱变化的差异。由此提出对700 nm~750 nm区间内的反射光谱做主成分变换,并利用第2主成分与LAI建立反演模型( 即主成分分析法),取得了较好效果,表明在植被指数趋近于饱和以至于无法区分二者相关性时,主成分分析法可以作为一种简单 而有效提取水稻叶面积指数信息的补充手段。  相似文献   

5.
Canopy water content (CWC) is important for mapping and monitoring the condition of the terrestrial ecosystem. Spectral information related to the water absorption features at 970 nm and 1200 nm offers possibilities for deriving information on CWC. In this study, we compare the use of derivative spectra, spectral indices and continuum removal techniques for these regions. Hyperspectral reflectance data representing a range of canopies were simulated using the combined PROSPECT + SAILH model. Best results in estimating CWC were obtained by using spectral derivatives at the slopes of the 970 nm and 1200 nm water absorption features. Real data from two different test sites were analysed. Spectral information at both test sites was obtained with an ASD FieldSpec spectrometer, whereas at the second site HyMap airborne imaging spectrometer data were also acquired. Best results were obtained for the derivative spectra. In order to avoid the potential influence of atmospheric water vapour absorption bands the derivative of the reflectance on the right slope of the canopy water absorption feature at 970 nm can best be used for estimating CWC.  相似文献   

6.
The red edge position (REP) in the vegetation spectral reflectance is a surrogate measure of vegetation chlorophyll content, and hence can be used to monitor the health and function of vegetation. The Multi-Spectral Instrument (MSI) aboard the future ESA Sentinel-2 (S-2) satellite will provide the opportunity for estimation of the REP at much higher spatial resolution (20 m) than has been previously possible with spaceborne sensors such as Medium Resolution Imaging Spectrometer (MERIS) aboard ENVISAT. This study aims to evaluate the potential of S-2 MSI sensor for estimation of canopy chlorophyll content, leaf area index (LAI) and leaf chlorophyll concentration (LCC) using data from multiple field campaigns. Included in the assessed field campaigns are results from SEN3Exp in Barrax, Spain composed of 35 elementary sampling units (ESUs) of LCC and LAI which have been assessed for correlation with simulated MSI data using a CASI airborne imaging spectrometer. Analysis also presents results from SicilyS2EVAL, a campaign consisting of 25 ESUs in Sicily, Italy supported by a simultaneous Specim Aisa-Eagle data acquisition. In addition, these results were compared to outputs from the PROSAIL model for similar values of biophysical variables in the ESUs. The paper in turn assessed the scope of S-2 for retrieval of biophysical variables using these combined datasets through investigating the performance of the relevant Vegetation Indices (VIs) as well as presenting the novel Inverted Red-Edge Chlorophyll Index (IRECI) and Sentinel-2 Red-Edge Position (S2REP). Results indicated significant relationships between both canopy chlorophyll content and LAI for simulated MSI data using IRECI or the Normalised Difference Vegetation Index (NDVI) while S2REP and the MERIS Terrestrial Chlorophyll Index (MTCI) were found to have the strongest correlation for retrieval of LCC.  相似文献   

7.
This study describes the retrieval of state variables (LAI, canopy chlorophyll, water and dry matter contents) for summer barley from airborne HyMap data by means of a canopy reflectance model (PROSPECT + SAIL). Three different inversion techniques were applied to explore the impact of the employed method on estimation accuracies: numerical optimization (downhill simplex method), a look-up table (LUT) and an artificial neural network (ANN) approach. By numerical optimization (Num Opt), reliable estimates were obtained for LAI and canopy chlorophyll contents (LAI × Cab) with r2 of 0.85 and 0.94 and RDP values of 1.81 and 2.65, respectively. Accuracies dropped for canopy water (LAI × Cw) and dry matter contents (LAI × Cm). Nevertheless, the range of leaf water contents (Cw) was very narrow in the studied plant material. Prediction accuracies generally decreased in the order Num Opt > LUT > ANN. This decrease in accuracy mainly resulted from an increase in offset in the obtained values, as the retrievals from the different approaches were highly correlated. The same decreasing order in accuracy was found for the difference between the measured spectra and those reconstructed from the retrieved variable values. The parallel application of the different inversion techniques to one collective data set was helpful to identify modelling uncertainties, as shortcomings of the retrieval algorithms themselves could be separated from uncertainties in model structure and parameterisation schemes.  相似文献   

8.
In this paper, we carried out a laboratory experiment to study changes in canopy reflectance of Tamarugo plants under controlled water stress. Tamarugo (Prosopis tamarugo Phil.) is an endemic and endangered tree species adapted to the hyper-arid conditions of the Atacama Desert, Northern Chile. Observed variation in reflectance during the day (due to leaf movements) as well as changes over the experimental period (due to water stress) were successfully modelled by using the Soil-Leaf-Canopy (SLC) radiative transfer model. Empirical canopy reflectance changes were mostly explained by the parameters leaf area index (LAI), leaf inclination distribution function (LIDF) and equivalent water thickness (EWT) as shown by the SLC simulations. Diurnal leaf movements observed in Tamarugo plants (as adaptation to decrease direct solar irradiation at the hottest time of the day) had an important effect on canopy reflectance and were explained by the LIDF parameter. The results suggest that remote sensing based assessment of this desert tree should consider LAI and canopy water content (CWC) as water stress indicators. Consequently, we tested fifteen different vegetation indices and spectral absorption features proposed in literature for detecting changes of LAI and CWC, considering the effect of LIDF variations. A sensitivity analysis was carried out using SLC simulations with a broad range of LAI, LIDF and EWT values. The Water Index was the most sensitive remote sensing feature for estimating CWC for values less than 0.036 g/cm2, while the area under the curve for the spectral range 910–1070 nm was most sensitive for values higher than 0.036 g/cm2. The red-edge chlorophyll index (CIred-edge) performed the best for estimating LAI. Diurnal leaf movements had an effect on all remote sensing features tested, particularly on those for detecting changes in CWC.  相似文献   

9.
Computer simulation models have seldom been applied for estimating the structural and biophysical variables of forest canopy. In this study, an approach for the estimation of leaf area index (LAI) using the information contained in hyperspectral, multi-angle images and the inversion of a computer simulation model are explored. For this purpose, L-systems combined with forest growth model ZELIG were applied to render 3-D forest architectural scenarios. The Radiosity-graphics combined model (RGM) was used to estimate forest LAI from the Compact High-Resolution Imaging Spectrometer/Project for On-Board Autonomy (CHRIS/PROBA) data. LAI inversion was performed using the look-up table (LUT) method. The estimated LAI was evaluated against in situ LAI measurement and compared against the LAI predictions from CHRIS data obtained using the Li-Strahler geometric-optical canopy reflectance model (GOMS). The results indicated that the method used in this study can be efficient strategy to estimate LAI by RGM model inversion.  相似文献   

10.
When crops senescence, leaves remain until they fall off or are harvested. Hence, leaf area index (LAI) stays high even when chlorophyll content degrades to zero. Current LAI approaches from remote sensing techniques are not optimized for estimating LAI of senescent vegetation. In this paper a two-step approach has been proposed to realize simultaneous LAI mapping over green and senescent croplands. The first step separates green from brown LAI by means of a newly proposed index, ‘Green Brown Vegetation Index (GBVI)’. This index exploits two shortwave infrared (SWIR) spectral bands centred at 2100 and 2000 nm, which fall right in the dry matter absorption regions, thereby providing positive values for senescent vegetation and negative for green vegetation. The second step involves applying linear regression functions based on optimized vegetation indices to estimate green and brown LAI estimation respectively. While the green LAI index uses a band in the red and a band in the red-edge, the brown LAI index uses bands located in the same spectral region as GBVI, i.e. an absorption band located in the region of maximum absorption of cellulose and lignin at 2154 nm, and a reference band at 1635 nm where the absorption of both water and dry matter is low. The two-step approach was applied to a HyMap image acquired over an agroecosystem at the agricultural site Barrax, Spain.  相似文献   

11.
This study is aimed at demonstrating the feasibility of the large scale LAI inversion algorithms using red and near infrared reflectance obtained from high resolution satellite imagery. Radiances in digital counts were obtained in 10 m resolution acquired on cloud free day of August 23, 2007, by the SPOT 5 high resolution geometric (HRG) instrument on mostly temperate hardwood forest located in the Great Lakes – St. Lawrence forest in Southern Quebec. Normalized difference vegetation index (NDVI), scaled difference vegetation index (SDVI) and modified soil-adjusted vegetation index (MSAVI) were applied to calculate gap fractions. LAI was inverted from the gap fraction using the common Beer–Lambert's law of light extinction under forest canopy. The robustness of the algorithm was evaluated using the ground-based LAI measurements and by applying the methods for the independently simulated reflectance data using PROSPECT + SAIL coupled radiative transfer models. Furthermore, the high resolution LAI was compared with MODIS LAI product. The effects of atmospheric corrections and scales were investigated for all of the LAI retrieval methods. NDVI was found to be not suitable index for large scale LAI inversion due to the sensitivity to scale and atmospheric effects. SDVI was virtually scale and atmospheric correction invariant. MSAVI was also scale invariant. Considering all sensitivity analysis, MSAVI performed best followed by SDVI for robust LAI inversion from high resolution imagery.  相似文献   

12.
Leaf chlorophyll content is an important variable for agricultural remote sensing because of its close relationship to leaf nitrogen content. The triangular greenness index (TGI) was developed based on the area of a triangle surrounding the spectral features of chlorophyll with points at (670 nm, R670), (550 nm, R550), and (480 nm, R480), where Rλ is the spectral reflectance at wavelengths of 670, 550 and 480, respectively. The equation is TGI = −0.5[(670  480)(R670  R550)  (670  550)(R670  R480)]. In 1999, investigators funded by NASA's Earth Observations Commercialization and Applications Program collaborated on a nitrogen fertilization experiment with irrigated maize in Nebraska. Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) data and Landsat 5 Thematic Mapper (TM) data were acquired along with leaf chlorophyll meter and other data on three dates in July during late vegetative growth and early reproductive growth. TGI was consistently correlated with plot-averaged chlorophyll-meter values at the spectral resolutions of AVIRIS, Landsat TM, and digital cameras. Simulations using the Scattering by Arbitrarily Inclined Leaves (SAIL) canopy model indicate an interaction among TGI, leaf area index (LAI) and soil type at low crop LAI, whereas at high LAI and canopy closure, TGI was only affected by leaf chlorophyll content. Therefore, TGI may be the best spectral index to detect crop nitrogen requirements with low-cost digital cameras mounted on low-altitude airborne platforms.  相似文献   

13.
This paper reports a series of laboratory and field measurements of spectral reflectance under artificial and natural light conditions which demonstrate that effects of natural chlorophyll fluorescence are observable in the reflectance red edge spectral region. These are results from the progress made to link physiologically-based indicators to optical indices from hyperspectral remote sensing in the Bioindicators of Forest Sustainability Project. This study is carried out on twelve sites of Acer saccharum M. in the Algoma Region, Ontario (Canada), where field measurements, laboratory-simulation experiments, and hyperspectral CASI imagery have been carried out in 1997, 1998, 1999 and 2000 campaigns. Leaf samples from the study sites have been used for reflectance and transmittance measurements with the Li-Cor Model 1800 integrating sphere apparatus coupled to an Ocean Optics Model ST1000 fibre spectrometer in which the same leaves are illuminated alternatively with and without fluorescence-exciting radiation. A study of the diurnal change in leaf reflectance spectra, combined with fluorescence measurements with the PAM-2000 Fluorometer show that the difference spectra are consistent with observed diurnal changes in steady-state fluorescence. Small canopies of Acer saccharum M. have been used for laboratory measurements with the CASI hyperspectral sensor, and under natural light conditions with a fibre spectrometer in diurnal trials, in which the variation of measured reflectance is shown experimentally to be consistent with a fluorescence signature imposed on the inherent leaf reflectance signature. Such reflectance changes due to CF are measurable under natural illumination conditions, although airborne experiments with the CASI hyperspectral sensor produced promising but less convincing results in two diurnal experiments carried out in 1999 and 2000, where small variations of reflectance due to the effect of CF were observed.  相似文献   

14.
Hyperspectral remote sensing has demonstrated great potential for accurate retrieval of canopy water content (CWC). This CWC is defined by the product of the leaf equivalent water thickness (EWT) and the leaf area index (LAI). In this paper, in particular the spectral information provided by the canopy water absorption feature at 970 nm for estimating and predicting CWC was studied using a modelling approach and in situ spectroradiometric measurements. The relationship of the first derivative at the right slope of the 970 nm water absorption feature with CWC was investigated with the PROSAIL radiative transfer model and tested for field spectroradiometer measurements on two test sites. The first site was a heterogeneous floodplain with natural vegetation like grasses and various shrubs. The second site was an extensively grazed fen meadow.  相似文献   

15.
The results emerged out of the studies on spectral reflectance under normal and nitrogen and phosphorus stress condition in soybean (Glycine max L.) conducted at Marathwada Agricultural University experimental farm, Parbhani duringkharif 2004–05 showed that crop growth and bio-physiological parameters viz., Height, chlorophyll, leaf area index and total biomass influenced by pest and disease and nutrient stress resulted in detectable spectral reflectance variation. Poor crop growth, reduced canopy cover, chlorophyll content and biomass production are the effects observed in nutrient deficient crops. These above changes in soybean crop were related to spectral indices (RVI and NDVI) that are resulted in discrimination of stressed and normal (non-stressed) soybean crop.  相似文献   

16.
A time series of leaf area index (LAI) has been developed based on 16-day normalized difference vegetation index (NDVI) data from the Moderate Resolution Imaging Spectroradiometer (MODIS) at 250 m resolution (MOD250_LAI). The MOD250_LAI product uses a physical radiative transfer model which establishes a relationship between LAI, fraction of vegetation cover (FVC) and given patterns of surface reflectance, view-illumination conditions and optical properties of vegetation. In situ measurements of LAI and FVC made at 166 plots using hemispherical photography served for calibration of model parameters and validation of modelling results. Optical properties of vegetation cover, summarized by the light extinction coefficient, were computed at the local (pixel) level based on empirical models between ground-measured tree crown architecture at 85 sampling plots and spectral values in Landsat ETM+ bands. Influence of view-illumination conditions on optical properties of canopy was simulated by a view angle geometry model incorporating the solar zenith angle and the sensor viewing angle. The results revealed high compatibility of the produced MOD250_LAI data set with ground truth information and the 30 m resolution Landsat ETM+ LAI estimated using the similar algorithm. The produced MOD250_LAI was also compared with the global MODIS 1000-m LAI product (MOD15A2 LAI). Results show good consistency of the spatial distribution and temporal dynamics between the two LAI products. However, the results also showed that the annual LAI amplitude by the MOD15A2 product is significantly higher than by the MOD250_LAI. This higher amplitude is caused by a considerable underestimation of the tropical rainforest LAI by the MOD15A2 during the seasonal phases of low leaf production.  相似文献   

17.
Leaf to canopy upscaling approach affects the estimation of canopy traits   总被引:1,自引:0,他引:1  
In remote sensing applications, leaf traits are often upscaled to canopy level using sunlit leaf samples collected from the upper canopy. The implicit assumption is that the top of canopy foliage material dominates canopy reflectance and the variability in leaf traits across the canopy is very small. However, the effect of different approaches of upscaling leaf traits to canopy level on model performance and estimation accuracy remains poorly understood. This is especially important in short or sparse canopies where foliage material from the lower canopy potentially contributes to the canopy reflectance. The principal aim of this study is to examine the effect of different approaches when upscaling leaf traits to canopy level on model performance and estimation accuracy using spectral measurements (in-situ canopy hyperspectral and simulated Sentinel-2 data) in short woody vegetation. To achieve this, we measured foliar nitrogen (N), leaf mass per area (LMA), foliar chlorophyll and carbon together with leaf area index (LAI) at three vertical canopy layers (lower, middle and upper) along the plant stem in a controlled laboratory environment. We then upscaled the leaf traits to canopy level by multiplying leaf traits by LAI based on different combinations of the three canopy layers. Concurrently, in-situ canopy reflectance was measured using an ASD FieldSpec-3 Pro FR spectrometer, and the canopy traits were related to in-situ spectral measurements using partial least square regression (PLSR). The PLSR models were cross-validated based on repeated k-fold, and the normalized root mean square errors (nRMSEcv) obtained from each upscaling approach were compared using one-way analysis of variance (ANOVA) followed by Tukey’s post hoc test. Results of the study showed that leaf-to-canopy upscaling approaches that consider the contribution of leaf traits from the exposed upper canopy layer together with the shaded middle canopy layer yield significantly (p < 0.05) lower error (nRMSEcv < 0.2 for canopy N, LMA and carbon) as well as high explained variance (R2 > 0.71) for both in-situ hyperspectral and simulated Sentinel-2 data. The widely-used upscaling approach that considers only leaf traits from the upper illuminated canopy layer yielded a relatively high error (nRMSEcv>0.2) and lower explained variance (R2 < 0.71) for canopy N, LMA and carbon. In contrast, canopy chlorophyll upscaled based on leaf samples collected from the upper canopy and total canopy LAI exhibited a more accurate relationship with spectral measurements compared with other upscaling approaches. Results of this study demonstrate that leaf to canopy upscaling approaches have a profound effect on canopy traits estimation for both in-situ hyperspectral measurements and simulated Sentinel-2 data in short woody vegetation. These findings have implications for field sampling protocols of leaf traits measurement as well as upscaling leaf traits to canopy level especially in short and less foliated vegetation where leaves from the lower canopy contribute to the canopy reflectance.  相似文献   

18.
本文介绍了美国俄罗冈州西部黄松的叶面积指数(LAI)与小型航空光谱制图成像仪(CASI)获取的高光谱分辨率数据进行的相关分析。在试验场地上使用LAI-2000植物冠层分析仪测得8个LAI值(0.87—2.72)。对CASI数据进行一阶和二阶微分处理,以减少土壤背景光谱对森林光谱的影响。使用逐步回归分析方法探索LAI与CASI数据的关系。由回归分析产生多元线性方程和相应的拟合度(GOF)及估计LAI的标准误(SE)。结果说明光谱微分技术能够提高LAI和CASI数据的相关性,因而可以改善LAI的估计精度。如,对于单通道LAI预测的最高GOF值是0.681,SE是0.345,而经一阶和二阶光谱微分处理后,GOF被分别提高到0.904和0.898,SE被分别降低到0.189和0.195.  相似文献   

19.
水稻叶片不同光谱形式反演叶绿素含量的对比分析研究   总被引:7,自引:0,他引:7  
通过对常优1号和武粳15两个品种水稻叶片的反射率R、lg(1/R)、反射率一阶微分(FD)和反射率归一化(BN)等光谱形式的测量和计算,分析了叶片光谱不同变化形式与叶绿素含量的相关关系,建立了统计方程,并进行了比较与评价,同时,对反演方程的最佳波段选择进行了探讨。结果表明,叶绿素含量与反射率一阶微分光谱方程的相关性最强,而采用lg(1/R)的光谱形式能够提高遥感反演叶绿素含量的效果。经验证,两个水稻品种叶绿素含量的模拟值与实测值的复相关系数R2分别达到0.641和0. 818。  相似文献   

20.
The potential usefulness of spectral properties and vegetation indices in varietal discrimination of potato genotypes was studied in the field experiment. Spectral measurements were recorded in different bands in blue (450–520 nm), green (520–590 nm), red (620–680 nm) and infrared (770–860 nm) of the electromagnetic spectrum at different stages during crop growth period. A ground based hand held multiband radiometer (Model/041) was used for the purpose. The mean per cent green reflectance value among different genotypes was lowest in genotype MS/86-89, while it was observed highest in genotype JX-216. Significant difference among these genotypes was found at all growth stages except 6 week after planting. Consequent to variation in spectral reflectance the vegetation indices like, NDVI, RVI, TVI and DVI showed significant difference among genotypes at all growth stages except at 8th week after planting. The vegetation indices are good indicators of crop growth and condition. Similarly, fresh weight, dry weight, and leaf area index were also highest in MS/86-89, followed by KUFRI Bahar and KUFRI Sutlej while in case of leaf area index it was followed by Kufri Sutlej and Kufri Bahar. JX-23 was highest in chlorophyll content and tuber yield followed by MS/86-89 and JW-160, while lowest chlorophyll content was seen in MS/89-1095 and poorest tuber yield in MS/89-60. Most of the genotypes exhibited considerable variation in their spectral response and vegetation indices thereby indicating the possibility of their discrimination through remote sensing technique.  相似文献   

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

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