一种基于仿射标定的多方向莫尔层析方法
技术领域
本发明涉及一种基于仿射标定的多方向莫尔层析方法。
背景技术
在现代航空、航天及能源工程等领域,复杂流场的三维显示和关键参数的定量测量是一项基本而关键的技术。光学层析技术由于具有非接触、瞬态和全场三维测量的优点而成为复杂流场诊断的重要方法。其中,莫尔层析是一种以携带待测场偏折信息的多方向投影数据重建待测场三维分布的技术,具有装置简单,抗干扰能力强,动态测量范围大,适合于在有强震动的恶劣环境下的测试的特点,对于复杂流场的定性显示和定量测试都具有很高的实用价值。由于复杂的瞬态流场往往具有非轴对称结构的特点,因此,建立多方向的光学层析装置是莫尔层析技术的研究核心之一。
文献1(J.Doi,and S.Sato,"Three-dimensional modeling of theinstantaneous temperature distribution in a turbulent flame using amultidirectional interferometer,"Opt.Eng.46,015601-015607(2007))介绍了一种基于Twyman-Green干涉仪的八个投影方向的光学层析系统。该系统每个投影方向都有物光和参考光。系统体积庞大,装置复杂,不利于三维流场诊断的实际使用。相比于非共光路的光学层析系统,文献2(B.Zhang,Z.Wu,and M.Zhao,"Deflection tomographicreconstructions of a threedimensional flame structure and temperaturedistribution of premixed combustion,"Appl.Opt.54,1341-1349(2015))公开了一种具有6个投影方向的共光路的莫尔层析系统。该系统利用反射镜和分光镜将一束激光光束分为6路,经过待测场后,利用反射镜将6路光束集中到一个探测模块,系统装置非模块化,在装调过程中会带来偏差。
基于多方向莫尔层析技术的三维流场重建需要先将多方向投影重映射到一个统一的坐标系中,而目前现有的多方向光学层析光路都忽视了这个问题,这对后续的待测场三维重建质量产生影响。
因此,需要一种新的多方向莫尔层析技术以解决上述问题。
发明内容
本发明的目的是针对在现有技术的不足,提供一种基于仿射标定的多方向莫尔层析方法。
为实现上述发明目的,本发明的基于仿射标定的多方向莫尔层析方法可采用如下技术方案:
一种基于仿射标定的多方向莫尔层析装置,包括6路独立的单方向莫尔层析装置、点阵标定板和电机,所述电机的旋转轴与6路独立的单方向莫尔层析装置均垂直,所述点阵标定板设置在所述电机的旋转轴上,所述点阵标定板与所述单方向莫尔层析装置垂直,所述单方向莫尔层析装置包括发射模块和探测模块,所述探测模块包括CCD相机,同一所述单方向莫尔层析装置的发射模块和探测模块相对设置,相邻所述单方向莫尔层析装置之间的夹角为30°,所述待测场设置在所述发射模块和探测模块之间,所述发射模块和探测模块均设置在同一个待测场的四周,所述点阵标定板和电机均设置在所述待测场中。
更进一步的,所述发射模块包括激光器、第一凸透镜L1和第二凸透镜L2,所述第一凸透镜L1和第二凸透镜L2依次设置在所述激光器的一侧,所述探测模块包括第三凸透镜L3、第一空间滤波器SF1、第四凸透镜L4、介质过滤片IF、第一光栅G1、第二光栅G2、第五凸透镜L5、第二空间滤波器SF2和第六凸透镜L6,
所述激光器正对所述CCD相机设置,所述激光器和CCD相机之间依次设置有所述第一凸透镜L1、第二凸透镜L2、第三凸透镜L3、第一空间滤波器SF1、第四凸透镜L4、介质过滤片IF、第一光栅G1、第二光栅G2、第五凸透镜L5、第二空间滤波器SF2和第六凸透镜L6,所述第一凸透镜L1和第二凸透镜L2将激光扩束准直为平行光,所述第三凸透镜L3、第一空间滤波器SF1和第四凸透镜L4构成一个开普勒望远镜系统,所述第五凸透镜L5、第二空间滤波器SF2和第六凸透镜L6构成一个光学4F系统。
更进一步的,所述第一光栅G1和第二光栅G2均为周期为50线/mm的Ronchi光栅。
更进一步的,所述点阵标定板上设置有16个圆形标定点,所述圆形标定点的直径均为0.4cm。采集得到的莫尔条纹图效果更好。
更进一步的,所述电机为步进电机。
有益效果:本发明的基于仿射标定的多方向莫尔层析装置由6路相互独立的单方向莫尔层析光路组成。本发明的多方向莫尔层析装置包含有6个发射模块、6个探测模块以及待测场区,装置模块化,装调简单,多方向莫尔层析装置强有力的支持了三维流场诊断的工程应用。
更进一步的,利用如上所述的基于仿射标定的多方向莫尔层析装置,
包括以下步骤:
1)、利用电机带动点阵标定板旋转,进行标定图像采集;建立仿射投影模型并求解,对6路独立的单方向莫尔层析装置的CCD相机的空间位置及参数的进行标定;
2)、6路独立的单方向莫尔层析装置的探测模块同时采集多方向莫尔条纹图;
3)、通过相位提取得到多方向莫尔剪切投影相位;
4)、根据步骤1)得到的CCD相机的空间位置及参数将各投影方向的图像重映射到一个统一的坐标系中;
5)、通过莫尔层析技术重建待测场。
更进一步的,建立仿射投影模型并求解,对6路独立的单方向莫尔层析装置的CCD相机的空间位置及参数的进行标定,包括以下步骤:
1)、定义重建坐标系(xct,yct,zct)和重投影坐标系(xr,yr,zr),其中,重建坐标系(xct,yct,zct)的原点Oct是点阵标定板在电机的旋转轴方向上的中点,zct与电机的旋转轴同轴,xct轴与第i个单方向莫尔层析装置一致,其中,1≤i≤6;重投影坐标系(xr,yr,zr)的原点Or与Oct一致,yr轴与zct轴方向一致,zr轴平行于第i个单方向莫尔层析装置,其中,1≤i≤6;
2)、对世界坐标系(xw,yw,zw)、摄像机坐标系(xc,yc,zc)和离散图像坐标系(u,v)进行设定,并确定世界坐标系(xw,yw,zw)、摄像机坐标系(xc,yc,zc)和离散图像坐标系(u,v)之间的转换关系:
世界坐标系(xw,yw,zw)的xw轴和yw轴均位于点阵标定板上,世界坐标系的原点Ow与重建坐标系的原点Oct一致,zw轴与所述点阵标定板的法线方向一致,yw轴与yr轴方向一致;摄像机坐标系(xc,yc,zc)的原点Oc位于旋转轴上的主点位置,zc轴与第i个单方向莫尔层析装置一致,其中,1≤i≤6,yc轴与yr轴方向一致;离散图像坐标系(u,v)位于CCD靶面上,Oc点在CCD靶面的投影点为(u0,v0),Or点在CCD靶面的投影点为(um,vm);
从世界坐标系(xw,yw,zw)到重投影坐标系(xr,yr,zr)的变换为:
其中,R为3×3正交单位旋转矩阵,rj表示旋转矩阵R的第i行,其中,j=1,2,3,t为三维平移向量;
重投影坐标系(xr,yr,zr)到摄像机坐标系(xc,yc,zc)的变换为:
其中,d1和d2定义为Oc和Or的抵消因子;
摄像机坐标系(xc,yc,zc)到离散图像坐标系(u,v)的变换是一个仿射变换,两者关系为:
其中,αx和αy是缩放因子,(u0,v0)为Oc点在CCD靶面的投影点;
将重投影坐标系(xr,yr,zr)到摄像机坐标系(xc,yc,zc)的变换导入摄像机坐标系(xc,yc,zc)到离散图像坐标系(u,v)的变换,得到离散图像坐标系(u,v)与重投影坐标系(xr,yr,zr)之间的投影关系为:
则离散图像坐标系(u,v)与世界坐标系(xw,yw,zw)之间的投影关系为:
将离散图像坐标系(u,v)与世界坐标系(xw,yw,zw)之间的投影关系改写为:
式中,
根据离散图像坐标系(u,v)与世界坐标系(xw,yw,zw)之间的投影关系的两种表达式得到:
3)、根据步骤1)和步骤2可得:
其中,r1和r2是单位向量,根据上述两式计算得到αx和αy,根据r3=r1×r2求得r3;
4)、用电机控制点阵标定板旋转360°,每间隔1°利用CCD相机采集6幅与6路独立的单方向莫尔层析装置方向一致的投影图,Or点在CCD靶面的投影点(um,vm)通过对投影图中点阵标定板的主点位置做平均来得出。其中,仿射标定大大提高相机空间位置以及内部参数测量结果的精确度,
更进一步的,步骤3)中通过相位提取得到多方向莫尔剪切投影相位包括以下步骤:
(1)、获得无待测场的多方向相位;
(2)、加入待测场,然后得到多方向莫尔条纹图;
(3)、将步骤(2)的多方向莫尔条纹图与步骤(1)的多方向相位相减,获得多方向莫尔剪切投影相位。
有益效果:本发明的基于仿射标定的多方向莫尔层析方法通过建立多方向投影仿射标定模型及其求解方法,获得各方向之间的夹角,为每一路投影在层析重建坐标中定位,同时将各方向投影图统一到同一个层析重建坐标系中。该仿射标定模型不仅适用于莫尔层析光路,也适用于其它激光干涉层析技术。
附图说明
图1是本发明多方向莫尔层析系统装置图;
图2是本发明单方向莫尔层析装置图;
图3是本发明多方向莫尔层析系统工作步骤流程图;
图4是本发明仿射标定模型的重建坐标系和重投影坐标系设置图;
图5是本发明仿射标定模型的单方向坐标设置图;
图6是实施例中6个投影方向的标定板投影图;
图7是实施例中采集的莫尔条纹图和莫尔剪切投影相位
图8是实施例中丙烷火焰三维温度场重建结果。
具体实施方式
下面结合附图和具体实施例,进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
请参阅图1、图2、图3、图4和图5所示,本发明的基于仿射标定的多方向莫尔层析装置,包括6路独立的单方向莫尔层析装置、点阵标定板和电机。优选的,电机为步进电机。点阵标定板上设置有16个圆形标定点,圆形标定点的直径均为0.4cm。采集得到的莫尔条纹图效果更好。电机的旋转轴与6路独立的单方向莫尔层析装置均垂直,点阵标定板设置在电机的旋转轴上,点阵标定板与单方向莫尔层析装置垂直,单方向莫尔层析装置包括发射模块和探测模块,探测模块包括CCD相机,同一单方向莫尔层析装置的发射模块和探测模块相对设置,相邻单方向莫尔层析装置之间的夹角为30°,待测场设置在发射模块和探测模块之间,发射模块和探测模块均设置在同一个待测场的四周,点阵标定板和电机均设置在待测场中。
发射模块包括激光器、第一凸透镜L1和第二凸透镜L2,第一凸透镜L1和第二凸透镜L2依次设置在激光器的一侧,探测模块包括第三凸透镜L3、第一空间滤波器SF1、第四凸透镜L4、介质过滤片IF、第一光栅G1、第二光栅G2、第五凸透镜L5、第二空间滤波器SF2和第六凸透镜L6,其中,第一光栅G1和第二光栅G2均为周期为50线/mm的Ronchi光栅。
激光器正对CCD相机设置,激光器和CCD相机之间依次设置有第一凸透镜L1、第二凸透镜L2、第三凸透镜L3、第一空间滤波器SF1、第四凸透镜L4、介质过滤片IF、第一光栅G1、第二光栅G2、第五凸透镜L5、第二空间滤波器SF2和第六凸透镜L6,第一凸透镜L1和第二凸透镜L2将激光扩束准直为平行光,第三凸透镜L3、第一空间滤波器SF1和第四凸透镜L4构成一个开普勒望远镜系统,第五凸透镜L5、第二空间滤波器SF2和第六凸透镜L6构成一个光学4F系统。
本发明的基于仿射标定的多方向莫尔层析装置由6路相互独立的单方向莫尔层析光路组成。本发明的多方向莫尔层析装置包含有6个发射模块、6个探测模块以及待测场区,装置模块化,装调简单,多方向莫尔层析装置强有力的支持了三维流场诊断的工程应用。
利用如上的基于仿射标定的多方向莫尔层析装置,
包括以下步骤:
1)、利用电机带动点阵标定板旋转,进行标定图像采集;建立仿射投影模型并求解,对6路独立的单方向莫尔层析装置的CCD相机的空间位置及参数的进行标定;
2)、6路独立的单方向莫尔层析装置的探测模块同时采集多方向莫尔条纹图;
3)、通过相位提取得到多方向莫尔剪切投影相位;
4)、根据步骤1)得到的CCD相机的空间位置及参数将各投影方向的图像重映射到一个统一的坐标系中;
5)、通过莫尔层析技术重建待测场。
建立仿射投影模型并求解,对6路独立的单方向莫尔层析装置的CCD相机的空间位置及参数的进行标定,包括以下步骤:
1)、定义重建坐标系(xct,yct,zct)和重投影坐标系(xr,yr,zr),其中,重建坐标系(xct,yct,zct)的原点Oct是点阵标定板在电机的旋转轴方向上的中点,zct与电机的旋转轴同轴,xct轴与第i个单方向莫尔层析装置一致,其中,1≤i≤6;重投影坐标系(xr,yr,zr)的原点Or与Oct一致,yr轴与zct轴方向一致,zr轴平行于第i个单方向莫尔层析装置,其中,1≤i≤6;
2)、对世界坐标系(xw,yw,zw)、摄像机坐标系(xc,yc,zc)和离散图像坐标系(u,v)进行设定,并确定世界坐标系(xw,yw,zw)、摄像机坐标系(xc,yc,zc)和离散图像坐标系(u,v)之间的转换关系:
世界坐标系(xw,yw,zw)的xw轴和yw轴均位于点阵标定板上,世界坐标系的原点Ow与重建坐标系的原点Oct一致,zw轴与点阵标定板的法线方向一致,yw轴与yr轴方向一致;摄像机坐标系(xc,yc,zc)的原点Oc位于旋转轴上的主点位置,zc轴与第i个单方向莫尔层析装置一致,其中,1≤i≤6,yc轴与yr轴方向一致;离散图像坐标系(u,v)位于CCD靶面上,Oc点在CCD靶面的投影点为(u0,v0),Or点在CCD靶面的投影点为(um,vm);
从世界坐标系(xw,yw,zw)到重投影坐标系(xr,yr,zr)的变换为:
其中,R为3×3正交单位旋转矩阵,rj表示旋转矩阵R的第i行,其中,j=1,2,3,t为三维平移向量;
重投影坐标系(xr,yr,zr)到摄像机坐标系(xc,yc,zc)的变换为:
其中,d1和d2定义为Oc和Or的抵消因子;
摄像机坐标系(xc,yc,zc)到离散图像坐标系(u,v)的变换是一个仿射变换,两者关系为:
其中,αx和αy是缩放因子,(u0,v0)为Oc点在CCD靶面的投影点;
将重投影坐标系(xr,yr,zr)到摄像机坐标系(xc,yc,zc)的变换导入摄像机坐标系(xc,yc,zc)到离散图像坐标系(u,v)的变换,得到离散图像坐标系(u,v)与重投影坐标系(xr,yr,zr)之间的投影关系为:
则离散图像坐标系(u,v)与世界坐标系(xw,yw,zw)之间的投影关系为:
将离散图像坐标系(u,v)与世界坐标系(xw,yw,zw)之间的投影关系改写为:
式中,m1和m2是人为设定的向量。
根据离散图像坐标系(u,v)与世界坐标系(xw,yw,zw)之间的投影关系的两种表达式得到:
3)、根据步骤1)和步骤2可得:
其中,r1和r2是单位向量,根据上述两式计算得到αx和αy,根据r3=r1×r2求得r3;
4)、用电机控制点阵标定板旋转360°,每间隔1°利用CCD相机采集6幅与6路独立的单方向莫尔层析装置方向一致的投影图,Or点在CCD靶面的投影点(um,vm)通过对投影图中点阵标定板的主点位置做平均来得出。其中,仿射标定大大提高相机空间位置以及内部参数测量结果的精确度,
步骤3)中通过相位提取得到多方向莫尔剪切投影相位包括以下步骤:
(1)、获得无待测场的多方向相位;
(2)、加入待测场,然后得到多方向莫尔条纹图;
(3)、将步骤(2)的多方向莫尔条纹图与步骤(1)的多方向相位相减,获得多方向莫尔剪切投影相位。
本发明的基于仿射标定的多方向莫尔层析方法通过建立多方向投影仿射标定模型及其求解方法,获得各方向之间的夹角,为每一路投影在层析重建坐标中定位,同时将各方向投影图统一到同一个层析重建坐标系中。该仿射标定模型不仅适用于莫尔层析光路,也适用于其它激光干涉层析技术。
实施例1
请参阅图1,一种基于仿射标定的多方向莫尔层析系统。系统由6路相互独立的单方向莫尔层析光路组成,每两个投影方向间夹角为30°,每个方向发射端使用功率为200mW,波长为532nm的固体激光器。探测端使用德国AVT公司生产的型号为AVT Guppy Pro F-125B工业数字摄像机,分辨率为1292*964。标定过程中使用的步进电机为大恒光电公司生产的型号为GCD-012060M的步进电机。点阵标定板尺寸为12cm×12cm,上面有直径为0.4cm的16个标定点。系统工作时,探测端同时采集光线经过待测场后的莫尔条纹图,莫尔层析技术的投影图像记录了探测光线经过待测场后的偏折角信息,通过相位提取得到多方向的莫尔剪切投影相位,通过建立多方向仿射投影模型并求解,实现对多个CCD相机空间位置以及相应参数的标定,同时将各方向投影图重映射到一个统一的层析重建坐标系中,最后重建出三维待测场。
请参阅图2,单方向莫尔层析光路图。透镜L1和L2将激光扩束准直为直径为50mm的平行光。透镜L3和L4以及空间滤波器SF1构成了一个开普勒望远镜系统,透镜L5和L6以及空间滤波器SF2构成了一个4F系统。探测端光路使得到达CCD探测器的光束直径缩小到2.5mm。G1和G2是两片光栅周期为50线/mm的Ronchi光栅。IF是介质滤光片。
请参阅图3,多方向莫尔层析系统的工作流程包括以下步骤:
步骤1、待测场区放置的步进电机带动点阵标定板旋转,进行标定图像采集。通过建立多方向仿射投影模型并求解,实现对多个CCD相机空间位置以及相应参数的标定。
请参阅图4,对重建坐标系和重映射坐标进行设置,所有坐标系的设置都符合右手定则。定义多方向投影图的重建坐标系为(xct,yct,zct),zct与步进电机的旋转轴方向一致,xct与第一个投影方向的光轴方向一致,原点Oct是标定板在步进电机旋转轴方向上的中点。每一个投影方向都可以定义一个重投影坐标系(xr,yr,zr)。其原定Or与Oct一致,yr与zct方向一致,zr平行于各投影方向的光轴方向。
请参阅图5,分别对世界坐标系、摄像机坐标系、离散图像坐标系进行设置,坐标系的设置也都符合右手定则,确定坐标系之间的转换关系。具体为:
定义世界坐标系(xw,yw,zw)的平面位于点阵标定板上,其原定Ow定义在标定板的中点,zw是标定板的法线方向,yw与yr方向一致。定义摄像机坐标系(xc,yc,zc)的原点Oc位于旋转轴上的主点位置,Zc与光轴方向一致,yc与yr方向一致。定义离散图像坐标系(u,v)位于CCD靶面,坐标原点位于图像左上角。Oc在CCD靶面的投影点为(u0,v0),Or在CCD靶面的投影点为(um,vm)。
根据上述坐标系的设置,则有从世界坐标系(xw,yw,zw)到重投影坐标(xr,yr,zr)的变换为:
其中R为3×3正交单位旋转矩阵,可用ri(i=1,2,3)表示旋转矩阵R的第i行。t为三维平移向量。重投影坐标系(xr,yr,zr)到摄像机坐标系(xc,yc,zc)的变换为:
其中d1和d2定义为Oc和Or的抵消因子。从摄像机坐标系(xc,yc,zc)到图像坐标系(u,v)的变换是一个仿射变换,两者关系为:
其中αx和αy是缩放因子。将(2)式代入(3)式可得:
合并上述(2),(3),(4)三式可得离散图像坐标系(u,v)与世界坐标系(xw,yw,zw)之间的投影关系为:
将公式(5)改写如下:
根据公式(5)、(6)可得:
式中,m1和m2是人为设定的向量。
综上可得:
由于r1,r2是单位向量,所以根据公式(8)、(9)可以计算出αx和αy。由于R是单位正交矩阵,根据r3=r1×r2可以求得r3。调节步进电机旋转360°,每间隔1°采集一次标定板的投影图,Or在CCD靶面的投影点为(um,vm)可以通过对投影图中标定板主点位置做平均来得出。综上所述,基于仿射标定模型的相机参数可以求得。
步骤2、多方向莫尔层析光路的探测端同时采集携带待测场信息的莫尔条纹图。
步骤3、通过相位提取得到多方向莫尔剪切投影相位。首先获得无待测场的多方向相位,然后加入待测场,从激光器发出的光线经过待测场时会发生偏折,形成莫尔条纹图。通过与待测场的相位相减,可以获得待测场的多方向莫尔剪切投影相位。
步骤4、基于步骤1求出的相机标定参数,对各投影方向的图像进行重映射。
步骤5、通过莫尔层析技术重建待测场。
下面请参阅具体实施例对本发明做进一步说明。
实施例2
请参阅图1,建立一个具有6个投影采集方向的基于仿射标定的莫尔层析系统,在丙烷预混火焰燃烧从瞬态到稳态的燃烧过程中,对三个瞬态时刻以及一个稳态时刻的温度场进行三维重建。实施例中曝光时间为20μs,采集帧频为5fbs。图6是从不同投影方向上采集到的点阵投影图。为实现上述多个方向的投影匹配过程,用电机控制点阵标定板旋转360°,每间隔1°采集6幅不同方向的投影图。再根据步骤1描述的算法计算出不同投影方向上的主点坐标(um,vm),从而得出不同方向上的内参数矩阵。各投影方向摄像内参数矩阵如表1所示:
表1各方向内参数矩阵计算结果
|
主点坐标(u<sub>m</sub>,v<sub>m</sub>) |
缩放比例α<sub>x</sub> |
缩放比例α<sub>y</sub> |
方向1 |
(639.34,464.98) |
13.00 |
12.93 |
方向2 |
(622.11,456.41) |
12.80 |
12.90 |
方向3 |
(648.48,478.42) |
13.20 |
13.19 |
方向4 |
(630.99,493.53) |
13.00 |
12.90 |
方向5 |
(687.68,476.19) |
13.10 |
13.10 |
方向6 |
(660.13,466.82) |
13.20 |
13.20 |
图7是本系统6个探测模块采集的丙烷火焰燃烧的四个时刻的莫尔条纹图以及莫尔剪切投影相位。其中前三组为丙烷火焰燃烧过程中的瞬态时刻,第四组为丙烷火焰燃烧过程中的稳态时刻。
图8是实施例中对应图7燃烧状态的丙烷火焰的温度场三维重建结果。
虽然上面已经示出和描述了本发明的实施例,但不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。