专利名称：傅立叶域光学相干层析技术中可变插值间隔的插值方法
技术领域：
本发明涉及一种插值方法，特别是一种新型的、可用于傅立叶域光学相干层析等有 插值需要的仪器的可变插值间隔的插值方法。
背景技术：
傅立叶域光学相干层析系统是一种新型的非接触、高分辨率的光学探测系统，它通过光学 干涉的方法对目标进行纵向扫描，最后通过二维或三维重建，得出目标的结构信息、多普勒 信息以及偏振信息，因此它可广泛应用于医学成像、各种工业损伤探测之中。傅立叶域光学
相干层析技术中，参考光和信号光在光学分光器3中形成干涉，干涉信号由衍射光栅9进行 分光，由透镜10聚焦到线扫描CCD11上，由CCD11将模拟信号转换为数字信号，如图2， 光谱仪8由衍射光栅9，透镜10和线扫描CCD11构成；CCD11采集到的光栅出射波长按波 长线性分布，但数据重建要求波长信息在k空间中线性分布，需要对数据进行插值。傅立叶 域光学相干层析系统中的插值方法有多种，如离散傅立叶变换补零插值，B样条拟合，直接 线性插值等，而大多数的傅立叶域光学相干层析系统采用离散傅立叶变换补零插值和直接线 性插值相结合的方法，即N点数据通过离散傅立叶变换运算之后，产生N点的频率域数据， 在高频点补充N^N点的零值，产生MgN+N的数据，然后通过反傅立叶变换得到M*N+N点 的数据，其中，M为补零的倍数，最后采用线性插值的方法得到N点的插值数据。由傅立叶
域光学相干层析系统扫描采集的一组数据向量^1 = {^12,......, }，根据传统的离散傅立叶
变换补零插的步骤为
1) 将数据进行离散傅立叶变换，得到一组新数据
义i (0 =》i (") exP(-/ ;
2) 对新数据进行补零插值，得到M倍补零后的数据
0, 其他
X(/陽顧+N), 廳-* S / S層_ 1
3
3) 对M倍补零后的数据进行反向离散傅立叶变换，得到拓展了 M+l倍后的数据；
4) 对拓展后的数据根据线性分布K空间进行线性插值，可以得到插值后的数据。
这种方法简单，成熟，但具有运算量大，远远达不到实时处理的要求，插值间隔与插值精
度由补零倍数M固定，不能随意改变插值间隔等缺点，并且由于进行离散傅立叶变换补零插 值后再进行线性插值，因而插值精度也受到了影响，这些局限性严重制约了傅立叶域光学相 干层析系统的应用。
发明内容
本发明要解决的技术问题是克服现有傅立叶域光学相干层析技术的插值方法中插值 精度低、运算速度慢，插值间隔固定且不能任意改变等问题，本发明提供一种新型的插值 方法，该方法精度高，运算速度快，并具有可变的插值精度和插值间隔等特点，能有效得提 高傅立叶域光学相干层析系统的运算速度和插值精度。
本发明解决其技术问题所采用的的技术方案是用于傅立叶域光学相干层析技术中可
变插值间隔的插值方法，其特征在于以下步骤
(1 )衍射光栅9通过透镜10映射在像数点为N的线扫描CCD11上的波长采用光谱
仪准确标定，得到一组与CCD11中每个像素对应的均匀分布的波长向量 义4d.......，4}，波长差为A;i，波长向量在CCD11上的实际位置系数为
/"ifec 1 = {";" = 1,2，......, #};
(2) 由第一个波长^和最后一个波长^，利用公式k-2兀/X可以得到CCDll首尾两 像素点的波数&和kN，利用Id和kw组成一个长度为N的线性分布的波数向量
& = {、 = & —=……J}，利用公式；i = ^反算出相应的波长向量
义2="' = 4， = 1，2,......，iV}，从而利用A;i可以计算出每个波数&对应的^'在CCD11
中的虚拟位置系数/^^2 = {&=4^ + 1; = 1,2,……,w;
A义
(3) 由CCD11采集的数据为实数以及实信号的在离散傅立叶变换中厄米对称的特 点，在离散傅立叶变换的过程中，在高频点增加一些数据点，得到一个离散傅立叶变换补
零插值传递函数7^("A) = l + J>os(f《"-O)，将/mfel = {";" = 1,2，……，iV}和
/"&jc2 = {^=^^ + 1;" = 1,2,……，7V}中不同的位置 n 和sn依次代入 △A
7F(",0 = l + 2>os(^/0-; 》到中产生一个忖"^的加权系数矩阵//歸0^ ),这样，对
插值加权系数的处理完成；
(4)傅立叶域光学相干层析系统中的CCD11采集一组纵轴扫描数据向量 x-(;^，;^,......，^}，进行数据插值，可到插值后的数据;c'—^,x^......，^}。其计算公
式如下
可以用 一 切可利用的窗函数对加权过程进行截断，根据窗长和
2 = K = +1;" = 1,2,......,获得插值起始位置Min和结束位置Max ，利用公式
△A
对原始数据进行插值，其中W (n)为所需加窗的窗函数，从而提高新插值方法的运算速 度。
在傅立叶域光学相干层析数据处理时，可以用一切可利用的窗函数对加权系数进行 截断，将加权系数/^w(n,&)中的数据进行加窗截断，减少数据处理的数据处理长度和数
据处理量。，其计算公式如下
Max
公式中W(n)为一切可利用的加窗函数，插值运算起始位置Min和结束位置Max由加窗
窗长和虚拟位置系数//7^^2 = {^=^^ + 1; = 1,2,......,iVn寻到，从而提高可变插值间隔
A义
插值运算的处理速度，而且可以将加权系数存储在计算机中，方便运算时调用，避免重 复计算。
本发明与现有技术相比有如下优点
1、可在事前提取波长和波数信息，构造与K空间相对应的非线性分布波长向量以及 这个波长向量在CCD11像素点所对应的虚拟位置系数向量，从而由传递函数计算出加权
系数矩阵/^,;V("人)，传统离散傅立叶变换补零值插的精度由补零倍数M固定，只能达到
1/(M+1)的位置精度，因为虚拟位置系数Sn的的位置不是由传统离散傅立叶变换补零插
值算法中的补零倍数M固定的，可以为计算机数据精度内的任意一个实数，从而实现可 变的插值精度和插值间隔；
2.、可利用一切可利用的加窗函数对加权系数矩阵进行加窗截取，并存储到计算中， 方便计算时调用，从而避免重复计算，由传统离散傅立叶变换补零插值方法中有一次N
点和一次M*N+N点快速傅立叶变换，因此需要*log2(AO + (M;Iog2((M +1)* AO次复
数乘法，本方法需要N4L次实数乘法，其中，N为CCD11像素点数，L为加窗函数的窗 长，从而提高插值方法的处理运算速度，提高离散傅立叶域光学相干层析系统的实时处理 能力。
图1为傅立叶域光学相千层析系统数据插值流程图2为傅立叶域光学相干层析系统结构图，其中1为光源，2为光学隔离器，3光学 分光器，4为偏振控制器，5为PZT转换器，6为扫描控制器，7为采样目标，8为光谱 仪，9为衍射光栅，IO为透镜，ll为线扫描CCD;
图3为插值效果比较图4为二维图像重构。 具体实施方案
下面结合附图及具体实施方式
详细介绍本发明。本实施方案采用的是傅立叶域光学相干 层析系统，对其采集到的数据进行插值，其流程如图l所示；具体步骤如下
(1) 本实施方案从图2所示的光谱仪准确标定入射到CCD11的波长，中心波长为 849.72mn，光谱分辨率A义-0.0674nm，线扫描CCD11上像素点数N=2048, CCDU上 首尾两个像素点的波长分别为"-780.7024nm， lw=918.4006nm，各个波长在CCD11中
位置系数为/"lxl = = 1,2,……,AO :
(2) 由;t:^可得CCD11首尾两个像素点所对应的两个波数为H-^和^-^，
义 义l 义jv
计算一个在K空间线性分布的波数向量f = {^=&十%^*("-1);" = 1,2,……,}，由这 个在K空间线性分布的波数向量通过公式A =^得到 一 组非均匀的
"-"'A',……,H从而可以利用公式^^^S^计算波长义'-W，A',……,；'}对应在
CCD11中的虚拟位置系数/^^2 = ^ =^^ + 1;" = 1,2,......，W;
A义
(3) 由实际波长在CCD11上位置系数向量/"1jc1 = {";" = 1,2,......，A^和在CCD11的
虚拟位置系数向量/"^12 = { =4^ + 1;" = 1，2,......,W，依次抽取不同的n和sn由传递
A义
2 卩:W
函数IF("A)^l + J]c0S(苦/("-,得到加权系数矩阵^w,w("A);
,=1 W
(4) 由图2所示的傅立叶域光学相干层析系统的CCD11采集到的一组干涉信号数据 为^^k^,^2,……,^v}，通过加窗长度为L=ll的布莱克曼即Blackman窗函数
『(")=0.42 + 0.5*(^(，) + 0.08*003(，)截断加权系数，由插值公式得到插值之后的数
据，其计算公式如下
Sn由/丄2-^二^^+l;"-l^……，M给出，iWrox: = /"cfecO)+^:^,M/" = /raflfec(^) —
A/1 2 2
(5) 对傅立叶域光学相干层析系统中CCD11采集到数据重复步骤(3)进行插值， 由离散傅立叶变换对插值后的各组数据x'(s)进行傅立叶变换得到X，(s)，令对比度 Contrast-6，亮度偏差Brightness-—82，对X，(s)中的各个点进行对数运算可以得到图像的 灰度值Intensity,其计算公式如下
/她肌'z^ = CoWraW * (10 * log 10(X 'W + Sn'g/rf"ew) + 255 计算出灰度值还要进行截取，少于0的赋值为0，高于255的赋值为255，其目的是使其 灰度范围在
之间，以符合计算机图像的灰度输出范围,由扫描控制器6对采样目标 7重复线扫描，并对CCD11采集到的干涉信号数据中进行插值和映射，可以重构出二维 图像或都三维图像，图4为一张重构的二维图像。
为了对比，在实施本方案的同时，运用传统的离散傅立叶变换补零插值方法，在本次 实验中，取一组线扫描数据4倍抽取后再进行插值，插值后的数据如图3所示，与原始 数据相减后，新方法均值为0.1409，方差为0.2524，传统的离散傅立叶变换补零方法均值 为0.1448,方差为0.2564。可见均值和方差的对比均为基于可变间隔的插值方法的为佳。 用傅立叶域光学相干层析系统对一物体进行扫描，采集到2048*300的数据进行处理及重 构图像，重构后的图像如图4。在CPU为酷睿Q9300，内存为4GB时，运算时间由原来的9秒减少到目前400毫秒，大大地提高了处理速度。
虽然通过参照本方法的说明和具体实施方案，已经对本方法进行了图示和描述，但普 通的技术人员应该明白，可以在形式上和细节上对其作各种各样的改变，而不偏离所附 权利要求书所限定的本方法的精神和范围。
权利要求
1、傅立叶域光学相干层析系统中可变插值间隔的插值方法，其特征在于包括以下步骤(1)将衍射光栅(9)通过透镜(10)映射在像素点为N的线扫描CCD(11)上的波长准确标定，得到一组均匀分布的波长向量<maths id="math0001" num="0001" ><math><![CDATA[ <mrow><msub> <mover><mi>λ</mi><mo>→</mo> </mover> <mn>1</mn></msub><mo>=</mo><mo>{</mo><msub> <mi>λ</mi> <mn>1</mn></msub><mo>,</mo><msub> <mi>λ</mi> <mn>2</mn></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub> <mi>λ</mi> <mi>N</mi></msub><mo>}</mo><mo>,</mo> </mrow>]]></math> id="icf0001" file="A2008101191310002C1.tif" wi="37" he="6" top= "54" left = "95" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>波长差为Δλ，波长向量在CCD(11)上的实际位置系数为<maths id="math0002" num="0002" ><math><![CDATA[ <mrow><mover> <mi>Index</mi> <mo>→</mo></mover><mn>1</mn><mo>=</mo><mo>{</mo><mi>n</mi><mo>;</mo><mi>n</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>}</mo><mo>;</mo> </mrow>]]></math> id="icf0002" file="A2008101191310002C2.tif" wi="48" he="6" top= "66" left = "77" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>(2)利用波长和波数的关系，将在波数空间均匀分布的波数向量反算出一组非均匀分布的波长<maths id="math0003" num="0003" ><math><![CDATA[ <mrow><msub> <mover><mi>λ</mi><mo>→</mo> </mover> <mn>2</mn></msub><mo>=</mo><mo>{</mo><msub> <msup><mi>λ</mi><mo>′</mo> </msup> <mn>1</mn></msub><mo>,</mo><msup> <msub><mi>λ</mi><mn>2</mn> </msub> <mo>′</mo></msup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msup> <msub><mi>λ</mi><mi>N</mi> </msub> <mo>′</mo></msup><mo>}</mo><mo>,</mo> </mrow>]]></math> id="icf0003" file="A2008101191310002C3.tif" wi="42" he="6" top= "85" left = "42" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>并与<maths id="math0004" num="0004" ><math><![CDATA[ <mrow><msub> <mover><mi>λ</mi><mo>→</mo> </mover> <mn>1</mn></msub><mo>=</mo><mo>{</mo><msub> <mi>λ</mi> <mn>1</mn></msub><mo>,</mo><msub> <mi>λ</mi> <mn>2</mn></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub> <mi>λ</mi> <mi>N</mi></msub><mo>}</mo> </mrow>]]></math> id="icf0004" file="A2008101191310002C4.tif" wi="34" he="6" top= "86" left = "97" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>利用Δλ进行比较得到 id="icf0005" file="A2008101191310002C5.tif" wi="3" he="5" top= "85" left = "176" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/>各个波长在CCD(11)中的虚拟位置系数<maths id="math0005" num="0005" ><math><![CDATA[ <mrow><mover> <mi>Index</mi> <mo>→</mo></mover><mn>2</mn><mo>=</mo><mo>{</mo><msub> <mi>s</mi> <mi>n</mi></msub><mo>;</mo><mi>n</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>}</mo><mo>;</mo> </mrow>]]></math> id="icf0006" file="A2008101191310002C6.tif" wi="50" he="6" top= "97" left = "86" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>(3)由于CCD(11)采集到的干涉信号数据为实数以及实信号在离散傅立叶变换中厄米对称性的特点，得到一个离散傅立叶变换补零插值传递函数<maths id="math0006" num="0006" ><math><![CDATA[ <mrow><mi>TF</mi><mrow> <mo>(</mo> <mi>n</mi> <mo>,</mo> <msub><mi>s</mi><mi>n</mi> </msub> <mo>)</mo></mrow><mo>=</mo><mn>1</mn><mo>+</mo><munderover> <mi>Σ</mi> <mrow><mi>i</mi><mo>=</mo><mn>1</mn> </mrow> <mfrac><mi>N</mi><mn>2</mn> </mfrac></munderover><mi>cos</mi><mrow> <mo>(</mo> <mfrac><mrow> <mn>2</mn> <mi>π</mi></mrow><mi>N</mi> </mfrac> <mi>i</mi> <mrow><mo>(</mo><mi>n</mi><mo>-</mo><msub> <mi>s</mi> <mi>n</mi></msub><mo>)</mo> </mrow> <mo>)</mo></mrow><mo>,</mo> </mrow>]]></math> id="icf0007" file="A2008101191310002C7.tif" wi="60" he="13" top= "125" left = "19" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>由<maths id="math0007" num="0007" ><math><![CDATA[ <mrow><mover> <mi>Index</mi> <mo>→</mo></mover><mn>1</mn><mo>=</mo><mo>{</mo><mi>n</mi><mo>;</mo><mi>n</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>}</mo> </mrow>]]></math> id="icf0008" file="A2008101191310002C8.tif" wi="46" he="6" top= "129" left = "88" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>和<maths id="math0008" num="0008" ><math><![CDATA[ <mrow><mover> <mi>Index</mi> <mo>→</mo></mover><mn>2</mn><mo>=</mo><mo>{</mo><msub> <mi>s</mi> <mi>n</mi></msub><mo>;</mo><mi>n</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>}</mo> </mrow>]]></math> id="icf0009" file="A2008101191310002C9.tif" wi="48" he="6" top= "129" left = "140" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>在TF(n，sn)中产生一个N*N的加权系数HN*N(n，sn)，这样，对插值加权系数的处理完成；(4)CCD(11)采集一组干涉信号扫描数据向量x＝{xn1，xn2，......，xnN}，利用加权系数加窗截断后对数据进行插值运算，可得到插值后的数据x′＝{xs1，xs2，......，xsN}。
2、根据权利要求1中所述的傅立叶域光学相干层析系统中可变插值间隔的插值方法， 其特征在于在傅立叶域光学相干层析数据处理时，可以用一切可利用的窗函数对加权 系数进行截断，将加权系数/^w("A)中的数据进行加窗截断，减少数据处理的数据处理 长度和数据处理量，其计算公式如下公式中W(n)为一切可利用的加窗函数，插值运算起始位置Min和结束位置Max由加窗 窗长和虚拟位置系数//^"2 = ";"-1，2，......,A^得到，从而提高插值运算的处理速度。
全文摘要
傅立叶域光学相干层析技术中可变插值间隔的插值方法，其特征在于事先通过傅立叶域光学相干层析系统特性对波长信息进行提取，由在波数k空间均匀分布的波数得到一组在波长空间非均匀分布的波长向量，从而得到这个波长向量在CCD中的虚似位置系数，通过对离散傅立叶补零插值方法传递函数产生加权系数矩阵，在系统运行时由加权系数矩阵和采集到的数据进行插值或者对数据矩阵进行加窗截断后进行插值，得到所要求的插值数据。本插值方法简单、易于实现，得到一个更高精度，更快运算速度的插值方法，从而提高傅立叶域光学相干层析数据处理的精度和速度，提高了傅立叶域光学相干层析系统的实时处理能力。
文档编号G01J3/28GK101358879SQ20081011913
公开日2009年2月4日 申请日期2008年8月27日 优先权日2008年8月27日
发明者史国华, 张雨东, 李喜琪, 凌 魏 申请人:中国科学院光电技术研究所