主权项 |
一种基于点集匹配的图像配准方法,用于将浮动图像与参考图像进行配准以得到配准图像,其特征在于,包括如下步骤:步骤S1:对于所述浮动图像和所述参考图像进行特征点的提取,使该浮动图像和参考图像表示为浮动图像特征点集和参考图像特征点集;步骤S2:构造12个四元数表征的空间映射信息的集合:<maths num="0001"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>Ω</mi><mo>=</mo><mo>{</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>,</mo><mn>0</mn><mo>,</mo><mn>0</mn><mo>,</mo><mn>0</mn><mo>)</mo></mrow><mi>T</mi></msup><mo>,</mo><msup><mrow><mo>(</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>,</mo><mn>0</mn><mo>,</mo><mn>0</mn><mo>)</mo></mrow><mi>T</mi></msup><mo>,</mo><msup><mrow><mo>(</mo><mn>0</mn><mo>,</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>,</mo><mn>0</mn><mo>)</mo></mrow><mi>T</mi></msup><mo>,</mo><msup><mrow><mo>(</mo><mn>0</mn><mo>,</mo><mn>0</mn><mo>,</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>)</mo></mrow><mi>T</mi></msup><mo>,</mo></mrow></mtd></mtr><mtr><mtd><mrow><msup><mrow><mo>(</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>)</mo></mrow><mi>T</mi></msup><mo>,</mo><msup><mrow><mo>(</mo><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>)</mo></mrow><mi>T</mi></msup><mo>,</mo><msup><mrow><mo>(</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>)</mo></mrow><mi>T</mi></msup><mo>,</mo><msup><mrow><mo>(</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>)</mo></mrow><mi>T</mi></msup><mo>,</mo></mrow></mtd></mtr><mtr><mtd><mrow><msup><mrow><mo>(</mo><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>)</mo></mrow><mi>T</mi></msup><mo>,</mo><msup><mrow><mo>(</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>)</mo></mrow><mi>T</mi></msup><mo>,</mo><msup><mrow><mo>(</mo><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>)</mo></mrow><mi>T</mi></msup><mo>,</mo><msup><mrow><mo>(</mo><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>)</mo></mrow><mi>T</mi></msup><mo>}</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001191627200000011.GIF" wi="1566" he="367" /></maths>绕轴<img file="FDA0001191627200000012.GIF" wi="454" he="78" />旋转α角度用四元数表示为:<img file="FDA0001191627200000013.GIF" wi="965" he="96" />ω<sub>x</sub>,ω<sub>y</sub>,ω<sub>z</sub>表示轴<img file="FDA0001191627200000014.GIF" wi="48" he="43" />坐标;对于所述浮动图像特征点集分别作用这12个空间映射,得到12个中间图像特征点集,该12个中间图像特征点集中至少存在一个点集与所述参考图像的特征点集的绕轴旋转角度小于等于90度;步骤S3:对于每个中间图像特征点集,与所述参考图像特征点集进行匹配,得到12个几何变换,该12个几何变换分别与步骤S2对应的空间映射作用后,求得所述浮动图像特征点集到参考图像特征点集的12个最终空间映射信息;步骤S4:将所述12个最终空间映射信息分别作用于所述浮动图像,比较其与参考图像的像素点位置误差,挑选最小的误差所对应的最终空间映射为满足图像配准要求的最优解。 |