首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 62 毫秒
1.
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.  相似文献   

2.
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.  相似文献   

3.
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.  相似文献   

4.
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.  相似文献   

5.
Advanced space-borne thermal emission and reflection radiometer (ASTER) data were evaluated for the hydrothermal alteration mapping of the Nimu porphyry copper deposit, southern Tibet. According to the metallogenic model for porphyry copper, we chose hydrothermal alteration mineral association and then established a remote sensing model. Relative absorption-band depth (RBD) and simple band combination methods in visible near infrared (VNIR) and shortwave infrared (SWIR) bands were used to retrieve information about lithological distributions. Principal component analyses (PCA) were applied to extract the prospecting information based on the spectral information of argillization mineral association, propylitization mineral association, and iron stained (limonite) effects. Results of the study match up well with known copper occurrences and a circular structure in study area, while circular structures usually have close relationship with mineralization. Combined with field validation, ASTER data are proved able to characterise the alteration zone of these porphyry deposits. With VNIR and SWIR bands, argillic, propylitic and ferritization alterations in a single deposit can be effectively discriminated.  相似文献   

6.
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.  相似文献   

7.
The Sanandaj-Sirjan Zone (SSZ) is considered as an important region for gold exploration in the western sector of Iran. Its mountainous topography and unpaved routes make its study challenging for researchers and raise the costs for mining companies strating new exploration plans. Gold mineralization mainly occurs as irregular to lenticular sulfide-bearing quartz veins along shear zones in deformed mafic to intermediate metavolcanic and metasedimentary rocks. In this investigation, ASTER data are used for mapping hydrothermal alteration minerals and to better discriminate geological structural features associated with orogenic gold occurrences in the area. Image transformation techniques such as specialized band ratioing and Principal Component Analysis are used to delineate lithological units and alteration minerals. Supervised classification techniques, namely Spectral Angle Mapper (SAM) and Spectral Information Divergence (SID) are applied to detect subtle differences between indicator alteration minerals associated with ground-truth gold locations in the area. The directional filtering technique is applied to help in tracing along the strike the different linear structures. Results demonstrate that the integration of image transformation techniques and supervised classification of ASTER data with fieldwork and geochemical exploration studies has a great efficiency in targeting new prospects of gold mineralization in the SSZ. The approach used in this research provides a fast, cost-efficient means to start a comprehensive geological and geochemical exploration programs in the study area and elsewhere in similar regions.  相似文献   

8.
In the process of exploring pre-earthquake thermal anomalies using satellite infrared data, Blackett et al. (2011) found that the previously reported anomalies before the 2001 Mw 7.7 Gujarat earthquake, in India, were related to positive biases caused by data gaps due to cloud cover and mosaicing of neighboring orbits of MODIS satellite data. They supposed that such effects could also be responsible for other cases. We noted a strip-shaped TIR anomaly on March 17th, 2010, 28?days before the Ms. 7.1 Yushu earthquake (Qin et al., 2011). Here we again investigate multi-year infrared satellite data in different bands to discriminate whether the anomaly is associated with the earthquake, or is only bias caused by the data gaps. From the water vapor images, we find lots of clouds that have TIR anomalies. However, on the cloudiness background, there is an obvious strip-shaped gap matching the tectonic faults almost perfectly. In particular, the animation loops of hourly water vapor images show that the cloud kept moving from west to east, while they never covered the strip-shaped gap. We consider that the cloud with this special spatial pattern should have implied the abnormal signals associated with the seismogenic process. Based on current physical models, the satellite IR anomalies both on TIR and water vapor bands can qualitatively be explained using synthetic mechanisms.  相似文献   

9.
Multispectral data covering an area near Marysvale, Utah, collected with the airborne National Aeronautics and Space Administration (NASA) 24-channel Bendix multispectral scanner, were analyzed to detect areas of hydrothermally altered, potentially mineralized rocks. Spectral bands were selected for analysis that approximate those of the Landsat 4 Thematic Mapper and which are diagnostic of the presence of hydrothermally derived products. Hydrothermally altered rocks, particularly volcanic rocks affected by solutions rich in sulfuric acid, are commonly characterized by concentrations of argillic minerals such as alunite and kaolinite. These minerals are important for identifying hydrothermally altered rocks in multispectral images because they have intense absorption bands centered near a wavelength of 2.2 μm. Unaltered volcanic rocks commonly do not contain these minerals and hence do not have the absorption bands.A color-composite image was constructed using the following spectral band ratios: 1.6μm/2.2μm, 1.6μm/0.48μm, and 0.67μm/1.0μm. The particular bands were chosen to emphasize the spectral contrasts that exist for argillic versus non-argillic rocks, limonitic versus nonlimonitic rocks, and rocks versus vegetation, respectively. The color-ratio composite successfully distinguished most types of altered rocks from unaltered rocks. Some previously unrecognized areas of hydrothermal alteration were mapped. The altered rocks included those having high alunite and/or kaolinite content, siliceous rocks containing some kaolinite, and ash-fall tuffs containing zeolitic minerals. The color-ratio-composite image allowed further division of these rocks into limonitic and nonlimonitic phases. The image did not allow separation of highly siliceous or hematitically altered rocks containing no clays or alunite from unaltered rocks.A color-coded density slice image of the 1.6μm/2.2μm band ratio allowed further discrimination among the altered units. Areas containing zeolites and some ash-fall tuffs containing montmorillonite were readily recognized on the color-coded density slice as having less intense 2.2-μm absorption than areas of highly altered rocks. The areas of most intense absorption, as depicted in the color-coded density slice, are dominated by highly altered rocks containing large amounts of alunite and kaolinite. These areas form an annulus, approximately 10 km in diameter, which surrounds a quartz monzonite intrusive body of Miocene age. The patterns of most intense alteration are interpreted as the remnants of paleohydrothermal convective cells set into motion during the emplacement of the central intrusive body.  相似文献   

10.
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.  相似文献   

11.
In this work historical investigations and modern results of classification of the Krasnoyarsk Reservoir are presented. The paper presents results of studying the dynamics of phytopigments and other optically active components, using multispectral satellite data. Several approaches to interpreting satellite data for optically complex inland water bodies are offered. Based on results of historical investigations it is shown that the spatial distribution of phytoplankton in the reservoir stems back to the time of its formation. Color index in the red spectral region (CIR) is introduced. A relationship between the color index and chlorophyll concentration is investigated. The CIR, derived from the AVHRR data, has been found to be related to chlorophyll concentration. Based on MODIS data, the waters of the Krasnoyarsk Reservoir have been classified in accordance with their optical spectral variability, using the technique of unsupervised IsoData classification. An empirical relationship between multispectral MODIS data and the ground-truth measurements of chlorophyll concentration has been found.  相似文献   

12.
Current satellite imaging systems offer a trade-off between high spatial and high spectral resolution providing panchromatic images at a higher spatial resolution and multispectral images at a lower spatial resolution but rich in spectral information while a wide range of applications need the highest level of this information, simultaneously. Image fusion techniques as means of enhancing the information content of initial panchromatic and multispectral images produce new images, titled pan-sharpened, which inherent the advantages of the initial images. Considering the impact of fusion accuracy on the quality of corresponding applications, it is necessary to evaluate the quality of these processed images. During the last decade, a lot of quality evaluation metrics have been proposed which are mostly inspired by traditional image quality metrics. These methods are mostly based on applying quality metrics at the pixel level and evaluating final quality value based on averaging of obtained metric values through the whole image. However, obtained results clearly show that the behaviour of image fusion quality is inconsistent amongst different image objects. In this article, by applying image fusion quality metrics (IFQMs) to image objects, an object-level strategy for quality assessment of the image fusion process is proposed. The proposed strategy is applied to different satellite imagery covering residential and agricultural areas. Experimental results show higher capabilities of object-level quality assessment strategy in the quality assessment of the fusion process. Evaluating fusion quality at the object level provides the potential of fusion quality assessment for each individual image object in compliance with different parameters such as the type of objects and the effective size of objects in data set.  相似文献   

13.
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.  相似文献   

14.
星载多光谱遥感器太阳定标技术的进展   总被引:1,自引:0,他引:1  
星载多光谱遥感器的太阳定标器一般选择太阳作为基准光源 ,通过它将太阳辐射引入星载遥感器并调节到星载遥感器的动态范围内 ,对星载遥感器进行绝对辐射定标 ,也可对星载遥感器性能变化进行监测和校正。文章介绍一些最具代表性的星载多光谱遥感器的太阳定标器 ,并进行了分析 ,以反映太阳定标技术的现状与发展。目前采用太阳漫射器的星上定标方法可以实现全视场、全孔径、端点到端点的定标。这一方法的严重缺点是漫射器反射比随时间变化。为了解决这一问题 ,设计一种比辐射计或反射比标定装置来监测漫射器的辐亮度、反射比以及太阳常数。因此这种定标方法是很有前途的 ,在现代一些先进的星载多光谱遥感器上获得应用。通过对此方法的分析 ,提出了在太阳漫射器研制中的一些关键技术。  相似文献   

15.
Flood forecasting has long been a major topic of hydrologic research. Recent events and studies indicate that the success of flood forecasting in Taiwan depends heavily on the accuracy of real-time rainfall forecasting. In this study, we demonstrate a multi-spectral spatial convolution approach for real-time rainfall forecasting using geostationary weather satellite images. The approach incorporates cloud-top temperatures of three infrared channels in a spatial convolution context. It not only characterizes the input–output relationship between cloud-top temperature and rainfall at the ground level, but also is more consistent with physical and remote sensing principles than single-pixel matches. Point rainfall measurements at raingauge sites are up-scaled to pixel-average-rainfall by block kriging, then related to multi-spectral cloud-top temperatures derived from Geostationary Meteorological Satellite images by spatial convolution. The kernel function of the multispectral spatial convolution equation is solved by the least squares method. Through a cross-validation procedure, we demonstrate that the proposed approach is capable of achieving high accuracy for 1- to 3-h-lead pixel-average-rainfall forecasting.  相似文献   

16.
温场可编程控制的小型空间晶体炉设计与测试   总被引:1,自引:0,他引:1  
详细介绍了一台温场可编程控制的小型空间晶体炉。其工作原理是,在计算机控制下,炉膛温场发生有序地变化,完成晶体生长过程。由于没有机械运动部件,从而彻底消除了机械振动对微重力环境的干扰。此炉是为返回式卫星搭载试验而设计,配备以Intel80c31为核心的计算机温度控制及数据采集系统。该系统采用数字PID方式独立调节加热器中7个结构相同的加热单元的温度,可以编程实现多种温场变化模式。  相似文献   

17.
Since the 2001 eruption, the volcanic fumes from the Oyama of Miyakejima Island, located approximately 200 km southwest from the mainland of Japan, have affected the vegetation health and biomass of the Tokyo metropolitan area. In this study, we evaluated the potential for measuring/mapping relative forest damage/recovery in forests in the Bosoh Peninsula, Chiba Prefecture, Japan, using LANDSAT Enhanced Thematic Mapper data and ASTER Level-1B data. The simple ratio which is derived from the near-infrared and red reflective response was found to correlate well with ground-based measurements of forest damage caused by continuous SO2 contained in the volcanic fumes. Accumulated data by the synchronized field campaign observation of the Japanese cedar field spectra and by satellite analysis indicates the following interesting features: (1) regional differences of SO2 levels in the forestal area can be estimated by the simple spectrum ratio of the near-infrared and red bands of multispectral satellite data; (2) time-series images using the simple ratios indicated the location of still damaged and recovered forests between the 2001 eruption and 2003; (3) if the vegetations under study are homogeneous single-species with a similar leaf area index and a similar age, the simple ratio was useful to evaluate vegetation damage or recovery among individual forests throughout many of the forests in the Bosoh Peninsula; (4) the damage level of the Tokyo Bay side was found to be consistently higher than the Pacific Ocean side of the peninsula. This obvious difference is likely caused, not only by SO2 from Miyakejima Island, but also by long term and chronic discharged gas from heavy industrial complexes.  相似文献   

18.
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.  相似文献   

19.
The Arabian-Nubian Shield (ANS) underwent a regional carbonation event with CO2 fluxes derived from the mantle and circulated along post-accretionary shear zones, thus affecting the redistribution of gold mineralization. Formerly, the analyses of remotely-sensed data have suggested regions of gold mineralization linked to the post-accretionary structured framework, but the results were insufficient for detailed prospecting in the ANS. In this research, aeromagnetic data and Landsat-8 imagery were integrated for delineating new high potential zones of gold mineralization in the Allaqi-Heiani Suture (AHS) zone, the ANS, South Eastern Desert (SED) of Egypt. Aeromagnetic data were enhanced using the Center for exploration targeting (CET) grid enhancement technique to detect the main structures that control hydrothermal alterations in the study area. Principal Component Analysis (PCA) technique and ratios of spectral bands were applied to Landsat-8 data for mapping zones of hydrothermal alteration associated with gold mineralization. The traced structural elements from RTP, CET and Landsat-8 maps were statistically plotted and interpreted using rose diagrams. These diagrams indicated that the NW and NE trends are the most dominant in the RTP map and that the intersection zones obtained from the CET analysis are aligned along the NW-oriented trend. While the NNW, WNW and NE trends were the most dominant ones in Landsat maps. The integration of the results indicates several new high potential zones of gold and Cu/Ni mineralization types, which are mostly concentrated at the western part of the AHS zone. The hydrothermal alteration zones and associated gold mineralization are strongly linked to NW, NE, E-W and WNW trending lineaments. The approach used in this study can be applied to other parts of the ANS and other regions with similar geological conditions around the world.  相似文献   

20.
遥感图像融合的目的是融合高光谱分辨率、低空间分辨率的多光谱(MS)图像和高空间分辨率、低光谱分辨率的全色(PAN)图像,得到高光谱分辨率与高空间分辨率的融合图像。遥感图像的注入模型中如何确定注入细节及注入系数是该技术研究的关键。针对注入细节优化,先通过模拟MS传感器的特性来定义一种多尺度高斯滤波器,再用该滤波器卷积PAN图像以提取细节,得到与MS图像高度相关的细节。针对注入系数优化,综合考虑光谱信息与细节信息提出一种自适应的注入量系数。为更好地保留边缘信息,提出一种新的边缘保持权重矩阵,实现光谱信息与空间的双保真。将优化后的注入系数与注入细节相乘注入到上采样后的MS图像中,得到融合结果。对所提方法进行性能分析,并在各卫星数据集上进行大量测试,与一些先进的遥感图像融合方法进行对比,实验结果表明,所提方法在主观与综合客观指标上都能达到最优。   相似文献   

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

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