用卫星高光谱数据提取江西德兴铜矿矿山废水的pH值污染指标_图

1、第22卷第1112期地质通报Vol22.N 0811。122003 年 12 月 GEOLOGICAL BULLETIN OF CHINADec.2003用卫星高光谱数据提取江西德兴铜矿矿山废水的pH值污染指标刘圣伟1,甘甫平z一 ,王润生S(1中国地质太学地球科学与资源学院，北京100083;2中国国土责掉航空物探遥 感中心，北京1000833北京走学地球与空间科学学院.北京100871摘要:在对德兴铜矿矿山废水的光谱特征深入分析研究的基础上.总结了不同类型水体(酸性水、碱性水以及河流 水的特征光谱.并利用地物谱特征开展矿山废水 pH值污染指标提取研究。针对亦体光谱反射率低、特征光潜不明显的

2、特点,采用矿区卫星H师面on高光谱数据，应用Isa算法和掩膜技术识别出水体分布并进一步与 MNF变换有效 结合,根据波段散点图进行不属pH值水体的有效分割。为矿山唆水 污染的诊断和监侧提供了新技木和理论支撵。关键词:高光谱;矿山废水;Hyperion;光谱特征；pH值;穗兴铜矿中圈分娄号：P627;X141;X.52文献标识码:A文章编号：1671 2552t200311 121013 08中国是世界上的矿业大国之一，矿业开发总规模居世界第3位。但是随着矿产 资源的开发，越来越 突出的矿山生态环境问题不仅威胁到人民的生命 安全，而且严重 制约了国民经济的发展。矿山废水 排放是矿山生态环境最严重

3、的污染源之一。矿 山废水主要包括采矿酸性废水和选矿碱性废水，并含大量重金属及有毒、有害元素 如果不进行治理，直接排放，不仅会对矿区及其周围环境造成污染和危害，而且会造 成矿山资源的浪费。近10多年来，在矿山废水 污染的治理方面，中国虽然已经取得 较大的进展，从单项治理发展到全面规划、综合治理，但是全国由矿 山废水污染引起 的环境恶化趋势还未得到有效的遏 制。中国目前每年因采矿产生的废水、废液的 排放总 量约占全国工业废水排放总量的100,6以上，处理率仅 为4.23%,大量未经处 理的废水排入江河湖海，污染严重，因此必须加大对矿山废水的治理力度。为了对污染的生态环境进行有效的治理和修复，需要及

4、时了解其受污染的程度及其变化趋势，以便给生态环境的健康状况 确诊”然后对症下药”。污染生态诊断 可以采用包括对水体、土壤、动 物和植被直接采样进行化学分析诊断、物理方法 诊断以及生物指示诊断等各种方法 “1尽管这些方法能够对不同地点污染的强度 及变化进行准确的测 定，但是工作程序通常繁杂而昂贵。不仅如此，而且 对水体、土 壤和动植物直接采样的技术无法提供适用于大范围污染地区制图的区域信息。然 而遥感技术却可以通过提供污染地区和周边区域的航空或卫星影像，仅使用直接采样技术进行区域调查所花代价的一小部分就能够监测区域乃至全球性大气、水 体、土壤、植被的污染情况，掌握污染源的位置、污染 物的性质及其

5、扩散的动态变 化，及时了解污染物对 生态系统的影响，从而采取积极的防护措施。在这种 情况下， 采用先进的高光谱技术对矿山废水污染进行直接识别尤为必要。高光谱地物识别 基于地物谱的指纹效应，利用高光谱数据的重建光谱与标准光谱或实测光谱的定量 对比分析，实现利用宏观手段（遥感技术进行微观地物（如水体、污染信息等的直接 识别,并且已在诸如地质、农业以及生态等领域 开展了深入研究。但矿山污染物因 光谱反射率低、光 谱混合（复合等现象明显而极难识别，需要做专门收稿日期：20030626;惨订日期：20030926基金项目：国家自然科学基金项目（40201034国土资源部 十五”重点科研项目（200220

6、6共同资助作者简介:刘圣伟（1972 一，男。博士 .从事高光谱遥感拄术及其应用方面的研究工作。E-mail:fianshanmoun HYPERLINK mailto: 万 方数据1014地质通报GEOLOGICAL BULLETIN OF CHD/A 一一一的研究。目前对矿区环境污染这个世界性难题 ，欧河，东 边为富家坞河。美一些发达国家已经相继利用高光谱技术展开矿山污染调查由于数据获取以及其他客观原因，中国在这方面的研究较少或刚刚起步，在未来几年内将陆续有 一系列搭载有高光谱仪的卫星上大，如ARIES 一 1卫星搭载的ARES 1以及 PRISM卫星都将搭载高光谱仪等。这些光谱仪的光谱分

7、辨率一般都 在几纳米到 20m之间,空间分辨率在几米到50m之间，个别还具有全色波段。其数据特点完全适 用于资源环境与污染调查以及动态监测等，并基本进入 商业运营。中国正在规划中 的资源环境小卫星（1+2星计划也主要是采用高光谱成像仪。这些将开创高光谱技术的航天时代标志着高光谱应用的高潮即将到来，并可能取代多光谱数据而成为下 一步应用的主流。德兴铜矿从20世纪70年代就开始设点进行废水监测工作，主要是采用堰流法 进行流量测定；80年代以来，环保部门对德兴铜矿各排污口的废水流量和 水质情况 进行了测定，其pH值、铜铁等金属离子含 量超标比较严重日。近2年来，国土资源 部先后组织 开展了以江西省在

8、内的全国5省为试点的矿山地质 环境调查工作.初步 建立了矿山环境调查治理数据 库系统，为矿山环境保护规划的制定和实施提供了 依 据。目前，以现时相的高分辨率多光谱遥感数据 为主，辅以以往时相的航空遥感图像， 采用现行的 遥感数据处理和信息提取技术，可以完成矿山尾矿 及固体废料排放现状 与动态变化、土地利用现状与 动态变化、次生地质灾害和地质环境等项调查工作,但对矿山环境污染状况进行定量分析还没有比较成熟的遥感支撑技术。因此，本文以江西省德兴铜矿矿山为实验区，在对德兴铜矿区水体光谱特征详细分析研究的基础上，利用卫星高光谱Hyperion数据从高光谱的图像识别和光谱识别 2方面进行了 水体pH值污

9、染指标的提取研究。1研究区概况德兴铜矿是一个世界级的大型斑岩铜矿，位于江西省东北部德兴县境内，已探明 铜储量800 x104t以上,年铜产量约占全国的1/4。矿区属中低山丘陵,地势东南高、 西北低，河流水系发育。矿区以北的乐安河是矿区主要的牛活用水来源，以南的德兴 河是德兴市以及下游的主要生活用水水源。以西为大坞德兴铜矿矿山废水主要由采矿场产生的酸性废水和选矿厂产生的碱性废水组成。当降水或地下涌 水流过废石场或露天采矿场的硫化矿石或废石时，硫化矿物的 氧化作用产生酸性废水：4FeS2+1502+14H20=8H2S04+4Fe(OH3碱性废水均来源于选矿厂，由于选矿作业是在 碱性条件下进行的，

10、在浮选过程中 需填加大量的石 灰和其他药剂。主要的碱性废水为大量尾矿和精矿脱水工序中产生的高碱性废水。德兴铜矿每年产生的酸性废水(pH=2.4-2.7约0.1x10 %3并含大量的Cu、Fe和Al等金属离子：每日产生碱性尾矿废水(pH=11.312.3和精矿废水 (pH=11.81224分另I约为22x1041113和2X104m3,所有废水由4个尾矿库以及 若 干酸性废水调节库进行处理和调节。德兴铜矿从 20世纪70年代初期就开始研究治 理矿山废水，废水治理经历了点源治理一集中治理一清洁生产加末端治理3个阶段,初步形成了目前的废水综合治理格局。但是，由于历史和现实的多方面原因，矿山环 保还滞

11、后于生产发展，存在诸如露天采矿场酸性废水未经处理、村民人为干扰引起酸性废水流失而直接排入河流、环保设施运行效率低、尾矿库排洪水和坝前渗水水质不稳定、pH值严重超标、悬浮物和金属离子含量也时有超标等许多问题，废水外排现象时有发生工酸、碱性废水破坏水体的自然缓冲作用，消灭和抑制细菌等微生物的生长，妨碍 水体的自净功能，腐蚀管道和机械等设施。酸碱污染不仅能改变水体的 pH值，而且大大增加水中的无机盐类和水的硬度，因此污 水的pH值被列为检验污水水质的重 要指标之一 01。2水体光谱特征采用ASD公司生产的FieldSpecPro FR便携式分光辐射光谱仪，选择矿区内 不同类型的水体分 别进行了野外光

12、谱测量，测试研究矿区碱性废水、酸 性废水以及 矿区流域内受矿山废水扰动的河流水系的光谱特征。酸性水的光谱特征德兴铜矿的杨桃坞废石场、祝家废石场、露天采 矿场、堆浸场均产生酸性废 水，输入酸性废水调节 库。酸性水pH值为2.5左右”水体浑浊，呈橙红色（图1 一 Ao其波谱曲线在640mn附近的橙红光谱区 内存在一个较强的且波形对称的反射 峰，峰值在深万方数据第22卷第11 1期刘圣伟等：州卫星高光潜数据提取汀西德兴铜矿J发水pH污染指标普毒啸1015罔1酸性水（A及其波谱曲线（B Fig1Acid water（Aand its spectra（B水处达28%,近岸处约为10%,随后反射率迅速降低

13、,最后在红外波段趋于平缓。近岸与远岸处水体的光谱曲线波形基本相似，近岸水体波谱640nm处的波峰反射率较高，推断是受到浅水处水体底质反射影响的缘故（图1 一 Bi碱性水的光谱特征尾矿库内澄清的尾矿碱性水 pH值为11312.3171。水质较为清澈，色发蓝（图2 A,其中含大量的选冶药剂。其波谱曲线与酸性水体明显不同（图2B,从380rim处反射率开始逐渐上升，在蓝绿光波段565nm附近达到峰值，然后迅速下降，形成左右不对称的反射峰。与酸性水体丰 u反该反射峰值随着测试点距岸距离的增加而升高，可能是悬浮物的A B浓度随之增大的反映。离岸340m处的波谱形态比 较相似，仪是反射峰币同。 离岸05m

14、处波谱曲线在580700i-lln波段出现较高的反射平台，推断也是受 到了浅 水处水体底质反射的影响.2.3流域河水光谱特征矿区东南高，西北低。流经矿区的水系主要有北 边的乐安河，西边的大坞河。大 坞河流向自南向北，并在西北汇人乐安河。本次主要对大坞河（图3 A和乐安河 （图3 B交汇处的河水进行光谱测试。由于受矿山酸性废水直接排人的影响，大坞河水体浑浊,蜀黄绿色。其波谱反射率值从 350nlrl开始迅速卜开，在580800nm之 间区域达到峰值并形成一高反射率平台，显巨岸由远到近反射峰值由低到高皑目2碱性水(A及其渡谱曲线(BFig.2 Alkaline water(Aand t氐 spec

15、tra(B第22卷第117：第22卷第111： B 万方数据哥善倒地质通报 GEOLOGICA 1.BULLETIN(IF CHINA 哥盎蜡图3河水(A,B搜其光谱曲线(c,DFig31Liver water(A,Band its spectra(C,D变化，同时波峰位置向长波方向偏移(图3 Co推断是由水体中悬浮物的浓度 和成分变化引起的，可能I兑明距岸较远的水域悬浮物浓度相应较高。曲线在450480nm之自J存在弱吸收带，在680nm、760l-lnl以及980nm附近存在明最的 强吸收，推断与水中悬浮 物所含的Fe、Cu等金属离子密切相关。乐安河水比大坞河清澈，水质比较洁净，足矿区生活

16、用水来源。其波谱反射率曲 线与大坞河水相比有 较大差异（图3 D,波潜反射率值明显降低，从350 nm处反射 率缓慢升高，在600850llril区间高反射率平台处存在一系列明显的吸收特征。可能 是受到水下 摹底沉积物的影响，距岸较近处反射率值相应较高。总之，矿DC内酸性废水以及受酸性废水污染的 夫坞河水的浑浊度最高，早橙红色,其波谱曲线在 橙红波段 具有强反射峰，波形对称分布，存在类似于Fe、Fe2+年lcu离予等金属离子跃迁 所产生的光谱 特征;而碱性尾矿库水则比较清澈，浑浊度低,呈蓝绿色,波谱曲线在蓝 绿光区具有特征的高反射峰，呈雎对称性分布。3高光谱数据处理及信息提取3.1Hyperi

17、on数据大气校正及去噪处理目前遥感影像已被广泛地应用于污染生态诊断及生产管理等领域，受到各国重 视。然而，由卫星获取的大量图像大部分被大气粒子吸收和散射来自地表的辐射引起的大气粒子效应所污染，大气校正的目的就是通过消除大气影响而从遥感影像中 恢复反映地物特性的表面反射率，从而使遥感影像分类的精确性得到较大的改善。 这、问题已经受到遥感研究人员的重视，并已经设计出许多解决方法。大气校正的 算法基本由2个主要步骤组成。苗先，利用地物表 面的特殊特性或大气组成的n接 测量以及利用理论模型对大气光学特征进行平估,根据大气的光学特性1016地须埴报可以通过辐射传输算法汁算 H1与人气校正有关的各万方数据

18、第22卷第1卜12期州圣伟等:用p星高光谱数据提取江西德兴铜矿蟆水 pH 弓染指标I【117种量值。然后，通过倒转过程得出地表反射率.使遥感影像得到校正，从而保证了卫星影像的有效性。采用江两彳惠兴铜矿地区2003年1月的Hyperion数据覆盖地面42kmX75km,、数据经系统校正和辐射定标处理之后，采HACORN v 4软件作大气校止，将辐射值数据转换成地而视反射率数据。ACORN是hnSpec公司基于 MODTRAN 4开发的大气校正和波潜重建软件。其原理是，用成像光谱数据940nm和i14011111附近水蒸气通道的吸收深度，逐个像元地反演水蒸气住的含量，代人大气传输模型MODTtCA

19、N 4计算大气中气体、分子和气溶胶的光学效应，对数据作大气校正.将其转换成地间视反射率数据。Hyperion数据除有少最云覆盖外，质量较女？,图像清晰。射图像进行去条带处理和运用空问域低通滤波对数据作平滑化处理以消除和压制条带和噪声后，图像质量有较大改善（图4,水体波谱曲线的特征谱带表现明显（图5。32水体信息及pH值污染指标提取321方法将本次野外文测的水体波谱与已有污染水体图4人气校一后的HyPenon影像Fig 4Hyperion image after atmospheric correctiorl波谱”进行综合分析，发现污染水体普遍县有以下 光谱特征:光谱反射率很低, 一般小于5%,

20、在浅水区水底质的影响p町达20%;由于悬浮物的存在 和不同酸碱 度水色的影响，往往在绿光和红光谱段具有特殊的光谱特征，但特征非常微弱。已有 研究证 叫，当目标地物的诊断特征在整个光谱范围内反射率都很低时，Tricorder和光谱特征匹配等光谱分类 方法无法取得很好的效果4-111。由于矿山废水的反射卒很低许且缺乏明确的光谱特征，为r对水体信息进行有效的提取，采用Ribeiro等提出的方法。该方法基于以下2个步骤:通过积分光 谱分析（ISA:Intergral Spectral Analysis对罔像进行增强处理，从而对低反射率区域进 行识别区分;使用最大噪声组分变换（MNF:Maxi mum

21、Noise Fraction Transformation和散点图分 析进行低反射率地物的分割。ISA算法是对各波段的反射率强度与它们各自的半极值宽度（FWHM:FullWidth of Half Maxi nltlll:的乘积进行求和。其计算公式为SI-2（RbFWHMb式L1,s,指光谱积分；R指亮度值。ISA算法的原理基于地物的反射率特性，即低反射率地物如水体等的光谱积分 值与其他地物卡|jI比要低得多。但是，由于该方法会引起对阴影或其他低反 射率地 物的错误指示，在这种情况下，要对相关因素 进行综合考虑，例如地形坡度可以用来对 阴影区域进行指示一。图5影像数据大气校正前后水体渡谱曲线对

22、比（蓝色校正前波谱，红色校正后波谱Fig 5Spectral contrast of waterbefore and after atmospheric correction(blue.before Correction;red,after第22堂第1种量以然后: 影像得到校正 采用江。第22卷第1H：种量值然后，通 星：的旦至I廿；H数据，覆盖地 辐射定标处卫 校正,将辐3correction 万方数据琏矮通报GEOLOGICAL BULLETIN OF CIIIB4 对低反射率区域的识别，目的 是进行掩膜处理并施用最大噪声组分变换（MNF。最大噪声绢分变换（MNF是利用图像的噪声组分矩阵（

23、01的特征向量对图彳t进行变换，使按特征值由大到小排序的变换分量所包含的噪声成分逐渐减小，而蚓像质 量顺玖提高。上为图像的总办方差矩阵.M为图像噪 声的价方篮矩阵。MNF相当于所有波段噪卢方差都 J【I等时的主成分分析，因此可分为2步实现，第一步 先将图像变换到个新的坐标系统，使变换后图 像噪声的协方差矩阵为单位阵；第二步再对变换 后的罔像施行丰成分变换。此改进的算法称为噪声凋节主成分变换（NAPC ”最大噪声组分变 换（MNF肯2个重要性质.一是对图像的任何波段 作比伤J扩展，变换结束不变；电变换使幽像矢量、 信息分量和加性噪声分量与.相垂直.乘性噪声可通 过对数变换转换为加性噪声。变换后可

24、钊对性地照 各分量图像进行去噪，或禽弃噪声占优势的分量， 最大噪声组分变换（MNF是由主H分分析（PCA 衍生出来的，它们之间的主耍【（别在于MNll变换考 虑到了主组分变换和噪音的数据犒篇、就幽像质量 方面而言，最人噪声组分变换（MNF被认为是一种 高光浒数据处理的响效斤法”j 200笄322结果和结论进彳T ISA增强和掩膜处弹后从网像中直接分割 出水体（图6,册掩膜图像以 MNF变换进行增强和 噪声消除处理 然后利用MNF变换不同波段的散点 图进行分 割。根据MNF波段1与MNFi0:段6的散点图（图7,町以将矿区水体夫致分为3类 （图9右。利用陵散点 图口阻区分出：酸性水（红色，位于祝

25、家废石场及其 附近的酸 性废水调节库和受酸性废水影响的大坞 河;碱性水（绿色，位于2号及4号尾矿库；以 及酸碱 度基本偏中性的水域（蓝色，位于尾矿库边缘澄清 程度较高的区域和矿区生 活水源地（图9右图左上角 的蓝色图斑。其光谱特征（图8与野外实测的不J司酸碱 度水的光谱特征类似（囤上幽3,反映出不同水 质在600nm附近的反射存在明显的差 异。根据德兴铜矿矿山酸碱Tt废水水质分析结果，米源于祝家废石场、露天采矿场 和堆浸场的酸性废水的pH值为24-27,酸性废水输人酸性废水调节库，部分直接排 人大坞河。来源于大山选矿场和泗洲选矿场的尾矿溢流水和精矿溢流水的 PH值分 别为11P12.3和11.

26、81224_排人2号及4号尾矿库。通过 与该结果的对比验证，发 现分类结果与德兴铜矿酸 碱性废水实际分布情况基本吻合，证明了此方法的 有效性 和适用性。由于水体光谱反射率低，特征光谱不明显，同时ffl6水体信息（红色的光谱识别 （左及其影像光谱（右Fig 6discrhalination of waSpectml00 r information(red,leftand its image1018对低 并施用最 换（MNF） 征向量对 的变换分 F=L 七在 a H fspectm（nghtloirt对低反射率 并施用最大噪声 换（MNF）是利用 征向量对图像进第22卷第1112期刘圣伟等：用卫

27、星高光谱数据提取江西德兴铜矿腹水pH污染指标 判7MNF波段1与MNF波段6散点罔Fig 7Scatterplot bctVI ,een MNF B1and MNF B6Fig8Image spectra ot water bodyi中不同颜色所表示的影像九皓与闭9地物对应图9Hyperion影像（左与水体酸碱性信息提取（右Fig9Hyperion image（1eftand relative pH for various kinds of water（right,red:relative low pH;blue:relativemiddle pH;gxeen:relative high pH

28、一些小水体由于卒问分辨率以及水体深浅的影响，其光谱特征困混合而难以提 取出存在信息。采用积分光 谱分析，即IsA算法增强和抢膜处理技术，可以对水体 进 行识别分割。水的酸碱性从光谱特征上雌以区分 ，但通过MNF变换后利用散点图 可以进行相对划分。册于因铜矿进行的酸性废水和尾矿碱性废水的酸碱巾和治理.蹦及尾矿砂的自然分割造成尾矿庠内水域第22卷第1172期刘圣伟等：用卫星7第22卷 第1132期 刘圣伟等：用卫专酸碱性的差异，需要进一步的野外验证，以修改技术1020地质通报GEOLOGICAL 在矿山废水信息提取和矿方法，深化高光谱遥感技术山废水污染诊断和监测工作中的应用。参考文献1孙铁琦，周启

29、星，李培军污染生态学Mj北京:科学H版社，200112j Mars J C,Crowiey j K Mapping rmne w*tes and analyzing3rea affected by selenium rich waterrunoff dieast ldahousing AVIRIS imagery and 山目出 elevation damUP.emote Seining ofEnvi nt,2003,84:4224363Swayze G A.Smith K S,Clark R N.et al Using i 蚂 ging3pecII uw。Py to mapacidic mi

30、ne wa=steJEnvironmental Science and Technology,2000,34:47544Minitorfog and assessing the environment impact of mining inEu pe using advance earth obsemnontechniquesfMINEOM/OL http:/www brgm ff/rahlco5德兴铜矿忐编委会德兴铜科志M19936朱训，黄祟轲,芮宗瑶，等.德兴斑岩铜矿M北京:地质出版 社.19817昊飞德兴铜矿矿山废水治理现状及其前景U铜业工程.2000,1:27298赵冬覃，刘玉机中国污

31、染水体光谱特征M北京:海洋出版 社,20019Clark R N,Galla 曲 el-A J,Swayze G A Material absorptionBULLETIN OF CHINA 2003 年 band depth “pping of imaging spectrometer data using a complete bal d shape leassquares fit w m hbrary reference spectraAIn:Summaries of the 2nd Annual JPL Airborne Geoscience Worlcshop CI JPL Pnhh

32、cafion,1990.54:176 186101Clark R N,Swayze G A Mapping minerals,amorphou s nutefials,envlro- entalmaterials,vegemdon,water.ice and snow,and other materials:The USGS 6corder algorithmA In:Suedes of the Fifth JPL Airborne Earth Science WothshopCJPL Pubhcadon.1995.1:3ch4011Ribeiro M N C.Carvalho J o A,G

33、uimar=ies E M,et al Integral spectral analysis(ISAapplied to AVIRISdata for manganesemineralized larefitcs inSio Jogo山 Alianga/GO Brazil f c Ninth JPL Aithome Visible/lnfra md Imaging Spectrometer(AVIRISWorkshop.JPL Publication,2000I“1812】Green A A,B M,SvAtzer P. al A transformafionfor orderingmuMsp

34、ectraldata in cof image qualityLm implicationsfor noise removal H IEEE Tram Geosci.Re mONb Sens,1988,26(1:657413Lee J B,Woodyart S,Bennan M Enhancement of higll spectral resolution ieH lOt sensing data byam e Adjust pfincipal components 密 afformJIEEE Trans.GecBei Re- mote Sens.1990.28(3:295-304Using

35、 hyperion data to extract pH information of mine waste water in the Dexing copper mine,Jiangxi Province,ChinaLlU Shengweil.GAN Fupin92 ,WANG Kunshen91 2(1School ofEcwth Sciences and Resources,China Univemity ofGeosciences.Beqing 100083,China;2.China Aero Geop女口 ieal Survey and Remote Sensing Center

36、for Land and Resources,& “/ng0083t China;3.School of Earth and Spoce Sciences,Peking Unwenizy,Belling 100871,ChinaAbstract:On the basis of an in depth analysis of山 e spectral features of mine waste water in the Dexing copper mine.the characteristic spectra of different typesof waters(acid water,alka

37、line water and river water were summed up and the pH information extraction frommine waste water was studied using the hyperspectral imaging data.Due to the lowreflectance and lack of specific absorption features of waters, the ISA algorithm and mask were used to dlscriminate the waters Then the seg

38、tnentation of waters with different pH values was accomplished by using MNF and scatter-plot analysis.This smdy can supply and help thewhole process of diagnosing and monitoring rfline waste water pollution.Key words:hyperspecrtal imaging;minewaste water;hyperion;spectral features;pH;Dexing copper m

39、ine用卫星高光谱数据提取江西德兴铜矿矿山废水的pH值污染指标作者:刘圣伟，甘甫平,王润生作者单位：刘圣伟（中国地质大学地球科学与资源学院，北京,100083 ,甘甫平（中国国土资源 航空物探 遥感中心，北京,100083;北京大学地球与空间科学学院，北京,100871，王润 生（中国地质大 学地球科学与资源学院，北京,100083;中国国土资源航空物探遥感中 心，北京,100083刊名:地质通报英文刊名：GEOLOGICAL BULLETIN OF CHINA 年，卷（期:2003,22（11 被引用次 数：7次参考文献（13条1.吴飞德兴铜矿矿山废水治理现状及其前景期刊论文卜铜业工程200

40、0(012.朱训;黄崇轲;芮宗瑶德兴斑岩铜矿19813.德兴铜矿志编辑委员会 德兴铜矿志1993Minitoring and assessing the environment impact of mining in Europe using advance earthobservation techniques(MINEOSwayze G A;Smith K S;Clark R N Using imaging spectroscopy to map acidic mine waste外文期刊20006. Mars J C;Crowley J K Mapping mine wastes and

41、 analyzing areas affected by selenium-rich waterrunoff southeast Idaho using AVIRIS imagery and digital elevation data阶文期干U 2003(3Clark R N;Gallagher A J;Swayze G A Material absorption band depth mapping of imaging spectrometerdata using a complete band shape least-squares fit with library reference

42、 spectra 1990驱冬至；刘玉机 中国污染水体光谱特征 2001LEE J B;Woodyatt S;Berman M Enhancement of high spectral resolution remotesensing data by anoise-Adjust principal components transform卜文期干U 1990(03Green A A;Berman M;Switzer P A transformation for ordering multispectral data in terms of imagequality with implications for noise removal外文期刊1988(01Ribe