发明名称 股骨三维模型可视化方法
摘要 股骨三维模型可视化方法,涉及图像处理领域,它为了解决现有采用直接数字X光图像DR无法获取三维空间信息,而采用CT由于体积庞大不适于在突发情况下的骨折影像显示的问题,它包括以下步骤:步骤一、建立通用股骨三维模型;步骤二、对待建模的股骨正位DR图像和股骨侧位DR图像进行预处理;步骤三、确定正位图像的配准变换关系,确定侧位图像的配准变换关系;步骤四、实现待建模的股骨的三维姿态估计,完成股骨三维模型姿态的可视化。本发明利用双平面DR图像实现股骨的三维模型可视化,成像迅速。适用于股骨三维影像显示。
申请公布号 CN102651145B 申请公布日期 2014.11.05
申请号 CN201210099633.3 申请日期 2012.04.06
申请人 哈尔滨工业大学;哈尔滨理工大学 发明人 张斌;孙金玮;孙绍斌;孟尚
分类号 G06T17/10(2006.01)I;G06T19/20(2011.01)I 主分类号 G06T17/10(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 张果瑞
主权项 股骨三维模型可视化方法,其特征在于,它包括以下步骤:步骤一、建立通用股骨三维模型,通过通用股骨三维模型的正位投影获取通用股骨的正位股骨三维模型投影轮廓,通过通用股骨三维模型的侧位投影获取通用股骨的侧位股骨三维模型投影轮廓;步骤二、对待建模的股骨正位DR图像和股骨侧位DR图像进行预处理,提取股骨正位DR图像中的股骨边缘轮廓获得正位股骨边缘轮廓,提取股骨侧位DR图像中的股骨边缘轮廓获得侧位股骨边缘轮廓;步骤三、对步骤二提取的待建模的股骨正位DR图像的边缘轮廓与步骤一所述的通用正位股骨三维模型投影轮廓进行仿射配准,确定正位图像的配准变换关系;对步骤二提取的待建模的股骨侧位DR图像的边缘轮廓与步骤一所述的通用侧位股骨三维模型投影轮廓进行仿射配准,确定侧位图像的配准变换关系;正位图像的配准变换关系和侧位图像的配准变换关系包含了八个自由度的配准参数{α,β,t<sub>x</sub>,t<sub>y</sub>,t<sub>z</sub>,s<sub>x</sub>,s<sub>y</sub>,s<sub>z</sub>},其中α为绕X轴的旋转角度、β为绕Y轴的旋转角度、t<sub>x</sub>为沿X轴的平移量、t<sub>y</sub>为沿Y轴的平移量、t<sub>z</sub>为沿Z轴的平移量、s<sub>x</sub>为沿X轴的比例系数、s<sub>y</sub>为沿Y轴的比例系数和s<sub>z</sub>为沿Z轴的比例系数;变换关系的获取过程为:步骤三一一、由待建模的股骨边缘轮廓上的点构成参考图像点集P,P={p<sub>1</sub>,p<sub>2</sub>,...,p<sub>i</sub>},i=1,2,...n,n为正整数;由通用股骨三维模型投影轮廓上的点构成浮动图像点集Q,Q={q<sub>1</sub>,q<sub>2</sub>,...,q<sub>i</sub>},并对两个点集进行采样,保证数据量相等;步骤三一二、将参考图像点集P与浮动图像点集Q代入公式(15)中,依据最小二乘准则进行参数优化,获得α,β,t<sub>x</sub>,t<sub>y</sub>,t<sub>z</sub>,s<sub>x</sub>,s<sub>y</sub>,s<sub>z</sub>共八个自由度的配准变换参数;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mrow><mo>(</mo><mi>&alpha;</mi><mo>,</mo><mi>&beta;</mi><mo>,</mo><msub><mi>t</mi><mi>x</mi></msub><mo>,</mo><msub><mi>t</mi><mi>y</mi></msub><mo>,</mo><msub><mi>t</mi><mi>z</mi></msub><mo>,</mo><msub><mi>s</mi><mi>x</mi></msub><mo>,</mo><msub><mi>s</mi><mi>y</mi></msub><mo>,</mo><msub><mi>s</mi><mi>z</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>arg</mi><mi></mi><mi>min</mi><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mrow><mo>|</mo><mo>|</mo><msub><mi>p</mi><mi>i</mi></msub><mo>-</mo><mrow><mo>(</mo><msub><mi>S</mi><mrow><mn>2</mn><mo>&times;</mo><mn>2</mn></mrow></msub><msub><mi>R</mi><mrow><mn>2</mn><mo>&times;</mo><mn>2</mn></mrow></msub><msub><mi>q</mi><mi>i</mi></msub><mo>+</mo><msub><mi>T</mi><mrow><mn>2</mn><mo>&times;</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>|</mo><mo>|</mo></mrow><mn>2</mn><mn>2</mn></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000533106130000011.GIF" wi="1612" he="174" /></maths>其中S<sub>2×2</sub>为二维比例变换矩阵,R<sub>2×2</sub>为二维旋转变换矩阵,T<sub>2×1</sub>为二维平移向量,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>q</mi><mi>i</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mi>i</mi></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000533106130000012.GIF" wi="198" he="167" /></maths>为浮动点集中某点的二维坐标向量,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>p</mi><mi>i</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>a</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><msub><mi>b</mi><mi>i</mi></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000533106130000013.GIF" wi="197" he="160" /></maths>为参考点集中某点的二维坐标向量;对于股骨的正位图像,得到正位图像的配准变换关系的具体表达式为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>S</mi><mi>AP</mi></msub><msub><mi>R</mi><mi>AP</mi></msub><msub><mi>q</mi><mi>i</mi></msub><mo>+</mo><msub><mi>T</mi><mi>AP</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>s</mi><mi>y</mi></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>s</mi><mi>z</mi></msub></mtd></mtr></mtable></mfenced><mfenced open='[' close=']'><mtable><mtr><mtd><mi>cos</mi><mi>&alpha;</mi></mtd><mtd><mo>-</mo><mi>sin</mi><mi>&alpha;</mi></mtd></mtr><mtr><mtd><mi>sin</mi><mi>&alpha;</mi></mtd><mtd><mi>cos</mi><mi>&alpha;</mi></mtd></mtr></mtable></mfenced><msub><mi>q</mi><mi>i</mi></msub><mo>+</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>t</mi><mi>y</mi></msub></mtd></mtr><mtr><mtd><msub><mi>t</mi><mi>z</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>16</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000533106130000021.GIF" wi="1344" he="202" /></maths>对于股骨的侧位图像,得到侧位图像的配准变换关系的具体表达式为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>S</mi><mi>LAT</mi></msub><msub><mi>R</mi><mi>LAT</mi></msub><msub><mi>q</mi><mi>i</mi></msub><mo>+</mo><msub><mi>T</mi><mi>LAT</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>s</mi><mi>x</mi></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>s</mi><mi>z</mi></msub></mtd></mtr></mtable></mfenced><mfenced open='[' close=']'><mtable><mtr><mtd><mi>cos</mi><mi>&beta;</mi></mtd><mtd><mi>sin</mi><mi>&beta;</mi></mtd></mtr><mtr><mtd><mo>-</mo><mi>sin</mi><mi>&beta;</mi></mtd><mtd><mi>cos</mi><mi>&beta;</mi></mtd></mtr></mtable></mfenced><msub><mi>q</mi><mi>i</mi></msub><mo>+</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>t</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><msub><mi>t</mi><mi>z</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>17</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000533106130000022.GIF" wi="1337" he="205" /></maths>步骤四、根据步骤三获得的正位图像的配准变换关系和侧位图像的配准变换关系确定股骨三维模型的姿态,实现股骨的三维模型的可视化。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号