理论酮究 遥感信息 一种基于颜色恒常性理论的城市高分辨率 遥感影像阴影消除方法 徐秋红①,叶勤①,② (①同济大学测量与国土信息系,上海200092;②同济大学遥感与空间信息技术研究中心,上海200092) 摘要:阴影是遥感影像中普遍存在的一种现象,它的存在影响对被遮挡区域信息的获取,并且对影像的目视效 果造成不好的影响。因此,对阴影的去除在遥感影像信息提取中是一个需研究的问题,特别是对于具有大量高层建 筑物的特大城市区域。本文采用颜色恒常性的理论,来实现对高分辨率遥感影像上阴影的去除,利用上海市陆家嘴 区域的QuiekBird高分辨率影像对算法进行了实验,结果表明该方法对阴影的去除具有较好的效果,并且计算简便。 关键词:颜色恒常性;阴影检测;阴影去除 doi:10.3969/j.issn.1000—3177.2010.04.003 中图分类号:TP79 文献标识码:A 文章编号:1000—3177(2010)11O一0013一O4 A Method of Shadow Elimination from City High Resolution Remote Sensing Images Based on Colour Constancy XU Qiu-hong ̄,YE Qin①,0 (①Department of Surveying and Geo—informatics,Tongji University,Shanghai 200092; ②Research Center of Remote Sensing&Saptial Imformation Technology,Tongji University,Shanghai 200092) Abstract:The shadow is a widespread phenomenon in remote sensing images,and it brings about serious problems to the in— formation gathering in the oeeluded areas and has a harmful influence on visual effect of image.Therefore,the removal of the shadow from the images in remote sensing information extraction is a problem to be studied,especially in the large urban areas with a large number of high-rise buildings.This paper is based on the theory of color constancy to remove the shadow from high—resolution remote sensing images,and uses the algorithm to conduct experiments with high-resolution Quickbird image in the region of Lujiazui in Shanghai.The results show that this method has a good effect of eliminating the shadow and the calcu— lation is simple. Key words:colour constancy;shadow detection;shadow removal 1 前 言 阴影的产生是源于光线被物体的遮挡,一般阴 影区域的灰度值与对比度比非阴影区域要小。这影 响了影像的质量,使得阴影区信息量减少,难以判 读,并对影像的后期处理工作如图像匹配、模式识别 如何去除阴影是遥感影像处理领域一直所关注的问 题,特别是在特大城市中心区域,高层建筑物的阴影 使高分辨率遥感影像的质量受到很大影响。 阴影的去除实际上就是将阴影区域中像素点的 亮度和颜色,恢复到非阴影区域光照的效果,从而使 整幅图像具有视觉上的光照一致性。由于在遥感影 和地物提取等带来困难,甚至会导致错误的结果。 像上建筑物阴影只存在于图像的部分区域,因此阴 收稿日期:2009--07—31 修订日期:2009--09—21 基金项目:863项目(编号:2OO7AA12z178),“十一五”国家科技支撑计划项目“村镇信息快速采集与处理关键技术研究”课题支持 (课题编号:2OO6BAJO9Bo1)。 作者简介:徐秋红(1985 ̄),女,硕士研究生,主要研究方向:遥感与数字摄影测量、图像处理。 E-mail:xuaured@sina.com 13 遥感信息 理论研究 的变化,然而人对该表面的色彩的认知在一定范围 内保持恒定 ]。由于阴影区的光照条件与非阴影区 差异很大,是一非标准光照下的场景图像,阴影消除 的目的是将其变换到标准光照(即非阴影区光照)条 件下,这个过程中的关键是估计光源颜色。 基于这一理论,在光源颜色估计上,众多学者根 据采用的假设条件和前提的异同提出了多种颜色恒 影去除包括两方面的工作:阴影区域的检测以及阴 影去除。 关于阴影区域的检测,这些年来,人们提出了多 种阴影检测方法,这些方法大致可以分为两类:基于 物体几何特性和基于阴影性质。具体有基于直方 图、基于同态系统阴影检测技术、归一化处理、基于 阴影强度与几何特性以及基于彩色不变模型等方 法。在文献E61中利用阴影区域在亮度和色度上的 常计算算法。这些算法可归为两类:(1)依据图像底 特性,采用光谱比和Otsu(最大类问方差法)阈值分 割技术,在ⅥQ(美国国家电视标准委员会(N )支 持的一种彩色电视传输的标准模式)空间进行阈值 分割,再利用数学形态学来去除草地、树木等对建筑 物阴影的影响;并且文中给出这种方法与在HSV 和YCbCr空间进行建筑物阴影提取效果的比较,认 为在YIQ空间最为理想。因此本文在针对城市高 分辨率遥感影像的阴影提取时,也采用这种方法。 而对于阴影的去除,因为遥感影像上每个像元的 亮度或色彩都是太阳光照函数和地物反射函数的复 合函数,这一复合函数相当复杂,从理论上来讲,要完 全消除影像上的阴影、恢复阴影区域中地物的本来面 目几乎是不可能的。传统的图像增强方法,如直方图 均衡、同态滤波、归一化处理等,对改善影像的阴影都 有一定的作用,但处理后的影像阴影仍然明显。 针对这一问题,本文重点研究基于颜色恒常性 的方法来实现对高分辨率遥感影像中的阴影去除, 特别是特大城市中心区的建筑物阴影去除。在处理 中,对阴影区域利用颜色恒常性算法进行阴影消除 处理,而对非阴影区域不做任何处理。实验表明这 种处理方法不仅具有较好的阴影去除效果,而且阴影 去除是在阴影检测的基础上的,故对非阴影区域不做 处理,这保证了非阴影区域的信息不受损失,弥补了 常用的阴影去除算法处理结果的缺点。例如,基于直 方图均衡的方法来讲,消除了直方图处理结果对于大 范围、复杂地形影像不适用的缺点;基于区域弥补的 阴影去除算法,该算法用临近非阴影区域的像素来填 充阴影区域,因此必将导致阴影区图像细节的丢失, 特别是当阴影区域是复杂纹理表面或者阴影区域很 大时,这种算法将导致比较严重的图像失真。 2颜色恒常性理论 颜色恒常是指在照明发生变化的条件下人们对 物体表面的知觉趋于稳定的心理倾向。同一表面在 不同的照明条件下会产生不同的反射谱,人的颜色 机制确实能分辨这种由照明变化导致的表面反射谱 14一 层特征的算法进行光源估计;目前主要有Gray- World算法,Max-RGB算法和Gray-Edge算法等。 这类算法计算量低,计算速度快捷有效,但是每种算 法也要基于一定的成像假设,只有对符合条件的图 像结合才能够取得较精确的估计。(2)依据统计学 的方法进行光源估计;典型的代表是Color by Cor— relation算法,该算法精确度不高,且需要大量的先 验知识,但具有很好的通用性。 在上述研究基础上,G Finlayson和Trezzi提出 了一种颜色恒常算法Shades of Gray算法,这种方 法假设图像的阴影区域和非阴影区域的场景图像部 分都应满足Minkowski范式的颜色恒常算法口 。 此外,G Finlayson和Trezzi指出Gray-World和Max- RGB算法都包含在明可夫斯基范式(Minkows ̄ norm)内,明可夫斯基范式如式(1): …(厂(z, )) dxdy I 一 J 一定1一忌f【 1——丽 t= — 1 ( , 1J l(1l ) 式中e表示当前场景中的光源颜色,厂表示图 像各通道的灰度值,忌为比例系数,P可以取[1, o。)内任意整数。当P=1时,即为Gray-World算 法,当P—o。时,为Max-RGB算法。Minkowski 范式参数P决定了估计光源所使用到的各个算法的 权重。文献ElsJ中用此方法对普通场景图像做阴影 去除,实验效果表明p一6效果最好,并且相对而 言,该算法计算简便。 与普通场景图像相比,高分辨率遥感图像具有 成像范围大、涉及的地物类型复杂、光照条件复杂等 特点,特别在大城市的中心区域,其建筑物阴影会造 成相邻地物间极大的光照条件差异。而针对普通场 景图像适用的光源估计算法,在高分辨率遥感图像 上的适应性还需研究。因此本文将采用该算法对城 市中心区高分辨率遥感图像进行阴影去除,分析比 较该算法对高分辨率遥感图像阴影去除的效果。 圈匿Ⅱ日 ④得到去除阴影后的图像为 遥感信息 3 基于颜色恒常性理论的阴影去除方法 基于颜色恒常性的阴影去除方法是建立在阴影 检测的基础上的,阴影检测的结果直接影响阴影去 除效果的好坏。 news(i,_『)===g—h(i, )+h(i, ) (6) ⑤将处理后的每幅灰度图像news ( ,.『),ne— wSG( ,-f),newsB( ,J)按R,G,B中输出。 图1为算法流程图: 3.1阴影检测步骤 本文选用前言中提到的方法进行了高分辨率遥 感影像上的阴影检测,得到阴影区域为基于颜色恒 常性的阴影去除做准备,其检测过程、步骤如下: ①将基于RGB的彩色图像转换到YIQ空间, 并提取出其中表示亮度I和色调H的分量,YIQ中 的Y表示亮度,Q表示色调。 ②利用光谱比技术,得到(H+1)/(I+1)比值 图像。(H+1)/(I+1)比值图像能增强亮度低的阴 影区域,使得色调值增加。在(H+1)/(I+1)比值 图像中,阴影区域像素值比非阴影区域内的像素值 要高。 ③利用Otsu方法对比值图像进行分割。Otsu 法利用比值图像中的直方图自动确定分割阈值,像 素值大于阈值部分即为阴影。 ④基于建筑物阴影与绿化树及其阴影在形态上 的差异,对图像进行形态学处理。首先对分割后图 像进行中值滤波处理,消除噪声影响;然后选择结构 元对图像进行腐蚀运算,断开图像中细长部分;再利 用开启运算,将连通区域面积较小阴影部分去掉;最 后对图像进行膨胀,以得到建筑物阴影。 3.2阴影消除步骤 将原始彩色航空图像分解为S ( , ),s。( , ), 图1算法流程图 4实验结果与分析 本文在Windows XP操作系统下,利用Matlab 的编程平台做的实验。选择上海2005年上海埔东 Se( , )三幅灰度图像,利用颜色恒常性理论分别进 行阴影去除处理。具体步骤如下: ①根据阴影检测的成果将灰度图像分成阴影区 域和非阴影区域两幅的图像g(i,J)以及h(i, ) s(i,.『)一g(i, )+h( , ) (2) 陆家嘴地区QuickBird多光谱影像作为实验数据, 图象分辨率为2.4m、大小为246×384像素,上面具 有明显的建筑物阴影特征如图2。图3为利用3.1 中方法进行阴影区域检测的结果,其中白色区域为 阴影区域,黑色为非阴影区域。 ②对阴影区域和非阴影区域分别利用明可夫斯 基范式(Minkowski norm),即式(1),进行颜色恒常 性计算来求出光源颜色 1和 2。其中P值需根据 实际图像来给定。 ③将阴影区域变换到标准光照条件下,即为 g—b(i, )一g(i, )・(e1) (3) 同样非阴影区变换到标准光照下为 h b(i, )一h(i, ,)・(P2) (4) g—b(i, )和h—b(i,.『)场景光照是一致的,当阴 影区域变换到非阴影区域光照条件下时,即为 g—h( , )一g(i, )・(e1)-1(P2) (5) 图2输入图像 图3阴影检测结果 15 遥感信息 圈回圃圈 Contrast:= (8) 同时也对原始图象上非阴影区的地物影像进行 了统计,计算结果如表1: 表1亮度与对比度的比较 户 原图 2 4 6 8 lO 对比度 拉伸图 Brightness 69.94 81.O0 83.25 84.41 84.】3 83.28 66.42 9.17 图4 一2时结果 图5 p一4时结果 Contrast 25.93 7.8O 8.08 8.27 8.2O 8.17 从表1及实验结果图像可以分析得出: ①对于本文所用的方法,随着P取值的增大, 亮度与对比度增大,当P取6时,亮度、对比度值达 到峰值,随后户值增大,两者又开始变小。因此,当 P取6时去除阴影后的高分辨率遥感影像具有最好 的效果,这与普通场景图像的实验结果一致。 ②从目视效果来判断,采用颜色恒常算法阴影 去除,阴影区的颜色与非阴影区接近,没有出现大的 图6 p一6时结果 图7户一8时结果 色调偏移,同时阴影区内的细节地物也得到了好的 体现;P一2时阴影区处理后的效果最接近非阴影 区,而其他几幅图像阴影区与非阴影区,感觉差异较 大,这主要是因为,P取其他值时,两者亮度上的差 异相对较大(非阴影区的平均亮度为70,而处理后 的阴影区达到了81甚至以上)上。 ③在阴影与非阴影的边界出现较明显的痕迹,这 一方面是由于两者之间的灰度差异,另一方面也和阴 影检测算法检测出的边界有微小移位而造成的。 ④对阴影区域做对比度拉伸的图像虽然对比度 得到了提高(9.17),并且比其他方法的对比度值要 图8户一10时结果 图9对比度拉伸结果 大,但亮度比原图非阴影区低,而且从目视效果看, 出现了较大色调偏差。 本文利用3.2中的方法,对Minkowski范式参 数P取值为2、4、6、8和1o来实现阴影区域的去除, 实验结果如图4~图9所示。本文还对该图像的阴 影区域利用对比度反差线性拉伸方法来消除阴影区 域的影响,结果如图9。 为了对阴影消除结果做客观质量评估,本文对 阴影去除后各实验结果上阴影区的亮度(bright— 5 结束语 本文探讨了利用颜色恒常理论来实现对高分辨 率遥感影像阴影的去除,主要是高大建筑物阴影的 去除。采用光谱比和Otsu阈值分割技术在YIQ空 间进行阔值分割、再利用数学形态学去除细小阴影 来提取建筑物的阴影,然后基于颜色恒常性算法对 阴影区进行阴影去除。实验证明对于高分辨率遥感 ness)和对比度(contrast)进行了统计。 亮度反映了图像的明暗程度。其计算方法如式(7)。 卅 3 影像,该方法具有较好的阴影去除效果、不影响阴影 区域的细节特征、阴影区的色调与非阴影区基本保 持一致,弥补了常用去阴影方法的缺点。由于只对 阴影区进行处理,因此非阴影区域的信息不受任何 ∑∑∑S(i, ,愚) Bright舢 一旦 一 (7) 对比度反映图像的黑白反差,当对比度越大,图像 中的黑白反差也就越明显。式(8)为对比度计算方法。 影响。就算法本身而言,它与同态滤波等方法相比, 又具有较低的计算复杂度。 (下转第111页) l6 一 2010.4 GIS技斑 遥感信息 6 结束语 本文在对空间信息服务集成技术的现状与存在 问题的分析基础上,提出基于服务总线的空间信息 集成框架,并结合相关实例,进行了验证说明。该方 及服务路由动态选择机制,提高了服务集成的灵活 性、稳定性,对空间信息服务集成技术研究有着重要 的参考价值。在今后的研究工作中,我们将对服务 集成的安全性、事务性、服务节点间数据高效传输等 关键性问题进行进一步研究与探讨,以便更好满足 空间信息的普适化应用发展趋势要求。 法解决了现阶段空间信息服务集成中点对点模式的 弊端,支持多种异构平台的服务接入,采用异步通信 参考文献 1-13 王方雄,边馥苓.从GISystem到GIService:GISystem发展的必然趋势[J].华中师范大学学报:自然科学版,2004,38 (4):528—532. [23 乐鹏.语义支持的空间信息智能服务关键技术研究[D].武汉:武汉大学。2007。 a1.The enterprise service bus:Making service-oriented architecture [3] SCHMIDT M T,HU TCHISON B,LAMBROS P,etreal[EB/OL].http://www.research.ibm.corn/journal/sj/444/schmidt.html,2005—12—01. NS0N R.Understand enterprise service bus scenarios and solutions in service-oriented architecture[EB/OL].ht— [4] R0BItp://www.ibm.com/developerworks/webservices/library/ws—esb-scen2.html,2OO4—12—31. [5] I.Foster,C.Kesselman,and S.Tuecke.The anatomy of the grid:Enabling sealable virtual organizations.International[]]. Supercomputer Applications,2001,15(3). [6] 徐志伟,冯百明,李伟.网格计算技术[M].北京:电子工业出版社,2004. obus Allinace.http://www.globus.org [7] The Gl[8] 袁占亭,张秋余,翟志万.基于软件总线技术的软件工程方法[J].计算机工程与应用,2003(15):137—139. [9] 王凌云,李琦,等.Web Service与地理信息互操作[J].测绘科学,2004,29(1):38—41. (上接第16页) 参考文献 [1]Highnam R,Brady M.Model—based image enhancement of far infrared images[J].IEEE Transactions on Pattern Analysis and Machine Intelligence,1997,19(4):410—415. [2]Voicu L I,Myler H R,Weeks A R.Practical considerations on color image enhancement using homomorphic filtering[J]. Journal of Electronic Imaging,1997,6(1):108—113. [3]虢建宏,田庆久.一种IKONOS影像阴影检测与去除的方法口].遥感信息,2004(4):32—34. [4]Jiang C,Ward M O M.Shadow identification.computer vision and pattern recognition[C]//IEEE Computer Society Con— ference on.1992(6):606—612. [5] Salvador E,Cavallaro A,Ebrahimi T.Shadow identification and classification using invariant color models[C]//IEEE In— ternational Conference.2001,3(7—11):1545—1548. [6]喻红艳,李利军.彩色航空影像中的建筑物阴影提取[j].科学技术与工程,2007,7(7):1389—1393. [7]陈硕,沈模卫.颜色恒常理论及模型探索[J].心理科学,2003,26(2):215—218. [8]G Buchsbaum.A spatial processor model for object color perception[J].Franck Inst,1980(310):1--26. [9]E Land,J McCann.Lightness and retinex theory[J].Journal of the Optical ocSiety of America,1971,61(1):1—11. [1O]J van de Weijer,T Gevers,A Gijsen ̄.Edge-based color constancy[J].IEEE Trans on Image Processing,2007,16(9): 2207—2214. [11] G Finlayson,S Hordley,P Hube1.Color by correlation:A simple,unifying framework for color constancy[J].IEEE Trans.Pattern Ane1.Machine Intell,2001,23(11):1209—1221. [12]K Barnard,L Martin,B Funt,A Coath.A data set for color research[J].Color Research and Application,2002,27(3): 147—151. [13]王宁.图像的阴影检测与去除算法研究[D].北京:北京交通大学,2008. [14]江兴方,金龙,何隆华,陶纯堪.基于Retinex彩色图像增强及其闭合区域提取的研究[J].光子学报,2007,36(3):565—567.