CN102038519A - 超声诊断的空间复合成像方法 - Google Patents
超声诊断的空间复合成像方法 Download PDFInfo
- Publication number
- CN102038519A CN102038519A CN 201010616777 CN201010616777A CN102038519A CN 102038519 A CN102038519 A CN 102038519A CN 201010616777 CN201010616777 CN 201010616777 CN 201010616777 A CN201010616777 A CN 201010616777A CN 102038519 A CN102038519 A CN 102038519A
- Authority
- CN
- China
- Prior art keywords
- window function
- image data
- different directions
- diagnosis
- ultrasound image
- 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)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了超声诊断的空间复合成像方法,其中对N个不同方向超声扫查中发射和接受通道的波束信号分别通过窗函数处理,然后对获得的N个不同方向的原始超声图像数据通过直方图均衡算法或匹配算法进行处理,再对处理后的超声图像数据进行空间复合成像,得到最终的超声图像,本领域的技术人员需要作一些调试,选用合适的窗函数处理,则加强了目标组织中诊断部位的波束信号强度,合适的直方图均衡算法或直方图匹配算法则进一步提升诊断部位的图像数据,增强最终的超声图像中诊断部位的信息,这样最终的超声图像的对比分辨率和组织结构的层次感较好,显示了诊断部位的更多有用信息。
Description
技术领域
本发明涉及超声诊断的空间复合成像方法。
背景技术
在超声诊断仪中,现有的空间复合成像方法,包括:对目标组织进行多个不同方向的超声扫查,在图1中,以向左斜扫得到超声图像A,以正扫得到超声图像B,以向右斜扫得到超声图像C,然后根据超声图像B的空间坐标对相同坐标的三幅图像点幅值进行线性平均,得到新的像点组成的超声图像D,超声图像D的纹理结构细腻,空间组织结构丰富,斑点噪声较少,但是这种对相同空间不同时间帧和不同方向扫查的图像做线性平均的方法,会降低超声图像D的对比分辨率和组织结构的层次感,使得超声图像D中和诊断部位有关的有用信息不多。
发明内容
本发明的目的是提供先对不同方向的超声发射和接收波束信号进行窗函数处理再对不同方向的超声图像数据进行直方图均衡或匹配算法处理的空间复合成像方法,该方法可以提高超声图像的对比分辨率和组织结构的层次感,显示诊断部位的更多信息。
为达到上述目的,本发明采用的技术方案是:超声诊断的空间复合成像方法,包括:
对目标组织进行N个不同方向的超声扫查,N>=2,对所述N个不同方向的超声发射通道和接受通道中的波束信号分别进行窗函数处理,然后获得所述N个不同方向的超声图像数据D’i,1=<i<=N,接着通过直方图均衡算法或直方图匹配算法对所述超声图像数据D’i进行处理,获得处理后的超声图像数据Di,最后对N个不同方向的所述超声图像数据Di中相同位置的图像数据取平均,得到最终的超声图像,
所述的窗函数处理用于加强所述目标组织中诊断部位的波束信号强度,所述直方图均衡算法或直方图匹配算法用于进一步提升所述诊断部位的图像数据,增强所述最终的超声图像中所述诊断部位的信息。
所述窗函数采用高斯窗函数、汉宁窗函数或者海明窗函数,基于所述加强诊断部位的波束信号强度的目的,对不同方向的波束信号采用合适的窗函数进行信号处理。
由于上述技术方案的运用,本发明与现有技术相比具有下列优点:对N个不同方向超声扫查中发射和接受通道的波束信号分别通过窗函数处理,然后对获得的N个不同方向的原始超声图像数据通过直方图均衡算法或匹配算法进行处理,再对处理后的超声图像数据进行空间复合成像,得到最终的超声图像,本领域的技术人员需要作一些调试,选用合适的窗函数处理,则加强了目标组织中诊断部位的波束信号强度,合适的直方图均衡算法或直方图匹配算法则进一步提升诊断部位的图像数据,增强最终的超声图像中诊断部位的信息,这样最终的超声图像的对比分辨率和组织结构的层次感较好,显示了诊断部位的更多有用信息。
附图说明
附图1为解释现有的空间复合成像方法的示意图;
附图2为解释本发明中空间复合成像方法的示意图;
附图3为本发明中窗函数处理的示意图。
具体实施方式
下面结合附图来进一步阐述本发明的结构。
在图2-3中,超声诊断的空间复合成像方法,通过对多个不同方向扫查得到的超声图像进行空间复合叠加从而得到高质量的超声图像。在对目标组织的不同方向扫查时,由于扫查角度的不同和扫查时间的略微差异,则不同扫查方向得到的实际超声回波信号是不相同,这就需要我们对相同空间不同角度和时间的图像帧进行不同的信号和图像处理,才能获得最佳的空间复合成像。具体来说,先通过对N(N>=2)个不同方向超声扫查中发射和接受通道的波束信号分别进行变迹窗控制,然后对N个不同方向的原始超声图像数据进行自适应的灰阶映射,再对映射后的超声图像数据进行空间复合叠加成像,如图2所示,这样就能得到对比分辨率和组织结构的层次感都较好的超声图像。
在图3中,通过三个方向向左斜扫、正扫、向右斜扫分别扫查目标组织,在向左斜扫时,选用合适的窗函数,对向左斜扫时发射和接收通道中的波束信号进行处理,得到向左斜扫的超声图像数据D’1,同样,在正扫和向右斜扫时,也需要选用合适的窗函数,来对正扫和向右斜扫时发射和接收通道中的波束信号进行处理,得到正扫和向右斜扫的超声图像数据D’2和D’3,窗函数处理用于控制波束的宽度,目的在于加强目标组织中诊断部位的波束信号强度,窗函数可以选用高斯函数或者汉宁窗函数(α=0.5),或者海明窗函数(α=0.54),x为目标组织的横向位置。窗函数处理时,要根据具体的扫查情况,本领域的技术人员进行一些调试,才能选用到合适的窗函数以及合适的权重系数,合适的窗函数处理才能达到加强有用的波束信号强度的目的。譬如,如果选用高斯函数作为窗函数,可以通过扫查超声仿体模块,观察三个方向向左斜扫、正扫、向右斜扫所得的超声图像中的固定深度的靶点形状(一般以圆形为最佳)和亮度,选取对左斜扫、正扫、向右斜的超声图像数据的最佳的σ值来调试得到与之对应的窗函数,达到加强有用的波束信号强度的目的。
在得到不同方向(向左斜扫、正扫、向右斜扫)的超声图像数据D’1、D’2和D’3后,根据直方图分析分别进行自适应的灰阶映射,即采用直方图均衡算法或者直方图匹配算法对超声图像数据D’1、D’2和D’3分别进行处理,得到处理后的超声图像数据D1、D2和D3。
直方图均衡算法是把原始图像的灰度直方图从比较集中的某个灰度区间变成在全部灰度范围内的均匀分布,从而增加象素灰度值的动态范围,达到增强图像整体对比度的效果。具体的算法描述可参见百度检索http://baike.baidu.com/view/1164383.htm。在不同方向扫查得到的超声图像数据D’1、D’2和D’3中,图像的灰度集中区间不相同,采用直方图均衡算法可分别将超声图像数据D’1、D’2和D’3的灰度集中区间提取出来,即提取出目标组织中和诊断部位有关的有用数据,这样可以将诊断部位的信息增强。
直方图匹配算法,就是将图像直方图以标准图像的直方图为标准作变换,使两图像的直方图相同和近似,从而使两幅图像具有类似的色调和反差。譬如,假如超声图像D2包含了较多的用于显示诊断部位的有用信息,则可以把超声图像D2作为标准图像,然后把超声图像D’1和超声图像D’3的直方图以超声图像D2的直方图为标准作变换,从而使得处理后的超声图像D1和超声图像D3与超声图像D2具有类似的图像色调、反差以及灰阶分布,达到较好的空间复合成像效果,具体的变换方法可以参见百度检索“直方图匹配算法”。直方图匹配算法同样提升了有用图像数据,增强显示诊断部位的信息。
在进行图像数据处理时,可以将直方图均衡算法和匹配算法相结合,譬如说,超声图像D’2通过直方图均衡算法处理,得到新的图像D2,然后以新的图像D2作为标准图像,对超声图像D’1和超声图像D’3进行直方图匹配算法处理。
最后对多个不同方向的超声图像数据D1、D2和D3中相同位置的图像数据取平均,得到最终的超声图像D,该超声图像D的对比分辨率和组织结构的层次感较好,显示了诊断部位的更多有用信息。
Claims (2)
1.超声诊断的空间复合成像方法,其特征在于:包括:
对目标组织进行N个不同方向的超声扫查,N>=2,对所述N个不同方向的超声发射通道和接受通道中的波束信号分别进行窗函数处理,然后获得所述N个不同方向的超声图像数据D’i,1=<i<=N,接着通过直方图均衡算法或直方图匹配算法对所述超声图像数据D’i进行处理,获得处理后的超声图像数据Di,最后对N个不同方向的所述超声图像数据Di中相同位置的图像数据取平均,得到最终的超声图像,
所述的窗函数处理用于加强所述目标组织中诊断部位的波束信号强度,所述直方图均衡算法或直方图匹配算法用于进一步提升所述诊断部位的图像数据,增强所述最终的超声图像中所述诊断部位的信息。
2.根据权利要求1所述的空间复合成像方法,其特征在于:所述窗函数采用高斯窗函数、汉宁窗函数或者海明窗函数,基于所述加强诊断部位的波束信号强度的目的,对不同方向的波束信号采用合适的窗函数进行信号处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010106167772A CN102038519B (zh) | 2010-12-31 | 2010-12-31 | 超声诊断的空间复合成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010106167772A CN102038519B (zh) | 2010-12-31 | 2010-12-31 | 超声诊断的空间复合成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102038519A true CN102038519A (zh) | 2011-05-04 |
CN102038519B CN102038519B (zh) | 2012-11-28 |
Family
ID=43905304
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010106167772A Expired - Fee Related CN102038519B (zh) | 2010-12-31 | 2010-12-31 | 超声诊断的空间复合成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102038519B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102551808A (zh) * | 2012-02-29 | 2012-07-11 | 飞依诺科技(苏州)有限公司 | 超声诊断中图形扩展的成像方法 |
CN102727255A (zh) * | 2012-07-13 | 2012-10-17 | 深圳市理邦精密仪器股份有限公司 | 一种超声图像空间复合成像的方法及装置 |
CN106056567A (zh) * | 2016-07-22 | 2016-10-26 | 珠海医凯电子科技有限公司 | 一种超声波空间复合成像方法 |
CN106373103A (zh) * | 2016-09-08 | 2017-02-01 | 飞依诺科技(苏州)有限公司 | 一种超声数据的复合方法及装置 |
CN106780349A (zh) * | 2016-11-21 | 2017-05-31 | 深圳开立生物医疗科技股份有限公司 | 基于fpga的空间复合成像方法以及装置 |
CN109549669A (zh) * | 2018-12-05 | 2019-04-02 | 余姚市华耀工具科技有限公司 | 定向检测式b超成像仪 |
CN111839588A (zh) * | 2020-07-29 | 2020-10-30 | 深圳开立生物医疗科技股份有限公司 | 一种超声图像生成方法、装置、设备及存储介质 |
TWI717589B (zh) * | 2016-04-05 | 2021-02-01 | 日商島津製作所股份有限公司 | 診斷畫像系統 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6126599A (en) * | 1998-10-01 | 2000-10-03 | Atl Ultrasound, Inc. | Ultrasonic diagnostic imaging system with real time spatial compounding processor |
CN101199430A (zh) * | 2006-12-15 | 2008-06-18 | 深圳迈瑞生物医疗电子股份有限公司 | 空间复合成像方法、设备及其超声成像系统 |
CN101617946A (zh) * | 2008-07-04 | 2010-01-06 | 深圳迈瑞生物医疗电子股份有限公司 | 用于超声成像的空间复合方法与装置以及超声成像系统 |
-
2010
- 2010-12-31 CN CN2010106167772A patent/CN102038519B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6126599A (en) * | 1998-10-01 | 2000-10-03 | Atl Ultrasound, Inc. | Ultrasonic diagnostic imaging system with real time spatial compounding processor |
CN101199430A (zh) * | 2006-12-15 | 2008-06-18 | 深圳迈瑞生物医疗电子股份有限公司 | 空间复合成像方法、设备及其超声成像系统 |
CN101617946A (zh) * | 2008-07-04 | 2010-01-06 | 深圳迈瑞生物医疗电子股份有限公司 | 用于超声成像的空间复合方法与装置以及超声成像系统 |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102551808A (zh) * | 2012-02-29 | 2012-07-11 | 飞依诺科技(苏州)有限公司 | 超声诊断中图形扩展的成像方法 |
CN102727255A (zh) * | 2012-07-13 | 2012-10-17 | 深圳市理邦精密仪器股份有限公司 | 一种超声图像空间复合成像的方法及装置 |
CN102727255B (zh) * | 2012-07-13 | 2014-04-30 | 深圳市理邦精密仪器股份有限公司 | 一种超声图像空间复合成像的方法及装置 |
TWI717589B (zh) * | 2016-04-05 | 2021-02-01 | 日商島津製作所股份有限公司 | 診斷畫像系統 |
CN106056567A (zh) * | 2016-07-22 | 2016-10-26 | 珠海医凯电子科技有限公司 | 一种超声波空间复合成像方法 |
CN106373103A (zh) * | 2016-09-08 | 2017-02-01 | 飞依诺科技(苏州)有限公司 | 一种超声数据的复合方法及装置 |
CN106780349A (zh) * | 2016-11-21 | 2017-05-31 | 深圳开立生物医疗科技股份有限公司 | 基于fpga的空间复合成像方法以及装置 |
WO2018090654A1 (zh) * | 2016-11-21 | 2018-05-24 | 深圳开立生物医疗科技股份有限公司 | 基于fpga的空间复合成像方法以及装置 |
CN106780349B (zh) * | 2016-11-21 | 2020-02-14 | 深圳开立生物医疗科技股份有限公司 | 基于fpga的空间复合成像方法以及装置 |
CN109549669A (zh) * | 2018-12-05 | 2019-04-02 | 余姚市华耀工具科技有限公司 | 定向检测式b超成像仪 |
CN111839588A (zh) * | 2020-07-29 | 2020-10-30 | 深圳开立生物医疗科技股份有限公司 | 一种超声图像生成方法、装置、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN102038519B (zh) | 2012-11-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102038519B (zh) | 超声诊断的空间复合成像方法 | |
US9008383B2 (en) | Enhancing quality of ultrasound image in ultrasound system | |
US20210132223A1 (en) | Method and Apparatus for Ultrasound Imaging with Improved Beamforming | |
CN101683274B (zh) | 偏转凸阵数字扫描变换和凸阵空间复合的方法及成像系统 | |
CN105074498B (zh) | 具有对梯形的扇区的空间复合的超声诊断成像系统 | |
CN103181779A (zh) | 超声波能量转换装置以及超声波成像系统与方法 | |
CN101617946B (zh) | 用于超声成像的空间复合方法与装置以及超声成像系统 | |
CN1553785A (zh) | 超声波摄像装置及其摄像数据处理方法以及处理程序 | |
CN104013424B (zh) | 一种基于深度信息的超声宽景成像方法 | |
US8152722B2 (en) | Method and apparatus for quickly determining an imaging region in an ultrasonic imaging system | |
CN105496515A (zh) | 穿刺增强方法及系统 | |
US8724880B2 (en) | Ultrasonic diagnostic apparatus and medical image processing apparatus | |
CN104490418A (zh) | 一种基于信号统计分析的超声图像自动优化方法 | |
CN102539532A (zh) | 一种基于二维邻域合成孔径聚焦的超声c扫描成像方法 | |
CN108354627A (zh) | 一种提高帧频的超声波束形成方法 | |
CN110477947A (zh) | 基于深度学习的平面波波束合成方法、系统、存储介质、设备 | |
CN108022209A (zh) | 颜色值的获取方法及装置、空间复合成像的方法及设备 | |
CN107997783A (zh) | 一种基于超声波方向性的自适应超声波束合成方法和系统 | |
US9081084B2 (en) | Ultrasound system and method for processing beam-forming based on sampling data | |
JP2000149015A (ja) | 画像のエッジ強調方法及びイメ―ジング装置 | |
CN116990824A (zh) | 集群侧扫系统的图形地理信息编码和融合方法 | |
AU2013273647B2 (en) | Improvements in and relating to ophthalmoscopes | |
US20110245676A1 (en) | Method and apparatus for ultrasound signal acquisition and processing | |
CN101390756B (zh) | 超声成像系统中的扫描线插线处理方法 | |
CN102423265B (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 | ||
C53 | Correction of patent for invention or patent application | ||
CB03 | Change of inventor or designer information |
Inventor after: Tian Yuan Inventor before: Tian Lin |
|
COR | Change of bibliographic data |
Free format text: CORRECT: INVENTOR; FROM: TIAN LIN TO: TIAN YUAN |
|
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: 20121128 Termination date: 20201231 |