首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Land surface temperature (LST) is an important factor in global change studies, heat balance and as control for climate change. A comparative study of LST over parts of the Singhbhum Shear Zone in India was undertaken using various emissivity and temperature retrieval algorithms applied on visible and near infrared (VNIR), and thermal infrared (TIR) bands of high resolution Landsat-7 ETM+ imagery. LST results obtained from satellite data of October 26, 2001 and November 2, 2001 through various algorithms were validated with ground measurements collected during satellite overpass. In addition, LST products of MODIS and ASTER were compared with Landsat-7 ETM+ and ground truth data to explore the possibility of using multi-sensor approach in LST monitoring. An image-based dark object subtraction (DOS3) algorithm, which is yet to be tested for LST retrieval, was applied on VNIR bands to obtain atmospheric corrected surface reflectance images. Normalized difference vegetation index (NDVI) was estimated from VNIR reflectance image. Various surface emissivity retrieval algorithms based on NDVI and vegetation proportion were applied to ascertain emissivities of the various land cover categories in the study area in the spectral range of 10.4–12.5 μm. A minimum emissivity value of about 0.95 was observed over the reflective rock body with a maximum of about 0.99 over dense forest. A strong correlation was established between Landsat ETM+ reflectance band 3 and emissivity. Single channel based algorithms were adopted for surface radiance and brightness temperature. Finally, emissivity correction was applied on ‘brightness temperature’ to obtain LST. Estimated LST values obtained from various algorithms were compared with field ground measurements for different land cover categories. LST values obtained after using Valor’s emissivity and single channel equations were best correlated with ground truth temperature. Minimum LST is observed over dense forest as about 26 °C and maximum LST is observed over rock body of about 38 °C. The estimated LST showed that rock bodies, bare soils and built-up areas exhibit higher surface temperatures, while water bodies, agricultural croplands and dense vegetations have lower surface temperatures during the daytime. The accuracy of the estimated LST was within ±2 °C. LST comparison of ASTER and MODIS with Landsat has a maximum difference of 2 °C. Strong correlation was found between LST and spectral radiance of band 6 of Landsat-7 ETM+. Result corroborates the fact that surface temperatures over land use/land cover types are greatly influenced by the amount of vegetation present.  相似文献   

2.
The present study is an assessment and identification of urban heat island (UHI) in the environment of one of the fastest urbanizing city of India, Delhi Metropolis, employing satellite image of ASTER and Landsat 7 ETM+ in the thermal infrared region 3–14 μm. Temporal (2001 and 2005) ASTER datasets were used to analyze the spatial structure of the thermal urban environment subsequently urban heat island (UHI) in relation to the urban surface characteristics and land use/land cover (LULC). The study involves derivation of parameters governing the surface heat fluxes, constructing statistics of ASTER thermal infrared images along with validation through intensive in situ measurements. The average images reveal spatial and temporal variations of land surface temperature (LST) of night-time and distinct microclimatic patterns. Central Business District (CBD) of Delhi, (Connaught Place, a high density built up area), and commercial/industrial areas display heat islands condition with a temperature greater than 4 °C compared to the suburbs. The small increase in surface temperature at city level is mainly attributed to cumulative impact of human activities, changes in LULC pattern and vegetation density. In this study the methodology takes into account spatially-relative surface temperatures and impervious surface fraction value to measure surface UHI intensity between the urban land cover and rural surroundings. Both the spatial and temporal variation in surface temperature associated with impervious surface area (ISA) has been evaluated to assess the effect of urbanization on the local climate.  相似文献   

3.
The eastern part of the Rich area consists of the massive Paleozoic and Meso-Cenozoic cover formations that present the geodynamic development of the study area, where is characterized by various carbonate facies of Jurassic age. The geographical characteristic of the study area leaves the zone difficult to map by conventional methods. The objective of this work focuses on the mapping of the constituent lithological units of the study area using multispectral data of Landsat OLI, ASTER, and Sentinel 2A MSI. The processing of these data is based on a precise methodology that distinguishs and highlights the limits of the different lithological units that have an approximate similarity of spectral signature. Three techniques were used to enhance the image including Principal Component Analysis (PCA), Minimum Noise Fraction (MNF), and Independent Component Analysis (ICA). Lithological mapping was performed using two types of supervised classification : Maximum likelihood classifier (MLC) and Support Vector Machine (SVM).The results of processing data show the effectiveness of Sentinel 2A data in mapping of lithological units than the ASTER and Landsat OLI data. The classification evaluation of two methods of the Sentinel 2A MSI image showed that the SVM method give a better classification with an overall accuracy of 93,93% and a Kappa coefficient of 0.93, while the MLC method present an overall accuracy of 82,86% and a Kappa coefficient of 0.80. The results of mapping obtained show a good correlation with the geological map of the study area as well as the efficiency of remote sensing in identification of different lithological units in the Central High Atlas.  相似文献   

4.
Object-based rice mapping using time-series and phenological data   总被引:1,自引:0,他引:1  
Remote sensing techniques are often used in mapping rice, but high quality time-series remote sensing data are difficult to obtain due to the cloudy weather of rice growing areas and long satellite revisit interval. As such, rice mapping is usually based on mono-temporal Landsat TM/ETM+ data, which have large uncertainties due to the spectral similarity of different vegetation types. Moreover, conventional pixel-based classification method is unable to meet the required accuracy for rice mapping. Therefore, this study proposes a new strategy for mapping rice in cloud-prone areas using fused data of Landsat-8 OLI time-series and phenological parameters, based on the object-based method. We determine the critical growth stages of paddy rice from observed phenological data and MODIS-NDVI time-series data. The spatial and temporal adaptive reflectance fusion model (STARFM) is used to blend the MODIS and Landsat data to obtain a multi-temporal Landsat-like dataset for classification. Finally, an object-oriented algorithm is used to extract rice paddies from the Landsat-like, time-series dataset. The validation experiments show that the proposed method can provide high accuracy rice mapping, with an overall accuracy of 92.38% and a kappa coefficient of 0.85.  相似文献   

5.
Spectral transformation methods, including correlation coefficient (CC) and Optimum Index Factor (OIF), band ratio (BR) and principal component analysis (PCA) were applied to ASTER and Landsat TM bands for lithological mapping of Soghan ophiolitic complex in south of Iran. The results indicated that the methods used evidently showed superior outputs for detecting lithological units in ophiolitic complexes. CC and OIF methods were used to establish enhanced Red–Green–Blue (RGB) color combination bands for discriminating lithological units. A specialized band ratio (4/1, 4/5, 4/7 in RGB) was developed using ASTER bands to differentiate lithological units in ophiolitic complexes. The band ratio effectively detected serpentinite dunite as host rock of chromite ore deposits from surrounding lithological units in the study area. Principal component images derived from first three bands of ASTER and Landsat TM produced well results for lithological mapping applications. ASTER bands contain improved spectral characteristics and higher spatial resolution for detecting serpentinite dunite in ophiolitic complexes. The developed approach used in this study offers great potential for lithological mapping using ASTER and Landsat TM bands, which contributes in economic geology for prospecting chromite ore deposits associated with ophiolitic complexes.  相似文献   

6.
The concerns over land use/land cover (LULC) change have emerged on the global stage due to the realisation that changes occurring on the land surface also influence climate, ecosystem and its services. As a result, the importance of accurate mapping of LULC and its changes over time is on the increase. Landsat satellite is a major data source for regional to global LULC analysis. The main objective of this study focuses on the comparison of three classification tools for Landsat images, which are maximum likelihood classification (MLC), support vector machine and artificial neural network (ANN), in order to select the best method among them. The classifiers algorithms are well optimized for the gamma, penalty, degree of polynomial in case of SVM, while for ANN minimum output activation threshold and RMSE are taken into account. The overall analysis shows that the ANN is superior to the kernel based SVM (linear, radial based, sigmoid and polynomial) and MLC. The best tool (ANN) is then applied on detecting the LULC change over part of Walnut Creek, Iowa. The change analysis of the multi temporal images indicates an increase in urban areas and a major shift in the agricultural practices.  相似文献   

7.
The use of satellite images for mineral exploration has been very successful in pointing out the presence of minerals such as smectite and kaolinite which are important in the identification of hydrothermal alterations. Shortwave infrared (SWIR) bands from Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) with the wavelength of ASTER SWIR bands between 1.65 and 2.43 μm has a good potential for mapping a hydrothermal alteration minerals such as alunite, pyrophyllite, kaolinite, illite–muscovite–sericite, and carbonate. In this range, hydroxide minerals which have been produced by hydrothermal alteration exhibit good absorption compared to shorter or longer wavelengths. In this research which aims to remove atmospheric and topographic effects from ASTER SWIR data, the authors used the log-residual method (LRM) with the minimum noise fraction (MNF) transformation to create a pixel purity index (PPI) which was used to extract the most spectrally pure pixels from multispectral images. Spectral analyses of the clay mineralogy of the study area (east Zanjan, in northern Iran) were obtained by matching the unknown spectra of the purest pixels to the U.S. Geological Survey (USGS) mineral library. Three methods, spectral feature fitting (SFF), spectral angle mapping (SAM), and binary encoding (BE) were used to generate a score between 0 and 1, where a value of 1 indicates a perfect match showing the exact mineral type. In this way, it was possible to identify certain mineral classes, including chlorite, carbonate, calcite–dolomite–magnesite, kaolinite–smectite, alunite, and illite. In this research, two main propylitic and phyllic–argillic zones could be separated using their compositions of these minerals. These two alteration zones are important for porphyry copper deposits and gold mineralization in this part of Iran.  相似文献   

8.
This study investigates the application of spectral image processing methods to ASTER data for mapping hydrothermal alteration zones associated with porphyry copper mineralization and related host rock. The study area is located in the southeastern segment of the Urumieh–Dokhtar Volcanic Belt of Iran. This area has been selected because it is a potential zone for exploration of new porphyry copper deposits. Spectral transform approaches, namely principal component analysis, band ratio and minimum noise fraction were used for mapping hydrothermally altered rocks and lithological units at regional scale. Spectral mapping methods, including spectral angle mapper, linear spectral unmixing, matched filtering and mixture tuned matched filtering were applied to differentiate hydrothermal alteration zones associated with porphyry copper mineralization such as phyllic, argillic and propylitic mineral assemblages.  相似文献   

9.
Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) and Phased Array L-band Synthetic Aperture Radar (PALSAR) images covering the Um Rus area in the Central Eastern Desert of Egypt were evaluated for mapping geologic structure, lithology, and gold-related alteration zones. The study area is covered by Pan-African basement rocks including gabbro and granodiorite intruded into a variable mixture of metavolcanics and metasediments. The first three principal component analyses (PCA1, PCA2, PCA3) in a Red-Green-Blue (RGB) of the visible through shortwave-infrared (VNIR + SWIR) ASTER bands enabled the discrimination between lithological units. The results show that ASTER band ratios ((2 + 4)/3, (5 + 7)/6, (7 + 9)/8) in RGB identifies the lithological units and discriminates the granodiorite very well from the adjacent rock units.The granodiorites are dissected by gold-bearing quartz veins surrounded by alteration zones. The microscopic examination of samples collected from the alteration zones shows sericitic and argillic alteration zones. The Spectral Angle Mapper (SAM) and Spectral Information Divergence (SID) supervised classification methods were applied using the reference spectra of the USGS spectral library. The results show that these classification methods are capable of mapping the alteration zones as indicated by field verification work. The PALSAR image was enhanced for fracture mapping using the second moment co-occurrence filter. Overlying extracted faults and alteration zone classification images show that the N30E and N-S fractures represent potential zones for gold exploration. It is concluded that the proposed methods can be used as a powerful tool for ore deposit exploration.  相似文献   

10.
This paper proposed a quantitative method of land use change trajectory, which means the succession among different land use types across time, to examine the effects of urbanization on an urban heat island (UHI). To accomplish this, multi-temporal images from Landsat 5 Thematic Mapper (TM) and Landsat 7 Enhanced Thematic Mapper Plus (ETM+) of Xiamen City in China from 1987 to 2007 were selected. First, the land use change trajectories were extracted based on the classified images from different years. Then the land surface temperatures (LST) were retrieved and the magnitudes of the UHI were evaluated using the UHI intensity (UHII) indicator. Finally, the indices of the contribution to UHI intensity (CUHII) were constructed and calculated to quantify the effects of each land use change trajectory on the UHI during urbanization. The results demonstrated that the land use change trajectories and CUHII are effective and useful in quantifying the effects of urbanization on UHI. In Xiamen City, a total of 2218 land use change trajectories were identified and 530 of them were the existing urban or urbanization trajectories. The UHII presents a trend of continuous increase from 0.83 °C in 1987 to 2.14 °C in 2007. With respect to the effects of urban growth on UHI, the contribution of existing urban area to UHI decreased during urbanization. Prior to 2007, the existing urban area of trajectory NO. 44444 had the most significant effect on UHI with the greatest CUHII, while the value has decreased from 55.00% in 1987 to 13.03% in 2007 because of the addition of new urbanized area. In 2007, the greatest CUHII was replaced by a trajectory from farmland to built-up area (NO. 22224) with the CUHII of 21.98%, followed by the existing urban area of trajectory NO. 44444 with the CUHII of 13.03%. These results provide not only a new methodology to assess the environmental effects of urbanization, but also decision-supports for the planning and management of cities.  相似文献   

11.
This paper discusses an approach for river mapping and flood evaluation based on multi-temporal time series analysis of satellite images utilizing pixel spectral information for image classification and region-based segmentation for extracting water-covered regions. Analysis of MODIS satellite images is applied in three stages: before flood, during flood and after flood. Water regions are extracted from the MODIS images using image classification (based on spectral information) and image segmentation (based on spatial information). Multi-temporal MODIS images from “normal” (non-flood) and flood time-periods are processed in two steps. In the first step, image classifiers such as Support Vector Machines (SVM) and Artificial Neural Networks (ANN) separate the image pixels into water and non-water groups based on their spectral features. The classified image is then segmented using spatial features of the water pixels to remove the misclassified water. From the results obtained, we evaluate the performance of the method and conclude that the use of image classification (SVM and ANN) and region-based image segmentation is an accurate and reliable approach for the extraction of water-covered regions.  相似文献   

12.
The purpose of this study was twofold: to develop a methodology for the estimation of land surface temperature for non-urban areas and to analyze the sensitivity of the methodology. The key element of the methodology was the development of emissivity maps based on CORINE Land Cover and the ASTER spectral library. Land surface temperatures were derived from NOAA/AVHRR data and the methodology was applied at a national scale in Greece, with emphasis given to non-urban areas. A sensitivity analysis was performed in order to determine the variables that mainly affect the estimation of land surface temperature. A varying propagation error was identified depending on the temperature and humidity of the atmosphere, as well as the land cover type. The methodology was applied to a series of 25 AVHRR images and the results were compared to in-situ measurements from representative stations.  相似文献   

13.
COVID-19 pandemic has had a major impact on our society, environment and public health, in both positive and negative ways. The main aim of this study is to monitor the effect of COVID-19 pandemic lockdowns on urban cooling. To do so, satellite images of Landsat 8 for Milan and Rome in Italy, and Wuhan in China were used to look at pre-lockdown and during the lockdown. First, the surface biophysical characteristics for the pre-lockdown and within-lockdown dates of COVID-19 were calculated. Then, the land surface temperature (LST) retrieved from Landsat thermal data was normalized based on cold pixels LST and statistical parameters of normalized LST (NLST) were calculated. Thereafter, the correlation coefficient (r) between the NLST and index-based built-up index (IBI) was estimated. Finally, the surface urban heat island intensity (SUHII) of different cities on the lockdown and pre-lockdown periods was compared with each other. The mean NLST of built-up lands in Milan (from 7.71 °C to 2.32 °C), Rome (from 5.05 °C to 3.54 °C) and Wuhan (from 3.57 °C to 1.77 °C) decreased during the lockdown dates compared to pre-lockdown dates. The r (absolute value) between NLST and IBI for Milan, Rome and Wuhan decreased from 0.43, 0.41 and 0.16 in the pre-lockdown dates to 0.25, 0.24, and 0.12 during lockdown dates respectively, which shows a large decrease for all cities. Analysis of SUHI for these cities showed that SUHII during the lockdown dates compared to pre-lockdown dates decreased by 0.89 °C, 1.78 °C, and 1.07 °C respectively. The results indicated a high and substantial impact of anthropogenic activities and anthropogenic heat flux (AHF) on the SUHI due to the substantial reduction of huge anthropogenic pressure in cities. Our conclusions draw attention to the contribution of COVID-19 lockdowns (reducing the anthropogenic activities) to creating cooler cities.  相似文献   

14.
Variable Endmember Constrained Least Square (VECLS) technique is proposed to account endmember variability in the linear mixture model by incorporating the variance for each class, the signals of which varies from pixel to pixel due to change in urban land cover (LC) structures. VECLS is first tested with a computer simulated three class endmember considering four bands having small, medium and large variability with three different spatial resolutions. The technique is next validated with real datasets of IKONOS, Landsat ETM+ and MODIS. The results show that correlation between actual and estimated proportion is higher by an average of 0.25 for the artificial datasets compared to a situation where variability is not considered. With IKONOS, Landsat ETM+ and MODIS data, the average correlation increased by 0.15 for 2 and 3 classes and by 0.19 for 4 classes, when compared to single endmember per class.  相似文献   

15.
Remote sensing satellite data plays a vital role and capable in detecting minerals and discriminating rock types for explorations of mineral resources and geological studies. Study of spectral absorption characters of remotely sensed data are under consideration by the exploration and mining companies, and demonstrating the spectral absorption characters of carbonates on the cost-effective multispectral image (rather than the hyperspectral, Lidar image) for easy understanding of all geologists and exploration communities of carbonates is very much important. The present work is an integrated study and an outcome of recently published works on the economic important carbonate rocks, includes limestone, marl, listwaenites and carbonatites occurred in parts of the Sultanate of Oman. It demonstrates the spectral sensitivity of such rocks for simple interpretation over satellite data and describes and distinguishes them based on the absorptions of carbonate minerals in the spectral bands of advanced spaceborne thermal emission and reflection radiometer (ASTER) for mapping and exploration studies. The study results that the ASTER spectral band 8 discriminates the carbonate rocks due to the presence of predominantly occurred carbonate minerals; the ASTER band 5 distinguishes the limestones and marls (more hydroxyl clay minerals) from listwaenite (hydrothermally altered rock) due to the presence of altered minerals and the ASTER band 4 detects carbonatites (ultramafic intrusive alkaline rocks) which contain relatively more silicates. The study on the intensity of the total absorptions against the reflections of these rocks shows that the limestones and marls have low intensity in absorptions (and high reflection values) due to the presence of carbonate minerals (calcite and dolomite) occurred in different proportions. The listwaenites and carbonatites have high intensity of absorptions (low reflection values) due to the occurrence of Mn-oxide in listwaenites and carbonates in carbonatites apart the influence of major carbonate minerals that occurred predominantly in these rocks. The study of ASTER thermal infrared (TIR) spectral bands distinguished the marls have low emissivity of energy due to the presence of hydroxyl bearing alumina-silicate minerals from the other rocks such as limestones, listwaenites and carbonatites which have high emissivity due to the absence of hydroxyl bearing alumina-silicate minerals and the presence of carbonate minerals and carbonates.  相似文献   

16.
The aim of this study was to identify landslide-related factors using only remotely sensed data and to present landslide susceptibility maps using a geographic information system, data-mining models, an artificial neural network (ANN), and an adaptive neuro-fuzzy interface system (ANFIS). Landslide-related factors were identified in Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) satellite imagery. The slope, aspect, and curvature of topographic features were calculated from a digital elevation model that was made using the ASTER imagery. Lineaments, land-cover, and normalized difference vegetative index layers were also extracted from the imagery. Landslide-susceptible areas were analyzed and mapped based on occurrence factors using the ANN and ANFIS. The generalized bell-shaped built-in membership function of the ANFIS was applied to landslide susceptibility mapping. Analytical results were validated using landslide test location data. In the validation results, the ANN model showed 80.42% prediction accuracy and the ANFIS model showed 86.55% prediction accuracy. These results suggest that the ANFIS model has a better performance than does the ANN in predicting landslide susceptibility.  相似文献   

17.
The structural transformation of Polish agriculture after 1989 has been leading to significant changes in land use. As a result a large part of formerly ploughed fields lies abandoned and they occur across considerable variability in soil typological units. Accurate information about soil conditions within the abandoned fields facilitates proper management in the new socio-economic situation. Such information can be collected using satellite images since the structure and condition of the vegetation growing on the abandoned fields reflects soil properties. The objective of this study is to evaluate the relationships between physical and chemical attributes of soil within the abandoned fields and spectral reflectance patterns recorded by ASTER sensors onboard Terra Satellite. Soil samples were collected at five abandoned fields which have not been ploughed since 2000 and analyzed in the laboratory to determine their physical and chemical properties. Nine ASTER nearly cloud-free pictures were used for this study in order to derive the remote-sensing attributes of the abandoned fields. In order to evaluate the relationships between soil fertility and remotely sensed data, partial least-squares (PLS) and a multiple linear regression (MLR) analysis between these two datasets were carried out. In the regression analysis, only soil TEB (total exchangeable bases) stock in the whole profile displayed the highest correlation with remotely sensed data acquired in April and May and the best predictors were NDVI and LSWI vegetation indices.  相似文献   

18.
Precise glacier information is important for assessing climate change in remote mountain areas. To obtain more accurate glacier mapping, rough set theory, which can deal with vague and uncertainty information, was introduced to obtain optimal knowledge rules for glacier mapping. Optical images, thermal infrared band data, texture information and morphometric parameters were combined to build a decision table used in our proposed rough set theory method. After discretizing the real value attributes, decision rules were calculated through the decision rule generation algorithm for glacier mapping. A decision classifier based on the generated rules classified the multispectral image into glacier and non-glacier areas. The result of maximum likelihood classification (MLC) was used to compare with the result of the classification based on the rough set theory. Confusion matrix and visual interpretation were used to evaluate the overall accuracy of the results of the two methods. The accuracies of the rough set method and maximum likelihood classification were compared, yielding overall accuracies of 94.15% and 93.88%, respectively. It showed the area difference based on rough set was smaller by comparing the glacier areas of the rough set method and MLC with visual interpreter, respectively. The high accuracy for glacier mapping and the small area difference for glacier based on rough set theory demonstrated that this method was effective and promising for glacier mapping.  相似文献   

19.
In this paper a method based on in situ measured reflectance spectra was developed for accurately mapping salt fields in the Taibei Salt Farm in Lianyungang City of East China. After radiometric correction and geometric rectification, six multispectral ETM+ bands were fused with the panchromatic band via principal components analysis (PCA). The fused data were used to map salt fields in the farm. Salt fields were mapped 91.95% correctly in comparison with ground statistics from raw multispectral bands. This accuracy level rose to 96.4% with the use of the panchromatic band (15 m resolution). However, PCA-fused data produced the highest accuracy of 98.8%. At the zone level, coarse resolution data resulted in an accuracy of 98.93% for crystallization ponds, but only 82.68% for condensation ponds. Use of the panchromatic band alone improved the accuracy for condensation ponds but decreased the accuracy for crystallization and evaporation ponds. Results derived from PCA-fused data are highly consistent with the ground statistics at a discrepancy between 0.35% for evaporation zone and 1.80% for condensation zone. It is concluded that a fine spectral resolution is important to separation of ponds of a varying salt content while a fine spatial resolution is important to their mapped area. Ponds of a lower salt content can be mapped more accurately from fused data than from raw spectral bands.  相似文献   

20.
This research focuses on the application of HyMap airborne hyperspectral data and ASTER satellite multispectral data to mineral exploration and lithologic mapping in the Arctic regions of central East Greenland. The study area is the Kap Simpson complex in central East Greenland. The Kap Simpson complex is one of the largest exposed Palaeogene felsic complexes of East Greenland. It has been the target of several mineral exploration projects. The analysis of the HyMap data produced a detailed picture of the spatial distribution of the alteration minerals in the Kap Simpson complex, unavailable from field-based studies in the area. The analysis of the ASTER data produced mineral maps which due to the moderate spatial and spectral resolution of the ASTER imagery can be useful for reconnaissance level mineral exploration. Colour composites of the mean normalized ASTER thermal bands display lithological information and detected a large felsic igneous intrusion that has not been shown on the recently compiled geological maps of the area. The results of this research have considerable potential to evaluate the use of hyperspectral and multispectral remote sensing for geological purposes in the Arctic regions of central East Greenland.  相似文献   

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

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