CN113507887B - 血管造影成像***的时间校准 - Google Patents
血管造影成像***的时间校准 Download PDFInfo
- Publication number
- CN113507887B CN113507887B CN202080017712.8A CN202080017712A CN113507887B CN 113507887 B CN113507887 B CN 113507887B CN 202080017712 A CN202080017712 A CN 202080017712A CN 113507887 B CN113507887 B CN 113507887B
- Authority
- CN
- China
- Prior art keywords
- watermark
- timestamp
- image
- computer
- radiopaque
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 26
- 239000003550 marker Substances 0.000 claims abstract description 53
- 238000000034 method Methods 0.000 claims description 56
- 238000012545 processing Methods 0.000 claims description 16
- 238000012937 correction Methods 0.000 claims description 14
- 230000001360 synchronised effect Effects 0.000 claims description 13
- 230000005855 radiation Effects 0.000 claims description 10
- 230000000747 cardiac effect Effects 0.000 claims description 9
- 238000002583 angiography Methods 0.000 claims description 7
- 230000001902 propagating effect Effects 0.000 claims description 2
- 238000001514 detection method Methods 0.000 claims 1
- 239000002872 contrast media Substances 0.000 abstract description 5
- 239000000126 substance Substances 0.000 abstract description 4
- 210000004204 blood vessel Anatomy 0.000 abstract description 2
- 230000015654 memory Effects 0.000 description 13
- 238000004891 communication Methods 0.000 description 9
- 239000002184 metal Substances 0.000 description 8
- 238000010801 machine learning Methods 0.000 description 6
- 210000003484 anatomy Anatomy 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000010365 information processing Effects 0.000 description 4
- 230000007246 mechanism Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000002591 computed tomography Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000000877 morphologic effect Effects 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 230000002792 vascular Effects 0.000 description 2
- 239000013598 vector Substances 0.000 description 2
- 240000007651 Rubus glaucus Species 0.000 description 1
- 235000011034 Rubus glaucus Nutrition 0.000 description 1
- 235000009122 Rubus idaeus Nutrition 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 239000011248 coating agent Substances 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000002085 persistent effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B90/00—Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges
- A61B90/39—Markers, e.g. radio-opaque or breast lesions markers
-
- 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/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/504—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of blood vessels, e.g. by angiography
-
- 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/46—Arrangements for interfacing with the operator or the patient
- A61B6/467—Arrangements for interfacing with the operator or the patient characterised by special input means
- A61B6/468—Arrangements for interfacing with the operator or the patient characterised by special input means allowing annotation or message recording
-
- 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/48—Diagnostic techniques
- A61B6/486—Diagnostic techniques involving generating temporal series of image data
-
- 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/48—Diagnostic techniques
- A61B6/486—Diagnostic techniques involving generating temporal series of image data
- A61B6/487—Diagnostic techniques involving generating temporal series of image data involving fluoroscopy
-
- 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
-
- 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/5288—Devices using data or image processing specially adapted for radiation diagnosis involving retrospective matching to a physiological signal
-
- 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/5294—Devices using data or image processing specially adapted for radiation diagnosis involving using additional data, e.g. patient information, image labeling, acquisition parameters
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T1/00—General purpose image data processing
- G06T1/0021—Image watermarking
- G06T1/0085—Time domain based watermarking, e.g. watermarks spread over several images
-
- 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/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
- G06T7/0014—Biomedical image inspection using an image reference approach
- G06T7/0016—Biomedical image inspection using an image reference approach involving temporal comparison
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/22—Image preprocessing by selection of a specific region containing or referencing a pattern; Locating or processing of specific regions to guide the detection or recognition
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/20—ICT specially adapted for the handling or processing of medical images for handling medical images, e.g. DICOM, HL7 or PACS
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B90/00—Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges
- A61B90/39—Markers, e.g. radio-opaque or breast lesions markers
- A61B2090/3966—Radiopaque markers visible in an X-ray image
-
- 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/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/404—Angiography
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V2201/00—Indexing scheme relating to image or video recognition or understanding
- G06V2201/02—Recognising information on displays, dials, clocks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V2201/00—Indexing scheme relating to image or video recognition or understanding
- G06V2201/03—Recognition of patterns in medical or anatomical images
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Radiology & Medical Imaging (AREA)
- Surgery (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Pathology (AREA)
- Optics & Photonics (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Multimedia (AREA)
- Vascular Medicine (AREA)
- Dentistry (AREA)
- Human Computer Interaction (AREA)
- Quality & Reliability (AREA)
- Physiology (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
通过向血管内注射化学造影剂并对造影剂因变于时间的通路进行成像来获得血管造影数据,从而生成图像的序列。为了校正嵌在图像元数据中的未校准时间戳的误差,使用不透射线标志物在获得的图像中生成嵌入时间戳数据的水印。不透射线标志物以动态水印的形式在x射线图像上造成不透明,该动态水印对时间戳进行编码。然后,对标志物在水印(从不透射线标志物投射)中的位置进行处理和分析,以生成图像的精确时间戳。通过生成精确时间戳,提供了该图像与其他数据源同步的计算。
Description
相关申请的交叉引用
本申请要求于2019年3月8日提交的美国临时申请第62/815,476号的优先权,该临时申请的内容通过引用整体并入本文。
技术领域
该领域涉及生成时间戳的血管造影或其他x射线成像***,并且特别地,涉及生成具有亚秒精度的精确时间戳以允许与外部信号同步的血管造影或其他x射线成像***。
背景技术
血管造影采集装备的制造商没有预见到需要在血管造影图像上包括精确时间戳。在血管造影图像中作为元数据存储的时间戳示出了从几秒到几分钟的误差范围。当前的血管造影成像***不包括用于重新校准***时间戳的机构,诸如应用编程接口。
血管造影图像的通用文件格式使用DICOM标准,其中一系列血管造影图像中的每个图像包括具有DICOM元数据的报头。虽然DICOM元数据包括用于存储时间戳信息的标签,但存储的时间戳已经被确定为不正确,并且不指示采集相对应的血管造影图像的正确时间。
附加地,现有血管造影装备不包括用于与外部信号(诸如心脏信号)同步的硬件输入端口。此外,此类装备不具有用于统一存储外部心脏信号和血管造影图像的内部数据库或文件结构。
虽然使用计算机断层扫描和MRI获得的心脏门控图像提供了计算机断层扫描或MRI数据与外部信号的定时的事后关联,但该装备已经被制造为包括允许计算机断层扫描或MRI硬件的内部时钟与外部信号同步的设计规范。然而,对于缺乏用于内部时钟同步的硬件或软件的血管造影***,这些类型的***不会起作用。
因此,在缺乏内部时钟同步的血管造影和其他x射线成像***中,需要用于获得精确时间戳的技术。
发明内容
本技术涉及用于在x射线成像***(诸如具有不同步到时间标准(例如,诸如全球或国际时间标准、美国国家标准与技术研究所(NIST)时间标准,或任何其他合适的时间标准)的内部时钟的血管造影成像***))中生成具有亚秒精度的精确时间戳的***、方法和计算机可读介质,从而允许将成像***获得的图像同步到外部信号。
多个可移动不透射线标志物被定位在辐射发射器和辐射检测器阵列之间。使得所述多个不透射线标志物中的每个不透射线标志物以各自的频率移动。当所述不透射线标志物正在移动时,获得一系列x射线图像,每个图像包括由所述多个移动不透射线标志物投射的动态生成的水印,其中每个不透射线标志物的位置被示出在所述水印上。一系列图像被处理以测量水印的多个标记中的每个标记的旋转位置。基于所述标记的所测量的旋转位置和所述多个不透射线标志物中的每个不透射线标志物的所述各自的频率,确定时间戳。
在各方面中,水印对具有亚秒精度的时间戳进行编码。在各方面中,多个不透射线标志物中的每个不透射线标志物以由同步到时间标准的控制器确定的速率旋转。
在各方面中,基于由水印编码的时间戳校正元数据时间戳。
在各方面中,多个步进电机中的每个步进电机以不同的指定频率旋转,以便使得相关联的不透射线标志物以相同的频率旋转。在各方面中,至少两个不透射线标志物用于获得亚秒精度。
在各方面中,在用于根据以比心率更快的速率采集的血管造影图像对心脏频率现象时空重建的技术中,如本文所提供的生成的时间戳可以用于使图像与外部信号同步(参见,于2018年11月13日公告的美国专利第10,123,761号,其内容通过引用整体并入本文)。
应当理解,该发明内容并不旨在标识本公开的实施例的关键或必要特征,也不旨在用于限制本公开的范围。通过下文的描述,本公开的其他特征将变得容易理解。
附图说明
通常,各个附图中的相似附图标记用于指定相似部件。
图1示出了根据本文所提供的技术的用于基于x射线图像的动态水印生成精确时间戳的血管造影***。
图2是根据本文所提供的技术的用于使用图像处理技术动态地处理带水印的x射线图像的流程图。
图3是根据本文所提供的技术的用于使用结合图像处理技术的机器学习动态地处理带水印的x射线图像的流程图。
图4是根据本文所提供的技术的用于根据基于经动态处理的带水印的x射线图像校正元数据时间戳来生成时间戳的流程图。
图5是根据本文所提供的技术的用于依据经动态处理的带水印的x射线图像来生成时间戳的流程图。
图6是根据本文所提供的技术的用于基于x射线图像的动态水印生成具有亚秒精度的时间戳的操作的高级流程图。
图7是根据本文所提供的技术的可以与血管造影***一起使用以存储和分析包括水印的血管造影图像的计算机***或信息处理设备的框图。
图8A和图8B示出了根据本文所提供的技术的不同类型的刻度盘的示例。
具体实施方式
血管造影成像***包括硬件和软件,该硬件包括x射线发射器、x射线检测器阵列和电子硬件,该软件将由x射线检测器阵列捕获的数据转化成图像。为了获得血管造影图,向患者血管内注射化学造影剂的团注,并获得x射线的序列或时间系列。血管***的解剖结构的二维投影被捕获为在x射线投影路径中通过血管***的阻止x射线通路的化学造影剂。
化学造影剂的团注通路定时可以提供医学信息。例如,在团注通路中相对较早采集的图像更可能反映动脉解剖结构,而相对较晚采集的图像更可能反映静脉解剖结构。根据采集的时间排序的这些图像的聚合包括血管造影图。
每个图像都具有时间戳,时间戳与由x射线信号检测器阵列产生的信号被聚合成数字化图像的时刻相对应。该时间戳可以存储在元数据中。x射线和其他医学图像的通用图像文件格式使用DICOM标准。该标准提供了被包括作为元数据的时间戳标签,以指示x射线图像的采集时刻。
然而,血管造影成像***是典型的闭合***,其不允许内部时钟的校准。因此,从图像元数据中提取的x射线图像时间戳对于其他应用(诸如为了计算目的将图像与其他时间序列数据同步)通常是不令人满意的。
在一些情况下,对时间戳误差的测量已显示误差高达20分钟或更长。附加地,DICOM标准没有提供足够的时间分辨率(例如,亚秒分辨率),以允许血管造影图像与外部信号(例如,心脏信号)精确同步。
本发明技术提供基于同步到NIST时间标准的用于形成水印的多个不透射线标志物来生成精确时间戳,该水印由所述多个不透射线标志物投射。在一些方面中,基于由不透射线标志物提供的信息,DICOM时间戳经历校正过程。在其他方面中,在不利用DICOM时间戳的情况下,仅从由不透射线标志物提供的信息产生时间戳。
利用精确时间戳,血管造影数据可以与外部信号(诸如由手指脉搏血氧仪或由心电图产生的心脏信号)同步。在实施例中,心脏信号可以与小波算法相结合(在一些情况下为互相关的),以从血管造影图像序列产生单个移动血管脉搏波的时空重建,例如,如上文通过引用并入本文的美国专利第10,123,761号中所公开的。附加地,利用小波算法的外部信号增加血管造影图像序列的信噪比。
如下文所描述的,提供了用于将动态生成的水印***血管造影图的像素中的装置、方法和计算机可读介质,该水印对可以被转换为时间戳的信息进行编码(例如,通过自动化软件)。
参考图1,示出了根据本文所提供的技术的用于基于x射线图像的动态水印生成精确时间戳的***/设备。在该示例中,x射线源3和x射线检测器阵列10是血管造影成像***的一部分,该血管造影成像***不提供直接校准其内部时钟或以其他方式生成具有精确时间戳的图像的能力。
***包括计算机(在图7中以计算***80示出)、控制器1和电机11(例如,步进电机等)、不透射线刻度盘2、x射线源3和x射线检测器阵列10。在该示例中,x射线检测器10示出了来自x射线投影到不透射线刻度盘2上的嵌入时间水印4。在该示例中,不透射线标志物是刻度盘的形式。
在操作中,该***生成图像的校准系列或序列5,图像的校准系列或序列5可以经过处理以将图像内的水印6转换成用于图像的时间戳校正7。时间戳校正7可以包括用于校正元数据的误差信号,或用于替换元数据的时间戳。下文进一步描述这些部件。
计算机可以包括任何合适的计算***,计算***包括独立计算机、服务器、平板计算机、膝上型计算机等。例如,计算机可以包括不在x射线源3和用于生成一系列x射线图像的x射线检测器阵列10之间的视线中的小型便携式计算机。
计算机可被配置为使电机旋转,该电机控制不透射线标志物以指定频率的旋转或移动。例如,便携式计算机可以附连到使用python库编程的控制器以用于控制电机。python库被编程为以预定编程速率转动每个电机。
在如图1所示的不透射线参考标志物***中,提供了多个可移动不透射线标志物(例如,刻度盘2),该多个可移动不透射线标志物在x射线图像上的相对位置对x射线图像采集时间戳进行编码。在各方面中,不透射线参考标志物***包括能够控制不透射线标志物的旋转的任何***。通过示例的方式,可移动不透射线标志物可以是由控制器同步到时间标准的刻度盘的形式。不透射线标志物的移动可以参考任何合适的坐标系(例如,角度坐标系、x-y坐标系等)。不透射线标志物相对于不透射线参考系的相对位置被投射为x射线图像上的动态水印。例如,不透射线标志物可以被投射为图像中的相对应的标记(作为水印的一部分)。
在如图8A所示的一个实施例中,可移动不透射线标志物可以是绕圆形框架810的中心旋转的指针(hand)820的形式。在该方面中,刻度盘的框架是固定的,并且指针的旋转由控制器进行控制。可以以这种方式实现多个刻度盘,以达到适当的时间分辨率。
在如图8B所示的另一个实施例中,可移动不透射线标志物可以包括刻度盘本身的盘面,该盘面可以被配置为绕该盘面的中心或某个其他点旋转。在该实施例中,不透射线标志物840被定位在旋转盘面830上(例如,以刻度线的形式或其他合适的形式等)。
在每种情况下,一个或多个电机可以用于(例如,通过旋转可移动指针或可移动盘面)旋转刻度盘的可移动不透射线标志物。在示例中,刻度盘形式的不透射线标志物的旋转由单独的电机驱动。在其他示例中,单个电机可以经由产生不同旋转频率的齿轮耦合到多个刻度盘。在一方面,电机可以直接连接到刻度盘。在另一个方面,电机可以经由一个或多个连接部件连接到刻度盘。
相应地,尽管关于刻度盘的旋转描述了包括如通篇所提供的各种形式的设备,但将理解,定时由不透射线标志物进行编码。因此,本实施例涵盖(例如,由控制器控制的)可移动不透射线标志物对时间信息进行编码的任何合适的形式。不透射线标志物可以包括刻度盘的形式。
由控制器控制的电机可以控制刻度盘2以咔哒声或平滑扫掠或任何其他合适的方式移动,以在获得x射线图像的序列的同时以指定频率旋转。刻度盘被放置于x射线源3和x射线检测器阵列10之间的视线内。在各方面中,刻度盘可以被放置在金属安装支架中,并且具有附连到可旋转轴的金属指针。金属指针在x射线图像中是不透射线的,并且显示为例如,标记。在其他方面中,刻度盘的盘面可以旋转,并且不透射线标志物可以以一定图案布置在刻度盘的盘面上,其中该图案允许基于盘面的旋转对定时信息进行编码。在其他方面中,不透射线标志物可以由其他类型的不透射线材料制成,或由覆盖或涂覆有不透射线涂料的非不透射线材料制成。与从x射线管3传播到x射线检测器阵列10上的x射线相对应的像素被不透明化,从而在图像的检测器阵列10中形成水印。本文所提供的技术涵盖将刻度盘(例如,金属标志物或不透射线标志物)直接附连到电机,以及通过一个或多个连接设备将刻度盘连接到电机。
利用电机及其各自以所编程速率旋转的刻度盘,获得了血管造影图像的序列。在一些情况下,在与血管造影图像系列通常使用的相同的历时和图像帧速率的情况下获得图像。刻度盘上的不透射线标志物对x射线不透明,并且当在刻度盘在x射线源路径中的情况下收集图像时,在图像上形成包括与不透射线标志物的位置相对应的标记的水印。
因此,水印由诸如一个或多个刻度盘/(多个)可旋转不透射线标志物之类的设备产生,每个不透射线标志物以指定频率旋转并且直接或间接地同步到控制器。在一些方面中,控制器可以被同步到NIST时间戳标准(例如,通过NIST同步控制器)。然而,本技术可以用于同步到任何合适的时间标准。来自不透射线标志物的水印被嵌入到血管造影图的像素中,并且可以被转换成可以与图像元数据中的时间戳进行比较的形式。这可以用于确定时间戳校正,时间戳校正可以被应用于图像序列的元数据(例如,DICOM时间戳),或者可以独立于图像元数据被转换为时间戳。
可以使用任何合适数量的不透射线标志物或刻度盘,其中多个刻度盘中的每个刻度盘以任何合适的频率旋转不透射线标志物,以产生对时间信息进行编码的水印。特别地,可以将刻度盘的数量和不透射线标志物的旋转的频率配置为以指定的分辨率对时间戳进行编码。可以在指定频率(振荡的时间周期)下采用更多数量的电机和刻度盘,以提供小时、分钟、秒和亚秒分辨率。可以在指定频率下采用至少两个电机,以提供分钟和秒分辨率。
通过非限制性示例的方式,不在x射线源和x射线检测器阵列之间的视线内的小型便携式计算机(例如,Raspberry Pi 3(树莓派3),https://www.raspberrypi.org,Linux操作***)被附连到驱动数字步进电机的硬件单元(例如,Pi板电机板-R1.0(PiPlatesMotorplate-R1.0)https://pi-plates.com/motorr1/)。一组两个步进电机(例如,Polulu零件号SY20SSH30-0604A https://www.pololu.com/product/1204)被放置在x射线源和x射线检测器阵列之间的视线内。步进电机被放置在金属安装支架中,并且具有附连到轴的充当定时指针的金属指针。
步进电机控制单元使用python库控制步进电机(https://pi-plates.com/code/#Object_Oriented_MOTORplate_Demo_Using_Tkinter),并且被编程为以所编程速率转动步进电机刻度盘。步进电机刻度盘中的一个刻度盘被编程为每10秒完成一个循环,另一个刻度盘被编程定为每秒完成一个循环。金属指针对x射线不透明,并在x射线检测器阵列上投射阴影,从而允许在图像上留下两个刻度盘指针的图像印记。
在一些方面中,当人类未被血管造影成像时,获得血管造影图像(包括水印)。在该情况下,时间戳校正被存储在计算机存储器或计算机数据库中,以用于随后校正人类临床血管造影图的时间戳。可以在对人类患者进行血管造影研究之前和/或之后获得一组校准时间戳,以验证时间戳校正在一天或一周或更长时间内尚没有显著的漂移。校准可以在人类临床血管造影之前和/或之后执行。
在其他方面中,在为人类获得血管造影图像的同时获得血管造影图像系列,但是刻度盘的位置被放置成使得所得水印不会模糊与血管造影相关的信息。因此,在人类临床血管造影期间,可以通过将刻度盘装置放置在被研究的相关人类解剖结构的x射线投影的线之外来校准图像。
因此,这种方法独立于血管造影成像***而运行。将校准x射线图像的序列5导出以供进一步的分析,诸如图2或图3中所描述的用于基于被投射以形成水印(例如,包括与不透射线标志物相对应的一系列标记)的不透射线标志物的位置(例如,通过图像处理技术或通过机器学习***)来生成时间戳的过程(包括计算机中运行的软件程序)。
图2和图3中提供的示例参考刻度盘,该刻度盘包括具有旋转刻度盘指针的非可移动框架。在该示例中,不透射线标志物是刻度盘的形式。然而,这些技术可扩展到其他刻度盘形式,诸如本文所提供的刻度盘的旋转盘面,或者其中不透射线标志物由控制器控制并且定时通过由多个移动不透射线标志物投射到x射线图像上的水印进行编码的任何其他形式。
图2是根据本文所提供的技术的用于使用图像处理技术动态地处理带水印的x射线图像以获得时间戳的流程图。可以标识水印,并且可以基于刻度盘指针投影的位置和刻度盘指针的旋转的频率对经编码的时间戳进行解码。在该示例中,不透射线标志物包括刻度盘指针。不透射线刻度盘指针被投射到x射线图像上,作为标记/投影,在这种情况下,该标记/投影具有刻度盘指针的形状。
在操作210处,获得一系列的x射线图像,每个图像包括通过一个或多个不透射线刻度盘指针的投射而生成的水印。水印对与时间戳相对应的信息进行编码。在该示例中,可以将刻度盘定位成使得水印被生成在图像的角落附近,或者生成在不包括相关医疗信息的区域中。如该示例中所示,通过金属的合适的放置,刻度盘的边框也可以任选地位于水印中。
在操作220处,针对每个图像,获得图像内水印的位置,包括每个刻度盘投影的中心。例如,图像分析软件可以被配置为标识圆形形状,并且标识每个所标识的圆的中心。在其他方面中,水印(一条或多条线条)的位置可以被提供给图像分析软件。参考标记也可以被包括以指示图像的适当取向(例如,当参考标记被定位在图像的右上角时,该图像可以被分析为指示整个图像未经旋转或转置)。在一些方面中,刻度盘投影可以以预定尺寸出现在图像中,并且图像分析***可以被配置为标识预定半径的圆。
在操作230处,针对每个刻度盘,执行形态分析以提取与每个刻度盘指针投影相对应的线条。形态分析可以用于标识与所投影的刻度盘指针相对应的线条。在一些方面中,所标识的线条可以作为向量存储在存储器中。任何合适的功能均可用于提取。
在操作240处,基于该提取确定每个刻度盘指针投影/标记的旋转位置或取向。例如,刻度盘指针投影的旋转位置可以通过刻度盘指针投影的当前位置和参考物的已知与360°圆的旋转的起始(0°)相对应的位置之间的角度来确定。在一些方面中,向量的取向可以按弧度计算。
在操作250处,基于刻度盘指针的旋转的频率和每个刻度盘指针投影的旋转位置来确定时间戳。例如,旋转位置(按弧度)可以基于已知的旋转频率和转换因子转换为时间戳,如下所示。
例如,对于以1/10Hz(每10秒1圈)旋转的第一刻度盘指针,可以基于以下方式确定以秒为单位的定时:10.*刻度盘1位置1/2.Pi
对于以1Hz(每1秒1圈)旋转的第二刻度盘指针,可以基于以下方式确定以秒为单位的定时:
1.*刻度盘2位置2/2.Pi
以秒为单位的最终时间可以通过组合这两个时间来确定:10.*刻度盘1位置1/2.Pi+1.*刻度盘2位置2/2.Pi
在其他方面中,可以基于从查找表(例如,将每个刻度盘指针的旋转位置与时间戳值相关联的查找表)获得的值来确定时间戳,其中查找表是特定于刻度盘指针的频率。
通常,可以使用任意数量的刻度盘来实现期望的时间分辨率。通过示例的方式,图2和图3涉及具有两个刻度盘的***。然而,在其他情况下,可以使用三个、四个、五个、六个、七个、八个或更多个刻度盘。每个刻度盘可以被配置为以预定的旋转频率旋转。通过示例的方式,第一刻度盘可以被配置为以每秒10转来旋转,第二刻度盘可以被配置为以每秒1转来旋转,第三刻度盘可以被配置为以每10秒1转来旋转。此类配置可以允许以高精度获得分钟和亚分钟尺度的精确度。
图3是根据本文所提供的技术的用于使用结合图像处理技术的机器学习来处理x射线图像的动态时间戳水印的流程图。机器学习方法可以用于自动化对图像内的水印的标识(对时间戳进行编码)。一旦被标识,经编码的时间戳可以基于刻度盘投影的位置和刻度盘的旋转的频率来进行解码。
在操作310处,获得一系列x射线图像,每个图像包括水印,水印包括多个刻度盘投影。在该示例中,可以将刻度盘定位成使得水印被生成在图像的角落附近,或者生成在不包括相关医疗信息的区域中。参考标记也可以被包括以指示图像的适当取向(例如,当参考标记被定位在图像的右上角时,该图像可以被分析为指示整个图像未经旋转或转置)。在操作320处,使用机器学习(ML)***在每个图像内标识水印的位置。在操作330处,使用霍夫变换或用于检测线条的其他合适的变换来检测每个刻度盘指针投影。在操作340处,确定每个刻度盘指针投影(例如,来自霍夫变换的线条)的旋转位置。例如,刻度盘指针投影的旋转位置可以通过刻度盘指针投影的当前位置和参考物的已知与360°圆的旋转的起始(0°)相对应的位置之间的角度来确定。在操作350处,基于刻度盘投影的旋转位置和每个刻度盘指针的旋转的频率来确定时间戳。在其他方面中,可以基于从查找表获得的值来确定时间戳,该查找表例如将每个刻度盘指针投影的旋转位置与时间戳值相关联的查找表,其中查找表是特定于刻度盘的频率。
图4是根据本文所提供的技术的用于通过基于带水印的图像校正元数据时间戳来生成时间戳的流程图。
可以从根据图1的血管造影研究中获得一系列图像,其中每个图像包括动态生成的水印,并且与包括与元数据时间戳(包括至少小时和分钟)相对应的标签的报头相关联。在该示例中,可以DICOM格式提供的时间戳元数据410通过现有的血管造影成像***嵌入到图像元数据中。DICOM格式包括指示每个图像的采集时刻的元数据标签(时间戳)。然而,DICOM时间戳不正确。
通过确定误差、并以该误差校正每个DICOM时间戳图像410的时间戳来生成经校正的DICOM时间戳图像420。例如,从一系列图像的图像中提取DICOM元数据时间戳。将动态水印时间戳和DICOM元数据时间戳进行比较,并且动态水印时间戳和DICOM元数据时间戳的算术平均值差(mean difference)或平均值差(average difference)(跨图像集合)被用于确定针对血管造影研究的时间戳校正。DICOM时间戳可以通过误差进行校正,以生成具有第二分辨率的经校正的DICOM时间戳图像420。
然而,DICOM时间戳具有第二分辨率,并且需要亚秒分辨率。为了更新时间戳以包括亚秒分辨率,可以分析水印以获得亚秒定时信息。
因此,通过基于投射到水印上的旋转刻度盘/不透射线标志物的位置将经校正的DICOM时间戳420与亚秒定时相结合来生成经校正的时间戳430。时间戳被更新以包括基于DICOM时间戳未提供的亚秒分辨率。
该校正可以被存储以在随后利用人类临床血管造影图校正DICOM元数据时间戳时使用。可以在人类临床血管造影之前和之后确定误差,以说明DICOM元数据因变于时间(例如,在几天、几周或甚至几个月内)的漂移。在一些方面中,经校正的时间戳430可以替换DICOM时间戳,而在其他方面中,经校正的时间戳430可以与DICOM时间戳分开存储在不同的标签中。
时间戳校正信息(例如,平均误差)可以以可计算的可读格式保存,以允许随后与其他相关数据源(诸如包括可以由与血管造影***并行的单独硬件***同时获得的心电图的生理信号)进行同步的计算。这些方法可以用于使来自没有实时连接到血管造影采集装备的***的数据与由血管造影装备生成的数据同步。因此,如美国专利第10,123,761号(“761专利”)所教导的,该方法提供外部信号(例如,心脏信号)与血管造影图像的事后集成,以准备应用提供对心脏频率活动的时空描述的重建算法。
图5是根据本文所提供的技术的用于依据经动态水印标记的x射线图像来生成时间戳的流程图。在该示例中,虽然DICOM时间戳可能存在于元数据中,但它未被利用。相反,仅使用来自多个时间同步的不透射线标志物/刻度盘的水印来生成时间戳。合适数量的不透射线标志物/刻度盘被编程为以预定义频率旋转,以获得期望的时间戳范围和分辨率(例如,按照小时、分钟、秒、亚秒等)。
获得图像510,其中水印对来自多个不透射线标志物/刻度盘的定时信息进行编码,每个不透射线标志物/刻度盘以不同频率旋转。对水印进行处理(例如,根据图2或图3),使得每个不透射线标志物/刻度盘的旋转位置被确定。基于每个不透射线标志物/刻度盘的旋转位置和旋转的频率,确定每个不透射线标志物/刻度盘的时间,并且将各个时间组合成具有亚秒精度的时间戳,以从时间同步的刻度盘生成时间戳520。
图6是根据本文所提供的技术的用于基于x射线图像的动态时间戳水印生成具有亚秒精度的时间戳的操作的高级流程图。
在操作610处,将多个可移动不透射线标志物定位在辐射发射器和辐射检测器阵列之间。使得多个不透射线标志物中的每个标志物以各自的频率移动(例如,通过控制器)。
在操作620处,当不透射线标志物正在移动时,获得一系列x射线图像,每个图像包括由多个移动不透射线标志物投射的动态生成的水印,其中每个不透射线标志物的位置被示出在水印上(例如,被示出为标记)。水印对时间戳进行编码,时间戳可以基于频率和位置进行解码。
在操作630处,该系列被处理以测量投射到图像上的多个标记中的每个标记的旋转位置。在操作640处,基于标记的所测量的旋转位置和多个不透射线标志物中的每个标志物的各自的频率,确定具有亚秒精度的时间戳。
本文所提供的技术的优点包括生成精确时间戳,这允许外部信号与血管造影图像同步。因此,通过生成正确的时间戳,外部心脏信号或其他信号可以与获得的图像同步,以增加获得的血管造影图像的信噪比(例如,结合互相关小波算法(cross-correlatedwavelet algorithm),该算法依据血管造影图像序列产生单个移动血管脉搏波的时空重建,如‘761专利中所公开的)。
该同步可以作为事后关联过程的一部分来执行。根据需要,可以使用插值来微调图像的时间戳与外部信号的同步。
图7是可以与本发明的实施例一起使用的计算机***或信息处理设备的框图。计算机***或信息处理设备80被示出为可以与图1-图6和图8的***一起使用以获得一系列带水印的血管造影图像。一旦获得,***可以根据本文所提供的技术中的任何一种或多种技术来处理血管造影图像以生成精确时间戳。
图7是根据本公开中的技术所编程的通用计算机***80或用于本文所提供的实施例的专用信息处理设备的图示,并且并不旨在限制本文所公开的主题的范围。本领域的普通技术人员将认识到仍在本公开的范围和等效物内的计算机***80的其他变型、修改和替代方案。
在一个实施例中,计算机***80包括监视器82、计算机84(其包括(多个)处理器86、总线子***88、存储器子***90和盘子***92)、用户输出设备94、用户输入设备96和通信接口98。监视器82可以包括被配置为生成信息的可视表示或显示的硬件和/或软件元件。监视器82的一些示例可以包括熟悉的显示设备,诸如电视监视器、阴极射线管(CRT)、液晶显示器(LCD)等。在一些实施例中,监视器82可以提供输入接口,诸如结合触摸屏技术。
计算机84可以包括熟悉的计算机部件,诸如一个或多个中央处理单元(CPU)、存储器或存储设备、图形处理单元(GPU)、通信***、接口卡等。如图7所示,计算机84可以包括一个或多个处理器86,该一个或多个处理器86经由总线子***88与多个***设备通信。(多个)处理器86可以包括市售的中央处理单元等。总线子***88可以包括用于使计算机84的各种部件和子***彼此按照预期通信的机构。尽管总线子***88示意性地示出为单个总线,但总线子***的替代实施例可以利用多个总线子***。与(多个)处理器86通信的***设备可以包括存储器子***90、盘子***92、用户输出设备94、用户输入设备96、通信接口98等。
存储器子***90和盘子***92是被配置为存储数据的物理存储介质的示例。存储器子***90可以包括多个存储器,该多个存储器包括用于在程序执行期间对程序代码、指令和数据进行易失性存储的随机存取存储器(RAM),以及存储固定程序代码、指令和数据的只读存储器(ROM)。盘子***92可以包括为程序和数据提供永久(非易失性)存储的多个文件存储***。其他类型的物理存储介质包括软盘、移动硬盘、光学存储介质(诸如CD-ROM、DVD和条形码)、半导体存储器,诸如闪存、只读存储器(ROM)、电池支持的易失性存储器、网络存储设备等。
存储器子***90和盘子***92可以被配置为存储提供本文所讨论的技术的功能或特征的编程和数据构造。当由(多个)处理器86执行时实现或以其他方式提供功能的软件代码模块和/或处理器指令可以存储在存储器子***90和盘子***92中。
用户输入设备94可以包括被配置为接收来自用户的输入以供计算机***80的部件进行处理的硬件和/或软件元件。用户输入设备可以包括用于向计算机***84输入信息的所有可能类型的设备和机构。这些可以包括键盘、小键盘、触摸屏、并入显示器的触摸界面、音频输入设备(诸如麦克风和语音识别***)以及其他类型的输入设备。在各种实施例中,用户输入设备94可以被实现为计算机鼠标、轨迹球、跟踪板、操纵杆、无线遥控器、绘图板、语音命令***、眼睛跟踪***等。在一些实施例中,用户输入设备94被配置为允许用户经由命令、动作或手势(诸如点击按钮等)选择监视器82上可能出现的对象、图标、文本等或以其他方式与监视器82上可能出现的对象、图标、文本等交互。
用户输出设备96可以包括被配置为从计算机***80的部件向用户输出信息的硬件和/或软件元件。用户输出设备可以包括用于从计算机84输出信息的所有可能类型的设备和机构。这些可以包括显示器(例如,监视器82)、打印机、触摸或力反馈设备、音频输出设备等。
通信接口98可以包括被配置为提供与其他设备的单向或双向通信的硬件和/或软件元件。例如,通信接口98可以(诸如经由互联网连接)提供计算机84与其他通信网络和设备之间的接口。
图7表示能够体现本发明的实施例的计算机***。对于本领域的普通技术人员显而易见的是,许多其他硬件和软件配置适合与本发明一起使用。例如,计算机可以是台式、便携式、机架安装或平板计算机配置。附加地,计算机可以是一系列联网计算机。在其他实施例中,上述技术可以在芯片或辅助处理板上实现。
由此可见,在从前面的描述中显而易见的那些对象中,上文阐述的对象是有效获得的并且因为在不脱离本发明的精神和范围的情况下,可以在实施上述方法和(多个)构造中进行某些改变,旨在将上述描述中包含的以及附图中示出的所有内容解释为说明性的,而不是限制性的。
还应当理解,以下权利要求旨在涵盖本文所描述的本发明的所有通用和特定特征并且本发明的范围的在语言上的所有陈述可以说介于两者之间。
已经根据优选实施例描述了这些实施例,并且认识到,除了明确说明的那些以外,等效例、替代例和修改例是可能的,并且这些等效例、替代例和修改例在附加权利要求的范围内。
Claims (28)
1.一种对x射线图像进行动态时间戳水印标记的方法,包括:
将多个可移动不透射线标志物定位在成像设备的辐射发射器和辐射检测器阵列之间;
使得所述多个不透射线标志物中的每个不透射线标志物以各自的频率移动;
当所述不透射线标志物正在移动时,获得一系列x射线图像,每个图像包括由所述多个移动不透射线标志物投射的动态生成的水印,其中每个不透射线标志物的位置被示出在所述水印上;
对所述水印进行解码以生成时间戳;以及
使用通过对所述水印进行解码生成的所述时间戳来校正由所述成像设备生成的时间戳。
2.如权利要求1所述的方法,进一步包括:
使得所述多个不透射线标志物中的每个不透射线标志物以各自的频率旋转;以及
处理所述一系列x射线图像以测量所述多个标志物中的每个标志物的旋转位置。
3.如权利要求2所述的方法,其特征在于,对所述水印进行解码以生成时间戳包括:基于所述水印的所述多个标记中的每个标记的所测量的旋转位置和已知的频率,确定所述时间戳。
4.如权利要求1所述的方法,其特征在于,所述水印对具有亚秒精度的所述时间戳进行编码。
5.如权利要求1所述的方法,其特征在于,所述不透射线标志物中的每个不透射线标志物是刻度盘的形式,并且多个刻度盘中的每个刻度盘以由同步到国际时钟标准的控制器确定的速率旋转。
6.如权利要求1所述的方法,其特征在于,每个不透射线标志物是刻度盘的形式,并且其中每个刻度盘连接到多个电机中的一个电机,所述方法进一步包括:
对于每个刻度盘,以与同其他刻度盘连接的电机不同的指定频率旋转连接到所述刻度盘的电机,以便使得每个刻度盘以与所述刻度盘连接到的电机相同的频率旋转。
7.如权利要求1所述的方法,其特征在于,所述不透射线标志物中的每个不透射线标志物是刻度盘的形式,并且以不同频率旋转的至少两个刻度盘用于获得亚秒精度。
8.如权利要求1所述的方法,其特征在于,所述不透射线标志物被配置为旋转。
9.如权利要求1所述的方法,其特征在于,每个不透射线标志物包括可旋转刻度盘。
10.如权利要求1所述的方法,其特征在于,与所述多个不透射线标志物相关联的所述各自的频率彼此不同。
11.如权利要求1所述的方法,其特征在于,由所述成像设备生成的所述时间戳通过以下操作来校正:基于通过对所述水印进行解码生成的时间戳与由所述成像设备生成的时间戳的比较来确定时间戳校正。
12.如权利要求11所述的方法,其特征在于,所述时间戳校正基于跨图像集合的通过对所述水印进行解码生成的时间戳与由所述成像设备生成的时间戳之间的算术平均值差或平均值差。
13.如权利要求11所述的方法,进一步包括将所述时间戳校正用于由所述成像设备生成的时间戳的后续校正。
14.如权利要求1所述的方法,进一步包括使用经校正的时间戳来将图像集合与外部信号同步。
15.如权利要求14所述的方法,其特征在于,所述外部信号是心脏信号。
16.如权利要求1所述的方法,其特征在于,由所述成像设备生成的时间戳嵌入在图像元数据中。
17.一种对x射线图像进行动态时间戳水印标记的***,所述***包括:
多个电机;
不透射线标志物,所述不透射线标志物直接或间接地连接到所述多个电机,所述不透射线标志物被配置为响应于所述多个电机的操作而选择性地旋转;以及
控制器,所述控制器被配置为控制所述多个电机中的每个电机的操作,以使每个电机以各自的频率旋转,从而使得连接到所述电机的不透射线标志物以所述各自的频率旋转;
一个或多个计算机处理器;
一个或多个计算机可读存储介质;
程序指令,所述程序指令存储在所述一个或多个计算机可读存储介质上以供所述一个或多个计算机处理器中的至少一个计算机处理器执行,所述程序指令包括用于以下内容的指令:
处理来自成像设备的一系列x射线图像,每个图像包括由所述不透射线标志物投射的动态生成的水印,其中每个不透射线标志物的位置被示出在所述水印上;
对所述水印进行解码以生成时间戳;以及
使用通过对所述水印进行解码生成的所述时间戳来校正由所述成像设备生成的时间戳。
18.如权利要求17所述的***,其特征在于,每个不透射线标志物使像素相对从x射线源传播到x射线检测阵列上的x射线不透明,以在水印的像素阵列中形成图像。
19.如权利要求17所述的***,其特征在于,所述控制器被同步到国际时钟标准。
20.如权利要求17所述的***,其特征在于,至少两个不透射线标志物用于获得亚秒精度。
21.如权利要求17所述的***,其特征在于,至少三个不透射线标志物用于对所述时间戳进行编码,并且其中所述程序指令进一步包括指令以:
在不使用与所述图像相对应的元数据时间戳的情况下,将所述时间戳转换为具有亚秒精度的时间戳。
22.一种对x射线图像进行动态时间戳水印标记的***,所述***包括:
多个电机;
不透射线标志物,所述不透射线标志物直接或间接地连接到电机,所述不透射线标志物被配置为响应于所述多个电机的操作而选择性地旋转;以及
控制器,所述控制器被配置为控制所述多个电机中的每个电机的操作,以使每个电机以各自的频率旋转,从而使得相对应的不透射线标志物以所述各自的频率旋转;
一个或多个计算机处理器;
一个或多个计算机可读存储介质;
程序指令,所述程序指令存储在所述一个或多个计算机可读存储介质上以供所述一个或多个计算机处理器中的至少一个计算机处理器执行,所述程序指令包括用于以下内容的指令:
处理一系列x射线图像以生成时间戳,其中每个图像包括对与所述时间戳相对应的信息进行编码的动态生成的水印,其中所述动态生成的水印由所述多个移动不透射线标志物投射;
将从所述水印生成的所述时间戳与被包含在图像元数据中的时间戳进行比较;以及
基于从所述水印生成的所述时间戳,校正在人类临床血管造影期间获得的与所述图像相关联的元数据时间戳。
23.一种或多种计算机可读存储介质,所述一种或多种计算机可读存储介质共同具有包含在其中的程序指令,所述程序指令能由计算机执行以使得所述计算机:
获得一系列x射线图像,每个图像包括与定位在成像设备的辐射发射器和辐射检测器阵列之间的多个可移动不透射线标志物相对应的动态生成的水印,其中每个不透射线标志物的旋转位置被投射到所述水印上作为标记,以及其中每个不透射线标志物与各自的频率相关联;对所述水印进行解码以生成时间戳;以及
使用通过对所述水印进行解码生成的所述时间戳来校正由所述成像设备生成的时间戳。
24.如权利要求23所述的计算机可读存储介质,进一步包括测量所述多个标记中的每一个标记的旋转位置,并且其中通过对所述水印进行解码生成的所述时间戳是基于所述水印的所述多个标记中的每个标记的所测量的旋转位置和已知的频率。
25.如权利要求24所述的计算机可读存储介质,其特征在于,通过对所述水印进行解码生成的所述时间戳具有亚秒精度。
26.如权利要求23所述的计算机可读存储介质,其特征在于,所述水印与多个不透射线标志物相对应,所述多个不透射线标志物被配置为根据国际时钟标准旋转。
27.如权利要求23所述的计算机可读存储介质,其特征在于,每个不透射线标志物由电机旋转,所述电机被配置为使得所述不透射线标志物以指定频率旋转。
28.一种或多种计算机可读存储介质,所述一种或多种计算机可读存储介质共同具有包含在其中的程序指令,所述程序指令能由计算机执行以使得所述计算机:
获得一系列x射线图像,每个图像包括与定位在辐射发射器和辐射检测器阵列之间的多个可移动不透射线标志物相对应的动态生成的水印,其中每个不透射线标志物的旋转位置被投射到所述水印上作为标记,以及其中每个不透射线标志物与各自的频率相关联;
确定所述多个标记中的每个标记的旋转位置;
对所述水印进行解码以生成时间戳;以及
基于经解码的时间戳校正元数据时间戳。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201962815476P | 2019-03-08 | 2019-03-08 | |
US62/815,476 | 2019-03-08 | ||
PCT/US2020/021754 WO2020185706A1 (en) | 2019-03-08 | 2020-03-09 | Temporal calibration of an angiographic imaging system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113507887A CN113507887A (zh) | 2021-10-15 |
CN113507887B true CN113507887B (zh) | 2023-12-05 |
Family
ID=72335341
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202080017712.8A Active CN113507887B (zh) | 2019-03-08 | 2020-03-09 | 血管造影成像***的时间校准 |
Country Status (9)
Country | Link |
---|---|
US (1) | US11282202B2 (zh) |
EP (1) | EP3911239A4 (zh) |
JP (1) | JP7360470B2 (zh) |
KR (1) | KR20210124380A (zh) |
CN (1) | CN113507887B (zh) |
AU (1) | AU2020234781B2 (zh) |
CA (1) | CA3130150C (zh) |
GB (1) | GB2595166B (zh) |
WO (1) | WO2020185706A1 (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA3126986A1 (en) | 2019-02-06 | 2020-08-13 | William E. Butler | Spatiotemporal reconstruction of a moving vascular pulse wave from a plurality of lower dimensional angiographic projections |
US11291422B2 (en) | 2019-03-27 | 2022-04-05 | William E. Butler | Reconstructing cardiac frequency phenomena in angiographic data |
WO2020206430A1 (en) | 2019-04-04 | 2020-10-08 | Butler William E | Intrinsic contrast optical cross-correlated wavelet angiography |
WO2023076244A1 (en) * | 2021-10-25 | 2023-05-04 | Lightlab Imaging, Inc. | Temporal data generation with scattered x-ray radiation |
US20230389892A1 (en) * | 2022-06-02 | 2023-12-07 | Canon U.S.A., Inc. | Devices, systems, and methods for automated delay detection between medical-imaging devices |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS53123188A (en) * | 1977-03-30 | 1978-10-27 | American Science & Eng Inc | Xxray image regenerator |
DE2948109A1 (de) * | 1979-11-29 | 1981-06-04 | Decker, Kurt, Prof. Dr.med., 8000 München | Einrichtung zum kennzeichnen der reihenfolge bei einer serie von roentgenaufnahmen |
CN1726871A (zh) * | 2004-07-29 | 2006-02-01 | 西门子公司 | 可视化血管、尤其是心脏冠状血管内沉积物的方法和装置 |
WO2008139167A2 (en) * | 2007-05-11 | 2008-11-20 | Cameron Carpenter | Radiological quality assurance phantom |
EP2587450A1 (en) * | 2011-10-27 | 2013-05-01 | Nordson Corporation | Method and apparatus for generating a three-dimensional model of a region of interest using an imaging system |
CN103260519A (zh) * | 2010-10-20 | 2013-08-21 | 美敦力导航公司 | 用于重建图像投影的方法和设备 |
WO2014162273A1 (en) * | 2013-04-03 | 2014-10-09 | Koninklijke Philips N.V. | Vessel segmentation |
CN108078576A (zh) * | 2016-11-21 | 2018-05-29 | 株式会社东芝 | 医学图像处理的装置、方法、计算机可读程序以及移动对象跟踪装置和放射治疗*** |
Family Cites Families (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4896344A (en) * | 1984-10-15 | 1990-01-23 | Grady John K | X-ray video system |
US4698836A (en) * | 1984-10-24 | 1987-10-06 | Minasian John L | Means and techniques useful in marking X-ray film |
US5090042A (en) * | 1990-12-24 | 1992-02-18 | Bejjani Fadi J | Videofluoroscopy system for in vivo motion analysis |
US6267502B1 (en) * | 1998-04-10 | 2001-07-31 | Minrad Inc. | Alignment verification device and method of using the same with a visual light beam and an x-ray |
DE19936364A1 (de) * | 1999-08-03 | 2001-02-15 | Siemens Ag | Verfahren zur Identifizierung und Lokalisierung von Marken in einem 3D-Volumendatensatz |
JP2003522576A (ja) * | 2000-02-18 | 2003-07-29 | ウィリアム・ボーモント・ホスピタル | 平坦なパネル画像装置を有するコーンビームコンピュータ断層撮像装置 |
CA2487044C (fr) * | 2002-05-23 | 2010-05-11 | Centre National De La Recherche Scientifique (C.N.R.S.) | Dispositif de stereoradiographie et procedure d'utilisation |
US7697974B2 (en) * | 2003-10-10 | 2010-04-13 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus for analysis of angiographic and other cyclical images |
US8768026B2 (en) * | 2003-11-26 | 2014-07-01 | Hologic, Inc. | X-ray imaging with x-ray markers that provide adjunct information but preserve image quality |
JP4758355B2 (ja) * | 2003-12-22 | 2011-08-24 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 患者の体内に医療用機器を案内するためのシステム |
EP1715800A2 (en) * | 2004-02-10 | 2006-11-02 | Koninklijke Philips Electronics N.V. | A method, a system for generating a spatial roadmap for an interventional device and a quality control system for guarding the spatial accuracy thereof |
EP1744675B1 (en) * | 2004-04-29 | 2013-10-30 | Koninklijke Philips N.V. | Viewing system for control of ptca angiograms |
FR2882245B1 (fr) * | 2005-02-21 | 2007-05-18 | Gen Electric | Procede de determination du deplacement 3d d'un patient positionne sur une table d'un dispositif d'imagerie |
US7356120B2 (en) * | 2005-09-23 | 2008-04-08 | Accuray Incorporated | Integrated quality assurance for in image guided radiation treatment delivery system |
US20080025475A1 (en) * | 2006-07-25 | 2008-01-31 | Synarc, Inc. | Apparatus for determining the position and orientation in medical imaging |
US7863897B2 (en) | 2007-09-07 | 2011-01-04 | The General Hospital Corporation | Method and apparatus for characterizing the temporal resolution of an imaging device |
US7643615B2 (en) * | 2007-09-27 | 2010-01-05 | Carestream Health, Inc. | Radiographic marker for inclination angle of imaging receiver |
ES2646819T3 (es) * | 2010-05-12 | 2017-12-18 | Trophy | Aparato de alineación para radiografía intrabucal dental |
US8582847B2 (en) | 2010-10-19 | 2013-11-12 | Siemens Medical Solutions Usa, Inc. | System for verifying data integrity in an X-ray imaging system |
JP5859232B2 (ja) * | 2011-07-08 | 2016-02-10 | 国立大学法人東北大学 | 評価用補助具および評価用装置 |
US9084888B2 (en) * | 2011-11-10 | 2015-07-21 | Varian Medical Systems, Inc. | Systems and methods for segmentation of radiopaque structures in images |
US9936896B2 (en) * | 2012-01-12 | 2018-04-10 | Siemens Medical Solutions Usa, Inc. | Active system and method for imaging with an intra-patient probe |
WO2014124447A1 (en) * | 2013-02-11 | 2014-08-14 | Angiometrix Corporation | Systems for detecting and tracking of objects and co-registration |
US10433911B2 (en) * | 2013-09-18 | 2019-10-08 | iMIRGE Medical INC. | Optical targeting and visualization of trajectories |
US10499813B2 (en) * | 2014-09-12 | 2019-12-10 | Lightlab Imaging, Inc. | Methods, systems and apparatus for temporal calibration of an intravascular imaging system |
WO2016139347A1 (en) | 2015-03-04 | 2016-09-09 | 3Shape A/S | A system and method for increasing the accuracy of a medical imaging device |
US10123761B2 (en) | 2015-07-01 | 2018-11-13 | William E. Butler | Device and method for spatiotemporal reconstruction of a moving vascular pulse wave in the brain and other organs |
US10517561B2 (en) * | 2017-05-17 | 2019-12-31 | Carestream Health, Inc. | Self-calibrating technique for x-ray imaging scanners |
-
2020
- 2020-03-09 AU AU2020234781A patent/AU2020234781B2/en active Active
- 2020-03-09 CN CN202080017712.8A patent/CN113507887B/zh active Active
- 2020-03-09 EP EP20769736.8A patent/EP3911239A4/en active Pending
- 2020-03-09 KR KR1020217028376A patent/KR20210124380A/ko not_active Application Discontinuation
- 2020-03-09 US US16/813,513 patent/US11282202B2/en active Active
- 2020-03-09 CA CA3130150A patent/CA3130150C/en active Active
- 2020-03-09 WO PCT/US2020/021754 patent/WO2020185706A1/en unknown
- 2020-03-09 GB GB2111722.1A patent/GB2595166B/en active Active
- 2020-03-09 JP JP2021549509A patent/JP7360470B2/ja active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS53123188A (en) * | 1977-03-30 | 1978-10-27 | American Science & Eng Inc | Xxray image regenerator |
DE2948109A1 (de) * | 1979-11-29 | 1981-06-04 | Decker, Kurt, Prof. Dr.med., 8000 München | Einrichtung zum kennzeichnen der reihenfolge bei einer serie von roentgenaufnahmen |
CN1726871A (zh) * | 2004-07-29 | 2006-02-01 | 西门子公司 | 可视化血管、尤其是心脏冠状血管内沉积物的方法和装置 |
WO2008139167A2 (en) * | 2007-05-11 | 2008-11-20 | Cameron Carpenter | Radiological quality assurance phantom |
CN103260519A (zh) * | 2010-10-20 | 2013-08-21 | 美敦力导航公司 | 用于重建图像投影的方法和设备 |
EP2587450A1 (en) * | 2011-10-27 | 2013-05-01 | Nordson Corporation | Method and apparatus for generating a three-dimensional model of a region of interest using an imaging system |
WO2014162273A1 (en) * | 2013-04-03 | 2014-10-09 | Koninklijke Philips N.V. | Vessel segmentation |
CN108078576A (zh) * | 2016-11-21 | 2018-05-29 | 株式会社东芝 | 医学图像处理的装置、方法、计算机可读程序以及移动对象跟踪装置和放射治疗*** |
Non-Patent Citations (1)
Title |
---|
MSCT扫描方式对辐射剂量及图像质量影响的研究;兰永树 等;临床放射学杂志;第29卷(第7期);第960-963页 * |
Also Published As
Publication number | Publication date |
---|---|
US11282202B2 (en) | 2022-03-22 |
GB202111722D0 (en) | 2021-09-29 |
AU2020234781B2 (en) | 2023-03-30 |
JP7360470B2 (ja) | 2023-10-12 |
EP3911239A4 (en) | 2022-10-05 |
CN113507887A (zh) | 2021-10-15 |
KR20210124380A (ko) | 2021-10-14 |
GB2595166A (en) | 2021-11-17 |
JP2022523199A (ja) | 2022-04-21 |
CA3130150C (en) | 2024-06-25 |
GB2595166B (en) | 2023-07-12 |
US20200286237A1 (en) | 2020-09-10 |
AU2020234781A1 (en) | 2021-09-09 |
EP3911239A1 (en) | 2021-11-24 |
CA3130150A1 (en) | 2021-09-17 |
WO2020185706A1 (en) | 2020-09-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113507887B (zh) | 血管造影成像***的时间校准 | |
EP3501601B1 (en) | Medical apparatus and computer program. | |
CN102067176B (zh) | 结合局部运动监测、校正和评估的辐射成像 | |
US10376217B2 (en) | Positioning unit for positioning a patent, imaging device and method for the optical generation of a positioning aid | |
US8582856B2 (en) | Image processing apparatus, image processing method, and program | |
US11158046B2 (en) | Estimating measurements of craniofacial structures in dental radiographs | |
US20150297120A1 (en) | Method For Tracking Motion of Subject in Real Time and for Correcting Medical Image | |
US20140371578A1 (en) | Real-time display of vasculature views for optimal device navigation | |
EP2024935A1 (en) | Method and apparatus for reconstructing an image | |
CN104854471A (zh) | 针对磁共振弥散加权成像(dwi)的运动检测与校正方法 | |
CN1895170A (zh) | X射线诊断装置、摄影角度决定装置、程序存储介质和方法 | |
US10827999B2 (en) | Dynamic analysis apparatus and system for measuring temporal changes in blood vessels | |
EP3629295A1 (en) | Tomographic image generation apparatus, method, and program | |
Schaller et al. | Inverse C-arm positioning for interventional procedures using real-time body part detection | |
Konchak et al. | A Pascal computer program for digitizing lateral cephalometric radiographs | |
US11380006B2 (en) | Size measurement using angle-constrained radiographic imaging | |
WO2021123364A1 (en) | Apparatus, system, method and computer program for providing a nuclear image of a region of interest of a patient | |
US10709403B2 (en) | Processing of interventional radiology images by ECG analysis | |
JP2015105841A (ja) | 画像診断支援装置に画像処理を実行させるためのコンピュータプログラム、装置及び方法 | |
WO2007054876A2 (en) | Signal processing unit for producing images | |
US20220343567A1 (en) | Systems and methods for rendering objects translucent in x-ray images | |
CN113781593B (zh) | 四维ct图像的生成方法、装置、终端设备及存储介质 | |
RU2732735C1 (ru) | Способ контроля положения и дыхания пациента с применением набора инфракрасных датчиков глубины | |
CN110546684A (zh) | 时变数据的定量评估 | |
US20230380791A1 (en) | Method for planning a recording of an x-ray image of a patient by an x-ray system, and overall system |
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 |