CN1364400A - X射线成象***中自动确定自动亮度***采样窗口大小和位置的方法和装置 - Google Patents

X射线成象***中自动确定自动亮度***采样窗口大小和位置的方法和装置 Download PDF

Info

Publication number
CN1364400A
CN1364400A CN00808555A CN00808555A CN1364400A CN 1364400 A CN1364400 A CN 1364400A CN 00808555 A CN00808555 A CN 00808555A CN 00808555 A CN00808555 A CN 00808555A CN 1364400 A CN1364400 A CN 1364400A
Authority
CN
China
Prior art keywords
sampling window
abs
image
subregion
ray
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN00808555A
Other languages
English (en)
Inventor
维蒙·T·詹森
R·拉里·安德顿
巴里·K·哈诺维尔
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.)
OEC Medical Systems Inc
Original Assignee
OEC Medical Systems Inc
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
Family has litigation
First worldwide family litigation filed litigation Critical https://patents.darts-ip.com/?family=23187405&utm_source=***_patent&utm_medium=platform_link&utm_campaign=public_patent_search&patent=CN1364400(A) "Global patent litigation dataset” by Darts-ip is licensed under a Creative Commons Attribution 4.0 International License.
Application filed by OEC Medical Systems Inc filed Critical OEC Medical Systems Inc
Publication of CN1364400A publication Critical patent/CN1364400A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • H—ELECTRICITY
    • H05—ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
    • H05G—X-RAY TECHNIQUE
    • H05G1/00—X-ray apparatus involving X-ray tubes; Circuits therefor
    • H05G1/08—Electrical details
    • H05G1/26—Measuring, controlling or protecting
    • H05G1/30—Controlling
    • H05G1/36—Temperature of anode; Brightness of image power

Landscapes

  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Toxicology (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • X-Ray Techniques (AREA)
  • Image Analysis (AREA)

Abstract

用于在使用自动亮度***(ABS)采样窗口的闭环X射线成象***中提供自动亮度控制的方法。根据从与X射线***和被处理的图象有关的数据得到的包含空间灰度分布数据(24)的统计信息,调整ABS采样窗口的位置(28)、尺寸(30)和形状(32),由此根据被修改的ABS采样窗口的统计数据,使自动亮度控制能进行亮度和功率调节。

Description

X射线成象***中自动确定自动亮度***采样窗口 大小和位置的方法和装置
背景
1.发明领域
本发明一般涉及X射线成象装置。更具体地，本发明涉及在使用自动亮度控制采样窗口的闭环成象***中的自动亮度控制。
2.本领域现状
X射线成象***是众所周知的医疗诊断和***工具。在高压发生器向X射线管提供电源时产生X射线。一个电路通过加热该管的丝极而为该管的X射线辐照作准备。第二电路产生高压电势，用于在X射线管内加速从丝极(或阴极)向阳极发射的X射线。
X射线管内的丝极是盘绕的钨丝，当被电流加热时钨丝发射电子。这是低电压电路。加热丝极只需要相对较小的功率，并且丝极电流中小的变化导致X射线管电流大的变化。
从丝极发射的电子被聚焦在钨靶(阳极)中的点上。当电子因突然减速而在钨阳极内相互作用时产生X射线光子。调整钨靶表面的角度，使X射线向着X射线管输出窗口的方向反射。
为了理解高压发生器控制与其对X射线束穿透率及随后的诊断图象质量的影响之间的关系，关键在于记住X射线束的强度随着施加到管上的千伏电势(kVp)和管电流(mA)而变化。所产生的X射线光子的质量和强度几乎完全取决于kVp。
由于曝露在辐射下是有害的，因此所有实际使用的方法都降低辐照到患者的X射线。一种方法是用会部分或完全吸收X射线的材料使X射线束准直。正确的束准直通过禁止不衰减(不受阻碍)的X射线光子到达图象增强器也有助于视频成象***。
不同的临床程序具有各自的束准直要求。最普通的***使用固定的、可调整的叶片和可调整的光圈(iris)组合准直。***还可使用具有独立的电动控制和位置反馈的光圈和叶片准直仪。
在检查患者的过程中，当X射线光子穿过患者并接着通过使用图象增强器管或一些其它的X射线转换器件而被转变成可见光光子(lightphoton)时形成图象。X射线光子在撞击图象增强器的输入表面之前穿过各种质量和成分的组织和材料。X射线会穿透或会被任何位于X射线束路径上的东西吸收。当X射线撞到图象增强器的输入屏时，X射线光子被转变成可见光光子。在光阴极内，可见光光子被转变成电子并被静电透镜聚焦到输出磷光体上。输出磷光体再次把电子转变成可见光光子，在图象增强器管的输出窗口处图象是可见的。
在此阶段，在图象中计算物体的实际大小或尺寸需要以下知识：1)图象增强器的大小；2)图象增强器静电透镜的放大率；以及3)所述物体因其到图象增强器管表面的距离而引起的放大率。
当图象在图象增强器的输出上显示时，视频照相机捕捉此图象。自动亮度***(ABS)控制应用程序基于图象亮度统计动态地确定照相机增益、kVp和mA。在ABS采样窗口内的峰值亮度和区域平均亮度在设置ABS控制参数(kVp、mA和照相机增益)时用作常规因数。还应提到，图象增强器也可以是一些其它类型的图象接收器，如平板X射线图象接收器或一些其它的扫描图象X射线接收器。
当患者在检查过程中重新位于X射线束之下时，由于当X射线束穿过不同厚度和密度的身体组织和骨头时X射线束衰减发生变化，因此视频图象的亮度改变。为了补偿这些图象亮度的变化，已发明各种自动亮度补偿***。
例如，在美国专利4573183中，从视频信号图象的平均亮度得到自动亮度控制。此平均亮度信息用于产生控制照相机的视频增益信号并且产生控制X射线管电源的亮度反馈信号。为了响应亮度反馈信号，X射线管电源产生用于X射线管的偏压和丝极电流，由此调节图象增强管形成的X射线图象的亮度。通过改变视频照相机的增益系数或光圈，也可控制在监视器上得到的图象亮度。相应地，用于自动亮度控制的反馈改变影响图象亮度的X射线发射和照相机增益。此反馈控制可被校正，以产生与患者X射线吸收变化无关且一致的图象显示亮度。从视频增益和来自共有平均亮度值、峰值亮度或图象数值组合的亮度反馈信号得到标准的自动亮度控制。然而，共有亮度反馈信号对这些影响监视器上图象形成的***参数的控制不够理想。
本领域现状的另一实例在美国专利4703496中描述。此专利清楚地叙述，在每个视频图象区域中图象元素的辉度被取平均，以产生具有与平均图象亮度成正比的电压的信号。
平均亮度值用作反馈信号，以控制X射线管的激励和视频增益，由此维持最佳水平基本不变的视频图象亮度。然而，此***相对复杂，它利用三个单独的环路来调节管电流、偏压和视频增益。
上述***复杂且难以维护。相应地，在本领域现状中提供一种便于操作的自动亮度控制***是有益的，并且此***应能补偿患者不正确位置、患者体形差别以及躯干到肢端成象迅速变化的影响。
本发明目的及概述
本发明的一个目的是提供一种用于在闭环成象***中进行自动亮度***(ABS)控制的方法。
另一目的是提供一种用于ABS的方法，此ABS使用采样窗口确定ABS采样区域边界，此方法能补偿当最大的X射线衰减不对准图象中心的情况。
另一目的是提供一种用于使用ABS采样窗口的ABS的方法，此方法能补偿不衰减的X射线光子面积在成象区域内的情况。
另一目的是提供一种用于使用ABS采样窗口的ABS的方法，其中，作为图象内空间灰度分布函数的采样窗口位置是可调整的。
另一目的是提供一种用于使用ABS采样窗口的ABS的方法，其中，作为图象增强器尺寸、图象增强器放大模式以及图象kVp值的函数的采样窗口尺寸是可调整的。
另一目的是提供一种用于使用ABS采样窗口的ABS的方法，其中，作为准直函数的采样窗口的形状是可调整的。
还有一目的是提供一种用于ABS的方法，其中，在图象处理***中以动态函数执行此方法，由此尽管患者位置不正确、患者体形变化并且躯干到肢端成象迅速变化，也提供最佳的成象。
本发明的优选实施例是用于在使用ABS采样窗口的闭环成象***中提供自动亮度控制的方法。根据从感兴趣的图象区域内的数据得到的统计信息，调整ABS采样窗口的位置和几何形状。
在本发明的第一方面中，ABS采样窗口根据图象内空间灰度分布函数而移动。
在本发明的第二方面中，图象被划分成多个分区，其中，分析每个分区以确定对于ABS最有用的采样区域而不是确定图象的中心。
在本发明的第三方面中，考虑准直对图象的侵占。
在本发明的第四方面中，根据kVp、图象增强器的尺寸和放大模式调整ABS采样窗口的尺寸。
在本发明的第五方面中，通过在算法中增加瞬时取平均而减缓ABS采样窗口的移动。
在本发明的第六方面中，为灰度值设置阀值由此限制极值对算法的影响。
从以下结合附图而进行的详细描述中，本发明的这些和其它目的、特征、优点和替代方面对于本领域中技术人员是显而易见的。
附图简述
图1是一个图象的正视图，包括在圆形空白区域内建立的感兴趣区域(ROI)。
图2示出常规ABS采样窗口，此窗口直接对准诊断X射线图象的中心。
图3示出根据本发明的ABS采样窗口在已移动和确定尺寸后的状况。
图4示出被划分成两个横跨X坐标轴的分区的感兴趣区域。
图5示出被划分成两个横跨Y坐标轴的分区的感兴趣区域。
图6是图1图象的正视图，其中ROI现在被划分。
图7示出划分成象限的感兴趣区域，这些象限基于准直而被剪切。
图8是综述本发明优选实施例步骤的流程图。
图9示出在临床图象的模拟测试中***是如何体现的。
图10为示出代入公式1、2和3中的值与所得到的偏移X和Y的表格。
图11为X射线***的示意图，此***具有改变ABS采样窗口尺寸、位置和形状的调整单元。
本发明详述
现在参照附图，对本发明的各种元件分配数字编号，并且讨论本发明以使本领域中技术人员能实施和使用本发明。应理解，以下仅仅是本发明原理的示例性描述并且不应视为是对后附权利要求的限制。
在对优选实施例进行详细描述之前，对本发明作个综述是有帮助的。相应地观察到，本发明有利地提供一种利用ABS采样窗口的自动亮度***(ABS)。在图象中定义ABS采样窗口的两个特性为：1)窗口位置；和2)窗口几何形状。
通过比较感兴趣分区的图象区域的平均灰度值来确定ABS采样窗口的位置。基于kVp和图象增强器的尺寸和放大率调整ABS采样窗口的尺寸。进一步地，ABS采样窗口的形状也基于ABS采样窗口的准直来调整。
本发明的优选实施例提供用于闭环成象***中的ABS。本方法提供操作ABS采样窗口的动态ABS功能，由此提高成象***提供最佳成象的能力。通过使本***自动化，所产生的图象具有更好的诊断信息。
图1是视频图象10的正视图。感兴趣区域的尺寸(ROISize)在圆形空白区域12内建立。在此图中，示出的ROI 14为正方形区域，但如以后将解释地，ROI 14也可以是圆形或其它任意的几何形状。圆形空白区域12外部的区域16包括与将要进行的计算无关的信息。所有有用的信息在圆形空白区域12内示出。ROISize定义一些被检查的象素。
必须了解与图象有关的一些特性，以便理解在本发明中用于确定ABS窗口位置和尺寸的方法。首先，在图象中更暗的区域代表图象中最厚的或最密的组织。在此定义一区域作为大得足以代表明显的解剖学性质或特征的图象区域。
应该理解，在图象中最暗象素的位置对于ABS采样窗口而言不可能是最佳的中心点位置，因为最暗的象素可以是金属修补物或相似的材料。进而，所有无阻碍(非准直)的图象区域的灰度组成提供用于判断ABS采样窗口位置的有用信息。
图2示出临床X射线图象，其中，在腿38中被X射线照射的组织密度最大的区域(如髋关节34)不在图象30的中央。同时示出的ABS采样窗口32使用常规方法在目标图象30内使ABS采样窗口32居中。图2中的ABS采样窗口还包括具有不衰减的X射线光子的区域40。在平均或峰值灰度值计算中包括不衰减的X射线光子产生令人误解的输出。结果，增益或kVp不会增加到足够高的水平，因为ABS***会相信图象已足够亮了。
图3示出使用本发明方法进行定位的ABS采样窗口36。通过在图象30内在密实的患者组织34上使ABS采样窗口36居中，对图象的灰度范围进行优化。应该指出，根据本发明方法定位的ABS采样窗口36避免捕捉不衰减的X射线光子。重新确定ABS窗口的位置和尺寸也捕捉更密实些的组织区域，以确保采用正确的kVp、mA和增益在正确的动态范围内观察受验组织。换句话说，由于因重新定位而使不衰减的X射线光子面积最小化，因此组织图象不会太暗，并且在确定X射线强度时会考虑患者组织更密实的区域。
现在参照图4，在本发明中用于移动ABS采样窗口的方法比较ABS采样窗口14左半部60的平均灰度值和右半部62的平均灰度值，以确定水平(X轴)偏移。基于此描述还显而易见，左半部或右半部还可与整个图象12相比较，但比较ABS采样窗口的左右两半是更有效的。不考虑全部图象10中先前空白的外部区域16。
图5示出在ABS采样窗口14的上半部64和下半部66之间所作的相同类型的比较，以确定垂直(Y轴)偏移。
为便于描述和执行，本发明的方法计算图象中各象限的平均灰度值，然后结合适当的象限形成一个半区，此半区接着与所有象限平均灰度的总和相比较。在本质上，这与前述一样，比较一个半区和另一个半区。图6是图1中图象10的正视图，其中图象ROI 14被划分成多个象限。被划分的区域Q1、Q2、Q3和Q4定义成其统计信息被用作ABS采样窗口调整的各个区域。重要的是应注意，对于准直而言剪切得是否正确将影响分区形状。然而，分区(或象限)的数量、位置和形状根据具体的执行而改变。
尽管示出的分区Q1、Q2、Q3和Q4在图6中是邻接的，但分区不必是互相接触的。分区可以是互相分隔开，或者甚至是重叠的。比较图象四个象限的平均灰度值提供足够的信息，以对ABS采样窗口作相应的位置调整。如前面所提到的，每个象限不必具有相同的形状或面积。比较每个象限的灰度平均值使它们的统计权重相同，这与每个象限采样的象素比例无关。
在划分后，下一步是确定图象每个分区的平均灰度值。计算图象中任何可选择部分的平均灰度值的方法对于本领域技术人员是熟知的。例如，灰度值可在零(0)到255之间变化，零对应于完全黑的象素而255对应于完全白的象素。表示不衰减的X射线束的象素一般会被分配灰度值255，不衰减的X射线束是因没有患者组织、骨头或其它材料阻碍而得到的。相反，当X射线不能穿透物体时，象素一般被分配给值0。相应地，表示组织、骨头或其它已至少被部分X射线光子穿透的材料的象素被分配给在0-255数值之间的绝对值。本领域中技术人员应该清楚，如果需要可把灰度值扩展到12位、16位或更多。
在确定所有分区的平均灰度值之后，下一步是计算ABS采样窗口的X和Y轴向偏移量。这些偏移量被确定为在分区中的空间灰度分布的函数。
在图6中，ROI 14被划分成四个相同的分区。相应地，本优选方法所用的公式见下式(1)和(2)：公式1
  Xoff＝(ROISize/2)-(ROISize×(Q2％GS+Q4％GS))其中：Xoff＝ABS采样窗口的X偏移量
  ROISize＝ROI中横跨X轴的总象素数
  Q2％GS＝在第二分区中灰度分布的百分比
  Q4％GS＝在第四分区中灰度分布的百分比公式2
  Yoff＝(ROISize/2)-(ROISize×(Q3％GS+Q4％GS))其中：Yoff＝作用到ABS采样窗口的Y偏移量
ROISize＝ROI中横跨Y轴的总象素数
Q3％GS＝在第三分区中灰度分布的百分比
Q4％GS＝在第四分区中灰度分布的百分比
为计算公式1和2的结果，首先需要计算在分区中灰度分布的百分比。此结果用公式3决定。公式3
Q(n)％GS＝Q(n)AverageGS/SumOfQAverages其中：Q(n)％GS＝在分区n中灰度分布的百分比
Q(n)AverageGS＝在分区n中的平均灰度分布
SumOfQAverages＝在1到n的所有分区中平均灰度分布的总和
当已确定X偏移量和Y偏移量后，根据偏移量移动ABS采样窗口。基于此描述可看出，也可用Q1得到等效的公式。
在前述实施例中，应注意，图6中的分区Q1、Q2、Q3和Q4都在图象ROI 14的中点汇合。相应地，X偏移和Y偏移可以认为是从图象ROI14的中点或者从沿着图象ROI 14边缘的任意点作用的偏移。得到的结果是相同的。ABS采样窗口会正确地移动。
基于发生器kVp值、图象增强器的尺寸和放大模式确定ABS采样窗口的尺寸。在本发明中已研究kVp和组织结构成象之间的近似关系：公式4
kVp＝40+(2×以cm为单位的组织厚度)
根据在本发明中进行的观察，即大多数组织部分在横向和轴向上具有大致相同的横截面尺寸，确定此近似关系。此关联用于创建ABS采样窗口和kVp之间的关系，在此通过求解下面公式得到ABS采样窗口的尺寸：公式5
ABS窗口直径(cm)＝(kVp-40)÷2ABS采样窗口直径(如投影到图象接收器或图象增强器表面上的)也与图象增强器的尺寸和放大模式有关，因此图象增强器的模式必须标称化。这可通过在每种情况下把从以上公式计算得到的ABS采样窗口直径除以图象增强器视域直径而实现。
例如，假设kVp值是60kVp，对应于10cm厚的组织厚度并因此对应于组织的近似宽度或高度。另外，对投影到图象增强器上的采样窗口直径计算为10cm或最初为23cm视域图象直径的43％。那么，如果23cm图象增强器(或相似的图象接收器)以14cm放大模式操作，全屏幕的采样窗口百分比就是10/14或71％。因而ABS采样直径应为当前视域图象直径的71％。根据此描述应该认识到，在公式4和5中所用的常数可根据本发明的应用或ABS窗口和kVp之间所希望的的尺寸关系而修改。
由于kVp值用在ABS采样窗口中得到的数据修改，因此，直径范围必须有上限和下限以使ABS采样窗口不会变得不稳定。例如，如果在X射线***中观察较小的组织部分如患者的手指(如1cm厚)，那么由于尺寸的改变，ABS采样窗口有可能变得如此之小，以致于定位方法实际上会移动ABS窗口掠过手指，因为ABS窗口不知道手指在变暗的分区中。如果ABS采样窗口错过被观察的组织，它就只会捕捉不衰减的X射线光子，因此手指图象会变得更暗和不能使用。上限存在是因为一旦患者的组织超过图象尺寸的一定百分比就没有理由增加ABS采样窗口的尺寸，因为通过增加此尺寸得不到任何好处。例如，如果患者的躯干占满图象，那么增加ABS采样直径超过一定量时就不再会明显地影响kVp。
直径百分比的上限和下限可通过实验确定。优选的直径百分比是下限为图象直径的大约40％并且上限为约70％。此直径限制可根据图象增强器尺寸、图象增强器放大模式和具体的X射线应用而改变。
由于本发明方法的优选实施例比较每个分区内的平均灰度值和所有象限灰度值平均的总和，因此分区不需要具有相同的面积或一致的尺寸。当因准直仪侵占而必需修改ROI(感兴趣区域)时这是有好处的。
光圈和叶片准直仪提供X射线屏障以减少对患者的照射。当准直仪进入图象区域时它们产生黑暗的边界。如果可确定准直仪的位置，在采集平均灰度数据时，应把这些黑暗区域排除在外。图7说明因准直而引起地如何修改分区尺寸。
图8示出本优选实施例的步骤。第一步骤20选择图象ROI和相应的ABS采样窗口。下一步22把图象ROI划分成所需数量、尺寸和形状的分区。步骤24从公式3确定图象ROI每个分区的平均灰度值。步骤26根据在公式1、2和3确定的数值，找到ABS采样窗口的X和Y轴偏移量。步骤28根据在步骤26确定的偏移量移动ABS采样窗口。
步骤30基于发生器kVp值、图象增强器的尺寸和放大模式确定ABS采样窗口的大小。接着在步骤32中，修改采样窗口以计算准直仪对采样区域的侵占。最后在步骤34中，成象***在ABS采样窗口中进行必要的测量，由此确定用于调节图象亮度的ABS值。这些测量包括在ABS采样窗口内收集统计的图象象素数据。包括峰值和平均象素值的信息用作输入，以确定将在ABS控制逻辑模块内进行的调整。对ABS采样窗口位置、尺寸和形状确定的改进在由X射线***生成的组织图象内提供更好的灰度值动态范围。由本发明方法对ABS采样窗口所作的调整能从X射线图象获得更好的医疗和诊断信息。
图9提供在对临床图象的模拟测试中***如何操作的示例。具体地，临床图象40在计算机显示器的窗口内示出。应记住，从ABS采样窗口42得到的统计信息输入到用于成象***闭环控制的ABS控制***中。
在示范窗口中，在程序窗口左下角中有可选择的数值，这些数值输入以模拟一般由尺寸确定和象素采样方法提供的数值。这些数值包括尺寸44、偏移量46、粒度48、暗度限值(dark clip)50、直径52和偏移系数54。
尺寸44指采样窗口以象素为单位的尺寸。偏移量46指作用到ABS采样窗口的以象素为单位的偏移量。粒度48指象素采样水平。例如，如果象素被跳过并且在象限或分区中只有第2、4、8等象素被采样，那么采样就可以加速。暗度限值50指灰度阀值，这在以后解释。直径52指ABS采样窗口以象素为单位的直径，此直径由上述的kVp值、图象增强器尺寸和放大率控制。偏移系数54指减少ABS采样窗口从其中心原点移动的量。
上述的暗度限值50指阀值。对于某些灰度值(如矫形器具或不衰减的X射线束)有可能剧烈地偏离对分区的阅读。暗度限值50使***能忽略落在设定的灰度值范围之外的灰度值。在本发明优选实施例中，象素灰度值的有效范围可设定上界和下界。例如，范围可建立在50到200之间。此范围建立暗度限值和亮度阀值，这排除非常亮和非常暗的值。限值的实际范围可通过各种因素如自适应空间采样、黑暗掩蔽采样或直接实验来确定。
图10提供用于确定公式1、2和3结果的计算的结果页。给定的采样数量是对图9中临床图象40的分析结果。此结果表示获得X偏移量-42和Y偏移量-39。注意，ROISize值是以象素为单位的ROI正方形边的长度或ROI圆的直径。
四个分区Q1、Q2、Q3和Q4的百分比灰度值Q(n)包括相对较大范围的数值。这些数值是在组成图象ROI中整个图象的每个分区内的亮象素的百分比。注意，对于Q1平均灰度值是56.93，这意味着分区相当暗(在0-255的范围内)。
图11示出与X射线***一起使用的本发明。示出产生X射线104的X射线管100，X射线104穿过光圈和叶片准直仪102。X射线104接着穿过患者106并进入增强器管108(图象增强器)，增强器管108又产生可被视频照相机110检测到的图象。视频信号111成为图象处理器的输入，图象处理器又产生数字图象数据124，ABS采样窗口处理器114按照以上详细描述处理数字图象数据124。本发明的ABS采样窗口处理器114基于发生器kVp值130、图象增强器的尺寸和放大模式122确定ABS采样窗口的尺寸。
按照以上结合公式1-3所描述的，在114中计算ABS采样窗口的位置。最后，如果需要的话，基于从准直仪102得到的准直信息128剪切ABS采样窗口。应该提到，ABS采样窗口处理器可以是单一用途特定芯片(ASIC)、多用途ASIC或可编程通用处理器。
修改的ABS采样窗口(ROI定义)信息接着传递到图象数据分析仪116，图象数据分析仪116产生关于ABS采样窗口内图象的具体统计信息，此统计数据134传递到ABS控制118。随后，增益信号126基于接收到的图象统计从ABS控制118发送到照相机110。然后，基于接收到的统计信息134，ABS控制118发送关于kVp(千伏)138和mA(毫安)136的信号到高压(HV)发生器以控制***图象亮度。
本发明还有几个应该研究的替代实施例。本发明的一个替代实施例提供另一个确定ABS采样窗口尺寸的准则。整个图象平均灰度值与象限平均灰度值之比可用于确定ABS采样窗口的尺寸。因此，如果整个图象是不同寻常地亮，那么这意味着在图象内的组织较小并且采样窗口尺寸应减小。
在本发明的另一实施例中，确定ABS采样窗口尺寸的因数基于组织选择技术。如果ABS应用知道在图象区域中的组织例如为侧向的腕关节，那么狭窄的矩形采样窗口比圆形采样窗口提供更有意义的统计信息。
偏移量46(在图9中示出)也是替代实施例。具体地，本发明可配置成从图象中心移动ABS采样窗口的距离比确定的X偏移量和Y偏移量更小或更大。当本发明方法对于监视器上的图象不会自动地实现最佳ABS时，这有助于细微调整。
另一替代实施例提供一种减缓ABS采样窗口移动的方法。可通过为本方法增加瞬时取平均来获得此减缓效果。换句话说，为了提供动态(接近于实时)补偿，可对每帧或以高达每秒30次的速度进行取平均计算。
应该理解，上述设置仅仅是本发明原理应用的示例。只要不偏离本发明的精神和范围，本领域中技术人员可创造大量的变更和替代的设置。后附权利要求用于涵盖这些变更和设置。

Claims (27)

1.一种在X射线图象中确定自动亮度***(ABS)采样窗口的尺寸、位置和形状以及在闭环X射线成象***中调节图象亮度的方法，所述方法包括以下步骤：
(a)基于发生器kVp值的函数，确定ABS采样窗口的尺寸；
(b)基于在感兴趣的图象区域内的空间灰度分布，确定ABS采样窗口的位置；
(c)基于***中X射线的准直，确定ABS采样窗口的几何形状；以及
(d)根据从按步骤(a)-(c)中修改的ABS采样窗口得到的数据，调节图象的亮度控制，由此补偿患者不正确的位置、患者形体的变化以及躯干到肢端成象的迅速变化所带来的不利效果。
2.如权利要求1所述的方法，其中，确定自动亮度***(ABS)采样窗口尺寸的步骤进一步包括以下步骤：
使用来自高压发生器的千伏电势值，寻找ABS采样窗口以厘米为单位的直径；
所述直径除以以厘米为单位的图象增强器模式，得到ABS采样窗口相对于全图象直径的百分比；以及
设置ABS采样窗口直径为全图象直径的ABS采样窗口百分比。
3.如权利要求2所述的方法，其中，使用来自高压发生器的千伏电势值寻找自动亮度***(ABS)采样窗口以厘米为单位的直径的步骤进一步包括：从千伏电势值减去40并把结果除以2。
4.如权利要求1所述的方法，其中，基于发生器kVp值的函数确定自动亮度***(ABS)采样窗口尺寸的步骤进一步包括：基于发生器kVp值和图象接收器模式的函数，确定ABS采样窗口的尺寸。
5.如权利要求1所述的方法，其中，确定自动亮度***(ABS)采样窗口尺寸的步骤进一步包括：基于图象增强器的尺寸和图象增强器的放大模式确定增强器的模式。
6.如权利要求1所述的方法，其中，确定自动亮度***(ABS)采样窗口尺寸的步骤进一步包括：选择限制全图象视域直径最大值和最小值的ABS采样窗口尺寸。
7.如权利要求6所述的方法，其中，自动亮度***(ABS)采样窗口被限制在全图象视域直径的40％-70％。
8.如权利要求1所述的方法，其中，基于***中X射线的准直确定自动亮度***(ABS)采样窗口形状的步骤进一步包括：修改采样窗口以计算准直仪对ABS采样窗口的侵占。
9.如权利要求1所述的方法，其中，确定自动亮度***(ABS)采样窗口位置的步骤进一步包括以下步骤：
(1)在被处理的图象内选择感兴趣的区域；
(2)把感兴趣区域划分成多个分区；以及
(3)通过比较分区内相应的灰度值，基于图象内的灰度分布确定ABS采样窗口的位置。
10.如权利要求9所述的方法，其中，在显示的图象内选择感兴趣区域的步骤进一步包括在图象内生成圆形空白区域的步骤，其中感兴趣区域从圆形空白区域内选择。
11.如权利要求9所述的方法，其中，在图象内选择感兴趣区域的步骤进一步包括在感兴趣区域内建立自动亮度***(ABS)采样窗口的步骤，其中ABS采样窗口内的统计灰度分布数据用来确定对图象亮度的正确调节。
12.如权利要求9所述的方法，其中，把感兴趣区域划分成多个分区的步骤进一步包括选择自动亮度***(ABS)采样窗口以包含图象圆形区域的步骤。
13.如权利要求12所述的方法，其中，把感兴趣的图象区域(ROI)采样窗口划分成多个象限的步骤进一步包括把感兴趣的图象区域(ROI)划分成相等大小象限的步骤。
14.如权利要求9所述的方法，其中，通过比较分区内相应的灰度值，基于图象内的灰度分布确定ABS采样窗口位置的步骤进一步包括移动自动亮度***(ABS)采样窗口的步骤，以便ABS采样窗口内统计数据的测量会使***对亮度控制***进行调节，从而在监视器上获得改进的图象。
15.如权利要求14所述的方法，其中，移动自动亮度***(ABS)采样窗口的步骤进一步包括首先确定ABS采样窗口的X偏移量和Y偏移量的步骤。
16.如权利要求15所述的方法，其中，确定用于自动亮度***(ABS)采样窗口的X偏移量和Y偏移量的步骤进一步包括以下步骤：
(1)确定感兴趣的图象区域(ROI)中每个分区的平均灰度值；以及
(2)利用每个分区的平均灰度值和感兴趣区域的尺寸确定X偏移量和Y偏移量。
17.如权利要求16所述的方法，其中，确定感兴趣的采样窗口图象区域(ROI)中每个分区的平均灰度值的步骤进一步包括为每个分区计算以下公式的步骤：
     Q(n)％GS＝Q(n)AverageGS/SumOfQAverages
其中：Q(n)％GS＝在分区n中灰度分布的百分比
    Q(n)AverageGS＝在分区n中的平均灰度分布
    SumOfQAverages＝在1到n的所有分区中平均灰度分布的总和。
18.如权利要求16所述的方法，其中，利用四个分区中每一个的平均灰度值确定X偏移量和Y偏移量的步骤进一步包括利用下列公式确定X偏移量的步骤：
    Xoff＝(ROISize/2)-(ROISize×(Q2％GS+Q4％GS))
其中：Xoff＝自动亮度***采样窗口的X偏移量
    ROISize＝感兴趣区域中横跨X轴的总象素数
    Q2％GS＝在第二分区中灰度分布的百分比
    Q4％GS＝在第四分区中灰度分布的百分比。
19.如权利要求16所述的方法，其中，利用四个分区中每一个的平均灰度值确定X偏移量和Y偏移量的步骤进一步包括利用下列公式确定Y偏移量的步骤：
    Yoff＝(ROISize/2)-(ROISize×(Q3％GS+Q4％GS))
其中：Yoff＝作用到自动亮度***采样窗口的Y偏移量
    ROISize＝感兴趣区域中横跨Y轴的总象素数
    Q3％GS＝在第三分区中灰度分布的百分比
    Q4％GS＝在第四分区中灰度分布的百分比。
20.如权利要求16所述的方法，其中，此方法进一步包括后述步骤：根据被X射线辐照且在监视器上显示的患者组织厚度，使用整个感兴趣的自动亮度***图象区域(ROI)的平均灰度分布与分区灰度之比，由此调整ABS采样窗口的定位。
21.如权利要求15所述的方法，其中，此方法进一步包括通过利用瞬时取平均来减缓自动亮度***(ABS)采样窗口移动的步骤。
22.如权利要求9所述的方法，其中，通过比较分区内相应的灰度值，基于图象内的灰度分布确定ABS采样窗口位置的步骤进一步包括建立用于忽略灰度值的上限和下限由此避免使灰度分布值偏离的步骤。
23.如权利要求22所述的方法，其中，此方法进一步包括后述步骤：如果部分X射线辐照在所述上限和下限之外，从用于确定空间灰度分布的计算中删除一部分感兴趣的图象区域(ROI)采样窗口。
24.如权利要求9所述的方法，其中，把感兴趣的图象区域(ROI)采样窗口划分成多个分区的步骤进一步包括后述步骤：相对于感兴趣的图象区域(ROI)采样窗口原点使分区居中，以便多个分区互相不接触并与原点相等地分隔开。
25.一种在闭环X射线成象***中对在监视器上显示的图象的亮度进行自动调节的方法，所述方法包括以下步骤：
(1)在图象内选择感兴趣的区域；
(2)在感兴趣区域内建立自动亮度***(ABS)采样窗口；
(3)通过在ABS采样窗口内进行统计灰度分布的测量而调节ABS采样窗口的位置；以及
(4)根据ABS采样窗口内灰度统计调节图象亮度，由此补偿患者不正确的位置、患者形体的变化以及躯干到肢端成象的迅速变化所带来的不利效果。
26.如权利要求25所述的方法，其中调节自动亮度***(ABS)采样窗口位置的步骤进一步包括后述步骤：在感兴趣的图象区域(ROI)内基于与至少一个其它具有空间灰度分布的分区相比较的分区的空间灰度分布，调节自动(ABS)采样窗口的位置。
27.用于在X射线图象内确定自动亮度***(ABS)采样窗口的尺寸、位置和形状以及根据在被调整的ABS采样窗口内的灰度数据调节闭环X射线成象***中图象亮度的装置，其中包括：
耦合到X射线成象***的位置处理器，基于空间灰度分布，在从X射线成象***接收到的X射线图象中确定ABS采样窗口的位置，其中所述选择的分区在感兴趣区域中具有所需的统计灰度分布；
耦合到位置处理器的尺寸处理器，基于发生器kVp值的函数确定ABS采样窗口的尺寸；
耦合到尺寸处理器的形状处理器，基于X射线成象***中X射线准直确定ABS采样窗口的形状；
耦合到形状处理器的统计处理器，用于接收ABS采样窗口数据并基于ABS采样窗口数据产生统计信息；以及
耦合到统计处理器的自动亮度控制，根据从统计处理器收到的关于被修改的ABS采样窗口的统计信息调节图象亮度，由此补偿患者不正确的位置、患者形体的变化以及躯干到肢端成象的迅速变化所带来的不利效果。
CN00808555A 1999-05-07 2000-05-05 X射线成象***中自动确定自动亮度***采样窗口大小和位置的方法和装置 Pending CN1364400A (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US09/306,908 US6175614B1 (en) 1999-05-07 1999-05-07 Method and apparatus for automatic sizing and positioning of ABS sampling window in an x-ray imaging system
US09/306,908 1999-05-07

Related Child Applications (1)

Application Number Title Priority Date Filing Date
CN2008102133619A Division CN101365288B (zh) 1999-05-07 2000-05-05 确定自动亮度***采样窗口尺寸位置和形状的方法和装置

Publications (1)

Publication Number Publication Date
CN1364400A true CN1364400A (zh) 2002-08-14

Family

ID=23187405

Family Applications (2)

Application Number Title Priority Date Filing Date
CN2008102133619A Expired - Lifetime CN101365288B (zh) 1999-05-07 2000-05-05 确定自动亮度***采样窗口尺寸位置和形状的方法和装置
CN00808555A Pending CN1364400A (zh) 1999-05-07 2000-05-05 X射线成象***中自动确定自动亮度***采样窗口大小和位置的方法和装置

Family Applications Before (1)

Application Number Title Priority Date Filing Date
CN2008102133619A Expired - Lifetime CN101365288B (zh) 1999-05-07 2000-05-05 确定自动亮度***采样窗口尺寸位置和形状的方法和装置

Country Status (8)

Country Link
US (1) US6175614B1 (zh)
EP (1) EP1181846B1 (zh)
JP (1) JP4606593B2 (zh)
CN (2) CN101365288B (zh)
AT (1) ATE495651T1 (zh)
AU (1) AU4826800A (zh)
DE (1) DE60045509D1 (zh)
WO (1) WO2000069228A1 (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100405355C (zh) * 2003-05-30 2008-07-23 西门子公司 对数字x射线图像进行图像处理的自学习方法及其设备
CN102413769A (zh) * 2009-05-05 2012-04-11 皇家飞利浦电子股份有限公司 采集x射线图像的方法和包括自动楔定位的x射线图像采集装置
CN101721219B (zh) * 2008-10-17 2012-08-22 佳能株式会社 控制装置及控制方法
CN102113889B (zh) * 2009-12-30 2013-01-23 航天恒星空间技术应用有限公司 自动亮度透视控制***控制信号的双窗口拾取方法
CN103177183A (zh) * 2011-12-21 2013-06-26 三星电子株式会社 医学设备及使用所述医学设备的图像显示方法
CN103211607A (zh) * 2013-04-28 2013-07-24 北京东方惠尔图像技术有限公司 X 光机及其控制方法和装置
US10383594B2 (en) 2015-12-01 2019-08-20 Shenyang Neusoft Medical Systems Co., Ltd. Image brightness adjustment
CN111789609A (zh) * 2020-08-10 2020-10-20 上海联影医疗科技有限公司 一种成像***和方法

Families Citing this family (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1053659B1 (en) * 1998-12-08 2008-09-24 Koninklijke Philips Electronics N.V. An x-ray examination apparatus having an object absorption dependent brightness control
CN1331022C (zh) * 2000-03-31 2007-08-08 皇家菲利浦电子有限公司 用于获取辐射图像的***和方法
US6327336B1 (en) * 2000-06-05 2001-12-04 Direct Radiography Corp. Radiogram showing location of automatic exposure control sensor
US7062714B1 (en) * 2000-07-28 2006-06-13 Ge Medical Systems Global Technology Company, Llc Imaging system having preset processing parameters adapted to user preferences
JP2002072998A (ja) * 2000-08-25 2002-03-12 Internatl Business Mach Corp <Ibm> 輝度制御装置、輝度調節システム、コンピュータシステム、液晶ディスプレイ装置、輝度制御方法、コンピュータプログラム及び記憶媒体
US6501819B2 (en) * 2000-12-18 2002-12-31 Ge Medical Systems Global Technology Company, Llc Medical diagnostic method and apparatus to control dual energy exposure techniques based on image information
DE60228555D1 (de) * 2001-03-29 2008-10-09 Toshiba Kk Röntgendiagnosevorrichtung
DE10132816A1 (de) 2001-05-31 2002-12-05 Philips Corp Intellectual Pty Vorrichtung und Verfahren zur Anpassung der Strahlungsdosis einer Röntgenstrahlungsquelle
US6614877B2 (en) * 2001-11-21 2003-09-02 Ge Medical Systems Global Technology Company Llc Method and apparatus for enhancing the contrast of a medical diagnostic image acquired using collimation
DE10163215B4 (de) * 2001-12-21 2020-02-20 Philips Gmbh System und Verfahren mit automatisch optimierter Bilderzeugung
US6795526B2 (en) * 2002-03-04 2004-09-21 Ge Medical Systems Global Technology Co., Llc Automatic exposure control for a digital image acquisition system
US7180981B2 (en) * 2002-04-08 2007-02-20 Nanodynamics-88, Inc. High quantum energy efficiency X-ray tube and targets
US20040196958A1 (en) * 2002-11-29 2004-10-07 Werner Beck Operating device for a diagnostic imaging unit
EP1651112B1 (en) * 2003-07-30 2010-10-27 Philips Intellectual Property & Standards GmbH X-ray unit having an automatically adjustable collimator
US7734119B2 (en) * 2004-09-21 2010-06-08 General Electric Company Method and system for progressive multi-resolution three-dimensional image reconstruction using region of interest information
JP4325552B2 (ja) * 2004-12-24 2009-09-02 セイコーエプソン株式会社 画像処理装置、画像処理方法及び画像処理プログラム
DE102005052993B4 (de) * 2005-11-07 2014-08-21 Siemens Aktiengesellschaft Verfahren zur automatisierten Auswertung eines dreidimensionalen Abbildes eines seitensymmetrischen Organsystems
JP5032082B2 (ja) * 2006-09-29 2012-09-26 富士フイルム株式会社 放射線画像撮像装置及びその制御方法、並びに、放射線量検出器調整装置及びその調整方法
JP5464799B2 (ja) * 2007-11-16 2014-04-09 キヤノン株式会社 画像処理装置、画像処理方法及びプログラム
US8228570B2 (en) * 2008-01-28 2012-07-24 Air Techniques, Inc. Process and apparatus for scanning imaged storage plates and having automatic gain adjustment
JP5305747B2 (ja) * 2008-06-17 2013-10-02 キヤノン株式会社 放射線画像撮影装置及びその駆動方法、並びに、プログラム
JP5491836B2 (ja) * 2009-01-30 2014-05-14 株式会社東芝 超音波診断装置、超音波画像処理装置、医用画像診断装置及び医用画像処理装置
US8694075B2 (en) * 2009-12-21 2014-04-08 General Electric Company Intra-operative registration for navigated surgical procedures
DE102012216269A1 (de) * 2012-09-13 2014-03-13 Siemens Aktiengesellschaft Röntgensystem und Verfahren zur Erzeugung von Bilddaten
US9594033B2 (en) * 2014-07-22 2017-03-14 The Boeing Company Visible X-ray indication and detection system for X-ray backscatter applications
JP6441015B2 (ja) * 2014-10-06 2018-12-19 キヤノンメディカルシステムズ株式会社 Ｘ線診断装置及びｘ線管制御方法
CN108924432B (zh) * 2018-06-12 2021-01-05 东软医疗***股份有限公司 一种信号偏差的校正方法、装置及设备
US11490872B2 (en) 2020-08-21 2022-11-08 GE Precision Healthcare LLC C-arm imaging system and method
CN113347369B (zh) * 2021-06-01 2022-08-19 中国科学院光电技术研究所 一种深空探测相机曝光调节方法、调节***及其调节装置

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5973898A (ja) * 1982-10-19 1984-04-26 Hitachi Medical Corp Ｘ線自動露出装置
US4573183A (en) 1984-06-29 1986-02-25 General Electric Company X-Ray image brightness control
NL8502569A (nl) * 1985-09-20 1987-04-16 Philips Nv Roentgenonderzoekapparaat met een locaal opgedeelde hulpdetector.
DE3732634A1 (de) * 1987-09-28 1989-04-06 Siemens Ag Roentgendiagnostikeinrichtung
EP0362427B1 (de) * 1988-10-05 1993-01-27 Siemens Aktiengesellschaft Röntgendiagnostikanlage mit einem Detektor für die mittlere Bildhelligkeit
EP0383963B1 (de) * 1989-02-20 1993-09-08 Siemens Aktiengesellschaft Röntgendiagnostikeinrichtung
US5003572A (en) 1990-04-06 1991-03-26 General Electric Company Automatic brightness compensation for x-ray imaging systems
US5461658A (en) * 1993-05-21 1995-10-24 U.S. Philips Corporation X-ray examination apparatus
DE4330787A1 (de) * 1993-09-10 1995-03-23 Siemens Ag Verfahren zum Betrieb eines Röntgenbelichtungsautomaten
JPH08102396A (ja) * 1994-09-29 1996-04-16 Shimadzu Corp Ｘ線透視撮影装置
US5574764A (en) * 1995-06-06 1996-11-12 General Electric Company Digital brightness detector
US5675624A (en) 1995-12-14 1997-10-07 General Electric Company Adaptive x-ray brightness and display control for a medical imaging system
DE69815252T2 (de) * 1997-04-24 2004-04-29 Koninklijke Philips Electronics N.V. Belichtungssteuerung auf basis von einem bedeutenden teil eines röntgenstrahlbildes
EP0909527B1 (en) * 1997-04-24 2008-12-31 Koninklijke Philips Electronics N.V. X-ray examination apparatus including an exposure control system

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100405355C (zh) * 2003-05-30 2008-07-23 西门子公司 对数字x射线图像进行图像处理的自学习方法及其设备
CN101721219B (zh) * 2008-10-17 2012-08-22 佳能株式会社 控制装置及控制方法
CN102413769A (zh) * 2009-05-05 2012-04-11 皇家飞利浦电子股份有限公司 采集x射线图像的方法和包括自动楔定位的x射线图像采集装置
CN102413769B (zh) * 2009-05-05 2014-11-26 皇家飞利浦电子股份有限公司 采集x射线图像的方法和包括自动楔定位的x射线图像采集装置
CN102113889B (zh) * 2009-12-30 2013-01-23 航天恒星空间技术应用有限公司 自动亮度透视控制***控制信号的双窗口拾取方法
CN103177183A (zh) * 2011-12-21 2013-06-26 三星电子株式会社 医学设备及使用所述医学设备的图像显示方法
CN103211607A (zh) * 2013-04-28 2013-07-24 北京东方惠尔图像技术有限公司 X 光机及其控制方法和装置
CN103211607B (zh) * 2013-04-28 2015-06-03 北京东方惠尔图像技术有限公司 X光机及其控制方法和装置
US10383594B2 (en) 2015-12-01 2019-08-20 Shenyang Neusoft Medical Systems Co., Ltd. Image brightness adjustment
CN111789609A (zh) * 2020-08-10 2020-10-20 上海联影医疗科技有限公司 一种成像***和方法
WO2022032455A1 (en) * 2020-08-10 2022-02-17 Shanghai United Imaging Healthcare Co., Ltd. Imaging systems and methods
CN111789609B (zh) * 2020-08-10 2023-12-08 上海联影医疗科技股份有限公司 一种成像***和方法

Also Published As

Publication number Publication date
EP1181846B1 (en) 2011-01-12
AU4826800A (en) 2000-11-21
EP1181846A4 (en) 2009-05-06
US6175614B1 (en) 2001-01-16
WO2000069228A1 (en) 2000-11-16
ATE495651T1 (de) 2011-01-15
CN101365288B (zh) 2011-12-07
DE60045509D1 (de) 2011-02-24
EP1181846A1 (en) 2002-02-27
JP4606593B2 (ja) 2011-01-05
CN101365288A (zh) 2009-02-11
JP2003522370A (ja) 2003-07-22

Similar Documents

Publication Publication Date Title
CN1364400A (zh) X射线成象***中自动确定自动亮度***采样窗口大小和位置的方法和装置
EP1443857B1 (en) X-ray image enhancement
CN103854941B (zh) 在带电粒子显微镜中执行样本的层析成像的方法
WO2010113415A1 (ja) Ｘ線撮影装置及びその制御方法、コンピュータプログラム
US10132761B2 (en) Method of constructing 3D image, image processor, and electron microscope
US20110174972A1 (en) Apparatus and methods for controlling electron microscope stages
WO2011093058A1 (en) Radiation imaging apparatus, and control method and program of the apparatus
DE102011054114A1 (de) Verfahren und System für den Betrieb eines Elektronenstrahlsystems
US9782138B2 (en) Imaging apparatus and method
CZ2017837A3 (cs) Zařízení pro vyzařování svazků paprsků s nabitými částicemi a způsob ovládání tohoto zařízení
CN107847198A (zh) 辐射图像处理方法及射线照相***
CN104323787B (zh) X射线断层扫描方法和***
US6437330B1 (en) Method and apparatus for adjusting a charged particle beam of a beam optical system
JP2005207827A (ja) Ｘ線検査装置、ｘ線検査方法およびｘ線検査装置の制御プログラム
JPH11296680A (ja) コントラスト改善機能を有する画像処理装置及び荷電粒子線装置
US11676796B2 (en) Charged particle beam device
JP2020027101A (ja) Ｘ線撮像装置
JP2011098009A (ja) Ｘ線診断装置
EP2383767A1 (en) Method of imaging an object
US20230115486A1 (en) Charged Particle Beam System and Control Method Therefor
US6839405B2 (en) System and method for electronic shaping of X-ray beams
CN114531767B (zh) 一种手持式x光机可视化x射线定位方法及***
JP6218365B2 (ja) 医用画像診断装置
KR101497105B1 (ko) 전자 현미경 및 전자 현미경을 이용한 영상 처리 방법
US10114209B2 (en) Microscope apparatus

Legal Events

Date Code Title Description
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C06 Publication
PB01 Publication
C12 Rejection of a patent application after its publication
RJ01 Rejection of invention patent application after publication