CN104809460B - 晶体中心位置图生成方法及晶体像素查找表生成方法 - Google Patents

晶体中心位置图生成方法及晶体像素查找表生成方法 Download PDF

Info

Publication number
CN104809460B
CN104809460B CN201510214482.5A CN201510214482A CN104809460B CN 104809460 B CN104809460 B CN 104809460B CN 201510214482 A CN201510214482 A CN 201510214482A CN 104809460 B CN104809460 B CN 104809460B
Authority
CN
China
Prior art keywords
connected domain
generating method
event
photo
germ nucleus
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510214482.5A
Other languages
English (en)
Other versions
CN104809460A (zh
Inventor
周鑫
李丽平
李强
吕新宇
安少辉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai United Imaging Healthcare Co Ltd
Original Assignee
Shanghai United Imaging Healthcare Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai United Imaging Healthcare Co Ltd filed Critical Shanghai United Imaging Healthcare Co Ltd
Priority to CN201510214482.5A priority Critical patent/CN104809460B/zh
Priority to CN201611266541.4A priority patent/CN106611173B/zh
Publication of CN104809460A publication Critical patent/CN104809460A/zh
Priority to US15/142,295 priority patent/US9928437B2/en
Application granted granted Critical
Publication of CN104809460B publication Critical patent/CN104809460B/zh
Priority to US15/897,110 priority patent/US10176393B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16—Measuring radiation intensity
    • G01T1/161—Applications in the field of nuclear medicine, e.g. in vivo counting
    • G01T1/164—Scintigraphy
    • G01T1/1641—Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras
    • G01T1/1647—Processing of scintigraphic data
    • 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/24—Aligning, centring, orientation detection or correction of the image
    • 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/24—Aligning, centring, orientation detection or correction of the image
    • G06V10/245—Aligning, centring, orientation detection or correction of the image by locating a pattern; Special marks for positioning
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16—Measuring radiation intensity
    • G01T1/161—Applications in the field of nuclear medicine, e.g. in vivo counting
    • G01T1/164—Scintigraphy
    • G01T1/1641—Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras
    • G01T1/1644—Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras using an array of optically separate scintillation elements permitting direct location of scintillations
    • 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/28—Quantising the image, e.g. histogram thresholding for discrimination between background and foreground patterns
    • 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/40—Extraction of image or video features
    • 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/40—Extraction of image or video features
    • G06V10/46—Descriptors for shape, contour or point-related descriptors, e.g. scale invariant feature transform [SIFT] or bags of words [BoW]; Salient regional features
    • G06V10/462—Salient features, e.g. scale invariant feature transforms [SIFT]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Multimedia (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Image Analysis (AREA)
  • Measurement Of Radiation (AREA)
  • Image Processing (AREA)

Abstract

本发明提供的晶体中心位置图生成方法，采用连续非极大值衰退的方法进行晶体中心位置定位，本方法本质上是对图像中的非波峰部分进行连续衰减，从而突出波峰部分，对于一维或者二维寻峰算法而言，极大地提高了寻峰的准确度。同时，本发明还提供了晶体像素查找表生成方法。

Description

晶体中心位置图生成方法及晶体像素查找表生成方法
技术领域
本发明涉及正电子发射断层成像领域，尤其是涉及一种正电子发射断层成像系统探测器晶体中心位置图生成方法及晶体像素查找表生成方法。
背景技术
正电子湮灭事例经PET探测器像素阵列解码之后产生事例解码之后产生事例的位置分布，统计之后形成光子事件二维直方图，由于探测器实现中编码与解码的非线性，生成的光子事件二维直方图会形成不规则形变。以一个探测器block为例，晶体阵列为16*16的等距规则排列像素阵列，解码之后会发生桶形、蝶形、旋转、压缩、扩张等各种形式的形变。若不做任何处理，数据获取时，在线处理无法判断晶体像素的行列坐标，直接导致图像分辨率下降，严重的可以导致图像信息错误。如果能确定之后的行列坐标和解码前的硬件晶体之间的对应关系，就能判断晶体像素的行列坐标。目前业界广泛使用的做法是使用晶体像素查找表(Crystal Lookup Table，CLT)来描述解码前后的晶体对应关系。
由于PET探测器的晶体漂移随时间推移会不停恶化，使CLT每隔一段时间就要被重新绘制。现阶段CLT靠手工绘制，但是一台PET设备往往有几万个晶体，手工绘制过程耗时耗力，成为PET设备维护中的主要时间瓶颈。为了既保证矫正效果，又缩短维护时间，目前最常用的方法就是半自动绘制CLT，即计算机自动生成一个中间结果：往往不是分割好的CLT，而是晶体中心位置图CCPM(Crystal Central Position Map)。自动生成的CCPM质量越好，技术检查和修正结果花费的时间越少。
现有技术中，已经公开的自动定位各晶体中心位置的算法，大致可以分为基于一维投影和直接在二维图像域处理两大类。基于一维投影的方法先把光子事件二维直方图投影到两个方向，最常见的是将图像分别投影在X方向和Y方向分别投影，则二维检测问题被简化为一维寻峰问题。求解可以使用一维的极值检测，一维求导定位法，甚至更加复杂的Gaussian fitting。但是投影算法过度依赖光子事件二维直方图中晶体行列的规整性，当晶体行列出现扭曲时，一维寻峰会失败。二维图像域处理包括特征检测域聚类算法两大类。特征检测的实际性能往往取决于图像质量和预处理，容易受图像噪声、灰度、不均匀性等影响出现漏检和错检。聚类算法缺少空间约束，更容易受伪影的影响，即使独立的伪影也能使算法结果出错，极大地影响探测结果，并导致后续的晶体映射错位，严重的会导致整个晶体阵列的映射完全变形。综上所述，现有众多的晶体中心位置检测算法，稳定性差，出现误检和漏检的可能性非常大。因此需要提出一种稳定性较好的晶体中心位置检测方法。
发明内容
为了解决现有技术中晶体中心位置检测方法错检率高，性能不稳定的问题，本发明提供了一种晶体中心位置图生成方法。
一种晶体中心位置图生成方法，包括如下步骤：
获取光子事件二维直方图I，将所述光子事件二维直方图I归一化至[0，1]之间，获得图像I0；对所述图像I0进行迭代处理得到图像Ii；使用自动阈值对所述图像Ii进行二值化处理，获得二值图像集SB；提取所述二值图像集SB中的每幅二值图像的连通域及其对应的中心位置；判断所述连通域彼此是否有重叠；保留有重叠的连通域中迭代次数最大的连通域；生成晶体中心位置图，所述晶体中心位置图包括有重叠的连通域中迭代次数最大的连通域对应的中心位置以及没有重叠的连通域对应的中心位置。
优选地，所述判断所述连通域彼此是否有重叠的步骤具体是：
将一个连通域内的任意位置与另一个连通域内的任意位置进行比较，如果相同即表示这两个连通域彼此有重叠。
优选地，所述保留有重叠的连通域中迭代次数最大的连通域的步骤具体为：比较所述有重叠的连通域的对应图像Ii的迭代次数i，剔除对应迭代次数较小的连通域。
优选地，所述有重叠的连通域中仅保留迭代次数最大的连通域的步骤具体为：比较所述彼此有重叠的连通域的面积，剔除面积较大的连通域。
优选地，还包括设置面积阈值，剔除连通域面积大于所述面积阈值的连通域。
优选地，所述将光子事件二维直方图I归一化至[0，1]之间采用线性收放的方式进行，即max(I)是所述光子事件二维直方图中的最大值，min(I)是所述光子事件二维直方图中的最小值。
优选地，所述光子事件二维直方图最大值max(I)按照比例ε获得，即取所述光子事件二维直方图的灰度直方图面积的p处的灰度为最大值max(I)，p＝1-ε，图中超出最大值的max(I)的值强制转换成max(I)。
优选地，所述比例其中n是所述光子事件二维直方图中实际包含的晶体数量，Np是光子事件二维直方图的像素数量。
优选地，所述自动阈值的选取基于光子事件二维直方图的灰度直方图，所述阈值将所述光子事件二维直方图的灰度直方图分成两部分，所述两部分的质心到所述自动阈值的距离相等。
本发明还提供了一种晶体像素查找表生成方法，采用上述的方法生成的晶体中心位置图，根据所述晶体中心位置图划分各相邻晶体条的分界线，生成晶体像素查找表。
与现有技术相比，本发明提供的晶体中心位置图生成方法，采用连续非极大值衰退的方法进行晶体中心位置定位，本方法本质上是对图像中的非波峰部分进行连续衰减，从而突出波峰部分，对于一维或者二维寻峰算法而言，极大地提高了寻峰的准确度。同时，本发明还提供了晶体像素查找表生成方法。
附图说明
图1为本发明晶体中心位置图生成方法的流程示意图；
图2为本发明的连续非极大值衰退的示意图；
图3为本发明生成的晶体像素查找表。
具体实施方式
在下面的描述中阐述了很多具体细节以便于充分理解本发明。但是本发明能够以很多不同于在此描述的其它方式来实施，本领域技术人员可以在不违背本发明内涵的情况下做类似推广，因此本发明不受下面公开的具体实施的限制。其次，本发明利用示意图进行详细描述，在详述本发明实施例时，为便于说明，所述示意图只是实例，其在此不应限制本发明保护的范围。
如背景技术所述，现有技术中的晶体中心位置检测方法，不管是一维投影法，还是直接在二维图像域处理法，其稳定性和精确性都有待提高。因此本发明提出一种稳定性较好的晶体中心位置检测方法。
如图1所示，一种晶体中心位置图生成方法，包括如下步骤：
步骤S10，获取光子事件二维直方图I，将所述光子事件二维直方图I归一化至0到1之间，获得图像I0。
对入射光子事件位置进行二维分布统计，获取光子事件二维直方图I，将所述光子事件二维直方图进行归一化处理至[0，1]之间。在一个实施例中，所述将光子事件二维直方图归一化至[0，1]之间采用线性收放的方式进行，即max(I)是所述光子事件二维直方图中的最大值，min(I)是所述光子事件二维直方图中的最小值。为了避免极端值对归一化的影响，对图像面积的ε处做截断处理，所述比例ε取决于晶体阵列中晶体的数量n和光子事件二维直方图的像素数Np。在本发明中，采用的是16*16的晶体阵列，则每个阵列有256个晶体，而光子事件二维直方图的分辨率(即像素数Np)是256*256。因此，ε＝256/(256*256)＝0.4％，那么图像直方图面积的99.6％处的值为max(I)，图像中大于max(I)的值都强制转成max(I)，即全图最大值就是max(I)。
在一个实施例中，在对图像进行处理归一化处理之前，对所述光子事件位置直方图进行均衡化矫正。本技术方案采用CLAHE算法对图像进行均衡化矫正，但是其他矫正的方法均可使用，例如用原图减去或除以图像中缓慢变化的背景内容。图像中缓慢变化的背景内容可以通过在图像域中使用大尺度平滑，在频域中压制高频部分，或使用曲面拟合等方法获得。
步骤S20，对所述图像I0进行迭代处理得到图像Ii。
所述迭代处理即为连续非极大值衰退。所述图像I0经过多次迭代处理，图像中的非波峰部分进了连续的衰减，突显出波峰部分，突显的波峰部分更易于被探测。
迭代之后的图像Ii＝Ii-1×I0，其中i∈{1…n}，i为迭代的次数，n为最大迭代次数，即终止条件。在一个实施例中，对图像I0进行50次迭代处理，如图2所示，分别为迭代1次、5次、10次、15次、20次、25次、30次、35次、40次的中间图像。从图中可以看出，经历迭代的次数越多，非波峰部分信号越弱，波峰部分越明显。
步骤S30，使用自动阈值对所述图像Ii进行二值化处理，获得二值图像Bi。
自动阈值算法有许多种，比如著名的otsu分割，最大熵分割，isodata等等。在一个实施例中采用isodata自动阈值法，但是其他自动阈值方法均可使用。具体实现方式如下：
先提取图像Ii中的非零最小值为阈值t的，开始迭代：
a)计算图像Ii中灰度大于t的所有像素的平均值msup；
b)计算图像Ii中灰度小于t的所有像素的平均值minf；
c)检验t是否等于msup和minf的平均值，如是，则结束，如否，则增加t的值(t＝t+1)，从回步骤a)。
迭代结束后t的值，正好是图像Ii的灰度直方图中大于t部分的质心和小于t部分的质心的平均值，即所述两部分的质心到所述自动阈值的距离相等。此时的t即为最后的阈值。以t对图像Ii进行二值化，大于t的位置赋为1，小于t的位置赋为0，得到图像Bi。
步骤S40，对所述二值图像集SB中的每幅二值图像Bi提取连通域及其对应的中心位置。
在一个实施例中，通过如下方式提取连通域及其对应的中心位置：建立一个空的堆栈T，遍历Bi上的像素位置(x，y)，当其值为1时，进行如下步骤：
1)将当前位置设为(x，y)；2)将(x，y)压入堆栈T；3)遍历Bi上(x，y)的四邻域((x+1，y)，(x-1，y)，(x，y+1)，(x，y-1))，如果其中的某个位置值为1，针对此位置，重复步骤1)-3)；4)将T中所有位置取出，并将Bi上这些位置对应的值改为0。取出的所有位置即为一个连通域，位置的个数即为此连通域面积，所有位置的平均值既为此连通域的中心。
步骤S50，判断所述连通域彼此是否有重叠。
将一个连通域内的任意位置与另一个连通域内的任意位置进行比较，如果相同即表示这两个连通域彼此有重叠。
在一个实施例中，如果一个连通域的中心位置与另一个连通域的任意位置相同，即表示这两个连通域彼此有重叠。
在一个实施例中，将新提取的连通域内的所有位置与已确定保留的所有连通域的中心位置进行一一比较，如果相同，即表示这新提取的连通域与已确定保留的某一个连通域有重叠。
步骤S60，保留有重叠的连通域中迭代次数最大的连通域。
由于最终的晶体中心位置图是由连通域集组成，重叠的连通域实质上是对同一个晶体重复记录，因此需要有选择性地过滤有重叠的连通域。连通域来源于二值图像Bi。来源于相同二值图像的连通域彼此之间是没有重叠的，但不同二值图像所产生的连通域会出现重叠。故相同迭代次数的连通域彼此没有重叠，彼此有重叠的连通域迭代次数必不相同。可以直接通过比较重叠连通域的迭代次数，保留迭代次数多的连通域来过滤连通域中的重叠。也可以通过比较连通域面积，保留面积较小的连通域达到同样的目的，因为迭代次数多的连通域的面积必然小于迭代次数少的连通域的面积。
实际应用中，为了减少重复运算，可以将步骤S40，S50和S60同时进行(即在提取连通域时就决定此连通域是保留还是丢弃)。同时可以通过设置连通域面积阈值c进一步加快运行速度。具体的，需要如下步骤：
1)建立确定保留的连通域集合Sc以及其对应的中心位置的集合Sp，初始使Sc和Sp为空集；
2)二值图像Bi＝n得到的连通域彼此之间不会有重叠。且Bi＝n是所有二值图像中迭代次数最大的，因此如果与其他二值图像得到的连通域有重叠，都应该保留二值图像Bi＝n得到的连通域。因此将所有二值图像Bi＝n得到的连通域加入Sc中，其中心位置加入Sp；
3)按从i＝n-1，i＝n-2，…直到i＝1的顺序从二值图像Bi中提取连通域。提取每个连通域时，如果其面积小于面积阈值c，且此连通域不包括任何Sp已有的中心位置，则将此连通域加入Sc中，其中心位置加入Sp，否则丢弃。得到的Sc即为筛选后的连通域集，Sp即为晶体中心位置集。
步骤S70，生成晶体中心位置图，所述晶体中心位置图包括有重叠的连通域中迭代次数最大的连通域对应的中心位置以及没有重叠的连通域对应的中心位置。在一个实施例中，将晶体中心位置集Sp中所有位置在二维图像上标出，即生成晶体中心位置图。
本发明还提供了晶体像素查找表生成方法，采用上述生成的晶体中心位置图，根据所述晶体中心位置图划分各相邻晶体条的分界线，生成晶体像素查找表。
根据晶体中心位置图划分各相邻晶体条的分界线，利用得到的分界线完成对晶体条区域的划分填充，将填充数据回朔到图像存储中，完成晶体像素查找表的生成。
对于实际的晶体越是到中间能量越大，即像素值越大；而越到周边，其能量越小。因此，分界线所在位置像素值是全局最小，以及分界线本身光滑连续的这两个条件，这两个条件互相制约，形成了传统的内能和外能。本发明使用更精准的约束条件：由于晶体实际上是以行列形式排列的，因此一条连续的分界线应该是相邻两行或两列晶体之间的分界线。由于在建立网格模板时使用样条把晶体中心点的行列都串了起来，这等于为分界线确立了范围。而且，每条分界线从模型上来说就是一条最优化路径。
本实施例使用动态规划确定各晶体条的分界线，对于一条相邻两行之间的分界线可以分为w段，其中w为晶体中心图像的宽度。动态规划中所需要的能量函数为En(x,x′)＝a*d(x,x′)+b*I(x)。其中，En(x,x′)为相邻两段之间的能量递增,x和x'分别是相邻两段的位置，d为取距离差函数。由于我们要求分界线尽可能平滑，因此相邻两段之间的距离差被限制在2个像素以下。I(x)为0-1均衡化后的图像像素值，a和b都是能量函数的权重，其数值大小根据实际需要可以调节，在此具体实施例中取a＝0.1，b＝1.5，由于图像中心信息比较准确，我们采用从中间向左右两边的顺序寻找最短路径。
具体地包括：(1)在正确的晶体中心位置图上以两条串起相邻两行晶体中心点的样条为边界线并将上述边界线分为w段，令边界线所围的图像中心位置为阶段0，记阶段0以左为阶段1至阶段e，记阶段0以右为阶段e+1至阶段w-1，记Ei(x)为阶段i的能量，1≤i≤w-1，令E0(x)＝1.5I(x)。(2)获取阶段1至阶段e所有阶段累积能量最小值，并记录其回朔位置x1’，其中每一阶段的累积能量Ei(x)＝Ei-1(x′)+En(x,x′)。同样阶段e+1至阶段w-1所有阶段累积能量最小值，并记录其回朔位置x2’。(3)根据起始位置进行左右路径连接的所有组合，计算总累积能量。(4)总累积能量最小的路径即为分界线。通过上述操作即可完成晶体分割，生成如图3所示的最终的晶体像素查找表。
本发明涉及的晶体中心位置图生成方法，采用连续非极大值衰退的方法进行晶体中心位置定位，本方法本质上是对图像中的非波峰部分进行连续衰减，从而突出波峰部分，对于一维或者二维寻峰算法而言，极大地提高了寻峰的准确度。同时，本发明还提供了晶体像素查找表生成方法。
虽然本发明披露如上，但本发明并非限定于此。任何本领域技术人员，在不脱离本发明的精神和范围内，均可作各种更动与修改，因此本发明的保护范围应当以权利要求所限定的范围为准。

Claims (10)

1.一种晶体中心位置图生成方法，其特征在于，包括如下步骤：
获取光子事件二维直方图I，将所述光子事件二维直方图I归一化至[0，1]之间，获得图像I0；
对所述图像I0进行迭代处理得到图像Ii；
使用自动阈值对所述图像Ii进行二值化处理，获得二值图像集SB；
提取所述二值图像集SB中的每幅二值图像的连通域及其对应的中心位置；
判断所述连通域彼此是否有重叠；
保留有重叠的连通域中迭代次数最大的连通域；
生成晶体中心位置图，所述晶体中心位置图包括有重叠的连通域中迭代次数最大的连通域对应的中心位置以及没有重叠的连通域对应的中心位置。
2.如权利要求1所述的晶体中心位置图生成方法，其特征在于，判断所述连通域彼此是否有重叠的具体步骤为：
将一个连通域内的任意位置与另一个连通域内的任意位置进行比较，如果相同即表示这两个连通域彼此有重叠。
3.如权利要求1所述的晶体中心位置图生成方法，其特征在于，所述保留有重叠的连通域中迭代次数最大的连通域的具体步骤为：
比较所述有重叠的连通域的对应图像Ii的迭代次数i，剔除对应迭代次数较小的连通域。
4.如权利要求1所述的晶体中心位置图生成方法，其特征在于，所述有重叠的连通域中仅保留迭代次数最大的连通域的具体步骤为：
比较所述彼此有重叠的连通域的面积，剔除面积较大的连通域。
5.如权利要求1所述的晶体中心位置图生成方法，其特征在于，还包括设置面积阈值，剔除连通域面积大于所述面积阈值的连通域。
6.如权利要求1所述的晶体中心位置图生成方法，其特征在于，所述将光子事件二维直方图I归一化至[0，1]之间采用线性收放的方式进行，即 max(I)是所述光子事件二维直方图中的最大值，min(I)是所述光子事件二维直方图中的最小值。
7.如权利要求6所述的晶体中心位置图生成方法，其特征在于，所述光子事件二维直方图最大值max(I)按照比例ε获得，即取所述光子事件二维直方图的灰度直方图面积的p处的灰度为最大值max(I)，p＝1-ε，图中超出最大值的max(I)的值强制转换成max(I)。
8.如权利要求7所述的晶体中心位置图生成方法，其特征在于，所述比例其中n是所述光子事件二维直方图中实际包含的晶体数量，Np是光子事件二维直方图的像素数量。
9.如权利要求1所述的晶体中心位置图生成方法，其特征在于，所述自动阈值的选取基于光子事件二维直方图的灰度直方图，所述阈值将所述光子事件二维直方图的灰度直方图分成两部分，所述两部分的质心到所述自动阈值的距离相等。
10.一种晶体像素查找表生成方法，其特征在于，采用如权利要求1-9任一项所述的方法生成的晶体中心位置图，根据所述晶体中心位置图划分各相邻晶体条的分界线，生成晶体像素查找表。
CN201510214482.5A 2015-04-29 2015-04-29 晶体中心位置图生成方法及晶体像素查找表生成方法 Active CN104809460B (zh)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN201510214482.5A CN104809460B (zh) 2015-04-29 2015-04-29 晶体中心位置图生成方法及晶体像素查找表生成方法
CN201611266541.4A CN106611173B (zh) 2015-04-29 2015-04-29 晶体中心位置图生成方法
US15/142,295 US9928437B2 (en) 2015-04-29 2016-04-29 Method and system for crystal identification
US15/897,110 US10176393B2 (en) 2015-04-29 2018-02-14 Method and system for crystal identification

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510214482.5A CN104809460B (zh) 2015-04-29 2015-04-29 晶体中心位置图生成方法及晶体像素查找表生成方法

Related Child Applications (1)

Application Number Title Priority Date Filing Date
CN201611266541.4A Division CN106611173B (zh) 2015-04-29 2015-04-29 晶体中心位置图生成方法

Publications (2)

Publication Number Publication Date
CN104809460A CN104809460A (zh) 2015-07-29
CN104809460B true CN104809460B (zh) 2017-09-26

Family

ID=53694271

Family Applications (2)

Application Number Title Priority Date Filing Date
CN201510214482.5A Active CN104809460B (zh) 2015-04-29 2015-04-29 晶体中心位置图生成方法及晶体像素查找表生成方法
CN201611266541.4A Active CN106611173B (zh) 2015-04-29 2015-04-29 晶体中心位置图生成方法

Family Applications After (1)

Application Number Title Priority Date Filing Date
CN201611266541.4A Active CN106611173B (zh) 2015-04-29 2015-04-29 晶体中心位置图生成方法

Country Status (1)

Country Link
CN (2) CN104809460B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106611173A (zh) * 2015-04-29 2017-05-03 上海联影医疗科技有限公司 晶体中心位置图生成方法

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106526651B (zh) * 2015-09-14 2019-08-23 上海联影医疗科技有限公司 一种探测器晶体位置表的建立方法及系统
CN109884682B (zh) * 2019-02-01 2020-12-01 东软医疗系统股份有限公司 一种晶体位置查找表的生成方法、装置、设备及介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101978288A (zh) * 2008-03-21 2011-02-16 株式会社岛津制作所 二维位置图校准方法以及放射线检测装置
CN103914860A (zh) * 2013-01-05 2014-07-09 苏州瑞派宁科技有限公司 晶体条位置查找表生成方法及装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI398664B (zh) * 2009-10-30 2013-06-11 Iner Aec Executive Yuan 閃爍偵檢器之晶體位置感測訊號識別與排序方法
CN104809460B (zh) * 2015-04-29 2017-09-26 上海联影医疗科技有限公司 晶体中心位置图生成方法及晶体像素查找表生成方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101978288A (zh) * 2008-03-21 2011-02-16 株式会社岛津制作所 二维位置图校准方法以及放射线检测装置
CN103914860A (zh) * 2013-01-05 2014-07-09 苏州瑞派宁科技有限公司 晶体条位置查找表生成方法及装置

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106611173A (zh) * 2015-04-29 2017-05-03 上海联影医疗科技有限公司 晶体中心位置图生成方法
CN106611173B (zh) * 2015-04-29 2020-06-02 上海联影医疗科技有限公司 晶体中心位置图生成方法

Also Published As

Publication number Publication date
CN106611173B (zh) 2020-06-02
CN104809460A (zh) 2015-07-29
CN106611173A (zh) 2017-05-03

Similar Documents

Publication Publication Date Title
CN109378052B (zh) 图像标注的预处理方法及系统
CN104778688B (zh) 点云数据的配准方法及装置
CN102567300B (zh) 图片文档的处理方法及装置
CN108921166A (zh) 基于深度神经网络的医疗票据类文本检测识别方法及系统
CN106683106B (zh) 一种晶体像素查找表生成的方法
CN104200444B (zh) 基于对称样本块的图像修复方法
CN112200117B (zh) 表格识别方法及装置
CN109242858B (zh) 基于自适应模板匹配的织物印花循环图案基元分割方法
CN110619333B (zh) 一种文本行分割方法、文本行分割装置及电子设备
CN104809460B (zh) 晶体中心位置图生成方法及晶体像素查找表生成方法
CN109087276B (zh) 基于smt料盘的x射线图像的元器件自动计数和定位方法
CN104680531B (zh) 一种连通量统计信息提取方法及vlsi结构
CN112990183B (zh) 离线手写汉字同名笔画提取方法、系统、装置
WO2021168703A1 (zh) 字符处理及字符识别方法、存储介质和终端设备
CN114612490B (zh) 一种基于显微镜图像的栅藻类细胞统计方法
CN111724354A (zh) 一种基于图像处理的多株小麦穗长与小穗数的测量方法
CN106526651A (zh) 一种探测器晶体位置表的建立方法及系统
CN112329641A (zh) 一种表格识别方法、装置、设备及可读存储介质
JP2017500662A (ja) 投影ひずみを補正するための方法及びシステム
CN110751732A (zh) 一种2d图像转3d图像的方法
CN107085855A (zh) 一种快速圆弧检测方法
CN104112135B (zh) 文本图像提取装置以及方法
CN102072707A (zh) 数字图像中圆的中心与半径快速检测方法
CN111968115A (zh) 基于栅格化图像处理方法的骨科耗材检测方法及系统
CN105427332B (zh) 一种用于水平集图像分割的快速符号距离函数计算方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CP01 Change in the name or title of a patent holder

Address after: 201807 No. 2258 Chengbei Road, Jiading Industrial Zone, Jiading District, Shanghai.

Patentee after: Shanghai Lianying Medical Technology Co., Ltd

Address before: 201807 No. 2258 Chengbei Road, Jiading Industrial Zone, Jiading District, Shanghai.

Patentee before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd.

CP01 Change in the name or title of a patent holder