发明名称 超声图像三维重建方法与系统
摘要 本发明提供一种超声图像三维重建方法,其包括以下步骤:投射步骤,用于确定光线的起始位置及采样的起始点及采样点的坐标并给灰度值、阻光度值赋初值;查表步骤,用于根据当前采样点的灰度值,判断当前采样点是否满足合成条件,对于满足合成条件的采样点,查表得到当前采样点的阻光度以及颜色值;合成步骤,对满足合成条件的采样点,计算当前采样点到所投射光线起始点以及投射光线结束到开始的距离,以及基于该距离计算累积的灰度和累积的阻光度,并根据所述的查表步骤中得到的当前采样点的阻光度以及颜色值,进行组合计算得到合成光线颜色值。
申请公布号 CN101996419B 申请公布日期 2012.11.07
申请号 CN201010546441.3 申请日期 2010.11.16
申请人 深圳市理邦精密仪器股份有限公司 发明人 凌锋;齐保春
分类号 G06T17/00(2006.01)I 主分类号 G06T17/00(2006.01)I
代理机构 深圳市科吉华烽知识产权事务所 44248 代理人 胡吉科;许建
主权项 1.一种超声图像三维重建方法,其包括以下步骤:投射步骤,用于确定光线的起始位置及采样的起始点及采样点的坐标并给灰度值、阻光度值赋初值;查表步骤,用于根据当前采样点的灰度值,判断当前采样点是否满足合成条件,对于满足合成条件的采样点,查表得到当前采样点的阻光度以及颜色值;合成步骤,对满足合成条件的采样点,计算当前采样点到所投射光线起始点以及投射光线结束到开始的距离,以及基于该距离计算累积的灰度和累积的阻光度,并根据所述的查表步骤中得到的当前采样点的阻光度以及颜色值,进行组合计算得到合成光线颜色值,具体包括:对于满足合成条件的第i采样点,计算采样点到所投射光线起始处的距离D<sub>pi</sub>,计算投射光线结束到开始的长度D,并计算两者的比值D/D<sub>pi</sub>;根据第i个采样点的颜色值C<sub>pi</sub>、阻光度值<img file="FDA00001680164700011.GIF" wi="78" he="54" />以及根据距离因素,得到的颜色值C′<sub>pi</sub>、阻光度值<img file="FDA00001680164700012.GIF" wi="78" he="56" />计算合成的I<sub>i</sub>,I′<sub>i</sub>,<img file="FDA00001680164700013.GIF" wi="66" he="49" /><maths num="0001"><![CDATA[<math><mrow><msub><mi>I</mi><mi>i</mi></msub><mo>=</mo><msub><mi>I</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mo>&PartialD;</mo><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mo>&PartialD;</mo><mi>pi</mi></msub><mo>&CenterDot;</mo><msub><mi>C</mi><mi>pi</mi></msub></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msubsup><mi>I</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>=</mo><msubsup><mi>I</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow><mo>&prime;</mo></msubsup><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mo>&PartialD;</mo><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><msubsup><mo>&PartialD;</mo><mi>pi</mi><mo>&prime;</mo></msubsup><mo>&CenterDot;</mo><msubsup><mi>C</mi><mi>pi</mi><mo>&prime;</mo></msubsup><mo>;</mo></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msub><mo>&PartialD;</mo><mi>i</mi></msub><mo>=</mo><msub><mo>&PartialD;</mo><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mo>&PartialD;</mo><mi>pi</mi></msub><mo>&CenterDot;</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mo>&PartialD;</mo><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msubsup><mo>&PartialD;</mo><mi>pi</mi><mo>&prime;</mo></msubsup><mo>=</mo><mrow><mo>(</mo><mi>D</mi><mo>/</mo><msub><mi>D</mi><mi>pi</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mo>&PartialD;</mo><mi>pi</mi></msub></mrow></math>]]></maths>C′<sub>pi</sub>=C<sub>tf</sub>(D<sub>pi</sub>)其中,I<sub>i</sub>与I<sub>i-1</sub>分别为前i个采样点、i-1个采样点合成的灰度值,I′<sub>i</sub>与I′<sub>i-1</sub>分别为前i个采样点、i-1个采样点根据距离因素合成的灰度值,<img file="FDA00001680164700018.GIF" wi="66" he="49" /><img file="FDA00001680164700019.GIF" wi="61" he="49" />分别为前i个采样点、i-1个采样点而合成的阻光度,D为当前光线的长度,D<sub>pi</sub>为第i个采样点到光线开始的距离,C<sub>tf</sub>为灰度-颜色值表,C′<sub>pi</sub>为第i个采样点根据D<sub>pi</sub>查找灰度-颜色值表所得到的颜色值大小,计算合成后的光线颜色I<sub>end</sub>=a·I<sub>i</sub>+b·I′<sub>i</sub>,其中,a,b为调节系数。
地址 518000 广东省深圳市南山区蛇口南海大道1019号南山医疗器械园B栋三楼