《基于GPU实现的超声成像扫描变换方法.pdf》由会员分享,可在线阅读,更多相关《基于GPU实现的超声成像扫描变换方法.pdf(5页珍藏版)》请在专利查询网上搜索。
1、(10)授权公告号 CN 102631222 B (45)授权公告日 2014.01.01 CN 102631222 B (21)申请号 201210127194.2 (22)申请日 2012.04.26 A61B 8/00(2006.01) (73)专利权人 珠海医凯电子科技有限公司 地址 519041 广东省珠海市金湾区机场东路 288 号医疗器械产业园 D 栋 4 楼 (72)发明人 丁波 朱逸斐 王建和 (74)专利代理机构 广州嘉权专利商标事务所有 限公司 44205 代理人 谭志强 US 2008/0161690 A1,2008.07.03, US 2005/0122333 A1,。
2、2005.06.09, EP 1424721 A2,2004.06.02, 全文 . CN 101606854 A,2009.12.23, CN 101543412 A,2009.09.30, CN 101396281 A,2009.04.01, CN 101125089 A,2008.02.20, 谈宜育等 . 一种 B 超数字扫描变换器 (DSC) 的 FPGA 实现 .中国医疗器械杂志 .2001, 第 25 卷 ( 第 4 期 ), 第 213-214 页 . 彭龙飞等 . 一种高效实现 B 超坐标变换的方 法 .声学技术 .2009, 第 28 卷 ( 第 4 期 ), 第 260-。
3、263 页 . (54) 发明名称 基于 GPU 实现的超声成像扫描变换方法 (57) 摘要 本发明公开了基于 GPU 实现的超声成像扫描 变换方法, 其获取扫描数据, 将扫描数据划分形成 n*m 份的矩形, 将获得的矩形分成三角形, 完成坐 标变换, 再对每个三角形依次进行绘制, 得到扫描 变换过后的图像。本发明基于 GPU 在不增加系统 成本和复杂度的前提下, 完成扫描变换, 提供更好 的实时性和图像质量, 本发明的方法可以使用最 廉价的硬件完成, 不但降低成本, 而且硬件结构简 单, 更适合小型化, 同时由于硬件结构简单, 整体 功耗也会比以往的系统小很多, 符合绿色环保的 理念。 (5。
4、1)Int.Cl. (56)对比文件 审查员 许流芳 权利要求书 1 页 说明书 2 页 附图 1 页 (19)中华人民共和国国家知识产权局 (12)发明专利 权利要求书1页 说明书2页 附图1页 (10)授权公告号 CN 102631222 B CN 102631222 B 1/1 页 2 1. 基于 GPU 实现的超声成像扫描变换方法, 其特征在于包括以下步骤 : 获取扫描数据 s(, d), 为探头扫描方向, d 为深度方向 ; 将扫描数据 s(, d) 在 方向分为 n 等份, 在 d 方向分为 m 等份, 形成 n*m 份的矩 形, n 和 m 的值为不小于 10 的整数 ; 将获得。
5、的矩形分成三角形, 矩形划分为三角形时按照矩形对角线划分成两个三角形, 计算与s(, d)平面的三角形对应的在xy平面上的三角形顶点坐标, 坐标变换关系为 : x= (d0+d)sin(), y= (d0+d)cos() ; 把 s(, d) 作为 GPU 的纹理, 方向为纹理的 u 坐标, d 方向为纹理的 v 坐标, xy 平 面三角形的坐标作为 GPU 的顶点坐标, xy 平面三角形对应的 u 坐标和 v 坐标为纹理坐标, 对每个三角形依次进行绘制, 得到扫描变换过后的图像, 获取扫描数据 s(, d) 后先创建 好对应大小的 texture, 扫描数据 s(, d) 划分为三角形后计算。
6、好每个三角形的顶点和纹 理坐标, 顶点和纹理坐标存在内存中, 绘制时打开 GPU 的双线性插值开关, 将接收到的扫描 数据上传到预先创建好的 texture, 绑定该 texture, 用预先计算好的顶点和纹理坐标逐个 绘制三角形。 权 利 要 求 书 CN 102631222 B 2 1/2 页 3 基于 GPU 实现的超声成像扫描变换方法 技术领域 0001 本发明涉及医用超声成像技术, 特别是基于 GPU(Graphics Process Unit) 图形 处理技术实现的超声成像扫描变换方法。 背景技术 0002 大部分的腹部 B 超探头的换能器是扇形的, 所以从换能器所得到的扫描数据是。
7、 极坐标系的。在屏幕上显示的时候我们需要把极坐标转换成直角坐标才能显示。扫描变 换本质上是一个从极坐标系变换到直角坐标系的过程。如图 1 和 2 所示, 探头扫描方向为 , 深度方向为 d, 扫描所得的数据是 s(, d), 通过坐标变换, 将数据变换到 xy 平面下的 i(x,y), 具体变换公式为 : 0003 x= (d0+d)sin() 0004 y= (d0+d)cos() 0005 坐标变换是一个非常耗时的操作。 为了能够达到实时性的要求, 即使采用查表法, 对 CPU 的负担也是很重的。通常需要非常高性能的 CPU, 不采用双线性插值, 或者采用 DSP 或者 FPGA 辅助。 。
8、0006 但是这些方法都有缺点 : 1) 高性能的 CPU 通常成本很高, 而且无法使设备小型化 ; 2) 不采用双线性插值能够省去一部分计算量, 但是也使图像质量有明显下降 ; 3) 采用 DSP 或者 FPGA 不但增加成本, 同时也增加了系统的复杂度。 发明内容 0007 针对上述问题, 本发明提供一种实现快速、 低成本、 高性能的基于 GPU 实现的超声 成像扫描变换方法。 0008 本发明为解决其技术问题所采用的技术方案是 : 0009 基于 GPU 实现的超声成像扫描变换方法, 包括以下步骤 : 获取扫描数据 s(, d) ; 将扫描数据 s(, d) 在 方向分为 n 等份, 在。
9、 d 方向分为 m 等份, 形成 n*m 份的矩形 ; 将 获得的矩形分成三角形, 计算与 s(, d) 平面的三角形对应的在 xy 平面上的三角形顶点 坐标, 坐标变换关系为 : x= (d0+d)sin(), y= (d0+d)cos() ; 把 s(, d) 作为 GPU 的纹 理, 方向为纹理的 u 坐标, d 方向为纹理的 v 坐标, xy 平面三角形的坐标作为 GPU 的顶点 坐标, xy 平面三角形对应的 u 坐标和 v 坐标为纹理坐标, 对每个三角形依次进行绘制, 得到 扫描变换过后的图像。 0010 进一步, n 和 m 的值为不小于 10 的整数。 0011 进一步, 矩形。
10、划分为三角形时按照矩形对角线划分成两个三角形。 0012 进一步, 获取扫描数据s(, d)后先创建好对应大小的texture, 扫描数据s(, d) 划分为三角形后计算好每个三角形的顶点和纹理坐标, 顶点和纹理坐标存在内存中, 绘 制时打开 GPU 的双线性插值开关, 将接收到的扫描数据上传到预先创建好的 texture, 绑定 该 texture, 用预先计算好的顶点和纹理坐标逐个绘制三角形。 说 明 书 CN 102631222 B 3 2/2 页 4 0013 本发明的有益效果是 : 本发明基于 GPU 在不增加系统成本和复杂度的前提下, 完 成扫描变换, 提供更好的实时性和图像质量,。
11、 本发明的方法可以使用最廉价的硬件完成, 不 但降低成本, 而且硬件结构简单, 更适合小型化, 同时由于硬件结构简单, 整体功耗也会比 以往的系统小很多, 符合绿色环保的理念。 附图说明 0014 下面结合附图和具体实施方式进行进一步的说明 : 0015 图 1 为变换前的极坐标系扫描数据图 ; 0016 图 2 为变换后的直角坐标系扫描数据图 ; 0017 图 3 为变换前的极坐标系划分成三角形的扫描数据图 ; 0018 图 4 为变换后的直角坐标系划分成三角形的扫描数据图。 具体实施方式 0019 参照图 1、 2、 3 和 4, 本发明的基于 GPU 实现的超声成像扫描变换方法, 用于将。
12、从换 能器所得到的极坐标系下的扫描数据转换成可在屏幕上显示的直角坐标系数据。 下面对本 发明的详细步骤进行描述。 0020 获取扫描数据s(, d)的步骤, 如图1所示, 由于此数据由扇形换能器生成的, 因 此此获取到的扫描数据s(, d)是极坐标系下的, 这些数据必须变换成直角坐标系下才能 显示。 0021 将扫描数据划分的步骤 : 将扫描数据 s(, d) 在 方向分为 n 等份, 在 d 方向 分为 m 等份, 形成 n*m 份的矩形。一般地, 划分的数量需要足够大, 才能形成 n*m 份的矩形, 在一般的实施方式中, n、 m 的数据优选的不小于 10 的整数。 0022 划分三角形与。
13、计算三角形顶点坐标的步骤 : 将获得的矩形分成三角形, 计算 与 s(, d) 平面的三角形对应的在 xy 平面上的三角形顶点坐标, 坐标变换关系为 : x= (d0+d)sin(), y= (d0+d)cos()。4) 当 m,n 足够大时, 我们可以认为在三角形内, 每个 点的坐标变换是线性分布的。之所以划分为三角形, 因为 GPU 一般只能对三角形进行绘制。 作为优选的实施方式, 矩形划分为三角形时按照矩形对角线划分成两个三角形。划分成的 示意图如图 3 所示。 0023 绘制的步骤 : 把 s(, d) 作为 GPU 的纹理, 方向为纹理的 u 坐标, d 方向为纹 理的 v 坐标, 。
14、xy 平面三角形的坐标作为 GPU 的顶点坐标, xy 平面三角形对应的 u 坐标和 v 坐标为纹理坐标, 对每个三角形依次进行绘制, 得到扫描变换过后的图像, 如图 4 所示。 0024 在具体实施过程中, 可选用 Direct3D 作为 GPU 加速的 API。当选用 Direct3D 作 为 GPU 加速的 API 来实现本方法时, 一般采用以下方式进行 : 获取扫描数据 s(, d) 后先 创建好对应大小的 texture, 扫描数据 s(, d) 划分为三角形后计算好每个三角形的顶点 和纹理坐标, 顶点和纹理坐标存在内存中, 绘制时打开 GPU 的双线性插值开关, 将接收到的 扫描数据上传到预先创建好的 texture, 绑定该 texture, 用预先计算好的顶点和纹理坐标 逐个绘制三角形。 0025 本发明的实施方式并不受上述实施例的限制, 只要其以基本相同的手段达到本发 明的技术效果, 都应属于本发明的保护范围。 说 明 书 CN 102631222 B 4 1/1 页 5 图 1 图 2 图 3 图 4 说 明 书 附 图 CN 102631222 B 5 。