CN102283675A - 一种医学超声宽景成像中的旋转判断及误差纠正方法 - Google Patents
一种医学超声宽景成像中的旋转判断及误差纠正方法 Download PDFInfo
- Publication number
- CN102283675A CN102283675A CN 201110139790 CN201110139790A CN102283675A CN 102283675 A CN102283675 A CN 102283675A CN 201110139790 CN201110139790 CN 201110139790 CN 201110139790 A CN201110139790 A CN 201110139790A CN 102283675 A CN102283675 A CN 102283675A
- Authority
- CN
- China
- Prior art keywords
- registration
- rotation
- template
- standard deviation
- anglec
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明公开了一种医学超声宽景成像中的旋转判断及误差纠正方法,主要包括模板配准结果标准差的统计,基于该标准差的图像旋转判断,根据旋转角度大小降低图像序列配准的间隔帧数。具体是采用多模板的方式配准两张超声图像的同时,统计每个模板配准结果的标准差;配准图像之间的相对旋转角度越大,该标准差就越大,故可以用该标准差来作为图像之间旋转角度大小的判断;对旋转角度较大的时候,通过降低配准序列配准的间隔帧数来提高配准的准确度。本发明创新性的提出基于模板配准结果标准差的旋转判断方法,提高了实时超声宽景成像中图像序列帧之间存在较大角度的旋转时的配准准确度。
Description
技术领域
本发明涉及计算机辅助医学成像技术,尤其涉及一种医学超声宽景成像中的旋转判断及误差纠正方法。
背景技术
医学超声诊断是现代影像医学的重要组成部分,在人身体的很多组织(比如心脏、肝脏、神经血管等)的病理诊断中有很大的价值。然而由于受超声探头宽度和扫描角度的限制,实时超声扫描视野有限,不能反映被观察组织的整体图像,这影响了医生诊断的效果。针对这问题,人们提出了超声宽景成像技术。它是通过探头的移动获取一系列二维切面图像,利用计算机重建的方法把这一系列二维图像拼接为一幅连续超宽视野的切面图像。现在的宽景成像技术主要包括图像配准、图像几何变换、全景图重建和图像显示四部分。图像配准是宽景成像的核心。传统的超声宽景成像技术在配准两张图像时,一般用多模板的方法完成图像的配准;另外,由于宽景成像中探头的扫面速度是很慢的,图像序列相邻帧之间的变化很小,故一般采用相隔多帧配准一次的方法配准整个图像序列。模板配准中,一般在水平和垂直方向上搜索模板配准位置。这种方法只适合图像之间的旋转角度较小的时候。当该旋转角度太大的时候,这种方法的准确率不高。特别是在对整个图像序列相隔多帧配准一次的时候:相隔帧数越多,图像之间的相对旋转角度就可能越大,配准的准确度可能越低。故对配准图像之间的旋转角度的估计是很重要的。
发明内容
本发明的目的在于克服现有技术的缺点和不足,提供一种医学超声宽景成像中的旋转判断及误差纠正方法;克服现有宽景成像时配准图像之间相对旋转角度较大的时候配准准确率不高的问题,在估计出图像之间旋转角度的基础上,对图像序列相隔帧数进行修正,提高整个图像序列的配准准确度。
本发明通过下述技术方案实现:
一种医学超声宽景成像中的旋转判断及误差纠正方法,其步骤如下:
(1)、统计模板配准结果的标准差;
(2)、根据模板配准结果标准差来判断配准图像之间旋转角度的大小;
(3)、根据旋转角度的大小,相应的修正图像序列的间隔帧数。
所述步骤(1)采用多模板的方式配准两帧超声图像,并统计模板配准结果的标准差,具体步骤有:在待配准图像中选择多个模板,在水平和垂直方向搜索每个模板在固定图像中的匹配位置,得到每个模板的配准结果(Rb(i),i=1...Nb,Nb是模板的个数)。这种搜索方式的搜索范围中没有加入旋转角度,故只适合图像之间相对旋转角度比较小的情况。当旋转角度较大时,每个模板的配准结果准确度会有所降低,而且每个模板配准结果之间的差异也会比较大。故可以用模板配准结果的标准差来作为旋转角度大小的度量。所有模板配准结果的标准差(SDr)的计算如下:
所述步骤(2)根据模板配准结果标准差来判断配准图像之间旋转角度的大小。配准图像之间的相对旋转越大的时候,模板配准结果标准差就越大。当模板配准结果标准差大于阈值TSD的时候,就确定配准图像之间的旋转角度过大。本发明中,通过实验来确定阈值TSD的具体值,过程如下:
根据不同的参数值(Δx,Δy,θ)在一张大图中截取两张大小为360×360的小图(如图2中矩形框所示)。然后用基于多模板的配准方法配准这两张图,得到配准结果为(Δx′,Δy′,θ′)。由于两图之间的实际变换参数(Δx,Δy,θ)是已知的,则配准结果中水平方向、垂直方向、旋转角度的误差分别为:
εx=Δx′-Δx
εy=Δy′-Δy
εθ=θ′-θ
实验中,Δx分别取0、10、20、30、40、50;Δy一直取0;θ取的范围为0~10°,每0.01°取一个值。Δx取不同的值,得到的结果很类似。这里以Δx=0为例子进行说明,图3是Δx=0时候的配准结果误差和标准差图。从图3(a)、(b)中可见,当θ小于一定角度(5~6°,设为θ0)时,εx、εy、εθ的值都比较小,当θ大于θ0时,εx、εy、εθ的值都突然变大,意味着配准结果突然变差,甚至是错误的。从图3(c)可见,SDblock的值在θ大于θ0时也突然变大,变大前的SDblock值大致为10,故这里阈值TSD的值设为10。
所述步骤(3)根据旋转角度的大小,相应的修正图像序列的间隔帧数。本发明中,在配准整个图像序列的时候采用的相隔多帧配准一次的方法,以减小计算量,提高成像速度。当图像序列对应的组织旋转比较明显时,相隔帧数越大,配准图像之间的相对旋转角度越大,配准的准确度越低。本发明中,在完成一次配准之前,首先根据模板配准结果的标准差来判断图像之间的旋转角度是否过大,假如角度较小则该次配准结果有效,继续下一次配准;假如角度过大,则降低相隔帧数,重新进行配准。具体为:设配准失败的两帧图像序号分别为m和n,则要重新进行两次配准,配准的图像帧分别为m和(m+n)/2、(m+n)/2和n。假如这两次配准依旧失败,则继续折半配准,直到配准成功。
经过上述步骤,完成对整个超声图像序列配准的旋转判断,优化配准过程,并得到一张大的宽景图像。
本发明首次提出了基于模板配准结果标准差的旋转判断的新技术。根据旋转较大的时候,模板配准结果准确性不高及结果标准差较大两个特点,用该标准差作为旋转大小的一个度量。当旋转角度较大的时候,就降低图像序列配准的相隔帧数,从而降低误差,提高了实时超声宽景成像中图像序列帧之间存在较大角度的旋转时的配准准确度。
附图说明
图1是本发明医学超声宽景成像中的旋转判断及误差纠正方法的示意框图;其中:图(a)为两帧图像的配准结果的旋转判断示意框图;图(b)为多帧图像序列的配准的误差纠正示意框图。
图2是本发明医学超声宽景成像中的旋转判断及误差纠正方法,从大图中截取两个小图的示意图。
图3是本发明发明医学超声宽景成像中的旋转判断及误差纠正方法中,对不同的旋转角度的曲线图,横坐标是旋转角度,其中纵坐标:图(a)为εx、εy;图(b)为εθ;图(c)为SDblock。
图4是用于本发明的实验装置示意图。
图5是本发明方法改进前的宽景图。
图6是本发明方法改进后的宽景图。
具体实施方式
下面结合具体实施例对本发明作进一步具体详细描述,但本发明的实施方式不限于此,对于未特别注明的工艺参数,可参照常规技术进行。
实施本发明需要超声设备仪采集到的AVI超声图像视频文件,用计算机从AVI文件中提取出图像数据进行宽景成像,用纯平型显示器显示用户图形界面,可采用C语言编制各类处理程序,便能较好地实施本发明。
实施例
如图1所示,配准两张图的时候,采用多模板的配准方法,并统计每个模板配准结果的标准差,并用该标准差作为配准图像之间旋转角度大小的度量;在配准整个图像序列的时候,本发明采用相隔多帧配准一次的方法,并根据图像之间旋转角度的大小来调整相隔帧数,从而提高宽景成像的准确度。
本发明医学超声宽景成像中的旋转判断及误差纠正方法,通过下述步骤实现:
(1)、统计模板配准结果的标准差;
(2)、根据模板配准结果标准差来判断配准图像之间旋转角度的大小;
(3)、根据旋转角度的大小,相应的修正图像序列配准的间隔帧数。
如图1(a)所示。上述步骤(1)采用多模板的方式配准两帧超声图像,并统计模板配准结果的标准差,具体操作为:在待配准图像中选择多个模板,在水平和垂直方向搜索每个模板在固定图像中的匹配位置,得到每个模板的配准结果:Rb(i),i=1...Nb,Nb是模板的个数,然后计算所有模板配准结果的标准差(SDr):
上述步骤(2)根据模板配准结果标准差来判断配准图像之间旋转角度的大小,当模板配准结果标准差大于阈值TSD的时候,就确定配准图像之间的旋转角度过大;阈值TSD的选择过程如下:
如图2所示,根据不同的参数值(Δx,Δy,θ)在一张大图中截取两张大小为360×360的小图,然后用基于多模板的配准方法配准这两张图,得到配准结果为(Δx′,Δy′,θ′)。由于两图之间的实际变换参数(Δx,Δy,θ)是已知的,则配准结果中水平方向、垂直方向、旋转角度的误差分别为:
εx=Δx′-Δx
εy=Δy′-Δy
εθ=θ′-θ
实验中,Δx分别取0、10、20、30、40、50;Δy一直取0;θ取的范围为0~10°,每0.01°取一个值。Δx取不同的值,得到的结果很类似。这里以Δx=0为例子进行说明,图3是Δx=0时候的配准结果误差和标准差图。从图3(a)、(b)中可见,当θ小于一定角度(5°~6°,设为θ0)时,εx、εy、εθ的值都比较小,当θ大于θ0时,εx、εy、εθ的值都突然变大,意味着配准结果突然变差,甚至是错误的。从图3(c)可见,SDblock的值在θ大于θ0时也突然变大,变大前的SDblock值大致为10,故这里阈值TSD的值设为10。
上述步骤(3)根据旋转角度的大小,相应的修正图像序列的间隔帧数,具体实施方法是:如图1(b)所示,在配准整个图像序列的时候采用的相隔多帧配准一次的方法,以减小计算量,提高计算速度。当图像序列对应的组织旋转比较明显的时候,相隔帧数越大,配准图像之间的相对旋转角度越大,配准的准确度越低。本发明中,在完成一次配准之前,首先根据模板配准结果的标准差来判断图像之间的旋转角度是否过大,假如角度较小则该次配准结果有效,继续下一次配准;假如角度过大,则降低相隔帧数,重新进行配准。具体为:设配准失败的两帧图像序号分别为m和n,则要重新进行两次配准,配准的图像帧分别为m和(m+n)/2、(m+n)/2和n。假如这两次配准依旧失败,则继续折半配准,直到配准成功。
下面的实验是采用本发明所述的方法前后的宽景成像效果,用于验证本发明方法的有效性。
为了测量本发明方法的成像准确度,我们做了一个实验装置。如图4所示,在一个塑料盒子中固定三个塑料棒子。塑料棒子的直径为10mm,且棒子之间的相互距离是是已知的(分别为112.5mm、111.0mm、200.0mm,上面两个塑料棒子之间夹角为127.0°),可以作为距离和角度的测量标准,来测量宽景成像的准确性。实验中,三个塑料棒子嵌入到一块猪肉中,用超声探头在猪肉上表面采集到一系列图像数据,用于宽景成像。
图5、图6分别是没有采用和采用本发明进行改进前后的宽景图像。从图5中可见,改进前的宽景图中有两个比较明显的误差区域。图6是改进后的宽景图,基本上跟图4装置图中的猪肉形状相符。
下表是宽景图中各个塑料棒子距离的测量距离和角度及其误差率。
从表中可见,用本发明方法改进后的宽景图中,除了1-2棒的距离准确度稍微降低一点,其他的距离的准确度都有所提高;特别是1-2-3棒子的角度,误差率从14.25%提高到7.32%,改进效果明显。证明了本发明方法对宽景成像中旋转角度所造成的误差的纠正效果。
如上所述便可较好地实现本发明。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。
Claims (4)
1.一种医学超声宽景成像中的旋转判断及误差纠正方法,其特征在于包括如下步骤:
(1)、统计模板配准结果的标准差;
(2)、根据模板配准结果标准差来判断配准图像之间旋转角度的大小;
(3)、根据旋转角度的大小,相应的修正图像序列配准的间隔帧数。
2.根据权利要求1所述医学超声宽景成像中的旋转判断及误差纠正方法,其特征在于所述步骤(1)采用多模板的方式配准两帧超声图像,并统计模板配准结果的标准差,具体操作为:在待配准图像中选择多个模板,在水平和垂直方向搜索每个模板在固定图像中的匹配位置,得到每个模板的配准结果:Rb(i),i=1...Nb,Nb是模板的个数,然后计算所有模板配准结果的标准差(SDr):
3.根据权利要求1或2所述的医学超声宽景成像中的旋转判断及误差纠正方法,其特征在于所述步骤(2)根据模板配准结果标准差来判断配准图像之间旋转角度的大小,当模板配准结果标准差大于阈值TSD的时候,就确定配准图像之间的旋转角度过大。
4.根据权利要求3所述医学超声宽景成像中的旋转判断及误差纠正方法,其特征在于所述步骤(3)根据旋转角度的大小,相应的修正图像序列的间隔帧数,具体如下:
在配准整个图像序列的时候采用的相隔多帧配准一次的方法,在完成一次配准之前,首先根据模板配准结果的标准差来判断图像之间的旋转角度是否过大,若角度较小则该次配准结果有效,继续下一次配准;若角度过大,则降低相隔帧数,重新进行配准,具体为:设配准失败的两帧图像序号分别为m和n,则要重新进行两次配准,配准的图像帧分别为m和(m+n)/2、(m+n)/2和n,若这两次配准依旧失败,则继续折半配准,直到配准成功。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110139790 CN102283675B (zh) | 2011-05-27 | 2011-05-27 | 一种医学超声宽景成像中的旋转判断及误差纠正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110139790 CN102283675B (zh) | 2011-05-27 | 2011-05-27 | 一种医学超声宽景成像中的旋转判断及误差纠正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102283675A true CN102283675A (zh) | 2011-12-21 |
CN102283675B CN102283675B (zh) | 2013-04-17 |
Family
ID=45330630
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110139790 Expired - Fee Related CN102283675B (zh) | 2011-05-27 | 2011-05-27 | 一种医学超声宽景成像中的旋转判断及误差纠正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102283675B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106999733A (zh) * | 2014-12-10 | 2017-08-01 | 因赛泰克有限公司 | 用于优化经颅声治疗的***和方法 |
CN107106144A (zh) * | 2015-02-26 | 2017-08-29 | 株式会社日立制作所 | 超声波成像装置以及图像处理装置 |
CN108742549A (zh) * | 2018-06-26 | 2018-11-06 | 京东方科技集团股份有限公司 | 一种图像信息生成方法和脉搏波测量*** |
CN111310851A (zh) * | 2020-03-03 | 2020-06-19 | 四川大学华西第二医院 | 一种人工智能超声辅助***及其应用 |
CN111968222A (zh) * | 2020-08-04 | 2020-11-20 | 西北工业大学 | 一种人体组织非静止状态下三维超声重建方法 |
CN114373004A (zh) * | 2022-01-13 | 2022-04-19 | 强联智创(北京)科技有限公司 | 一种基于动态级联网络的无监督三维图像刚性配准方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6159152A (en) * | 1998-10-26 | 2000-12-12 | Acuson Corporation | Medical diagnostic ultrasound system and method for multiple image registration |
US20030036701A1 (en) * | 2001-08-10 | 2003-02-20 | Dong Fang F. | Method and apparatus for rotation registration of extended field of view ultrasound images |
CN1839760A (zh) * | 2005-03-29 | 2006-10-04 | 深圳迈瑞生物医疗电子股份有限公司 | 宽景成像处理方法 |
CN101455576A (zh) * | 2007-12-12 | 2009-06-17 | 深圳迈瑞生物医疗电子股份有限公司 | 超声宽景成像的方法、装置与*** |
CN101612048A (zh) * | 2008-06-27 | 2009-12-30 | 西门子公司 | 一种超声波探头和超声波扫描装置 |
CN101901481A (zh) * | 2010-08-11 | 2010-12-01 | 深圳市蓝韵实业有限公司 | 一种图像拼接方法 |
-
2011
- 2011-05-27 CN CN 201110139790 patent/CN102283675B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6159152A (en) * | 1998-10-26 | 2000-12-12 | Acuson Corporation | Medical diagnostic ultrasound system and method for multiple image registration |
US20030036701A1 (en) * | 2001-08-10 | 2003-02-20 | Dong Fang F. | Method and apparatus for rotation registration of extended field of view ultrasound images |
CN1839760A (zh) * | 2005-03-29 | 2006-10-04 | 深圳迈瑞生物医疗电子股份有限公司 | 宽景成像处理方法 |
CN101455576A (zh) * | 2007-12-12 | 2009-06-17 | 深圳迈瑞生物医疗电子股份有限公司 | 超声宽景成像的方法、装置与*** |
CN101612048A (zh) * | 2008-06-27 | 2009-12-30 | 西门子公司 | 一种超声波探头和超声波扫描装置 |
CN101901481A (zh) * | 2010-08-11 | 2010-12-01 | 深圳市蓝韵实业有限公司 | 一种图像拼接方法 |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106999733A (zh) * | 2014-12-10 | 2017-08-01 | 因赛泰克有限公司 | 用于优化经颅声治疗的***和方法 |
CN106999733B (zh) * | 2014-12-10 | 2019-08-30 | 因赛泰克有限公司 | 用于优化经颅声治疗的***和方法 |
CN107106144A (zh) * | 2015-02-26 | 2017-08-29 | 株式会社日立制作所 | 超声波成像装置以及图像处理装置 |
CN107106144B (zh) * | 2015-02-26 | 2020-08-11 | 株式会社日立制作所 | 超声波成像装置以及图像处理装置 |
CN108742549A (zh) * | 2018-06-26 | 2018-11-06 | 京东方科技集团股份有限公司 | 一种图像信息生成方法和脉搏波测量*** |
US11497408B2 (en) | 2018-06-26 | 2022-11-15 | Beijing Boe Technology Development Co., Ltd. | Image information generation method, pulse wave measurement system and electronic device |
CN111310851A (zh) * | 2020-03-03 | 2020-06-19 | 四川大学华西第二医院 | 一种人工智能超声辅助***及其应用 |
CN111968222A (zh) * | 2020-08-04 | 2020-11-20 | 西北工业大学 | 一种人体组织非静止状态下三维超声重建方法 |
CN114373004A (zh) * | 2022-01-13 | 2022-04-19 | 强联智创(北京)科技有限公司 | 一种基于动态级联网络的无监督三维图像刚性配准方法 |
CN114373004B (zh) * | 2022-01-13 | 2024-04-02 | 强联智创(北京)科技有限公司 | 一种动态图像配准方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102283675B (zh) | 2013-04-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102283675B (zh) | 一种医学超声宽景成像中的旋转判断及误差纠正方法 | |
CN107665497B (zh) | 一种医学图像中计算心胸比的方法 | |
US8265363B2 (en) | Method and apparatus for automatically identifying image views in a 3D dataset | |
CN100455266C (zh) | 宽景成像处理方法 | |
CN107123137B (zh) | 医学图像处理方法及设备 | |
US8036436B2 (en) | System and method for segmenting a region in a medical image | |
US9607405B2 (en) | Method and device for detecting displacement in elastography | |
CN105809702B (zh) | 一种基于Tsai算法的改进位姿估计方法 | |
CN100561518C (zh) | 基于感兴趣区域的自适应医学序列图像插值方法 | |
CN107481228A (zh) | 基于计算机视觉的人体背部脊柱侧弯角度测量方法 | |
WO2017036231A1 (zh) | 从光学相干断层图像获取视网膜结构的方法及其*** | |
JPH08206117A (ja) | 超音波診断装置 | |
CN101909165B (zh) | 一种基于混合测度的视频数据宽景成像方法 | |
CN101996416B (zh) | 3d人脸捕获方法和设备 | |
CN106934821A (zh) | 一种基于icp算法和b样条的锥形束ct和ct图像配准方法 | |
CN108682025A (zh) | 一种图像配准方法及装置 | |
CN110458874B (zh) | 一种图像非刚性配准方法和*** | |
WO2005046473A1 (en) | A method for processing slice images | |
CN103295209A (zh) | Dr图像的拼接方法及*** | |
CN109242947A (zh) | 三维超声图像显示方法 | |
CN103544695A (zh) | 一种高效的基于博弈框架的医学图像分割方法 | |
CN103295262A (zh) | 管状腔体组织的旋转多角度曲面重建方法及装置 | |
CN112329699A (zh) | 一种像素级精度的人眼注视点定位方法 | |
CN106846338A (zh) | 基于混合模型的视网膜oct图像视***结构分割方法 | |
CN110599501B (zh) | 一种胃肠结构的真实尺度三维重建及可视化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130417 Termination date: 20200527 |