一种基于连续太赫兹共焦扫描偏振成像的图像融合方法
【技术领域】
[0001] 本发明涉及基于连续太赫兹共焦扫描偏振成像的图像融合方法。
【背景技术】
[0002] 太赫兹(Terahertz，简称THz)辐射通常指的是频率在O.l-lOTHz范围内的电磁辐 射，此波段内的电磁辐射可以穿透大多数非金属非极性物质。
[0003] 偏振成像(polarization imaging)是一种通过与物质发生相互作用后的光束的 偏振信息来获得物体的偏振特性图像的技术。它可以提取出物体表面各向异性特性的信 息，从而在一定程度上可以提高普通成像的质量。
[0004] 图像融合（image fusion)是一种将多源信道所采集到的关于同一目标的图像数 据综合起来的技术。它可以最大限度的提取各自信道中的有利信息，以提高图像信息的利 用率。
[0005] 截至目前，已有可见光等波段偏振成像图像融合方面的研究，而连续太赫兹共焦 扫描偏振成像图像融合方法则未见报导。由于太赫兹辐射的特殊性，图像融合的常规算法 诸如空间滤波、频域滤波等在太赫兹波段难以奏效;并且上述常规算法往往需要基于经验 的考虑来选取参数，难以实现图像的自动化处理，有一定的局限性。另外连续太赫兹共焦扫 描偏振成像效率较低，为获得足够好的图像质量，需要采取大量的点，耗时较长，因而它的 应用受到一定的限制。
【发明内容】
[0006] 本发明是为了解决目前连续太赫兹共焦扫描偏振成像中成像效率较低以及成像 质量较差的问题，而提出的一种基于连续太赫兹共焦扫描偏振成像的图像融合方法。
[0007] -种基于连续太赫兹共焦扫描偏振成像的图像融合方法按以下步骤实现：
[0008] 步骤一:对连续太赫兹共焦扫描偏振成像中得到的太赫兹光的偏振度P进行对比 度拉伸操作，得到Pi，Pi为对比度拉伸操作后的P;
[0009] 步骤二:对步骤一中得到的PiW及连续太赫兹共焦扫描偏振成像中得到的太赫兹 光的总强度I进行插值操作，得到内和12，？ 2为插值操作后的?1，12为为插值操作后的1;
[0010] 步骤三:对步骤二中得到的内和12进行校正操作，得到p3和13,所述p 3为校正操作后 的P2,13为校正操作后的12;
[0011 ] 步骤四：对步骤三中得到的P3和13进行加权融合，得到IP，所述IP为目标扫描图。
[0012] 发明效果：
[0013] 本发明实现了自动选取参数对太赫兹光的偏振度P进行对比度拉伸操作，太赫兹 光的偏振度P的直方图分布地更为均匀，有效地提高了它的对比度。对偏振分量I、P的插值 操作弥补了太赫兹共焦扫描成像效率较低的缺点，可以通过减少扫描点数来减少扫描时 间，再通过插值操作来提高图像的质量。对I、p进行的加权融合，充分利用了物体表面偏振 情况和太赫兹波本身的信息。
【附图说明】
[0014] 图1为本发明方法的流程图；
[0015] 图2为本发明方法中太赫兹光的偏振度P的原始图像；
[0016] 图3为本发明方法中对太赫兹光的偏振度P进行对比度拉伸后得到的丹的图像； [0017]图4为本发明方法中对朽进行插值操作后得到的内的图像；
[0018] 图5为本发明方法中太赫兹光的总强度I的原始图像；
[0019] 图6为本发明方法中对太赫兹光的总强度I进行插值操作后得到的12的图像；
[0020] 图7为本发明方法中对内进行校正操作后得到的P3的图像；
[0021]图8为本发明方法中对12进行校正操作后得到的13的图像；
[0022] 图9为本发明方法中对13和P3进行加权融合后得到的IP的图像。
【具体实施方式】
【具体实施方式】 [0023] 一:如图1所示，一种基于连续太赫兹共焦扫描偏振成像的图像融合 方法包括以下步骤：
[0024] 步骤一:对连续太赫兹共焦扫描偏振成像中得到的太赫兹光的偏振度P进行对比 度拉伸操作，得到Pi，Pi为对比度拉伸操作后的P;
[0025] 步骤二:对步骤一中得到的Pi以及连续太赫兹共焦扫描偏振成像中得到的太赫兹 光的总强度I进行插值操作，得到P2和1 2,所述I为太赫兹光的总强度，P2为插值操作后的丹， 12为为插值操作后的I;
[0026] 步骤三:对步骤二中得到的内和12进行校正操作，得到P3和13,所述P 3为校正操作后 的P2,13为校正操作后的12;
[0027] 步骤四：对步骤三中得到的P3和13进行加权融合，得到IP，所述IP为（最终得到的） 目标扫描图，它结合了 I和P分量各自的优点。
【具体实施方式】 [0028] 二:本实施方式与一不同的是所述步骤一中得到 具体过程为：
[0029] 将连续太赫兹共焦扫描偏振成像中得到的太赫兹光的偏振度P的灰度级归一化至 [0，1]区间内，求取P的归一化直方图P(r k)，然后根据p(rk)求取归一化直方图累加和s(rk)， 计算公式如下：
[0032]其中所述k为灰度级的级数，k=l，2,…，L，L为[0，1]区间含有的总灰度级数，rk为 第k级亮度，η为P含有的像素数，nk为灰度级为Λ的图像P中的像素数；
[0033]求取s(rk) = 0.^trk的值，对偏振分量图像P进行对比度拉伸操作，将[a，l]区间线 性映射到[0，1 ]区间，计算公式如下：
[0035] 其中（i，j)为像素点的坐标，P(i，j)为坐标为（i，j)时的太赫兹光的偏振度，PKi， j)为对比度拉伸操作后的P(i，j)。通过此操作便提高了偏振分量图像P的对比度。
[0036] 对于一般的物体，连续太赫兹波照射在其表面反射后的光的退偏化程度较低，这 会导致偏振分量图像P大部分区域的值接近1，对比度比较低。通过这一步所述的算法，可以 舍弃40%的低灰度值部分，拉伸高灰度值部分，来提高图像的对比度。
【具体实施方式】 [0037] 三:本实施方式与一或二不同的是:所述步骤二中得 到P2和12的具体过程为：
[0038] 根据需要选取放大倍数为n，1 < η < 10,对连续太赫兹共焦扫描偏振成像中得到的 偏振分量图像I以及步骤一中得到的Ρ!进行双三值插值操作，得到Ρ2和1 2。
[0039]由于太赫兹共焦扫描偏振成像效率较低，为获得足够好的图像质量，需要采取大 量的点，耗时较长。可以在成像时采取少量的点，减少成像时间，再对图像进行双三值插值 操作，提高成像质量，这样便充分提高了效率。
【具体实施方式】 [0040] 四：本实施方式与一至三之一不同的是:所述步骤三 中得到Ρ3和13的具体过程为：
[0041] 利用Sobel检测器对步骤二中得到的12进行边缘检测，根据边缘检测的结果计算 图像需要旋转的角度，对步骤二中得到的P2