CN117338315A - 一种光子计数ct的成像方法和系统 - Google Patents
一种光子计数ct的成像方法和系统 Download PDFInfo
- Publication number
- CN117338315A CN117338315A CN202210737557.8A CN202210737557A CN117338315A CN 117338315 A CN117338315 A CN 117338315A CN 202210737557 A CN202210737557 A CN 202210737557A CN 117338315 A CN117338315 A CN 117338315A
- Authority
- CN
- China
- Prior art keywords
- target object
- current
- parameters
- historical
- scanning
- 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 58
- 238000013170 computed tomography imaging Methods 0.000 title claims abstract description 9
- 238000012545 processing Methods 0.000 claims abstract description 150
- 238000003384 imaging method Methods 0.000 claims description 91
- 239000000463 material Substances 0.000 claims description 45
- 238000000354 decomposition reaction Methods 0.000 claims description 22
- 230000009467 reduction Effects 0.000 claims description 18
- 238000004422 calculation algorithm Methods 0.000 claims description 11
- 230000006870 function Effects 0.000 claims description 8
- 230000010354 integration Effects 0.000 claims description 8
- 230000008569 process Effects 0.000 description 21
- 238000002591 computed tomography Methods 0.000 description 18
- 230000008859 change Effects 0.000 description 15
- 238000012986 modification Methods 0.000 description 15
- 230000004048 modification Effects 0.000 description 15
- 239000000126 substance Substances 0.000 description 15
- 239000002872 contrast media Substances 0.000 description 11
- 230000005855 radiation Effects 0.000 description 8
- 210000000577 adipose tissue Anatomy 0.000 description 7
- 238000003745 diagnosis Methods 0.000 description 7
- 210000001519 tissue Anatomy 0.000 description 7
- 230000008901 benefit Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 5
- 210000000056 organ Anatomy 0.000 description 5
- 238000002059 diagnostic imaging Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000006872 improvement Effects 0.000 description 4
- 238000002595 magnetic resonance imaging Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 238000002604 ultrasonography Methods 0.000 description 4
- 230000005684 electric field Effects 0.000 description 3
- 238000004497 NIR spectroscopy Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 229910052791 calcium Inorganic materials 0.000 description 2
- 239000011575 calcium Substances 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 239000010931 gold Substances 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 229910052740 iodine Inorganic materials 0.000 description 2
- 239000011630 iodine Substances 0.000 description 2
- 238000012014 optical coherence tomography Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000002600 positron emission tomography Methods 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 238000004445 quantitative analysis Methods 0.000 description 2
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- LEHOTFFKMJEONL-UHFFFAOYSA-N Uric Acid Chemical compound N1C(=O)NC(=O)C2=C1NC(=O)N2 LEHOTFFKMJEONL-UHFFFAOYSA-N 0.000 description 1
- TVWHNULVHGKJHS-UHFFFAOYSA-N Uric acid Natural products N1C(=O)NC(=O)C2NC(=O)NC21 TVWHNULVHGKJHS-UHFFFAOYSA-N 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000000712 assembly Effects 0.000 description 1
- 238000000429 assembly Methods 0.000 description 1
- 230000003190 augmentative effect Effects 0.000 description 1
- 229910052788 barium Inorganic materials 0.000 description 1
- DSAJWYNOEDNPEQ-UHFFFAOYSA-N barium atom Chemical compound [Ba] DSAJWYNOEDNPEQ-UHFFFAOYSA-N 0.000 description 1
- 238000012984 biological imaging Methods 0.000 description 1
- 230000037396 body weight Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 229940039231 contrast media Drugs 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 239000002784 hot electron Substances 0.000 description 1
- 238000002608 intravascular ultrasound Methods 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000003607 modifier Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- -1 rolling (Gd) Inorganic materials 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 239000000725 suspension Substances 0.000 description 1
- 239000013077 target material Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000005979 thermal decomposition reaction Methods 0.000 description 1
- 229940116269 uric acid Drugs 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/42—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4208—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
- A61B6/4241—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with 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 for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computerised tomographs
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/482—Diagnostic techniques involving multiple energy imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/54—Control of apparatus or devices for radiation diagnosis
Abstract
本说明书实施例提供一种光子计数CT的成像方法和系统。所述方法包括:获取目标对象的当前基本信息;获取所述目标对象的参考基本信息、历史扫描参数和历史处理参数;基于所述当前基本信息、所述参考基本信息以及所述历史扫描参数,确定所述目标对象的当前扫描参数;和/或,基于所述当前基本信息、所述参考基本信息以及所述历史处理参数,确定所述目标对象的当前处理参数。所述方法还包括:基于所述当前扫描参数或预设扫描参数对所述目标对象进行扫描,获取探测器输出结果;以及基于所述当前处理参数,对所述探测器输出结果进行处理,获取所述目标对象的目标图像。
Description
技术领域
本说明书涉及医学成像领域,特别涉及一种成像方法和系统。
背景技术
光子计数CT(Photon-Counting Computed Tomography,PCCT)因具有可实现材料成分分析、降低患者辐射剂量、提高CT定量分析准确性和实现超高空间分辨率等优点,已被广泛应用于医学成像。在PCCT成像中,操作者可通过断层扫描图像获得目标对象的物质定量结果,基于定量结果进行有理可依的定性判断。但是,物质分解定量结果的精确度同时取决于扫描系统、扫描参数、对数据的后续处理过程。
因此,希望提出一种光子计数CT的成像方法和系统,可以确定目标对象的当次的扫描参数和/或后处理参数。
发明内容
本说明书一个方面提供一种光子计数CT的成像方法。所述方法包括:获取目标对象的当前基本信息;获取所述目标对象的参考基本信息和历史参数,所述历史参数包括历史扫描参数和历史处理参数中的至少一个;基于所述当前基本信息、所述参考基本信息以及所述历史参数,确定所述目标对象的当前扫描参数和/或当前处理参数。
在一些实施例中,所述方法进一步包括:基于所述当前扫描参数或预设扫描参数对所述目标对象进行扫描,获取探测器输出结果;以及,基于所述当前处理参数,对所述探测器输出结果进行处理,获取所述目标对象的所述目标图像。
在一些实施例中,所述参考基本信息和所述历史参数来自所述目标对象的个人数据库,所述方法进一步包括:从所述数据库中获取所述目标对象的历史图像;将所述目标对象的所述目标图像与所述历史图像进行对比,并输出对比结果。
在一些实施例中,所述目标图像包括所述目标对象的物质分解图像,所述目标图像与所述历史图像具有相同或相近的置信度,所述置信度反映物质分解图像中所示定量结果的准确性。
在一些实施例中,所述对比结果包括以下中的至少一种:所述目标对象的当前基物质图像与历史基物质图像的对比结果、当前各基物质按照第一权重组合获得的当前组合图像与历史各基物质按照第二权重组合获得的历史组合图像的对比结果,其中,所述第一权重组合使得所述当前组合图像与所述历史组合图像中的非目标区域的数值结果保持一致。
在一些实施例中,基于所述当前基本信息、所述参考基本信息以及所述历史参数,确定所述目标对象的当前扫描参数和/或当前处理参数,包括:将所述目标对象的所述当前基本信息与所述参考基本信息进行比较,获取信息变化量,所述信息变化量与所述目标对象的衰减信息相关;基于所述信息变化量和所述历史参数,确定所述目标对象的当前扫描参数和/或当前处理参数。
在一些实施例中,所述基于所述信息变化量和所述历史参数,确定所述目标对象的当前扫描参数和/或当前处理参数,包括:将所述信息变化量转化为扫描参数变化量,并基于所述扫描参数变化量和所述历史扫描参数,确定所述目标对象的当前扫描参数;和/或,将所述信息变化量转化为处理参数变化量,并基于所述处理参数变化量和所述历史处理参数,确定所述目标对象的当前处理参数。
在一些实施例中,所述扫描参数包括以下中的至少一种:管电压、管电流、能量仓的数量、各能量仓的阈值、单视野下的数据对应的积分时间、机架转速和扫描层厚;所述处理参数包括以下中的至少一种:各能量仓的加权因子、降噪算法、降噪等级、滤波函数和重建矩阵。
本说明书另一个方面提供一种光子计数CT的成像系统。所述系统包括:获取模块,用于获取目标对象的当前基本信息,以及所述目标对象的参考基本信息和历史参数,所述历史参数包括历史扫描参数和历史处理参数中的至少一个;确定模块,用于基于所述当前基本信息、所述参考基本信息以及所述历史参数,确定所述目标对象的当前扫描参数和/或当前处理参数;扫描模块,用于基于所述当前扫描参数或预设扫描参数,对所述目标对象进行扫描,以获取探测器输出结果;数据处理模块,用于基于所述当前处理参数,对所述探测器输出结果进行处理,以获取所述目标对象的目标图像。
本说明书另一个方面提供一种计算机可读存储介质,所述存储介质存储计算机指令,当计算机读取所述计算机指令时,所述计算机执行如前所述的成像方法。
本说明书另一个方面提供一种光子计数CT的成像装置,该装置包括处理器以及存储器,所述存储器用于存储指令,当所述处理器执行所述指令时实现如前所述的成像方法。
附图说明
本说明书将以示例性实施例的方式进一步说明,这些示例性实施例将通过附图进行详细描述。这些实施例并非限制性的,在这些实施例中,相同的编号表示相同的结构,其中:
图1是根据本说明书一些实施例所示的示例性成像系统的应用场景示意图;
图2是根据本说明书一些实施例所示的示例性成像系统的模块示意图;
图3是根据本说明书一些实施例所示的示例性成像方法的流程示意图;
图4是根据本说明书另一些实施例所示的示例性成像方法的流程示意图;
图5是根据本说明书另一些实施例所示的示例性成像方法的流程示意图。
具体实施方式
为了更清楚地说明本说明书实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单的介绍。显而易见地,下面描述中的附图仅仅是本说明书的一些示例或实施例,对于本领域的普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图将本说明书应用于其它类似情景。除非从语言环境中显而易见或另做说明,图中相同标号代表相同结构或操作。
应当理解,本文使用的“系统”、“装置”、“单元”和/或“模块”是用于区分不同级别的不同组件、元件、部件、部分或装配的一种方法。然而,如果其他词语可实现相同的目的,则可通过其他表达来替换所述词语。
如本说明书和权利要求书中所示,除非上下文明确提示例外情形,“一”、“一个”、“一种”和/或“该”等词并非特指单数,也可包括复数。一般说来,术语“包括”与“包含”仅提示包括已明确标识的步骤和元素,而这些步骤和元素不构成一个排它性的罗列,方法或者设备也可能包含其它的步骤或元素。
本说明书中使用了流程图用来说明根据本说明书的实施例的系统所执行的操作,相关描述是为帮助更好地理解医学成像方法和/或系统。应当理解的是,前面或后面操作不一定按照顺序来精确地执行。相反,可以按照倒序或同时处理各个步骤。同时,也可以将其他操作添加到这些过程中,或从这些过程移除某一步或数步操作。
图1是根据本说明书一些实施例所示的示例性成像系统的应用场景示意图。
如图1所示,成像系统100可以包括成像设备110、处理设备120、一个或以上终端130、存储设备140以及网络150。成像系统100中的组件可以各种方式中的一种或以上方式连接。仅作为示例,如图1中所示,成像设备110可以通过网络150连接到处理设备120。作为另一示例,成像设备110可以直接连接到处理设备120,如成像设备110和处理设备120可以如图中虚线双向箭头所指示进行连接。作为又一示例,存储设备140可以直接与处理设备120连接(图1中未示出)或通过网络150连接。作为再一示例,一个或以上终端130可以直接与处理设备120连接(如由连接终端130和处理设备120的虚线双向箭头所示)或通过网络150连接。
成像设备110可以用于对检测区域内的目标对象进行扫描,得到该目标对象的扫描数据(例如,历史图像、目标图像等)。在一些实施例中,目标对象可以包括生物对象和/或非生物对象。例如,目标对象可以包括身体的特定部分,例如头部、胸部、腹部等,或其组合。又例如,目标对象可以是有生命或无生命的有机和/或无机物质的人造成分。在一些实施例中,与目标对象有关的扫描数据可以包括目标对象的投影数据、一个或以上扫描图像等。
在一些实施例中,成像设备110可以包括用于疾病诊断或研究目的的非侵入性生物成像装置。例如,成像设备110可以包括单模态扫描仪和/或多模态扫描仪。单模态扫描仪可以包括例如超声波扫描仪、X射线扫描仪、计算机断层扫描(CT)扫描仪、磁共振成像(MRI)扫描仪、超声检查仪、正电子发射断层扫描(PET)扫描仪、光学相干断层扫描(OCT)扫描仪、超声(US)扫描仪、血管内超声(IVUS)扫描仪、近红外光谱(NIRS)扫描仪、远红外(FIR)扫描仪等。多模态扫描仪可以包括例如X射线成像-磁共振成像(X射线-MRI)扫描仪、正电子发射断层扫描-X射线成像(PET-X射线)扫描仪、单光子发射计算机断层扫描-磁共振成像(SPECT-MRI)扫描仪、正电子发射断层扫描-计算机断层摄影(PET-CT)扫描仪、数字减影血管造影-磁共振成像(DSA-MRI)扫描仪等。上面提供的扫描仪仅用于说明目的,而无意限制本说明书的范围。如本文所用,术语“成像模态”或“模态”广泛地是指收集、生成、处理和/或分析目标对象的成像信息的成像方法或技术。
在一些实施例中,成像设备110可以包括用于执行成像和/或相关分析的模块和/或组件。例如,成像设备110可以包括射线发生装置、附属装置和成像装置。其中,射线发生装置是指完成射线(例如,X射线)产生并对其进行控制的装置。射线附属装置是指为满足临床诊疗的需要而设计的各种与射线发生装置配套的设施,例如,可以包括机械设备如检查床、诊断床、导管床、摄影床等、各种支撑、悬吊装置、制动装置、滤线栅、保持装置、遮线器等。在一些实施例中,射线成像装置可以有多种形式,例如,数字成像装置可以包括探测器、计算机系统及图像处理软件等;其他成像装置可以包括荧光屏、胶片暗匣、影像增强器、影像电视等。
在一些实施例中,由成像设备110获取的数据(例如,目标对象的目标图像等)可以被传送到处理设备120以供进一步分析。附加地或替代地,由成像设备110获取的数据可以被发送到终端设备(例如,终端130)用于显示和/或存储设备(例如,存储设备140)用于存储。
处理设备120可以处理从成像设备110、终端130、存储设备140和/或其他存储设备获取和/或提取的数据和/或信息。例如,处理设备120可以从终端130或存储设备140获取目标对象的当前基本信息、参考基本信息、历史扫描参数和历史处理参数,并基于这些信息确定目标对象的当前扫描参数和/或当前处理参数。又如,处理设备120可以从成像设备110获取目标对象的目标图像,将目标图像与该目标对象的历史图像进行对比,并输出对比结果。
在一些实施例中,处理设备120可以是单个服务器或服务器组。服务器组可以是集中式或分布式的。在一些实施例中,处理设备120可以是本地的或远程的。在一些实施例中,处理设备120可以被实现在云平台上。仅作为示例,云平台可以包括私有云、公共云、混合云、社区云、分布云、内部云、多层云等或其任意组合。
在一些实施例中,处理设备120可以在计算设备上实现。在一些实施例中,处理设备120可以在终端(例如,终端130)上实现。在一些实施例中,处理设备120可以在成像设备(例如,成像设备110)上实现。例如,处理设备120可以集成在终端130和/或成像设备110中。
终端130可以与成像设备110和/或处理设备120连接,用于输入/输出信息和/或数据。例如,用户可以通过终端130与成像设备110进行交互,以控制成像设备110的一个或多个部件(例如,输入患者信息等)。又例如,成像设备110可以将生成的医学图像(例如,目标图像)输出至终端130,以展示给用户。
在一些实施例中,终端130可以包括移动设备131、平板计算机132、膝上型计算机133等,或其任意组合。在一些实施例中,移动设备131可以包括智能家居设备、可穿戴设备、智能移动设备、虚拟现实设备、增强现实设备等,或其任意组合。
在一些实施例中,一个或以上终端130可以远程操作成像设备110。在一些实施例中,终端130可以经由无线连接来操作成像设备110。在一些实施例中,一个或以上终端130可以是处理设备120的一部分。在一些实施例中,可以省略终端130。
存储设备140可以存储数据和/或指令。在一些实施例中,存储设备140可以存储从终端130和/或处理设备120获取的数据。例如,存储设备140可以存储能量仓的目标数量、扫描协议等。在一些实施例中,存储设备140可以存储数据和/或指令,处理设备120可以执行或使用该数据和/或指令来执行本申请中描述的示例性方法。
在一些实施例中,存储设备140可以包括大容量存储设备、可移动存储设备、易失性读写内存、只读内存(ROM)等或其任意组合。示例性大容量存储器可以包括磁盘、光盘、固态驱动器等。示例性可移动存储器可以包括闪存驱动器、软盘、光盘、内存卡、压缩盘、磁带等。示例性易失性读写内存可以包括随机存取内存(RAM)。在一些实施例中,存储设备140可以被实现在云平台上。在一些实施例中,存储设备140可以是处理设备120的一部分。
网络150可以包括可以促进成像系统100的信息和/或数据的交换的任何合适的网络。在一些实施例中,成像系统100的一个或以上组件(例如,成像设备110、一个或以上终端130、处理设备120或存储设备140)可以与成像系统100的一个或以上其他组件通信以传输信息和/或数据。在一些实施例中,网络150可以是任何类型的有线或无线网络或其组合。例如,网络150可以是和/或包括公共网络(例如,互联网)、私有网络(例如,局域网(LAN)、广域网(WAN)等)、有线网络(例如,以太网)、无线网络(例如,802.11网络、Wi-Fi网络等)、蜂窝网络(例如,长期演进(LTE)网络)、帧中继网络、虚拟专用网络(“VPN”)、卫星网络、电话网络、路由器、集线器、交换机、服务器计算机和/或其任何组合。在一些实施例中,网络150可以包括一个或以上网络接入点。
应当注意,对成像系统100的以上描述仅出于说明的目的,而无意于限制本说明书的范围。对于本领域普通技术人员而言,可以根据本说明书进行各种变型和修改。然而,这些变化和修改不脱离本说明书的范围。例如,成像设备110、处理设备120与终端130可以共用一个存储设备140,也可以有各自的存储设备。
光子计数探测器为CT扫描带来了真正的物质分辨能力。操作者可以通过断层扫描图像获得具有一定精确度的目标对象的物质分解定量结果。基于该定量结果,操作者可以在特定的临床场景(例如,肿瘤的性质、患者康复进展等)中进行有理可依的定性判断,从而可以减轻对医师临床经验的需求。然而,物质分解定量结果的精确度同时取决于扫描系统、扫描参数、对数据的后续处理过程。
本说明书实施例中提供一种光子计数CT的成像方法和系统,从每个目标对象的个人数据库中获取参考基本信息、历史扫描参数和历史处理参数,结合目标对象的当前基本信息,确定目标对象的当前扫描参数和/或当前处理参数,并基于当前扫描参数和当前处理参数获得目标对象的目标图像,以期使得目标对象的一系列扫描中每次扫描的扫描结果具有定量可对比性,便于追踪目标对象病情的发展,帮助提高诊断效率和诊断结果准确性。
图2是根据本说明书一些实施例所示的示例性成像系统的模块示意图。
如图2所示,在一些实施例中,成像系统200可以包括获取模块210、确定模块220、扫描模块230以及数据处理模块240。在一些实施例中,成像系统200可以集成在成像设备110或处理设备120中。
获取模块210可以用于获取目标对象的相关信息。在一些实施例中,获取模块210可以用于获取目标对象的当前基本信息。在一些实施例中,获取模块210可以用于获取目标对象的参考基本信息、历史扫描参数和历史处理参数。在一些实施例中,基本信息可以包括目标对象的3D模型、体脂率。在一些实施例中,参考基本信息、历史扫描参数和历史处理参数来自目标对象的个人数据库。
在一些实施例中,获取模块210可以用于从数据库中获取目标对象的历史图像。在一些实施例中,目标图像包括目标对象的物质分解图像。在一些实施例中,目标对象的目标图像与历史图像具有相同或相近的置信度,置信度反映物质分解图像中所示定量结果的准确性。确定模块220可以用于确定目标对象的当前扫描参数和/或当前处理参数。
在一些实施例中,确定模块220可以用于基于当前基本信息、参考基本信息以及历史扫描参数,确定目标对象的当前扫描参数。在一些实施例中,确定模块220可以用于将目标对象的当前基本信息与参考基本信息进行比较,获取信息变化量,并基于信息变化量和历史扫描参数,确定目标对象的当前扫描参数。其中,信息变化量与目标对象的衰减信息相关。在一些实施例中,确定模块220可以将信息变化量转化为扫描参数变化量,并基于扫描参数变化量和历史扫描参数,确定目标对象的当前扫描参数。
在一些实施例中,确定模块220可以用于基于当前基本信息、参考基本信息以及历史处理参数,确定目标对象的当前处理参数。在一些实施例中,确定模块220可以用于基于信息变化量和历史处理参数,确定目标对象的当前处理参数。在一些实施例中,确定模块220可以将信息变化量转化为处理参数变化量,并基于处理参数变化量和历史处理参数,确定目标对象的当前扫描参数。
在一些实施例中,扫描参数可以包括以下中的至少一种:造影剂类别、管电压、管电流、能量仓的数量、各能量仓的阈值、单视野下的数据对应的积分时间、机架转速和扫描层厚等。在一些实施例中,处理参数可以包括以下中的至少一种:各能量仓的加权因子、降噪算法、降噪等级、滤波函数和重建矩阵等。
扫描模块230可以用于对目标对象进行扫描。在一些实施例中,扫描模块230可以用于基于目标对象的当前扫描参数或预设扫描参数,对目标对象进行扫描,以获取探测器输出结果。
数据处理模块240可以用于对探测器输出结果进行处理。在一些实施例中,数据处理模块240可以用于基于目标对象的当前处理参数,对探测器输出结果进行处理,以获取目标对象的目标图像。
在一些实施例中,数据处理模块240可以用于将目标对象的目标图像与历史图像进行对比,并输出对比结果。在一些实施例中,对比结果可以包括以下中的至少一种:目标对象的当前基物质图像与历史基物质图像的对比结果、当前各基物质按照第一权重组合获得的当前组合图像与历史各基物质按照第二权重组合获得的历史组合图像的对比结果。其中,第一权重组合使得当前组合图像与历史组合图像中的非目标区域的数值结果保持一致。
需要注意的是,以上对于成像系统200及其模块的描述,仅为描述方便,并不能把本说明书限制在所举实施例范围之内。可以理解,对于本领域的技术人员来说,在了解该系统的原理后,可能在不背离这一原理的情况下,对各个模块进行任意组合,或者构成子系统与其他模块连接。然而,这些修正和改变仍在本说明书的范围之内。
图3是根据本说明书一些实施例所示的示例性成像方法的流程示意图。
在一些实施例中,成像方法300可以由成像设备110或处理设备120执行。例如,成像方法300可以以程序或指令的形式存储在存储装置(如存储设备140)中,当成像设备110或处理设备120执行该程序或指令时,可以实现成像方法300。在一些实施例中,成像方法300可以由成像系统200执行。
步骤310,获取目标对象的当前基本信息。在一些实施例中,步骤310可以由处理设备120或获取模块210执行。
目标对象可以指需要进行扫描的对象。在一些实施例中,目标对象可以包括需要进行扫描的患者,或患者的扫描部位。
当前信息可以指目标对象当次扫描的相关信息。相应地,当前基本信息可以指目标对象当次扫描对应的基本信息。
在一些实施例中,基本信息可以包括目标对象的3D模型。3D模型可以反映人体或扫描部位包含的器官/组织,器官/组织的结构、密度、体积、位置、厚度等信息。在一些实施例中,可以通过摄像装置(例如,深度摄像头)、红外成像仪、多角度的3D扫描(例如,CT扫描、PET扫描、MR扫描、DR扫描等)等方式,获取目标对象的当前3D模型。例如,可以基于红外成像仪获得的目标对象的3D图像,其能够通过热分解显示体内得器官分布、位置和大小等信息。
在一些实施例中,基本信息还可以包括目标对象的年龄、性别、身高、体重、体脂率等或其任意组合。在一些实施例中,可以从存储设备(例如,存储设备140)中获取目标对象的基本信息。在一些实施例中,可以从终端设备(例如,终端130)或成像设备(例如,成像设备110)获取目标对象的基本信息。在一些实施例中,可以从医疗系统获取目标对象的基本信息。
步骤320,获取目标对象的参考基本信息、历史扫描参数和历史处理参数。在一些实施例中,步骤320可以由处理设备120或获取模块210执行。
扫描参数可以反映对患者进行扫描时使用的参数的类型和/或相应的数值。在一些实施例中,扫描参数可以包括造影剂类别、管电压、管电流、能量仓的数量、各能量仓的阈值、单视野下的数据对应的积分时间、机架转速和扫描层厚等,或其任意组合。
造影剂(又称对比剂,contrast media)是为增强影像观察效果而注入(或服用)到人体组织或器官的化学制品。例如,造影剂类别可以包括碘(I)、金(Au)、轧(Gd)、钡(Ba)等。
医学成像时,射线(例如,X射线)是从射线管球中发出的,射线管球有阴极侧(灯丝)和阳极侧(靶面)。
管电压可以指加载到射线管的阴极侧和阳极侧之间的电压,以用来形成高压电场,使灯丝发出的热电子能在高压电场的加速下高速轰击靶面,激发射线。射线管球产生的射线束的最大光子能量等于高速电子流的最大能量,而高速电子流的最大能量取决于管电压的峰值,通过改变管电压可改变光子最大能量及射线谱。
管电流是指灯丝加热产生的电子,在阴阳两级高压电场作用下,向阳极高速运动,形成的电流。管电流和时间的乘积决定射线的量(发射光子数量)。
光子计数型探测器(Photon Counting Detector,PCD)可以通过设置不同的阈值,对比该设定阈值能量高的光子进行计数,从而将每个入射的光子进行分区间计数,让原始的CT数据带有能量维度的信息。其中,划分的能量区间称为能量仓,设定的阈值称为能量仓(Bin)阈值,可以反映每个能量仓的区间大小。能量仓的数量可以指需要划分的能量区间的个数。
单视野(view)下的数据对应的积分时间可以指CT成像时机架旋转一圈,采集单组数据(即对应单个视野的数据)所用的时间。例如,机架旋转一圈用时1秒钟,采集9800组数据,即9800个视野,则单视野下的数据对应的积分时间为1/9800。
机架转速即成像设备(例如,CT成像设备)的机架的旋转速度。
扫描层厚可以指成像设备(例如,CT成像设备)扫描一个周期经过的层厚。层厚可以指计算机体层成像断面图像在垂直于被检体(例如,目标对象)或检查床方向上覆盖的长度。扫描层厚与机器排数、探测器数量有关。机器越好,其排数越高,扫描层厚相对越薄,获得的重建图像越清晰。
处理参数可以反映扫描完成后,对扫描数据进行处理和/或基于扫描数据进行图像重建时使用的处理参数。在一些实施例中,处理参数可以包括各能量仓的加权因子、降噪算法、降噪等级、滤波函数和重建矩阵等或其任意组合。其中,各能量仓的加权因子可以反映对探测器的输出结果进行处理时,探测器的每个能量仓的输出数据的权重占比。各能量仓的加权因子、降噪算法、降噪等级、滤波函数和重建矩阵等参数可以影响物质分解定量结果。
在一些实施例中,处理参数还可以包括图像CT值、窗宽、窗位、伪影校正算法等或其任意组合。图像CT值、窗宽、窗位、伪影校正算法等可以影响物质分解图像的图像质量。
参考基本信息、历史扫描参数和历史处理参数等可以反映目标对象的历史扫描数据。例如,参考基本信息反映患者先前扫描对应的3D模型和/或体脂率等身体状态信息,历史扫描参数反映患者先前扫描使用的扫描参数,历史处理参数反映患者先前进行后处理时使用的处理参数。
在一些实施例中,参考基本信息、历史扫描参数和历史处理参数可以来自目标对象的个人数据库。例如,可以为每个患者建立个人数据库,当每次完成对患者的扫描和/或诊断后,将该次扫描对应的患者的体脂率、体重、3D模型等基本信息,患者扫描时使用的扫描参数、扫描完成后使用的后处理参数、扫描图像、影像检查结果等数据归档,存入与该患者对应的个人数据库中。
在一些实施例中,可以基于身份标识从数据库中获取目标对象的参考基本信息、历史扫描参数和历史处理参数。例如,可以根据患者的身份证号码、电话号码或就诊编号,匹配该患者的数据库,从该数据库中获取其参考基本信息、历史扫描参数和历史处理参数。
步骤330,基于当前基本信息、参考基本信息以及历史扫描参数,确定目标对象的当前扫描参数。在一些实施例中,步骤330可以由处理设备120或确定模块220执行。
在一些实施例中,可以将目标对象的当前基本信息与参考基本信息进行比较,获取信息变化量。其中,信息变化量与目标对象的衰减信息相关。衰减信息可以反映目标对象对射线的衰减程度。例如,衰减程度=衰减系数*衰减长度。
在一些实施例中,可以基于信息变化量和历史扫描参数,确定目标对象的当前扫描参数。在一些实施例中,可以基于信息变化量,确定是否调整历史扫描参数,以确定当前扫描参数。例如,可以当信息变化量小于或等于预设阈值时,直接将历史扫描参数确定为目标对象的当前扫描参数;当信息变化量大于预设阈值时,对历史扫描参数进行调整,以确定目标对象的当前扫描参数。
在一些实施例中,可以将信息变化量转化为扫描参数变化量,并基于扫描参数变化量和历史扫描参数,确定目标对象的当前扫描参数。更多相关内容可以参见图4(例如,步骤440)及其相关描述,此处不再赘述。
在一些实施例中,无法确定目标对象的信息变化量和/或扫描参数变化量时,可以基于目标对象的当前基本信息确定预设扫描参数。例如,患者第一次进行CT扫描时,由于尚不存在历史诊断数据,可以基于患者的扫描部位、3D模型、身高、年龄、性别、体重等中的一种或多种,确定该患者的预设扫描参数。又如,目标对象当前扫描的造影剂条件发生变化,导致扫描的Bin数、Bin阈值、管电流等参数的变化,从而无法确定目标对象的扫描参数变化量,此时,可以基于目标对象当前的造影剂类别、3D模型、体脂率等,确定该目标对象的预设扫描参数。
步骤340,基于当前基本信息、参考基本信息以及历史处理参数,确定目标对象的当前处理参数。在一些实施例中,步骤340可以由处理设备120或确定模块220执行。
在一些实施例中,可以基于信息变化量和历史处理参数,确定目标对象的当前处理参数。在一些实施例中,可以基于信息变化量,判断是否调整历史处理参数,以确定目标对象的当前处理参数。例如,可以当信息变化量小于或等于预设阈值时,直接将历史处理参数确定为目标对象的当前扫描参数;当信息变化量大于预设阈值时,对历史处理参数进行调整,以确定目标对象的当前处理参数。
在一些实施例中,当信息变化量超过预设阈值时,可以对历史扫描参数和历史处理参数同时进行调整,或仅调整其中的一项,以确定目标对象的当前扫描参数和当前处理参数。例如,可以将目标对象的历史扫描参数确定为当前扫描参数,基于信息变化量对目标对象的历史处理参数进行调整,确定当前处理参数。又如,可以基于信息变化量对目标对象的历史扫描参数进行调整,确定当前扫描参数,将目标对象的历史处理参数确定为当前处理参数。再如,可以基于目标对象的信息变化量,同时对历史扫描参数和历史处理参数进行调整,确定目标对象的当前扫描参数和当前处理参数。
当信息变化量大于预设阈值时,表示患者当前的衰减信息相较于上次扫描时的衰减信息变化较大,通过调整历史扫描参数和/或历史处理参数,可以保证患者前后扫描结果具有定量可比性的同时,使得当前扫描获取的图像满足诊断需求。
在一些实施例中,可以基于当前扫描参数对目标对象进行扫描,获取探测器输出结果;以及,基于当前处理参数,对探测器输出结果进行处理,获取目标对象的目标图像。在一些实施例中,可以将目标对象的目标图像与历史图像进行对比,并输出对比结果。更多相关内容可以参见图5及其相关描述,此处不再赘述。
应当注意的是,上述有关方法300的描述仅仅是为了示例和说明,而不限定本说明书的适用范围。对于本领域技术人员来说,在本说明书的指导下可以对方法300进行各种修正和改变。例如,可以先执行步骤340后执行步骤330,或同时执行步骤330和步骤340,或仅执行步骤330和步骤340中的一个。然而,这些修正和改变仍在本说明书的范围之内。
图4是根据本说明书另一些实施例所示的示例性成像方法的流程示意图。
在一些实施例中,成像方法400可以由成像设备110或处理设备120执行。例如,成像方法400可以以程序或指令的形式存储在存储装置(如存储设备140)中,当成像设备110或处理设备120执行该程序或指令时,可以实现成像方法400。在一些实施例中,成像方法400可以由成像系统200执行。
其中,步骤410:获取目标对象的当前基本信息,步骤420:获取目标对象的参考基本信息、历史扫描参数和历史处理参数,分别与方法300中步骤310和步骤320类似,更多内容可以参见图3及其相关描述,此处不再赘述。
步骤430,将目标对象的当前基本信息与参考基本信息进行比较,获取信息变化量。在一些实施例中,步骤430可以由处理设备120或确定模块220执行。
信息变化量可以反映目标对象的基本信息的变化带来的衰减程度的改变。
在一些实施例中,信息变化量可以基于衰减系数和衰减长度确定。不同材料(例如,人体的不同组织、器官、部位,不同类别的造影剂等)对不同能量的X射线的衰减系数不同。衰减长度取决于衰减模型(例如,目标对象的3D模型)在X射线光路中的长度,也可以理解为患者或扫描部位的厚度。
在一些实施例中,可以基于目标对象的当前基本信息的衰减系数和衰减长度,以及参考基本信息的衰减系数和衰减长度,确定信息变化量。在一些实施例中,可以基于目标对象的当前基本信息和参考基本信息之间的衰减长度差值以及平均衰减系数差值,确定信息变化量。
仅作为示例,用ΔB表示信息变化量,ΔB=(L+ΔL)×(α+Δα)-L×α。其中,L表示目标对象的参考基本信息B1对应的衰减长度;ΔL表示目标对象的当前基本信息B2和参考基本信息B1之间的衰减长度差值;α表示B1对应的平均衰减系数;Δα表示B1和B2之间的平均衰减系数差值。
在一些实施例中,可以基于目标对象的年龄、性别、身高、体重、体脂率以及3D模型等中的至少一种,确定衰减系数和/或衰减长度。例如,可以基于目标对象的参考基本信息B1对应的年龄、性别、身高、体重、体脂率以及3D模型,提取出患者在相应扫描体位和扫描部位下暴露在X射线光路中的衰减长度L,和/或,提取出患者在相应扫描体位和扫描部位下暴露在X射线光路中人体组织的平均衰减系数α。又如,可以基于目标对象的当前基本信息B2对应的年龄、性别、身高、体重,将目标对象简化为3D桶状模型或带入标准3D模型,以确定相应的衰减长度和衰减系数。再如,可以基于目标对象的当前基本信息B2对应的3D模型,将目标对象简化为3D桶状模型或带入标准3D模型,以确定相应的衰减长度和衰减系数。
可以理解,上述关于信息变化量的确定方式仅作为示例,在一些实施例中,可以通过其他可行的方式确定信息变化量,本说明书对此不作限制。
步骤440,基于信息变化量和历史扫描参数,确定当前扫描参数。在一些实施例中,步骤440可以由处理设备120或确定模块220执行。
在一些实施例中,当信息变化量小于或等于预设阈值时,可以将目标对象的历史扫描参数确定为当前扫描参数。预设阈值可以为任意合理的数值,本说明书对此不作限制。
在一些实施例中,当信息变化量大于预设阈值时,可以将信息变化量转化为扫描参数变化量,并基于扫描参数变化量和历史扫描参数,确定目标对象的当前扫描参数。扫描参数变化量可以反映各扫描参数需要调整的数值,或要调整为的目标数值。
在一些实施例中,扫描参数变化量可以包括造影剂类别、管电压、管电流、能量仓的数量、各能量仓的阈值、单视野下的数据对应的积分时间、机架转速和扫描层厚等参数中,一个参数或多个参数的变化量的组合。
在一些实施例中,可以基于造影剂的使用条件(例如,是否使用造影剂、使用的造影剂类别等),确定扫描参数变化量。例如,在不改变造影剂的使用条件的扫描场景下,将管电压、能量仓的数量、能量仓阈值的设定保持不变(即使用历史数值),扫描参数变化量ΔS可以为管电流的变化量ΔmA、单view下的数据对应的积分时间的变化量Δt、扫描层厚变化量等中的一个或多个的组合。又如,在需要改变造影剂使用条件的扫描场景下,可以通过迭代方式,按预设的参数调节优先级进行计算,确定最佳扫描参数调节组合为扫描参数变化量ΔS。
在一些实施例中,可以通过保证探测器各能量仓的输出数据特性与上次扫描结果可比拟,确定扫描参数变化量。在一些实施例中,可以在保证探测器的各能量仓的输出结果一致的前提下,获得衰减信息与扫描参数之间的对应关系。例如,可以通过在保证探测器的能量仓的输出结果一致的前提下,对多个不同衰减能力的人体的3D模型进行仿真模拟扫描测试,确定衰减信息与扫描参数之间的对应关系(例如,3D模型的厚度与管电流之间的对应关系)。在一些实施例中,可以根据衰减信息与扫描参数之间的对应关系,将信息变化量转化为扫描参数变化量。
在一些实施例中,可以通过迭代计算,将信息变化量转化为扫描参数变化量。例如,可以将目标对象的信息变化量作为输入,进行扫描参数组迭代计算,以确定扫描参数变化量。
在一些实施例中,可以基于扫描参数变化量,调整目标对象的历史扫描参数,以确定目标对象的当前扫描参数。例如,可以基于扫描参数变化量,调整相应类型的扫描参数的数值,以确定目标对象的当前扫描参数。
步骤450,基于信息变化量和历史处理参数,确定当前处理参数。在一些实施例中,步骤450可以由处理设备120或确定模块220执行。
在一些实施例中,当信息变化量小于或等于预设阈值时,可以将目标对象的历史处理参数确定为当前处理参数。
在一些实施例中,当信息变化量大于预设阈值时,可以通过调整历史处理参数,确定目标对象的当前处理参数。在一些实施例中,可以将信息变化量转化为处理参数变化量,基于处理参数变化量和历史处理参数确定当前处理参数。
在一些实施例中,当前处理参数可以包括各能量仓的加权因子、降噪算法、降噪等级、滤波函数和重建矩阵等中的一个或多个。
在一些实施例中,可以基于探测器各能量仓的输出数据特性,确定当前处理参数。例如,当前扫描对应的探测器的某能量仓的输出数据的信噪比,相较于与历史处理参数对应的上一次扫描的该能量仓的输出数据的信噪比较差,则可以降低该能量仓对应的加权因子。
在一些实施例中,无法确定目标对象的信息变化量时,可以基于目标对象的当前基本信息和/或当前扫描的探测器输出结果,确定当前处理参数。例如,患者第一次进行CT扫描时,可以基于当前扫描的探测器输出结果确定该患者的当前处理参数。
影响物质定量分解结果准确性的主要因素为扫描参数中的Bin数和各Bin数据的信噪比程度。其它类别的扫描参数以及处理参数主要通过影响Bin数和各Bin数据的信噪比程度,影响物质定量分解结果的准确性。例如,Bin阈值与管电压、管电流组合决定分布在各Bin内光子数的多少,从而影响信噪比。又如,重建矩阵的大小决定了对相邻像素数据加和处理程度的大小,从而影响信噪比;降噪算法、降噪等级和滤波函数,可以在算法层面改变信噪比。因此,通过将使得各Bin数据的信噪比与上次扫描结果可比拟作为条件,确定当前扫描参数和/或当前处理参数,可以提高物质定量分解结果准确性。
应当注意的是,上述有关方法400的描述仅仅是为了示例和说明,而不限定本说明书的适用范围。对于本领域技术人员来说,在本说明书的指导下可以对方法400进行各种修正和改变。例如,可以先执行步骤450再执行步骤440,或同时执行步骤440和步骤450,或仅执行步骤440和步骤450中的其中一个。然而,这些修正和改变仍在本说明书的范围之内。
图5是根据本说明书另一些实施例所示的示例性成像方法的流程示意图。
如图5所示,在一些实施例中,可以基于目标对象的当前扫描参数和当前处理参数确定目标图像,并基于目标图像和历史图像确定对比结果。在一些实施例中,成像方法500可以由成像设备110或处理设备120执行。例如,成像方法500可以以程序或指令的形式存储在存储装置(如存储设备140)中,当成像设备110或处理设备120执行该程序或指令时,可以实现成像方法500。在一些实施例中,成像方法500可以由成像系统200执行。
步骤510,基于当前扫描参数或预设扫描参数,对目标对象进行扫描,获取探测器输出结果。在一些实施例中,步骤510可以由处理设备120或扫描模块230执行。
在一些实施例中,可以根据当前扫描参数或预设扫描参数,进行参数设置,以对目标对象进行扫描,获取成像设备的探测器输出结果。在一些实施例中,可以由系统自动设置扫描参数,或由用户手动调节扫描参数。例如,系统100可以自动控制成像设备110将相应参数设置为与当前扫描参数匹配的数值,或由医护人员通过终端130或成像设备110的输入装置手动输入相应扫描参数的数值。
在一些实施例中,探测器输出结果可以包括各能量仓的光子计数值。
步骤520,基于当前处理参数,对探测器输出结果进行处理,获取目标对象的目标图像。在一些实施例中,步骤520可以由处理设备120或数据处理模块240执行。
在一些实施例中,可以根据当前处理参数中的一种或多种,对探测器输出结果进行处理,获取目标对象的目标图像。例如,可以根据确定的当前各能量仓的加权因子,对探测器输出的各能量仓的输出数据进行加权计算,并通过重建矩阵进行重建计算,以获取目标对象的目标图像。又如,可以根据确定的当前各能量仓的加权因子对探测器输出的各能量仓的输出数据进行加权计算,同时根据降噪等级利用确定的降噪算法进行降噪计算,通过滤波函数、重建矩阵进行计算,以获取目标对象的目标图像。
在一些实施例中,目标图像可以包括目标对象的物质分解图像。通过根据在不同能量下获得的材料衰减系数的差异,对探测器的输出数据进行基材料分解获得的图像,即为物质分解图像,也称基物质图像。例如,可以分解获得“水-碘”图、“水-钙”图、“水-碘-钙”图、“水-钙-金属”图、脂肪图、尿酸图等。
步骤530,从数据库中获取目标对象的历史图像。在一些实施例中,步骤530可以由处理设备120或获取模块210执行。
在一些实施例中,目标对象的目标图像与历史图像具有相同或相近的置信度。在一些实施例中,目标对象的目标物质分解图像结果与历史物质分解图像结果具有相同的置信度。其中,置信度可以反映物质分解图像中所示定量结果的准确性。例如,目标图像与历史图像在背景区域有着相同或相近的图像质量,此时可认为目标图像与历史图像具有相同的置信度,其所对应的物质分解定量结果也有着相同的准确性。
在一些实施例中,可以从目标对象的个人数据库中获取历史图像。在一些实施例中,历史图像可以包括目标对象上次扫描获得的重建图像。在一些实施例中,历史图像可以包括目标对象先前多次扫描获得的多张重建图像。
步骤540,将目标对象的目标图像与历史图像进行对比,并输出对比结果。在一些实施例中,步骤540可以由处理设备120或数据处理模块240执行。
对比结果可以反映当前图像与历史图像之间对比产生的物质分解定量分析结果。在一些实施例中,对比结果可以包括以下中的至少一种:目标对象的当前基物质图像与历史基物质图像的对比结果、当前各基物质按照第一权重组合获得的当前组合图像与历史各基物质按照第二权重组合获得的历史组合图像的对比结果。其中,第一权重组合使得当前组合图像与历史组合图像中的非目标区域的数值结果保持一致。
在一些实施例中,可以通过对目标图像和历史图像进行基物质分解,确定相应的基物质图像,进而获得当前基物质图像与历史基物质图像的对比结果。在一些实施例中,可以从数据库中获取目标对象的历史基物质图像。在一些实施例中,基物质图像之间的对比结果可以包括,当前基物质图像的数值(例如,含量值)与历史基物质图像的数值(例如,含量值),和/或两组数值之间的差异。
第一权重组合可以反映当前基物质图像对应的各基物质之间的加权因子。相应地,第二权重组合可以反映历史基物质图像对应的各基物质之间的加权因子。在一些实施例中,第一权重组合和第二权重组合可以相同或不同。
在一些实施例中,可以按照第一权重组合,对当前的各基物质进行计算,获得当前组合图像。在一些实施例中,可以按照第二权重组合,对历史的各基物质进行组合,获得历史组合图像。在一些实施例中,当前组合图像与历史组合图像的对比结果可以包括图像、数值、差异等,或其任意组合。例如,可以将当前组合图像与历史组合图像输出展示给医护人员,或将当前组合图像与历史组合图像之间的差异输出展示给医护人员。
非目标区域可以指不需要关注或不需要重点关注的组织区域,例如,正常组织区域、背景区域等。在一些实施例中,当前组合图像与历史组合图像中的非目标区域的数值结果保持一致可以指,当前组合图像与历史组合图像中的非目标区域的CT值或其他等效的数值之间保持一致。通过保证当前组合图像与历史组合图像中的非目标区域的数值结果一致,可以凸显患者的目标区域(例如,病灶区域、目标对象区域)的变化,从而提高诊断效率和诊断结果准确性。
应当注意的是,上述有关方法500的描述仅仅是为了示例和说明,而不限定本说明书的适用范围。对于本领域技术人员来说,在本说明书的指导下可以对方法500进行各种修正和改变。然而,这些修正和改变仍在本说明书的范围之内。
本说明书实施例可能带来的有益效果包括但不限于:(1)通过基于目标对象的当前基本信息和参考基本信息、历史扫描参数、历史处理参数等,确定目标对象的当前扫描参数和/或当前处理参数,可以使得同一患者每次扫描的扫描参数、处理参数之间具有一定的关联性;(2)通过基于目标对象的历史扫描数据确定当前扫描数据,可以帮助患者在一系列扫描中的每次临床结果具有定量的可对比性,从而为医生判断患者病情状态提供可信度更高的数据支撑;(3)为每个患者建立存储历次诊断相关数据的个人数据库,可以帮助对患者在每次扫描时的扫描状态进行比对,以及确定更优的扫描参数和处理参数,以保证患者前后扫描结果具有定量可比性;(4)将目标对象的目标图像和历史图像进行对比,并输出比较结果,可以帮助医生快速了解患者前后扫描之间的状态变化,帮助对病情的诊断和治疗。
需要说明的是,不同实施例可能产生的有益效果不同,在不同的实施例里,可能产生的有益效果可以是以上任意一种或几种的组合,也可以是其他任何可能获得的有益效果。
上文已对基本概念做了描述,显然,对于本领域技术人员来说,上述详细披露仅仅作为示例,而并不构成对本说明书的限定。虽然此处并没有明确说明,本领域技术人员可能会对本说明书进行各种修改、改进和修正。该类修改、改进和修正在本说明书中被建议,所以该类修改、改进、修正仍属于本说明书示范实施例的精神和范围。
同时,本说明书使用了特定词语来描述本说明书的实施例。如“一个实施例”、“一实施例”、和/或“一些实施例”意指与本说明书至少一个实施例相关的某一特征、结构或特点。因此,应强调并注意的是,本说明书中在不同位置两次或多次提及的“一实施例”或“一个实施例”或“一个替代性实施例”并不一定是指同一实施例。此外,本说明书的一个或多个实施例中的某些特征、结构或特点可以进行适当的组合。
此外,除非权利要求中明确说明,本说明书所述处理元素和序列的顺序、数字字母的使用、或其他名称的使用,并非用于限定本说明书流程和方法的顺序。尽管上述披露中通过各种示例讨论了一些目前认为有用的发明实施例,但应当理解的是,该类细节仅起到说明的目的,附加的权利要求并不仅限于披露的实施例,相反,权利要求旨在覆盖所有符合本说明书实施例实质和范围的修正和等价组合。例如,虽然以上所描述的系统组件可以通过硬件设备实现,但是也可以只通过软件的解决方案得以实现,如在现有的服务器或移动设备上安装所描述的系统。
同理,应当注意的是,为了简化本说明书披露的表述,从而帮助对一个或多个发明实施例的理解,前文对本说明书实施例的描述中,有时会将多种特征归并至一个实施例、附图或对其的描述中。但是,这种披露方法并不意味着本说明书对象所需要的特征比权利要求中提及的特征多。实际上,实施例的特征要少于上述披露的单个实施例的全部特征。
一些实施例中使用了描述成分、属性数量的数字,应当理解的是,此类用于实施例描述的数字,在一些示例中使用了修饰词“大约”、“近似”或“大体上”来修饰。除非另外说明,“大约”、“近似”或“大体上”表明所述数字允许有±20%的变化。相应地,在一些实施例中,说明书和权利要求中使用的数值参数均为近似值,该近似值根据个别实施例所需特点可以发生改变。在一些实施例中,数值参数应考虑规定的有效数位并采用一般位数保留的方法。尽管本说明书一些实施例中用于确认其范围广度的数值域和参数为近似值,在具体实施例中,此类数值的设定在可行范围内尽可能精确。
针对本说明书引用的每个专利、专利申请、专利申请公开物和其他材料,如文章、书籍、说明书、出版物、文档等,特此将其全部内容并入本说明书作为参考。与本说明书内容不一致或产生冲突的申请历史文件除外,对本说明书权利要求最广范围有限制的文件(当前或之后附加于本说明书中的)也除外。需要说明的是,如果本说明书附属材料中的描述、定义、和/或术语的使用与本说明书所述内容有不一致或冲突的地方,以本说明书的描述、定义和/或术语的使用为准。
最后,应当理解的是,本说明书中所述实施例仅用以说明本说明书实施例的原则。其他的变形也可能属于本说明书的范围。因此,作为示例而非限制,本说明书实施例的替代配置可视为与本说明书的教导一致。相应地,本说明书的实施例不仅限于本说明书明确介绍和描述的实施例。
Claims (10)
1.一种光子计数CT的成像方法,其特征在于,包括:
获取目标对象的当前基本信息;
获取所述目标对象的参考基本信息和历史参数,所述历史参数包括历史扫描参数和历史处理参数中的至少一个;
基于所述当前基本信息、所述参考基本信息以及所述历史参数,确定所述目标对象的当前扫描参数和/或当前处理参数。
2.根据权利要求1所述的方法,其特征在于,所述方法进一步包括:
基于所述当前扫描参数或预设扫描参数,对所述目标对象进行扫描,获取探测器输出结果;
基于所述当前处理参数,对探测器输出结果进行处理,获取所述目标对象的目标图像。
3.根据权利要求2所述的方法,其特征在于,所述参考基本信息和所述历史参数来自所述目标对象的个人数据库,所述方法进一步包括:
从所述数据库中获取所述目标对象的历史图像;
将所述目标对象的所述目标图像与所述历史图像进行对比,并输出对比结果。
4.根据权利要求3所述的方法,其特征在于,所述目标图像包括所述目标对象的物质分解图像,所述目标图像与所述历史图像具有相同或相近的置信度,所述置信度反映物质分解图像中所示定量结果的准确性。
5.根据权利要求3所述的方法,其特征在于,所述对比结果包括以下中的至少一种:所述目标对象的当前基物质图像与历史基物质图像的对比结果、当前各基物质按照第一权重组合获得的当前组合图像与历史各基物质按照第二权重组合获得的历史组合图像的对比结果,其中,所述第一权重组合使得所述当前组合图像与所述历史组合图像中的非目标区域的数值结果保持一致。
6.根据权利要求1所述的方法,其特征在于,基于所述当前基本信息、所述参考基本信息以及所述历史参数,确定所述目标对象的当前扫描参数和/或当前处理参数,包括:
将所述目标对象的所述当前基本信息与所述参考基本信息进行比较,获取信息变化量,所述信息变化量与所述目标对象的衰减信息相关;
基于所述信息变化量和所述历史参数,确定所述目标对象的当前扫描参数和/或当前处理参数。
7.根据权利要求6所述的方法,其特征在于,所述基于所述信息变化量和所述历史参数,确定所述目标对象的当前扫描参数和/或当前处理参数,包括:
将所述信息变化量转化为扫描参数变化量,并基于所述扫描参数变化量和所述历史扫描参数,确定所述目标对象的当前扫描参数;和/或
将所述信息变化量转化为处理参数变化量,并基于所述处理参数变化量和所述历史处理参数,确定所述目标对象的当前处理参数。
8.根据权利要求1所述的方法,其特征在于,所述扫描参数包括以下中的至少一种:管电压、管电流、能量仓的数量、各能量仓的阈值、单视野下的数据对应的积分时间、机架转速和扫描层厚;所述处理参数包括以下中的至少一种:各能量仓的加权因子、降噪算法、降噪等级、滤波函数和重建矩阵。
9.一种光子计数CT的成像系统,其特征在于,包括:
获取模块,用于获取目标对象的当前基本信息,以及获取所述目标对象的参考基本信息和历史参数,所述历史参数包括历史扫描参数和历史处理参数中的至少一个;
确定模块,用于基于所述当前基本信息、所述参考基本信息以及所述历史参数,确定所述目标对象的当前扫描参数和/或当前处理参数;
扫描模块,用于基于所述当前扫描参数或预设扫描参数,对所述目标对象进行扫描,以获取探测器输出结果;
数据处理模块,用于基于所述当前处理参数,对所述探测器输出结果进行处理,以获取所述目标对象的目标图像。
10.一种计算机可读存储介质,所述存储介质存储计算机指令,当计算机读取所述计算机指令时,所述计算机执行如权利要求1-8中任一项所述的方法。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210737557.8A CN117338315A (zh) | 2022-06-27 | 2022-06-27 | 一种光子计数ct的成像方法和系统 |
PCT/CN2023/102968 WO2024002109A1 (en) | 2022-06-27 | 2023-06-27 | Methods and systems for imaging |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210737557.8A CN117338315A (zh) | 2022-06-27 | 2022-06-27 | 一种光子计数ct的成像方法和系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117338315A true CN117338315A (zh) | 2024-01-05 |
Family
ID=89359931
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210737557.8A Pending CN117338315A (zh) | 2022-06-27 | 2022-06-27 | 一种光子计数ct的成像方法和系统 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN117338315A (zh) |
WO (1) | WO2024002109A1 (zh) |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120271840A1 (en) * | 2011-04-25 | 2012-10-25 | General Electric Company | Systems and methods for storing and providing scan protocol information |
WO2013050897A1 (en) * | 2011-10-04 | 2013-04-11 | Koninklijke Philips Electronics N.V. | Adaptive dual-pass targeted reconstruction and acquisition |
US9795356B2 (en) * | 2013-12-18 | 2017-10-24 | General Electric Company | System and method of X-ray dose distribution for computed tomography based on simulation |
US9869641B2 (en) * | 2016-04-08 | 2018-01-16 | Ellumen, Inc. | Microwave imaging device |
US11497459B2 (en) * | 2018-01-26 | 2022-11-15 | General Electric Company | Methods and system for optimizing an imaging scan based on a prior scan |
CN110464326B (zh) * | 2019-08-19 | 2022-05-10 | 上海联影医疗科技股份有限公司 | 一种扫描参数推荐方法、系统、装置及存储介质 |
EP3889969A1 (en) * | 2020-04-02 | 2021-10-06 | Koninklijke Philips N.V. | Medical imaging system |
CN112690810B (zh) * | 2020-12-22 | 2023-08-15 | 上海联影医疗科技股份有限公司 | 基于先验信息的扫描方法和医学扫描系统 |
-
2022
- 2022-06-27 CN CN202210737557.8A patent/CN117338315A/zh active Pending
-
2023
- 2023-06-27 WO PCT/CN2023/102968 patent/WO2024002109A1/en unknown
Also Published As
Publication number | Publication date |
---|---|
WO2024002109A1 (en) | 2024-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3143935B1 (en) | Tomography apparatus and method of reconstructing tomography images | |
US10628972B2 (en) | Diagnostic imaging method and apparatus, and recording medium thereof | |
Ning et al. | Flat panel detector-based cone-beam volume CT angiography imaging: system evaluation | |
CN110751702B (zh) | 图像重建方法、系统、装置及存储介质 | |
US10638993B2 (en) | Computed tomography apparatus and control method for the same | |
JP4344191B2 (ja) | 画像形成システムの低線量画像シミュレーションのための方法及びシステム | |
US10143433B2 (en) | Computed tomography apparatus and method of reconstructing a computed tomography image by the computed tomography apparatus | |
KR20140019631A (ko) | 방사선 영상 촬영 장치 및 방법 | |
US11497459B2 (en) | Methods and system for optimizing an imaging scan based on a prior scan | |
US10052077B2 (en) | Tomography imaging apparatus and method | |
CN111374690A (zh) | 医学成像方法及系统 | |
US20140072096A1 (en) | Contrast-enhanced imaging of objects | |
CN103860185B (zh) | 多能量图像的确定 | |
CN112842370A (zh) | 用于x射线成像中的参数噪声调制的方法和系统 | |
CN111938681A (zh) | 用于双能量成像中的自动管电势选择的系统和方法 | |
KR20160066941A (ko) | 의료 영상 처리 장치 및 그에 따른 의료 영상 처리 방법 | |
US9858688B2 (en) | Methods and systems for computed tomography motion compensation | |
US10552992B2 (en) | Poly-energetic reconstruction method for metal artifacts reduction | |
US10383589B2 (en) | Direct monochromatic image generation for spectral computed tomography | |
CN109924998A (zh) | 医学成像方法及光子计数能谱ct成像设备 | |
CN117338315A (zh) | 一种光子计数ct的成像方法和系统 | |
Qiu | The Principle and State-of-art Applications for CT Detector | |
KR101590874B1 (ko) | 의료 영상 장치의 영상 재구성 방법 및 그 의료 영상 장치 | |
WO2024002125A1 (en) | Methods and systems for determining scanning parameters | |
CN110730977B (zh) | 低剂量成像方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |