GISdevelopment.net ---> AARS ---> ACRS 1990 ---> Geology/Disaster 2

An automatic technique and its effect by using Landsat TM data to extract hydrothermal alteration information in subtropical volcanic rock area

Zhao Tuanhong, Zhang Fuxianf, Chen Nanfeng
Remote Sensing Division, Dept. of Earth Sciences,
Zhejiang University Hangzhou, China


Abstract
A new automatic technique for using Landsat TM data to extract hydrothermal alteration zones n subtropical volcanic rock area covered by vegetation is developed here, and its applications is very successful in Bamao Au and Ag multimetal mineral area, Xinchang country, Eastern Zhejiang province. Not only are known alteration zones in the area displayed effectively, but also several unknown hydrothermal alteration zone, silicified quarts viens, ore-controlling structures and volcanic apparatus, which are very useful for prospecting, are revealed here.

Introduction
Landsat TM and 7(2.08-2.35mm), in hydroxyl-bearing clay minerals exhibit strong absorption (2.20mm)1,3,4, is very useful for geological exploration. With the developing of Landsat-4,5, a series of studies for mapping altered rocks have been done by many schlars4,5,6,7, and have got some valuable results. But the recent literatures show that there are still two problems in this field:
  1. Most people select arid or semiarid area as study are. The effective studies of technique and application are still lacked under range vegetation, especially in the subtropical area.

  2. Enhancing the alteration information, people generally delineate the altered rocks by visual interpretation instead of automatic technique. Such a case with a bad accuracy and low working efficiency, cannot meet the needs of geological exploration.

    Therefore, it is in an interested problem how to use Landsat TM data for altered rocks enhancement and automatic extraction in vegetated area.
Alteration Zone Automatic Extraction
For improving the application effects, a new effective technique, which can use Landsat TM data for mapping alteration zone automatically under range vegetation, is designed here; its flow chart is show in Fig 1.

Its working principle as following :
  1. Alteration Information Enhancement

    TM5/7 is considered as an effective band ratio for alteration information enhancement in arid or semiarid area. But in vegetated area, vegetation impedes the alteration information discrimination of TM5/7 ratio image, because it is both widely distributed in the surficial environment and can be spectrally similar to clays when sampled by TM8,6. Analysis of the study area show that both green vegetation and clay minerals contribute to the lighter pixels on the TM5/7 image, and most of vegetation has a larger TM5/7 ratio than clays. Whereas we must find a new method to remove the effects of vegetation, so that we can enhance the alteration information effectively in such a area covered by vegetation.

    It is well known that TM4/3 band ratio serves as an excellent vegetation index image. The analysis of the study area shows that patches of vegetation correspond to the lightest pixels and far lighter that others do on the TM4/3 ratio image.


    Fig. 1. Flow chart of alteration zone automatic extraction

    The principle component transformation is used for analyzing the TM5/7 and TM4/3 band ratio image here. Once computed, the first principal component (pc#1) contain the contribution of main information in both ratios-that is vegetation; the second principal component (pc#2) contains the contribution of altered minerals. Such a method is called as "directed principal component analysis (DPCA)" by S.J Fraser and A.A Green. They have used this method for anlaysing ATM data in tropical savannah woodland with 50to 70 percent total vegetation cover, and the clay and vegetation have been identified successfully, but few people have done such a study with TM data in subtropical area.

    This method has been developed for using Landsat TM data enhance the alteration information here, and its application is very successful in Bamao area.

  2. Alteration zone automatic Extraction

    On the PC#2 image, the lightest pixels correspond to altered rocks, and the lightness is vegetally increased with altered rocks. For extracting the alteration zones exactly, following steps are taken here:

    1. Density slicing. A bilateral selection method, combining both sample training of known alteration zones and histogram training of PC#2 is used to select a threshold that can describe the real digital number (DN) range of altered pixels.

      Supposed that DN range of altered pixels is (D_e, 255), where D+e is lower limit, D is an approximate lower limit, and e is a correction factor (e<=0). An approximate mean value (M) of known altered pixels is calculated by statistics of alteration zone. Considering that the M is often decreased due to other alteration free pixels inside the "Statistics window", and relatively speaking, weaker alteration zone has little significance for geological exploration, we let D=M. The histogram of PC#2 can display the distribution sate of DN over all images directly. According to the signature of histogram near D, a best "e" can be found to determine the threshold, so that we can slice off the alteration information. Letting altered pixels be 1, other be 0, a(0,1) binary image (ALTA) is acquired.

    2. Random noise removing. On the ALTA image, random noise, coming from TM3, 4,5,7 band data, is one disturber. For removing such a disturbance, we have made further analysis to determine the statistical signatures of "1" pixels on raw TM data, and a "discrimination window" of alteration information is selected. Then such "1" pixels, being outside the "window" and being considered as random noise, can e removed, an a new image (ALTB) is made here.

    3. Alteration zone delineating. By further analysis of ALTB image, we find that there are two types of information displayed on it, one is contributed by altered rocks with small scale and weak alteration, which is almost no significance for geological exploration, and another is contribute by the altered rocks with, large scale and strong alteration, which had important significance for geological exploration. We must find a method to delete the first one and remain the second one, but alteration zone is displayed as discrete cluster on the image, some other information, such as water bodies, residential points inside the alteration zone, have disturbed it. So it is difficult to filter such information using ordinary methods.

      A new filter algorithm - called as contexture filter, is presented here: a first, we use a "close" algorithms of "Dilation-erosion" to connect all these discrete "1" pixels one alteration zone: and then statistical signature of "1" pixels in spectrum and space is computed with scanning and tracing method" At last, sample training is made to determine the principles for filtering, so that we can extract the valuable alteration information.

      Using this algorithms, not only can the weakness of ordinary filter methods, with which discrete points inside the valuable alteration zone will be deleted, be overcome, but also the edge effects that will be induced by contexture or texture methods can be avoided. Its effects are very satisfactory.

    4. Displaying. Using "close' method to connect eh discrete clusters into consistent clusters; adding TM1 images a background image on it ; Letting different colors represent varieties of alteration grades respectively. Then an alteration zone displaying image, on which had topographic information, is made well.
Application Results
Bamao is situated in subtropical climate zone with 70-80 per cent total cover in Eastern Zhejiang province. Upper Jurassic and Lower Cretaceous volcanic rocks are the predominant rock units exposed at the surface. T shi range is chosen as the study area because of the presence of extensive hydrothermal alternation, excellent mineralized conditions. Nine mineral occurrences have been discovered here, and all their orebodies occur insilicified quarts veins, which are enclosed by clay zones.

The Ladnsat-5 TM data selected for analysis were acquired on February 11, 1987. For distinguishing the intensity of clay alteration, an alteration zone-displaying image on which two grades of alteration are displayed is made successfully, and our field checking has shown that the displaying image has very good effects.

Fig.2. Displaying image of alteration zone, Red color--Strong alteration, Green color--Weak alteration


Fig.3. Schematic interpretation map

  1. All nine known mineralized alteration zones are displayed exactly, except Changzhao (A7) (with near 100% vegetation cover).

  2. Distribution shape of alteration zones are expressed very well. for example, Qianguanling (A4) is a strong alteration district with a wide range; the linear alteration zone trending SN in Houan (A5) corresponds to a silicified quarts vein that trend SN; and huanggongkeng (A9) is contributed by an epidotization zone trending NW etc.

  3. There are alteration zoning in some known mineral occurrences, it is displayed as the strong alteration zone surrounding by weak alteration. Such as Qianguanling, Houan, Xiaolontinshan (A1), Caopenghsna (A2) etc.

  4. We have checked seven previously unknown alteration zones revealed on the image, and find it all be displayed exactly. Such a Wangjianshan (B5), eastern Houan (B40 are strong sericitization; Bainu (B6) is a chloritization zone; and Xianjiancun (B2) corresponds to a silicified quarts vein. Especially, there is a crossed alterationzone displayed in Western Lanyan (B3), the one trending NNE corresponds to a silicified quarts vein, and another one trending NWW is controlled by a fault by which there is a strong sericitization zone.

  5. The distribution of alteration zones corresponds to known ore-controlling structures, volcanic apparatus. For example, Xiaolongtingshan, Caopengshan, Xishan (A3), Qianguanling lie in the Bamao volcanicn basin; furthermore, by the fault passed 1. Xinchang country seat; 2. Xinchang river; 3. Boundary of training area; 4. Changzhao reservoir; 5. Known mineral occurrence and its serial number; 6. Checked unknown alteration zoen and its srail numerb; 7. Inferred volcanic vent; 8. Inferred linear structure; 9. Inferred silicified quarts vein; 10. Inferred mineral potential district and its serial number through Lanyan, Huanggongkeng, Liumenshan and Yanwo, a series of alteration zones are revealed. The alteration zone Wangjiangshan, Baibu, Huangputing is situated in a fault zone trending NE, and all the three zones lie in the edge of volcanic vents predicted here.
  6. Basing on the displaying image, we have delineated three strong alterationdistricts (C1, C2, C3) a regional tectonics just passed by it. Especially, the district C1 is displayed as very strong alteration area. Out field checking has shown it is mainly contributed by strong sericitization minerals, and silicified quarts vein also has been found there, perhaps it is a very valuable district for multimetal deposit exploration, and has been come into notice.

    In general, this study is successful. It has important significance for large scale mineral prospecting, and can supply can effective method for alteration information extraction, geological structure feature analysis and integrated mineral prospecting in subtropical volcanic rock area in well vegetation.
References
  1. Ailicia.W.Seigrists: Optimum Spectral Bands for Rock Discrimination. Photogrammetric Engineering and Remote Sensing, Vol. XL. VI, No.9 1980, 1027.

  2. Segal.D.B: Use of Landsat Multispectral Data for Definition of Limonitic Exposure in Heavily Vegetated AREas. Economic Geology, 78, 711, 1983.

  3. F.A. Kruse: Munsell Color Analysis of Landsat Color-Ratio Composite Images of Limonite Areas of Southwest New Mexico. Proceedings of International Symposium on remote Sensing of Environment, 3rd Thematic Conference: "Remote Sensing for Exploration Geology", 1984, 761-763

  4. D.R. Gladwell, P. Lawrence, M.Ancziger: The Application of Rapid, Semi- Quantitative Clay Mineral Determination at the Cortez Gold Mine, Nevada. Proceedings of international symposium on Remote Sensing of Environment. Fourth thematic Conference: "remote Sensing for Exploration Geology", April 1-4, 1985, 403.

  5. Nevin L. miller and Hristopher D. Elvidge: The iron Absorption Index : A comparison of Ratio-Based and Basetine-Based Techniques for the Mapping of Iron Oxides. Proceedings o Inernational Symposium on Remote Sensing of Environment, Fourth Thematic Conference: "Remote Sensing for Exploration Geology", April 1-4, 1985, 405-415.

  6. Elvidge.C.D, and Lyon R.J.P: Estimation of the Vegetation Contribution of The 1.65mm/2.22mm Ratio in Airborne Thematic Mapper imagery of The Virginia Range, Nevada. International Journal of Remote Sensing, No. 6, 1985,75.

  7. S.J. Fraser and A.A. Green: A Software Defoliant for Geological Analysis of Band Ratios. International Journal of Remote Sensing, Vol. 8, No. 3, 1987, 525-532.