CN106023278A - 一种图像重建的方法和装置 - Google Patents
一种图像重建的方法和装置 Download PDFInfo
- Publication number
- CN106023278A CN106023278A CN201610355835.8A CN201610355835A CN106023278A CN 106023278 A CN106023278 A CN 106023278A CN 201610355835 A CN201610355835 A CN 201610355835A CN 106023278 A CN106023278 A CN 106023278A
- Authority
- CN
- China
- Prior art keywords
- crystal
- virtual
- actual
- single photon
- crystals
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 82
- 239000013078 crystal Substances 0.000 claims abstract description 496
- 230000004044 response Effects 0.000 claims abstract description 65
- 230000003111 delayed effect Effects 0.000 claims description 54
- 230000008569 process Effects 0.000 claims description 45
- 238000002600 positron emission tomography Methods 0.000 description 47
- 238000001514 detection method Methods 0.000 description 15
- 238000010586 diagram Methods 0.000 description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 9
- 230000000694 effects Effects 0.000 description 8
- 239000003814 drug Substances 0.000 description 5
- 229940079593 drug Drugs 0.000 description 5
- 230000033001 locomotion Effects 0.000 description 5
- 239000000523 sample Substances 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000004323 axial length Effects 0.000 description 4
- 238000009825 accumulation Methods 0.000 description 3
- 239000008710 crystal-8 Substances 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 238000002591 computed tomography Methods 0.000 description 2
- 238000002595 magnetic resonance imaging Methods 0.000 description 2
- 241001497337 Euscorpius gamma Species 0.000 description 1
- WQZGKKKJIJFFOK-GASJEMHNSA-N Glucose Natural products OC[C@H]1OC(O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-GASJEMHNSA-N 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 239000008103 glucose Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000002503 metabolic effect Effects 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 102000004169 proteins and genes Human genes 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/037—Emission tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/04—Positioning of patients; Tiltable beds or the like
- A61B6/0487—Motor-assisted positioning
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
- A61B6/4241—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector using energy resolving detectors, e.g. photon counting
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10104—Positron emission tomography [PET]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10116—X-ray image
- G06T2207/10128—Scintigraphy
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Nuclear Medicine (AREA)
- Measurement Of Radiation (AREA)
Abstract
本申请提供一种图像重建的方法和装置,该方法包括:获得连续进床扫描模式下的扫描数据;将实际晶体虚拟对应一个以上虚拟晶体,根据所述扫描数据获得所述虚拟晶体所在响应线的延迟随机符合数据;利用实际晶体的晶体接收效率对所述延迟随机符合数据进行去噪声处理,得到去噪声处理后的随机符合数据;利用去噪声处理后的随机符合数据进行图像重建。本申请实施例中,可以在图像重建时使用满足图像重建要求的随机符合数据,实现了对延迟随机符合数据去噪声的效果,从而使得参与图像重建的数据更加准确,图像重建的质量得到提高。
Description
技术领域
本申请涉及医疗技术领域,特别涉及一种图像重建的方法和装置。
背景技术
临床检查影像技术可以包括多种技术,其中一种技术为:在被检体(例如,患者)的体内注入某种物质(一般是生物生命代谢中必须的物质,如:葡萄糖、蛋白质等),并在该物质上标记短寿命的放射性核素(如18F,11C等),然后让被检体躺在扫描床上,并在探测装置中接受扫描。在扫描过程中,被检体中的放射性核素可以产生正电子湮灭事件,该正电子湮灭事件可以产生γ光子,探测装置需要接收该γ光子,并根据γ光子的接收信息(如接收时间、接收位置等)进行分析,以获得来自于同一正电子湮灭事件的γ光子,此为探测到符合事件。探测装置可以根据符合事件的信息重建出被检体的体内图像,该图像可以反映生命代谢活动的情况,进而达到诊断被检体的目的。
由于探测装置的分析存在误差,图像重建所依据的数据可能包括噪声数据,这些噪声数据可能将非符合事件确定为符合事件,因此会影响图像重建的质量。
发明内容
本申请实施例提供一种图像重建的方法,应用于包括多个实际晶体的医疗设备,该方法包括:
获得连续进床扫描模式下的扫描数据;
将实际晶体虚拟对应一个以上虚拟晶体,根据所述扫描数据获得所述虚拟晶体所在响应线的延迟随机符合数据;
利用实际晶体的晶体接收效率对所述延迟随机符合数据进行去噪声处理,得到去噪声处理后的随机符合数据;
利用去噪声处理后的随机符合数据进行图像重建。
本申请实施例提供一种图像重建的装置,应用于包括多个实际晶体的医疗设备,该装置包括:
获得模块,用于获得连续进床扫描模式下的扫描数据;将实际晶体虚拟对应一个以上虚拟晶体,根据所述扫描数据获得所述虚拟晶体所在响应线的延迟随机符合数据;
处理模块,用于利用实际晶体的晶体接收效率对所述延迟随机符合数据进行去噪声处理,得到去噪声处理后的随机符合数据;
重建模块,用于利用去噪声处理后的随机符合数据进行图像重建。
基于上述技术方案,本申请实施例中,针对不满足图像重建要求的延迟随机符合数据,可以使用晶体接收效率对延迟随机符合数据进行去噪声处理,继而使得去噪声处理后的随机符合数据可以满足图像重建要求,从而可以在图像重建时使用满足图像重建要求的随机符合数据,实现对延迟随机符合数据去噪声的效果,从而使得参与图像重建的数据更加准确,图像重建的质量得到提高。
附图说明
为了更加清楚地说明本发明实施例或者现有技术中的技术方案,下面将对本发明实施例或者现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明中记载的一些实施例,对于本领域普通技术人员来讲,还可以根据这些附图获得其它的附图。
图1是本申请一示例性实施例示出的一种弦采集扫描数据场景示意图;
图2(a)和2(b)是本申请一示例性实施例示出的一种即时弦图和延迟弦图示意图;
图3是本申请一示例性实施例示出的一种PET连续进床扫描场景示意图;
图4是本申请一示例性实施例示出的一种PET系统的立体示意图;
图5是本申请一示例性实施例示出的一种虚拟PET系统的示意图;
图6是本申请一示例性实施例示出的一种晶体环的示意图;
图7是本申请一示例性实施例示出的一种单光子累加示意图;
图8是本申请一示例性实施例示出的一种扫描过程示意图;
图9是本申请一示例性实施例示出的一种响应线扫描时间示意图;
图10是本申请一示例性实施例示出的一种响应线扫描晶体示意图;
图11是本申请一示例性实施例示出的两个连续晶体块的示意图;
图12是本申请一示例性实施例示出的一种图像重建的方法流程图;
图13是本申请一示例性实施例示出的一种扫描场景示意图;
图14是本申请一示例性实施例示出的一种医疗设备的结构图;
图15是本申请一示例性实施例示出的一种图像重建的装置的结构示意图。
具体实施方式
在本发明使用的术语仅仅是出于描述特定实施例的目的,而非限制本发明。本发明和权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式,除非上下文清楚地表示其它含义。还应当理解,本文中使用的术语“和/或”是指包含一个或多个相关联的列出项目的任何或所有可能组合。
应当理解,尽管在本发明可能采用术语第一、第二、第三等来描述各种信息,但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如,在不脱离本发明范围的情况下,第一信息也可以被称为第二信息,类似地,第二信息也可以被称为第一信息。取决于语境,此外,所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”。
本申请实施例中提供了一种图像重建的方法,该方法可以用于对被检体进行连续进床扫描模式中的图像重建的过程,且该方法可以应用于医疗设备上。该医疗设备可以如CT(Computed Tomography,电子计算机断层扫描)系统、PET(Positron Emission ComputedTomography,正电子发射型计算机断层显像)系统、DR(Digital Radiography,数字放射显影系统)、MRI(Magnetic Resonance Imaging,磁共振成像系统)等。如下以PET系统为例来描述该方法,可以理解的是,该方法也可以应用于其它系统。其中,PET系统的扫描模式包括分床位扫描模式以及连续进床扫描模式。分床位扫描模式是指:被检体躺在扫描床上,当扫描床运动到某一位置,进行一段时间的扫描,扫描完成后,扫描床运动到下一个位置再进行一段时间的扫描。连续进床扫描模式是指:被检体躺在扫描床上,且扫描床按照一个固定速率运动,并进行持续扫描。
在根据被检体的扫描数据进行图像重建的过程中,扫描数据可能包括噪声数据,而该噪声数据中可能包含随机符合数据。对于随机符合数据,在图像重建中需要使用,而对于其它噪声数据,在图像重建中不需要使用。因此,为了使参与图像重建的数据更加准确,需要得到准确的随机符合数据。
在一个例子中,可以从扫描数据中获得符合数据。如图1所示,如果在一个符合时间窗内,实际晶体A和实际晶体B(后续过程对实际晶体进行详细说明)分别接收到正电子湮灭产生的光子γ1和γ2(从扫描数据中获知该接收情况),被称作发生一个符合事件,而每次发生符合事件就将符合数据加1。为了获得出实际晶体A和实际晶体B的符合数据,假设视野中心O到实际晶体A和实际晶体B连线的距离r是径向坐标,而过O点与两个实际晶体A和实际晶体B连线的垂线OC与x轴的夹角φ为角度坐标,(r,φ)唯一地表示一对实际晶体对应的位置。将所有位于(r,φ)的所发生的符合事件的次数累计的和作为矩阵中(r,φ)处的值,而矩阵中(r,φ)处的值就是实际晶体A和实际晶体B的符合数据。
在实际应用中,符合数据可以为即时符合数据或者随机符合数据。
针对即时符合数据,基于实际晶体A和实际晶体B的扫描数据(接收到正电子湮灭产生的光子γ1和γ2),将实际晶体A和实际晶体B接收到光子γ1和γ2的事件按时间顺序分别排列成两列,如图2(a)所示。上面一列是实际晶体A接收到光子γ1的事件,下面一列是实际晶体B接收到光子γ2的事件,如果实际晶体A接收到光子γ1的事件与实际晶体B接收到光子γ2的事件之间的时间差在一个符合时间窗内,则认为发生一个符合事件,将即时符合数据加1。如图2(a)所示,可以认为一共发生3个符合事件,因此即时符合数据为3。
针对随机符合数据,可以采用延迟符合窗法获得随机符合数据,例如,基于图2(a)的场景,在将实际晶体A和实际晶体B接收到光子γ1和γ2的事件按时间顺序分别排列成两列之后,可以将其中任意一列事件序列延迟一定时间,如图2(b)所示,为将实际晶体A接收到光子γ1的事件延迟一定时间后的示意图,而实际晶体B接收到光子γ2的事件的时间保持不变。基于此,如果实际晶体A接收到光子γ1的事件与实际晶体B接收到光子γ2的事件之间的时间差在一个符合时间窗内,则认为发生一个符合事件,将随机符合数据加1。如图2(b)所示,可以认为一共发生2个符合事件,因此随机符合数据为2。
考虑到采用延迟符合窗法获得的随机符合数据的噪声水平较大,需要对随机符合数据进行去噪声处理。而在连续进床扫描模式中,目前并没有有效的对随机符合数据进行去噪声处理的方式。基于此,在本申请实施例中,将提供一种连续进床扫描模式中,对随机符合数据进行去噪声处理的方式,使得参与重建的数据更加准确,进而提高图像重建的质量。在本申请实施例中,将通过延迟符合窗法获得的随机符合数据称为延迟随机符合数据,在对延迟随机符合数据进行去噪声处理后,可以得到去噪声处理后的随机符合数据。
针对本申请实施例中提供的图像重建的去噪声方法,参见图3所示,图3示例了一个被检体在PET系统中执行连续进床扫描的场景示意图,图中箭头所示为被检体的进床方向(例如,被检体躺在扫描床上,扫描床沿该箭头方向移动),可以是按照一个固定速率运动。图3所示的PET系统,其实是该PET系统的一个截面,该PET系统可以包括很多探测环,探测环的示例可以参见图4。图4所示是PET系统的立体示意图,其中示例了四个探测环21(实际的探测环数量会更多,探测环也可以称为Block(探测器)环),每个探测环由多个Block22拼装而成,而每个Block22包括若干晶体23。其中,多个探测环构成一个内部空间24,被检体就位于该内部空间24中,在这个内部空间发生的正电子湮灭事件所产生的单光子(即γ光子),可以被探测环中的晶体23接收到。图3所示的PET系统,可以是按照图4中的虚箭头所示方向的剖面得到,由图4可以看到,可以将PET系统理解为包括很多晶体环,各晶体环是沿着探测环21相同方向的一个环状晶体串,其中可以包括很多单个的晶体,比如一环100个晶体。
本申请实施例中,虚拟了一个PET系统,实际上,通过将实际晶体虚拟对应一个以上的虚拟晶体,继而形成虚拟PET的系统,而实际晶体与虚拟晶体的对应关系是根据实际晶体位置和床码信息确定的。可以参见图5的虚拟的PET系统示意图,该虚拟的PET系统10包括很多的虚拟晶体,例如,虚拟晶体11的大小可以认为与实际的PET系统的实际晶体一致。为了区分方便,将虚拟的PET系统中的晶体称为虚拟晶体,将实际的PET系统中的晶体称为实际晶体。
进一步的,该虚拟的PET系统10也包括很多探测环,探测环形式上与图4所示的实际的PET系统是相同的,图5是该虚拟的PET系统的一个剖面图。虚拟的PET系统与实际的PET系统的区分在于:虚拟的PET系统10比实际的PET系统要大,主要是探测环的数量会增加(比如,实际的PET系统包括20个探测环,而虚拟的PET系统包括40个探测环,此处数量仅是示例)。虚拟的PET系统10比实际的PET系统要大的原因是:由图3可以看到,实际的PET系统在进床方向(该方向后续可以称为Z轴)上的长度通常比被检体的长度要小,而由图5可以看到,虚拟的PET系统在Z轴方向上的长度比被检体的长度要大,这样,就可以通过虚拟的PET系统10,将连续进床扫描虚拟成被检体在PET系统中进行固定扫描(类似于分床扫描),使得被检体在PET系统中固定不动就可以得到全身各个部位的扫描。
图5所示的虚拟的PET系统10,还有一个特点是,其中的每一个虚拟晶体都固定对应被检体的一个扫描位置。比如,图5中所示的虚拟晶体11对应被检体的扫描位置13,由于虚拟晶体11的长度可以认为与实际的PET系统的实际晶体一致,所以与一个虚拟晶体对应的扫描位置也可以认为是一个实际晶体的长度。比如,假设x是被检体上的一个点位置,(x-δ,x+δ]的范围可以认为是上述扫描位置13,δ为实际晶体大小的一半。对于扫描位置15,可以是与虚拟晶体14对应,可以看到,每一个虚拟晶体与被检体12的一个扫描位置对应。
进一步的,上述的每一个虚拟晶体与一个扫描位置对应,其中的“对应”可以这样理解:假设被检体12在虚拟的PET系统10中固定不动的扫描,则被检体12的上述扫描位置13与虚拟晶体11在位置上正对,可以类似于传统的分床扫描中的扫描床固定不动时的实际晶体与扫描位置的一一对应。
此外,如前面提到的,图5只是示意了虚拟的PET系统的一个剖面,其实该虚拟的PET系统也是一个环形的探测装置,类似于图4所示的立体系统,当被检体放置于虚拟的PET系统的内部空间时,对应某一个扫描位置的虚拟晶体的数量很多,其实是一个晶体环。可以参见图6的示例,以扫描位置13为例,虚拟晶体11、虚拟晶体41、虚拟晶体42等一个环的很多虚拟晶体,都是与该扫描位置对应的虚拟晶体,也就是说,与某个扫描位置对应的虚拟晶体是一个环的多个虚拟晶体,这些虚拟晶体只是相对于扫描位置的角度不同。可以想象,当被检体躺在扫描床上,位于图4所示的内部空间24中时,其实就是位于一个环状探测装置的内部,被检体的任一个位置都对应着一个环的虚拟晶体。
经过上述描述可以看到,该虚拟的PET系统10,即为虚拟成一个类似于分床扫描的PET系统,尽管被检体在PET系统中以连续进床方式扫描,扫描床是不断运动的,但是,通过将PET系统10假想成一个相当于扫描床不动时扫描的PET系统,相当于将被检体放置在图5的虚拟的PET系统10中进行固定不动的扫描,虚拟的PET系统10中的各个虚拟晶体固定对应被检体的一个扫描位置。
本申请实施例中,如果要利用图5所示的虚拟的PET系统10,对响应线上的延迟随机符合数据进行去噪声处理,则可能执行如下过程:获得该响应线上的虚拟晶体的单光子计数率;获得该虚拟晶体在响应线上对应的各个实际晶体的晶体接收效率;利用该晶体接收效率对延迟随机符合数据进行去噪声处理。
如下提供一种获得虚拟晶体的单光子计数率的例子,在具体实施方式中并不局限于下面列举的获得方式。为了获得虚拟晶体的单光子计数率,首先获得虚拟晶体的单光子计数,然后,根据得到的虚拟晶体的单光子计数,以及虚拟晶体对应的扫描位置的扫描时间,获得虚拟晶体的单光子计数率,比如,通过虚拟晶体的单光子计数除以该扫描时间得到虚拟晶体的单光子计数率。
在一个例子中,可以统计虚拟的PET系统中的每一个虚拟晶体在被检体的扫描过程中,所获取到的单光子计数,而且该单光子计数就是在整个扫描过程中,该虚拟晶体接收到的单光子(即正电子湮灭事件中产生的γ光子)的数量,该虚拟晶体接收到的单光子的数量就是虚拟晶体的单光子计数。
虚拟晶体的单光子计数的获得:
每一个虚拟晶体的单光子计数都是实际的PET系统中的多个实际晶体的单光子计数的累加,具体是将哪些实际晶体接收到的哪部分单光子计数,统计到与这些实际晶体对应的某个虚拟晶体,可以结合图5和图7来说明。
如图5所示,仍以被检体12的其中一个扫描位置为例,比如,以扫描位置13为例,该扫描位置13对应的是虚拟晶体11。在连续进床扫描模式中,被检体将随着扫描床一起移动,而被检体上的扫描位置13也将随之一起移动,图7中示出了扫描位置13被扫描过程的两个端点,从13-1的位置开始接受扫描,并沿着箭头51的方向(这也是扫描床的进床方向)随被检体一起运动,直至运动到13-2的位置停止接受扫描(比如,可以是移出了PET系统)。此外,图7中还示出了这两个端点位置所对应的晶体环,其中的虚线52表示这两个晶体环之间是连续的,还有很多晶体环,可以想象,这是一个立体的实际的PET系统,只是在图7中仅示出了该PET系统两个端侧的两个晶体环。
请继续参见图7所示,实际晶体53、实际晶体54、实际晶体55直至实际晶体58(在图7中仅示例了一部分实际晶体,实际数量可能更多),组成了一个“晶体串”,该晶体串的方向与箭头51所示的行进方向平行,并且,在扫描位置13从13-1的位置到13-2的位置这个过程的扫描中,上述晶体串中的每一个实际晶体都会有一个与扫描位置13“对应”的时间。比如,扫描位置13位于图7所示的最左端的晶体环中时,则实际晶体53处于与扫描位置13“对应”的位置,当扫描位置13位于图7所示的最右端的晶体环中时,则实际晶体58处于与扫描位置13“对应”的位置,而且,中间的各个晶体,比如实际晶体54、实际晶体55等,都会有一个与扫描位置13“对应”的位置。
图8示例了某个实际晶体与扫描位置13位置对应的运动过程示意图,结合图8所示,以实际的PET系统的其中一个实际晶体55为例,图8示意了连续进床扫描过程中的五个时刻,分别为T1、T2、T3、T4和T5,从T1至T5的顺序,被检体12沿着箭头所示的进床方向在移动,可以是按照一定速率匀速进床。其中,实际晶体55的位置固定,随着被检体12的进床移动,实际晶体55的中心线S所对应的被检体12上的位置也在不断移动。在T3时刻,扫描位置13的左侧“x-δ”(x是扫描位置13的中心线,δ是实际晶体长度的一半)到达实际晶体55的中心线S,在T4时刻,扫描位置13的右侧“x+δ”到达实际晶体55的中心线S,在T3至T4的时间段内,可以认为扫描位置13与实际晶体55对应。
按照上述的原理,被检体12在连续进床的过程中,对于图7中所示的实际晶体53至实际晶体58之间的这个晶体串来说,每个实际晶体都会有一个类似图8所示的对应扫描位置13的时间段,在该时间段内该实际晶体接收到的单光子数,即属于扫描位置13对应的虚拟晶体11接收的单光子计数的一部分。即,可以将实际晶体53至实际晶体58之间的这个晶体串中的每个实际晶体在对应扫描位置13时接收到的光子数累加,得到图6中所示的对应扫描位置13的虚拟晶体11的单光子计数。其中,虚拟晶体11与上述的晶体串的关系,由图6和图7可以看到,如果假设扫描位置13处于晶体环的环心,那么,虚拟晶体11与上述晶体串中的每个实际晶体在这个晶体环中的位置是相同的,即实际晶体与所述虚拟晶体在对应所述扫描位置时,与扫描位置的相对位置关系相同(如果以晶体环360度来看,则两者在环上的角度位置相同)。
同理,在图7中,另一个晶体串包括实际晶体60、实际晶体61、实际晶体62直至实际晶体64,将这些实际晶体在对应扫描位置13时接收到的单光子数,可以均累加到图6所示的虚拟晶体42中。其它虚拟晶体的单光子计数的获得不再详述,原理与上述相同,都是由与该虚拟晶体位于相同位置的一个晶体串的各个实际晶体的单光子数的累加。如果将一个实际晶体在对应某个扫描位置时接收到的单光子计数称为一个单光子计数分量,那么实际的PET系统的一个晶体串的多个实际晶体,将获取到对应该扫描位置的多个单光子计数分量,将这些单光子计数分量叠加,即得到对应该扫描位置的虚拟晶体的单光子计数。
通过上述方法,可以对虚拟的PET系统中的每个虚拟晶体,得到该虚拟晶体在扫描过程中的单光子计数,即对于某个虚拟晶体来说,确定该虚拟晶体对应的被检体的扫描位置,还确定该虚拟晶体所对应的实际PET系统中的某个晶体串,将该扫描位置在整个扫描过程中与上述晶体串中的各个实际晶体对应时的单光子计数累加到虚拟晶体中,继而可以得到该虚拟晶体的单光子计数。
在一个例子中,为了实现上述的将实际晶体得到的单光子计数向对应的虚拟晶体的累加,可以按照如下方式处理:实际晶体在接收到γ光子以后,可以将该γ光子的接收信息传输至后端的数据处理系统,该接收信息中可以包括γ光子的接收时间,数据处理系统可以根据被检体的进床运动信息以及该接收时间,确定所述的γ光子对应的被检体的扫描位置,并将该光子的计数累加到对应上述扫描位置的虚拟晶体中。例如,被检体的进床运动信息可以包括被检体的初始位置信息和被检体的连续进床速度。比如,在扫描开始时,被检体中的某个扫描位置对应的是实际晶体C1,被检体的连续进床速度是每秒钟s1的距离,那么根据γ光子的接收时间以及上述连续进床速度,就可以知道在所述“接收时间”的时刻,实际晶体C1对应的被检体中的扫描位置。
经过上述的描述获知,可以构造出一个虚拟的PET系统10,并且将连续进床扫描模式下的扫描过程,虚拟为被检体在该虚拟的PET系统10中,固定扫描一段时间的扫描过程,而且还可以通过上述的方法,获得该扫描过程中,虚拟的PET系统10中的各个虚拟晶体所获得的单光子计数。
如下提供一种获得实际晶体的晶体接收效率的例子,在具体实施方式中并不局限于下面列举的获得方式。其中,当在单位时间内有N个单光子发射到实际晶体上时,如果被识别到M个单光子,则实际晶体的晶体接收效率η=M/N。
实际晶体的晶体接收效率的获得:
在一个例子中,PET系统中的每个实际晶体的晶体接收效率可以使用一个近似值,在PET系统中心放置一个圆柱型水模,该圆柱型水模的药物活度保证PET系统内的剂量和扫描待检体时接近,在扫描床不动的情况下扫描圆柱型水模,获得所有实际晶体的单光子计数率,并使用如下公式获得实际晶体的晶体接收效率:其中,i表示实际晶体的晶体号,Si表示实际晶体i上的单光子计数率,表示所有实际晶体的单光子计数率的均值。基于上述分析可以看出,在同一剂量下,可以获得实际晶体的单光子计数率与所有实际晶体的单光子计数率的均值之间的比值,并将所述比值作为所述实际晶体的晶体接收效率。
但实际上,在不同的药物活度下,实际晶体的晶体接收效率不同,为了更准确的获得实际晶体的晶体接收效率,在另一个例子中,可以分别获得不同剂量(不同剂量对应不同的药物活度)下的每个实际晶体的晶体接收效率,并根据不同剂量下的数据建立一个实际晶体上的单光子计数率与晶体接收效率的函数关系,以使得在得到实际晶体的单光子计数率后,可以根据函数关系获得对应的晶体接收效率,且晶体接收效率是根据不同剂量的测试确定,较为准确。基于上述分析可以看出,在不同剂量下,可以根据实际晶体的单光子计数率,以及单光子计数率与晶体接收效率之间的函数关系,得到所述实际晶体的晶体接收效率,所述函数关系是根据不同剂量下获得的晶体接收效率得到。
例如,为了获得某一药物活度下的PET系统中的每个实际晶体的晶体接收效率,可以在PET系统中心放置轴向长度超过PET系统的轴长的一个圆柱型水模,该圆柱型水模的药物剂量为D0,在扫描床不动的情况下扫描圆柱型水模,获得所有实际晶体的单光子计数率,并使用如下公式获得实际晶体的晶体接收效率:其中,i表示实际晶体的晶体号,Si(0)表示实际晶体i上的单光子计数率,表示所有实际晶体的单光子计数率的均值。
此外,为了获得不同剂量下的PET系统中的每个实际晶体的晶体接收效率,如上所述,在PET系统中心放置轴向长度超过PET系统的轴长的一个圆柱型水模,该圆柱型水模的药物剂量为Dt,在扫描床不动的情况下扫描圆柱型水模,获得所有实际晶体的单光子计数率,并使用如下公式获得实际晶体的晶体接收效率:通过对不同的Dt做如上测试,获得ηi(Si(t))。
基于此,可以根据不同剂量的测试结果建立函数关系ηi(S(t)),该函数关系用于表示实际晶体上的单光子计数率S(t)与实际晶体的晶体接收效率之间的关系,使得只要知道实际晶体上的单光子计数率S(t),即可获得晶体接收效率。
可以按照上述的方法,获得实际的PET系统中的其中一个实际晶体的晶体接收效率。此外,当获得出一个实际晶体的晶体接收效率之后,可以很方便快捷的获取PET系统中的其它实际晶体的晶体接收效率。具体的,在临床活度范围内,不同活度下的每个实际晶体的单光子计数率不同,而且随着活度增大,所有实际晶体的单光子计数率都是增大的。但是,在一个固定剂量Dt下,不同实际晶体之间的单光子计数率Si(t)之间的比例关系是固定的,因此,只要知道一个实际晶体上的单光子计数率,就可以获得所有实际晶体上的单光子计数率,而且知道几个实际晶体上的单光子计数率的均值,也可以获得所有实际晶体上的单光子计数率。即只要知道一个实际晶体上的单光子计数率或者几个实际晶体上的单光子计数率的均值,就可以获得所有实际晶体上的单光子计数率,再根据上述的函数关系,就可以得到对应的实际晶体的晶体接收效率。
例如,实际晶体1与实际晶体2的比例关系是1:2,实际晶体2与实际晶体3的比例关系是1:2,实际晶体3与实际晶体4的比例关系是1:2。基于此,假设知道实际晶体2上的单光子计数率为A,则可以获得实际晶体1上的单光子计数率为A/2,获得实际晶体3上的单光子计数率为A*2,获得实际晶体4上的单光子计数率为A*4。假设上述4个实际晶体上的单光子计数率的均值为B,则4个实际晶体上的单光子计数率的总和是4B,因此,可以获得实际晶体1上的单光子计数率为4B/15,获得实际晶体2上的单光子计数率为8B/15,获得实际晶体3上的单光子计数率为16B/15,获得实际晶体4上的单光子计数率为32B/15。
为了实现上述利用实际晶体上的单光子计数率获得晶体接收效率的过程,首先要获得实际晶体上的单光子计数率。在一个例子中,对于每个实际晶体来说,可以统计出该实际晶体接收到的单光子计数,并根据单光子计数以及单光子计数的接收时间,获得所述实际晶体的单光子计数率,如通过实际晶体接收到的单光子计数除以该接收时间,得到所述实际晶体的单光子计数率。在另一个例子中,可以获得对应一个扫描位置的多个单光子计数分量,根据多个单光子计数分量,得到与所述扫描位置对应的虚拟晶体的单光子计数率,利用所述虚拟晶体的单光子计数率获得所述虚拟晶体对应的实际晶体的单光子计数率。
其中,针对第一种方式,对于每个实际晶体来说,可以统计出该实际晶体在扫描过程中所获取到的单光子计数,该单光子计数就是该实际晶体接收到的单光子(即正电子湮灭事件中产生的γ光子)的数量,而该实际晶体的单光子计数率可以通过实际晶体的单光子计数除以单光子计数的接收时间得到。
其中,针对第二种方式,针对利用虚拟晶体的单光子计数率获得实际晶体的单光子计数率的过程,在上面的内容中已经指出,将实际晶体53至实际晶体58之间的这个晶体串中的每个实际晶体在对应扫描位置13时接收到的光子数累加,得到虚拟晶体11的单光子计数,因此一个虚拟晶体的单光子计数是多个实际晶体的单光子计数的和,且一个虚拟晶体的单光子计数率是多个实际晶体的单光子计数率的和。假设虚拟晶体的单光子计数率为C,则说明实际晶体53至实际晶体58这多个实际晶体上的单光子计数率的总和是C,基于实际晶体53至实际晶体58这多个实际晶体的比例关系,就可获得实际晶体的单光子计数率。
针对延迟随机符合数据的去噪声处理过程,由于延迟随机符合数据的噪声水平较大,不符合图像重建的质量要求,因此需要对延迟随机符合数据进行去噪声处理,以使去噪声处理后的随机符合数据符合图像重建的质量要求。在介绍对延迟随机符合数据进行去噪声处理之前,首先介绍符合图像重建的质量要求的去噪声处理后的随机符合数据,这一随机符合数据就是本申请中需要得到的随机符合数据,其是通过对延迟随机符合数据进行去噪声处理后得到。
其中,该去噪声处理后的随机符合数据是响应线上的随机符合数据,在一个例子中,响应线上的随机符合数据=响应线的随机符合计数率*扫描时间。假设某条响应线对应虚拟的PET系统中的虚拟晶体i和虚拟晶体j,在上述公式中,随机符合计数率即该虚拟晶体i和虚拟晶体j对应的响应线上单位时间内得到的随机符合数据,而扫描时间是这条响应线在实际的PET系统中的扫描时间,可以称为响应线的扫描时间。例如,可以结合图9来说明上述“响应线的扫描时间”,在图9中,虚拟的PET系统中的与被检体中的某个扫描位置x对应的响应线Y,该响应线Y在实际的PET系统中,由开始扫描到结束扫描的时间段,即图9中的被检体由时间Ts运动至时间Te的时间段。如下提供一种获得响应线的随机符合计数率的例子,在具体方式中并不局限于下面列举的获得方式。
响应线的随机符合计数率的获得:在一个例子中,对于虚拟晶体i和虚拟晶体j对应的响应线上的随机符合计数率的获得,可以按照如下方式:对于每个虚拟晶体,根据虚拟晶体的单光子计数率,以及所述虚拟晶体在响应线上对应的各个实际晶体的晶体接收效率,得到所述响应线的随机符合计数率。
如下的公式(1),即为虚拟晶体i和虚拟晶体j对应的响应线的随机符合计数率的获得公式,下面将说明对于该公式的推导过程以及该公式的含义,以对上面提到的响应线上的随机符合计数率的获得方法进行说明。
其中,响应线在图9所示的实际的PET系统中扫描时,所经过的实际晶体的对数是m对,结合图10所示,图10中示例了响应线的扫描过程,实线的响应线表示该响应线的起始扫描点,两条虚线的响应线表示在扫描过程中的该响应线的两个位置,整个扫描过程经过了m对实际晶体对。比如,按照响应线的方向,实际晶体81和实际晶体82是一个晶体对,实际晶体83和实际晶体84是一个晶体对,共有m对实际晶体对。并且上述响应线的扫描过程中经过的实际晶体对中的实际晶体,可称为虚拟晶体在响应线上对应的各个实际晶体,比如虚拟晶体i在响应线上对应的各个实际晶体包括实际晶体81、实际晶体83等。
在公式(1)中,i和j是响应线上的两个虚拟晶体的晶体号,x1、x2、…xm是虚拟晶体i在响应线上对应的各个实际晶体,y1、y2、…ym是虚拟晶体j在响应线上对应的各个实际晶体,ηix1、ηix2、ηixm是实际晶体x1、x2、…xm对应的晶体接收效率,ηjy1、ηjy2、ηjym是实际晶体y1、y2、…ym对应的晶体接收效率,是虚拟的PET系统中由被检体发射到虚拟晶体i的单光子计数率,是虚拟的PET系统由被检体发射到虚拟晶体j的单光子计数率,2τ是符合时间窗。
在上面的过程中已经介绍,每个虚拟晶体可以对应多个实际晶体,如虚拟晶体1对应实际晶体1、实际晶体2和实际晶体3,虚拟晶体6对应实际晶体7、实际晶体8和实际晶体9,基于此,假设在第一个扫描位置,虚拟晶体1对应实际晶体1,虚拟晶体6对应实际晶体7,在第二个扫描位置,虚拟晶体1对应实际晶体2,虚拟晶体6对应实际晶体8,在第三个扫描位置,虚拟晶体1对应实际晶体3,虚拟晶体6对应实际晶体9。针对公式(1),当i为1,j为6时,则实际晶体对的数量m为3,x1为实际晶体1、x2为实际晶体2、x3为实际晶体3、y1为实际晶体7、y2为实际晶体8、y3为实际晶体9。
在公式(1)中,可以另αij是虚拟晶体i和虚拟晶体j的晶体对接收效率,基于此,公式(1)可以简化为在得到响应线的随机符合计数率rij之后,由于响应线上的随机符合数据为随机符合计数率*扫描时间,因此随机符合数据为T为对应的扫描时间。
基于上述分析,随机符合数据就是本申请中需要得到的去噪声处理后的随机符合数据,即对延迟随机符合数据进行去噪声处理的结果是该随机符合数据在本申请实施例中,可以另Rij为该随机符合数据
在一个例子中,可以获得多个延迟随机符合数据,对于延迟随机符合数据的获得方式,已经在上面的过程中介绍,即采用延迟符合窗法获得延迟随机符合数据,在此不再重复赘述,可以直接使用这多个延迟随机符合数据进行处理。延迟随机符合数据Ril是虚拟晶体i和虚拟晶体l对应的延迟随机符合数据,l∈B,且B是包含虚拟晶体j的一组虚拟晶体的集合。延迟随机符合数据Rjk是虚拟晶体j和虚拟晶体k对应的延迟随机符合数据,k∈A,且A是包含虚拟晶体i的一组虚拟晶体的集合。Rlk是虚拟晶体l和虚拟晶体k对应的延迟随机符合数据。其中,如图11所示,A和B可以是探测环中相对的两个连续晶体块。
为了确定出虚拟晶体的集合A和虚拟晶体的集合B,在一个例子中,可以给出集合A内的虚拟晶体的数量,并给出集合B内的虚拟晶体的数量。通常情况下,集合A内的虚拟晶体的数量与集合B内的虚拟晶体的数量相同,假设数量为H个,如图11所示,该数量H为10个。而且,集合A内的虚拟晶体是连续的,集合B内的虚拟晶体也是连续的。在此基础上,只要集合A内包含虚拟晶体i,集合B内包含虚拟晶体j即可。因此,集合A内包含虚拟晶体i在内的连续的H个虚拟晶体,集合B内包含虚拟晶体j在内的连续的H个虚拟晶体。
在一个例子中,虚拟晶体i可以位于连续的H个虚拟晶体的中间,例如,如果H/2为整数,则虚拟晶体i是集合A中的第H/2个虚拟晶体,如果H/2为小数,则虚拟晶体i是集合A中的对H/2向上取整个虚拟晶体。同理,虚拟晶体j是集合B中的第H/2个虚拟晶体,或者对H/2向上取整个虚拟晶体。
当然,虚拟晶体i和虚拟晶体j也可以是其它位置的虚拟晶体,如虚拟晶体i是集合A中的第一个虚拟晶体,或者最后一个虚拟晶体等,虚拟晶体j是集合B中的第一个虚拟晶体,或者最后一个虚拟晶体等,本申请实施例中对于虚拟晶体i和虚拟晶体j的位置不再说明,只要得到集合A和集合B即可。
需要说明的是,此处使用的是两个虚拟晶体对应的延迟随机符合数据,而上述延迟符合窗法获得的是两个实际晶体对应的延迟随机符合数据,因此,需要利用两个实际晶体对应的延迟随机符合数据,获得两个虚拟晶体对应的延迟随机符合数据。在一个例子中,对于这两个虚拟晶体,可以获知每个虚拟晶体对应的所有实际晶体,并利用这些实际晶体对应的延迟随机符合数据,获得这两个虚拟晶体对应的延迟随机符合数据。例如,虚拟晶体1对应实际晶体1、实际晶体2和实际晶体3,虚拟晶体2对应实际晶体4、实际晶体5和实际晶体6,则虚拟晶体1和虚拟晶体2对应的延迟随机符合数据可以为:实际晶体1与实际晶体4对应的延迟随机符合数据,加上实际晶体2与实际晶体5对应的延迟随机符合数据,加上实际晶体3与实际晶体6对应的延迟随机符合数据。
例如,当i为11,j为16时,B可以是虚拟晶体15、虚拟晶体16、虚拟晶体17组成的一组虚拟晶体的集合,A可以是虚拟晶体10、虚拟晶体11、虚拟晶体12组成的一组虚拟晶体的集合。基于此,k的取值可以为10、11、12,l的取值可以为15、16、17。当k的取值为12,l的取值为17时,则Ril是虚拟晶体11和虚拟晶体17对应的延迟随机符合数据,Rjk是虚拟晶体16和虚拟晶体12对应的延迟随机符合数据,Rlk是虚拟晶体17和虚拟晶体12对应的延迟随机符合数据。对于其它k和l的取值组合,均有对应的延迟随机符合数据Ril、Rjk、Rlk,在此不再重复赘述。
在一个例子中,延迟随机符合数据Ril可以通过如下公式获得:延迟随机符合数据Rjk可以通过如下公式获得:延迟随机符合数据Rlk可以通过如下公式获得:进一步的,由于l∈B,k∈A,因此,RiB可以通过如下公式获得:Σl∈BRil,因此,同理,RjA可以通过如下公式获得:Σk∈ARjk,因此,
基于此,通过整理公式(2),得出由于Rij为并且Rlk为因此,可以通过整理公式(3),得出进一步的,通过整理公式(4),就可以得出之后,将上述RiB和RjA代入公式(5)中,可以得出而Rij就是去噪声处理后的随机符合数据。
在上述公式中,αil是虚拟晶体i和虚拟晶体l的晶体对接收效率,αjk是虚拟晶体j和虚拟晶体k的晶体对接收效率,αlk是虚拟晶体l和虚拟晶体k的晶体对接收效率。例如,当i为11,j为16,k的取值为12,l的取值为17时,则αil是虚拟晶体11和虚拟晶体17的晶体对接收效率,αjk是虚拟晶体16和虚拟晶体12的晶体对接收效率,αlk是虚拟晶体17和虚拟晶体12的晶体对接收效率。
其中,αil的获得公式与上述αij的获得公式类似,只是将αij中的j替换成l,如αjk的获得公式与上述αij的获得公式类似,只是将αij中的j替换成k,将αij中的i替换成j,αlk的获得公式与上述αij的获得公式类似,只是将αij中的j替换成k,将αij中的i替换成l,在此不再赘述。
由上述描述,本申请实施例的图像重建的方法,应用于包括多个实际晶体的医疗设备,并可以包括图12中所示的如下处理步骤进行处理:
在步骤1201中,获得连续进床扫描模式下的扫描数据。
其中,在连续进床扫描模式下的扫描过程中,被检体中的放射性核素可以产生正电子湮灭事件,该正电子湮灭事件可以产生γ光子,实际晶体需要接收该γ光子,而这些γ光子的接收信息(如接收时间、接收位置等)就是扫描数据。当然,扫描数据还可以包括其它信息,所有在连续进床扫描模式下的扫描过程中,得到的数据均可以认为是扫描数据,在此不再详加赘述。
在步骤1202中,将实际晶体虚拟对应一个以上虚拟晶体,根据扫描数据获得虚拟晶体所在响应线的延迟随机符合数据。
其中,在上面的过程中已经介绍到,可以通过将实际晶体虚拟对应一个以上的虚拟晶体,继而形成虚拟PET的系统,而实际晶体与虚拟晶体的对应关系是根据实际晶体位置和床码信息确定的,在此不再重复赘述。
其中,在上面的过程中已经介绍到,基于扫描数据,可以采用延迟符合窗法获得两个实际晶体对应的延迟随机符合数据。此外,针对两个虚拟晶体来说,可以根据这两个虚拟晶体对应的所有实际晶体的延迟随机符合数据,获得这两个虚拟晶体对应的延迟随机符合数据,在此不再重复赘述。
在步骤1203中,利用实际晶体的晶体接收效率对延迟随机符合数据进行去噪声处理,得到去噪声处理后的随机符合数据。
在一个例子中,在步骤1203之前,还可以获得实际晶体的单光子计数率,并利用实际晶体的单光子计数率得到该实际晶体的晶体接收效率。
所述获得实际晶体的单光子计数率的过程,可以包括如下方式:
方式一、统计所述实际晶体接收到的单光子计数,并根据所述单光子计数以及所述单光子计数的接收时间,得到所述实际晶体的单光子计数率。
方式二、获得对应一个扫描位置的多个单光子计数分量,根据所述多个单光子计数分量,得到与所述扫描位置对应的虚拟晶体的单光子计数率;利用所述虚拟晶体的单光子计数率,得到所述虚拟晶体对应的实际晶体的单光子计数率;其中,一个单光子计数分量是所述扫描位置运动到对应的一个实际晶体时,所述实际晶体接收到的单光子计数,且所述实际晶体与所述虚拟晶体在对应所述扫描位置时与所述扫描位置的相对位置关系相同。
其中,虚拟晶体固定对应被检体的某个扫描位置,并且向该虚拟晶体中累加单光子计数的各个实际晶体,在对应所述扫描位置时,实际晶体和虚拟晶体与扫描位置的相对位置关系相同,且每一个实际晶体累加的计数是该扫描位置运动到对应该晶体时得到,可以称为一个单光子计数分量。
针对方式二、在根据所述多个单光子计数分量,得到与所述扫描位置对应的虚拟晶体的单光子计数率的过程中,根据被检体的进床运动信息、以及接收到的单光子的接收时间,确定所述单光子对应的被检体的扫描位置,并将所述多个单光子计数分量累加到与所述扫描位置对应的虚拟晶体中,得到所述虚拟晶体的单光子计数;根据所述虚拟晶体的单光子计数、以及所述虚拟晶体对应的扫描位置在扫描过程中的扫描时间,得到所述虚拟晶体的单光子计数率。
在一个例子中,利用实际晶体的单光子计数率得到实际晶体的晶体接收效率的过程,具体可以包括如下方式:方式一、在不同剂量下,根据实际晶体的单光子计数率,以及单光子计数率与晶体接收效率之间的函数关系,得到所述实际晶体的晶体接收效率,所述函数关系是根据不同剂量下获得的晶体接收效率得到。方式二、在同一剂量下,获得实际晶体的单光子计数率与所有实际晶体的单光子计数率的均值之间的比值,并将所述比值作为所述晶体接收效率。
在一个例子中,利用晶体接收效率对延迟随机符合数据进行去噪声处理,得到去噪声处理后的随机符合数据的过程,具体可以包括如下方式:确定响应线上的两个虚拟晶体、以及所述两个虚拟晶体在响应线上对应的各个实际晶体,并获取各个实际晶体的晶体接收效率;根据所述各个实际晶体的晶体接收效率,得到所述两个虚拟晶体的晶体对接收效率,并利用所述晶体对接收效率对所述延迟随机符合数据进行去噪声处理,得到去噪声处理后的随机符合数据。
本发明实施例中,根据所述各个实际晶体的晶体接收效率,得到所述两个虚拟晶体的晶体对接收效率的过程,具体可以包括:针对所述响应线上的虚拟晶体i和虚拟晶体j,根据所述虚拟晶体i在所述响应线上对应的m个实际晶体的晶体接收效率、所述虚拟晶体j在所述响应线上对应的m个实际晶体的晶体接收效率,得到所述虚拟晶体i和所述虚拟晶体j的晶体对接收效率。
在一个例子中,可以利用如下公式获得所述两个虚拟晶体的晶体对接收效率:该公式的内容参见上面的实施例。
本发明实施例中,利用所述晶体对接收效率对所述延迟随机符合数据进行去噪声处理,得到去噪声处理后的随机符合数据的过程,具体可以包括:根据包含虚拟晶体j的一组虚拟晶体的集合A、包含虚拟晶体i的一组虚拟晶体的集合B、虚拟晶体i和属于集合B的虚拟晶体l的晶体对接收效率、虚拟晶体j和属于集合A的虚拟晶体k的晶体对接收效率、虚拟晶体l和虚拟晶体k的晶体对接收效率,对虚拟晶体i和虚拟晶体l对应的延迟随机符合数据、虚拟晶体j和虚拟晶体k对应的延迟随机符合数据、虚拟晶体l和虚拟晶体k对应的延迟随机符合数据,进行去噪声处理,得到去噪声处理后的虚拟晶体i和虚拟晶体j的随机符合数据。
在一个例子中,可以利用如下公式对所述延迟随机符合数据进行去噪声处理:该公式的内容参见上面的实施例。
在步骤1204中,利用去噪声处理后的随机符合数据进行图像重建。
例如,可以在图像重建的数据中,对于每条响应线都使用该响应线上的去噪声处理后的随机符合数据,从而使得图像重建的数据更加准确。
基于上述技术方案,本申请实施例中,针对不满足图像重建要求的延迟随机符合数据,可以使用晶体接收效率对延迟随机符合数据进行去噪声处理,继而使得去噪声处理后的随机符合数据可以满足图像重建要求,从而在图像重建时使用满足图像重建要求的随机符合数据,实现了对延迟随机符合数据去噪声的效果,从而使得参与图像重建的数据更加准确,图像重建的质量得到提高。
图13示例了一个对被检体进行扫描的应用场景,例如,可以是通过连续进床扫描模式对被检体进行扫描。工作人员1301可以通过控制台1302操作探测装置(该装置包括多个实际晶体)1303,以对扫描床1304上的被检体1305进行连续进床扫描。在扫描过程中,探测装置1303中的各个实际晶体可以将接收到的单光子计数传送到后端的处理装置1306,该处理装置1306可以是用于对探测装置接收到的数据进行处理并进行图像重建的装置。该处理装置1306就可以执行上面方法实施例中描述的图像重建的方法。其中,控制台1302、探测装置1303、处理装置1306就可以是医疗设备内部的几个功能模块。
本申请实施例还提供了一种图像重建的装置,该装置可以应用于包括多个实际晶体的医疗设备,图14是根据一示例性实施例示出的一种医疗设备1400的框图。参照图14,医疗设备1400可以包括处理组件1401,其进一步包括一个或者多个处理器,以及由存储器1402所代表的存储器资源,用于存储可由处理组件1401执行的指令,例如应用程序。存储器1402中存储的应用程序可以包括一个或者一个以上的每一个对应于一组指令的模块。
本申请的实施例中,图像重建的装置可以位于所述的存储器1402中,并且处理组件1401可以通过该图像重建的装置执行本申请实施例的图像重建的方法,以对被检体进行连续进床扫描模式时的图像重建。
医疗设备1400还可以包括一个电源组件1403,该电源组件1403被配置为执行医疗设备1400的电源管理。一个有线或者无线网络接口1404被配置为将医疗设备1400连接到网络,以及一个输入输出(I/O)接口1405。
本申请实施例的图像重建的装置,该装置用于对被检体进行连续进床扫描模式中的图像重建;该装置可以从逻辑上划分为多个模块,比如参见图15所示,该装置可以包括:获得模块1501、处理模块1502和重建模块1503;其中:
获得模块1501,用于获得连续进床扫描模式下的扫描数据;将实际晶体虚拟对应一个以上虚拟晶体,根据所述扫描数据获得所述虚拟晶体所在响应线的延迟随机符合数据;
处理模块1502,用于利用实际晶体的晶体接收效率对所述延迟随机符合数据进行去噪声处理,得到去噪声处理后的随机符合数据;
重建模块1503,用于利用去噪声处理后的随机符合数据进行图像重建。
在一个例子中,所述获得模块1501,还用于获得实际晶体的单光子计数率,利用所述实际晶体的单光子计数率得到所述实际晶体的晶体接收效率。
在一个例子中,所述获得模块1501,具体用于在获得实际晶体的单光子计数率的过程中,统计所述实际晶体接收到的单光子计数,并根据所述单光子计数以及所述单光子计数的接收时间,得到所述实际晶体的单光子计数率;或者,
获得对应一个扫描位置的多个单光子计数分量,根据所述多个单光子计数分量,得到与所述扫描位置对应的虚拟晶体的单光子计数率;利用所述虚拟晶体的单光子计数率,得到所述虚拟晶体对应的实际晶体的单光子计数率;
其中,一个单光子计数分量是所述扫描位置运动到对应的一个实际晶体时,所述实际晶体接收到的单光子计数,且所述实际晶体与所述虚拟晶体在对应所述扫描位置时与所述扫描位置的相对位置关系相同。
在一个例子中,所述获得模块1501,具体用于在利用所述实际晶体的单光子计数率得到所述实际晶体的晶体接收效率的过程中,在不同剂量下,根据所述实际晶体的单光子计数率,以及单光子计数率与晶体接收效率之间的函数关系,得到所述实际晶体的晶体接收效率,所述函数关系是根据不同剂量下计算的晶体接收效率得到;或者,在同一剂量下,获得所述实际晶体的单光子计数率与所有实际晶体的单光子计数率的均值之间的比值,并将所述比值作为所述晶体接收效率。
在一个例子中,所述处理模块1502,具体用于在利用实际晶体的晶体接收效率对所述延迟随机符合数据进行去噪声处理,得到去噪声处理后的随机符合数据的过程中,确定响应线上的两个虚拟晶体、以及所述两个虚拟晶体在所述响应线上对应的各个实际晶体,并获取所述各个实际晶体的晶体接收效率;根据所述各个实际晶体的晶体接收效率,得到所述两个虚拟晶体的晶体对接收效率,并利用所述晶体对接收效率对所述延迟随机符合数据进行去噪声处理,得到去噪声处理后的随机符合数据。
在一个例子中,所述处理模块1502,具体用于在根据所述各个实际晶体的晶体接收效率,得到所述两个虚拟晶体的晶体对接收效率的过程中,针对所述响应线上的虚拟晶体i和虚拟晶体j,根据所述虚拟晶体i在所述响应线上对应的m个实际晶体的晶体接收效率、所述虚拟晶体j在所述响应线上对应的m个实际晶体的晶体接收效率,得到所述虚拟晶体i和所述虚拟晶体j的晶体对接收效率。
在一个例子中,所述处理模块1502,具体用于在利用所述晶体对接收效率对所述延迟随机符合数据进行去噪声处理,得到去噪声处理后的随机符合数据的过程中,根据包含虚拟晶体j的一组虚拟晶体的集合A、包含虚拟晶体i的一组虚拟晶体的集合B、虚拟晶体i和属于集合B的虚拟晶体l的晶体对接收效率、虚拟晶体j和属于集合A的虚拟晶体k的晶体对接收效率、虚拟晶体l和虚拟晶体k的晶体对接收效率,对虚拟晶体i和虚拟晶体l对应的延迟随机符合数据、虚拟晶体j和虚拟晶体k对应的延迟随机符合数据、虚拟晶体l和虚拟晶体k对应的延迟随机符合数据,进行去噪声处理,得到去噪声处理后的虚拟晶体i和虚拟晶体j的随机符合数据。
本申请实施例的图像重建的功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台图像重建设备执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,RandomAccess Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述仅为本申请的较佳实施例而已,并不用以限制本申请,凡在本申请的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本申请保护的范围之内。
Claims (14)
1.一种图像重建的方法,应用于包括多个实际晶体的医疗设备,其特征在于,该方法包括:
获得连续进床扫描模式下的扫描数据;
将实际晶体虚拟对应一个以上虚拟晶体,根据所述扫描数据获得所述虚拟晶体所在响应线的延迟随机符合数据;
利用实际晶体的晶体接收效率对所述延迟随机符合数据进行去噪声处理,得到去噪声处理后的随机符合数据;
利用去噪声处理后的随机符合数据进行图像重建。
2.根据权利要求1所述的方法,其特征在于,该方法还包括:
获得实际晶体的单光子计数率,利用所述实际晶体的单光子计数率得到所述实际晶体的晶体接收效率。
3.根据权利要求2所述的方法,其特征在于,所述获得实际晶体的单光子计数率的过程,包括:
统计所述实际晶体接收到的单光子计数,并根据所述单光子计数以及所述单光子计数的接收时间,得到所述实际晶体的单光子计数率;或者,
获得对应一个扫描位置的多个单光子计数分量,根据所述多个单光子计数分量,得到与所述扫描位置对应的虚拟晶体的单光子计数率;利用所述虚拟晶体的单光子计数率,得到所述虚拟晶体对应的实际晶体的单光子计数率;
其中,一个单光子计数分量是所述扫描位置运动到对应的一个实际晶体时,所述实际晶体接收到的单光子计数,且所述实际晶体与所述虚拟晶体在对应所述扫描位置时与所述扫描位置的相对位置关系相同。
4.根据权利要求2所述的方法,其特征在于,利用所述实际晶体的单光子计数率得到所述实际晶体的晶体接收效率的过程,包括:
在不同剂量下,根据所述实际晶体的单光子计数率,以及单光子计数率与晶体接收效率之间的函数关系,得到所述实际晶体的晶体接收效率,所述函数关系是根据不同剂量下计算的晶体接收效率得到;或者,
在同一剂量下,获得所述实际晶体的单光子计数率与所有实际晶体的单光子计数率的均值之间的比值,并将所述比值作为所述晶体接收效率。
5.根据权利要求1所述的方法,其特征在于,所述利用实际晶体的晶体接收效率对所述延迟随机符合数据进行去噪声处理,得到去噪声处理后的随机符合数据的过程,包括:
确定响应线上的两个虚拟晶体、以及所述两个虚拟晶体在所述响应线上对应的各个实际晶体,并获取所述各个实际晶体的晶体接收效率;
根据所述各个实际晶体的晶体接收效率,得到所述两个虚拟晶体的晶体对接收效率,并利用所述晶体对接收效率对所述延迟随机符合数据进行去噪声处理,得到去噪声处理后的随机符合数据。
6.根据权利要求5所述的方法,其特征在于,根据所述各个实际晶体的晶体接收效率,得到所述两个虚拟晶体的晶体对接收效率的过程,包括:
针对所述响应线上的虚拟晶体i和虚拟晶体j,根据所述虚拟晶体i在所述响应线上对应的m个实际晶体的晶体接收效率、所述虚拟晶体j在所述响应线上对应的m个实际晶体的晶体接收效率,得到所述虚拟晶体i和所述虚拟晶体j的晶体对接收效率。
7.根据权利要求6所述的方法,其特征在于,利用所述晶体对接收效率对所述延迟随机符合数据进行去噪声处理,得到去噪声处理后的随机符合数据的过程,包括:
根据包含虚拟晶体j的一组虚拟晶体的集合A、包含虚拟晶体i的一组虚拟晶体的集合B、虚拟晶体i和属于集合B的虚拟晶体l的晶体对接收效率、虚拟晶体j和属于集合A的虚拟晶体k的晶体对接收效率、虚拟晶体l和虚拟晶体k的晶体对接收效率,对虚拟晶体i和虚拟晶体l对应的延迟随机符合数据、虚拟晶体j和虚拟晶体k对应的延迟随机符合数据、虚拟晶体l和虚拟晶体k对应的延迟随机符合数据,进行去噪声处理,得到去噪声处理后的虚拟晶体i和虚拟晶体j的随机符合数据。
8.一种图像重建的装置,应用于包括多个实际晶体的医疗设备,其特征在于,该装置包括:
获得模块,用于获得连续进床扫描模式下的扫描数据;将实际晶体虚拟对应一个以上虚拟晶体,根据所述扫描数据获得所述虚拟晶体所在响应线的延迟随机符合数据;
处理模块,用于利用实际晶体的晶体接收效率对所述延迟随机符合数据进行去噪声处理,得到去噪声处理后的随机符合数据;
重建模块,用于利用去噪声处理后的随机符合数据进行图像重建。
9.根据权利要求8所述的装置,其特征在于,
所述获得模块,还用于获得实际晶体的单光子计数率,利用所述实际晶体的单光子计数率得到所述实际晶体的晶体接收效率。
10.根据权利要求9所述的装置,其特征在于,
所述获得模块,具体用于在获得实际晶体的单光子计数率的过程中,统计所述实际晶体接收到的单光子计数,并根据所述单光子计数以及所述单光子计数的接收时间,得到所述实际晶体的单光子计数率;或者,
获得对应一个扫描位置的多个单光子计数分量,根据所述多个单光子计数分量,得到与所述扫描位置对应的虚拟晶体的单光子计数率;利用所述虚拟晶体的单光子计数率,得到所述虚拟晶体对应的实际晶体的单光子计数率;
其中,一个单光子计数分量是所述扫描位置运动到对应的一个实际晶体时,所述实际晶体接收到的单光子计数,且所述实际晶体与所述虚拟晶体在对应所述扫描位置时与所述扫描位置的相对位置关系相同。
11.根据权利要求9所述的装置,其特征在于,
所述获得模块,具体用于在利用所述实际晶体的单光子计数率得到所述实际晶体的晶体接收效率的过程中,在不同剂量下,根据所述实际晶体的单光子计数率,以及单光子计数率与晶体接收效率之间的函数关系,得到所述实际晶体的晶体接收效率,所述函数关系是根据不同剂量下计算的晶体接收效率得到;或者,在同一剂量下,获得所述实际晶体的单光子计数率与所有实际晶体的单光子计数率的均值之间的比值,并将所述比值作为所述晶体接收效率。
12.根据权利要求8所述的装置,其特征在于,
所述处理模块,具体用于在利用实际晶体的晶体接收效率对所述延迟随机符合数据进行去噪声处理,得到去噪声处理后的随机符合数据的过程中,确定响应线上的两个虚拟晶体、以及所述两个虚拟晶体在所述响应线上对应的各个实际晶体,并获取所述各个实际晶体的晶体接收效率;根据所述各个实际晶体的晶体接收效率,得到所述两个虚拟晶体的晶体对接收效率,并利用所述晶体对接收效率对所述延迟随机符合数据进行去噪声处理,得到去噪声处理后的随机符合数据。
13.根据权利要求12所述的装置,其特征在于,
所述处理模块,具体用于在根据所述各个实际晶体的晶体接收效率,得到所述两个虚拟晶体的晶体对接收效率的过程中,针对所述响应线上的虚拟晶体i和虚拟晶体j,根据所述虚拟晶体i在所述响应线上对应的m个实际晶体的晶体接收效率、所述虚拟晶体j在所述响应线上对应的m个实际晶体的晶体接收效率,得到所述虚拟晶体i和所述虚拟晶体j的晶体对接收效率。
14.根据权利要求13所述的装置,其特征在于,
所述处理模块,具体用于在利用所述晶体对接收效率对所述延迟随机符合数据进行去噪声处理,得到去噪声处理后的随机符合数据的过程中,根据包含虚拟晶体j的一组虚拟晶体的集合A、包含虚拟晶体i的一组虚拟晶体的集合B、虚拟晶体i和属于集合B的虚拟晶体l的晶体对接收效率、虚拟晶体j和属于集合A的虚拟晶体k的晶体对接收效率、虚拟晶体l和虚拟晶体k的晶体对接收效率,对虚拟晶体i和虚拟晶体l对应的延迟随机符合数据、虚拟晶体j和虚拟晶体k对应的延迟随机符合数据、虚拟晶体l和虚拟晶体k对应的延迟随机符合数据,进行去噪声处理,得到去噪声处理后的虚拟晶体i和虚拟晶体j的随机符合数据。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610355835.8A CN106023278A (zh) | 2016-05-25 | 2016-05-25 | 一种图像重建的方法和装置 |
CN201710249045.6A CN107133998A (zh) | 2016-05-25 | 2017-04-17 | 图像重建的方法和装置 |
US15/604,859 US10535165B2 (en) | 2016-05-25 | 2017-05-25 | Reconstructing image |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610355835.8A CN106023278A (zh) | 2016-05-25 | 2016-05-25 | 一种图像重建的方法和装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106023278A true CN106023278A (zh) | 2016-10-12 |
Family
ID=57094275
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610355835.8A Pending CN106023278A (zh) | 2016-05-25 | 2016-05-25 | 一种图像重建的方法和装置 |
CN201710249045.6A Pending CN107133998A (zh) | 2016-05-25 | 2017-04-17 | 图像重建的方法和装置 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710249045.6A Pending CN107133998A (zh) | 2016-05-25 | 2017-04-17 | 图像重建的方法和装置 |
Country Status (2)
Country | Link |
---|---|
US (1) | US10535165B2 (zh) |
CN (2) | CN106023278A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106725573A (zh) * | 2017-01-16 | 2017-05-31 | 沈阳东软医疗系统有限公司 | 一种连续进床pet系统的计数丢失校正方法和装置 |
CN107464270A (zh) * | 2017-07-17 | 2017-12-12 | 沈阳东软医疗系统有限公司 | 一种图像重建方法和装置 |
CN110097611A (zh) * | 2019-04-28 | 2019-08-06 | 上海联影智能医疗科技有限公司 | 图像重建方法、装置、设备及存储介质 |
CN110384480A (zh) * | 2018-04-18 | 2019-10-29 | 佳能株式会社 | 被检体信息取得装置、被检体信息处理方法和存储介质 |
CN112971824A (zh) * | 2021-02-08 | 2021-06-18 | 上海联影医疗科技股份有限公司 | Pet动态图像扫描方法、装置和计算机设备 |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109833055B (zh) * | 2019-01-07 | 2023-07-04 | 东软医疗系统股份有限公司 | 图像重建方法和装置 |
CN109965897B (zh) * | 2019-05-10 | 2022-03-22 | 上海联影医疗科技股份有限公司 | Pet扫描仪校正方法、装置、计算机设备和可读存储介质 |
CN111080734B (zh) * | 2019-11-25 | 2023-10-20 | 中国科学院深圳先进技术研究院 | 一种处理正电子发射断层扫描pet数据的方法及终端 |
US11249206B2 (en) * | 2020-01-06 | 2022-02-15 | Canon Medical Systems Corporation | Method and system for PET detector efficiency normalization |
US20220350037A1 (en) * | 2021-04-30 | 2022-11-03 | Siemens Medical Solutions Usa, Inc. | Mean randoms estimation from list mode data |
CN115100315B (zh) * | 2022-08-24 | 2022-11-15 | 威海市博华医疗设备有限公司 | 基于ct的高清成像方法及装置 |
Family Cites Families (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6373059B1 (en) * | 2000-10-31 | 2002-04-16 | Ge Medical Systems Global Technology Company, Llc | PET scanner septa |
US7414246B2 (en) * | 2006-01-03 | 2008-08-19 | Koninklijke Philips Electronics N.V. | Achieving accurate time-of-flight calibrations with a stationary coincidence point source |
US8487264B2 (en) * | 2008-07-31 | 2013-07-16 | Shimadzu Corporation | Radiation tomography apparatus |
EP2433261B1 (en) * | 2009-05-20 | 2014-03-05 | Koninklijke Philips N.V. | Continuous time-of-flight scatter simulation method |
CN102028485A (zh) * | 2009-10-01 | 2011-04-27 | 株式会社东芝 | Pet扫描仪系统以及其中的图像重建方法 |
CN202092697U (zh) * | 2010-12-29 | 2011-12-28 | 兰州海默科技股份有限公司 | 多相流的正电子断层成像装置 |
CN102789510B (zh) * | 2011-05-18 | 2015-02-18 | 上海生物医学工程研究中心 | 一种获取pet系统几何校正参数的方法 |
CN102274041B (zh) * | 2011-05-23 | 2013-09-04 | 江苏中惠医疗科技股份有限公司 | 正电子发射断层成像电子学系统的全数字化处理装置 |
CN103099637B (zh) * | 2013-02-21 | 2014-12-10 | 清华大学 | 一种用于双平板pet探测器的图像重建方法 |
CN103099639B (zh) * | 2013-03-01 | 2014-12-24 | 江苏中惠医疗科技股份有限公司 | Pet成像系统的环形拓扑结构及其实现方法 |
CN103417235B (zh) * | 2013-07-31 | 2015-09-30 | 沈阳东软医疗系统有限公司 | 一种随机噪声的校正方法和装置 |
CN103559391B (zh) * | 2013-10-23 | 2016-10-05 | 沈阳东软医疗系统有限公司 | 判定随机符合事件和判定真符合事件的方法及装置 |
CN103996213A (zh) * | 2014-04-01 | 2014-08-20 | 平生医疗科技(昆山)有限公司 | 一种pet图像重建方法及系统 |
CN104146726B (zh) * | 2014-08-19 | 2016-09-14 | 中国科学院高能物理研究所 | Pet系统符合探测响应的生成方法 |
CN104408763B (zh) * | 2014-10-29 | 2017-08-11 | 沈阳东软医疗系统有限公司 | 一种图像重建方法及装置 |
CN105125231B (zh) * | 2015-09-18 | 2018-02-16 | 沈阳东软医疗系统有限公司 | 一种pet图像环状伪影的去除方法和装置 |
-
2016
- 2016-05-25 CN CN201610355835.8A patent/CN106023278A/zh active Pending
-
2017
- 2017-04-17 CN CN201710249045.6A patent/CN107133998A/zh active Pending
- 2017-05-25 US US15/604,859 patent/US10535165B2/en active Active
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106725573A (zh) * | 2017-01-16 | 2017-05-31 | 沈阳东软医疗系统有限公司 | 一种连续进床pet系统的计数丢失校正方法和装置 |
CN107464270A (zh) * | 2017-07-17 | 2017-12-12 | 沈阳东软医疗系统有限公司 | 一种图像重建方法和装置 |
CN107464270B (zh) * | 2017-07-17 | 2020-08-11 | 东软医疗系统股份有限公司 | 一种图像重建方法和装置 |
CN110384480A (zh) * | 2018-04-18 | 2019-10-29 | 佳能株式会社 | 被检体信息取得装置、被检体信息处理方法和存储介质 |
CN110097611A (zh) * | 2019-04-28 | 2019-08-06 | 上海联影智能医疗科技有限公司 | 图像重建方法、装置、设备及存储介质 |
CN110097611B (zh) * | 2019-04-28 | 2023-09-22 | 上海联影智能医疗科技有限公司 | 图像重建方法、装置、设备及存储介质 |
CN112971824A (zh) * | 2021-02-08 | 2021-06-18 | 上海联影医疗科技股份有限公司 | Pet动态图像扫描方法、装置和计算机设备 |
WO2022166988A1 (zh) * | 2021-02-08 | 2022-08-11 | 上海联影医疗科技股份有限公司 | Pet动态扫描方法、装置和计算机设备 |
Also Published As
Publication number | Publication date |
---|---|
CN107133998A (zh) | 2017-09-05 |
US10535165B2 (en) | 2020-01-14 |
US20170345189A1 (en) | 2017-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106023278A (zh) | 一种图像重建的方法和装置 | |
CN106725573B (zh) | 一种连续进床pet系统的计数丢失校正方法和装置 | |
EP3224801B1 (en) | Multi-modality imaging system and method | |
CN106974671B (zh) | 正电子发射断层成像系统及其图像重建方法 | |
CN105844599B (zh) | 一种图像重建的去噪声方法和装置 | |
CN105125231B (zh) | 一种pet图像环状伪影的去除方法和装置 | |
JP2012118079A (ja) | Petデータ処理方法、petデータ処理装置、コンピュータ読み取り可能な記録媒体、及びデータ処理方法 | |
US10482634B2 (en) | Systems and methods for imaging with anisotropic voxels | |
CN104408756A (zh) | 一种pet图像重建方法及装置 | |
Krzemien et al. | Overview of the software architecture and data flow for the J-PET tomography device | |
US10222490B2 (en) | PET scanner with emission and transmission structures in a checkerboard configuration | |
JP6037048B2 (ja) | 輪郭画像生成装置および核医学診断装置 | |
Lee et al. | Impact of system design parameters on image figures of merit for a mouse PET scanner | |
CN106859686B (zh) | 成像方法和成像系统 | |
KR102283454B1 (ko) | 양전자 방출 단층 촬영기의 비행시간 정보를 이용하여 관측된 데이터에 포함된 반응깊이 정보를 보정하여 영상을 재구성하는 영상재구성방법 | |
CN101524278A (zh) | 一种图像深度效应校正方法 | |
El Bitar et al. | Acceleration of fully 3D Monte Carlo based system matrix computation for image reconstruction in small animal SPECT | |
Gundlich et al. | From 2D PET to 3D PET: issues of data representation and image reconstruction | |
JP5710166B2 (ja) | 画像解析装置及び方法 | |
JP2011002306A (ja) | Pet装置の逐次近似画像再構成法 | |
CN111080737A (zh) | 图像重建方法、装置及pet扫描系统 | |
JP4317124B2 (ja) | 画像処理システムおよび画像処理方法 | |
Leadbeater | The development of positron imaging systems for applications in industrial process tomography | |
WO2021153555A1 (ja) | 減弱分布画像作成装置、画像処理装置、放射線断層撮影システム、減弱分布画像作成方法、画像処理方法 | |
US20230206516A1 (en) | Scatter estimation for pet from image-based convolutional neural network |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20161012 |
|
WD01 | Invention patent application deemed withdrawn after publication |