CN102169575B - 一种非基色单色光图像数据的处理方法及其处理装置 - Google Patents
一种非基色单色光图像数据的处理方法及其处理装置 Download PDFInfo
- Publication number
- CN102169575B CN102169575B CN 201110101912 CN201110101912A CN102169575B CN 102169575 B CN102169575 B CN 102169575B CN 201110101912 CN201110101912 CN 201110101912 CN 201110101912 A CN201110101912 A CN 201110101912A CN 102169575 B CN102169575 B CN 102169575B
- Authority
- CN
- China
- Prior art keywords
- image
- view data
- data
- value
- primary colours
- 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.)
- Expired - Fee Related
Links
- 238000012545 processing Methods 0.000 title claims abstract description 20
- 238000000034 method Methods 0.000 title abstract description 7
- 239000003086 colorant Substances 0.000 claims description 24
- 238000012937 correction Methods 0.000 claims description 5
- 238000004040 coloring Methods 0.000 claims description 4
- LTXREWYXXSTFRX-QGZVFWFLSA-N Linagliptin Chemical compound N=1C=2N(C)C(=O)N(CC=3N=C4C=CC=CC4=C(C)N=3)C(=O)C=2N(CC#CC)C=1N1CCC[C@@H](N)C1 LTXREWYXXSTFRX-QGZVFWFLSA-N 0.000 claims 2
- 241000233855 Orchidaceae Species 0.000 description 7
- 230000004297 night vision Effects 0.000 description 7
- 230000005540 biological transmission Effects 0.000 description 3
- 230000007123 defense Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 102220032806 rs367543154 Human genes 0.000 description 3
- 102220074140 rs45479698 Human genes 0.000 description 3
- 102220109513 rs760091696 Human genes 0.000 description 3
- 102220065847 rs794726907 Human genes 0.000 description 3
- 239000004744 fabric Substances 0.000 description 2
- 238000002594 fluoroscopy Methods 0.000 description 2
- 238000005286 illumination Methods 0.000 description 2
- 210000000936 intestine Anatomy 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 210000002784 stomach Anatomy 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000001174 ascending effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000004438 eyesight Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 230000008521 reorganization Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
Images
Landscapes
- Image Processing (AREA)
Abstract
一种非基色单色光图像数据的处理方法及其处理装置,利用同一像素的色光分量的不同,运用数据处理单元获取高动态范围图像数据,降低图像的噪声,大大提高图像质量;使原来图像中由于超过该色分量传感器信号最大值而过曝光造成数据丢失的数据,可以从欠曝光数据中恢复;而原来由于欠曝光分量中暗的部分的图像数据的噪声大,可以从过曝光相同像素的低噪声图像数据中恢复,提高图像质量;其它部分的图像数据也可以用过曝光,正常曝光,欠曝光的数据平均获得,降低图像的噪声,增加图像的动态范围,提高图像质量。
Description
技术领域:本发明涉及图像处理领域,特别涉及非基色单色光图像数据的处理方法及其处理装置。
背景技术:目前非基色单色光大量运用于医疗设备、物体的无损伤性探查,夜视等领域,无论是医疗设备、物体的无损伤性探查,夜视等均是采集荧光图像,并将其转换为电信号获得数码图像。但由于光电转换设备的制约,图像的灰阶总是不能超过设备的采集阈值,所以,必然使设备的采集阈值确定了图像与物的差距。
图像的灰阶是有局限的,怎样获得灰阶高动态范围图像,是使图像更接近被摄物,也是提高图像质量的重要课题,从中国专利,或其它国家专利和见到的记载中,有通过多次曝光进行处理后合成一幅高动态图像,但存在不同时刻的曝光可能存在的差异的重要缺点,也有中国专利公开号为CN101917551A,名称为:一种单次曝光的高动态范围图像获取方法,运用电荷耦合器件输出的模拟电信号通过串行信号线,分别输入高增益模拟数字转换器和低增益模拟信号转换器进行采样,由高增益模拟信号转换器通过数字并行总线输出一幅亮图,由低增益信号转换器通过数字并行总线输出一幅暗图,通过处理后合成一幅高动态图像。但存在电荷耦合元件在采集信号时受其电荷耦合元件本身的动态范围限制。
发明内容:为了克服或缩小物与图像的差距,本发明的目的在于提出一种利用非基色单色光图像数据的同一像素的色光分量的不同,运用数据处理单元获取高动态范围图像数据的方法及其处理装置。
本发明的技术方案是:一种非基色单色光图像数据的处理方法,其特征在于,至少包括下列方法对图像数据进行处理:
A.通过二色及其以上的处理设备,采集到二色以上的图像数据F12…i(x,y,z),其中,x为图像的水平坐标,y为图像的垂直坐标,z=1,2,……,i(i≥2)
分离各色分量为:
F1(x,y)=F12…i(x,y,1)
F2(x,y)=F12…i(x,y,2)
……
Fi(x,y)=F12…i(x,y,i)
……
其中,M为图像的水平最大坐标,N为图像的垂直最大坐标;
B.将F1(x,y),F2(x,y),……Fi(x,y)传输到数据处理单元,并进行:
B.1选择F1(x,y),F2(x,y),……Fi(x,y)中的n个相同像素点(n为1至M×N之间),求出F1(x,y),F2(x,y),……Fi(x,y)的n个相同像素点和:
……
选择S1,S2,……Si的其中之一为Sj(j为1至i);
B.2利用计算单元计算,每个S1,S2,……Si与Sj比例的关系:
……
……
在K1j,K2j,……,Kij中,选择Kmj为同一性校正的标准比例,m为1至i;
B.3利用计算单元,对F1(x,y),F2(x,y),……,Fi(x,y)的各像素值进行同一性校正,
……
Dj(x,y)=Fj(x,y)×Kjj/Kmj=Fj(x,y)/Kmj
……
……
B.4利用计算单元从D1(x,y),D2(x,y),……Di(x,y)中再造F(x,y),
其中,F(x,y)是D1(x,y),D2(x,y),……Di(x,y)同一像素的最小值至最大值之间的数值。
所述的选择F1(x,y),F2(x,y),……Fi(x,y)中的n个相同像素点时,处理设备的采集最高值即最大信号幅度数值L,设定系数a,a为0.7-1之间,设定系数b,b为0-0.3之间,设置为最高阈值TH=a×L和最低阈值TL=b×L,对图像中每一个象素点,去除所有分量亮度值超过TH的象素点和所有分量亮度值低于TL的象素点,将剩下合格的象素点的1,2,……,i分量分别求和,得出S1,S2,……,Si。
所述的一种非基色单色光图像数据的处理方法,其特征在于:所述的a为0.97-1之间,所述的b为0-0.1之间。
所述的一种非基色单色光图像数据的处理方法,其特征在于:所述的S1,S2,……Si选最大为Sj。
所述的一种非基色单色光图像数据的处理方法,其特征在于:所述在K1j,K2j,……,Kij中,选最大为Kmj,所述的选择Kmj为同一性校正的标准比例。
所述的一种非基色单色光图像数据的处理方法,其特征在于:所述的B4中的再造F(x,y)对D1(x,y),D2(x,y),……,Di(x,y)中选取,遵循S1,S2,…..,Si中数值大的所对应的D1(x,y),D2(x,y),……,Di(x,y)对F(x,y)中的小数值区域贡献大,S1,S2,…..,Si中数值小的所对应的D1(x,y),D2(x,y),……,Di(x,y)对F(x,y)中的大数值区域贡献大的原则。
所述的利用计算单元从D1(x,y),D2(x,y),……Di(x,y)中再造F(x,y),
设定c为0.7-1之间,d为0-0.6之间,L为最大信号幅度数值,设置为最高阈值TH1=c×L和最低阈值TL1=d×L,
按照S1,S2,……Si由小到大的顺序排列D1(x,y),D2(x,y),……,Di(x,y),得到D1S(x,y),D2S(x,y),……DiS(x,y),其对应的K1S,K2S,……KiS是按相反顺序排列,即由大到小的顺序排列,
对图像中每一个象素点,进行如下处理,得到再造新图像,
如果最小D1S(x,y)≤TL1×K1S/Kmj
F(x,y)=DiS(x,y),即用最大;
如果最小D1S(x,y)>TL1×K1S/Kmj和DiS(x,y)≤TH1×KiS/Kmj
F(x,y)=(D1S(x,y)+D2S(x,y)+……+DiS(x,y))/i;
如果DiS(x,y)>TH1×KiS/Kmj和Di-1S(x,y)≤TH1×Ki-1S/Kmj
F(x,y)=(D1S(x,y)+D2S(x,y)+……+Di-1S(x,y))/(i-1);
如果Di-1S(x,y)>TH1×Ki-1S/Kmj和Di-2S(x,y)≤TH1×Ki-2S/Kmj
F(x,y)=(D1S(x,y)+D2S(x,y)+……+Di-2S(x,y))/(i-2);
……
如果D3S(x,y)>TH1×K3S/Kmj和D2S(x,y)≤TH1×K2S/Kmj
F(x,y)=(D1S(x,y)+D2S(x,y))/2;
如果D2S(x,y)>TH1×K2S/Kmj
F(x,y)=D1S(x,y),即用最小。
一种非基色单色光图像数据的处理装置,它由一个非基色单色光的被摄物和被摄物图像数据的处理设备,进行运算和处理数据的单元组成,其特征在于:被摄物图像数据的处理设备采集二色以上的图像数据F12…i(x,y,z),其中,x为图像的水平坐标,y为图像的垂直坐标,z=1,2,……,i(i≥2),并将F12…i(x,y,z)图像数据分离为F1(x、y)、F2(x、y)……Fi(x、y),i≥2,并将F1(x、y)、F2(x、y)……Fi(x、y)传输到运算和处理数据单元。
所述的一种非基色单色光图像数据的处理装置,其特征在于所述被摄物图像数据的处理设备和运算和处理数据单元是至少包括能够对下列图像数据进行处理的单元:
A.通过二色及其以上的处理设备,采集到二色以上的图像数据,F12…i(x,y,z),其中,
z=1,2,……,i(i≥2)
F1(x,y)=F12…i(x,y,1)
F2(x,y)=F12…i(x,y,2)
……
Fi(x,y)=F12…i(x,y,i)
……
B.将F1(x,y),F2(x,y),……Fi(x,y)传输到数据处理单元,并进行:
B.1选择F1(x,y),F2(x,y),……Fi(x,y)中的n个相同像素点,n为1至M×N之间,求出F1(x,y),F2(x,y),……Fi(x,y)之间的比例关系:
……
选择S1,S2,……Si的其中之一为Sj,
B.2计算每个S1,S2,……Si与Sj比例的关系:
……
在K1j,K2j,……,Kij中,选择Kmj为同一性校正的标准比例,式中m为1至i,
B.3对F1(x,y),F2(x,y),…….Fi(x,y)进行各像素值进行同一性校正。
……
Dj(x,y)=Fj(x,y)/Kmj
……
……
B.4从D1(x,y),D2(x,y),……Di(x,y)中再造F(x,y),
其中,F(x,y)是D1(x,y),D2(x,y),……Di(x,y)同一像素的最大值或最小值之间的数值的。
根据国际照明委员会(CIE)于1931年创立CIE 1931色彩空间,即CIE 1931 XYZ色彩空间可知,一种非基色单色光均可分解为红R,绿G,兰B三色分量,为此,通过R,G,B三色采集设备,必可以采集到F1(x,y)、F2(x,y)、F3(x,y)三个数据,除极端情况下,F1(x,y)≠F2(x,y)≠F3(x,y),有大小之分,且相互之间成比例关系。映射到图像上,即有图像亮暗之分。为此,在适当曝光的情况下,对于F1(x,y)、F2(x,y)、F3(x,y)必有“过曝光,适中曝光,欠曝光”之分。正是利用“过曝光”对应的数据中,对较暗像素的信息记载的较为充分;“欠曝光”对应的数据中,对较亮像素的信息记载的较为充分的原理,对一次曝光同时采集到的F1(x,y)、F2(x,y)、F3(x,y)进行处理后,合成一幅含有高动态图像信息的数据F(x,y),提高图像的信噪比。
在本发明中运用数字图像理论,一个图像对应一个矩阵数据F12…i(x,y,z),其中,x为图像的水平坐标,y为图像的垂直坐标,z=1,2,……,i(i≥2)的二色以上图像数据,将F12…i(x,y,z)图像数据分离为F1(x、y)、F2(x、y)……Fi(x、y),i≥2,并将F1(x、y)、F2(x、y)……Fi(x、y)传输到运算和处理数据单元进行运算和处理。
根据色彩理论,任何一个非基色单色光是由二色以上的基色光组成。所以一个非基色单色光图像对应一个数据矩阵F12…i(x,y,z),可将F12…i(x,y,z)图像数据分离为F1(x、y)、F2(x、y)……Fi(x、y),i≥2。由于在一个F12…i(x,y,z)图像数据中,分离的F1(x、y)、F2(x、y)……Fi(x、y)同一像素对应的F1(x、y)、F2(x、y)……Fi(x、y)数据是有一定比例的,并且该比例对于该图像同一像素对应的F1(x、y)、F2(x、y)……Fi(x、y)是相同的。为此可利用每个S1,S2,……Si与Sj比例进行各像素值进行同一性校正。从理论上而言,被掩盖的数据可以用F1(x、y)、F2(x、y)……Fi(x、y)数据代替,重组F(x,y)。
本发明的优点在于:利用同一像素的色光分量的不同,运用数据处理单元获取高动态范围图像数据,降低图像的噪声,大大提高图像质量。使原来图像中由于超过该色分量传感器信号最大值而过曝光造成数据丢失的数据,可以从欠曝光数据中恢复;而原来由于欠曝光分量中暗的部分的图像数据的噪声大,可以从过曝光相同像素的低噪声图像数据中恢复,提高图像质量;其它部分的图像数据也可以用过曝光,正常曝光,欠曝光的数据平均获得,降低图像的噪声,增加图像的动态范围,提高图像质量。
附图说明:
图1-1为本发明采用专用灰度采集设备采集的灰度图像。
图1-2为本发明采用彩色采集设备采集的彩色图像。
图1-3为本发明1(红)分量图像。
图1-4为本发明2(绿)分量图像。
图1-5为本发明3(兰)分量图像。
图1-6为本发明再造图像。
图1-7为本发明原始图像的A部分图像。
图1-8为本发明再造图像A相同部分的图像。
图1-9为本发明原始图像的B部分图像。
图1-10为本发明再造图像B相同部分的图像。
图2-1为本发明采用专用灰度采集设备采集灰度图像。
图2-2为本发明采用彩色采集设备采集彩色图像。
图2-3为本发明1(红)分量图像。
图2-4为本发明2(绿)分量图像。
图2-5为本发明3(兰)分量图像。
图2-6为本发明再造的新图像。
图2-7为本发明原始图像的A部分图像。
图2-8为本发明再造图像A相同部分的图像。
图2-9为本发明原始图像的B部分图像。
图2-10为本发明再造图像B相同部分的图像。
图3-1为本发明夜视图像。
图3-2为本发明1(红)分量图像。
图3-3为本发明2(绿)分量图像。
图3-4为本发明3(兰)分量图像。
图3-5为本发明灰色图像。
图3-6为本发明再造图像。
图3-7为本发明原始图像的A部分图像。
图3-8为本发明再造图像A相同部分的图像。
具体实施方式:本发明的技术方案通过以下实施例作详细说明,但不限于下列实施例。
例1、X光灰阶分辨率及清晰度分辨率卡的照片图像的获取及处理
一般情况下,X光信息是被照射物体的在照射方向上密度信息的叠加形成的,因此,只有强度信息,并无色彩信息,一般采用专用灰度采集设备采集其强度信号,见图1-1。
本发明是用医用X光胃肠机对灰阶分辨率测试卡及清晰度分辨率测试卡进行X光拍摄的彩色图像采集装置,对该X光设备的影像增强器输出荧光屏采集到的彩色照片图像数据,见图1-2。上述图像数据的采集照相装置是佳能EOS450D照相机,将该照相装置获取的三色图像数据传输到处理运算单元,本例中用计算机代替。
1.计算机读入包含红R、绿G、兰B图像分量的彩色uint16(16位整数)图像。图像R分量为1分量,图像G分量为2分量,图像B分量为3分量。
F123(x,y,z)为彩色图像,x为图像水平坐标,y为图像垂直坐标,z为1,2,3图像分量坐标。z=1为1分量图像,z=2为2分量图像,z=3为3分量图像。
F1(x,y)为1分量图像,x为图像水平坐标,y为图像垂直坐标。
F2(x,y)为2分量图像,x为图像水平坐标,y为图像垂直坐标。
F3(x,y)为3分量图像,x为图像水平坐标,y为图像垂直坐标。
F1(x,y)=F123(x,y,1)
F2(x,y)=F123(x,y,2)
F3(x,y)=F123(x,y,3)
2.观察可以发现,1,2,3图像分量图像中,2分量最亮,1分量次之,3分量最暗。
1分量图像如图1-3。
F1(x,y)=F123(x,y,1)
M=2480,N=2480
2分量图像如图1-4
F2(x,y)=F123(x,y,2)
M=2480,N=2480
3分量图像如图1-5
F3(x,y)=F123(x,y,3)
M=2480,N=2480
3.处理运算单元对其进行如下程序的处理:
3.1.设置最高阈值TH和最低阈值TL
比如,16为无符号整数的最大为65535,
TH=65535*0.987=64684
TL=65535*0.037=2425
3.2.对图像中每一个象素点,去除所有分量亮度值超过TH的象素点和所有分量低于TL的象素点,将剩下合格的象素点的1,2,3分量分别求和,得出S1,S2,S3。
实际图像的水平像素数为2480,图像的垂直像素数为2480。
求出S1=3.7641e+010,S2=6.7229e+010,
S3=2.6009e+010,
选择最大的S2为Sj,j=2,S2=6.7229e+010,
3.3.计算分量2与分量1,3间的比例关系:
K12=S2/S1=1.7861,
K32=S2/S3=2.5848,
找出最大比例数为K32,确定为同一性校正的标准比例。
3.4.为了减少数据精度对图像质量的影响,将Uint16整数彩色图像数据转换为Double双精度数据彩色图像。
3.5.对F1(x,y),F3(x,y)各像素值进行同一性校正。
D2(x,y)=F2(x,y)/K32;
3.6.从D1(x,y),D2(x,y),D3(x,y)中再造F(x,y)。
3.6.1.计算新的阈值:比如:双精度最高亮度值为1,T31=1*0.987,T41=1*0.037
3.6.2.对于每一个象素点,进行如下处理:
再造的新图像如图1-6。
3.7.结果比较:
3.7.1.动态范围比较:
将原始双精度图像数据RGB彩色图像转换为灰色图像,图像的最大亮度为MAX=0.3869,最小亮度为MIN=6.3573e-004,MAX-MIN=0.3862,即原始图像的动态范围是0.3862。
我们再造的双精度图像的最大亮度为MAX=1,最小亮度为MIN=8.0285e-004,MAX-MIN=0.9992,即图像的动态范围是0.9992,显然,动态范围是原来动态范围的0.9992/0.3862=2.5872倍,有了明显扩展。
3.7.2.图像效果比较
图1-7是处理前原始图像的部分图像,图1-8是再造图像相同部分的图像;图1-9是处理前原始图像的部分图像,图1-10是再造图像相同部分的图像。从图中可以明显看出,由于再造图像的动态范围比原始的动态范围大,原始图像中由于过曝光丢失的信息可以在再造图像中清晰的显示出来。
例2、X光胸部透视视频图像的获取及处理
采用专用灰度采集设备采集其强度图像。见图2-1。
采用医用X光胃肠机对人体胸部透视图像进行X光摄像的彩色图像采集装置,对该X光设备的影像增强器输出荧光屏用录像方式采集到的一组彩色图像数据,选取其中一张图像,去除无用部分数据,再将其进行处理,见图2-2。其它图像按相同方法处理,就可以得到一组高动态图像,将其连续播放,就可看到高动态范围的视频影像。上述图像数据的采集装置是索尼NEX-5数码照相机,将该照相装置获取的三色图像数据传输到处理运算单元,本例中用计算机代替。
1.计算机读入包含R(红)G(绿)B(兰)图像分量的彩色uint16(16位整数)图像。图像R(红)分量为1分量,图像G(绿)分量为2分量,图像B(兰)分量为3分量。
F123(x,y,z)为彩色图像,x为图像水平坐标,y为图像垂直坐标,z为1,2,3图像分量坐标。z=1为1分量图像,z=2为2分量图像,z=3为3分量图像,
F1(x,y)为1分量图像,x为图像水平坐标,y为图像垂直坐标。
F2(x,y)为2分量图像,x为图像水平坐标,y为图像垂直坐标。
F3(x,y)为3分量图像,x为图像水平坐标,y为图像垂直坐标。
F1(x,y)=F123(x,y,1)
F2(x,y)=F123(x,y,2)
F3(x,y)=F123(x,y,3)
2.观察可以发现,1,2,3图像分量图像中,2分量最亮,1分量次之,3分量最暗。
1分量图像如图2-3。
F1(x,y)=F123(x,y,1)
M=1080,N=1080
2分量图像如图2-4
F2(x,y)=F123(x,y,2)
M=1080,N=1080
3分量图像如图2-5
F3(x,y)=F123(x,y,3)
M=1080,N=1080
3.处理运算单元对其进行如下程序的处理:
3.1.设置最高阈值TH和最低阈值TL
比如,16为无符号整数的最大为65535,
TH=65535*0.987=64684
TL=65535*0.037=2425
3.2.对图像中每一个象素点,去除最亮2分量亮度值超过TH的象素点和最暗3分量低于TL的象素点,将剩下合格的象素点的1,2,3分量分别求,得出S1,S2,S3。
实际图像的水平像素数为1080,图像的垂直像素数为1080。
求出S1=1.9382e+010,S2=3.5353e+010,
S3=1.3919e+010,
选择最大的S2为Sj,j=2,S2=3.5353e+010,
3.3.计算分量2与分量1,3间的比例关系:
K12=S2/S1=1.8240,
K32=S2/S3=2.5397
找出最大比例数为K32,确定为同一性校正的标准比例。
3.4.为了减少数据精度对图像质量的影响,将Uint16整数彩色图像数据转换为Double(双精度)数据彩色图像。
3.5.对F1(x,y),F3(x,y)各像素值的进行同一性校正。
D2(x,y)=F2(x,y)/K32;
3.6.从D1(x,y),D2(x,y),……Di(x,y)中再造F(x,y)。
3.6.1.计算新的阈值:比如:双精度最高亮度值为1,T31=1*0.987,T41=1*0.037
3.6.2.对于每一个象素点,进行如下处理:
再造的新图像如图2-6。
3.7.结果比较:
3.7.1.动态范围比较:
原始图像的最大亮度为MAX=0.3684,最小亮度为MIN=7.1865e-004,MAX-MIN=0.3676,即原始图像的动态范围是0.3676。
我们再造的双精度图像的最大亮度为MAX=0.6514,最小亮度为MIN=9.0122e-004,MAX-MIN=0.6505,即图像的动态范围是0.6505,显然,动态范围是原来动态范围的0.6505/0.3676=1.7695倍,有了明显扩展。
3.7.2.图像效果比较
图2-7是处理前原始图像的A部分图像,图2-8是再造图像A相同部分的图像;图2-9是处理前原始图像的B部分图像,图2-10是再造图像B相同部分的图像。从图中可以明显看出,由于再造图像的动态范围比原始的动态范围大,原始图像中由于过曝光丢失的信息可以在再造图像中清晰的显示出来。
例3、一张夜视成像仪图像获取及处理
从美国国防部(U.S.DEPARTMET OF DEFENSE)网站(http://www.defense.gov/news/newsarticle.aspx?id=1617)下载的夜视图像hires_Night Vision.jpg来进行测试,图像水平为1222,垂直为1875,数据位数为uint8。夜视图像见图3-1。
将该三色图像数据传输到处理运算单元,本例中用计算机代替。
1.计算机读入包含R(红)G(绿)B(兰)图像分量的彩色uint8(8位整数)图像。
图像R(红)分量为1分量,图像G(绿)分量为2分量,图像B(兰)分量为3分量。
F123(x,y,z)为彩色图像,x为图像水平坐标,y为图像垂直坐标,z为1,2,3图像分量坐标。z=1为1分量图像,z=2为2分量图像,z=3为3分量图像,
F1(x,y)为1分量图像,x为图像水平坐标,y为图像垂直坐标。
F2(x,y)为2分量图像,x为图像水平坐标,y为图像垂直坐标。
F3(x,y)为3分量图像,x为图像水平坐标,y为图像垂直坐标。
F1(x,y)=F123(x,y,1)
F2(x,y)=F123(x,y,2)
F3(x,y)=F123(x,y,3)
2.观察可以发现,1,2,3图像分量图像中,2分量最亮,1分量次之,3分量最暗。
1分量图像如图3-2。
F1(x,y)=F123(x,y,1)
M=1222,N=1875
2分量图像如图3-3
F2(x,y)=F123(x,y,2)
M=1222,N=1875
3分量图像如图3-4
F3(x,y)=F123(x,y,3)
M=1222,N=1875
3.处理运算单元对其进行如下程序的处理:
3.1.设置最高阈值TH和最低阈值TL
比如,8为无符号整数的最大为255,
TH=255*0.987=252
TL=255*0.037=9
3.2.对图像中每一个象素点,去除最亮2分量亮度值超过TH的象素点和最暗3分量低于TL的象素点,将剩下合格的象素点的1,2,3分量分别求,得出S1,S2,S3。
实际图像的水平像素数为1222,图像的垂直像素数为1875。
求出S1=59061635,S2=190827744,S3=49007209,
选择最大的S2为Sj,j=2,S2=190827744,
3.3.计算分量2与分量1,3间的比例关系:
K12=S2/S1=3.2310,
K32=S2/S3=3.8939,
找出最大比例数为K32,确定为同一性校正的标准比例。
3.4.为了减少数据精度对图像质量的影响,将Uint8整数彩色图像数据转换为Double(双精度)数据彩色图像.
3.5.对F1(x,y),F3(x,y)各像素值进行同一性校正。
D2(x,y)=F2(x,y)/K32;
3.6.从D1(x,y),D2(x,y),……Di(x,y)中再造F(x,y)。
3.6.1.计算新的阈值:比如:双精度最高亮度值为1,T31=1*0.987,T41=1*0.037
3.6.2.对于每一个象素点,进行如下处理:
再造的新图像如图3-6。
3.7.结果比较:
3.7.1.动态范围比较:
将原始夜视RGB彩色图像转换为灰色图像,见3-5。图像的最大亮度为MAX=0.2568,最小亮度为MIN=0,MAX-MIN=0.2568,即原始图像的动态范围是0.2568。
我们再造的双精度图像的最大亮度为MAX=1,最小亮度为MIN=0,MAX-MIN=1,即图像的动态范围是1,显然,动态范围是原来动态范围的1/0.2568=3.894倍,有了明显扩展。
3.7.2.图像效果比较
图3-7是处理前原始图像的A部分图像,图3-8是再造图像A相同部分的图像。从图中可以明显看出,由于再造图像的动态范围比原始的动态范围大,原始图像中由于过曝光丢失的信息可以在再造图像中清晰的显示出来,并且噪声也明显减少。
Claims (8)
1.一种非基色单色光图像数据的处理方法,其特征在于,至少包括下列步骤对图像数据进行处理:
A.通过二色及其以上的处理设备,采集到二色以上的图像数据F12…i(x,y,z),其中,x为图像的水平坐标,y为图像的垂直坐标,z=1,2,……,i,i≥2,
分离各色分量为:
F1(x,y)=F12…i(x,y,1)
F2(x,y)=F12…i(x,y,2)
……
Fi(x,y)=F12…i(x,y,i)
……
其中,M为图像的水平最大坐标,N为图像的垂直最大坐标;
B.将F1(x,y), F2(x,y),……Fi(x,y)传输到数据处理单元,并进行:
B.1选择F1(x,y),F2(x,y),……Fi(x,y)中的n个相同像素点,n为l至M×N之间;求出F1(x,y),F2(x,y),……Fi(x,y)的n个相同像素点和:
……
选择S1,S2,……Si的其中之一为Sj,j为1至i;
B.2利用计算单元计算,每个S1,S2,……Si与Sj比例的关系:
……
……
在K1j,K2j,……,Kij中,选择Kmj为同一性校正的标准比例,m为1至i;
B.3利用计算单元,对F1(x,y),F2(x,y),……,Fi(x,y)的各像素值进行同一性校正,
……
Dj(x,y)=Fj(x,y)×Kjj/Kmj=Fj(x,y)/Kmj
……
……
B.4利用计算单元从D1(x,y),D2(x,y),……Di(x,y)中再造F(x,y),
其中,F(x,y)是D1(x,y),D2(x,y),……Di(x,y)同一像素的最小值至最大值之间的数值。
2.根据权利要求1所述的一种非基色单色光图像数据的处理方法,其特征在于:所述的选择F1(x,y),F2(x,y),……Fi(x,y)中的n个相同像素点时,处理设备的采集最高值即最大信号幅度数值L,设定系数a,a为0.7-1之间,设定系数b,b为00.3之间,设置为最高阈值TH=a×L和最低阈值TL=b×L,对图像中每一个像素点,去除所有分量亮度值超过TH的像素点和所有分量亮度值低于TL的像素点,将剩下合格的像素点的1,2,……,i分量分别求和,得出S1,S2,……,Si。
3.根据权利要求2所述的一种非基色单色光图像数据的处理方法,其特征在于:所述的a为0.97-1之间,所述的b为0-0.1之间。
4.根据权利要求1所述的一种非基色单色光图像数据的处理方法,其特征在于:所述的S1,S2,……Si选最大为Sj。
5.根据权利要求1所述的一种非基色单色光图像数据的处理方法,其特征在于:所述在K1j,K2j,……,Kij中,选最大为Kmj,所述的Kmj为同一性校正的标准比例。
6.根据权利要求1所述的一种非基色单色光图像数据的处理方法,其特征在于:所述的B.4中的再造F(x,y)对D1(x,y),D2(x,y),……,Di(x,y)的选取,遵循S1,S2,……,Si中数值大的所对应的D1(x,y),D2(x,y),……,Di(x,y)对F(x,y)中的小数值区域贡献大,S1,S2,……,Si中数值小的所对应的D1(x,y),D2(x,y),……,Di(x,y)对F(x,y)中的大数值区域贡献大的原则。
7.根据权利要求1或6所述的一种非基色单色光图像数据的处理方法,其特征在于:所述的利用计算单元从D1(x,y),D2(x,y),……Di(x,y)中再造F(x,y),
设定c为0.7-1之间,d为0-0.6之间,L为最大信号幅度数值,设置为最高阈值TH1=c×L和最低阈值TL1=d×L,
按照S1,S2,……Si由小到大的顺序排列D1(x,y),D2(x,y),……,Di(x,y),得到D1S(x,y),D2S(x,y),……DiS(x,y),其对应的K1S,K2S,……KiS是按相反顺序排列,即由大到小的顺序排列,
对图像中每一个像素点,进行如下处理,得到再造新图像,
如果最小D1S(x,y)≤TL1×K1S/Kmj
F(x,y)=DiS(x,y),即用最大;
如果最小D1S(x,y)>TL1×K1S/Kmj和DiS(x,y)≤TH1×KiS/Kmj
F(x,y)=(D1S(x,y)+D2S(x,y)+……+DiS(x,y))/i;
如果DiS(x,y)>TH1×KiS/Kmj和Di-1S(x,y)≤TH1×Ki-1S/Kmj
F(x,y)=(D1S(x,y)+D2S(x,y)+……+Di-1S(x,y))/(i-1);
如果Di-1S(x,y)>TH1×Ki-1S/Kmj和Di-2S(x,y)≤TH1×Ki-2S/Kmj
F(x,y)=(D1S(x,y)+D2S(x,y)+……+Di-2S(x,y))/(i-2);
……
如果D3S(x,y)>TH1×K3S/Kmj和D2S(x,y)≤TH1×K2S/Kmj
F(x,y)=(D1S(x,y)+D2S(x,y))/2;
如果D2S(x,y)>TH1×K2S/Kmj
F(x,y)=D1S(x,y),即用最小。
8.一种非基色单色光图像数据的处理装置,它由一个非基色单色光的被摄物和被摄物图像数据的处理设备,进行运算和处理数据的单元组成,其特征在于:被摄物图像数据的处理设备采集二色以上的图像数据F12…i(x,y,z),其中,x为图像的水平坐标,y为图像的垂直坐标,z=1,2,……,i,i≥2,并将F12…i(x,y,z)图像数据分离为F1(x、y)、F2(x、y)……Fi(x、y),i≥2,并将F1(x、y)、F2(x、y)……Fi(x、y)传输到运算和处理数据单元;
所述被摄物图像数据的处理设备及运算和处理数据的单元是至少包括能够对下列图像数据进行处理的单元:
A.通过二色及其以上的处理设备,采集到二色以上的图像数据,F12…i(x,y,z),其中,z=1,2,……,i,i≥2,
F1(x,y)=F12…i(x,y,1)
F2(x,y)=F12…i(x,y,2)
……
Fi(x,y)=F12…i(x,y,i)
……
其中,M为图像的水平最大坐标,N为图像的垂直最大坐标;
B.将F1(x,y),F2(x,y),……Fi(x,y)传输到数据处理单元,并进行:
B.1选择F1(x,y),F2(x,y),……Fi(x,y)中的n个相同像素点,n为1至M×N之间;求出F1(x,y),F2(x,y),……Fi(x,y)的n个相同像素点和:
……
选择S1,S2,……Si的其中之一为Sj,j为1至i;
B.2利用计算单元计算,每个S1,S2,……Si与Sj比例的关系:
……
……
在K1j,K2j,……,Kij中,选择Kmj为同一性校正的标准比例,m为1至i;
B.3利用计算单元,对F1(x,y),F2(x,y),……,Fi(x,y)的各像素值进行同一性校正,
……
Dj(x,y)=Fj(x,y)×Kjj/Kmj=Fj(x,y)/Kmj
……
……
B.4利用计算单元从D1(x,y),D2(x,y),……Di(x,y)中再造F(x,y),
其中,F(x,y)是D1(x,y),D2(x,y),……Di(x,y)同一像素的最小值至最大值之间的数值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110101912 CN102169575B (zh) | 2011-04-22 | 2011-04-22 | 一种非基色单色光图像数据的处理方法及其处理装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110101912 CN102169575B (zh) | 2011-04-22 | 2011-04-22 | 一种非基色单色光图像数据的处理方法及其处理装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102169575A CN102169575A (zh) | 2011-08-31 |
CN102169575B true CN102169575B (zh) | 2012-12-05 |
Family
ID=44490728
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110101912 Expired - Fee Related CN102169575B (zh) | 2011-04-22 | 2011-04-22 | 一种非基色单色光图像数据的处理方法及其处理装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102169575B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110784645B (zh) * | 2019-10-09 | 2021-01-01 | 中国航空工业集团公司洛阳电光设备研究所 | 基于单色双通道sCMOS摄像机的灰度调制图像融合方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5818975A (en) * | 1996-10-28 | 1998-10-06 | Eastman Kodak Company | Method and apparatus for area selective exposure adjustment |
US6636646B1 (en) * | 2000-07-20 | 2003-10-21 | Eastman Kodak Company | Digital image processing method and for brightness adjustment of digital images |
CN1878245A (zh) * | 2005-06-07 | 2006-12-13 | 明基电通股份有限公司 | 校正数字图像的曝光度的方法 |
CN101116324A (zh) * | 2005-02-03 | 2008-01-30 | 索尼爱立信移动通讯股份有限公司 | 用于从多次曝光中创建高动态范围图片的方法和装置 |
CN101534375A (zh) * | 2008-03-12 | 2009-09-16 | 汤姆森特许公司 | 用于校正色差的方法 |
-
2011
- 2011-04-22 CN CN 201110101912 patent/CN102169575B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5818975A (en) * | 1996-10-28 | 1998-10-06 | Eastman Kodak Company | Method and apparatus for area selective exposure adjustment |
US6636646B1 (en) * | 2000-07-20 | 2003-10-21 | Eastman Kodak Company | Digital image processing method and for brightness adjustment of digital images |
CN101116324A (zh) * | 2005-02-03 | 2008-01-30 | 索尼爱立信移动通讯股份有限公司 | 用于从多次曝光中创建高动态范围图片的方法和装置 |
CN1878245A (zh) * | 2005-06-07 | 2006-12-13 | 明基电通股份有限公司 | 校正数字图像的曝光度的方法 |
CN101534375A (zh) * | 2008-03-12 | 2009-09-16 | 汤姆森特许公司 | 用于校正色差的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102169575A (zh) | 2011-08-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
USRE47523E1 (en) | Device and method for extending dynamic range in an image sensor | |
CN105847772B (zh) | 具有透明滤波器像素的成像*** | |
Xiao et al. | High dynamic range imaging of natural scenes | |
US10511785B2 (en) | Temporally aligned exposure bracketing for high dynamic range imaging | |
US8619143B2 (en) | Image sensor including color and infrared pixels | |
TW469733B (en) | Method of dynamically modulating image brightness and image catching device using the method | |
DE102019114533A1 (de) | Erkennung von visible-light-communication-quellen über einen hohen dynamikbereich | |
TW200810566A (en) | Imaging apparatus | |
CN110365920A (zh) | 图像处理 | |
DE202017104747U1 (de) | Monochrom-Farb-Abbildung unter Verwendung eines monochromatischen Bildgebers und eines Farbabbildungssensors | |
KR102495754B1 (ko) | 광역 영상의 동적 범위를 압축하는 영상 처리 장치 및 방법 | |
CN105719274A (zh) | 边缘检测***及方法 | |
CN1288549A (zh) | 用于处理图像的方法和*** | |
CN102169575B (zh) | 一种非基色单色光图像数据的处理方法及其处理装置 | |
CN102202182B (zh) | 一种采用线阵ccd的大动态范围图像获取的装置及方法 | |
Meyer et al. | Electronic image analysis of crop residue cover on soil | |
RU2543985C1 (ru) | Способ формирования сигналов телевизионных изображений различных участков спектра | |
CN105516674B (zh) | 一种具有hdr显示功能的监视设备 | |
CN111294522A (zh) | Hdr图像成像方法、装置以及计算机存储介质 | |
CN112399150B (zh) | 一种优化监控摄像机成像画面的方法 | |
CN105913406B (zh) | 基于SoPC的新型红外图像增强处理方法 | |
KR20140057131A (ko) | 엑스선 영상 장치 및 엑스선 영상 생성 방법 | |
JP2016212412A (ja) | 被写体構造のコントラスト画像を生成する方法およびそれに関連する装置 | |
CN202059486U (zh) | 一种采用线阵ccd的大动态范围图像获取的装置 | |
CN110852981A (zh) | 一种单光子成像图像处理方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121205 |