CN107507151A - 一种多光谱遥感图像真实色彩还原方法及*** - Google Patents
一种多光谱遥感图像真实色彩还原方法及*** Download PDFInfo
- Publication number
- CN107507151A CN107507151A CN201710781623.0A CN201710781623A CN107507151A CN 107507151 A CN107507151 A CN 107507151A CN 201710781623 A CN201710781623 A CN 201710781623A CN 107507151 A CN107507151 A CN 107507151A
- Authority
- CN
- China
- Prior art keywords
- remote sensing
- sensing image
- pixel
- spectral remote
- pixel value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 58
- 238000001914 filtration Methods 0.000 claims description 16
- 230000011218 segmentation Effects 0.000 claims description 16
- 238000010606 normalization Methods 0.000 claims description 7
- 230000003595 spectral effect Effects 0.000 claims description 6
- 238000003707 image sharpening Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- VMXUWOKSQNHOCA-UKTHLTGXSA-N ranitidine Chemical compound [O-][N+](=O)\C=C(/NC)NCCSCC1=CC=C(CN(C)C)O1 VMXUWOKSQNHOCA-UKTHLTGXSA-N 0.000 abstract description 32
- 230000004927 fusion Effects 0.000 abstract description 8
- 230000000694 effects Effects 0.000 abstract description 5
- 230000000007 visual effect Effects 0.000 abstract description 5
- 239000003086 colorant Substances 0.000 abstract description 3
- 238000005286 illumination Methods 0.000 description 5
- 230000009467 reduction Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000001965 increasing effect Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011017 operating method Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/80—Geometric correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/188—Vegetation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
- G06T2207/10036—Multispectral image; Hyperspectral image
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Multimedia (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种多光谱遥感图像真实色彩还原方法及***。该方法包括:获取多光谱遥感图像;采用线性拉伸法对多光谱遥感图像进行拉伸,得到第一多光谱遥感图像;对多光谱遥感图像中各像素点绿色波段的像素值进行调整,采用线性拉伸法将调整后的多光谱遥感图像进行拉伸,得到第二多光谱遥感图像;根据多光谱遥感图像的各像素点的归一化植被指数,对第一多光谱遥感图像和第二多光谱遥感图像进行融合,对融合后的图像进行滤波、分段拉伸以及锐化,达到将多光谱遥感图像真实色彩还原的效果。本发明提供的多光谱遥感图像真实色彩还原方法及***能够改善遥感图像的视觉效果,使其更加接近地物的真实颜色以及更加接近人眼观察到的地物的颜色。
Description
技术领域
本发明涉及图像处理领域,特别是涉及一种多光谱遥感图像真实色彩还原方法及***。
背景技术
进入21世纪,高分辨率遥感卫星进入了蓬勃发展的时期,众多国家竞相发射高分辨率卫星,其中绝大部分是高空间高分辨率的多光谱卫星,且图像的深度都为16byte,也就是遥感图像的深度不再是8byte,这样的图像使得地物纹理更加丰富。但是,计算机显示图像都是8byte,这就需要将16byte图像转换为8byte的图像进行展示。通常的方法有线性拉伸、直方图均衡化、高斯拉伸、二次开方拉伸等等,这些方法在主流的遥感图像处理软件:PCI、ENVI、ERDAS等都有实现,这些拉伸方法造成地物纹理变少,虽然地物的亮度得到了一定的增强,但是地物亮度并不均匀,而且,色彩与人眼看到的地物本身的色彩相差甚远,因此,需要一种新的方法来增强图像亮度、保留图像纹理细节,最终使得图像中的地物接近人眼所看到的地物本来的色彩。
现有的图像拉伸方法虽然能在一定程度上改善遥感图像的视觉效果,但是难以达到与人眼观察地物的真实感相同的效果。
另外,还有一种匀光算法,因为通常情况下图像光照是比较均匀的,但是有时候由于多种因素的影响,图像光照并不均匀,为了改善光照均匀度,有人提出了改善光照均匀度的图像处理方法,这种方法有人称为匀光,例如GeoDogging软件,这种方法是针对光照的差异来改善图像,但并不具有使图像中的地物颜色接近人眼所看到的地物颜色的效果。
发明内容
本发明的目的是提供一种多光谱遥感图像真实色彩还原方法及***,能够改善遥感图像的视觉效果,使其更接近地物本身的实际颜色,更加接近人眼观察到的地物的颜色。
为实现上述目的,本发明提供了如下方案:
一种多光谱遥感图像真实色彩还原方法,所述方法包括:
获取多光谱遥感图像;
采用线性拉伸法将所述多光谱遥感图像由16byte图像转化为8byte图像,得到第一多光谱遥感图像;
根据公式Gafter=Gbefore·0.8+NIR·0.2对所述多光谱遥感图像中各像素点绿色波段的像素值进行调整,其中,Gafter为所述像素点调整后的绿光波段的像素值,Gbefore为所述像素点调整前的绿光波段的像素值,NIR为所述像素点的红外波段的像素值;
采用线性拉伸法将调整后的所述多光谱遥感图像由16byte图像转化为8byte图像,得到第二多光谱遥感图像;
计算所述多光谱遥感图像的各像素点的归一化植被指数;
判断所述多光谱遥感图像中像素点的归一化植被指数是否大于零;
如果是,则将所述像素点的像素值确定为所述第二多光谱遥感图像中所述像素点对应的像素值;
如果否,则将所述像素点的像素值确定为所述第一多光谱遥感图像中所述像素点对应的像素值。
可选的,在所述获取多光谱遥感图像之后,所述采用线性拉伸法将所述多光谱遥感图像由16byte图像转化为8byte图像,得到第一多光谱遥感图像之前,还包括:
采用快速大气校正算法对所述多光谱遥感图像进行大气校正。
可选的,在确定所述像素点的像素值之后,还包括:
对所述多光谱遥感图像进行高斯滤波,得到滤波后的所述多光谱遥感图像。
可选的,在所述对所述多光谱遥感图像进行高斯滤波之后,还包括:
对滤波后的所述多光谱遥感图像进行分段拉伸。
可选的,所述对滤波后的所述多光谱遥感图像进行分段拉伸,具体包括:
判断滤波后的所述多光谱遥感图像像素点的像素值是否大于等于0且小于等于90,得到第一判断结果;
当所述第一判断结果表示是时,则根据公式y=x·(90-30)/90+30对所述像素点的像素值调整;
判断滤波后的所述多光谱遥感图像像素点的像素值是否大于90且小于180,得到第二判断结果;
当所述第二判断结果表示是时,则保持所述像素点的像素值不变;
判断滤波后的所述多光谱遥感图像像素点的像素值是否大于等于180且小于等于255,得到第三判断结果;
当所述第三判断结果表示是时,则根据公式y=(x-180)·(245-180)/(255-180)+180对所述像素点的像素值调整;
其中,x为所述像素点拉伸前的像素值,y为所述像素点拉伸后的像素值。
可选的,所述方法还包括:
采用USM锐化算法对拉伸后的所述多光谱遥感图像的图像纹理和边缘清晰度进行锐化。
本发明还提供了一种多光谱遥感图像真实色彩还原***,所述***包括:
多光谱遥感图像获取模块,用于获取多光谱遥感图像;
第一拉伸模块,用于采用线性拉伸法将所述多光谱遥感图像由16byte图像转化为8byte图像,得到第一多光谱遥感图像;
像素调整模块,用于根据公式Gafter=Gbefore·0.8+NIR·0.2对所述多光谱遥感图像中各像素点绿色波段的像素值进行调整,其中,Gafter为所述像素点调整后的绿光波段的像素值,Gbefore为所述像素点调整前的绿光波段的像素值,NIR为所述像素点的红外波段的像素值;
第二拉伸模块,用于采用线性拉伸法将调整后的所述多光谱遥感图像由16byte图像转化为8byte图像,得到第二多光谱遥感图像;
归一化植被指数计算模块,用于计算所述多光谱遥感图像的各像素点的归一化植被指数;
归一化植被指数判断模块,用于判断所述多光谱遥感图像中像素点的归一化植被指数是否大于零;
像素值第一确定模块,用于当所述多光谱遥感图像中像素点的归一化植被指数大于零时,将所述像素点的像素值确定为所述第二多光谱遥感图像中所述像素点对应的像素值;
像素值第二确定模块,用于当所述多光谱遥感图像中像素点的归一化植被指数小于等于零时,将所述像素点的像素值确定为所述第一多光谱遥感图像中所述像素点对应的像素值。
可选的,所述***还包括:
大气校正模块,用于采用快速大气校正算法对所述多光谱遥感图像进行大气校正。
高斯滤波模块,用于对所述多光谱遥感图像进行高斯滤波,得到滤波后的所述多光谱遥感图像;
分段拉伸模块,用于对滤波后的所述多光谱遥感图像进行分段拉伸。
可选的,所述分段拉伸模块,具体包括:
第一判断单元,用于判断滤波后的所述多光谱遥感图像像素点的像素值是否大于等于0且小于等于90,得到第一判断结果;
第一拉伸单元,当所述第一判断结果表示是时,根据公式y=x·(90-30)/90+30对所述像素点的像素值调整;
第二判断单元,用于判断滤波后的所述多光谱遥感图像像素点的像素值是否大于90且小于180,得到第二判断结果;
第二拉伸单元,用于当所述第二判断结果表示是时,保持所述像素点的像素值不变;
第三判断单元,用于判断滤波后的所述多光谱遥感图像像素点的像素值是否大于等于180且小于等于255,得到第三判断结果;
第三拉伸单元,用于当所述第三判断结果表示是时,根据公式y=(x-180)·(245-180)/(255-180)+180对所述像素点的像素值调整;
其中,x为所述像素点拉伸前的像素值,y为所述像素点拉伸后的像素值。
可选的,所述***还包括:
图像锐化模块,用于采用USM锐化算法对拉伸后的所述多光谱遥感图像的图像纹理和边缘清晰度进行锐化。
根据本发明提供的具体实施例,本发明公开了以下技术效果:本发明提供的多光谱遥感图像真实色彩还原方法及***,首先将多光谱遥感图像进行拉伸,得到第一多光谱遥感图像,然后,对多光谱遥感图像中各像素点绿色波段的像素值进行调整,并将调整后的多光谱遥感图像进行拉伸,得到第二多光谱遥感图像,之后,根据各像素点的归一化植被指数对第一多光谱遥感图像和第二多光谱遥感图像进行融合,并对融合后的图像进行滤波、分段拉伸以及锐化,实现了对多光谱遥感图像色彩的还原,改善了遥感图像的视觉效果,使其更接近地物本身的实际颜色,更加接近人眼观察到的地物的颜色。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例一的多光谱遥感图像真实色彩还原方法的流程示意图;
图2为本发明实施例二的多光谱遥感图像真实色彩还原方法的流程示意图;
图3为多光谱遥感图像由16byte图像拉伸为8byte图像的图像显示图;
图4为采用本发明实施例提供的方法对多光谱遥感图像处理后的图像显示图;
图5为本发明实施例多光谱遥感图像真实色彩还原***结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种多光谱遥感图像真实色彩还原方法及***,能够改善遥感图像的视觉效果,使其更接近地物本身的实际颜色,更加接近人眼观察到的地物的颜色。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明实施例一的多光谱遥感图像真实色彩还原方法的流程示意图,如图1所示,本发明提供的多光谱遥感图像真实色彩还原方法的操作步骤如下:
步骤101:获取多光谱遥感图像;
步骤102:采用线性拉伸法将所述多光谱遥感图像由16byte图像转化为8byte图像,得到第一多光谱遥感图像;
步骤103:根据公式Gafter=Gbefore·0.8+NIR·0.2对所述多光谱遥感图像中各像素点绿色波段的像素值进行调整,其中,Gafter为所述像素点调整后的绿光波段的像素值,Gbefore为所述像素点调整前的绿光波段的像素值,NIR为所述像素点的红外波段的像素值;
步骤104:采用线性拉伸法将调整后的所述多光谱遥感图像由16byte图像转化为8byte图像,得到第二多光谱遥感图像;
步骤105:计算所述多光谱遥感图像的各像素点的归一化植被指数;
步骤106:判断所述多光谱遥感图像中像素点的归一化植被指数是否大于零;
步骤107:如果所述多光谱遥感图像中像素点的归一化植被指数大于零,则将所述像素点的像素值确定为所述第二多光谱遥感图像中所述像素点对应的像素值;
步骤108:如果所述多光谱遥感图像中像素点的归一化植被指数小于等于零,则将所述像素点的像素值确定为所述第一多光谱遥感图像中所述像素点对应的像素值。
优选的,在步骤101和步骤102之间,还包括:采用快速大气校正算法对所述多光谱遥感图像进行大气校正。大气校能够消除大气和光照等因素对地物反射的影响,而且快速大气校正算法不需要外界更多的参数,能够自动从图像中收集不同物质的波谱信息,具有较高的精度。
优选的,在步骤108之后还包括:对所述多光谱遥感图像进行高斯滤波,得到滤波后的所述多光谱遥感图像;对滤波后的所述多光谱遥感图像进行分段拉伸;采用USM锐化算法对拉伸后的所述多光谱遥感图像的图像纹理和边缘清晰度进行锐化。其中,滤波窗口可以设置为5×5大小,USM锐化算法的权重可以设置为0.8,图像高斯模糊化小窗口半径为3,阈值设置为2。
优选的,具体的分段拉伸方法为:判断滤波后的所述多光谱遥感图像像素点的像素值是否大于等于0且小于等于90,得到第一判断结果;
当所述第一判断结果表示是时,则根据公式y=x·(90-30)/90+30对所述像素点的像素值调整;
判断滤波后的所述多光谱遥感图像像素点的像素值是否大于90且小于180,得到第二判断结果;
当所述第二判断结果表示是时,则保持所述像素点的像素值不变;
判断滤波后的所述多光谱遥感图像像素点的像素值是否大于等于180且小于等于255,得到第三判断结果;
当所述第三判断结果表示是时,则根据公式y=(x-180)·(245-180)/(255-180)+180对所述像素点的像素值调整;
其中,x为所述像素点拉伸前的像素值,y为所述像素点拉伸后的像素值。
本发明提供的多光谱遥感图像真实色彩还原方法,首先将多光谱遥感图像进行拉伸,得到第一多光谱遥感图像,然后对多光谱遥感图像中各像素点绿色波段的像素值进行调整,并将调整后的多光谱遥感图像进行拉伸,得到第二多光谱遥感图像,之后根据各像素点的归一化植被指数对第一多光谱遥感图像和第二多光谱遥感图像进行融合,并对融合后的图像进行滤波、分段拉伸以及锐化,实现了对多光谱遥感图像色彩的还原,使其更接近人眼观察到的地物的颜色。
作为本发明的又一实施例,如图2所示,具体流程步骤如下:
步骤201:读入四个通道的多光谱遥感图像;四个通道为红、绿、蓝、红外通道。
步骤202:图像大气校正。因为大气对地物反射光具有散射、吸收等作用,引起地物在不同波段的反射光发生不同程度的衰减,因此需要对遥感图像进行大气校正,从而消除大气和光照等因素对地物反射的影响。我们采用快速大气校正算法(QUAC),因为快速大气校正算法不需要外界更多的参数,能够自动从图像中收集不同物质的波谱信息,从而获得经验值实现大气校正,而且,校正精度与校正效果较好。
步骤203:通过计算归一化植被指数(NDVI)来识别植被地物。利用大气校正后的图像计算NDVI,NDVI是用来判断植被与非植被信息的有效指标,经过快速大气校正后的图像得到的NDVI能够更加准确的反映地物自身的光谱特征,是区分植被与非植被的有效指标。NDVI的计算公式为:
其中,NIR为近红外波段,R为红外波段。
NDVI能够在很大程度上消除与太阳高度角、卫星观测角、地形、阴影,以及大气散射与吸收等对地物辐射变化的影响。NDVI结果限定在[-1,1]之间,NDVI是植被生长状态及植被覆盖度的最佳指示因子。当NDVI为负值时表示地面覆盖为云、水、雪等,当NDVI值为0表示有岩石或裸土等,NIR和R近似相等,当NDVI为正值时表示有植被覆盖,且随覆盖度增大而增大。在遥感图像中植被往往占很大的比例,城市中植被所占比例远远低于农村,但通常也占有40%甚至更高,通过NDVI来识别植被,并在后续算法中增强这部分信息。
步骤204:16byte图像转换成8byte图像。拉伸方法采用线性拉伸。目的是将16byte数据从范围为0到65535压缩到范围为0到255的8byte数据,其计算公式为:
y=255*(x-minVal)/(maxVal-minVal)
其中,x为16byte图像的像素,minVal为16byte图像中像素最小值,maxVal为16byte图像中像素最大值。
图3为多光谱遥感图像由16byte图像拉伸为8byte图像的图像显示图,如图3所示,这样,16byte图像就转换成8byte图像,从而能够显示在屏幕上。
步骤205:增强植被信息的波段运算。就是将红外通道按照一定的比例增加到绿色通道中,这样做的目的是增加图像中植被的绿色信息,使植被更加鲜艳,更加接近人眼所看到的植被的色彩。植被在红外波段比在绿波段有更强的反射特征,也就是说植被在红外波段图像中的像素值更亮,为了提高绿色波段的亮度,按照一定的比例将红外波段增加到绿色波段中:
Gafter=Gbefore·0.8+NIR·0.2
其中,Gafter为所述像素点调整后的绿光波段的像素值,Gbefore为所述像素点调整前的绿光波段的像素值,NIR为所述像素点的红外波段的像素值,其加权比例为0.8和0.2。
因为植被在红外波段反射强,所以通过波段运算后的植被在图像中表现为更为鲜艳的绿色,而不是暗绿色。但是其他地物的光谱也被改变了,其他地物在图像中偏紫色,消除这种现象需要在第七步进行处理。
步骤206:波段运算后,将16byte图像转换成8byte图像,如图4所示。拉伸目的和拉伸方法同步骤204。
步骤207:图像融合。也就是将步骤204和步骤206得到的拉伸图像进行融合,融合时要根据NDVI值的大小来判断,如果NDVI大于0值,图像的像素值采用步骤206图像中的像素,否则就采用步骤204图像中的像素。这样就得到了初步真实色彩还原结果。
步骤208:滤波消除伪边缘。步骤207中两个图像依据NDVI值的大小来选择不同图像的像素,从而实现图像融合,这个过程中必然会出现两个图像像素之间的突然变化,从而导致伪边缘。为了消除伪边缘,需要采用滤波方法进行滤波,采用高斯滤波,滤波窗口设置为5×5大小。
步骤209:对融合后的图像进行分段拉伸。第一个目的是调整亮度,即将亮度低的地物的亮度进行一定的提高,将亮度高的地物的亮度整体降低,从而使得整个图像的亮度均匀分布在一个合理的范围之内,第二个目的是保持相邻图像之间在像素值上的一致性。设定分段拉伸图像最小值为30,最大值为245,将融合后的图像像素值分成三段,第一段为0到90,第二段为90到180,第三段为180到255,分段拉伸公式为:
y=x·(90-30)/90+30,{0≤x≤90}
y=x,{90≤x≤180}
y=(x-180)·(245-180)/(255-180)+180,{180≤x≤255}
步骤210:对图像纹理与边缘清晰度的锐化。因为图像拉伸、融合和滤波等操作,使得图像细节信息部分丢失,为了弥补这样的损失,采用图像锐化算法对图像进行操作,使得图像边缘更加清晰,纹理细节更加丰富,从而得到最终的真实色彩还原图像。我们这里采用了USM锐化算法,这个算法的权重设置为0.8,图像高斯模糊化小窗口半径为3,阈值设置为2。
本发明提供的多光谱遥感图像真实色彩还原方法,首先将多光谱遥感图像进行拉伸,得到第一多光谱遥感图像,然后对多光谱遥感图像中各像素点绿色波段的像素值进行调整,并将调整后的多光谱遥感图像进行拉伸,得到第二多光谱遥感图像,之后根据各像素点的归一化植被指数对第一多光谱遥感图像和第二多光谱遥感图像进行融合,并对融合后的图像进行滤波、分段拉伸以及锐化,实现了对多光谱遥感图像色彩的还原,使其更加接近地物的真实颜色以及更加接近人眼观察到的地物的颜色。
本发明还提供了一种多光谱遥感图像真实色彩还原***,如图5所示,所述***包括:
多光谱遥感图像获取模块501,用于获取多光谱遥感图像;
第一拉伸模块502,用于采用线性拉伸法将所述多光谱遥感图像由16byte图像转化为8byte图像,得到第一多光谱遥感图像;
像素调整模块503,用于根据公式Gafter=Gbefore·0.8+NIR·0.2对所述多光谱遥感图像中各像素点绿色波段的像素值进行调整,其中,Gafter为所述像素点调整后的绿光波段的像素值,Gbefore为所述像素点调整前的绿光波段的像素值,NIR为所述像素点的红外波段的像素值;
第二拉伸模块504,用于采用线性拉伸法将调整后的所述多光谱遥感图像由16byte图像转化为8byte图像,得到第二多光谱遥感图像;
归一化植被指数计算模块505,用于计算所述多光谱遥感图像的各像素点的归一化植被指数;
归一化植被指数判断模块506,用于判断所述多光谱遥感图像中像素点的归一化植被指数是否大于零;
像素值第一确定模块507,用于当所述多光谱遥感图像中像素点的归一化植被指数大于零时,将所述像素点的像素值确定为所述第二多光谱遥感图像中所述像素点对应的像素值;
像素值第二确定模块508,用于当所述多光谱遥感图像中像素点的归一化植被指数小于等于零时,将所述像素点的像素值确定为所述第一多光谱遥感图像中所述像素点对应的像素值。
优选的,所述***还包括:
大气校正模块,用于采用快速大气校正算法对所述多光谱遥感图像进行大气校正。
高斯滤波模块,用于对所述多光谱遥感图像进行高斯滤波,得到滤波后的所述多光谱遥感图像;
分段拉伸模块,用于对滤波后的所述多光谱遥感图像进行分段拉伸。
所述分段拉伸模块,具体包括:
第一判断单元,用于判断滤波后的所述多光谱遥感图像像素点的像素值是否大于等于0且小于等于90,得到第一判断结果;
第一拉伸单元,用于当所述第一判断结果表示是时,根据公式y=x·(90-30)/90+30对所述像素点的像素值调整;
第二判断单元,用于判断滤波后的所述多光谱遥感图像像素点的像素值是否大于90且小于180,得到第二判断结果;
第二拉伸单元,用于当所述第二判断结果表示是时,保持所述像素点的像素值不变;
第三判断单元,用于判断滤波后的所述多光谱遥感图像像素点的像素值是否大于等于180且小于等于255,得到第三判断结果;
第三拉伸单元,用于当所述第三判断结果表示是时,根据公式y=(x-180)·(245-180)/(255-180)+180对所述像素点的像素值调整;
其中,x为所述像素点拉伸前的像素值,y为所述像素点拉伸后的像素值。
所述***还包括:
图像锐化模块,用于采用USM锐化算法对拉伸后的所述多光谱遥感图像的图像纹理和边缘清晰度进行锐化。
本发明提供的多光谱遥感图像真实色彩还原***,首先将多光谱遥感图像进行拉伸,得到第一多光谱遥感图像,然后对多光谱遥感图像中各像素点绿色波段的像素值进行调整,并将调整后的多光谱遥感图像进行拉伸,得到第二多光谱遥感图像,之后根据各像素点的归一化植被指数对第一多光谱遥感图像和第二多光谱遥感图像进行融合,并对融合后的图像进行滤波、分段拉伸以及锐化,实现了对多光谱遥感图像色彩的还原,使其更加接近地物的真实颜色以及更加接近人眼观察到的地物的颜色。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的***而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种多光谱遥感图像真实色彩还原方法,其特征在于,所述方法包括:
获取多光谱遥感图像;
采用线性拉伸法将所述多光谱遥感图像由16byte图像转化为8byte图像,得到第一多光谱遥感图像;
根据公式Gafter=Gbefore·0.8+NIR·0.2对所述多光谱遥感图像中各像素点绿色波段的像素值进行调整,其中,Gafter为所述像素点调整后的绿光波段的像素值,Gbefore为所述像素点调整前的绿光波段的像素值,NIR为所述像素点的红外波段的像素值;
采用线性拉伸法将调整后的所述多光谱遥感图像由16byte图像转化为8byte图像,得到第二多光谱遥感图像;
计算所述多光谱遥感图像的各像素点的归一化植被指数;
判断所述多光谱遥感图像中像素点的归一化植被指数是否大于零;
如果是,则将所述像素点的像素值确定为所述第二多光谱遥感图像中所述像素点对应的像素值;
如果否,则将所述像素点的像素值确定为所述第一多光谱遥感图像中所述像素点对应的像素值。
2.根据权利要求1所述的多光谱遥感图像真实色彩还原方法,其特征在于,在所述获取多光谱遥感图像之后,所述采用线性拉伸法将所述多光谱遥感图像由16byte图像转化为8byte图像,得到第一多光谱遥感图像之前,还包括:
采用快速大气校正算法对所述多光谱遥感图像进行大气校正。
3.根据权利要求1所述的多光谱遥感图像真实色彩还原方法,其特征在于,在确定所述像素点的像素值之后,还包括:
对所述多光谱遥感图像进行高斯滤波,得到滤波后的所述多光谱遥感图像。
4.根据权利要求3所述的多光谱遥感图像真实色彩还原方法,其特征在于,在所述对所述多光谱遥感图像进行高斯滤波之后,还包括:
对滤波后的所述多光谱遥感图像进行分段拉伸。
5.根据权利要求4所述的多光谱遥感图像真实色彩还原方法,其特征在于,所述对滤波后的所述多光谱遥感图像进行分段拉伸,具体包括:
判断滤波后的所述多光谱遥感图像像素点的像素值是否大于等于0且小于等于90,得到第一判断结果;
当所述第一判断结果表示是时,则根据公式y=x·(90-30)/90+30对所述像素点的像素值调整;
判断滤波后的所述多光谱遥感图像像素点的像素值是否大于90且小于180,得到第二判断结果;
当所述第二判断结果表示是时,则保持所述像素点的像素值不变;
判断滤波后的所述多光谱遥感图像像素点的像素值是否大于等于180且小于等于255,得到第三判断结果;
当所述第三判断结果表示是时,则根据公式y=(x-180)·(245-180)/(255-180)+180对所述像素点的像素值调整;
其中,x为所述像素点拉伸前的像素值,y为所述像素点拉伸后的像素值。
6.根据权利要求4或5所述的多光谱遥感图像真实色彩还原方法,其特征在于,所述方法还包括:
采用USM锐化算法对拉伸后的所述多光谱遥感图像的图像纹理和边缘清晰度进行锐化。
7.一种多光谱遥感图像真实色彩还原***,其特征在于,所述***包括:
多光谱遥感图像获取模块,用于获取多光谱遥感图像;
第一拉伸模块,用于采用线性拉伸法将所述多光谱遥感图像由16byte图像转化为8byte图像,得到第一多光谱遥感图像;
像素调整模块,用于根据公式Gafter=Gbefore·0.8+NIR·0.2对所述多光谱遥感图像中各像素点绿色波段的像素值进行调整,其中,Gafter为所述像素点调整后的绿光波段的像素值,Gbefore为所述像素点调整前的绿光波段的像素值,NIR为所述像素点的红外波段的像素值;
第二拉伸模块,用于采用线性拉伸法将调整后的所述多光谱遥感图像由16byte图像转化为8byte图像,得到第二多光谱遥感图像;
归一化植被指数计算模块,用于计算所述多光谱遥感图像的各像素点的归一化植被指数;
归一化植被指数判断模块,用于判断所述多光谱遥感图像中像素点的归一化植被指数是否大于零;
像素值第一确定模块,用于当所述多光谱遥感图像中像素点的归一化植被指数大于零时,将所述像素点的像素值确定为所述第二多光谱遥感图像中所述像素点对应的像素值;
像素值第二确定模块,用于当所述多光谱遥感图像中像素点的归一化植被指数小于等于零时,将所述像素点的像素值确定为所述第一多光谱遥感图像中所述像素点对应的像素值。
8.根据权利要求7所述的多光谱遥感图像真实色彩还原***,其特征在于,所述***还包括:
大气校正模块,用于采用快速大气校正算法对所述多光谱遥感图像进行大气校正。
高斯滤波模块,用于对所述多光谱遥感图像进行高斯滤波,得到滤波后的所述多光谱遥感图像;
分段拉伸模块,用于对滤波后的所述多光谱遥感图像进行分段拉伸。
9.根据权利要求8所述的多光谱遥感图像真实色彩还原***,其特征在于,所述分段拉伸模块,具体包括:
第一判断单元,用于判断滤波后的所述多光谱遥感图像像素点的像素值是否大于等于0且小于等于90,得到第一判断结果;
第一拉伸单元,当所述第一判断结果表示是时,根据公式y=x·(90-30)/90+30对所述像素点的像素值调整;
第二判断单元,用于判断滤波后的所述多光谱遥感图像像素点的像素值是否大于90且小于180,得到第二判断结果;
第二拉伸单元,用于当所述第二判断结果表示是时,保持所述像素点的像素值不变;
第三判断单元,用于判断滤波后的所述多光谱遥感图像像素点的像素值是否大于等于180且小于等于255,得到第三判断结果;
第三拉伸单元,用于当所述第三判断结果表示是时,根据公式y=(x-180)·(245-180)/(255-180)+180对所述像素点的像素值调整;
其中,x为所述像素点拉伸前的像素值,y为所述像素点拉伸后的像素值。
10.根据权利要求8或9所述的多光谱遥感图像真实色彩还原***,其特征在于,所述***还包括:
图像锐化模块,用于采用USM锐化算法对拉伸后的所述多光谱遥感图像的图像纹理和边缘清晰度进行锐化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710781623.0A CN107507151B (zh) | 2017-09-02 | 2017-09-02 | 一种多光谱遥感图像真实色彩还原方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710781623.0A CN107507151B (zh) | 2017-09-02 | 2017-09-02 | 一种多光谱遥感图像真实色彩还原方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107507151A true CN107507151A (zh) | 2017-12-22 |
CN107507151B CN107507151B (zh) | 2020-09-15 |
Family
ID=60694983
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710781623.0A Active CN107507151B (zh) | 2017-09-02 | 2017-09-02 | 一种多光谱遥感图像真实色彩还原方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107507151B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109472751A (zh) * | 2018-10-29 | 2019-03-15 | 中国资源卫星应用中心 | 一种增强对比度的遥感卫星图像处理方法 |
CN109508397A (zh) * | 2018-11-16 | 2019-03-22 | 中国铁路设计集团有限公司 | 一种基于样图库的遥感影像处理方法 |
CN113222054A (zh) * | 2021-05-28 | 2021-08-06 | 珠江水利委员会珠江水利科学研究院 | 基于特征比值指数的遥感影像融合方法、***、设备及介质 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1770201A (zh) * | 2004-11-05 | 2006-05-10 | 北京师范大学 | 基于小波变换的可调节遥感影像融合方法 |
US20060126959A1 (en) * | 2004-12-13 | 2006-06-15 | Digitalglobe, Inc. | Method and apparatus for enhancing a digital image |
CN101246594A (zh) * | 2008-02-22 | 2008-08-20 | 华南师范大学 | 一种基于梯度场的优化融合遥感图像处理方法 |
CN103198463A (zh) * | 2013-04-07 | 2013-07-10 | 北京航空航天大学 | 基于整体结构和空间细节信息融合的光谱图像全色锐化方法 |
CN103268598A (zh) * | 2013-06-13 | 2013-08-28 | 武汉大学 | 基于Retinex理论的低照度低空遥感影像增强方法 |
CN103558190A (zh) * | 2013-10-22 | 2014-02-05 | 李云梅 | 基于绿光波段的内陆浑浊水体多光谱数据大气校正方法 |
CN103942769A (zh) * | 2013-12-10 | 2014-07-23 | 珠江水利委员会珠江水利科学研究院 | 一种卫星遥感影像融合方法 |
CN105354865A (zh) * | 2015-10-27 | 2016-02-24 | 武汉大学 | 多光谱遥感卫星影像自动云检测方法及*** |
-
2017
- 2017-09-02 CN CN201710781623.0A patent/CN107507151B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1770201A (zh) * | 2004-11-05 | 2006-05-10 | 北京师范大学 | 基于小波变换的可调节遥感影像融合方法 |
US20060126959A1 (en) * | 2004-12-13 | 2006-06-15 | Digitalglobe, Inc. | Method and apparatus for enhancing a digital image |
CN101246594A (zh) * | 2008-02-22 | 2008-08-20 | 华南师范大学 | 一种基于梯度场的优化融合遥感图像处理方法 |
CN103198463A (zh) * | 2013-04-07 | 2013-07-10 | 北京航空航天大学 | 基于整体结构和空间细节信息融合的光谱图像全色锐化方法 |
CN103268598A (zh) * | 2013-06-13 | 2013-08-28 | 武汉大学 | 基于Retinex理论的低照度低空遥感影像增强方法 |
CN103558190A (zh) * | 2013-10-22 | 2014-02-05 | 李云梅 | 基于绿光波段的内陆浑浊水体多光谱数据大气校正方法 |
CN103942769A (zh) * | 2013-12-10 | 2014-07-23 | 珠江水利委员会珠江水利科学研究院 | 一种卫星遥感影像融合方法 |
CN105354865A (zh) * | 2015-10-27 | 2016-02-24 | 武汉大学 | 多光谱遥感卫星影像自动云检测方法及*** |
Non-Patent Citations (3)
Title |
---|
AYMAN M. T. AHMED: "《Color Restoration Techniques for Faded Colors of Old Photos, Printings and Paintings》", 《2009 IEEE INTERNATIONAL CONFERENCE ON ELECTRO/INFORMATION TECHNOLOGY》 * |
王书民等: "《一种基于全局和局部光照估计的Retinex图像增强算法》", 《科学技术与工程》 * |
郝贵青等: "《一种基于RGB颜色空间的色彩还原方法》", 《国外电子测量技术》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109472751A (zh) * | 2018-10-29 | 2019-03-15 | 中国资源卫星应用中心 | 一种增强对比度的遥感卫星图像处理方法 |
CN109472751B (zh) * | 2018-10-29 | 2021-04-02 | 中国资源卫星应用中心 | 一种增强对比度的遥感卫星图像处理方法 |
CN109508397A (zh) * | 2018-11-16 | 2019-03-22 | 中国铁路设计集团有限公司 | 一种基于样图库的遥感影像处理方法 |
CN113222054A (zh) * | 2021-05-28 | 2021-08-06 | 珠江水利委员会珠江水利科学研究院 | 基于特征比值指数的遥感影像融合方法、***、设备及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN107507151B (zh) | 2020-09-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104851113B (zh) | 多分辨率遥感影像的城市植被自动提取方法 | |
CN102682436B (zh) | 一种改进的多尺度Retinex理论的图像增强方法 | |
CN105354865B (zh) | 多光谱遥感卫星影像自动云检测方法及*** | |
CN109754377A (zh) | 一种多曝光图像融合方法 | |
CN108596103A (zh) | 基于最佳光谱指数选择的高分辨率卫星遥感影像建筑物提取方法 | |
CN103942769B (zh) | 一种卫星遥感影像融合方法 | |
CN107680054A (zh) | 雾霾环境下多源图像融合方法 | |
CN104537634B (zh) | 动态图像中去除雨滴影响的方法和*** | |
CN107507151A (zh) | 一种多光谱遥感图像真实色彩还原方法及*** | |
CN107358585A (zh) | 基于分数阶微分及暗原色先验的雾天图像增强方法 | |
CN105046701B (zh) | 一种基于构图线的多尺度显著目标检测方法 | |
CN104182949A (zh) | 基于直方图特征点配准的影像匀色与融合方法及*** | |
CN106683056A (zh) | 一种机载光电红外数字图像处理方法及装置 | |
CN102789635A (zh) | 一种图像增强方法和装置 | |
CN109544487A (zh) | 一种基于卷积神经网络的红外图像增强方法 | |
CN107680114A (zh) | 一种基于计算机图像处理的草地盖度测算方法 | |
CN108257094A (zh) | 基于暗通道的遥感影像薄雾快速去除方法 | |
CN108288258A (zh) | 一种针对恶劣天气条件下的低质图像增强方法 | |
WO2018166083A1 (zh) | 一种对比度增强和视觉显著度优化方法在高尔夫球场图中的应用 | |
CN109815894A (zh) | 一种针对哨兵2a影像的建筑物阴影提取处理方法 | |
CN109919859A (zh) | 一种户外场景图像去雾增强方法、计算设备及其存储介质 | |
CN104637036A (zh) | 一种中国古画增强方法 | |
CN107146264A (zh) | 一种提取材质表面几何和光照物理属性的方法 | |
CN102063700A (zh) | 卫星遥感图像生成方法及其*** | |
Wang et al. | An adaptive correction algorithm for non-uniform illumination panoramic images based on the improved bilateral gamma function |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |