淮南积水沉陷区水深遥感定量分析.pdf

第3 6 卷第4 期中国矿业大学学报 V 0 1 ．3 6N o ．4 2 0 0 7 年7 月J o u r n a lo fC h i n aU n i v e r s i t yo fM i n i n g8 LT e c h n o l o g yJ u l y2 0 0 7 文章编号1 0 0 0 1 9 6 4 2 0 0 7 0 4 0 5 3 7 0 5 淮南积水沉陷区水深遥感定量分析 武彦斌，彭苏萍，黄明，吴观茂 中国矿业大学煤炭资源与安全开采国家重点实验室，北京1 0 0 0 8 3 摘要为获取沉陷区积水深度，尝试采用L A N D S A T 5 T M 影像建立水深遥感定量关系．首先 对T M 影像进行几何校正、大气校正、沉陷区提取等步骤，然后输出像元反射率值，并与水深实 测控制点坐标匹配，使水深值与反射率值对应．计算T M 影像各波段及波段组合反射率与水深 值之间的相关性，其中T M 4 波段的相关性最大 0 ．5 6 1 ，以T M 4 波段为例建立水深与像元反射 率间的回归模型．模型以水深值2m 为阈值，包括水深实测值小于2m 和在2 ～6m 的2 个回归 方程．将模型应用到淮南潘一矿的沉陷区，计算反演水深值与实测值误差．结果表明水深值小于 2m 的区域平均绝对误差0 ．1 5 91m ，平均相对误差1 3 , 3 6 7 ％；水深值在2 ～6m 范围内的区域 平均绝对误差0 ．5 7 18m ，平均相对误差1 5 ．0 2 ％． 关键词沉陷区；水深；遥感；回归模型 中图分类号T P7 9 ；P2 3 4 ．3文献标识码A Q u a n t i t a t i v eA n a l y s i so fW a t e rD e p t hi n Hu a i n a nP o n d i n gS u b s i d e n c eB a s e do nR e m o t eS e n s l ‘n g W UY a n b i n ，P E N GS u p i n g ，H U A N GM i n g ，W UG u a n - m a o S t a t eK e yL a b o r a t o r yf o rC o a lR e s o u r c e sa n dS a f eM i n i n g ，C h i n aU n i v e r s i t yo fM i n i n g T e c h n o l o g y ， B e i j i n g1 0 0 0 8 3 ，C h i n a A b s t r a c t T om e a s u r et h ew a t e rd e p t hi np o n d i n gs u b s i d e n c e ，am o d e lb a s e do nT Mi m a g ew a s p r o p o s e d ．A f t e rg e o m e t r i cc a l i b r a t i o n ，a t m o s p h e r i cc o r r e c t i o na n dp o n d i n gs u b s i d e n c ee x t r a c t i o n ，t h ep i x e l sr e f l e c t i v i t yw a se x p o r t e d ．I no r d e rt of o u n dt h er e l a t i o nb e t w e e na c t u a lw a t e r d e p t ha n dp i x e l sr e f l e c t i v i t y ，t h ep i x e l sr e f l e c t i v i t yw e r em a t c h e dt oc o n t r o lp o i n t s ．T h ec o r r e l a t i o nc o e f f i c i e n t sb e t w e e nr e f l e c t a n c eo fb a n d s c o m b i n e db a n d s a n dw a t e rd e p t hw e r ec a l c u l a t e d ．T h ec o r r e l a t i o nc o e f f i c i e n tb e t w e e nr e f l e c t a n c eo fT M 4a n dw a t e rd e p t hi St h em a x i m a l 0 ．5 6 1 ．T h er e g r e s s i o nm o d e lw a se s t a b l i s h e db e t w e e nt h ea c t u a lm e a s u r e dw a t e rd e p t ha n d r e f l e c t a n c eo fT M 4 ．T h et h r e s h o l do ft h em o d e li s2m ，i n c l u d i n gt w or e g r e s s i o ne q u a t i o n s w h i c hc o r r e s p o n dt oa c t u a lm e a s u r e dw a t e rd e p t hl e s st h a n2ma n dw a t e rd e p t hf r o m2m t o6 m ．T h em o d e lw a sa p p l i e dt om e a s u r et h ed i f f e r e n c eb e t w e e nt h ea c t u a lw a t e rd e p t hr e t r i e v e d w a t e rd e p t hi np o n d i n gs u b s i d e n c eo fH u a i n a n ．T h er e s u l t ss h o wt h a tt h em e a na b s o l u t ee r r o r i s0 ．1 5 91ma n dt h em e a nr e l a t i v ee r r o ri S1 3 ．3 6 7 ％，w h e nt h ea c t u a lm e a s u r e dw a t e rd e p t hi S l e s st h a n2m ．T h em e a na b s o l u t ee r r o ri S0 ．5 7 18ma n dt h em e a nr e l a t i v ee r r o ri S1 5 ．0 2 ％， 收稿日期2 0 0 6 0 9 0 6 基金项目国家重点基础研究发展规划 9 7 3 项目 2 0 0 2 C B 2 1 1 7 0 7 ，2 0 0 6 C B 2 0 2 2 0 9 ；国家自然科学基金重大项目 5 0 4 9 0 2 7 1 ；“十五”国 家科技攻关重点项目 2 0 0 1 B A 8 0 3 8 0 4 0 3 ；国家创新研究群体基金 5 0 2 2 1 4 0 2 ；教育部“长江学者和创新团队发展计划” I R T 0 4 0 8 ；国家杰出青年基金项目 5 0 0 2 5 4 1 3 ；教育部科学技术研究重大项目 3 0 6 0 0 2 作者简介武彦斌 1 9 8 0 一 ，男，河北省藁城市人，博士研究生，从事遥感地质应用方面的研究． E m a i l w u y a n b i n 0 8 0 1 6 3 ．c o r n T e l 0 1 0 5 1 7 3 3 0 2 0 万方数据 5 3 8中国矿业大学学报第3 6 卷 w h e nt h ea c t u a lm e a s u r e dw a t e rd e p t hi sf r o m2mt o6m ． K e yw o r d s p o n d i n gs u b s i d e n c e ；w a t e rd e p t h ；r e m o t es e n s i n g ；r e g r e s s i o nm o d e l 淮南矿区地处淮河两岸，地下水位埋藏较浅， 矿区沉陷区大面积积水，利用遥感手段量测水深可 以发挥遥感“快速、大范围、准同步”的优势，为沉陷 区治理提供可靠的数据流．水深遥感国内外都有过 研究，主要集中在沿海及内陆淡水水域．例如 B i e r w i r t h [ 1 ] ，M g e n g e [ 2 3 曾分别利用T M 影像有效 地对近海岸水域进行了多时相水深测量；党福星利 用T M 影像数据对南海永暑礁景区水深反演模型 进行了详细的研究[ 3 1 ；黄家柱等利用T M 影像数 据建立了长江南通河段水深遥感模型[ 4 1 ；陈呜等在 对长江口水深的分段遥感模型研究基础上，还进一 步开展了利用该模型对长江口水域水沙特性的探 讨凹3 ；于瑞宏等利用多时相的T M 影像对干旱区 草型湖泊乌梁素海的水深遥感定量模型也进行了 有益的研究邸] ．但目前有关针对煤矿沉陷区积水水 体的波谱测试和波谱信息的研究还未见详细报道． 业已证明，从卫星遥感图像上得到的影像光谱数据 经处理后与水深存在很好的相关性，因此可以直接 利用其相关性把它转化成相应的遥感水深数据口] ， 因此，本文将重点讨论和建立一种能将遥感信息与 煤矿沉陷区水域水深相吻合的回归计算模型． 1 研究区概况 研究区域为淮南矿区潘一煤矿开采后的地表 积水沉陷区．潘一煤矿1 9 8 3 年1 2 月2 6 日正式建 成投产，1 9 9 6 年达到设计能力3 ．0M t ／a ，之后矿井 产量一直稳定在3 ．0M t ／a ．潘一煤矿1 9 8 7 年积水 沉陷区面积0 ．2 0 4k m 2 ，1 9 9 2 年增加到0 ．6 7 87 k m 2 ，1 9 9 8 年为1 ．8 5 53k m 2 ，其发展速度之快，令 人瞩目．积水沉陷区主要接受泥河 淮河支流 、地 下水和大气降水补给，以前二者为主． 2数据处理 2 ．1 数据获取 积水沉陷区的水深实地采集时间为2 0 0 5 年4 月5 日～4 月1 5 日，遥感数据采用2 0 0 5 年4 月6 日L A N D S A T 5 T M 影像数据 轨道号1 2 1 3 7 ， T M l ～T M 5 ，T M 7 ，当日淮南天气状况良好，云量 很小，得到的遥感图像资料和实测水深资料可以认 为是准同步的，将作为进一步分析的数据源． 2 ．2 遥感数据预处理 2 ．2 ．1 几何校正 利用潘一煤矿1 1 00 0 0 的实测地形图，坐标 系统为北京5 4 坐标系，通过在地形图和T M 影像 上选取同名地物控制点来对遥感影像进行几何精 校正．在E N V l 4 ．2 运行环境下，在T M 影像和实 测地形图上寻找同名不变地物点4 0 个，采用二次 多项式进行坐标转换，用立方卷积法进行灰度重采 样，误差控制在0 ．5 个像元以内． 2 ．2 ．2 大气校正 由于大气散射和吸收的影响，导致遥感成像过 程中辐射能量被削弱，为提高遥感信息的提取和参 数反演的精度，因此在利用遥感图像资料之前必须 对图像进行大气校正．本文采用文献[ 8 ] 的大气校 正模式． 定标及辐射校正是大气校正的前提，L A N D S A T 5 预先测出了各波段的辐射亮度值L 。和记录 值D N 之间的校正增量系数g a i n 和校正偏差量 o f f s e t ，测出它们之间呈线形关系，故其定标采取 线形公式进行转换 L d D N g a i n o f f s e t ， 1 式中g a i n ，o f f s e t 可以从影像的头文件中读取． 辐射校正旨在计算卫星所观测到的行星反射 率 未经过大气校正的反射率 ，假设条件是天空辐 照度各向同性，且地面是朗伯体反射，则行星反射 率R ， R 一讯蕊老专而， 2 式中E o A 为大气层外相应波长的太阳光谱辐照 度；W A 为各波段的波段宽度 参数值可从 L A N D S A T 5 相关资料中获取 ． 影像经过辐射校正后，根据成黑体减法 d a r k o b j e c ts u b t r a c t i o n [ 8 1 模型理论，程辐射反射率P 。l 是由影像直方图下限值决定的，该直方图的下限值 被假定为程辐射反射率和暗体反射率之和，程辐射 率I D 曲与程辐射L 。的关系如下 丌L 。 阳一爵丽赢’ 3 式中0 为太阳天顶角 影像头文件中读取 ． 大气校正假设1 地面朗伯面反射和天空辐 照度各向同性；2 天空晴朗无云，不考虑大气的折 射、湍流和偏振；3 忽略大气的多次散射以及邻近 像元漫射过程．地物表面反射辐照度与太阳辐照度 和地表反射率的函数关系为 r R g x [ E o A C O S 眈e - r ““0 z E d ] ，小 L n 一 丌 ，、‘士／ 万方数据 第4 期 武彦斌等淮南积水沉陷区水深遥感定量分析 式中R n 为地表反射率；E “为天空光漫散射到地 表的光谱辐照度；n 为大气光学厚度，由文献[ 9 ] 的方法计算得到． ‘ 大气衰减后的地表反射辐照度与大气本身的 程辐射L 。之和，为卫星传感器接收到的光谱辐射 亮度L i x L i 一L n e - r x “”％ L p ， 5 式中0 ，为卫星传感器观测角，L A N D S A T 5 的观 测角为兀／2 ． 由式 4 ， 5 得到地表反射率 ％一毒l 讯万焉南专赢j ’ 6 式中L e l x “”吒为大气的下行倾斜透过率； T 。。 e 、““气为大气沿观测方向的上行透过率． 如只考虑一次散射，则整层大气的散射率为 1 一L 。，不考虑大气的吸收，E d 、近似为 ． 1 一T 。x E o A ，“ ／2 , d 1 一丁一‘ u7 2 ．2 ．3 分离积水沉陷区水体 以淮南潘一煤矿东采区形成的积水沉陷区为 例，该处沉陷区与泥河 淮河支流 连成一体，积水 沉陷区的提取必须结合研究区域的地质和采矿资 料．在A r c G I S 平台中对校正后的T M 5 波段进行 水迹线提取，并对提取的水迹线建立掩膜 m a s k ， 并将该掩膜图像与T M 图像的其他波段相乘，以 去除其它各波段的陆上信息，最后与井上井下对照 图做叠加，提取出积水沉陷区．为了更好地与遥感 影像套合，需要将水深采样点与遥感影像叠加、采 样数据的加密与生成等值线图．图1 是潘一煤矿东 采区形成的积水沉陷区实测水深等值线图． 图1实测水深等值线 F i g ．1 C o n t o u r so fa c t u a l l ym e a s u r e dw a t e rd e p t h 2 ．3 水深值和遥感数据的分析 在E N V l 4 ．2 平台中，将截取并掩膜后的 各波段的反射率值输出为A S C I I x y z 文件格 式，利用E N V I 的二次开发语言I D L 编写程序，使 水深实测点的水深值与其像元反射率值匹配． 对实测点的水深值和各波段相应的像元反射 率值进行了相关分析，旨在选取较好的水深反演因 子．由于波段比值法可使水深差异在舍去其背景的 条件下提高对比度m ] ，因此还进行了波段组合运 算与水深值的相关性分析，结果见表1 ． 表1波段及波段组合的反射率与 水深值的相关系数 T a b l e1C o r r e l a t i o nc o e f f i c i e n t sb e t w e e nr e f l e c t a n c eo f b a n d s c o m b i n e db a n d s a n dr e l e v a n tw a t e rd e p t h s 3 定量模型建立 从表1 可见，各波段及波段组合的反射率与水 深值之间的相关性．研究中尝试用与水深值相关性 最好的T M 4 波段，建立水深一反射率回归模型．以 像元的地表反射率为自变量，水深值为因变量． 在纯水中、底质均一、大气条件较好时，可以利 用将遥感图像的灰度值与光在水中的衰减程度 水 深 之间的线性关系式，推算出水深值．但是，实际 中若光可以从太阳到水体内至水底，再到卫星接收 设备，则水体的辐射总量主要由大气反射、水面反 射与折射、水中悬浮物质反射、散射和水底反射组 成．假设影响水体辐射量因素 底质、水生生物、含 沙量、水面波浪状况、水体上方大气、卫星传感器正 常 在空间上均匀变化，则水深增大而水体辐射强 度减弱．回归模型的建立选用指数函数e 一来作为 回归函数的主体部分，再增加2 项幂函数的修正． 由于采用水深实测点的水深一反射率回归关系 代替整个区域，进行反演水深运算，以局部代替整 体，不能完全反映全部信息；力口上目前未见煤矿沉 陷区水体波谱测试的详细报道，因此，根据研究区 域的平均水深值，以实测水深值2m 为阈值，模型 分为浅水区域 水深值小于2m 和深水区域 水深 值2 ～6m ． 浅水区域反演模型 H 一0 ．6 5 9 2 R 卧 0 ．4 6 0 4 2 0 ．9 3 7 8 R 。 0 ．4 6 0 4 0 ．6 5 7 9 ， D 0 ．4 2 4 1 e [ 12 4 1 8 ‘H “。 H Ⅲ1 2 0 6 ] 一 ，Q 、 3 5 0 ．4 1 4 6 ． 深水区域反演模型为 H 一0 ．4 6 0 4 2 5 2 万方数据 5 4 0中国矿业大学学报第3 6 卷 H 一一0 ．6 2 3 1 R 。。 0 ．4 1 5 7 2 1 ．0 2 9 4 R n 0 ．4 1 5 7 0 ．0 6 3 2 ， D 一0 ．8 8 4 1 e [ _ 1 删H 2 5 舢1 0 0 5 ] 一 9 H 一0 ．1 4 5 7 3 H 一0 ．4 1 5 7 2 2 54 4 0 ．1 1 4 6 ． 式中R n 为像元的地表反射率；H 为中间变量；D 为反演水深值． 根据统计学原理，利用给定类型的曲线回归， 数据越集中，估算精度越高，因此增加中间变量 1 1 0 ．6 5 92 R 一0 ．4 6 04 2 ．0 ．9 3 78 x 2 5 2 ．0 螽 冀。．。 0 ．5 ’；。i i ； l 。4 ’2 7 5 2 _ 4 1 H 4 - 6 0 ＼D o ．4 2 4l e Ⅲ4 1m 5 即2 ⋯“． ．4 6 04 3 ／3 5 _ 嬲氛‘＼≮ 0L 』J ▲J L 0 ．3 50 ．4 00 ．4 50 ．5 00 ．5 50 6 00 ．6 5 T M 4 反射率 a 浅水区域 H ；式 8 ， 9 中，增加第一修正项为一幂函数形 式的非线性修正，利用幂函数幂次越高，函数值变 化越剧烈的特点可以在一定程度上降低由于以局 部区域的回归方程代替整体时水深值产生的偏差； 增加第二修正项的目的是保证回归方程增加第一 修正项后仍然为一单调曲线，使水深和反射率保持 一一对应关系． 图2 a ，b 分别反映相应水深值与T M 4 波段像 元的地表反射率之间对应关系，其中纵轴表示水 深值，单位m ；横轴表示T M 4 的地表反射率． H 0 ．6 2 3I R g 04 4 57 2 1 ．0 2 94 x l O0 ．1 50 ．2 00 ．2 50 ．3 00 ．3 50 ．4 0 T M 4 反射率 b 深水区域 图2水深值与T M 4 反射率关系 F i g ．2 C o r r e l a t i o nb e t w e e nr e f l e c t a n c eo fT M 4a n dw a t e rd e p t h 图3 a ，b 分别表示实测水深等值面和反演水深 值等值面． 为进一步说明反演精度，选取独立的浅水水深 检查点2 0 个、深水水深检查点2 5 个，验证反演公 式的精度．表2 数据显示水深在小于2m 时，平均 绝对误差0 ．1 5 91i n ，平均相对误差1 3 ．3 6 7 ％；表3 中数据显示水深在2 ～6m 范围内，平均绝对误差 0 ．5 7 18m ，平均相对误差1 5 ．0 2 ％．表2 ，3 分别是 浅水、深水区域检查点的模拟深水误差分析表． 图3水深等值面 F i g ．3 I s o s u r f a c eo fw a t e rd e p t h 表2浅水区域误差 T a b l e2E r r o ro fr e t r i e v e ds h a l l o ww a t e rd e p t h 点号实测值／m反演值／m绝对误差／m 相对误差／％点号实测值／m反演值／m绝对误差／m 相对误差／％ 10 ．4 0 10 ．3 4 240 ．0 5 861 4 ．6 11 11 ．3 5 31 ．5 4 220 ．1 8 921 3 ．9 9 20 ．5 0 70 ．4 3 150 ．0 7 551 4 ．8 91 21 ．4 0 81 ．2 1 600 ．1 9 20 1 3 ．6 4 30 ．5 2 60 ．4 5 250 ．0 7 351 3 ．9 8 1 3 1 ．5 0 21 ．6 4 53 0 ．1 4 331 0 ．5 4 4 0 ．6 5 40 ．7 5 480 ．1 0 08 1 5 ．4 11 41 ．5 4 21 ．6 9 860 ．1 5 661 0 ．1 5 50 ．7 5 10 ．8 5 450 ．1 0 351 3 ．7 81 51 ．6 0 31 ．8 1 670 ．2 1 371 3 ．3 3 60 ．8 5 10 ．7 3 480 ．1 1 621 3 ．6 61 61 ．7 0 11 ．4 5 940 ．2 4 161 4 ．2 0 70 ．8 6 20 ．9 9 040 ．1 2 841 4 ．9 01 71 ．8 1 41 ．5 7 340 ．2 4 061 3 ．2 6 80 ．9 5 11 ．0 7 380 ．1 2 281 2 ．9 l1 81 ．8 6 11 ．6 9 680 ．1 6 421 0 ．8 2 91 ．1 0 11 ．2 5 860 ．1 5 761 4 ．3 1 1 9 1 ．9 1 42 ．1 8 720 ．2 7 32 1 4 ．2 8 1 01 ．2 0 21 ．0 2 790 ．1 7 411 4 ．4 9 2 01 ．9 5 01 ．6 9 270 ．2 5 731 3 ．1 9 5 O 5 O 5 O 5 O 5 0 0 0 m ／遥聪簧 万方数据 第4 期 武彦斌等淮南积水沉陷区水深遥感定量分析5 4 1 表3深水区域误差 T a b l e3E r r o ro fr e t r i e v e dd e e pw a t e rd e p t h 点号实测值／m反演值／m绝对误差／m 相对误差／％点号实测值／m反演值／m绝对误差／m 相对误差／％ 1 2 ．2 0 02 ．5 1 690 ．3 1 691 4 ．4 1 1 4 3 ．9 0 13 ．2 9 790 ．6 0 311 5 ．4 6 22 ．2 5 32 ．5 8 050 ．3 2 751 4 ．5 41 54 ．0 0 53 ．4 1 470 ．5 9 031 4 ．7 4 32 ．3 6 42 ．0 1 970 ．3 4 431 4 ．5 61 64 ．2 0 24 ．8 2 580 ．6 2 381 4 ．8 5 4 2 ．6 0 12 ．1 4 810 ．4 5 291 7 ．4 1 】7 4 ．4 0 15 ．0 5 950 ．6 5 851 4 ．9 6 52 ．6 5 12 ．1 7 480 ．4 7 621 7 ．9 61 84 ．5 0 13 ．7 7 120 ．7 2 981 6 ．2 1 62 ．7 0 13 ．1 1 860 ．4 1 761 5 ．4 61 94 ．5 4 05 ．1 9 220 ．6 5 221 4 ．3 7 72 ．9 0 53 ．3 5 490 ．4 4 991 5 ．4 9 2 0 4 ．7 9 34 ．0 2 540 ．7 6 761 6 ．0 1 83 ．0 1 12 ．5 3 480 ．4 7 6215 ．8 12 1 5 ．0 0 8 5 ．7 4 330 ．7 3 531 4 ．6 8 93 ．1 5 43 ．5 9 160 ．4 3 761 3 ．8 72 25 ．2 3 44 ．4 7 220 ．7 6 181 4 ．5 5 1 03 ．2 0 23 ．6 9 420 ．4 9 221 5 ．3 72 35 ．3 0 95 ．9 8 760 ．6 7 861 2 ．7 8 1 l 3 ．4 1 43 ．8 4 37 0 ．4 2 97 1 2 ．5 9 2 4 5 ．8 5 】 5 ．0 1 82 0 ．8 3 281 4 ．2 3 1 23 ．7 1 23 ．1 5 870 ．5 5 331 4 ．9 02 56 ．0 1 55 ．1 0 010 ．9 1 491 5 ．2 1 1 33 ．8 1 24 ．3 8 320 ．5 7 121 4 ．9 8 遥感反演值和实测水深值之问存在误差的主 要原因有 1 部分水域含有大量悬浮泥沙信息，悬浮泥 沙对入射光的后向散射会影响水体的总辐射量； 2 由于T M 影像分辨率的限制，水陆结合处 的像元含有部分陆地信息； 3 淮南潘一矿附近的积水沉陷区内含有绿 藻、蓝藻等藻类，藻类物质中的叶绿素对水体反射 率曲线会产生影响； 4 实测水深图受测点数量、分布的限制，个别 小区域的水深值也会不准确． 4 结论 1 通过分析T M 影像各波段及波段组合的反 射率与水深值的相关系数，建立了积水沉陷区的 T M 影像反射率与水深值的反演关系式．以相关系 数最大的T M 4 波段所建立的水深遥感反演关系 式，以水深2i n 为阈值，水深在小于2r n 时，平均 绝对误差0 ．1 5 91m ，平均相对误差1 3 ．3 6 7 0o ；水 深在2 ～6m 范围内，平均绝对误差0 ．5 7 18m ，平 均相对误差1 5 ．0 2 ％．从反演值与实测值的误差值 来看，模型对研究水域的整体地形信息与遥感光谱 信息之间关系描述得较好，可以为矿区积水沉陷区 综合治理提供一定精度的数据。 2 水深已知数据的可靠性在很大程度上影响 水深计算，如果其能得到保证，则可以进一步提高 遥感水深探测的可靠性和实用性． 3 水体总辐射量的影响因素为多方面的，各 因素与水深值问的定量关系，需要做进一步研究． 参考文献 [ 1 ]B I E R W O R T HPN ，L E ET ，B U R N ER 。S h a l l o w s e a f l o o rr e f l e c t a n c ea n dw a t e rd e p t hd e r i v e dh yu n m i x i n gm u l t i s p e c t r a li m a g e r y [ J ] ．P h o t o g r a m m e t r i c E n g i n e e r i n ga n dR e m o t eS e n s i n g ，1 9 9 3 ，5 9 3 3 3 1 ～ 3 3 8 ． E 2 1M G E N G E I ，V ，S P I T Z E RRJ ．A p p l i c a t i o no fr e m o t e s e n s i n gd a t at om a p i n go fs h a l l o ws e a f l o o rn e a rb y N e t h e r l a n d s [ J ] ．I n t e r n a t i o n a lJ o u r n a lo fR e m o t e S e n s i n g ，1 9 9 1 ，5 7 5 4 7 3 4 7 9 ． [ 3 ]党福星，丁谦．利用多波段卫星数据进行浅海水深 反演方法研究[ J ] ．海洋通报，2 0 0 3 ，2 2 3 5 5 5 9 ． D A N GF u x i n g ，D I N GQ i a n ．At e c h n i q u ef o re x t r a c t i n gw a t e rd e p t hi n f o r m a t i o nf r o mm u l t is p e c t r a l s c a n n e rd a t ai nt h eS o u t hC h i n aS e a E J ] 。M a r i n eS e i e n c eB u l l e t i n ，2 0 0 3 ，2 2 3 5 55 9 ． E 4 3 黄家柱，尤玉明．长江南通河段卫星遥感水深探测试 验[ J ] ．水利科学进展，2 0 0 2 ，1 3 2 2 3 5 2 3 8 ． H U A N GJ i a z h u ，Y O UY u m i n g ．E x p e r i m e n to fw a t e rd e p t hs u r v e y i n gi nt h eN a n t o n gs e c t i o no { t h e Y a n g t z eR i v e r [ J ] ．A d v a n c e si nW a t e rS c i e n c e ，2 0 0 2 ， 1 3 2 2 3 5 ～2 3 8 ． C 5 ] 陈鸣，李士鸿，孔庆芬．卫星遥感长江口水域水深 [ J ] ．水利水运工程学报，2 0 0 3 2 6 1 6 4 ． C H E NM i n g ，L IS h i h o n g ，K NGQ i n gf e n ．F a t h o m i n gm e t h o db ys a t e l l i t er e m o t es e n s i n gf o rt h e Y a n g t z eR i v e re s t u a r yw a t e rb o d y [ J ] ．H y d r o S c i e n c e a n dE n g i n e e r i n g ，2 0 0 3 2 6 1 6 4 ． [ 6 ] 于瑞宏，刘廷玺，李畅游，等．干旱区草型湖泊悬浮固 体浓度及水深的遥感与分析[ J ] ．水利学报，2 0 0 5 ，3 6 7 8 5 38 6 2 ． Y UR u ih o n g ，I 。I UT i n g x i ，L IC h a n g y o u ，e ta 1 ．A n a l y s i so fs u s p e n d e ds e d i m e n tc o n c e n t r a t i o na n dw a t e rd e p t hi ng r a s st y p el a k e sb ym e a n so fr e m o t e s e n s i n g [ J ] ．J o u r n a lo fH y d r a u l i cE n g i n e e r i n g ，2 0 0 5 ， 3 6 7 8 5 3 8 6 2 ． 下转第5 5 3 页 万方数据 第4 期 孙爱军等基于时序的工业用水效率测算与耗水量预测 5 5 3 r 1 3 3B A T T E S EGE ，C O E L L ITJ ．Am o d e lf o rt e c h n i c a l i n e f f i c i e n c ye f f e c t si na s t o c h a s t i cf r o n t i e rp r o d u c t i o n f u n c t i o nf o rp a n e ld a t a [ J ] ．E m p i r i c a lE c o n o m i c s ， 1 9 9 5 ， 2 0 2 3 2 5 3 3 2 ． [ 1 4 3李子奈，叶阿忠．高等计量经济学[ M ] ．北京清华大 学出版社，2 0 0 0 ． [ 1 5 3 涂正革．我国大中型工业企业生产率与技术效率的 随机前沿模型分析[ D 3 ．武汉华中科技大学经济学[ 1 7 ] 院，2 0 0 5 ． [ 1 6 3 涂正革，肖耿．中国的工业生产力革命用随机前 上接第5 4 1 页 [ 7 3 [ 8 ] [ 9 ] 任明达，吕斯骅，张绪定．中国海岸卫星遥感解译 [ M ] ．北京海洋出版社，1 9 9 0 1 0 8 1 1 6 ． 刘小平，邓孺孺，彭晓鹃．基于T M 影象的快速大气 校正方法[ J ] ．地理科学，2 0 0 5 ，2 5 1 8 7 9 3 ． L I UX i a o - p i n g ，D E N GR u - r u ，P E N GX i a o j u a n ．A [ 1 0 3 f a s ta t m o s p h e r i cc o r r e c t i o nm e t h o db a s e do nT Mi m a g e r y [ J ] ．S c i e n t i aG e o g r a p h i c aS i n i c a ，2 0 0 5 ，2 5 1 8 7 9 3 ． G I L A B E R TMA ，C O N E S EC ，M A S E L L IF ．A n 沿生产模型对中国大中型工业企业全要素生产率增 长的分解及分析[ J ] ．经济研究，2 0 0 5 3 6 3 6 8 ． T UZ h e n g g e ，X I A OG e n g ．C h i n a ’Si n d u s t r i a lp r o d u c t i v i t yr e v o l u t i o n as t o c h a s t i cf r o n t i e rp r o d u c t i o n f u n c t i o na n a l y s i so fC h i n a ’Sl a r g ea n dm e d i u mi n d u s t r i a le n t e r p r i s e sd u r i