CN113411941B - 一种ct灯丝电流校准方法及系统 - Google Patents

一种ct灯丝电流校准方法及系统 Download PDF

Info

Publication number
CN113411941B
CN113411941B CN202110679696.5A CN202110679696A CN113411941B CN 113411941 B CN113411941 B CN 113411941B CN 202110679696 A CN202110679696 A CN 202110679696A CN 113411941 B CN113411941 B CN 113411941B
Authority
CN
China
Prior art keywords
current
exposure
filament
boost
interval
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
Application number
CN202110679696.5A
Other languages
English (en)
Other versions
CN113411941A (zh
Inventor
蒋唯
王斌
陈伟
朱炯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Minfound Medical Systems Co Ltd
Original Assignee
Minfound Medical Systems Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Minfound Medical Systems Co Ltd filed Critical Minfound Medical Systems Co Ltd
Priority to CN202110679696.5A priority Critical patent/CN113411941B/zh
Publication of CN113411941A publication Critical patent/CN113411941A/zh
Application granted granted Critical
Publication of CN113411941B publication Critical patent/CN113411941B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H05ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
    • H05GX-RAY TECHNIQUE
    • H05G1/00X-ray apparatus involving X-ray tubes; Circuits therefor
    • H05G1/08Electrical details
    • H05G1/26Measuring, controlling or protecting

Landscapes

  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Toxicology (AREA)
  • Electron Beam Exposure (AREA)
  • X-Ray Techniques (AREA)

Abstract

本发明提供了一种CT灯丝电流校准方法及系统,涉及医学影像设备领域,包括:获取球管电流区间,对球管电流区间进行分割,获得多个球管电流子区间;对于任一球管电流子区间:从预设数据库中获取初始曝光电流;曝光;判断曝光是否符合预设规则;若是,则生成目标灯丝电流区间;若否,当灯丝电流区间长度为零,预测第二曝光电流并更新初始曝光电流,再次曝光;当灯丝电流区间长度不为零,则获得目标曝光电流以更新初始曝光电流,再次曝光;迭代曝光,直至符合预设规则后生成目标灯丝电流区间;逐个对所有球管电流子区间处理生成目标灯丝电流区间后,结束校准过程;解决现有技术的灯丝电流校准依靠人工主观判断,稳定性较差且校准时间长的问题。

Description

一种CT灯丝电流校准方法及系统
技术领域
本发明涉及医学影像设备领域,尤其涉及一种CT灯丝电流校准方法及系统。
背景技术
电子计算机断层扫描(Computed Tomography,CT),以下简称CT扫描设备,是医学影像领域比较先进的临床检查影像技术。CT扫描设备是利用精确准直的X线束,与灵敏度极高的探测器一同围绕人体的某一部位作一个接一个的断面扫描。根据人体不同组织对X线束的吸收与透过率的不同,再通过探测器对透过人体的X线束进行测量,将获取的数据输入电子计算机,电子计算机对数据进行处理后,就可摄下人体被检查部位的断面或立体的图像。
CT球管中灯丝电流设置过大过小都会使得实际的灯丝电流的初始阶段的电流不稳,与预设值相差过大,在图像上产生伪影,现有校准是采用输入固定灯丝电流值后逐次逼近的方法来获得某个KV、MA条件下最优的灯丝电流值,灯丝电流测得后其值固定,相应的灯丝电流上升基本固定,则相应无效剂量也基本固定,其评价方法比较依靠人工的经验主观判断,而且大量耗费时间。
发明内容
为了克服上述技术缺陷,本发明的目的在于提供一种CT灯丝电流校准方法及系统,解决现有技术的灯丝电流校准依靠人工主观判断,稳定性较差且校准时间长的问题。
本发明公开了一种CT灯丝电流校准方法,应用于CT球管,包括以下:
获取球管电流区间,根据预设间距对所述球管电流区间进行分割,获得多个球管电流子区间;
对于任一所述球管电流子区间:
从预设数据库中获取所述球管电流子区间对应的初始曝光电流;
基于所述初始曝光电流进行曝光,计算灯丝电流曲线和实际电流曲线;
根据所述实际电流曲线判断曝光是否符合预设规则;
若是,则根据所述灯丝电流曲线生成目标灯丝电流区间;
若否,则根据所述灯丝电流曲线确定灯丝电流区间;当所述灯丝电流区间长度为零,则计算所述灯丝电流曲线积分,以预测第二曝光电流并更新所述初始曝光电流,再次曝光;当所述灯丝电流区间长度不为零,则采用二分法迭代获得目标曝光电流以更新所述初始曝光电流,再次曝光;迭代曝光,直至符合预设规则后生成目标灯丝电流区间;
逐个对所有球管电流子区间处理生成目标灯丝电流区间后,结束校准过程。
优选地,根据所述实际电流曲线判断曝光是否符合预设规则,包括以下:
根据所述实际电流曲线获取电流稳定到预设值的时间,判断是否超出预设标准时间;
若超出,则判断曝光不符合预设规则;
若未超出,则判断曝光符合预设规则。
优选地,当所述灯丝电流区间长度为零,则计算所述灯丝电流曲线积分,以预测第二曝光电流,包括以下:
根据下式计算灯丝电流曲线积分:
Figure BDA0003122405480000021
其中,Flagn为第n曝光电流的灯丝电流曲线积分,t0为灯丝电流曲线中毫安值的上升时间;tend为灯丝电流曲线中毫安值从波动到稳定的时间,
Figure BDA0003122405480000022
为灯丝电流的平均值;
根据下式预测第二曝光电流:
Figure BDA0003122405480000023
其中,boostn为第n+1曝光电流,t0为灯丝电流曲线中毫安值的上升时间;tend为灯丝电流曲线中毫安值从波动到稳定的时间,
Figure BDA0003122405480000024
为灯丝电流的平均值。
优选地,在基于所述第二曝光电流进行曝光后,包括以下:
对灯丝电流区间进行更新,当更新后灯丝电流区间长度为零,则再次计算第二曝光电流下的灯丝电流曲线积分,并预测第三曝光电流;
采用第二曝光电流更新初始曝光电流,采用第三曝光电流更新第二曝光电流,迭代曝光,直至灯丝电流区间长度不为零。
优选地,根据下式对灯丝电流区间进行更新:
Figure BDA0003122405480000031
其中,
Flag0=0
[boostmax,boostmin]为灯丝电流区间,boostn为第n曝光电流,Flagn为第n曝光电流的灯丝电流曲线积分。
优选地,所述采用二分法迭代获得目标曝光电流,包括以下:
根据下式获取目标曝光电流:
boostn+1=0.5*(boostmin+boostn);boostmin=boostn+1;Flagn<0
boostn+1=0.5*(boostmax+boostn);boostmax=boostn+1;Flagn>0
其中,boostn为第n曝光电流,为第n曝光电流的灯丝电流曲线积分,[boostmax,boostmin]为灯丝电流区间。
优选地,在结束校准过程后,还包括以下:记录球管曝光次数,当超出预设阈值时,发出校准警示信号。
优选地,所述初始曝光电流对应的灯丝电流区间的上限和下限均与所述初始曝光电流一致。
本发明还提供一种CT灯丝电流校准系统,包括球管,包括以下:
预处理模块,用于获取球管电流区间,根据预设间距对所述球管电流区间进行分割,获得多个球管电流子区间;
获取模块,用于从预设数据库中获取各个球管电流子区间对应的初始曝光电流;
曝光模块,用于基于所述初始曝光电流进行曝光,计算灯丝电流曲线和实际电流曲线;
判断模块,根据所述实际电流曲线判断曝光是否符合预设规则;
存储模块,用于当实际电流曲线判断曝光符合预设规则,则根据所述灯丝电流曲线生成目标灯丝电流区间;
处理模块,用于当实际电流曲线判断曝光不符合预设规则,则根据所述灯丝电流曲线确定灯丝电流区间;当所述灯丝电流区间长度为零,则计算所述灯丝电流曲线积分,以预测第二曝光电流并更新所述初始曝光电流,再次曝光;当所述灯丝电流区间长度不为零,则采用二分法迭代获得目标曝光电流以更新所述初始曝光电流,再次曝光;迭代曝光,直至符合预设规则后生成目标灯丝电流区间;
逐个对球管电流子区间处理生成目标灯丝电流区间后,结束校准过程。
采用了上述技术方案后,与现有技术相比,具有以下有益效果:
本方案中,对球管电流区间进行分区间校准,首先从数据库中获得预设的boostraw作为初始曝光条件,判断该初始曝光条件是否符合预设条件(合理),若否,根据灯丝电流区间长度选择基于灯丝电流曲线积分或二分法获得下一次曝光条件进行迭代曝光,直曝光合理获得灯丝电流区间,利用灯丝曝光实际值的范围(即上述灯丝电流区间)来预设灯丝值,根据实际的曝光值来预设下次的灯丝电流值(计算所述灯丝电流曲线积分,以预测第二曝光电流),实现自动校准,解决现有技术的灯丝电流校准依靠人工主观判断,稳定性较差且校准时间长的问题。
附图说明
图1为本发明所述一种CT灯丝电流校准方法及系统实施例一的流程图;
图2为本发明所述一种CT灯丝电流校准方法及系统实施例一中灯丝电流曲线的实际场景下参考示意图;
图3为本发明所述一种CT灯丝电流校准方法及系统实施例一中实际电流曲线的实际场景下参考示意图;
图4为本发明所述一种CT灯丝电流校准方法及系统实施例一中所述根据所述实际电流曲线判断曝光是否符合预设规则的流程图;
图5为本发明所述一种CT灯丝电流校准方法及系统实施例一中用于体现当曝光不符合预设规则的处理流程图;
图6为本发明所述一种CT灯丝电流校准方法及系统实施例二的模块示意图。
附图标记:
8-灯丝电流校准系统;80-球管;81-预处理模块;82-获取模块;83-曝光模块;84-判断模块;85-存储模块;86-处理模块。
具体实施方式
以下结合附图与具体实施例进一步阐述本发明的优点。
这里将详细地对示例性实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本公开相一致的所有实施方式。相反,它们仅是与如所附权利要求书中所详述的、本公开的一些方面相一致的装置和方法的例子。
在本公开使用的术语是仅仅出于描述特定实施例的目的,而非旨在限制本公开。在本公开和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式,除非上下文清楚地表示其他含义。还应当理解,本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。
应当理解,尽管在本公开可能采用术语第一、第二、第三等来描述各种信息,但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如,在不脱离本公开范围的情况下,第一信息也可以被称为第二信息,类似地,第二信息也可以被称为第一信息。取决于语境,如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”。
在本发明的描述中,需要理解的是,术语“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
在本发明的描述中,除非另有规定和限定,需要说明的是,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是机械连接或电连接,也可以是两个元件内部的连通,可以是直接相连,也可以通过中间媒介间接相连,对于本领域的普通技术人员而言,可以根据具体情况理解上述术语的具体含义。
在后续的描述中,使用用于表示元件的诸如“模块”、“部件”或“单元”的后缀仅为了有利于本发明的说明,其本身并没有特定的意义。因此,“模块”与“部件”可以混合地使用。
实施例一:本实施例提供一种CT灯丝电流校准方法,应用于CT球管,参阅图1,包括以下步骤:
S100:获取球管电流区间,根据预设间距对所述球管电流区间进行分割,获得多个球管电流子区间;
在上述步骤中,球管电流区间为0-maxMa毫安(由高压球馆决定),在本实施方式中,设置每隔50毫安进行校准,即10毫安,50毫安,100毫安到maxMa毫安,依次进行校准。
对于任一所述球管电流子区间,根据以下步骤进行校准:
S200:从预设数据库中获取所述球管电流子区间对应的初始曝光电流;
在上述步骤中,校准开始时会从数据库中获取每个毫安值对应一个可以曝光的初始值boostraw,所述初始曝光电流对应的灯丝电流区间的上限和下限均与所述初始曝光电流一致。灯丝电流区间即为[boostmax,boostmin],如上述初始状态为
boostmin=boostmax=boostraw
S300:基于所述初始曝光电流进行曝光,计算灯丝电流曲线和实际电流曲线;
具体的,灯丝电流曝光时,会产生灯丝电流曲线(filament曲线)(图2)和毫安曲线(即实际电流曲线)(图3),图2和图3的横坐标为曝光时各个角度投影的时间,纵坐标为灯丝电流和毫安值。作为补充的,曝光过程中除曝光电流外的其他曝光条件预先设置,按上述初始曝光电流进行曝光后,将得到球管内部的灯丝电流和实际的曝光的电流曲线,由于在给定的灯丝电流boostraw值后,高压会内部会运行内部的PID校准参数进行内部调整,计算内部灯丝电流的曲线的变化。
S400:根据所述实际电流曲线判断曝光是否符合预设规则;
在上述步骤中,判断曝光是否符合预设规则用于判断此次曝光是否合理。在本实施方式中,采用ma稳定到预设值来判断,具体的,所述根据所述实际电流曲线判断曝光是否符合预设规则,参阅图4,包括以下:
S410:根据所述实际电流曲线获取电流稳定到预设值的时间,判断是否超出预设标准时间;
具体的,通过计算毫安曲线,可以得到ma稳定到预设值的时间tused(即上述步骤中电流稳定到预设值的时间),可根据下式(1)进行判断:
tused<tset(boostraw可行,进行下一球管电流子区间校准)
tused>tset(boostraw不可行,进行下一步骤的校准) (1)
上述式(1)中tset即为预设标准时间,当电流稳定到预设值的时间未超出预设标准时间,则曝光合理,在可继续进行下一球管电流子区间的校准(即下述步骤S500),而电流稳定到预设值的时间超出预设标准时间,则曝光不合理,则需要进行校准(即下述步骤S600)。
S420:若超出,则判断曝光不符合预设规则,即执行步骤S600;
S430:若未超出,则判断曝光符合预设规则,即执行步骤S500。
S500:若是,则根据所述灯丝电流曲线生成目标灯丝电流区间;
根据上述,灯丝电流区间为[boostmax,boostmin],可根据所述灯丝电流曲线获取。
S600:若否,则根据所述灯丝电流曲线确定灯丝电流区间(S600-1)(即判断灯丝电流区间是否为零,即是否具有灯丝电流区间);当所述灯丝电流区间长度为零,则计算所述灯丝电流曲线积分,以预测第二曝光电流(S600-2),并基于所述第二曝光电流进行曝光;当所述灯丝电流区间长度不为零,则采用二分法迭代获得目标曝光电流(S600-3),基于所述目标曝光电流再次曝光;迭代曝光,直至符合预设规则后生成目标灯丝电流区间;
在上述步骤中,参阅图5,确认曝光不符合预设规则,即曝光不合理,则采用计算下一个曝光灯丝电流校准值boostright(不存在灯丝电流区间,即第二曝光电流),或者通过二分法迭代获取最佳的boost值(即boostright)(存在灯丝电流区间,即目标曝光电流),能够获取曝光的boost值和理想的正确值大小关系,将进行第二次曝光,曝光后继续迭代得到一个校准值boostright(即可采用第二曝光电流或目标曝光电流更新初始曝光条件,再次进行曝光,迭代)在这过程最后,也将得到一个灯丝电流的boost区间。
具体的,当所述灯丝电流区间长度为零,则计算所述灯丝电流曲线积分,以预测第二曝光电流,包括以下:
S610:根据下式(2)计算灯丝电流曲线积分:
Figure BDA0003122405480000071
其中,Flagn为第n曝光电流的灯丝电流曲线积分,t0为灯丝电流曲线中毫安值的上升时间;tend为灯丝电流曲线中毫安值从波动到稳定的时间,
Figure BDA0003122405480000072
为灯丝电流的平均值;
S620:根据下式(3)预测第二曝光电流:
Figure BDA0003122405480000073
其中,boostn为第n+1曝光电流,t0为灯丝电流曲线中毫安值的上升时间;tend为灯丝电流曲线中毫安值从波动到稳定的时间,
Figure BDA0003122405480000074
为灯丝电流的平均值。
具体的,上述基于所述第二曝光电流进行曝光后,包括以下:
S621:对灯丝电流区间进行更新,当更新后灯丝电流区间长度为零,则再次计算第二曝光电流下的灯丝电流曲线积分,并预测第三曝光电流;
具体的,根据下式(4)对灯丝电流区间进行更新:
Figure BDA0003122405480000075
其中,
Flag0=0
[boostmax,boostmin]为灯丝电流区间,boostn为第n曝光电流,Flagn为第n曝光电流的灯丝电流曲线积分。
S622:采用第二曝光电流更新初始曝光电流,采用第三曝光电流更新第二曝光电流,迭代曝光,直至灯丝电流区间长度不为零。
在上述步骤S611-S612中,即为迭代过程,在基于第二曝光电流(即上述boostright)曝光后,若曝光还是不合理,则基于第二次曝光产生的灯丝电流区间进行处理,若灯丝电流区间长度为零,则继续预测第三次曝光电流,并再次进行曝光,若灯丝电流区间不为零,则可执行步骤S630以获取目标曝光电流,即最佳的boost值。
更具体的,上述步骤S600中采用二分法迭代获得目标曝光电流,包括以下:
S630:根据下式(5)获取目标曝光电流:
boostn+1=0.5*(boostmin+boostn);boostmin=boostn+1;Flagn<0
boostn+1=0.5*(boostmax+boostn);boostmax=boostn+1;Flagn>0 (5)
其中,boostn为第n曝光电流,为第n曝光电流的灯丝电流曲线积分,[boostmax,boostmin]为灯丝电流区间。
在上述步骤中,获得目标曝光电流后还需要基于该目标曝光电流进行曝光,再次判断曝光是否合理,即循环上述步骤S200,直至判断曝光合理,并生成灯丝电流区间。
S700:逐个对所有球管电流子区间处理生成目标灯丝电流区间后,结束校准过程。
在执行步骤S700时,可以在对每一球管电流子区间执行上述步骤S200-S600之后进行判断是否还存在未执行校准的球管电流子区间,由于需要逐个判断,因此在所有球管电流子区间全部校准后才可结束过程。
在上述步骤S700结束校准过程后,还包括以下:记录球管曝光次数,当超出预设阈值时,发出校准警示信号。
由于考虑到在日常使用球管老化的情况下,boost值会随之变化,因此为了进一步确保上述步骤S200-S700自动校准过程的准确性,球管每曝光一万秒次,系统将会提示操作者进行灯丝电流校准,即重新设置初始曝光电流等。
在本实施方式中,利用灯丝曝光实际值的范围(即上述灯丝电流区间)来预设灯丝值,根据实际的曝光值来预设下次的灯丝电流值(计算所述灯丝电流曲线积分,以预测第二曝光电流),区别与现有技术中通过高压内部反馈的灯丝值,存在实际曝光的电流值和预设的电流值过大的情况,对球管电流区间进行分割,并对分割后的各个子区间进行逐个处理,通过上述步骤S200-S600的迭代,首先从数据库中获得与各个球管电流子区间对应的预设的boostraw作为初始曝光条件,判断该初始曝光条件是否符合预设条件(合理),若否,则判断是否存在曝光区间,即灯丝电流区间,当不存在曝光区间,则根据灯丝电流曲线积分预测下一次曝光条件并判断根据该曝光条件曝光是否符合条件,即步骤S600中S610-S630的迭代,当存在曝光区间,则采用二分法得到boostright,并根据该boostright生成下一次曝光条件并判断根据该曝光条件曝光是否符合条件,直至曝光合理获得灯丝电流区间,实现全条件灯丝电流均能在较短时间里稳定,且在环境变化时能在能自动微调调节灯丝电流,同时还可以降低辐射剂量,保护球管高压。
实施例二:本实施方式提供一种CT灯丝电流校准系统8,参阅图6,包括球管80,还包括以下:
预处理模块81,用于获取球管电流区间,根据预设间距对所述球管电流区间进行分割,获得多个球管电流子区间;
具体的,本方案中对多个球管电流子区间进行逐个校准,并在校准所有球管电流子区间后结束校准。
获取模块82,用于从预设数据库中获取各个球管电流子区间对应的初始曝光电流;
具体的,所述初始曝光电流对应的灯丝电流区间的上限和下限均与所述初始曝光电流一致。灯丝电流区间即为[boostmax,boostmin]。
曝光模块83,用于基于所述初始曝光电流进行曝光,计算灯丝电流曲线和实际电流曲线;
判断模块84,根据所述实际电流曲线判断曝光是否符合预设规则;
具体的,所述预设规则为根据所述实际电流曲线获取电流稳定到预设值的时间,判断是否超出预设标准时间,若超出,则不符合,即不合理,若未超出,则符合,即合理。
存储模块85,用于当实际电流曲线判断曝光符合预设规则,则根据所述灯丝电流曲线生成目标灯丝电流区间;
处理模块86,用于当实际电流曲线判断曝光不符合预设规则,则根据所述灯丝电流曲线确定灯丝电流区间;当所述灯丝电流区间长度为零,则计算所述灯丝电流曲线积分,以预测第二曝光电流,并基于所述第二曝光电流进行曝光;当所述灯丝电流区间长度不为零,则采用二分法迭代获得目标曝光电流,基于所述目标曝光电流再次曝光;迭代曝光,直至符合预设规则后生成目标灯丝电流区间;
具体的,当曝光不符合预设规则,判断是否存在灯丝电流区间,若存在,则采用二分法计算boostright,并根据boostright进行再次曝光,若不存在,则采用灯丝电流曲线积分预测出boostright,根据boostright进行再次曝光,迭代曝光,直至曝光符合预设规则后,得到灯丝电流区间。
逐个对球管电流子区间处理生成目标灯丝电流区间后,结束校准过程。
在上述实施方式中,采用预处理模块进行分割,并对分割后的球管电流子区间进行逐个校准,基于获取模块从数据库中获得与各个球管电流子区间对应的预设的boostraw作为初始曝光条件,并采用曝光模块曝光后判断该初始曝光条件是否符合预设条件(合理),若不符合,则采用处理模块进行迭代处理,即判断是否存在曝光区间,即灯丝电流区间,当不存在曝光区间,则根据灯丝电流曲线积分预测下一次曝光条件当存在曝光区间,则采用二分法得到boostright,并根据该boostright生成下一次曝光条件,再返回曝光模块进行曝光,直至灯丝电流区间,实现利用灯丝曝光实际值的范围来预设灯丝值,根据实际的曝光值来自动预设下次的灯丝电流值(boost),解决现有技术的灯丝电流校准依靠人工主观判断,稳定性较差且校准时间长的问题,实现全条件灯丝电流均能在较短时间里稳定。
应当注意的是,本发明的实施例有较佳的实施性,且并非对本发明作任何形式的限制,任何熟悉该领域的技术人员可能利用上述揭示的技术内容变更或修饰为等同的有效实施例,但凡未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何修改或等同变化及修饰,均仍属于本发明技术方案的范围内。

Claims (5)

1.一种CT灯丝电流校准方法,应用于CT球管,其特征在于,包括以下:
获取球管电流区间,根据预设间距对所述球管电流区间进行分割,获得多个球管电流子区间;
对于任一所述球管电流子区间:
从预设数据库中获取所述球管电流子区间对应的初始曝光电流;
基于所述初始曝光电流进行曝光,计算灯丝电流曲线和实际电流曲线;
根据所述实际电流曲线判断曝光是否符合预设规则;
若是,则根据所述灯丝电流曲线生成目标灯丝电流区间;
若否,则根据所述灯丝电流曲线确定灯丝电流区间;
当所述灯丝电流区间长度为零:
计算所述灯丝电流曲线积分,以预测第二曝光电流并更新所述初始曝光电流,再次曝光;
根据下式计算灯丝电流曲线积分:
Figure FDA0003574988070000011
其中,Flagn为第n曝光电流的灯丝电流曲线积分,t0为灯丝电流曲线中毫安值的上升时间;tend为灯丝电流曲线中毫安值从波动到稳定的时间,
Figure FDA0003574988070000012
为灯丝电流的平均值;
根据下式预测第二曝光电流:
Figure FDA0003574988070000013
其中,boostn+1为第n+1曝光电流,t0为灯丝电流曲线中毫安值的上升时间;tend为灯丝电流曲线中毫安值从波动到稳定的时间,
Figure FDA0003574988070000014
为灯丝电流的平均值;
在基于所述第二曝光电流进行曝光后,包括:对灯丝电流区间进行更新,当更新后灯丝电流区间长度为零,则再次计算第二曝光电流下的灯丝电流曲线积分,并预测第三曝光电流;采用第二曝光电流更新初始曝光电流,采用第三曝光电流更新第二曝光电流,迭代曝光,直至灯丝电流区间长度不为零;
其中,根据下式对灯丝电流区间进行更新:
Figure FDA0003574988070000021
其中,
Flag0=0
[boostmax,boostmin]为灯丝电流区间,boostn为第n曝光电流,Flagn为第n曝光电流的灯丝电流曲线积分;
当所述灯丝电流区间长度不为零:
则采用二分法迭代获得目标曝光电流以更新所述初始曝光电流,再次曝光;迭代曝光,直至符合预设规则后生成目标灯丝电流区间;
所述采用二分法迭代获得目标曝光电流,包括根据下式获取目标曝光电流:
boostn+1=0.5*(boostmin+boostn);boostmin=boostn+1;Flagn<0
boostn+1=0.5*(boostmax+boostn);boostmax=boostn+1;Flagn>0
其中,boostn为第n曝光电流,为第n曝光电流的灯丝电流曲线积分,[boostmax,boostmin]为灯丝电流区间;
逐个对所有球管电流子区间处理生成目标灯丝电流区间后,结束校准过程。
2.根据权利要求1所述的电流校准方法,其特征在于,根据所述实际电流曲线判断曝光是否符合预设规则,包括以下:
根据所述实际电流曲线获取电流稳定到预设值的时间,判断是否超出预设标准时间;
若超出,则判断曝光不符合预设规则;
若未超出,则判断曝光符合预设规则。
3.根据权利要求1所述的电流校准方法,其特征在于,在结束校准过程后,还包括以下:记录球管曝光次数,当超出预设阈值时,发出校准警示信号。
4.根据权利要求1所述的电流校准方法,其特征在于,包括以下:
所述初始曝光电流对应的灯丝电流区间的上限和下限均与所述初始曝光电流一致。
5.一种CT灯丝电流校准系统,包括球管,其特征在于,包括以下:
预处理模块,用于获取球管电流区间,根据预设间距对所述球管电流区间进行分割,获得多个球管电流子区间;
获取模块,用于从预设数据库中获取各个球管电流子区间对应的初始曝光电流;
曝光模块,用于基于所述初始曝光电流进行曝光,计算灯丝电流曲线和实际电流曲线;判断模块,根据所述实际电流曲线判断曝光是否符合预设规则;
存储模块,用于当实际电流曲线判断曝光符合预设规则,则根据所述灯丝电流曲线生成目标灯丝电流区间;
处理模块,用于当实际电流曲线判断曝光不符合预设规则,则根据所述灯丝电流曲线确定灯丝电流区间;
当所述灯丝电流区间长度为零:
计算所述灯丝电流曲线积分,以预测第二曝光电流并更新所述初始曝光电流,再次曝光;根据下式计算灯丝电流曲线积分:
Figure FDA0003574988070000031
其中,Flagn为第n曝光电流的灯丝电流曲线积分,t0为灯丝电流曲线中毫安值的上升时间;tend为灯丝电流曲线中毫安值从波动到稳定的时间,
Figure FDA0003574988070000032
为灯丝电流的平均值;
根据下式预测第二曝光电流:
Figure FDA0003574988070000033
其中,boostn+1为第n+1曝光电流,t0为灯丝电流曲线中毫安值的上升时间;tend为灯丝电流曲线中毫安值从波动到稳定的时间,
Figure FDA0003574988070000034
为灯丝电流的平均值;
在基于所述第二曝光电流进行曝光后,包括:对灯丝电流区间进行更新,当更新后灯丝电流区间长度为零,则再次计算第二曝光电流下的灯丝电流曲线积分,并预测第三曝光电流;采用第二曝光电流更新初始曝光电流,采用第三曝光电流更新第二曝光电流,迭代曝光,直至灯丝电流区间长度不为零;
其中,根据下式对灯丝电流区间进行更新:
Figure FDA0003574988070000035
其中,
Flag0=0
[boostmax,boostmin]为灯丝电流区间,boostn为第n曝光电流,Flagn为第n曝光电流的灯丝电流曲线积分;
当所述灯丝电流区间长度不为零:
则采用二分法迭代获得目标曝光电流以更新所述初始曝光电流,再次曝光;迭代曝光,直至符合预设规则后生成目标灯丝电流区间;
所述采用二分法迭代获得目标曝光电流,包括根据下式获取目标曝光电流:
boostn+1=0.5*(boostmin+boostn);boostmin=boostn+1;Flagn<0
boostn+1=0.5*(boostmax+boostn);boostmax=boostn+1;Flagn>0
其中,boostn为第n曝光电流,为第n曝光电流的灯丝电流曲线积分,[boostmax,boostmin]为灯丝电流区间;
逐个对所有球管电流子区间处理生成目标灯丝电流区间后,结束校准过程。
CN202110679696.5A 2021-06-18 2021-06-18 一种ct灯丝电流校准方法及系统 Active CN113411941B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110679696.5A CN113411941B (zh) 2021-06-18 2021-06-18 一种ct灯丝电流校准方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110679696.5A CN113411941B (zh) 2021-06-18 2021-06-18 一种ct灯丝电流校准方法及系统

Publications (2)

Publication Number Publication Date
CN113411941A CN113411941A (zh) 2021-09-17
CN113411941B true CN113411941B (zh) 2022-05-24

Family

ID=77681498

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110679696.5A Active CN113411941B (zh) 2021-06-18 2021-06-18 一种ct灯丝电流校准方法及系统

Country Status (1)

Country Link
CN (1) CN113411941B (zh)

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2776241B2 (ja) * 1994-03-31 1998-07-16 株式会社島津製作所 X線テレビジョン装置
US9338868B2 (en) * 2012-05-24 2016-05-10 Shimadzu Corporation Radiographic equipment
CN105430858B (zh) * 2015-11-06 2017-06-23 苏州博思得电气有限公司 一种x射线管的灯丝电流值校准方法及装置
CN105769232B (zh) * 2016-02-22 2018-01-12 上海联影医疗科技有限公司 Ct设备的x射线管灯丝预热方法和预热电路
EP3384847B1 (de) * 2017-04-05 2019-06-12 Siemens Healthcare GmbH Bestimmen eines referenzdosisparameters einer computertomographiebildgebung
CN109041392B (zh) * 2018-09-28 2020-01-24 苏州博思得电气有限公司 灯丝校正方法、装置及电子设备
CN211152311U (zh) * 2019-11-19 2020-07-31 明峰医疗系统股份有限公司 一种用于预防ct机探测器温控失控的保险系统
CN111050455B (zh) * 2019-12-04 2022-05-06 深圳市安健科技股份有限公司 球管管电流自动校准方法、装置、电子设备及存储介质
CN111568457B (zh) * 2020-05-25 2023-08-04 明峰医疗系统股份有限公司 一种分段控制灯丝电流快速软启动的调节方法
CN112816929A (zh) * 2021-01-04 2021-05-18 明峰医疗系统股份有限公司 一种ct设备的球管灯丝电流校准系统、方法及ct设备

Also Published As

Publication number Publication date
CN113411941A (zh) 2021-09-17

Similar Documents

Publication Publication Date Title
Hofheinz et al. Effects of cold sphere walls in PET phantom measurements on the volume reproducing threshold
CN100581462C (zh) 用于调整x射线设备的成像参数的装置和方法
US9020220B2 (en) X-ray computed tomography scanner, dose calculation method, and program
US8150110B2 (en) ROI-based rendering for diagnostic image consistency
CN109788926B (zh) 光谱计算机断层扫描(ct)的光谱校准
US6744849B2 (en) Image processing apparatus, image processing method, program, and storage medium
JP2006172457A (ja) 画像配置の最適化
US6148060A (en) Integrated automatic exposure control for portal imaging in radiotherapy
US20190059827A1 (en) System and method for imaging a subject
CN113411941B (zh) 一种ct灯丝电流校准方法及系统
EP1252550B1 (en) Automated calibration adjustment for film dosimetry
JP2003047606A (ja) X線源の放射線量を調節する装置及び方法
US6341172B1 (en) Acquisition scheme for an electron portal imaging system
CN112674785B (zh) 一种应用于ct设备的ma值调节方法、系统及计算机可读存储介质
US20170039735A1 (en) Computed tomography self-calibration without calibration targets
US20220313195A1 (en) Image processing apparatus and storage medium
CN110916703A (zh) 扫描剂量调制方法、装置、扫描设备及存储介质
US7111985B2 (en) Method and system for measuring table sag
CN111699508A (zh) 在预处置和后处置正电子发射断层摄影研究中校正标准化摄取值
US20200015767A1 (en) Sensitivity optimized patient positioning system for dark-field x-ray imaging
CN114972115A (zh) 一种图像降噪方法、装置、设备及介质
CN110678123B (zh) 改进的x射线图像中的几何测量
US6826252B2 (en) Computed tomography scanner and method for controlling the computed tomography scanner
CN113627492B (zh) 扫描对象尺寸的确定方法、扫描方法、装置和电子设备
US11759163B2 (en) Radiation imaging system, control apparatus, control method of radiation imaging system, and storage medium

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
GR01 Patent grant
GR01 Patent grant