CN1311786C - 增强包括外物的医疗诊断图像的对比度的方法和装置 - Google Patents
增强包括外物的医疗诊断图像的对比度的方法和装置 Download PDFInfo
- Publication number
- CN1311786C CN1311786C CNB02150590XA CN02150590A CN1311786C CN 1311786 C CN1311786 C CN 1311786C CN B02150590X A CNB02150590X A CN B02150590XA CN 02150590 A CN02150590 A CN 02150590A CN 1311786 C CN1311786 C CN 1311786C
- Authority
- CN
- China
- Prior art keywords
- data
- ray
- value
- anatomical data
- anatomical
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims description 46
- 238000003745 diagnosis Methods 0.000 title claims description 6
- 230000002708 enhancing effect Effects 0.000 title 1
- 238000012546 transfer Methods 0.000 claims abstract description 58
- 210000003484 anatomy Anatomy 0.000 claims abstract description 12
- 238000009826 distribution Methods 0.000 claims abstract description 7
- 230000005764 inhibitory process Effects 0.000 claims description 29
- 238000002224 dissection Methods 0.000 claims description 24
- 238000001514 detection method Methods 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000014509 gene expression Effects 0.000 claims description 8
- 235000012364 Peperomia pellucida Nutrition 0.000 claims 1
- 240000007711 Peperomia pellucida Species 0.000 claims 1
- 238000004422 calculation algorithm Methods 0.000 description 21
- 230000005855 radiation Effects 0.000 description 20
- 230000008569 process Effects 0.000 description 16
- 238000002594 fluoroscopy Methods 0.000 description 7
- 210000000988 bone and bone Anatomy 0.000 description 5
- 230000008859 change Effects 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 230000008676 import Effects 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 210000004872 soft tissue Anatomy 0.000 description 3
- 210000001519 tissue Anatomy 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 230000001915 proofreading effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/90—Dynamic range modification of images or parts thereof
- G06T5/94—Dynamic range modification of images or parts thereof based on local image properties, e.g. for local contrast enhancement
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/40—Image enhancement or restoration using histogram techniques
-
- 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
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Optics & Photonics (AREA)
- Heart & Thoracic Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
提供了一种X-射线系统(100)，该系统利用X-射线源(102)和检测器以在感兴趣区(174)中产生X-射线数据。该X-射线数据包括来自患者(104)的解剖部位的解剖数据、来自患者(104)之外的结构的非解剖数据。解剖和非解剖数据可能具有不相重叠的灰度分布。沿灰度分布识别抑制等级(210)。处理器模块基于抑制等级(210)选择传递函数(212)的最小值并自动地计算传递函数(212)以将X-射线数据转换为显示值。然后显示器基于显示值显示X-射线图像(172)。
Description
技术领域
本发明的某些实施例涉及医疗诊断系统，具体地说涉及调整显示的包括外物的诊断图像的对比度的技术和装置。
背景技术
大家都非常熟悉产生一系列患者的内部图像的X-射线系统，比如心脏学、放射学和荧光学系统。患者曝光在X-射线中，然后在X射线穿过患者之后检测它。发射脉冲辐射以产生实时显示在监视器上的连续序列图像。在X-射线经过患者时X-射线受到了衰减。X-射线所受到的衰减量通过显示的像素的灰度级表示在图像中。在灰度级之间的对比度代表衰减量。
骨头和不同类型的组织对X-射线的衰减量不同，由此能够检测并以不同的对比等级显示在图像监视器上。例如，骨头衰减X-射线的程度比肌肉衰减X-射线的程度更大，因此可能被比周围的人体组织更加暗地显示。仅包含软组织的人体可能比包含了软组织和骨头的人体的对比度范围更小。此外，散射辐射或在很大的患者图像中使用增加的kVp电平也可能降低对比度范围。
通过查询表或传递函数将该系统所检测的辐射级与所显示的图像的对比度关联起来。换句话说，该系统使用传递函数来将具体的辐射级指定到显示器的特定的灰度级。该系统通过改变传递函数的形状(例如斜率、偏移等)改变与灰度级的特定范围相关的所显示的图像的对比度范围。该系统可以具有表示不同的数学模型或形状的多个传递函数，从这些传递函数中选择一个传递函数来控制一个或多个显示图像的对比度。对于不同的过程或者在对不同的解剖体进行成像时选择特定的传递函数。
通过该系统可以设定在一个过程中用于显示多个图像的对比度范围，并且在整个过程中保持恒定。因此，对于一次扫描或一整系列扫描，给所检测的辐射级或特定的亮度级指定特定的灰度级。这并不是很有利，因为在整个相同的过程中，可以扫描在具有不同的对比度范围的患者体内的感兴趣区，因此对比度可能出现波动。一些图像可以具有较高的对比度等级，包括很暗的区域和很亮的区域，而另一些图像可能具有很低的对比度，呈浅或淡白色。操作员可以选择以调整所显示的图像的对比度以校正解剖体的对比度的变化，但是这比较费时、易于出错，并且在检查不同的组织时需要反复进行。
为在整个过程(荧光学过程、心脏学过程、放射学过程或其它的过程)中提供更加恒定的对比度，已经提出了自动对比度补偿算法。自动对比度补偿或对比度管理用于提供具有更好的诊断效用的更加满意的图像。在采集图像时检查它们的最大和最小亮度级。然后使用最大和最小的亮度级确定新的灰度传输函数以增强所显示的图像的对比度。因此，在相同的过程中放射科医师可以看到包含了骨头的图像和仅包含了软组织的图像而勿需手动调整对比度。
不幸的是，自动对比度补偿算法也对不透明材料或仪器敏感，这些不透明材料或仪器可能包含在患者体内或者可能由放射科医师、外科医师或其它人使用，并且可以出现在图像中。在外物比如卡钳或接骨螺钉进入视场中时，在自动对比度补偿启动时它能够识别新的最大和最小亮度级并调整对比度，以给该外物指定可用的灰度的低端或最小亮度级的某一范围。因此，降低了可用于显示感兴趣的解剖组织的灰度范围，所显示的解剖数据反映了对比度的突降。
操作员通过停止自动对比度选项可以补偿仪器的外观。然后通过将它移到左或右(实质上改变了最大和最小值)或者通过改变斜率可以改变传递函数，这种斜率可以增加或降低灰度范围。不幸的是，如果随后从视场中移去仪器，或者如果解剖数据的对比度变化了，则需要通过再次启动自动对比度选项或通过手动改变传递函数来再次调整该图像。在调整传递函数时，这可能比较费时，并且使患者和放射医师或技术人员受到过量的X-射线辐射。
因此，本行业中需要一种用于调整所显示的诊断图像的对比度以解决上述和先前所遇到的问题的方法和装置，该图像包括外物。
发明内容
根据至少一种实施例，提供了一种X-射线系统，该系统利用X-射线源和检测器以在感兴趣区中产生X-射线数据。该X-射线数据包括来自患者的解剖部位的解剖数据、来自患者之外的结构的非解剖数据。解剖和非解剖数据具有至少部分不相重叠的灰度分布。输入端沿灰度分布识别抑制等级(reject level)。处理器模块基于该抑制等级选择传递函数的最小值并自动地计算传递函数以将X-射线数据转换为显示值。然后显示器基于显示值显示X-射线图像。
根据至少一种实施例，提供了一种在医疗诊断装置中的自动对比度控制单元。该图像处理器用于产生在感兴趣区中的解剖和非解剖结构的采样数据。解剖和非解剖数据具有至少部分不相重叠的灰度分布。输入端沿采样数据的灰度分布识别抑制等级。条带图模块产生解剖和非解剖数据的灰度值的条带图。自动对比度模块部分基于该抑制等级选择传递函数的最小值并自动地计算传递函数以将采样数据转换为显示值。
根据至少一种实施例，提供了一种增强诊断图像的对比度的方法，将患者曝光在X-射线中，然后检测包括解剖和非解剖数据的X-射线数据。解剖和非解剖数据具有至少部分不相重叠的灰度分布。输入沿灰度分布的抑制等级。自动地计算基于该抑制等级的最小值。基于该最小值计算用于将X-射线数据转换为显示值的传递函数。基于该显示值显示X-射线图像。
附图说明
附图1所示为根据本发明的一种实施例利用自动对比度补偿的荧光X-射线系统。
附图2所示为根据本发明的一种实施例的感兴趣区和诊断图像。
附图3所示为根据本发明的一种实施例从感兴趣区中产生的传递函数和条带图。
附图4所示为根据本发明的一种实施例包括外物的感兴趣区和诊断图像。
附图5所示为根据一种本发明在包括外物的感兴趣区中从所检测的辐射中产生的传递函数和两个条带图。
附图6所示为根据本发明的一种实施例利用抑制等级的传递函数和两个条带图。
附图7所示为根据本发明的一种实施例的控制面板。
附图8所示为根据本发明的一种实施例使用自动对比度补偿和抑制等级增强所采集的诊断图像的对比度的方法。
在阅读本说明书时结合附图可以更好地理解前述的概述以及下文对本发明的实施例的详细描述。然而，应该理解的是，本发明并不限于在附图中所示的结构和手段。
具体实施方式
附图1所示为根据本发明的一种实施例利用自动对比度补偿的荧光X-射线系统。荧光检查系统100包括朝患者104集中X-射线束的X-射线源102。患者104可以是人、动物、测试模型或包括解剖结构的其它的目标。图像增强器106检测穿过患者104的辐射并产生图像。通过光学器件108将图像传递到CCD视频相机110。CCD视频相机110产生模拟视频信号，通过数字转换器(ADC)112将该模拟视频信号转换为数字信号。可任选择的是，ADC112可以包括在CCD视频相机110中。可替换的是，平板检测器可以替代图像增强器106、光学器件108和CCD视频相机110。然后将数字信号传送到图像处理器114。
荧光检查系统100包括感兴趣区(ROI)门116、条带图发生器118和自动对比度算法器120，这些部件协同工作以执行自动对比度补偿。为了清楚起见在附图1中分别显示了ROI门116、条带图发生器118和自动对比度算法器120，但是它们可以包括在图像处理器114或荧光检查系统100的其它部件中。操作员利用控制面板128来输入数据，比如选择感兴趣区和抑制等级，如下文所讨论。自动对比度算法器120的输出传输到显示器查询表(LUT)122。该显示器LUT122使用通过自动对比度算法器120所识别的参数并将来自图像处理器114的数字图像数据中所检测的辐射的等级相关到所选择的查询表或传递函数。数字到模拟转换器(DAC)124将显示器LUT122的输出从数字转换为模拟并将模拟信号传输到显示器126。然后显示器126显示患者104的诊断图像。
附图2所示为根据本发明的一种实施例的感兴趣区和诊断图像。监视器140或类似的装置显示诊断图像142比如通过显示器126所显示的图像。感兴趣区门116利用预定的参数或操作员通过控制面板128输入的数据识别感兴趣区144。感兴趣区144是诊断图像142的一部分，并且通常是诊断图像142的直径的百分之60至百分之100。例如，感兴趣区144可以固定在百分之90。条带图发生器118使用在感兴趣区144内的诊断数据以建立条带图。
附图3所示为根据本发明的一种实施例从感兴趣区中产生的传递函数和条带图。条带图发生器118利用从感兴趣区144的图像处理器114输出的在原始数据范围156上的数字图像数据以产生条带图150。自动对比度算法器120在数字图像数据的条带图150内确定所检测的最小值(检测的MIN)152和最大值(MAX)154。显示器LUT122利用所检测的MIN152和MAX154值以在所显示的对比度范围158上产生传递函数160。
仅仅作为举例，所显示的对比度范围158可以划分为256个离散的灰度级，其中0表示黑色，而255表示白色。在典型的荧光检查系统中，感兴趣区144的平均值可以显示在所显示的对比度范围158的上半部分中。例如，显示的图像的较亮的单元的值的范围为180至200，而所显示的最黑的单元的范围可以是50至100。
条带图发生器118可以包括存储在存储器(未示)中的值表，将数字图像数据与之进行比较。可替换的是，条带图发生器118可以使用其它的电子电路以检测在数字图像数据中的各种白色、黑色和灰度级。可以使用峰值白色(最大)和峰值黑色(最小)信号以产生驱动模拟电路的偏移和增益值。这种电路可以与例如模拟门电路或削波电路结合，这种模拟门电路或削波电路排除在操作员通过控制面板128所识别的某一等级之下的黑色值。
条带图发生器118在图表上绘制条带图150，水平轴代表在原始数据范围156上的数字图像数据在感兴趣区144内的所检测的辐射的等级，而垂直轴代表在显示器126上可得到的所显示的对比度范围158。所显示的对比度范围158可以是从黑色到白色的灰度。可替换的是，所显示的对比度范围158可以通过色标表示，在色标中每种色彩表示信号值或所检测的辐射的一种等级或等级带。条带图发生器118可以绘制在条带图150中的每个像素的值。可选择的是，数字图像数据可以划分为像素组，比如10个像素的方块，可以绘制每组像素的平均值。条带图发生器118可以绘制所采集的每单个图像或所采集的一定的百分比的图像的条带图150。
一旦产生了条带图150，自动对比度算法器120确定了在原始数据范围156内的所检测的MIN152和MAX154值。所检测的MIN152是在感兴趣区144内的最低的检测值或者是在感兴趣区144内具有最小的检测辐射的点。具有所检测的MIN152或更小的检测的辐射的像素可以通过在显示器126上的灰度(即，黑色或接近黑色)的最黑单元表示。MAX154是在感兴趣区144中的最高的检测值，并表示具有最大的检测辐射的像素。具有等于或高于MAX154的检测辐射的像素可以通过在显示器126上的最亮的单元(即，白色)表示。
此外，可以进一步对条带图150进行滤波以消除差异值。例如，为了识别所检测的MIN152，自动对比度算法器120可以设定必须满足的条件，比如在条带图150中绘制的行中具有大于预定数量的像素的10个值。可以在上端设置类似的限制以确定MAX154。因此，可以不考虑在条带图150的最小和最大端上存在的噪声。可替换的是，自动对比度算法器120可以设定如下的条件，使在条带图150接近并组合1％的幅值的阈值时条带图150的坡度或斜率可以在预定的水平之下。可选择的是，在识别所检测的MIN152和MAX154之前可以利用卷积核或窗口以通过样条或三次样条函数以对条带图150进行平滑。应该理解的是，存在识别在条带图150中的所检测的MIN152和MAX154值的进一步方法，并可以利用任何适合的其它方法。
在已经识别了所检测的MIN152和MAX154之后，显示器LUT122绘制从所检测的MIN152到MAX154的传递函数160并转换条带图150的辐射等级以占满显示器126的全部对比度范围。因此，在垂直轴上通常给所检测的MIN152分配零或接近零的值，而在垂直轴上给MAX154分配255的值。所示的传递函数160为非线性的，但是应该理解的是可以使用任何传递函数。例如，传递函数可以是线性(比如在附图5中所示的传递函数184)、对数、指数、S-形、双曲线切线或可以不拟合任何数学模型的另一形状。
附图4所示为根据本发明的一种实施例包括外物的感兴趣区和诊断图像。附图4并入了与在附图2中所示的部件类似的部件，比如监视器170、诊断图像172和感兴趣区174。此外，示出了外物176。外物176是不透明的，因此透射很少的辐射或不透射辐射。外物176可以是金属工具或零件，比如卡钳、铅手套、钢板、刀片或螺钉。外物176包括在诊断图像172中，因此诊断图像172包括解剖和非解剖数据。
附图5所示为根据本发明的一种实施例在包括外物的感兴趣区中从所检测的辐射中产生的传递函数和两个条带图。条带图发生器118利用来自感兴趣区174的原始数据范围194以建立条带图180和182。条带图180是外物176的非解剖数据的第一分布。所检测的MIN186是在原始数据范围194中所检测的最低的辐射级。条带图182是类似于在附图3的条带图150的患者104的解剖数据的第二分布。
在外物176引入到感兴趣区174中时，如附图4所示，自动对比度算法器120将所检测的MIN186识别为最低的可检测的辐射。显示器LUT122基于具有所显示的对比度范围192的检测的MIN186和MAX190产生传递函数。因此，所显示的解剖数据可以反映对比度的突降，因为可用于显示来自患者104的解剖数据(识别为在MIN188和MAX190之间的原始数据范围194)的对比度范围196小于所显示的对比度范围192。通过传递函数184的较低的等级表示非解剖数据，并使解剖数据具有传递函数184的较高的等级。
附图6所示为根据本发明的一种实施例的利用抑制等级的传递函数和两个条带图。所示的条带图200表示从外物176中所检测的非解剖数据，而条带图202表示从患者104中所检测的解剖数据，与附图5的条带图180和182类似。自动对比度算法器120确定对应于条带图202的MIN204和MAX206和对应于条带图200的检测的MIN208。正如下文所讨论，操作员通过控制面板128可以输入抑制等级210。
附图6阐述了点A和B和非解剖MAX218，它示出了外物176的最大的检测的辐射。下文结合点A和B和非解剖MAX218描述采集包括非解剖和解剖结构的图像的实例性方法。
通过单独利用自动对比度选项，可以产生附图5的传递函数184。然而，荧光检查系统100进一步包括用户控制的最小值，自动对比度算法器120不使用在原始数据范围216中的在该最小值之下的数据来确定最小和最大值。用户控制的最小值(抑制等级210)可以与自动对比度特征同时使用。
附图7所示为根据本发明的一种实施例的控制面板。控制面板220可以包括启动功能的软键，比如“启动自动对比度”222、“启动抑制等级”224、“增加抑制等级”226和“降低抑制等级”228。可替换的是，可以通过如下的方式选择通过软键222-228所访问的功能：在与荧光检查系统100相关的一个或多个附加控制面板上的一个或多个按键、跟踪球、开关或键盘。可替换的是，还可以通过硬件实现对“启动抑制等级”224、“增加抑制等级”226和“降低抑制等级”228的选择。例如，DAC124之后可以具有包括带有可变偏移和增益的放大器的电路。可以使用该电路来选择抑制等级210以改变在显示器126和其它的成像装置(比如打印机)上所显示的对比度范围214。
抑制等级210可以沿着原始数据范围216移动到任一点。可选择的是，抑制等级210可以移动到在非解剖MAX218和MIN204之间的值上以便自动对比度算法器120仅仅使用在抑制等级210之上的数字图像数据来确定最小和最大值。换句话说，通过使用抑制等级210来将非解剖数据和解剖数据的检测值分开，所显示的对比度范围214可以基于解剖数据的MIN204和MAX206。然而，如果抑制等级210移动到MIN204的右边，在患者104的解剖数据将受到影响，“剪切掉”部分图像或将其显示为黑色。如果抑制等级210移动到非解剖MAX218的左边，则使用一部分灰度级来将该部分的条带图部分200显示在抑制等级210的右边，由此降低了可用于解剖数据的灰度级。
虽然在引入了外物176或其它的不透明工具时可以建立可选择的一个抑制等级210，但是工具和解剖部位具有不同的密度级，因此透射不同的辐射级。开始可以将抑制等级210设定在检测的MIN208，然后操作员可以使用增加抑制等级226来调整抑制等级210，同时监测在显示器126上的诊断图像172。可替换的是，对大部分研究，开始可以将抑制等级210设定平均MIN204之下所确定的值上。
在显示诊断数据时使用黑色和白色的最大灰度级值不可取。因此，自动对比度算法器120可以使用适当的值来缓和由显示器LUT122所使用的MAX206和MIN204值。结果，在抑制等级210之下的值可以指定为黑色值或比黑色更黑的值，将在抑制等级210之下的数据从在抑制等级210之上的解剖数据中区别出来。
在外物不是非常明显地比患者104更加不透明的情况下，操作员可以选择通过抑制等级210切去一部分外物176。因此，虽然在灰影中可看见一部分外物176，并且不能如上文所讨论地那样以黑色整个地表示，但是可使用更大范围的显示对比度范围214来显示由条带图202所表示的解剖数据。
在采集荧光图像的同时，操作员可以选择“启动自动对比度”222以便自动调整所显示的图像的对比度。在将外物176引入到感兴趣区174的点上，操作员可以选择“启动抑制等级”224。然后操作员可以使用“增加抑制等级”226和/或“降低抑制等级”228以将抑制等级210调整到由条带图200所表示的非解剖数据之上的值。结果，解剖数据的最小值(MIN204)变为由自动对比度算法器120所确定的最小值，显示器LUT122使用该最小值来计算传递函数212。因此，在需要的时候(例如，在使用外物176时)可以引入抑制等级210，然后在不需要时(例如在外物176从感兴趣区176中取去时)关掉抑制等级210，尽管在自动对比度选项继续自动地调整所显示的图像时不是必需要停止抑制等级210。因此，即使外物176重复进入并离开感兴趣区174，所显示的解剖图像也不会由于外物176的缘故而改变对比度，因此不需要手动调整显示器对比度214。
附图8所示为根据本发明的一种实施例使用自动对比度补偿和抑制等级增强所采集的诊断图像的对比度的方法。对于附图8的方法，启动自动对比度选项。例如，通过选择软键“启动自动对比度”222或通过键盘或如前文所讨论的其它的用户输入可以启动该选项。预先定义显示器LUT122将其用于计算传递函数212的查询表。此外，通过用户输入或通过访问预先定义的值识别感兴趣区174。感兴趣区174识别由条带图发生器188所使用的原始数据范围216。
在步骤250中，X-射线源102将患者104曝光在X-射线辐射中。X-射线的水平可以由例如正成像的解剖部位、过程类型或患者102的尺寸所确定。
在步骤252中，使用图像增强器106、光学器件108、视频相机110、ADC112和图像处理器114采集可以作为诊断图像172显示的数字图像数据。数字图像数据可以存储在存储装置(未示)中以在后面进一步处理。
在步骤254中，条带图发生器188基于原始数据范围194和预先定义的感兴趣区174产生一个或多个条带图。例如，如果感兴趣区174包括如在附图4中所示的解剖和非解剖数据，在产生条带图180和182。如果感兴趣区174仅包括如在附图2中所示的解剖数据，则产生在附图3的条带图150。如前文所述，条带图发生器118可以为每个所采集的图像产生一个或多个条带图或者将一个以上的采集的图像组合以产生条带图。
在步骤256中，自动对比度算法器120确定在感兴趣区174内的原始数据范围194的最小和最大值。例如，如果在感兴趣区174中不存在外物176，则可以确定如在附图3中的检测的MIN152和MAX154。如果存在外物176，在可以确定如在附图5中的检测MIN186和MAX190。
接着，在步骤258中，显示器LUT122利用预先确定的查询表和步骤256的最大和最小值以计算传递函数184。传递函数184可以将条带图180和182与灰度级或所显示的对比度范围192关联起来。
在步骤260中，显示器126通过将传递函数184应用到来自图像处理器114的数字图像数据中来显示图像172。步骤256的最小的检测值可以显示为黑色，而最大的检测值可以显示为白色。如前文所讨论，还可以使用色标。如果外物176存在于感兴趣区174中，则可以以小于所显示的对比度范围192的对比度范围196显示解剖数据。因此，解剖数据的对比度可以降低而不会使图像模糊。
在步骤262中，该方法确定是否要采集更多的图像。如果是，则流程进行到步骤264中。如果否，则流程进行到步骤284，该方法结束。
在步骤264中，该方法确定是否启动抑制等级选项。如果否，该方法重复步骤250至262，使用感兴趣区174来产生新的条带图180和182，确定检测的MIN186和MAX190，并计算相关的传递函数184。如果操作员启动抑制等级选项，比如软键“启动抑制等级”224，则流程进行到步骤266。
在步骤266中，图像处理器144设定开始的抑制等级210。抑制等级210可以设定为等于或稍稍大于所检测的MIN208的等级。在步骤268中，与步骤250类似，X-射线源102将患者104曝光在X-射线辐射脉冲中。接着，在步骤270中，使用图像增强器106、光学器件108、视频相机110、ADC112和图像处理器114采集可以作为诊断图像172显示的数字图像数据。在步骤272中，与步骤254类似，条带图发生器188产生一个或多个条带图200和202。条带图发生器118再次利用先前定义的感兴趣区174。
在步骤274中，自动对比度算法器120识别抑制等级210。抑制等级210可以处于在步骤266中所设定的初始值。可替换的是，如先前所讨论，操作员可能已经通过使用软键“增加抑制等级”226和/或“降低抑制等级”228或通过使用键盘或跟踪球调整了抑制等级210。
在步骤276中，自动对比度算法器120确定MIN204和MAX206。自动对比度算法器120仅考虑在原始数据范围216中大于抑制等级210的数据。例如，如果抑制等级210设定在点A上，则最小值也设定在点A上。在点A之下与条带图200和202相关的任何数据都将被指定最低的灰度值，比如黑色或比黑色更黑的色。如果抑制等级210设定在点B上，则最小值将设定在MIN204上，因为自动对比度算法器120查找在抑制等级210之上的在原始数据范围216中的最低检测值。如前文所讨论，可以对条带图200和202进行滤波以消除差值。
在步骤278中，显示器LUT122利用在步骤276中所确定的最小和最大值并根据预先定义的查询表计算传递函数212。在步骤280中，显示器126显示图像172。如果由条带图200所代表的外物176位于感兴趣区174中，并且抑制等级210位于在非解剖MAX218和MIN204之间的任何位置，则可以将外物176显示为黑色。由条带图202所表示的解剖数据几乎利用所有的可用的显示对比度范围214。因此，使用与引入外物176之前的相同的显示对比度范围214来显示解剖数据。在步骤282中，该方法确定是否采集更多的图像。如果否，则流程进行到步骤284，该方法终止。如果是，则流程返回到步骤268。
通过使用前述的方法和装置来采集包括解剖和非解剖结构的诊断图像，可以采集具有改善的诊断用途的图像。因此，在引入了不透明物体时，利用抑制等级210连同自动对比度的荧光检查系统100能够使放射医师利用自动对比度的优点实施荧光检查程序，而不会造成用于显示解剖结构的对比度减低的缺陷。
虽然参考至少一种实施例已经描述了本发明，但是本领域的普通技术人员应该理解的是在不脱离本发明的范围的前提下可以作出各种变化和等效替代。此外，在本发明的教导下可以作出各种修改以适应特定的情况或材料而不脱离本发明的范围。因此，本发明并不限于所公开的具体的实施例，而是本发明包括落在附加的权利要求的范围中的所有的实施例。
Claims (29)
1.一种X-射线系统(100)，包括：
X-射线源(102)和检测器，该检测器产生具有在与感兴趣区(174)相关的X-射线数据范围上的灰度分布的X-射线数据，所说的X-射线数据包括表示在感兴趣区(174)中的患者(104)的解剖结构的解剖数据和包括在患者之外的感兴趣区(174)中的结构的非解剖数据，所说的解剖和非解剖数据具有至少部分不相重叠的灰度分布的灰度值；
沿所说的X-射线数据的所说的灰度分布识别抑制等级(210)的输入端；
处理器模块，该处理器模块自动地计算传递函数(212)以将X-射线数据转换为形成X-射线图像(172)的显示值，所说的处理器模块至少部分基于所说的抑制等级(210)选择所说的传递函数(212)的最小(204)值；以及
显示器(126)，该显示器基于所说的显示值显示所说的X-射线图像(172)。
2.权利要求1所述的X-射线系统(100)，其中传递函数(212)构成了在所说的解剖数据的灰度值和形成所说的X-射线图像(172)的显示值之间的一对一的相关。
3.权利要求1所述的X-射线系统(100)，其中解剖数据的灰度分布具有短于显示器范围的解剖灰度范围，所说的显示器值分布在显示器范围上，所说的传递函数(212)扩展所说的解剖数据以整个地跨越所说的显示器范围。
4.权利要求1所述的X-射线系统(100)，其中所说的处理器模块识别所说的解剖数据的最小和最大灰度值，并基于这些最小和最大灰度值计算所说的传递函数(212)。
5.权利要求1所述的X-射线系统(100)，其中所说的处理器模块产生所说的解剖和所说的非解剖数据的所说的灰度值的条带图(200，202)并至少部分地基于所说的条带图(200，202)自动地计算所说的传递函数(212)。
6.权利要求1所述的X-射线系统(100)，其中所说的传递函数(212)是线性的。
7.权利要求1所述的X-射线系统(100)，其中所说的传递函数(212)是非线性的。
8.权利要求1所述的X-射线系统(100)，其中所说的非解剖数据的所说的灰度分布与所说的解剖数据的所说的灰度分布分开并相区别。
9.权利要求1所述的X-射线系统(100)，其中所说的输入端将所说的抑制等级(210)识别在所说的解剖和非解剖数据的灰度分布之间。
10.权利要求1所述的X-射线系统(100)，其中，在所说的抑制等级(210)小于所说的解剖数据的最小灰度值时，所说的处理器模块基于独立于所说的抑制等级(210)的所说的解剖数据的所说的最小灰度值自动地设定所说的传递函数(212)的最小(204)值。
11.权利要求1所述的X-射线系统(100)，其中，在所说的抑制等级(210)大于所说的解剖数据的最小灰度值时，所说的处理器模块基于独立于该最小灰度值的抑制等级(210)自动地设定所说的传递函数(212)的最小(204)值。
12.一种在医疗诊断装置中的自动对比度控制单元，包括：
图像处理器(114)，该图像处理器(114)产生具有在与感兴趣区(174)相关的采样数据范围上的灰度分布的采样数据，所说的采样数据包括表示在感兴趣区(174)中的患者(104)的解剖结构的解剖数据和表示在患者(104)之外的感兴趣区(174)中的结构的非解剖数据，所说的解剖和非解剖数据具有至少部分不相重叠的灰度分布的灰度值；
沿所说的采样数据的所说的灰度分布识别抑制等级(210)的输入端；
产生所说的解剖和非解剖数据的灰度值的条带图(200，202)的条带图模块；以及
自动对比度模块，自动对比度模块自动地计算传递函数(212)以将采样数据转换为形成诊断图像(172)的显示值，所说的自动对比度模块至少部分基于所说的抑制等级(210)选择所说的传递函数(212)的最小(204)值。
13.权利要求12所述的自动对比度控制单元，其中该输入端进一步包括软键。
14.权利要求12所述的自动对比度控制单元，其中该输入端进一步包括手动拨号盘。
15.权利要求12所述的自动对比度控制单元，其中所说的医疗诊断装置是X-射线系统(100)，其中所说的采样数据是X-射线数据。
16.权利要求12所述的自动对比度控制单元，其中所说的非解剖数据的所说的灰度分布与所说的解剖数据的所说的灰度分布是分开的并相区别。
17.权利要求12所述的自动对比度控制单元，其中所说的输入端将所说的抑制等级(210)识别在所说的解剖和非解剖数据的灰度分布之间。
18.权利要求12所述的自动对比度控制单元，其中，在所说的抑制等级(210)小于所说的解剖数据的最小灰度值时，所说的处理器模块基于独立于所说的抑制等级(210)的所说的解剖数据的所说的最小灰度值自动地设定所说的传递函数(212)的最小(204)值。
19.权利要求12所述的自动对比度控制单元，其中，在所说的抑制等级(210)大于所说的解剖数据的最小灰度值时，所说的处理器模块基于独立于该最小灰度值的抑制等级(210)自动地设定所说的传递函数(212)的最小(204)值。
20.一种增强诊断图像(172)的对比度的方法，该方法包括：
将患者(104)曝光在X-射线中；
检测包括表示患者(104)的解剖结构的解剖数据和表示在患者(104)之外的结构的非解剖数据，所说的解剖和非解剖数据具有至少部分不相重叠的灰度分布的灰度值；
手动输入抑制等级(210)，所说的抑制等级(210)沿所说的灰度分布设定；
基于所说的抑制等级(210)自动计算沿着所说的至少部分不相重叠的灰度分布的最小(204)值；
基于所说的最小(204)值自动计算传递函数(212)以将所说的X-射线数据转换为显示值；以及
基于所说的显示值显示X-射线图像(172)。
21.权利要求20所述的方法，进一步包括基于所说的X-射线数据产生条带图(200，202)，所说的条带图(202，202)包括表示所说的非解剖数据的第一分布部分和表示所说的解剖数据的第二分布部分，所说的第一和第二分布部分全部分开并彼此不同。
22.权利要求20所述的方法，进一步包括：
基于所说的X-射线数据产生条带图(200，202)，所说的条带图(200，202)包括表示所说的非解剖数据的第一分布和表示所说的解剖数据的第二分布；以及
将第一分布部分作为单个灰度值显示。
23.权利要求20所述的方法，进一步包括以单个的灰度值显示在所说的抑制等级(210)之下的所说的X-射线数据。
24.权利要求20所述的方法，进一步包括将所说的抑制等级(210)设定为小于所说的最小(204)值。
25.权利要求20所述的方法，进一步包括将所说的抑制等级(210)设定为与所说的最小(204)值基本相同。
26.权利要求20所述的方法，其中在所说的患者(104)之外的所说的结构不透明。
27.权利要求20所述的方法，进一步包括检测包括解剖数据的X-射线数据，其中表示基于包括解剖和非解剖数据的所说的X-射线数据的解剖数据的所说的灰度分布基本类似于包括解剖数据的所说的X-射线数据的所说的灰度分布。
28.权利要求20所述的方法，其中所说的最小(204)值基于表示所说的解剖数据的分布。
29.权利要求20所述的方法，其中进一步包括计算最大(206)值，其中所说的最大(206)和最小(204)值基于在感兴趣区(174)中所说的X-射线数据。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US09/990,643 US6574307B1 (en) | 2001-11-21 | 2001-11-21 | Method and apparatus for enhancing the contrast of a medical diagnostic image that includes foreign objects |
| US09/990643 | 2001-11-21 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN1419892A CN1419892A (zh) | 2003-05-28 |
| CN1311786C true CN1311786C (zh) | 2007-04-25 |
Family
ID=25536371
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNB02150590XA Expired - Fee Related CN1311786C (zh) | 2001-11-21 | 2002-11-21 | 增强包括外物的医疗诊断图像的对比度的方法和装置 |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US6574307B1 (zh) |
| EP (1) | EP1324266B1 (zh) |
| JP (1) | JP2003198938A (zh) |
| KR (1) | KR100949225B1 (zh) |
| CN (1) | CN1311786C (zh) |
| DE (1) | DE60235293D1 (zh) |
Families Citing this family (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1298174C (zh) * | 2002-04-19 | 2007-01-31 | 汤姆森特许公司 | 阴极射线管的辅助线圈驱动电路 |
| JP4323770B2 (ja) * | 2002-10-18 | 2009-09-02 | キヤノン株式会社 | 画像処理装置、画像処理方法、プログラム及び記録媒体 |
| DE102004030833A1 (de) * | 2004-06-25 | 2006-01-26 | Siemens Ag | Röngtendiagnostikverfahren und zugehörige Vorrichtung |
| US20080144901A1 (en) * | 2006-10-25 | 2008-06-19 | General Electric Company | Cartoon-like exaggeration of medical images to emphasize abnormalities |
| US7912263B2 (en) * | 2007-08-06 | 2011-03-22 | Carestream Health, Inc. | Method for detecting clipped anatomy in medical images |
| JP4823996B2 (ja) * | 2007-11-02 | 2011-11-24 | 日本電産トーソク株式会社 | 輪郭検出方法及び輪郭検出装置 |
| US20110123074A1 (en) * | 2009-11-25 | 2011-05-26 | Fujifilm Corporation | Systems and methods for suppressing artificial objects in medical images |
| US20130085703A1 (en) * | 2011-09-30 | 2013-04-04 | Broadcom Corporation | Histogram-Based Linearization of Analog-to-Digital Converters |
| KR101313295B1 (ko) * | 2011-10-20 | 2013-09-30 | 주식회사 메도니카 | 엑스레이 영상 데이터 획득 방법 및 이를 위한 시스템 |
| KR101474768B1 (ko) * | 2011-12-21 | 2014-12-19 | 삼성전자 주식회사 | 의료기기 및 이를 이용한 영상표시방법 |
| JP5643889B2 (ja) * | 2013-10-24 | 2014-12-17 | 株式会社東芝 | 乳房ｘ線撮影装置及び画像処理装置 |
| KR102393294B1 (ko) * | 2014-09-26 | 2022-05-03 | 삼성전자주식회사 | 의료 영상 장치 및 의료 영상 장치의 제어 방법 |
| US11193898B1 (en) | 2020-06-01 | 2021-12-07 | American Science And Engineering, Inc. | Systems and methods for controlling image contrast in an X-ray system |
| EP4157093A4 (en) * | 2020-06-01 | 2024-01-24 | American Science & Eng Inc | SYSTEMS AND METHODS FOR CONTROLLING IMAGE CONTRAST IN AN X-RAY SYSTEM |
| CN116468729B (zh) * | 2023-06-20 | 2023-09-12 | 南昌江铃华翔汽车零部件有限公司 | 一种汽车底盘异物检测方法、系统及计算机 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5485501A (en) * | 1993-09-10 | 1996-01-16 | Siemens Aktiengesellschaft | Method for the operation of an automatic x-ray exposure unit |
| US5574764A (en) * | 1995-06-06 | 1996-11-12 | General Electric Company | Digital brightness detector |
| US6048311A (en) * | 1997-05-07 | 2000-04-11 | Washburn; Michael J. | Method and apparatus for ultrasound imaging using adaptive gray mapping |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5436829A (en) | 1992-11-05 | 1995-07-25 | General Electric Company | Method of achieving reduced dose X-ray fluoroscopy by employing transform-based estimation of Poisson noise |
| JPH0767864A (ja) * | 1993-09-07 | 1995-03-14 | Toshiba Medical Eng Co Ltd | Ｘ線診断装置 |
| JPH07239933A (ja) * | 1994-03-01 | 1995-09-12 | Toshiba Medical Eng Co Ltd | 画像処理装置 |
| JP3785817B2 (ja) | 1998-07-17 | 2006-06-14 | コニカミノルタホールディングス株式会社 | 放射線画像の画像処理装置 |
| US6633657B1 (en) | 1999-07-15 | 2003-10-14 | General Electric Company | Method and apparatus for controlling a dynamic range of a digital diagnostic image |
-
2001
- 2001-11-21 US US09/990,643 patent/US6574307B1/en not_active Expired - Lifetime
-
2002
- 2002-11-20 KR KR1020020072362A patent/KR100949225B1/ko not_active IP Right Cessation
- 2002-11-21 DE DE60235293T patent/DE60235293D1/de not_active Expired - Lifetime
- 2002-11-21 CN CNB02150590XA patent/CN1311786C/zh not_active Expired - Fee Related
- 2002-11-21 JP JP2002337731A patent/JP2003198938A/ja active Pending
- 2002-11-21 EP EP02258045A patent/EP1324266B1/en not_active Expired - Lifetime
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5485501A (en) * | 1993-09-10 | 1996-01-16 | Siemens Aktiengesellschaft | Method for the operation of an automatic x-ray exposure unit |
| US5574764A (en) * | 1995-06-06 | 1996-11-12 | General Electric Company | Digital brightness detector |
| US6048311A (en) * | 1997-05-07 | 2000-04-11 | Washburn; Michael J. | Method and apparatus for ultrasound imaging using adaptive gray mapping |
Also Published As
| Publication number | Publication date |
|---|---|
| EP1324266B1 (en) | 2010-02-10 |
| DE60235293D1 (de) | 2010-03-25 |
| KR20030041829A (ko) | 2003-05-27 |
| EP1324266A3 (en) | 2006-08-09 |
| EP1324266A2 (en) | 2003-07-02 |
| US20030095626A1 (en) | 2003-05-22 |
| CN1419892A (zh) | 2003-05-28 |
| KR100949225B1 (ko) | 2010-03-24 |
| JP2003198938A (ja) | 2003-07-11 |
| US6574307B1 (en) | 2003-06-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN1311786C (zh) | 增强包括外物的医疗诊断图像的对比度的方法和装置 | |
| DE69629445T2 (de) | Automatische Tonskalenabstimmung mittels Bildaktivitätsmessungen | |
| Zimmerman et al. | An evaluation of the effectiveness of adaptive histogram equalization for contrast enhancement | |
| US6614877B2 (en) | Method and apparatus for enhancing the contrast of a medical diagnostic image acquired using collimation | |
| US6415015B2 (en) | Method and system of compensation of thickness of an organ | |
| US7869637B2 (en) | Histogram calculation for auto-windowing of collimated X-ray image | |
| CN101930595A (zh) | 图像处理方法和图像处理设备 | |
| US6449390B1 (en) | Image processing apparatus and method therefor | |
| US7747058B2 (en) | Image processing method for windowing and/or dose control for medical diagnostic devices | |
| US6542628B1 (en) | Method for detection of elements of interest in a digital radiographic image | |
| EP2192508A1 (en) | Method and system for rendering of diagnostic images on a display | |
| US11213268B2 (en) | X-ray system with computer implemented methods for image processing | |
| US4891757A (en) | Medical imaging system and method having gray scale mapping means for indicating image saturation regions | |
| CN102169530A (zh) | 乳房造影图像的图像处理方法 | |
| CN106333698A (zh) | 运行用于检查对象的x射线系统 | |
| EP1253557A1 (en) | Apparatus for and method of generating an enhanced contrast information digital image | |
| US8442290B2 (en) | Simultaneous dual window/level settings for display of CT colonography images | |
| JP4314001B2 (ja) | 画像表示装置 | |
| US20030223628A1 (en) | Method and apparatus for generating a combined parameter map | |
| JP3378373B2 (ja) | 画像表示装置 | |
| Ivanov et al. | An adaptive system for quality improved digitalization of medical images | |
| CN101536034A (zh) | 用于放射图像的色调等级变换 | |
| DE102008040801B4 (de) | Verfahren zum Ableiten medizinischer Größen aus einer Bildfolge medizinischer Bilddaten, medizinische Vorrichtung und Analysesystern | |
| Ishigaki et al. | Verification of clinical test objects for diagnostic quality-guaranteed CRT monitor based on the visual evaluation |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20070425 Termination date: 20201121 |