[发明专利]一种PROPELLER技术中抑制非刚性运动的权重计算方法有效
申请号: | 201810155551.3 | 申请日: | 2018-02-23 |
公开(公告)号: | CN108577841B | 公开(公告)日: | 2021-09-10 |
发明(设计)人: | 罗海;朱高杰;周翔;王超;陈梅泞;刘霞;卞威;吴子岳 | 申请(专利权)人: | 奥泰医疗系统有限责任公司 |
主分类号: | A61B5/055 | 分类号: | A61B5/055 |
代理公司: | 成都高远知识产权代理事务所(普通合伙) 51222 | 代理人: | 谢一平;曾克 |
地址: | 610000 四川省*** | 国省代码: | 四川;51 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 propeller 技术 抑制 刚性 运动 权重 计算方法 | ||
1.一种PROPELLER技术中抑制非刚性运动的权重计算方法,其特征在于,包括:
采集一个以上不同角度的K空间数据带;
计算所述一个以上不同角度的K空间数据带中每两个K空间数据带之间的相关性,获取相关性矩阵;
根据所述相关性矩阵计算出每个所述K空间数据带的权重;
所述计算所述一个以上不同角度的K空间数据带中每两个K空间数据带之间的相关性,获取相关性矩阵的方法包括:
对所述一个以上不同角度的K空间数据带分别进行相同的数据处理,获取一个以上处理后的数据带;
将所述一个以上处理后的数据带分别变换到图像域,获取一个以上优化图像;
计算所述一个以上优化图像中每两个优化图像之间的灰度互信息;
根据所述灰度互信息构造所述相关性矩阵,所述相关性矩阵为对称矩阵;
所述数据处理包括:
对所述K空间数据带依次进行相位校正、旋转和平移校正,获取校正后数据带;
将所述校正后数据带网格化到笛卡尔坐标系下,获取所述K空间数据带中心圆形重叠区域数据和所述K空间数据带的网格矩阵;
将所述K空间数据带中心圆形重叠区域数据填充到所述网格矩阵的中心,获取填充后网格矩阵;
将所述填充后网格矩阵进行低通滤波,获取所述处理后的数据带;
所述将所述填充后网格矩阵进行低通滤波的方法为:
将所述填充后网格矩阵乘以一个二维低通窗函数;
所述计算所述一个以上不同角度的K空间数据带中每两个K空间数据带之间的相关性,获取相关性矩阵的方法还包括:
将所述一个以上优化图像的灰阶分别压缩到预定阶数,获取一个以上压缩后图像;
则所述计算所述一个以上优化图像中每两个优化图像之间的灰度互信息,为计算一个以上压缩后图像中每两个压缩后图像之间的灰度互信息;
所述计算所述一个以上优化图像中每两个优化图像之间的灰度互信息的方法为:
I(Mm,Mn)=H(Mm)+H(Mn)-H(Mm,Mn)
其中,H(Mm)为图像Mm的边际熵,H(Mn)为图像Mn的边际熵,H(Mm,Mn)为图像Mm和图像Mn的联合熵,I(Mm,Mn)为图像Mm和图像Mn之间的灰度互信息;
所述根据所述相关性矩阵计算出每个所述K空间数据带的权重的方法为:
计算出所述相关性矩阵的最大特征值;
计算出所述最大特征值所对应的特征向量;
根据所述最大特征值所对应的特征向量计算出每个所述K空间数据带的权重:
其中,γn为所述特征向量中第n个元素,γmin为所述特征向量中的最小元素,γmax为所述特征向量中的最大元素,a和p为常量,wn为第n个K空间数据带的权重。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于奥泰医疗系统有限责任公司,未经奥泰医疗系统有限责任公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201810155551.3/1.html,转载请声明来源钻瓜专利网。
- 上一篇:一种动态磁共振稳健PCA成像方法
- 下一篇:磁共振检查固定装置