首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 234 毫秒
1.
Discrete wavelet transform was used to eliminate the noise in the charge-coupled device near-infrared (CCD-NIR) spectra of apple. The influence of three parameters (wavelet function, decomposition level, and threshold) on the predictive ability of the calibration model was investigated. The result showed that the db, sym, and bior wavelet families performed well, while the coif, dmey, and haar wavelets were not able to denoise effectively. The best decomposition level was 2. The threshold selection rules of the default, Birge-Massart, and Penalty had good denoising results, while SURE, Sqtwolog, Heuristic SURE, and Minimax set all detailed coefficients to zero due to their high threshold values. The best denoising result was obtained with the combination of the bior3.3 wavelet function, two levels of decomposition, default threshold selection rule, and the soft thresholding method. The optimal model of soluble solids content was constructed. The relative standard deviation of prediction decreased from 7.79 to 5.82% after wavelet denoising.  相似文献   

2.
The representation of complex soil variation on wavelet packet bases   总被引:2,自引:0,他引:2  
The discrete wavelet packet transform (DWPT) is an advanced wavelet technique that has distinct advantages over the discrete wavelet transform (DWT) previously used in soil science. Because the DWT divides the spatial frequencies of the data into non‐uniform intervals, it may fail to resolve features of the spatial variation that a simpler spectral analysis might identify under stationary conditions. However, the DWPT allows us to retain the advantages of the DWT (particularly the lack of stationarity assumptions) while achieving much better resolution in the spatial frequency domain. This is at the cost of poorer spatial resolution at the higher frequencies. However, by selecting a best wavelet packet basis from among the many possible, according to some appropriate criterion, we can find a compromise between frequency and spatial resolution that is best suited to the representation of our particular data set. In this paper, I describe a number of analyses based on the DWPT and apply them to two data sets on the soil. The advantages of the DWPT with best basis selection are clearly illustrated. In particular, non‐stationary behaviour of a significant periodic component of variation was identified that could not be effectively resolved by the DWT. Improved resolution in the frequency domain allowed aspects of the non‐stationary variation of soil to be resolved in a multiresolution analysis (MRA) adapted to the variation of the data.  相似文献   

3.
含风机噪声的蛋鸡声音信号去噪方法比较   总被引:6,自引:3,他引:3  
蛋鸡声音可用来评价蛋鸡本身的福利状况,然而蛋鸡舍中往往存在着低频风机噪声干扰蛋鸡声音信号的时频特征。为了优化含风机噪声的蛋鸡声音信号,以海兰褐蛋鸡为例,利用数字化声音采集平台,采集了不同的蛋鸡叫声和风机噪声。采用LabVIEW软件,分析了蛋鸡声音和风机噪声的时频特征,同时比较了滤波器和小波去噪方法在去除风机噪声方面的应用效果。结果表明,蛋鸡产蛋期间的声音频率范围为400~2 500 Hz,而风机噪声的频率在600 Hz以内。在信噪比为-20~10 dB蛋鸡声音环境中,无限脉冲响应滤波器滤波后的均方根误差要比有限脉冲响应滤波器滤波后的均方根误差小,说明无限脉冲响应滤波器具有更好的滤波效果,与其他小波阈值去噪方法相比,以史坦无偏似然估计为阈值的小波去噪方法在去噪后的均方根误差最小,表明这种方法的去噪效果更好。该研究可为蛋鸡舍中风机噪声环境下的蛋鸡声音识别提供参考。  相似文献   

4.
The plausibility of the assumption that soil variation can be treated as a realization of a random spatial process that is stationary in the variance can break down in various ways. It is possible to test the assumption using methods based on the wavelet transform. To date these approaches have been applied using the discrete wavelet transform. A drawback of this approach is that it uses a partition of the spatial frequencies represented in the data into intervals (scales) that are arbitrarily defined in advance and are not necessarily suitable for the representation of the variation of the data in question. A solution to this problem is to identify the best basis for the data from a wavelet packet library. An interesting question is whether the structure of this best basis is in itself informative about the plausibility of the stationarity assumption. In this paper, I show that this is indeed the case. The best basis for a stationary random variable from some packet library is the basis on the maximum dilation of the mother wavelet, which gives the finest resolution in the frequency domain. I propose the ratio of the entropy cost functional for this basis to that of the empirical best basis as a measure of evidence against the null hypothesis of stationarity in the variance. Critical values of this statistic may be obtained by Monte Carlo methods. I demonstrate the method using data on the clay content of soil on a transect in central England. The null hypothesis of stationarity in the variance may be rejected. Tests for the uniformity of variance can then be applied to wavelet packets in the best basis. The dominant local feature that is reflected in this behaviour is the unique pattern of variation in alluvium around a drainage channel that crosses the transect. This variation contrasts with that seen at most positions on the transect, variation that arises from a more or less regular pattern of boundaries between contrasting Jurassic strata.  相似文献   

5.
Developing wheat kernels of three New Zealand wheat cultivars at one to nine weeks postanthesis (WPA) were analyzed for kernel weight, salt‐soluble protein, salt‐insoluble protein, ascorbate oxidase (AOX), protein disulfide isomerase (PDI), ascorbic acid (AA), and dehydroascorbic acid (DHA). Kernel weight and salt‐soluble protein increased from one to five or six WPA then declined. The amount of salt‐insoluble protein per kernel rapidly increased from two to five or six WPA, then remained constant. AOX activity increased from one to two or three WPA then declined to almost zero at seven WPA. PDI activity increased from one to three WPA, stayed constant to six WPA, then declined about 25%. DHA increased from one to two or three WPA then declined to zero at seven WPA. AA increased from one to three or five WPA then declined to zero at seven WPA. When DHA, at similar concentrations (≈1 mM) to that found in one to six WPA kernels, was mixed with reduced‐gluten proteins, it was reduced to AA. The results are interpreted as a scheme showing that, during the formation of protein‐disulfide bonds, electrons are transferred from nascent protein thiols, ultimately to oxygen, through a series of electron carriers including PDI, AA‐DHA, and AOX.  相似文献   

6.
基于小波变换的Wiener滤波方法在农产品图像去噪中的应用   总被引:2,自引:2,他引:0  
农产品图像的去噪是农产品图像处理中最基本、最重要的工作之一。为了更有效地去除农产品图像中的噪声。受二维离散Wiener滤波器计算方法的启发,提出了一种基于小波变换的Wiener滤波方法。该方法采用小波变换和Wiener滤波相结合的方法,具有稀疏性、多分辨率、去相关性、选基灵活性和在MSE意义上对图像进行最优估计的优点。该方法首先对含噪农产品图像ano做第一次小波变换得到低频图像a1和水平、垂直和对角三方向的高频图像 hd1、vd1及dd1;其次对低频图像a1做Wiener滤波得到a1w,再对3个高频图像分别做Wiener滤波并合成得到g1w ;接着对低频的a1w和高频g1w做小波逆变换,得到滤波图像“a1w+g1w”。同时,考虑到噪声主要在高频部分,所以直接把低频的a1和高频g1w做小波逆变换,得到滤波图像“a1+g1w”。这是对含噪图像ano做第1次小波变换的情况,其第2次、第3次及第4次变换的情况与此类似。这样可以得到许多滤波图像,然后根据图像信噪比PSNR和视觉效果,最终确定去噪效果最好的农产品图像。该方法应用于红枣、小麦杂草等农产品图像的去噪中,结果PSNR为158.23(视觉效果清晰),好于邻域平均法(PSNR 为154.14)、中值滤波法(PSNR 为155.82)、数学形态学(PSNR为154.07,视觉效果偏黑)、高斯滤波法(PSNR为153.79,视觉效果太黑)、直接维纳去噪(PSNR为154.14)和小波去噪(PSNR为158.18)等多种方法。试验结果表明,基于小波变换的Wiener滤波方法应用于农产品图像去噪具有信噪比高、视觉效果好等优点;基于小波变换的Wiener滤波方法用于农产品图像去噪是有效的、可行的。  相似文献   

7.
Impact assessments on river systems of the combined effect of bed and suspended sediment loads from quarries are difficult to find. In this study, bed and suspended loads were measured to determine the impact of a 20‐ha limestone quarry on the river system of its 5,000‐ha steep, diverse land use/land cover but mostly forested catchment. A network of hydrologic and sediment monitoring instruments was deployed over the catchment during two separate study periods when sediment loadings were measured from captured storms. Results showed that the quarry stood to make a disproportionately large contribution to the catchment's estimated 2·1 Mg ha−1 yr−1 suspended sediment load. Large storm events contributed most of the loadings with five events supplying 92% of total loadings at the outlet. A paired method approach to compare suspended sediment loads between two subcatchments showed that during eight storm events, the quarry yielded between 2 and 49·2 Mg ha−1 per event, whereas the forest never yielded more than 0·1 Mg ha−1. Furthermore, the contribution of sediments from the quarry to bed load was more than 75% at a section located 1·2 km downstream. Future management activities to reduce sediment and bed loads, not only from this catchment but also from all others with similar land use/land covers, should focus on improving quarry operations. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

8.
Diffuse reflectance spectroscopy using visible (vis), near‐infrared (NIR) and mid‐infrared (mid‐IR) energy can be a powerful tool to assess and monitor soil quality and function. Mathematical pre‐processing techniques and multivariate calibrations are commonly used to develop spectroscopic models to predict soil properties. These models contain many predictor variables that are collinear and redundant by nature. Partial least squares regression (PLSR) is often used for their analysis. Wavelets can be used to smooth signals and to reduce large data sets to parsimonious representations for more efficient data storage, computation and transmission. Our aim was to investigate their potential for the analyses of soil diffuse reflectance spectra. Specifically we wished to: (i) show how wavelets can be used to represent the multi‐scale nature of soil diffuse reflectance spectra, (ii) produce parsimonious representations of the spectra using selected wavelet coefficients and (iii) improve the regression analysis for prediction of soil organic carbon (SOC) and clay content. We decomposed soil vis‐NIR and mid‐IR spectra using the discrete wavelet transform (DWT) using a Daubechies’s wavelet with two vanishing moments. A multiresolution analysis (MRA) revealed their multi‐scale nature. The MRA identified local features in the spectra that contain information on soil composition. We illustrated a technique for the selection of wavelet coefficients, which were used to produce parsimonious multivariate calibrations for SOC and clay content. Both vis‐NIR and mid‐IR data were reduced to less than 7% of their original size. The selected coefficients were also back‐transformed. Multivariate calibrations were performed by PLSR, multiple linear regression (MLR) and MLR with quadratic polynomials (MLR‐QP) using the spectra, all wavelet coefficients, the selected coefficients and their back transformations. Calibrations by MLR‐QP using the selected wavelet coefficients produced the best predictions of SOC and clay content. MLR‐QP accounted for any nonlinearity in the data. Transforming soil spectra into the wavelet domain and producing a smaller representation of the data improved the efficiency of the calibrations. The models were computed with reduced, parsimonious data sets using simpler regressions.  相似文献   

9.
针对风电场并网点电压故障引起的风机大规模脱网问题,提出了基于柔性形态算子和经验模态分解(empirical mode decomposition,EMD)去噪技术的电网电压故障检测方法。首先,利用EMD对采样信号进行时频自适应预处理,从而确定噪声主导模态;然后,通过柔性形态学变换加阈值输出,有效放大信号奇异点,避免了因电网电压信号周期性变化和噪声引起的背景梯度对检测结果的影响,实现故障定位检测。通过对不同噪声强度的电压暂降故障信号进行检测对比分析发现,随着信号信噪比下降,标准形态学方法的检测误差进一步增大,当信噪比达到25db时,甚至出现了误检现象,而柔性形态EMD检测方法仍然可以有效检测故障扰动的起止时间,表明该方法与标准形态学和小波阈值方法相比,在简化运算过程的同时可以获得更高的检测精度。最后,对某风电场并网点故障电压的分析结果与实测数据的一致性,验证了该方法可以有效检测电网电压的瞬态故障信息,从而为风电场无功补偿装置的投切控制提供了依据。  相似文献   

10.
狭叶锦鸡儿叶切片显微图像在获取过程中不可避免的受到噪声污染,会对后续处理造成不良影响。针对现有噪声类型未知,去噪算法存在速度慢、效果不理想等问题,该文提出图像噪声类型估计-强度估计-去噪这一处理过程,实现对狭叶锦鸡儿叶切片显微图像降噪目的。首先采用平滑区直方图重构和拟合法确定噪声类型;然后在此基础上,应用基于图像块的SVD(singular valuable decomposition,SVD)域图像噪声强度估计法对噪声标准差进行估计;最后在确定噪声类型和强度基础上,采用几何均值滤波(geometric mean filtering,GMF)和三维块匹配滤波(block-matching and 3-D filtering,BM3D)对图像进行联合去噪。试验结果表明:该文噪声类型估计法估计出切片图像噪声类型为加性高斯噪声,高斯函数对随机选取的15幅狭叶锦鸡儿叶切片图像平滑区域直方图数据点拟合优度2R均值为0.996,平均均方根误差RMSE(root mean squared error,RMSE)为0.144 6;采用该文噪声强度估计法估计出的切片图像噪声标准差???[2.5,4.0],处理标准差较小噪声,该文算法处理精度、运行速度和稳定性等方面存在明显优势;GMF-BM3D算法在较好去除图像噪声同时,极大的保留了图像纹理、边缘和细节等信息,同时极大的提高了算法运行速度,处理后的图像BRISQUE(blind/referenceless image spatial quality evaluator,BRISQUE)值为10左右,相当于原图BRISQUE值的1/2左右。与传统BM3D算法相比,去噪效果相当,但耗时约相当于传统BM3D算法的1/9。与小波去噪算法(wavelet threshold,WT)算法相比,虽速度相对较慢,但去噪后图像BRISQUE值比使用WT法低4左右。因此,该算法较好实现了对狭叶锦鸡儿叶切片图像准确降噪,为其后续处理提供了可靠技术支持。  相似文献   

11.
In an earlier article, an intuitively appealing method for estimating the number of true null hypotheses in a multiple test situation was proposed. That article presented an iterative algorithm that relies on a histogram of observed p values to obtain the estimator. We characterize the limit of that iterative algorithm and show that the estimator can be computed directly without iteration. We compare the performance of the histogram-based estimator with other procedures for estimating the number of true null hypotheses from a collection of observed p values and find that the histogram-based estimator performs well in settings similar to those encountered in microarray data analysis. We demonstrate the approach using p values from a large microarray experiment aimed at uncovering molecular mechanisms of barley resistance to a fungal pathogen.  相似文献   

12.
离散小波变换和BP神经网络识别玉米种子纯度   总被引:3,自引:3,他引:0  
摘要:为快速有效地识别玉米种子纯度,针对玉米种子图像特征,对其图像处理方法和分类算法进行研究,提出一种基于离散小波变换和BP神经网络玉米种子纯度识别算法。该方法首先提取玉米种子冠部核心区域的RGB颜色模型特征参数,然后对三个色彩分量分别进行二层离散小波变换,提取各频带区域均值作为BP神经网络的输入样本,玉米种子的纯度分类作为神经网络的输出样本。实验结果表明该方法可准确识别玉米纯度并分类,正确识别率达94.5%。  相似文献   

13.
A model for predicting lime requirement (Bailey et al. , 1989) was tested using data from independent liming experiments on soils from Northern Ireland and New Zealand. Lime requirements predicted by the model were also compared with those obtained by the methods currently used by the Department of Agriculture for Northern Ireland and by the Ministry of Agriculture, Fisheries and Food in England and Wales.
The model predicted lime requirement precisely for the independent set of soils from Northern Ireland, and was as reliable at predicting lime requirement for New Zealand soils as the best of seven other quick-test methods. The model was more reliable at predicting lime requirement for soils from Northern Ireland than the method currently used by the Department of Agriculture for Northern Ireland, and appeared to be as reliable as the method used by the Ministry of Agriculture, Fisheries and Food in England and Wales.  相似文献   

14.
柑橘叶片叶绿素含量高光谱无损检测模型   总被引:13,自引:5,他引:13  
针对柑橘叶片叶绿素含量的传统化学检测,不仅耗时长且损伤柑橘叶片,还依赖检测者实操技术,无法集成于精细农业中变量喷施农机具的诸多弊端,该文探讨快速无损检测柑橘叶片叶绿素含量方法。以117棵园栽萝岗甜橙树为研究对象,选用ASD Field Spec 3光谱仪对萌芽期、稳果期、壮果促梢期、采果期共4个生长时期的柑橘叶片进行高光谱反射率采集,并同步采用分光光度法测得叶片的叶绿素含量;以原始光谱及其变换形式作为模型输入矢量,分别在主成分分析(principle component analysis,PCA)降维的基础上利用支持向量机回归(support vector regression,SVR)算法和在小波去噪的基础上利用偏最小二乘回归(partial least square regression,PLSR)算法对柑橘叶片叶绿素含量进行建模预测,全生长期整体建模的校正集和验证集最佳模型决定系数R2分别为0.8713和0.8670,均方根误差RMSE(root-mean-square error)分别为0.1517和0.1544,试验结果表明,高光谱可快速无损地对柑橘叶片叶绿素含量进行精确的定量检测,为柑橘不同生长期的营养监测提供理论依据。  相似文献   

15.
为了研究采用荧光光谱技术对生菜农药残留快速无损定性鉴别的可行性,该文通过采集180个生菜样品(3个浓度农药残留生菜,每个浓度农残生菜样本数为60,其中农药与水配比为1:500、1:1000、1:1200,即重度超标、轻微超标、标准农残)的荧光发射光谱,结合Savitzky-Golay(SG)、标准正态变量变换(standard normalized variable,SNV)、标准正态变量变换结合去趋势(standard normalized variable detrending,SNV detrending)、SG与SNV算法组合(SG-SNV)、SG与SNV detrending算法组合(SG-SNV detrending)对提取的荧光光谱进行预处理,基于全光谱、荧光特征峰值、小波特征建立支持向量机(Support Vector Machine, SVM)分类模型。其中,小波特征通过小波变换对原始光谱以及预处理后光谱进行特征选择获取,分别采用 db4、db6、sym5、sym7作为小波基函数。试验结果表明:基于小波特征、荧光特征峰值建立的 SVM 模型预测集识别率要高于基于全光谱建立的 SVM 模型。以 sym5作为小波基函数,基于 SG-SNV detrending预处理光谱选择的小波特征建立的SVM模型取得最优的预测集识别率93.33%,最佳小波分解层数为4。结果表明应用荧光光谱技术对生菜农药残留鉴别是可行的,为生菜农药残留快速、无损检测分析提供了参考。  相似文献   

16.
应用近红外光谱结合化学计量学方法对蜂蜜产地进行了判别分析。kennard-Stone法划分训练集和预测集。光谱用一阶导数加自归一化预处理后,再用小波变换(WT)进行压缩和滤噪。结合滤波后光谱信息,分别用径向基神经网络(RBFNN)和偏最小二乘-线性判别分析(PLS-LDA)建立了苹果蜜产地和油菜蜜产地判别模型。对不同小波基和分解尺度进行了详细讨论。对苹果蜜,WT-RBFNN模型和WT-PLS-LDA模型都是小波基为db1、分解尺度为2时的预测精度最好,都为96.2%。对油菜蜜:WT-RBFNN模型在小波基为db4和分解尺度为1时,预测精度最好;WT-PLS-LDA模型在小波基为db9、分解尺度也为1时,预测精度最好,为90.5%;预测精度WT-PLS-LDA模型优于WT-RBFNN模型。研究表明:WT结合线性的PLS-LDA建模比WT结合非线性的RBFNN建模更适于蜂蜜产地鉴别;近红外光谱结合WT-PLS-LDA可实现对蜂蜜产地的快速无损检测,为蜂蜜产地鉴别提供了一种新方法。  相似文献   

17.
小波变换耦合CARS算法提高土壤水分含量高光谱反演精度   总被引:4,自引:3,他引:1  
为实现干旱地区土壤水分含量(soil moisture content,SMC)的快速监测,该文以渭干河-库车河绿洲为靶区,采用小波变换(wavelet transform,WT)对反射光谱进行1~8层小波分解,通过相关性分析确定最大分解层数,再通过竞争性自适应重加权(competitive adaptive reweighted sampling,CARS)滤除冗余变量,筛选出与SMC相关性较好的波长变量,并叠加各层特征光谱的优选波长变量作为最优变量集,用偏最小二乘回归(partial least squares regression,PLSR)构建土壤水分含量预测模型并进行分析.结果显示:1)小波分解过程中,土壤反射率与SMC的相关性不断增强,到小波变换第6层分解(L6)处达到最高,因此小波变换最大分解层数为6层分解;2)通过对土样进行WT-CARS耦合算法筛选出变量,得出的最优变量集包括400~500、1 320~1 461、1 851~1 961、2 125~2 268 nm区域之间共131个波长变量;3)相对于全波段预测模型,各层特征光谱的CARS优选变量预测模型的精度均高,并且基于最优变量集的预测模型的精度最高,该模型的建模集均方根误差0.021、建模集决定系数0.721、预测集均方根误差0.028、预测集决定系数0.924、相对分析误差2.607.说明WT-CARS耦合算法使其在建立模型时尽可能少地损失光谱细节、较为彻底的去除噪声,同时还能对无信息变量进行有效去除,为该研究区SMC的预测提供新的思路.  相似文献   

18.
W.H. Lyford 《Geoderma》1974,11(3):195-208
Narrow bands of soils, 5 to 30 m wide, are common in southern New England immediately bordering wet areas, along drainageways, and where soil parent materials overlap or grade into one another. On the basis of detailed observations along 16- and 19-km trenches across representative landscapes, soil bodies less than 30 m wide occupy about 20% of the southern New England area. Yet on recently published soil survey maps such narrow soils seldom are delineated because of scale limitations, and the delineated areas less than 30 m wide occupy only about 1% of mapped areas. Narrow soils in forested areas are often more important than their width would indicate because of their influence on site quality, land use and natural-vegetation distribution.  相似文献   

19.
Populations of New England cottontails (Sylvilagus transitionalis) have declined substantially in recent decades in response to habitat loss and fragmentation. Among some remnant populations, cottontails occupy small patches of thicket habitat where they experience high mortality rates as a consequence of limited food during winter. This limitation causes rabbits to forage away from cover where they are exposed to predators. Although conservation efforts are emerging to reverse the decline of New England cottontails, most are directed toward improving long-term viability by increasing the abundance of suitable habitats. Such efforts do little to improve the short-term survival of remaining cottontails. To address this immediate need, we evaluated the use of supplemental food as an approach to improve overwinter survival rates. We speculated that by positioning feeders in close proximity to escape cover, rabbits would be less vulnerable to predators. We evaluated this approach using eastern cottontails (Sylvilagus floridanus) as a research surrogate because this species is readily available and has similar habitat requirements to New England cottontails. Transmitter-equipped eastern cottontails were randomly assigned to either a fed or unfed group. Remotely-triggered cameras were also used to gauge use of feeders by cottontails and visits by other species. Winter survival rates were substantially greater for fed rabbits (70%) than for unfed rabbits (32%). Cameras revealed that rabbits were the most frequent consumer and that there was only limited carnivore activity near feeders. We conclude that supplemental feeding may improve survival of remaining New England cottontails as efforts to increase habitat availability continue to develop.  相似文献   

20.
农田遥感图像在采集过程中会受到噪声影响,为得到准确的农田遥感图像数据,应对获取的农田遥感图像进行去噪预处理。农田遥感图像中的纹理承载了重要信息,在图像降噪的同时保持或增强图像纹理具有重要意义。由于纹理和噪声一样,在频域表现为高频信号,以分解和重构算法为基础的常见滤波(含小波变换)方法在降噪的同时,也会造成纹理清晰度的下降。该文结合农田遥感图像纹理呈现出来的直线特性,将剪切波(Shearlet)和变分理论相结合,提出了一种新的遥感农田图像保纹理降噪方法。该方法首先对较大的遥感图像分块进行shearlet变换,在降噪的同时识别不同图块图像的纹理含量;对细小纹理含量较少的平滑区域,采用保边降噪变分模型去除shearlet变换带来的人工伪影。为避免子图块边界带来的边界效应,该文基于中心仿射变换理论提出了一种新的图像延拓方法,有效提高了图像降噪的效果。试验结果表明,该文算法去噪后的峰值信噪比(peak signal to noise ratio,PSNR)平均值比全变分模型去噪算法大1 d B,该文算法去噪后的PSNR平均比曲线波去噪算法大2 d B。同基于Symmlet小波的Shearlet算法相比,该文算法处理后农田遥感图像中伪影减少,在高斯噪声标准偏差σ为10、20和30 d B时,峰值信噪比PSNR分别提高了13.99%、9.69%和7.75%。  相似文献   

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

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