基于MR图像的PET吸收校正图像生成方法以及计算机程序.pdf

上传人:n****g 文档编号:135327 上传时间:2018-01-28 格式:PDF 页数:31 大小:13.18MB
返回 下载 相关 举报
摘要
申请专利号:

CN201380017570.5

申请日:

2013.03.28

公开号:

CN104220898A

公开日:

2014.12.17

当前法律状态:

授权

有效性:

有权

法律详情:

授权|||专利申请权的转移IPC(主分类):G01T 1/161登记生效日:20160727变更事项:申请人变更前权利人:独立行政法人放射线医学综合研究所变更后权利人:国立研究开发法人放射线医学综合研究所变更事项:地址变更前权利人:日本千叶县变更后权利人:日本千叶县变更事项:申请人变更前权利人:株式会社岛津制作所变更后权利人:株式会社岛津制作所|||著录事项变更IPC(主分类):G01T 1/161变更事项:申请人变更前:国立研究开发法人放射线医学综合研究所变更后:国立研究开发法人量子科学技术研究开发机构变更事项:地址变更前:日本千叶县变更后:日本千叶县变更事项:申请人变更前:株式会社岛津制作所变更后:株式会社岛津制作所|||实质审查的生效IPC(主分类):G01T 1/161申请日:20130328|||公开

IPC分类号:

G01T1/161; A61B5/055; A61B6/03

主分类号:

G01T1/161

申请人:

独立行政法人放射线医学综合研究所; 株式会社岛津制作所

发明人:

山谷泰贺; 小畠隆行; 菅干生; 川口拓之; 山川善之

地址:

日本千叶县

优先权:

2012.03.28 JP 2012-074906

专利代理机构:

北京林达刘知识产权代理事务所(普通合伙) 11277

代理人:

刘新宇

PDF下载: PDF下载
内容摘要

当根据MR图像生成PET吸收校正图像时,根据像素值对通过MRI拍摄到的MR图像进行区域分割,参照既有的放射线吸收校正表来求出能够将放射线的射线衰减系数视为一致的区域的放射线吸收校正值,参照标准图像来求出包括放射线的衰减系数不同的多个组织的区域的放射线吸收校正值。由此,生成如下的PET吸收校正图像:虽然在MR图像中像素值相等,但能够区分放射线的衰减系数不同的组织,且还能够与个人差异、占位性病变(癌、脓肿等)、脏器缺损等病变部位对应。

权利要求书

1.  一种基于MR图像的PET吸收校正图像生成方法,其特征在于,
当根据MR图像即磁共振图像生成PET吸收校正图像即正电子放射断层扫描吸收校正图像时,
根据像素值对通过磁共振成像拍摄到的MR图像进行区域分割,参照既有的放射线吸收校正值表来求出能够将放射线的射线衰减系数视为一致的区域的放射线吸收校正值,参照标准图像来求出包括放射线的射线衰减系数不同的多个组织的区域的放射线吸收校正值。

2.
  根据权利要求1所述的基于MR图像的PET吸收校正图像生成方法,其特征在于,
当参照上述标准图像时,与被检者的MR图像相应地使标准图像变形。

3.
  根据权利要求2所述的基于MR图像的PET吸收校正图像生成方法,其特征在于,
根据生物体组织整体的图像或者生物体组织的一部分的图像来进行上述标准图像的参照。

4.
  根据权利要求1至3中的任一项所述的基于MR图像的PET吸收校正图像生成方法,其特征在于,
上述标准图像是PET吸收校正图像或者CT图像即计算机断层扫描图像或者UTE图像即超短回波时间图像。

5.
  根据权利要求1至4中的任一项所述的基于MR图像的PET吸收校正图像生成方法,其特征在于,
使上述MR图像与标准图像的分辨率一致。

6.
  根据权利要求1至5中的任一项所述的基于MR图像的PET吸收校正图像生成方法,其特征在于,
根据被检者的个人信息来改变上述放射线吸收校正值表和标准图像中的至少一方。

7.
  根据权利要求1至5中的任一项所述的基于MR图像的PET吸收校正图像生成方法,其特征在于,
上述放射线吸收校正值表和标准图像中的至少一方反复使用被检者本人的放射线吸收校正值、PET吸收校正图像、CT图像或者UTE图像。

8.
  一种用于使计算机执行根据权利要求1至7中的任一项所述的方法的计算机程序。

说明书

基于MR图像的PET吸收校正图像生成方法以及计算机程序
技术领域
本发明涉及基于MR(magnetic:resonance磁共振)图像的PET吸收校正图像生成方法以及计算机程序,特别是涉及在根据MR图像生成PET吸收校正图像时能够制作高精度的PET吸收校正图像的基于MR图像的PET吸收校正图像生成方法以及用于执行该方法的计算机程序。
背景技术
近年来,X射线CT(computed tomography:计算机断层扫描)、MRI(magnetic resonance imaging:磁共振成像)等能够微创地观察人体内部的检查装置广泛普及,对医疗诊断作出巨大贡献。这些装置主要提供将生物体的组织构造可视化为断层图像、体数据而得到的形态图像。与此相对地,以PET(positron emission tomography:正电子放射断层扫描)、SPECT(single photon emission CT:单光子放射断层扫描)为代表的核医学成像装置是提供机能图像的装置,该机能图像是将葡萄糖消耗量、局部血流量、氧气消耗量、神经递质的受体分布等生理学的机能信息定量地可视化而得到的,近年来随着癌、痴呆症、动脉硬化性疾患等疾病的增加,推进了核医学成像装置在研究、临床上的普及。并且,PET作为推进将生物体分子的行为图像化的分子成像研究的有力的方法而被关注。
近年来,开发了PET和X射线CT能够同时进行拍摄的PET/CT,在临床时得到普及。由此,将生物体机能与体内组织进行对照的诊断成为可能。例如,在利用PET进行的癌的诊断中,仅脓肿部分以高信号输出,因此有时难以获知脓肿位于哪个脏器。因此,与出色地观察脏器的形态的X射线CT相配合是有益的。
最近,不与X射线CT相配合而与MRI相配合地进行诊断的PET/MRI受到关注。MRI不仅能够以较高的空间分辨率观察生物体,还具有以下优点:与 X射线CT相比软组织的对比度优良、能够通过血循环动态的成像、MR分光法进行代谢产物浓度测量以获得机能图像。另外,PET/MRI还具有能够避免在PET/CT中成为问题的辐射照射等诸多优点,因此PET/MRI的实现备受期待(参照非专利文献1)。
PET装置根据从放射性药剂释放的湮灭放射线到达检测器来获得信息。在PET图像重构时,使用向180度方向释放的湮灭放射线检测数据,但该湮灭放射线在通过身体的各个组织而到达检测器的期间被吸收。因此,在被检体的深部,定量性发生较大紊乱。为了获得定量的药剂分布需要校正湮灭放射线的吸收。基于与PET的数据收集不同的如图1所例示那样一边使射线源12绕着被检体10旋转一边用检测器14检测的透射测量(称为透射扫描),来制作现有的PET图像重构时的吸收校正所需的放射线吸收图像(μ映射图)(参照专利文献1~3)。另外,关于PET/CT中的衰减校正,不进行上述透射扫描,而是日常地进行从X射线CT图像变换为μ映射图。
另一方面,MRI收集从组织的质子(氢原子核)获得的信号,因此不能直接获得各组织的放射线吸收率,在当前的PET/MRI的设计概念中不包括与X射线CT相同的透射线源,因此需要取代透射扫描的μ映射图的制作方法。
在此,作为迄今为止使用了MR图像的μ映射图的制作方法,提出了区域分割法和标准图像参照法,其中,该区域分割法是如图2所例示那样,按每个组织例如将MR图像m区域分割为高信号区域(软组织)、中信号区域(水)、低信号区域(空气和骨骼),代入组织固有的μ值来制作μ映射图(参照非专利文献2、3、4),该标准图像参照法是如图3所例示那样,利用仿射变换等使标准图像(标准MR图像ms、标准μ映射图μs)变形以与患者的MR图像m一致(参照非专利文献5、6)。
专利文献1:日本特开2005-283421号公报
专利文献2:日本特开2006-284346号公报
专利文献3:日本特开2008-22930号公报
非专利文献1:Bernd Pichler,Bernhard Scholkopf and Thomas Beyer  Matthias Hofmann,“Towards quantitative PET/MRI:a review of MR-based attenuation correction techniques,”EUROPEAN JOURNAL OF NUCLEAR MEDICINE AND MOLECULAR IMAGING,vol.36(Supplement1),pp.93-104,Mar.2009.
非专利文献2:E.RotaKops,P.Qin,M.Mueller-Veggian,and H.Herzog,“Attenuation Correction of PET Scanning Based on MR Images,”Nucl.Sci.Symp.&Med.Imag.Conference,10-112006.
非专利文献3:E.Rota Kops and H.Herzog,“Towards an MRI based attenuation correction for brain MR-PET,”2007BrainPET Conference,May2007.
非专利文献4:H.Zaidi,M.-L.Montandon,and D.O.Slosman,“Magnetic resonance imaging-guided attenuation and scatter corrections in threedimensional brain positron emission tomography,”Med.Phys.,vol.30,pp.937-948,2003.
非专利文献5:Marie-Louise Montandon and Habib Zaidi,“Atlas-guided non-uniform attenuation correction in cerebral3D PET imaging,”Neuroimage,vol.25,no.1,pp.278-286,Mar.2005.
非专利文献6:Marie-Louise Montandon and Habib Zaidi,“Quantitative analysis of template-based attenuation compensation in3D brain PET,”Computerized Medical Imaging and Graphics,vol.31,pp.28-38,Jan.2007.
发明内容
发明要解决的问题
然而,如图2所例示那样,前者的区域分割法存在以下问题:μ值大不相同的组织(在示例的情况下为骨骼和空气)在MR图像中均是低信号,因此不能进行区分。另外,具有与软组织相等的μ值的血液有时也不能与空气、骨骼 进行区分。
另一方面,如图3所例示那样,后者的标准图像参照法具有以下问题:不能与个人差异、占位性病变(癌、脓肿等)、脏器缺损等病变部位对应。
本发明是为了消除上述现有的问题而完成的,其目的在于生成如下的PET吸收校正图像:能够区分μ值大不相同的组织,且还能够与个人差异、占位性病变(癌、脓肿等)、脏器缺损等病变部位对应。
用于解决问题的方案
本发明通过以下方式来解决上述问题:当根据MR图像生成PET吸收校正图像时,根据像素值对通过MRI拍摄到的MR图像进行区域分割,参照既有的放射线吸收校正值表来求出能够将放射线的射线衰减系数视为一致的区域的放射线吸收校正值,参照标准图像来求出包括放射线的射线衰减系数不同的多个组织的区域的放射线吸收校正值。
在此,当参照上述标准图像时,能够与被检者的MR图像相应地使标准图像变形。
另外,能够根据生物体组织整体的图像或者生物体组织的一部分组织的图像来进行上述标准图像的参照。
另外,能够将上述标准图像设为PET吸收校正图像、CT图像或者UTE(Ultra short Echo Time:超短回波时间)图像。在此,UTE图像是指能够根据超短回波时间(Ultra short TE)内拍摄到的MRI图像制作的、骨骼呈高信号值的图像。
另外,能够使上述MR图像与标准图像的分辨率一致。
另外,能够根据被检者的个人信息(DNA、年龄、性别、身高、体重、出生地、居住地、病历)来改变上述放射线吸收校正值表和标准图像中的至少一方。
另外,上述放射线吸收校正值表和标准图像中的至少一方还能够反复使用被检者本人的放射线吸收校正值、PET吸收校正图像、CT图像或者UTE图像。
本发明还提供用于使计算机执行上述任一项所记载的方法的计算机程序。
发明的效果
根据本发明,对于能够通过参照MR图像的像素值来区分放射线吸收校正值(也称为μ值)的差异的区域(例如作为高信号区域的软组织和作为中信号区域的水),通过区域分割法来代入组织固有的μ值,对于不能利用MR图像的像素值进行区分的区域(例如包括骨骼和空气的低信号区域),通过标准图像参照法参照标准图像来求出μ值,因此能够在对不能利用MR图像的像素值进行区分的区域分配正确的μ值的同时,还与个人差异、脑的缺损部分等病变部位对应。
附图说明
图1是表示现有的透射法的概要的图。
图2是同样地表示区域分割法的图。
图3是同样地表示标准图像参照法的图。
图4是表示本发明的第一实施方式的处理过程的流程图。
图5是同样地表示图像的例子的图。
图6是表示本发明的第二实施方式的处理过程的流程图。
图7是表示本发明的第三实施方式的处理过程的流程图。
图8是表示本发明的第四实施方式的处理过程的流程图。
图9是表示本发明的第五实施方式的处理过程的流程图。
图10是表示本发明的第六实施方式的处理过程的流程图。
图11是根据第四实施方式~第六实施方式将考虑了分辨率的差时的μ映射图的例子进行比较并示出的图
图12是表示本发明的第七实施方式的处理过程的流程图。
图13是表示本发明的第八实施方式的处理过程的流程图。
图14是表示本发明的第九实施方式的处理过程的流程图。
图15是表示本发明的第十实施方式的处理过程的流程图。
图16是表示在第十实施方式中分类为骨骼和空气的区域的情形的图。
图17是表示本发明的第十一实施方式的处理过程的流程图。
图18是将现有方法和本发明方法对疑似鼻窦炎的处理结果进行比较并示出的图。
图19是表示本发明的第十二实施方式的处理过程的流程图。
图20是表示本发明的第十三实施方式的处理过程的流程图。
图21是表示本发明的第十四实施方式的处理过程的流程图。
图22是表示本发明的第十五实施方式的处理过程的流程图。
图23是表示本发明的第十六实施方式的处理过程的流程图。
图24是表示本发明的第十七实施方式的处理过程的流程图。
具体实施方式
下面,参照附图来详细地说明本发明的实施方式。
本发明的第一实施方式按如图4所示的过程来实施。
即,首先,在步骤S100中与现有的区域分割法同样地获取MR图像m。
接着,在步骤S110中,基于MR图像的像素值,如图5所例示那样例如分割为对应于软组织的高信号区域、对应于水的中信号区域、对应于空气和骨骼的低信号区域。
接着,在步骤S120中,针对各信号区域,判断各信号区域是否是包括射线衰减系数不同的多个组织的区域。当判断结果为“否”、例如判断为是对应于软组织的高信号区域或者对应于水的中信号区域时,对高信号区域分配软组织的射线衰减系数(例如μ=0.095),对中信号区域分配水的射线衰减系数(例如μ=0.097)。
另一方面,在步骤S200中,与图3所示的现有的标准图像参照法同样地,利用标准MR图像ms和标准射线衰减图像μs的数据库20求出用于进行在步骤S100中获取到的MR图像m的位置对准的使得m=F{ms}的变换F{}。
接着,在步骤S210中,同样地进行标准射线衰减图像的形状变换μ=F{μs}。
当步骤S120的判断结果为“是”、判断为是包括射线衰减系数不同的多个组织的区域、例如不能区分空气和骨骼的低信号区域时,进入步骤S300,将在步骤S210中进行形状变换后的标准射线衰减图像μs的像素值分配给低信号区域,由此对空气和骨骼分配不同的射线衰减系数。
此外,在第一实施方式中获取的MR图像与射线衰减图像的分辨率存在差异,因此直接进行合成有可能难以辨认。因此,在图6中示出解决了这种问题的第二实施方式。
在本实施方式中,在与第一实施方式相同的处理的过程中,在步骤S140中对分辨率高的MR图像施加模糊,以成为与射线衰减图像相同程度的分辨率。步骤S140的模糊如果在紧随步骤S130之后且应用最终的标准图像的μ值之前则能够利用。
根据本实施方式,步骤S130的结果(射线衰减图像)与标准图像(MR图像)的分辨率统一,因此能够显示正确的图像。但是,在原来的MR图像与射线衰减图像的分辨率的差大的情况下,强行施加模糊,如图6的示意图所示那样,在骨骼区域的周边存在射线源弱系数降低的情况。
此外,施加模糊的方法并不限定于第二实施方式,还能够如图7所示的第三实施方式那样,在即将结束前的步骤S150中施加模糊。
根据本实施方式,能够与PET检测器的分辨率相匹配,因此能够接近实际的射线衰减系数的测量。另外,如果不施加模糊,则能够削减工时且降低计算成本。
接着,在图8中示出利用另一方法吸收分辨率的差的本发明的第四实施方式。
关于本实施方式,在与图6所示的第二实施方式相同的处理中,当步骤S120的判断结果为“是”时,在步骤S290中针对低信号区域使区域膨胀,之后在步骤S300中分配形状变换后的标准射线衰减图像的像素值。
根据本实施方式,能够抑制骨骼区域周边的射线源弱系数降低,能够将射线衰减图像与MR图像更加顺畅地连接。
在本实施方式中,也能够与第三实施方式同样地在即将结束之前施加模糊。
接着,在图9中示出利用又一方法吸收分辨率的差的本发明的第五实施方式。
关于本实施方式,在与图6所示的第二实施方式相同的处理中,在步骤S140中进行模糊前后,通过步骤S280、S300来分配形状变换后的标准射线衰减图像的像素值。
根据本实施方式,能够一边保持骨骼的射线源弱系数(图像的高度H),一边进一步抑制骨骼区域周边的射线源弱系数降低。
由第二实施方式至第五实施方式的方法产生的效果还依赖于原来的MR图像与射线衰减图像的分辨率的差的程度,因此需要调整步骤S140的模糊的量、第四实施方式的S290的区域膨胀处理的量。
接着,在图10中示出将利用再一方法吸收分辨率的差的本发明的第六实施方式。
关于本实施方式,在与图9所示的第五实施方式相同的处理中,当在模糊处理S140之后进行分配时,使得进行与图8所示的第四实施方式相同的区域膨胀处理S290。与第二实施方式至第五实施方式的方法相比,第六实施方式的方法是能够稳定地校正MR图像与射线衰减图像的分辨率的差的方法。
在图11中将MRI图像、真正的μ映射图、由第四实施方式~第六实施方式得到的处理结果进行比较并示出。
此外,在上述实施方式中,均在进行MR图像与标准射线衰减图像的位置对准之后进行射线衰减图像的形状变换,但在图12中示出先提取骨骼并对该部分进行形状变换的本发明的第七实施方式。
在本实施方式中,首先在步骤S400中从标准MR图像ms和标准射线衰减图像μs的数据库20事先提取骨骼的部分μsb
接着,在与第一实施方式相同的步骤S200中进行MR图像的位置对准,因此求出使得m=F{ms}的变换F{}。
接着,在步骤S420中,针对在步骤S400中提取出的骨骼的部分变换射线源弱图像的形状,并与在步骤S130中求出的射线衰减系数进行合成。
在本实施方式中,先提取骨骼的部分,仅对骨骼的部分进行了形状变换,因此不需要步骤S120的判断并且是高精度的。
此外,在上述实施方式中,均利用了标准射线衰减图像μs的数据库20,但也能够如图13所示的第八实施方式那样取而代之地利用标准CT图像cs的数据库30。
在本实施方式中,在与第一实施方式相同的步骤S200中,进行与在步骤S100中获取到的MR图像m的位置对准,因此求出使得m=F{ms}的变换F{}。
接着,在步骤S520中,同样求出标准CT图像的形状变换c=F{cs}。
接着,在步骤S530中,求出从CT图像cs向射线衰减图像μ的变换。
而且,在步骤S120的判断结果为“是”的情况下,例如针对低信号区域,在步骤S300中分配形状变换后的标准射线衰减图像的像素值。
根据本实施方式,能够利用与MR图像之间分辨率的差小的CT图像来分配射线衰减系数,但即使是CT图像,也由于在标准图像中进行平均化处理而导致分辨率降低,因此如第三实施方式~第六实施方式那样的吸收分辨率的差的处理同样有效。
接着,在图14中示出同样利用了CT图像的本发明的第九实施方式。
关于本实施方式,在与图13所示的第八实施方式相同的处理过程中,针对在步骤S120中被认为包括射线衰减系数不同的多个组织的区域,在步骤S540中参照形状变换后的CT图像的像素值,在步骤S550中仅对包括多个组织的部分进行CT图像的区域分割,之后在步骤130中进行组织对放射线的射线衰减系数的分配。
具体地说,如图15所示的第十实施方式那样,在步骤S540结束后,进入步骤S560,判断在CT图像中是否为高信号。在判断结果为“是”、判断为是 骨骼的区域时,进入步骤S570,分配骨骼的射线衰减系数。
另一方面,在步骤S560的判断结果为“否”、判断为是空气的区域时,进入步骤S580,分配空气的射线衰减系数。
在本实施方式中,如图16所例示那样,在CT图像中识别了骨骼和空气,因此能够高精度地进行分配。
接着,在图17中示出同样利用CT图像来提取骨骼的本发明的第十一实施方式。
关于本实施方式,在与图12所示的第七实施方式相同的处理中,在步骤S600中从标准MR图像ms和标准CT图像cs的数据库30提取骨骼csb,接着在步骤S610中分配骨骼的射线衰减系数μsb,接着在步骤S200中进行MR图像的位置对准,在步骤S420中进行标准射线衰减图像的形状变换。
此外,变换F{}的估计并不完全,在产生了些许间隙、重叠的情况下,还能够根据需要而利用步骤S640通过放大缩小处理、最近相邻内插来进行校正的处理。
在图18中将通过现有方法和本发明方法对疑似鼻窦炎的被检者的图像进行处理而得到的结果进行比较并示出。
根据本发明方法,显然能够获得比区域分割法、标准图像参照法更优良的相关性。
图19示出在图14所示的第九实施方式中利用标准射线衰减图像来代替标准CT图像的第十二实施方式。与第九实施方式的不同之处在于,利用标准MR图像和标准射线衰减图像的数据库30作为数据库,在步骤S522、S542、S552中利用标准射线衰减图像来代替标准CT图像。其它处与第九实施方式相同,因此省略说明。
接着,对利用标准UTE图像来代替标准射线衰减图像或者标准CT图像的实施方式进行说明。
图20示出在图13所示的第八实施方式中利用标准UTE图像来代替标准CT图像的第十三实施方式。与第八实施方式的不同之处在于,利用标准MR 图像和标准UTE图像的数据库40作为数据库,在步骤S524、S534中利用标准UTE图像来代替标准CT图像。其它处与第八实施方式相同,因此省略说明。
图21示出在图14所示的第九实施方式中利用标准UTE图像来代替标准CT图像的第十四实施方式。与第九实施方式的不同之处在于,利用标准MR图像和标准UTE图像的数据库40作为数据库,在步骤S524、S544、S554中利用标准UTE图像来代替标准CT图像。其它处与第九实施方式相同,因此省略说明。
图22示出在图17所示的第十一实施方式中利用标准UTE图像来代替标准CT图像的第十五实施方式。与第十一实施方式的不同之处在于,利用了标准MR图像和标准UTE图像的数据库40作为数据库。其它处与第十一实施方式相同,因此省略说明。
并且,还能够通过如图23所示的第十六实施方式、图24所示的第十七实施方式那样基于被检者的个人信息(DNA、年龄、性别、身高、体重、出生地、居住地、病历等)来定制标准图像、射线衰减系数,来进一步提高精度。在图23、24中,S700是获取个人信息的步骤,S710是生成反映了个人信息的射线衰减系数的步骤,S729是生成反映了个人信息的标准MR图像和标准射线衰减图像的步骤。其它处与图4所示的第一实施方式相同,因此省略说明。
此外,与标准MR图像一起使用的标准图像并不限定于射线源弱图像,还能够通过按照数据库进行如下改写来应用于其它方式,在标准MR图像和标准CT图像的数据库30时,将S210处替换为图13的S520和S530,在标准MR图像和标准UTE图像的数据库40时,将S210处替换为图20的S524和S534即可。
另外,上述放射线吸收校正值表和标准图像中的至少一方还能够反复使用被检者本人的放射线吸收校正值、射线源弱图像、CT图像或者UTE图像。
此外,在上述实施方式中,本发明适用于脑图像,但本发明的适用对象并不限于此,区域分割数也并不限定于高信号区域、中信号区域、低信号区域这三个区域。
产业上的可利用性
本发明能够应用于使PET与MRI相结合的PET/MRI。
通过参照而将2012年3月28日申请的日本申请号2007-074906的说明书、附图以及权利要求书的公开内容整体引用到本说明书中。
附图标记说明
10:被检体;12:射线源;14:检测器;20:标准MR图像和标准射线衰减图像的数据库;30:标准MR图像和标准CT图像的数据库;40:标准MR图像和标准UTE图像的数据库。

基于MR图像的PET吸收校正图像生成方法以及计算机程序.pdf_第1页
第1页 / 共31页
基于MR图像的PET吸收校正图像生成方法以及计算机程序.pdf_第2页
第2页 / 共31页
基于MR图像的PET吸收校正图像生成方法以及计算机程序.pdf_第3页
第3页 / 共31页
点击查看更多>>
资源描述

《基于MR图像的PET吸收校正图像生成方法以及计算机程序.pdf》由会员分享,可在线阅读,更多相关《基于MR图像的PET吸收校正图像生成方法以及计算机程序.pdf(31页珍藏版)》请在专利查询网上搜索。

1、10申请公布号CN104220898A43申请公布日20141217CN104220898A21申请号201380017570522申请日20130328201207490620120328JPG01T1/161200601A61B5/055200601A61B6/0320060171申请人独立行政法人放射线医学综合研究所地址日本千叶县申请人株式会社岛津制作所72发明人山谷泰贺小畠隆行菅干生川口拓之山川善之74专利代理机构北京林达刘知识产权代理事务所普通合伙11277代理人刘新宇54发明名称基于MR图像的PET吸收校正图像生成方法以及计算机程序57摘要当根据MR图像生成PET吸收校正图像时,根。

2、据像素值对通过MRI拍摄到的MR图像进行区域分割,参照既有的放射线吸收校正表来求出能够将放射线的射线衰减系数视为一致的区域的放射线吸收校正值,参照标准图像来求出包括放射线的衰减系数不同的多个组织的区域的放射线吸收校正值。由此,生成如下的PET吸收校正图像虽然在MR图像中像素值相等,但能够区分放射线的衰减系数不同的组织,且还能够与个人差异、占位性病变癌、脓肿等、脏器缺损等病变部位对应。30优先权数据85PCT国际申请进入国家阶段日2014092886PCT国际申请的申请数据PCT/JP2013/0592352013032887PCT国际申请的公布数据WO2013/147013JA20131003。

3、51INTCL权利要求书1页说明书7页附图22页19中华人民共和国国家知识产权局12发明专利申请权利要求书1页说明书7页附图22页10申请公布号CN104220898ACN104220898A1/1页21一种基于MR图像的PET吸收校正图像生成方法,其特征在于,当根据MR图像即磁共振图像生成PET吸收校正图像即正电子放射断层扫描吸收校正图像时,根据像素值对通过磁共振成像拍摄到的MR图像进行区域分割,参照既有的放射线吸收校正值表来求出能够将放射线的射线衰减系数视为一致的区域的放射线吸收校正值,参照标准图像来求出包括放射线的射线衰减系数不同的多个组织的区域的放射线吸收校正值。2根据权利要求1所述的。

4、基于MR图像的PET吸收校正图像生成方法,其特征在于,当参照上述标准图像时,与被检者的MR图像相应地使标准图像变形。3根据权利要求2所述的基于MR图像的PET吸收校正图像生成方法,其特征在于,根据生物体组织整体的图像或者生物体组织的一部分的图像来进行上述标准图像的参照。4根据权利要求1至3中的任一项所述的基于MR图像的PET吸收校正图像生成方法,其特征在于,上述标准图像是PET吸收校正图像或者CT图像即计算机断层扫描图像或者UTE图像即超短回波时间图像。5根据权利要求1至4中的任一项所述的基于MR图像的PET吸收校正图像生成方法,其特征在于,使上述MR图像与标准图像的分辨率一致。6根据权利要求。

5、1至5中的任一项所述的基于MR图像的PET吸收校正图像生成方法,其特征在于,根据被检者的个人信息来改变上述放射线吸收校正值表和标准图像中的至少一方。7根据权利要求1至5中的任一项所述的基于MR图像的PET吸收校正图像生成方法,其特征在于,上述放射线吸收校正值表和标准图像中的至少一方反复使用被检者本人的放射线吸收校正值、PET吸收校正图像、CT图像或者UTE图像。8一种用于使计算机执行根据权利要求1至7中的任一项所述的方法的计算机程序。权利要求书CN104220898A1/7页3基于MR图像的PET吸收校正图像生成方法以及计算机程序技术领域0001本发明涉及基于MRMAGNETICRESONAN。

6、CE磁共振图像的PET吸收校正图像生成方法以及计算机程序,特别是涉及在根据MR图像生成PET吸收校正图像时能够制作高精度的PET吸收校正图像的基于MR图像的PET吸收校正图像生成方法以及用于执行该方法的计算机程序。背景技术0002近年来,X射线CTCOMPUTEDTOMOGRAPHY计算机断层扫描、MRIMAGNETICRESONANCEIMAGING磁共振成像等能够微创地观察人体内部的检查装置广泛普及,对医疗诊断作出巨大贡献。这些装置主要提供将生物体的组织构造可视化为断层图像、体数据而得到的形态图像。与此相对地,以PETPOSITRONEMISSIONTOMOGRAPHY正电子放射断层扫描、。

7、SPECTSINGLEPHOTONEMISSIONCT单光子放射断层扫描为代表的核医学成像装置是提供机能图像的装置,该机能图像是将葡萄糖消耗量、局部血流量、氧气消耗量、神经递质的受体分布等生理学的机能信息定量地可视化而得到的,近年来随着癌、痴呆症、动脉硬化性疾患等疾病的增加,推进了核医学成像装置在研究、临床上的普及。并且,PET作为推进将生物体分子的行为图像化的分子成像研究的有力的方法而被关注。0003近年来,开发了PET和X射线CT能够同时进行拍摄的PET/CT,在临床时得到普及。由此,将生物体机能与体内组织进行对照的诊断成为可能。例如,在利用PET进行的癌的诊断中,仅脓肿部分以高信号输出,。

8、因此有时难以获知脓肿位于哪个脏器。因此,与出色地观察脏器的形态的X射线CT相配合是有益的。0004最近,不与X射线CT相配合而与MRI相配合地进行诊断的PET/MRI受到关注。MRI不仅能够以较高的空间分辨率观察生物体,还具有以下优点与X射线CT相比软组织的对比度优良、能够通过血循环动态的成像、MR分光法进行代谢产物浓度测量以获得机能图像。另外,PET/MRI还具有能够避免在PET/CT中成为问题的辐射照射等诸多优点,因此PET/MRI的实现备受期待参照非专利文献1。0005PET装置根据从放射性药剂释放的湮灭放射线到达检测器来获得信息。在PET图像重构时,使用向180度方向释放的湮灭放射线检。

9、测数据,但该湮灭放射线在通过身体的各个组织而到达检测器的期间被吸收。因此,在被检体的深部,定量性发生较大紊乱。为了获得定量的药剂分布需要校正湮灭放射线的吸收。基于与PET的数据收集不同的如图1所例示那样一边使射线源12绕着被检体10旋转一边用检测器14检测的透射测量称为透射扫描,来制作现有的PET图像重构时的吸收校正所需的放射线吸收图像映射图参照专利文献13。另外,关于PET/CT中的衰减校正,不进行上述透射扫描,而是日常地进行从X射线CT图像变换为映射图。0006另一方面,MRI收集从组织的质子氢原子核获得的信号,因此不能直接获得各组织的放射线吸收率,在当前的PET/MRI的设计概念中不包括。

10、与X射线CT相同的透射线源,因此需要取代透射扫描的映射图的制作方法。说明书CN104220898A2/7页40007在此,作为迄今为止使用了MR图像的映射图的制作方法,提出了区域分割法和标准图像参照法,其中,该区域分割法是如图2所例示那样,按每个组织例如将MR图像M区域分割为高信号区域软组织、中信号区域水、低信号区域空气和骨骼,代入组织固有的值来制作映射图参照非专利文献2、3、4,该标准图像参照法是如图3所例示那样,利用仿射变换等使标准图像标准MR图像MS、标准映射图S变形以与患者的MR图像M一致参照非专利文献5、6。0008专利文献1日本特开2005283421号公报0009专利文献2日本特。

11、开2006284346号公报0010专利文献3日本特开200822930号公报0011非专利文献1BERNDPICHLER,BERNHARDSCHOLKOPFANDTHOMASBEYERMATTHIASHOFMANN,“TOWARDSQUANTITATIVEPET/MRIAREVIEWOFMRBASEDATTENUATIONCORRECTIONTECHNIQUES,”EUROPEANJOURNALOFNUCLEARMEDICINEANDMOLECULARIMAGING,VOL36SUPPLEMENT1,PP93104,MAR20090012非专利文献2EROTAKOPS,PQIN,MMUELL。

12、ERVEGGIAN,ANDHHERZOG,“ATTENUATIONCORRECTIONOFPETSCANNINGBASEDONMRIMAGES,”NUCLSCISYMPMEDIMAGCONFERENCE,101120060013非专利文献3EROTAKOPSANDHHERZOG,“TOWARDSANMRIBASEDATTENUATIONCORRECTIONFORBRAINMRPET,”2007BRAINPETCONFERENCE,MAY20070014非专利文献4HZAIDI,MLMONTANDON,ANDDOSLOSMAN,“MAGNETICRESONANCEIMAGINGGUIDEDATT。

13、ENUATIONANDSCATTERCORRECTIONSINTHREEDIMENSIONALBRAINPOSITRONEMISSIONTOMOGRAPHY,”MEDPHYS,VOL30,PP937948,20030015非专利文献5MARIELOUISEMONTANDONANDHABIBZAIDI,“ATLASGUIDEDNONUNIFORMATTENUATIONCORRECTIONINCEREBRAL3DPETIMAGING,”NEUROIMAGE,VOL25,NO1,PP278286,MAR20050016非专利文献6MARIELOUISEMONTANDONANDHABIBZAIDI,“。

14、QUANTITATIVEANALYSISOFTEMPLATEBASEDATTENUATIONCOMPENSATIONIN3DBRAINPET,”COMPUTERIZEDMEDICALIMAGINGANDGRAPHICS,VOL31,PP2838,JAN2007发明内容0017发明要解决的问题0018然而,如图2所例示那样,前者的区域分割法存在以下问题值大不相同的组织在示例的情况下为骨骼和空气在MR图像中均是低信号,因此不能进行区分。另外,具有与软组织相等的值的血液有时也不能与空气、骨骼进行区分。0019另一方面,如图3所例示那样,后者的标准图像参照法具有以下问题不能与个人差异、占位性病变癌、脓。

15、肿等、脏器缺损等病变部位对应。0020本发明是为了消除上述现有的问题而完成的,其目的在于生成如下的PET吸收校正图像能够区分值大不相同的组织,且还能够与个人差异、占位性病变癌、脓肿等、脏器缺损等病变部位对应。说明书CN104220898A3/7页50021用于解决问题的方案0022本发明通过以下方式来解决上述问题当根据MR图像生成PET吸收校正图像时,根据像素值对通过MRI拍摄到的MR图像进行区域分割,参照既有的放射线吸收校正值表来求出能够将放射线的射线衰减系数视为一致的区域的放射线吸收校正值,参照标准图像来求出包括放射线的射线衰减系数不同的多个组织的区域的放射线吸收校正值。0023在此,当参。

16、照上述标准图像时,能够与被检者的MR图像相应地使标准图像变形。0024另外,能够根据生物体组织整体的图像或者生物体组织的一部分组织的图像来进行上述标准图像的参照。0025另外,能够将上述标准图像设为PET吸收校正图像、CT图像或者UTEULTRASHORTECHOTIME超短回波时间图像。在此,UTE图像是指能够根据超短回波时间ULTRASHORTTE内拍摄到的MRI图像制作的、骨骼呈高信号值的图像。0026另外,能够使上述MR图像与标准图像的分辨率一致。0027另外,能够根据被检者的个人信息DNA、年龄、性别、身高、体重、出生地、居住地、病历来改变上述放射线吸收校正值表和标准图像中的至少一方。

17、。0028另外,上述放射线吸收校正值表和标准图像中的至少一方还能够反复使用被检者本人的放射线吸收校正值、PET吸收校正图像、CT图像或者UTE图像。0029本发明还提供用于使计算机执行上述任一项所记载的方法的计算机程序。0030发明的效果0031根据本发明,对于能够通过参照MR图像的像素值来区分放射线吸收校正值也称为值的差异的区域例如作为高信号区域的软组织和作为中信号区域的水,通过区域分割法来代入组织固有的值,对于不能利用MR图像的像素值进行区分的区域例如包括骨骼和空气的低信号区域,通过标准图像参照法参照标准图像来求出值,因此能够在对不能利用MR图像的像素值进行区分的区域分配正确的值的同时,还。

18、与个人差异、脑的缺损部分等病变部位对应。附图说明0032图1是表示现有的透射法的概要的图。0033图2是同样地表示区域分割法的图。0034图3是同样地表示标准图像参照法的图。0035图4是表示本发明的第一实施方式的处理过程的流程图。0036图5是同样地表示图像的例子的图。0037图6是表示本发明的第二实施方式的处理过程的流程图。0038图7是表示本发明的第三实施方式的处理过程的流程图。0039图8是表示本发明的第四实施方式的处理过程的流程图。0040图9是表示本发明的第五实施方式的处理过程的流程图。0041图10是表示本发明的第六实施方式的处理过程的流程图。0042图11是根据第四实施方式第六。

19、实施方式将考虑了分辨率的差时的映射图的例子进行比较并示出的图0043图12是表示本发明的第七实施方式的处理过程的流程图。说明书CN104220898A4/7页60044图13是表示本发明的第八实施方式的处理过程的流程图。0045图14是表示本发明的第九实施方式的处理过程的流程图。0046图15是表示本发明的第十实施方式的处理过程的流程图。0047图16是表示在第十实施方式中分类为骨骼和空气的区域的情形的图。0048图17是表示本发明的第十一实施方式的处理过程的流程图。0049图18是将现有方法和本发明方法对疑似鼻窦炎的处理结果进行比较并示出的图。0050图19是表示本发明的第十二实施方式的处理。

20、过程的流程图。0051图20是表示本发明的第十三实施方式的处理过程的流程图。0052图21是表示本发明的第十四实施方式的处理过程的流程图。0053图22是表示本发明的第十五实施方式的处理过程的流程图。0054图23是表示本发明的第十六实施方式的处理过程的流程图。0055图24是表示本发明的第十七实施方式的处理过程的流程图。具体实施方式0056下面,参照附图来详细地说明本发明的实施方式。0057本发明的第一实施方式按如图4所示的过程来实施。0058即,首先,在步骤S100中与现有的区域分割法同样地获取MR图像M。0059接着,在步骤S110中,基于MR图像的像素值,如图5所例示那样例如分割为对应。

21、于软组织的高信号区域、对应于水的中信号区域、对应于空气和骨骼的低信号区域。0060接着,在步骤S120中,针对各信号区域,判断各信号区域是否是包括射线衰减系数不同的多个组织的区域。当判断结果为“否”、例如判断为是对应于软组织的高信号区域或者对应于水的中信号区域时,对高信号区域分配软组织的射线衰减系数例如0095,对中信号区域分配水的射线衰减系数例如0097。0061另一方面,在步骤S200中,与图3所示的现有的标准图像参照法同样地,利用标准MR图像MS和标准射线衰减图像S的数据库20求出用于进行在步骤S100中获取到的MR图像M的位置对准的使得MFMS的变换F。0062接着,在步骤S210中,。

22、同样地进行标准射线衰减图像的形状变换FS。0063当步骤S120的判断结果为“是”、判断为是包括射线衰减系数不同的多个组织的区域、例如不能区分空气和骨骼的低信号区域时,进入步骤S300,将在步骤S210中进行形状变换后的标准射线衰减图像S的像素值分配给低信号区域,由此对空气和骨骼分配不同的射线衰减系数。0064此外,在第一实施方式中获取的MR图像与射线衰减图像的分辨率存在差异,因此直接进行合成有可能难以辨认。因此,在图6中示出解决了这种问题的第二实施方式。0065在本实施方式中,在与第一实施方式相同的处理的过程中,在步骤S140中对分辨率高的MR图像施加模糊,以成为与射线衰减图像相同程度的分辨。

23、率。步骤S140的模糊如果在紧随步骤S130之后且应用最终的标准图像的值之前则能够利用。0066根据本实施方式,步骤S130的结果射线衰减图像与标准图像MR图像的分辨率统一,因此能够显示正确的图像。但是,在原来的MR图像与射线衰减图像的分辨率的说明书CN104220898A5/7页7差大的情况下,强行施加模糊,如图6的示意图所示那样,在骨骼区域的周边存在射线源弱系数降低的情况。0067此外,施加模糊的方法并不限定于第二实施方式,还能够如图7所示的第三实施方式那样,在即将结束前的步骤S150中施加模糊。0068根据本实施方式,能够与PET检测器的分辨率相匹配,因此能够接近实际的射线衰减系数的测量。

24、。另外,如果不施加模糊,则能够削减工时且降低计算成本。0069接着,在图8中示出利用另一方法吸收分辨率的差的本发明的第四实施方式。0070关于本实施方式,在与图6所示的第二实施方式相同的处理中,当步骤S120的判断结果为“是”时,在步骤S290中针对低信号区域使区域膨胀,之后在步骤S300中分配形状变换后的标准射线衰减图像的像素值。0071根据本实施方式,能够抑制骨骼区域周边的射线源弱系数降低,能够将射线衰减图像与MR图像更加顺畅地连接。0072在本实施方式中,也能够与第三实施方式同样地在即将结束之前施加模糊。0073接着,在图9中示出利用又一方法吸收分辨率的差的本发明的第五实施方式。0074。

25、关于本实施方式,在与图6所示的第二实施方式相同的处理中,在步骤S140中进行模糊前后,通过步骤S280、S300来分配形状变换后的标准射线衰减图像的像素值。0075根据本实施方式,能够一边保持骨骼的射线源弱系数图像的高度H,一边进一步抑制骨骼区域周边的射线源弱系数降低。0076由第二实施方式至第五实施方式的方法产生的效果还依赖于原来的MR图像与射线衰减图像的分辨率的差的程度,因此需要调整步骤S140的模糊的量、第四实施方式的S290的区域膨胀处理的量。0077接着,在图10中示出将利用再一方法吸收分辨率的差的本发明的第六实施方式。0078关于本实施方式,在与图9所示的第五实施方式相同的处理中,。

26、当在模糊处理S140之后进行分配时,使得进行与图8所示的第四实施方式相同的区域膨胀处理S290。与第二实施方式至第五实施方式的方法相比,第六实施方式的方法是能够稳定地校正MR图像与射线衰减图像的分辨率的差的方法。0079在图11中将MRI图像、真正的映射图、由第四实施方式第六实施方式得到的处理结果进行比较并示出。0080此外,在上述实施方式中,均在进行MR图像与标准射线衰减图像的位置对准之后进行射线衰减图像的形状变换,但在图12中示出先提取骨骼并对该部分进行形状变换的本发明的第七实施方式。0081在本实施方式中,首先在步骤S400中从标准MR图像MS和标准射线衰减图像S的数据库20事先提取骨骼。

27、的部分SB。0082接着,在与第一实施方式相同的步骤S200中进行MR图像的位置对准,因此求出使得MFMS的变换F。0083接着,在步骤S420中,针对在步骤S400中提取出的骨骼的部分变换射线源弱图像的形状,并与在步骤S130中求出的射线衰减系数进行合成。0084在本实施方式中,先提取骨骼的部分,仅对骨骼的部分进行了形状变换,因此不需要步骤S120的判断并且是高精度的。说明书CN104220898A6/7页80085此外,在上述实施方式中,均利用了标准射线衰减图像S的数据库20,但也能够如图13所示的第八实施方式那样取而代之地利用标准CT图像CS的数据库30。0086在本实施方式中,在与第一。

28、实施方式相同的步骤S200中,进行与在步骤S100中获取到的MR图像M的位置对准,因此求出使得MFMS的变换F。0087接着,在步骤S520中,同样求出标准CT图像的形状变换CFCS。0088接着,在步骤S530中,求出从CT图像CS向射线衰减图像的变换。0089而且,在步骤S120的判断结果为“是”的情况下,例如针对低信号区域,在步骤S300中分配形状变换后的标准射线衰减图像的像素值。0090根据本实施方式,能够利用与MR图像之间分辨率的差小的CT图像来分配射线衰减系数,但即使是CT图像,也由于在标准图像中进行平均化处理而导致分辨率降低,因此如第三实施方式第六实施方式那样的吸收分辨率的差的处。

29、理同样有效。0091接着,在图14中示出同样利用了CT图像的本发明的第九实施方式。0092关于本实施方式,在与图13所示的第八实施方式相同的处理过程中,针对在步骤S120中被认为包括射线衰减系数不同的多个组织的区域,在步骤S540中参照形状变换后的CT图像的像素值,在步骤S550中仅对包括多个组织的部分进行CT图像的区域分割,之后在步骤130中进行组织对放射线的射线衰减系数的分配。0093具体地说,如图15所示的第十实施方式那样,在步骤S540结束后,进入步骤S560,判断在CT图像中是否为高信号。在判断结果为“是”、判断为是骨骼的区域时,进入步骤S570,分配骨骼的射线衰减系数。0094另一。

30、方面,在步骤S560的判断结果为“否”、判断为是空气的区域时,进入步骤S580,分配空气的射线衰减系数。0095在本实施方式中,如图16所例示那样,在CT图像中识别了骨骼和空气,因此能够高精度地进行分配。0096接着,在图17中示出同样利用CT图像来提取骨骼的本发明的第十一实施方式。0097关于本实施方式,在与图12所示的第七实施方式相同的处理中,在步骤S600中从标准MR图像MS和标准CT图像CS的数据库30提取骨骼CSB,接着在步骤S610中分配骨骼的射线衰减系数SB,接着在步骤S200中进行MR图像的位置对准,在步骤S420中进行标准射线衰减图像的形状变换。0098此外,变换F的估计并不。

31、完全,在产生了些许间隙、重叠的情况下,还能够根据需要而利用步骤S640通过放大缩小处理、最近相邻内插来进行校正的处理。0099在图18中将通过现有方法和本发明方法对疑似鼻窦炎的被检者的图像进行处理而得到的结果进行比较并示出。0100根据本发明方法,显然能够获得比区域分割法、标准图像参照法更优良的相关性。0101图19示出在图14所示的第九实施方式中利用标准射线衰减图像来代替标准CT图像的第十二实施方式。与第九实施方式的不同之处在于,利用标准MR图像和标准射线衰减图像的数据库30作为数据库,在步骤S522、S542、S552中利用标准射线衰减图像来代替标准CT图像。其它处与第九实施方式相同,因此。

32、省略说明。0102接着,对利用标准UTE图像来代替标准射线衰减图像或者标准CT图像的实施方式进行说明。说明书CN104220898A7/7页90103图20示出在图13所示的第八实施方式中利用标准UTE图像来代替标准CT图像的第十三实施方式。与第八实施方式的不同之处在于,利用标准MR图像和标准UTE图像的数据库40作为数据库,在步骤S524、S534中利用标准UTE图像来代替标准CT图像。其它处与第八实施方式相同,因此省略说明。0104图21示出在图14所示的第九实施方式中利用标准UTE图像来代替标准CT图像的第十四实施方式。与第九实施方式的不同之处在于,利用标准MR图像和标准UTE图像的数据。

33、库40作为数据库,在步骤S524、S544、S554中利用标准UTE图像来代替标准CT图像。其它处与第九实施方式相同,因此省略说明。0105图22示出在图17所示的第十一实施方式中利用标准UTE图像来代替标准CT图像的第十五实施方式。与第十一实施方式的不同之处在于,利用了标准MR图像和标准UTE图像的数据库40作为数据库。其它处与第十一实施方式相同,因此省略说明。0106并且,还能够通过如图23所示的第十六实施方式、图24所示的第十七实施方式那样基于被检者的个人信息DNA、年龄、性别、身高、体重、出生地、居住地、病历等来定制标准图像、射线衰减系数,来进一步提高精度。在图23、24中,S700是。

34、获取个人信息的步骤,S710是生成反映了个人信息的射线衰减系数的步骤,S729是生成反映了个人信息的标准MR图像和标准射线衰减图像的步骤。其它处与图4所示的第一实施方式相同,因此省略说明。0107此外,与标准MR图像一起使用的标准图像并不限定于射线源弱图像,还能够通过按照数据库进行如下改写来应用于其它方式,在标准MR图像和标准CT图像的数据库30时,将S210处替换为图13的S520和S530,在标准MR图像和标准UTE图像的数据库40时,将S210处替换为图20的S524和S534即可。0108另外,上述放射线吸收校正值表和标准图像中的至少一方还能够反复使用被检者本人的放射线吸收校正值、射线。

35、源弱图像、CT图像或者UTE图像。0109此外,在上述实施方式中,本发明适用于脑图像,但本发明的适用对象并不限于此,区域分割数也并不限定于高信号区域、中信号区域、低信号区域这三个区域。0110产业上的可利用性0111本发明能够应用于使PET与MRI相结合的PET/MRI。0112通过参照而将2012年3月28日申请的日本申请号2007074906的说明书、附图以及权利要求书的公开内容整体引用到本说明书中。0113附图标记说明011410被检体;12射线源;14检测器;20标准MR图像和标准射线衰减图像的数据库;30标准MR图像和标准CT图像的数据库;40标准MR图像和标准UTE图像的数据库。说。

36、明书CN104220898A1/22页10图1图2说明书附图CN104220898A102/22页11图3图4说明书附图CN104220898A113/22页12图5说明书附图CN104220898A124/22页13图6说明书附图CN104220898A135/22页14图7说明书附图CN104220898A146/22页15图8说明书附图CN104220898A157/22页16图9说明书附图CN104220898A168/22页17图10说明书附图CN104220898A179/22页18图11说明书附图CN104220898A1810/22页19图12说明书附图CN104220898A。

37、1911/22页20图13说明书附图CN104220898A2012/22页21图14说明书附图CN104220898A2113/22页22图15说明书附图CN104220898A2214/22页23图16说明书附图CN104220898A2315/22页24图17说明书附图CN104220898A2416/22页25图18说明书附图CN104220898A2517/22页26图19说明书附图CN104220898A2618/22页27图20说明书附图CN104220898A2719/22页28图21说明书附图CN104220898A2820/22页29图22说明书附图CN104220898A2921/22页30图23说明书附图CN104220898A3022/22页31图24说明书附图CN104220898A31。

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 物理 > 测量;测试


copyright@ 2017-2020 zhuanlichaxun.net网站版权所有
经营许可证编号:粤ICP备2021068784号-1