CN103705258B - 一种成像设备的自动曝光控制方法及装置 - Google Patents
一种成像设备的自动曝光控制方法及装置 Download PDFInfo
- Publication number
- CN103705258B CN103705258B CN201310585434.8A CN201310585434A CN103705258B CN 103705258 B CN103705258 B CN 103705258B CN 201310585434 A CN201310585434 A CN 201310585434A CN 103705258 B CN103705258 B CN 103705258B
- Authority
- CN
- China
- Prior art keywords
- exposure
- optical signal
- image
- parameter
- gray value
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000003384 imaging method Methods 0.000 title claims abstract description 28
- 230000008569 process Effects 0.000 claims abstract description 16
- 230000003287 optical effect Effects 0.000 claims description 144
- 230000002085 persistent effect Effects 0.000 claims description 29
- 238000001228 spectrum Methods 0.000 claims description 18
- 238000012937 correction Methods 0.000 claims description 13
- 239000000284 extract Substances 0.000 claims description 7
- 238000009499 grossing Methods 0.000 claims description 7
- 238000011946 reduction process Methods 0.000 claims description 7
- 230000001360 synchronised effect Effects 0.000 claims description 4
- 238000004148 unit process Methods 0.000 claims description 3
- 230000011664 signaling Effects 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 3
- 230000005855 radiation Effects 0.000 abstract description 2
- 230000001105 regulatory effect Effects 0.000 abstract 1
- 238000007689 inspection Methods 0.000 description 21
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000006698 induction Effects 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- PCTMTFRHKVHKIS-BMFZQQSSSA-N (1s,3r,4e,6e,8e,10e,12e,14e,16e,18s,19r,20r,21s,25r,27r,30r,31r,33s,35r,37s,38r)-3-[(2r,3s,4s,5s,6r)-4-amino-3,5-dihydroxy-6-methyloxan-2-yl]oxy-19,25,27,30,31,33,35,37-octahydroxy-18,20,21-trimethyl-23-oxo-22,39-dioxabicyclo[33.3.1]nonatriaconta-4,6,8,10 Chemical compound C1C=C2C[C@@H](OS(O)(=O)=O)CC[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@H]([C@H](C)CCCC(C)C)[C@@]1(C)CC2.O[C@H]1[C@@H](N)[C@H](O)[C@@H](C)O[C@H]1O[C@H]1/C=C/C=C/C=C/C=C/C=C/C=C/C=C/[C@H](C)[C@@H](O)[C@@H](C)[C@H](C)OC(=O)C[C@H](O)C[C@H](O)CC[C@@H](O)[C@H](O)C[C@H](O)C[C@](O)(C[C@H](O)[C@H]2C(O)=O)O[C@H]2C1 PCTMTFRHKVHKIS-BMFZQQSSSA-N 0.000 description 1
- 208000019155 Radiation injury Diseases 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000004980 dosimetry Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000008521 reorganization Effects 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- WFKWXMTUELFFGS-UHFFFAOYSA-N tungsten Chemical compound [W] WFKWXMTUELFFGS-UHFFFAOYSA-N 0.000 description 1
- 229910052721 tungsten Inorganic materials 0.000 description 1
- 239000010937 tungsten Substances 0.000 description 1
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种成像设备的自动曝光控制方法及装置。本发明采用了两次曝光过程,并对预曝光图像进行处理,对受检部位组织区域的图像灰度值进行统计,将灰度均值与推荐值的比值作为计算二次曝光参数(mA值或者曝光时间值)的主要依据,同时利用灰度最大值和最小值对调整后的图像灰度范围进行限定,确保二次曝光能够获得较高质量的图像,同时能够降低受检者所受的X射线辐射。而且本发明采取了比信噪比更直接的图像灰度值统计量进行曝光参数估计,对于各种图像均可以获得较好的处理结果,提高了自动曝光控制的准确性。
Description
技术领域
本发明涉及成像设备的自动曝光控制,尤其涉及一种基于两次曝光的控制方法及装置。
背景技术
自动曝光控制系统在基于X射线进行成像诊断的设备中,目的在于保证图像质量的同时,尽量降低X射线曝光的剂量,对受检者进行保护。在X射线源端,可以对曝光剂量进行控制的参数有:1)球管阴极与阳极之间的偏置电压;常用kV(千伏)值表示;2)球管的管电流,常用mA(毫安)值表示;3)一次曝光持续时间,常用时间值s(秒)表示。考虑到临床受检者的体型存在很大差异(如肥胖型、瘦小型、婴幼儿等),而且每次检查的身体部位也可能不同,所以单纯依靠临床物理师凭经验设置的曝光参数进行成像,其曝光剂量往往不是最优的。带来的影响是:如果曝光剂量过高,则受检者将承受过多的X射线照射,罹患重大疾病的风险增大;而如果曝光剂量过低,获得的图像质量较差,会严重影响医生的临床诊断。因此,实现智能化的曝光剂量控制对于X射线来说至关重要。
目前已经有至少两种实现自动曝光控制的方法。其中一种是通过在成像探测器表面设置电离室,可以使用一个或者多个电离室模块,通过电离室感应X射线直射以及穿过受检者特定部位后的强度变化,从而对采集图像的增益进行校正。另外一种不通过增加电离室等硬件模块,而单纯通过软件方式实现自动曝光剂量的控制,典型的如专利US6459765所示。该方法增加一次低剂量的预曝光过程,在计算机端对获取的预曝光图像进行处理,计算受检者身体部位的信噪比,通过将该计算值与系统推荐的理想信噪比进行比较,计算二者的比值,并将该比值作用于预曝光设定的mA值或者曝光时间值s上,实现对mA值或者曝光时间值s的调节,然后实施第二次曝光并采集图像。
上述两种自动曝光控制方法各有不足。采用电离室的方法缺点在于,这些电离室覆盖于成像探测器前表面,造成了X射线的额外衰减,会影响电离室覆盖区域的探测器单元对X射线的捕获。这就要求电离室模块应做的足够小,或者对X射线的衰减尽量小,以减少对探测器成像的影响。这样的要求会使得电离室的硬件制造成本增加,并且也会影响电离室对X射线探测的精度。而专利US6459765的方法缺点在于,采用信噪比作为评价标准,不仅要统计图像中组织区域的信号强度,还要统计非组织区域的噪声水平,从而实现对预曝光图像信噪比的评价。一旦组织区域的投影图像覆盖整个探测器表面,对噪声的估计将比较困难,会影响自动曝光控制的准确性。
上述问题也同样出现在其他相类似的成像设备的自动曝光控制系统中。
发明内容
本发明的目的在于提供一种成像设备的自动曝光控制方法及装置,以解决上述问题,提高自动曝光控制的准确性。
基于上述目的,本发明实施例提供了一种成像设备的自动曝光控制方法,该方法包括:
控制单元配置预曝光参数并将所述预曝光参数发送至光信号发生器;
所述光信号发生器根据所述预曝光参数产生预曝光光信号;
光信号探测器根据所述预曝光光信号形成预曝光图像并将所述预曝光图像发送至所述控制单元;
所述控制单元对所述预曝光图像进行处理,统计所述预曝光图像中组织区域的灰度值的最大值Imax、最小值Imin以及均值Imean,计算预设的组织区域的灰度值的均值Imean'与所述预曝光图像中组织区域的灰度值的均值Imean的比值R1,并根据所述预设的组织区域的灰度值的最大值Imax'与所述预曝光图像的灰度值的最大值Imax的比值和所述预设的组织区域的灰度值的最小值Imin'与所述预曝光图像中组织区域的灰度值的最小值Imin的比值修正所述比值R1为比值R2,根据所述比值R2确定二次曝光参数并发送所述二次曝光参数至光信号发生器;
所述光信号发生器根据所述二次曝光参数产生二次曝光光信号;
所述光信号探测器根据所述二次曝光光信号形成二次曝光图像并将所述二次曝光图像发送至所述控制单元;
所述预曝光参数和所述二次曝光参数均至少包括所述光信号发生器的偏置电压、所述光信号发生器的电流和所述光信号发生器的一次曝光持续时间。
优选的,根据所述预曝光参数得到的预曝光计量为正常曝光量的1/16~1/50。
优选的,所述对所述预曝光图像进行处理,统计所述预曝光图像的灰度值的最大值、最小值以及均值包括:
对所述预曝光图像进行平滑降噪处理;
提取所述预曝光图像的组织区域的边界轮廓;
统计所述边界轮廓内的图像的灰度值的最大值、最小值和均值。
优选的,所述根据所述预设的组织区域的灰度值的最大值Imax'与所述预曝光图像中组织区域的灰度值的最大值Imax的比值和所述预设的组织区域的灰度值的最小值Imin'与所述预曝光图像中组织区域的灰度值的最小值Imin的比值修正所述比值R1为比值R2包括:
若R1×Imax>Imax'或R1×Imin<Imin',则调整R1为R2,其中R2×Imax<Imax'且R2×Imin>Imin';
所述根据所述比值R2确定二次曝光参数包括根据所述比值R2确定所述光信号发生器的管电流和/或所述光信号发生器的一次曝光持续时间。
优选的,所述根据所述比值R2确定二次曝光参数包括:
调整所述光信号发生器的电流为预曝光参数中的所述光信号发生器的电流的R2倍,将调整后的所述光信号发生器的管电流、预曝光参数中的所述光信号发生器的一次曝光持续时间和预曝光参数中的所述光信号发生器的偏置电压一起确定为所述二次曝光参数;
或;
调整所述光信号发生器的一次曝光持续时间为预曝光参数中的所述光信号发生器的一次曝光持续时间的R2倍,将调整后的所述光信号发生器的一次曝光持续时间、预曝光参数中的所述光信号发生器的管电流和预曝光参数中的所述光信号发生器的偏置电压一起确定为所述二次曝光参数。
本发明实施例还提供了一种成像设备的自动曝光控制装置,该装置包括:光信号发生器、光信号探测器和控制单元;
所述控制单元,用于配置预曝光参数并将所述预曝光参数发送至所述光信号发生器;
所述控制单元,还用于对所述光信号探测器发送的预曝光图像进行处理,统计所述预曝光图像中组织区域的灰度值的最大值Imax、最小值Imin以及均值Imean,计算预设的组织区域的灰度值的均值Imean'与所述预曝光图像中组织区域的灰度值的均值Imean的比值R1,并根据所述预设的组织区域的灰度值的最大值Imax'与所述预曝光图像中组织区域的灰度值的最大值Imax的比值和所述预设的组织区域的灰度值的最小值Imin'与所述预曝光图像中组织区域的灰度值的最小值Imin的比值修正所述比值R1为比值R2,根据所述比值R2确定二次曝光参数并发送所述二次曝光参数至光信号发生器;
所述光信号发生器,用于根据所述预曝光参数产生预曝光光信号;
所述光信号发生器,还用于根据所述二次曝光参数产生二次曝光光信号;
所述光信号探测器,用于根据所述预曝光光信号形成预曝光图像并将所述预曝光图像发送至所述控制单元;
所述光信号探测器,还用于根据所述二次曝光光信号形成所述二次曝光图像并将所述二次曝光图像发送至所述控制单元;
所述预曝光参数和所述二次曝光参数均至少包括所述光信号发生器的偏置电压、所述光信号发生器的管电流和所述光信号发生器的一次曝光持续时间。
优选的,根据所述预曝光参数得到的所述预曝光计量为正常曝光量的1/16~1/50。
优选的,所述控制单元用于对所述预曝光图像进行平滑降噪处理,提取所述预曝光图像的组织区域的边界轮廓并统计所述边界轮廓内的图像的灰度值的最大值、最小值和均值。
优选的,所述控制单元用于在R1×Imax>Imax'或R1×Imin<Imin'时,调整R1为R2,并根据所述比值R2确定所述光信号发生器的电流和/或所述光信号发生器的一次曝光持续时间;其中R2×Imax<Imax'且R2×Imin>Imin'。
优选的,所述光信号发生器包括高压发生器和球管;所述控制单元包括曝光控制盒、曝光手闸和图像工作站;所述光信号探测器为平板X射线探测器;
其中,所述图像工作站,用于配置所述预曝光参数并将所述预曝光参数发送至所述曝光控制盒;
所述图像工作站,还用于对所述平板X射线探测器发送的预曝光图像进行处理,统计所述预曝光图像中组织区域的灰度值的最大值Imax、最小值Imin以及均值Imean,计算预设的组织区域的灰度值的均值Imean'与所述预曝光图像中组织区域的灰度值的均值Imean的比值R1,并根据所述预设的组织区域的灰度值的最大值Imax'与所述预曝光图像中组织区域的灰度值的最大值Imax的比值和所述预设的组织区域的灰度值的最小值Imin'与所述预曝光图像中组织区域的灰度值的最小值Imin的比值修正所述比值R1为比值R2,根据所述比值R2确定二次曝光参数并发送所述二次曝光参数至所述曝光控制盒;
所述曝光手闸,用于手动产生曝光使能信号并发送至所述曝光控制盒;
所述曝光控制盒,用于将所述预曝光参数、所述二次曝光参数和所述曝光使能信号发送至所述高压发生器并控制所述平板X射线探测器同步采集图像;
所述高压发生器,用于在接收到所述曝光使能信号时,根据所述预曝光参数或所述二次曝光参数为所述球管提供偏置电压以及管电流;
所述球管,用于根据所述预曝光参数或所述二次曝光参数产生X射线;
所述平板X射线探测器用于在接收到所述曝光使能信号时生成图像并将所述图像传输至所述图像工作站。
本发明的有益效果是:
本发明通过对预曝光图像进行处理,对受检部位组织区域的图像灰度值进行统计,将灰度均值与推荐值的比值作为计算二次曝光mA值或者曝光时间值的主要依据,同时利用灰度最大值和最小值对调整后的图像灰度范围进行限定,确保二次曝光能够获得较高质量的图像,同时能够降低受检者所受的X射线辐射。而且本发明采取了比信噪比更直接的图像灰度值统计量进行曝光参数估计,对于各种图像均可以获得较好的处理结果,解决了专利US6459765在组织区域覆盖整幅图像时不能准确估计噪声水平、造成信噪比计算可能出现较大误差的问题,提高了自动曝光控制的准确性。
附图说明
图1为基于X射线成像系统的自动曝光控制装置结构图;
图2为基于X射线成像系统的自动曝光控制方法流程图;
图3为图2中步骤S13的具体流程图;
图4为典型的X射线成像示意图;
图5提取图像组织区域边界示意图;
图6组织区域的边界轮廓提取效果示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例一提供了一种成像设备的自动曝光控制装置,该装置包括:光信号发生器、光信号探测器和控制单元;
控制单元,用于配置预曝光参数并将预曝光参数发送至光信号发生器。
控制单元,还用于对光信号探测器发送的预曝光图像进行处理,统计预曝光图像中组织区域的灰度值的最大值Imax、最小值Imin以及均值Imean,计算预设的组织区域的灰度值的均值Imean'与所述预曝光图像中组织区域的灰度值的均值Imean的比值R1,并根据所述预设的组织区域的灰度值的最大值Imax'与所述预曝光图像中组织区域的灰度值的最大值Imax的比值和所述预设的组织区域的灰度值的最小值Imin'与所述预曝光图像中组织区域的灰度值的最小值Imin的比值修正所述比值R1为比值R2,根据所述比值R2确定二次曝光参数并发送所述二次曝光参数至光信号发生器。
所述光信号发生器,用于根据所述预曝光参数产生预曝光光信号;
所述光信号发生器,还用于根据所述二次曝光参数产生二次曝光光信号;
所述光信号探测器,用于根据所述预曝光光信号形成预曝光图像并将所述预曝光图像发送至所述控制单元;
所述光信号探测器,还用于根据所述二次曝光光信号形成所述二次曝光图像并将所述二次曝光图像发送至所述控制单元;
所述预曝光参数和所述二次曝光参数均至少包括所述光信号发生器的偏置电压、所述光信号发生器的管电流和所述光信号发生器的一次曝光持续时间。
根据预曝光参数得到的预曝光剂量为正常曝光剂量的1/16~1/50,其中正常曝光量为现有的光信号发生器正常曝光时(只有一次曝光时)的曝光量
所述控制单元具体的用于对所述预曝光图像进行平滑降噪处理,提取所述预曝光图像的组织区域的边界轮廓并统计所述边界轮廓内的图像的灰度值的Imax、最小值Imin以及均值Imean,并在R1×Imax>Imax'或R1×Imin<Imin'时,调整R1为R2,根据比值R2确定所述光信号发生器的电流和/或所述光信号发生器的一次曝光持续时间;其中R2×Imax<Imax'且R2×Imin>Imin'。
其中根据所述比值R2确定二次曝光参数可具体包括:
调整所述光信号发生器的电流为预曝光参数中的所述光信号发生器的电流的R2倍,将调整后的所述光信号发生器的管电流、预曝光参数中的所述光信号发生器的一次曝光持续时间和预曝光参数中的所述光信号发生器的偏置电压一起确定为所述二次曝光参数;
或;
调整所述光信号发生器的一次曝光持续时间为预曝光参数中的所述光信号发生器的一次曝光持续时间的R2倍,将调整后的所述光信号发生器的一次曝光持续时间、预曝光参数中的所述光信号发生器的管电流和预曝光参数中的所述光信号发生器的偏置电压一起确定为所述二次曝光参数。
以下以X射线成像系统的自动曝光控制装置为例对本发明说明。
X射线成像系统的硬件组成结构如图1所示。其中,高压发生器11与球管12构成系统的X射线源组件,用于产生X射线。二者可以是分离器件,通过高压密封线缆连接,也可以整合在一起构成一体化器件。高压发生器11作为球管12的控制器件,通过与控制计算机或者控制器通信,为球管12提供偏置高压以及灯丝电流,并对球管12的工作状态进行监控。球管12作为被控器件,其内部有一真空管,在曝光状态下,管内阴极持续产生热电子,热电子在阴极与阳极之间的偏置高压作用下高速轰击阳极靶(靶面一般为钨材料),从而产生X射线。X射线照射受检者的待检部位,一部分射线被待检部位吸收,造成X射线的衰减,剩余的射线则到达X射线探测器13,实现对受检者待检部位的成像。X射线探测器13则感应穿过受检者的X射线,将其转换为电信号,并形成可被计算机处理的数字化灰度图像。曝光控制盒14用于接收来自图像工作站15以及曝光手闸16的控制指令,实现对X射线源与X射线探测器13的参数配置以及曝光过程控制。曝光手闸16仅用于手动产生曝光使能信号。图像工作站15可以为成像系统配置参数,还可以接收来自X射线探测器13的图像并进行处理。在自动曝光控制系统中,图像工作站15需要根据预曝光图像计算出二次曝光的优化剂量参数。
具体的,图像工作站15,用于配置预曝光参数并将预曝光参数发送至曝光控制盒14,还用于对X射线探测器13发送的预曝光图像进行处理,统计预曝光图像的灰度值的最大值Imax、最小值Imin以及均值Imean,计算预设的灰度值的均值Imean'与预曝光图像的灰度值的均值Imean的比值R1,并根据预设的灰度值的最大值Imax'与预曝光图像的灰度值的最大值Imax的比值和所述预设的灰度值的最小值Imin'与所述预曝光图像的灰度值的最小值Imin的比值修正所述比值R1为比值R2,根据比值R2确定二次曝光参数并发送二次曝光参数至曝光控制盒14。
所述曝光控制盒14,用于将预曝光参数、二次曝光参数和曝光手闸16产生的曝光使能信号发送至高压发生器11并控制X射线探测器13同步采集图像;
高压发生器11,用于在接收到曝光使能信号时,根据预曝光参数或所述二次曝光参数为球管12提供偏置电压以及电流,以便球管12根据预曝光参数或二次曝光参数产生X射线。
X射线探测器13用于在接收到曝光使能信号时按照预曝光参数生成预曝光图像,按照二次曝光参数生成二次曝光图像,并将图像传输至图像工作站15。
基于上述硬件组成结构,本发明所涉及的自动曝光控制方法包括两次曝光成像过程,流程图如图2所示。
S11、预曝光参数配置。在X射线成像系统准备就绪(即系统完成了开机启动及状态自检,球管预热,探测器的日常校正)的情况下,通过图像工作站(或者等效的具有人机交互界面的X射线成像系统控制平台),对受检者的待检部位进行选择。图像工作站基于选定的待检部位自动生成一组曝光参数,曝光参数应至少包括球管的焦点大小、kV值、mA值、曝光时间等。该组参数作为预曝光参数,因此要求mA值小,曝光时间短,产生的X射线剂量应相当于正常曝光量的1/16-1/50,以使受检者所受的辐射伤害尽量小。相应的曝光参数经由曝光控制盒发送给X射线源组件以及探测器,对二者的工作参数进行配置,并使二者处于等待曝光状态。
S12、启动预曝光。临床物理师按下曝光手闸,曝光使能信号传递至曝光控制盒,通过控制盒实现X射线组件与探测器的曝光及数据采集同步,启动X射线源曝光,并控制探测器采集受检者检测部位的投影图像,并将图像传输至图像工作站。
S13、预曝光图像处理。在图像工作站对预曝光图像进行处理,执行的操作包括图像去噪,组织区域边缘检测,组织区域图像灰度值的最大值、最小值和均值统计,并将统计值与临床推荐参数值进行比较,根据比值差异调整mA值或者曝光时间值s,使得二次曝光的X射线剂量最优化。
S14、二次曝光参数设置。根据第三步计算得到的mA值或者曝光时间s,重新生成一组曝光参数,并自动通过曝光控制盒对X射线源组件和探测器进行参数配置,使X射线源组件和探测器处理等待曝光状态。
S15、启动二次曝光。在不需要临床物理师干预的情况下,通过曝光控制盒自动产生曝光使能信号,控制X射线源组件进行二次曝光,并控制探测器采集受检者检测部位的二次投影图像,将改组投影图像传递至图像工作站,供后续处理和临床分析使用。
本发明还提供了一种步骤S13的具体实现方式,具体如图3所示,包括:
S21、对采集的预曝光图像进行校正和平滑降噪处理。其中降噪过程可以使用数字化的低通滤波器进行,目的是去除图像采集过程中引入的噪声,获得较为平滑的图像,方便后续的处理工作。
S22、在滤波的基础上,对图像中受检者组织区域的边界进行识别和提取。
典型的X射线成像结果应如图4a,4b,4c所示,具体说明如下:
图4a,4b,4c中ABCD所包围的矩形区域均为有效图像区域;图4a和图4b中,EFG或者EFGH所包围的规则或不规则区域(阴影部分)为检测部位的组织区域投影图像,X射线穿透受检者检测部位后,所成的图像灰度值相对较小,而且波动较大。图4a和图4b中,白色区域为X射线直射探测器的投影图像,其灰度值在整个图像中较大,但是波动较小,仅受噪声影响;
图4c比较特殊,但是临床上也大量存在这样的X射线图像,即检测部位的组织区域覆盖了整个探测器的成像区域,没有X射线直射到探测器上。
对于图4a和图4b所示的预曝光图像,可以通过一定的计算方法获得一个阈值b0,该阈值应以一定比例小于X射线直射探测器区域的图像的灰度均值。以图4a所示图像为例,对预曝光图像的每一行数据,分别从左右两侧比较其灰度值与设定阈值的大小关系,如图5a,5b所示。当遇到灰度值斜向下穿过阈值线时,对应的列序号即为组织区域边界所在位置。图5a中的I、J点为第r行图像中组织区域的边界,图5b中的I'、J'点为第r行图像中组织区域边界的估计值,由于噪声以及阈值设定的影响,其与I、J点可能存在一定的偏差,但是不会影响边界的变化趋势。所有图像行的I'、J'点分别连接起来,即构成了组织区域的边界。最终获得的组织区域的边界应可能如图6所示,其与理想的组织区域边界(图4a所示)会存在一定差异,但是该差异对大样本图像统计值的影响很小。另外,还可以采用梯度方法以及形态学方法,提高组织区域边界提取的精度,只是计算量会增加。对于图4c所示的预曝光图像,则无需进行组织区域边界的判别。
S23、计算组织区域内图像灰度值的最大值、最小值以及均值,并与系统推荐值比较。受检者各部位图像灰度均值的推荐值可通过临床试验获得,具有广泛适用性。设组织区域内图像灰度值得最大值、最小值和均值分别为Imax、Imin、Imean,系统推荐的图像灰度最大值、最小值和均值分别为Imax'、Imin'、Imean'。计算推荐的灰度均值与预曝光图像中组织区域灰度均值的比值:
R1=Imean'/Imean
如果同时存在R1×Imax<Imax',以及R1×Imin>Imin',则认为比值R1为合理值。否则可适当调整R1为R2,以同时满足R2×Imax<Imax',以及R2×Imin>Imin'。
其中R2优选满足R×Imax<Imax',以及R×Imin>Imin'的所有R值集合中,与R1的绝对值最小的R。
S24、通过比值R2确定mA值或者曝光时间的调整量。设预曝光时设定的mA值为MAps,曝光时间s为Tps,二次曝光的mA值为MAaec,曝光时间为Taec。我们选择调整mA值以及曝光时间的任意一个,则:
若调整mA值,则有:
MAaec=R×MAps
Taec=Tps若调整曝光时间值,则有:
Taec=R×Tps
MAaec=MAps
将MAaec、Taec值与预曝光之前设定的kV值一起,作为二次曝光的控制参数,即可完成自动曝光控制。
以上是对X射线成像系统自动曝光控制方法的介绍,上述方法可同样应用在其他类似的成像系统中。具体包括:
控制单元配置预曝光参数并将所述预曝光参数发送至光信号发生器;
所述光信号发生器根据所述预曝光参数产生预曝光光信号;
光信号探测器根据所述预曝光光信号形成预曝光图像并将所述预曝光图像发送至所述控制单元;
控制单元对所述预曝光图像进行处理,统计所述预曝光图像的灰度值的最大值Imax、最小值Imin以及均值Imean,计算预设的灰度值的均值Imean'与所述预曝光图像的灰度值的均值Imean的比值R1,并根据所述预设的灰度值的最大值Imax'与所述预曝光图像的灰度值的最大值Imax的比值和所述预设的灰度值的最小值Imin'与所述预曝光图像的灰度值的最小值Imin的比值修正所述比值R1为比值R2,根据所述比值R2确定二次曝光参数并发送所述二次曝光参数至光信号发生器;
所述光信号发生器根据所述二次曝光参数产生二次曝光光信号;
所述光信号探测器根据所述二次曝光光信号形成二次曝光图像并将所述二次曝光图像发送至所述控制单元;
所述预曝光参数和所述二次曝光参数均至少包括所述光信号发生器的偏置电压、所述光信号发生器的电流和所述光信号发生器的一次曝光持续时间。
其中,根据预曝光参数得到的预曝光剂量可限制为正常曝光剂量的1/16~1/50。
对预曝光图像进行处理,统计所述预曝光图像的灰度值的最大值、最小值以及均值的过程具体包括:
对所述预曝光图像进行平滑降噪处理;
提取所述预曝光图像的组织区域的边界轮廓;
统计所述边界轮廓内的图像的灰度值的最大值、最小值和均值。
而根据所述预设的灰度值的最大值Imax'与所述预曝光图像的灰度值的最大值Imax的比值和所述预设的灰度值的最小值Imin'与所述预曝光图像的灰度值的最小值Imin的比值修正所述比值R1为比值R2包括:
若R1×Imax>Imax'或R1×Imin<Imin',则调整R1为R2,其中R2×Imax<Imax'且R2×Imin>Imin';
所述根据所述比值R2确定二次曝光参数包括根据所述比值R2确定所述光信号发生器的电流和/或所述光信号发生器的一次曝光持续时间。
根据所述比值R2确定二次曝光参数包括:
调整所述光信号发生器的电流为预曝光参数中的所述光信号发生器的电流的R2倍,将调整后的所述光信号发生器的电流、预曝光参数中的所述光信号发生器的一次曝光持续时间和预曝光参数中的所述光信号发生器的偏置电压一起确定为所述二次曝光参数;
或;
调整所述光信号发生器的一次曝光持续时间为预曝光参数中的所述光信号发生器的一次曝光持续时间的R2倍,将调整后的所述光信号发生器的一次曝光持续时间、预曝光参数中的所述光信号发生器的电流和预曝光参数中的所述光信号发生器的偏置电压一起确定为所述二次曝光参数。
本发明由于不涉及噪声估计的问题,因此对图4c所示的组织区域覆盖整幅的情况同样可以获得较好的处理结果,避免了因组织区域覆盖整幅图像时不能准确估计噪声水平、造成信噪比计算可能出现较大误差,从而影响最终控制效果的问题,而且本发明采取了比信噪比更直接的图像灰度值统计量进行曝光参数估计,相比信噪比的方式处理时间短,因此缩短了预曝光与二次曝光之间的时间间隔,避免引起受检者不适,以及受检者体位移动造成图像质量的下降。确保二次曝光能够获得较高质量的图像,同时能够降低受检者所受的辐射。
需要说明的是,本发明的装置和方法实施例相对应,相关部分可互相参考。
以上的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应注意的是,以上仅为本发明的一个具体实施例而已,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (10)
1.一种成像设备的自动曝光控制方法,其特征在于,所述方法包括:
控制单元配置预曝光参数并将所述预曝光参数发送至光信号发生器;
所述光信号发生器根据所述预曝光参数产生预曝光光信号;
光信号探测器根据所述预曝光光信号形成预曝光图像并将所述预曝光图像发送至所述控制单元;
所述控制单元对所述预曝光图像进行处理,统计所述预曝光图像中组织区域的灰度值的最大值Imax、最小值Imin以及均值Imean,计算预设的组织区域的灰度值的均值Imean'与所述预曝光图像中组织区域的灰度值的均值Imean的比值R1,并根据所述预设的组织区域的灰度值的最大值Imax'与所述预曝光图像中组织区域的灰度值的最大值Imax的比值和所述预设的组织区域的灰度值的最小值Imin'与所述预曝光图像中组织区域的灰度值的最小值Imin的比值修正所述比值R1为比值R2,根据所述比值R2确定二次曝光参数并发送所述二次曝光参数至光信号发生器;
所述光信号发生器根据所述二次曝光参数产生二次曝光光信号;
所述光信号探测器根据所述二次曝光光信号形成二次曝光图像并将所述二次曝光图像发送至所述控制单元;
所述预曝光参数和所述二次曝光参数均至少包括所述光信号发生器的偏置电压、所述光信号发生器的管电流和所述光信号发生器的一次曝光持续时间;
所述根据所述预设的组织区域的灰度值的最大值Imax'与所述预曝光图像中组织区域的灰度值的最大值Imax的比值和所述预设的组织区域的灰度值的最小值Imin'与所述预曝光图像中组织区域的灰度值的最小值Imin的比值修正所述比值R1为比值R2包括:
若R1×Imax>Imax'或R1×Imin<Imin',则调整R1为R2,其中R2×Imax<Imax'且R2×Imin>Imin'。
2.如权利要求1所述的方法,其特征在于,根据所述预曝光参数得到的预曝光剂量为正常曝光剂量的1/16~1/50。
3.如权利要求1所述的方法,其特征在于,所述对所述预曝光图像进行处理,统计所述预曝光图像的灰度值的最大值、最小值以及均值包括:
对所述预曝光图像进行平滑降噪处理;
提取所述预曝光图像的组织区域的边界轮廓;
统计所述边界轮廓内的图像的灰度值的最大值、最小值和均值。
4.如权利要求1所述的方法,其特征在于,
所述根据所述比值R2确定二次曝光参数包括根据所述比值R2确定所述光信号发生器的管电流和/或所述光信号发生器的一次曝光持续时间。
5.如权利要求1所述的方法,其特征在于,所述根据所述比值R2确定二次曝光参数包括:
调整所述光信号发生器的管电流为预曝光参数中的所述光信号发生器的管电流的R2倍,将调整后的所述光信号发生器的管电流、预曝光参数中的所述光信号发生器的一次曝光持续时间和预曝光参数中的所述光信号发生器的偏置电压一起确定为所述二次曝光参数;
或;
调整所述光信号发生器的一次曝光持续时间为预曝光参数中的所述光信号发生器的一次曝光持续时间的R2倍,将调整后的所述光信号发生器的一次曝光持续时间、预曝光参数中的所述光信号发生器的管电流和预曝光参数中的所述光信号发生器的偏置电压一起确定为所述二次曝光参数。
6.一种成像设备的自动曝光控制装置,其特征在于,所述装置包括:光信号发生器、光信号探测器和控制单元;
所述控制单元,用于配置预曝光参数并将所述预曝光参数发送至所述光信号发生器;
所述控制单元,还用于对所述光信号探测器发送的预曝光图像进行处理,统计所述预曝光图像中组织区域的灰度值的最大值Imax、最小值Imin以及均值Imean,计算预设的组织区域的灰度值的均值Imean'与所述预曝光图像中组织区域的灰度值的均值Imean的比值R1,并根据所述预设的组织区域的灰度值的最大值Imax'与所述预曝光图像中组织区域的灰度值的最大值Imax的比值和所述预设的组织区域的灰度值的最小值Imin'与所述预曝光图像中组织区域的灰度值的最小值Imin的比值修正所述比值R1为比值R2,根据所述比值R2确定二次曝光参数并发送所述二次曝光参数至光信号发生器;
所述光信号发生器,用于根据所述预曝光参数产生预曝光光信号;
所述光信号发生器,还用于根据所述二次曝光参数产生二次曝光光信号;
所述光信号探测器,用于根据所述预曝光光信号形成预曝光图像并将所述预曝光图像发送至所述控制单元;
所述光信号探测器,还用于根据所述二次曝光光信号形成所述二次曝光图像并将所述二次曝光图像发送至所述控制单元;
所述预曝光参数和所述二次曝光参数均至少包括所述光信号发生器的偏置电压、所述光信号发生器的电流和所述光信号发生器的一次曝光持续时间;
所述控制单元用于在R1×Imax>Imax'或R1×Imin<Imin'时,调整R1为R2;其中R2×Imax<Imax'且R2×Imin>Imin'。
7.如权利要求6所述的装置,其特征在于,根据所述预曝光参数得到的所述预曝光剂量为正常曝光剂量的1/16~1/50。
8.如权利要求6所述的装置,其特征在于,
所述控制单元用于对所述预曝光图像进行平滑降噪处理,提取所述预曝光图像的组织区域的边界轮廓并统计所述边界轮廓内的图像的灰度值的最大值、最小值和均值。
9.如权利要求6所述的装置,其特征在于,所述控制单元根据所述比值R2确定所述光信号发生器的电流和/或所述光信号发生器的一次曝光持续时间。
10.如权利要求6所述的装置,其特征在于,所述光信号发生器包括高压发生器和球管;所述控制单元包括曝光控制盒、曝光手闸和图像工作站;所述光信号探测器为平板X射线探测器;
其中,所述图像工作站,用于配置所述预曝光参数并将所述预曝光参数发送至所述曝光控制盒;
所述图像工作站,还用于对所述平板X射线探测器发送的预曝光图像进行处理,统计所述预曝光图像中组织区域的灰度值的最大值Imax、最小值Imin以及均值Imean,计算预设的组织区域的灰度值的均值Imean'与所述预曝光图像中组织区域的灰度值的均值Imean的比值R1,并根据所述预设的组织区域的灰度值的最大值Imax'与所述预曝光图像中组织区域的灰度值的最大值Imax的比值和所述预设的组织区域的灰度值的最小值Imin'与所述预曝光图像的灰度值的最小值Imin的比值修正所述比值R1为比值R2,根据所述比值R2确定二次曝光参数并发送所述二次曝光参数至所述曝光控制盒;
所述曝光手闸,用于手动产生曝光使能信号并发送至所述曝光控制盒;
所述曝光控制盒,用于将所述预曝光参数、所述二次曝光参数和所述曝光使能信号发送至所述高压发生器并控制所述平板X射线探测器同步采集图像;
所述高压发生器,用于在接收到所述曝光使能信号时,根据所述预曝光参数或所述二次曝光参数为所述球管提供偏置电压以及灯丝电流;
所述球管,用于根据所述预曝光参数或所述二次曝光参数产生X射线;
所述平板X射线探测器用于在接收到所述曝光使能信号时生成图像并将所述图像传输至所述图像工作站。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310585434.8A CN103705258B (zh) | 2013-11-19 | 2013-11-19 | 一种成像设备的自动曝光控制方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310585434.8A CN103705258B (zh) | 2013-11-19 | 2013-11-19 | 一种成像设备的自动曝光控制方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103705258A CN103705258A (zh) | 2014-04-09 |
CN103705258B true CN103705258B (zh) | 2015-06-24 |
Family
ID=50398895
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310585434.8A Active CN103705258B (zh) | 2013-11-19 | 2013-11-19 | 一种成像设备的自动曝光控制方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103705258B (zh) |
Families Citing this family (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105078497B (zh) * | 2014-05-19 | 2017-12-08 | 锐珂(上海)医疗器材有限公司 | 数字化放射摄影装置与方法 |
CN106308828A (zh) * | 2015-07-01 | 2017-01-11 | 上海奕瑞光电子科技有限公司 | 数字图像探测装置及平板探测器无线外触发方法 |
CN106308827A (zh) * | 2015-07-01 | 2017-01-11 | 南京普爱医疗设备股份有限公司 | 一种三维c形臂自动剂量控制方法 |
CN105361898B (zh) * | 2015-12-11 | 2019-04-12 | 沈阳东软医疗系统有限公司 | 校正参数确定方法、装置及设备 |
CN105636327B (zh) * | 2015-12-28 | 2017-08-25 | 上海奕瑞光电子科技有限公司 | 一种基于透射可见光的曝光控制系统及方法 |
CN106954329B (zh) * | 2016-05-17 | 2018-07-13 | 江苏康众数字医疗设备有限公司 | 成像设备的自动曝光控制方法及曝光系统 |
CN106214171B (zh) * | 2016-09-07 | 2020-04-10 | 东软医疗系统股份有限公司 | 一种自动曝光控制方法和装置 |
CN106419944A (zh) * | 2016-11-09 | 2017-02-22 | 深海精密科技(深圳)有限公司 | C型臂x射线机自动曝光控制方法、控制系统及控制装置 |
CN108606804A (zh) * | 2016-12-13 | 2018-10-02 | 上海西门子医疗器械有限公司 | 自动曝光控制的自动校准系统及方法和医疗设备 |
CN108056784A (zh) * | 2017-12-01 | 2018-05-22 | 上海奕瑞光电子科技股份有限公司 | 一种数字x射线放射系统、自动曝光控制方法及系统 |
CN108294773A (zh) * | 2018-01-03 | 2018-07-20 | 沈阳东软医疗系统有限公司 | 一种曝光图像的获取方法及装置 |
CN110101399A (zh) * | 2018-02-01 | 2019-08-09 | 上海西门子医疗器械有限公司 | 用于x光机的自动曝光剂量调节方法和系统及存储介质 |
WO2019203193A1 (ja) * | 2018-04-17 | 2019-10-24 | 富士フイルム株式会社 | 撮影システム、撮影方法及びプログラム |
CN111281404A (zh) * | 2018-12-06 | 2020-06-16 | 上海西门子医疗器械有限公司 | X线影像设备的成像方法以及x线影像设备 |
CN110269635B (zh) * | 2019-07-09 | 2023-05-02 | 北京友通上昊科技有限公司 | X射线成像设备的自动曝光参数的计算方法及装置 |
WO2021114039A1 (zh) * | 2019-12-09 | 2021-06-17 | 深圳圣诺医疗设备股份有限公司 | 基于蒙版的自动曝光控制方法、装置、存储介质和电子设备 |
CN111568451A (zh) * | 2020-05-22 | 2020-08-25 | 上海联影医疗科技有限公司 | 一种曝光剂量调节方法和系统 |
CN112738391B (zh) * | 2020-12-23 | 2023-04-18 | 上海奕瑞光电子科技股份有限公司 | 自动曝光控制方法及系统 |
CN113014826B (zh) * | 2021-02-18 | 2023-05-05 | 科络克电子科技(上海)有限公司 | 一种图像感光强度参数调整方法、装置、设备及介质 |
CN115379130B (zh) * | 2022-08-25 | 2024-03-29 | 上海联影医疗科技股份有限公司 | 一种自动曝光控制系统、方法、装置和存储介质 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4910592A (en) * | 1988-01-13 | 1990-03-20 | Picker International, Inc. | Radiation imaging automatic gain control |
JP3958069B2 (ja) * | 2001-03-28 | 2007-08-15 | 株式会社東芝 | 放射線測定装置 |
JP5475574B2 (ja) * | 2010-07-02 | 2014-04-16 | 富士フイルム株式会社 | 放射線検出素子、及び放射線画像撮影装置 |
JP5766710B2 (ja) * | 2010-10-26 | 2015-08-19 | 富士フイルム株式会社 | 放射線画像撮影装置およびプログラム |
-
2013
- 2013-11-19 CN CN201310585434.8A patent/CN103705258B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN103705258A (zh) | 2014-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103705258B (zh) | 一种成像设备的自动曝光控制方法及装置 | |
WO2019105490A1 (zh) | 一种数字x射线放射系统、自动曝光控制方法及系统 | |
CN106413236B (zh) | 一种曝光参数调整方法和装置 | |
CN107753048B (zh) | X射线成像装置 | |
US6744846B2 (en) | Method and apparatus for automatic exposure control in CT scanning | |
CN105636327B (zh) | 一种基于透射可见光的曝光控制系统及方法 | |
JP5203946B2 (ja) | 自動露光制御のための方法および装置 | |
CN103699877A (zh) | 提高人脸识别效果的方法和系统 | |
US11166693B2 (en) | Radiation imaging system, control method of radiation imaging system and non-transitory computer-readable storage medium | |
JP7190344B2 (ja) | 画像処理装置、画像処理方法およびプログラム | |
CN112450955A (zh) | 一种ct成像自动剂量调整方法、ct成像方法、系统 | |
CN206453780U (zh) | C型臂x射线机自动曝光控制装置 | |
EP3993706B1 (en) | Narrow beam ct using a 3d fluence modulation and scatter shield system | |
KR20110051490A (ko) | 파노라마 영상 획득 방법 및 장치 | |
CN105662445A (zh) | 数字乳腺x射线成像系统的全自动曝光控制方法和装置 | |
US9820714B2 (en) | Radiographic imaging system, radiographic imaging method and non-transitory recording medium | |
US20140161221A1 (en) | X-ray computed tomography apparatus | |
US6987834B2 (en) | Optimized record technique selection in radiography and fluoroscopy applications | |
US9538977B2 (en) | X-ray diagnosis apparatus and dose distribution data generation method | |
CN106551701B (zh) | 一种无线平板探测器及其图像校正方法 | |
WO2021208374A1 (zh) | 一种自动曝光控制系统及图像校正方法 | |
JP7119015B2 (ja) | 画像データを取得する際の電磁干渉を軽減するためのシステムおよび方法 | |
CN103211607B (zh) | X光机及其控制方法和装置 | |
CN110575193A (zh) | 一种齿科用x射线机自动透视曝光远程控制方法 | |
CN106419944A (zh) | C型臂x射线机自动曝光控制方法、控制系统及控制装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CP03 | Change of name, title or address | ||
CP03 | Change of name, title or address |
Address after: 215000 Jiangsu Suzhou Industrial Park, 218 A2 bio Lake Park, building 501, B3 tower. Patentee after: Jiangsu Kang Zhong digital medical Polytron Technologies Inc Address before: 215123 the 5 floor of B3 5, bio nanotechnology Park, No. 218, Xing Hu Street, Jiangsu Industrial Park. Patentee before: Jiangsu Careray Digital Medical Equipment Co., Ltd. |