WO2015054946A1 - Ct投影旋转中心的确定方法和装置 - Google Patents

Ct投影旋转中心的确定方法和装置 Download PDF

Info

Publication number
WO2015054946A1
WO2015054946A1 PCT/CN2013/087534 CN2013087534W WO2015054946A1 WO 2015054946 A1 WO2015054946 A1 WO 2015054946A1 CN 2013087534 W CN2013087534 W CN 2013087534W WO 2015054946 A1 WO2015054946 A1 WO 2015054946A1
Authority
WO
WIPO (PCT)
Prior art keywords
sequence
projection
measured
cross
detector
Prior art date
Application number
PCT/CN2013/087534
Other languages
English (en)
French (fr)
Inventor
孟凡勇
李忠传
李静海
Original Assignee
中国科学院过程工程研究所
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 中国科学院过程工程研究所 filed Critical 中国科学院过程工程研究所
Publication of WO2015054946A1 publication Critical patent/WO2015054946A1/zh

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
    • G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
    • 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
    • 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/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A61B6/032—Transmission computed tomography [CT]
    • 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

Definitions

  • the present invention relates to the field of medical digital imaging and non-destructive testing technologies, and in particular, to a method and apparatus for determining a CT projection rotation center. Background technique
  • CT Computer Tomography
  • COR Center Of Rotation
  • the line connecting the focus of the ray source and the center of the projection rotation should be perpendicular to the detector and correspond to the central pixel of the detector, but it is difficult for the three to be completely and accurately adjusted during the installation of the CT, if not for the center of rotation of the projection Correction, the resulting reconstructed image will have artifacts, affecting the measurer to truly understand the cross-section, causing the measurer to misjudge.
  • the existing measurement correction methods of the projection rotation center include a direct measurement method, a model correction method, and a symmetric relationship method.
  • the direct measurement method is to determine the projection rotation center of the detector by directly measuring the source-detector of the CT and the rotation center of the rotary table. This method is limited by the measurement error and is difficult to implement in practical applications.
  • the model correction rule uses a special thin-walled cylinder model or a needle-shaped model. The calibration model is scanned first, the projection rotation center is determined by iterative calculation, and then the measured object is formally measured, using a dedicated model. It is cumbersome to measure the resulting projection center data to correct the official measurement data.
  • the symmetric relation method is to use the symmetry of the original data, and to determine the projection rotation center by adding or subtracting the calculation of the combination or the difference.
  • Such methods are limited by the object to be measured, projection blur, and projection angle, and cannot be accurately measured. Summary of the invention
  • the present invention provides a method and a device for determining a CT projection rotation center, which can be quickly performed by simply scanning the object to be measured and performing simple calculation on the collected raw data without using a dedicated calibration phantom. Accurately determine the projected rotation center of the CT to be tested.
  • an embodiment of the present invention provides a method for determining a CT projection rotation center, including: rotating a test object according to a predetermined incremental step size using a CT to be tested, and acquiring different rotation angles, each detection The projection value of the pixel point p ⁇ , ) ; where n, and N are integers, n [0,N], N> ⁇ , me[ ⁇ ,M], M is the total number of pixels of the detector, M is an integer greater than 1; each projection value of the detector pixel at the different rotation angles is obtained, and the projection values are divided into a first sequence p m (W) and a second sequence ⁇ + 180' ⁇ ) , where ⁇ is an integer, e [0, ⁇ ] ;
  • the rotating the scan of the measured object by using the CT to be tested according to the predetermined incremental step comprises:
  • the ray source of the CT to be tested is rotated according to a predetermined incremental step, and the rotating table of the object to be measured is fixed. Do not move, use the CT to be tested to scan the object under rotation; or
  • the rotating table on which the object to be measured is placed is rotated in a predetermined incremental step, and the source of the CT to be tested is fixed, and the object to be tested is rotated and scanned using the CT to be tested.
  • ⁇ (The number of correlations between ⁇ +180 specifically includes:
  • the method further includes: calibrating an original projection value of each detector pixel of the CT to be tested according to a determination result of a projection rotation center of the CT to be tested, and generating according to the calibrated projection data Corresponding reconstruction of the image.
  • an embodiment of the present invention provides a CT projection rotation center determining apparatus, including: a projection value acquiring unit, configured to perform a rotation scan on an object to be tested according to a predetermined incremental step size, and acquire The projection value p( , ) of each detector pixel at different rotation angles; where ⁇ and N are integers, " e [0,N], N>- 1 ⁇ -, me[l, ], M is The total number of detector pixels,
  • M is an integer greater than one
  • a classification sequence acquiring unit configured to acquire respective projection values of the detector pixel at the different rotation angles, and divide the projection values into a first sequence ⁇ ( ⁇ ) and a second sequence ⁇ + 180"), Where k
  • the cross-correlation sequence calculation unit is configured to separately calculate the correlation coefficient between the first sequence ⁇ ( ⁇ ) and the second sequence ⁇ ( ⁇ + 180") in each detector pixel to form a cross-correlation sequence? a central pixel point search unit for traversing the sequence of cross-correlation numbers (?), finding and ?() The detector pixel corresponding to the maximum value in the column, as the center pixel;
  • a rotation center determining unit configured to determine a projection rotation center of the CT to be tested according to a position of the central pixel point.
  • the projection value acquiring unit is specifically configured to:
  • the ray source of the CT to be tested rotates according to a predetermined incremental step size, and the rotating table of the object to be measured is fixed, and the object to be tested is rotated and scanned by using the CT to be tested, and the projection of the pixel points of each detector is obtained under different rotation angles.
  • the rotating table on which the object to be measured is placed is rotated according to a predetermined incremental step, and the radiation source of the CT to be tested is fixed, and the object to be tested is rotated and scanned by using the CT to be tested, and the projection of the pixel points of each detector is obtained under different rotation angles. Value W , ).
  • cross-correlation sequence calculation unit is specifically configured to:
  • the apparatus further includes: a calibration reconstruction unit:
  • the object to be tested is rotated and scanned, and the projection values of the detecting units in the detector array under different rotation angles are obtained, and the projection values of the same detecting unit at different rotation angles are 180°.
  • FIG. 1 is a schematic diagram showing characteristics of a CT scan projection in the prior art
  • FIG. 2 is a flow chart of a method for determining a CT projection rotation center according to a first embodiment of the present invention
  • FIG. 3 is a third embodiment of the present invention.
  • Figure 4 is a sequence diagram of the calculated cross-correlation coefficient of the second embodiment of the present invention.
  • FIG. 5 is a two-dimensional reconstruction diagram of a non-projection rotation center correction at the time of full-cycle sampling according to the second embodiment of the present invention.
  • FIG. 6 is a two-dimensional weight map corrected by a projection rotation center during full-cycle sampling according to a second embodiment of the present invention.
  • FIG. 7 is a projection sine of incomplete full-circumferential scan data with a perforated disk as a measurement object according to a third embodiment of the present invention.
  • Figure 8 is a sequence diagram of the calculated cross-correlation coefficient of the third embodiment of the present invention.
  • FIG. 9 is a structural diagram of a CT projection rotation center determining apparatus according to a fourth embodiment of the present invention, specifically J Jr
  • FIG. 1 A schematic diagram of the characteristics of a CT scan projection is shown in FIG. Among them, the CT machine in the third generation fan beam scanning mode is shown in FIG.
  • the three-generation CT includes a radiation source 110 and a detector array 120, wherein the emission signal of the radiation source 110 has a certain opening angle, and the typical angular angle range is 30 to 45 °, and the detector array 120 is disposed at the radiation source 110. Within the sector angle of the signal coverage, the detector array 120 includes a plurality of detector units, each of which is a detector pixel. At the time of scanning, the ray source 110 and the detector array 120 rotate around the axis of the object 130 to be measured.
  • the line connecting the focus of the ray source 110 with the center of the projection rotation should be perpendicular to the detector array 120, that is, the Y-axis position in FIG. 1, but in reality, it is difficult for the three to be completely precisely adjusted.
  • the central ray tends to have a certain offset from the Y-axis position.
  • the ray source 110 is offset somewhat, and therefore, the central ray 140 is offset from the Y-axis by an angle. At this time, if the image is reconstructed with the Y-axis as the center, there will be artifacts.
  • the same set of data for the two sets is repeated at 180 degrees. Therefore, by performing data processing on the projection data of the same detector pixel at different rotation angles, the detector pixel with the highest correlation is obtained. At point 150, the projected rotation center of the CT to be tested can be determined.
  • Step 110 Perform a rotation scan on the measured object according to a predetermined incremental step, using the CT to be tested, and obtain a projection value p( , ) of each detector pixel at different rotation angles. ; wherein, n, and N are integers, w e [0, N] , N> ⁇ , me [ ⁇ , M], M is a total number of detector pixels, M being greater than
  • the central processing unit of the CT to be tested controls the CT to be tested, and performs rotational scanning on the measured object according to a predetermined incremental step size.
  • the rotating scan of the measured object by using the CT to be tested according to the predetermined incremental step includes: the ray source of the CT to be tested is rotated according to a predetermined incremental step, and the rotating table of the measured object is fixed, and the test is performed.
  • CT scans the object under test; or
  • the rotating table on which the object to be measured is placed is rotated in a predetermined incremental step, and the source of the CT to be tested is fixed, and the object to be tested is rotated and scanned using the CT to be tested.
  • the incremental step size can be set according to the actual situation. As an example and not a limitation, it can be set to 0.05°, 0.1° or 0.2°.
  • ⁇ is not completely scanned all week ( ⁇ N ⁇ ).
  • N can also It is considered that an integer greater than ⁇ is not limited thereto.
  • the projection values of the pixel points of each detector form a matrix P under different rotation angles
  • Step 120 Acquire each projection value of the detector pixel at the different rotation angles, and divide the projection values into a first sequence ⁇ ( ⁇ ) and a second sequence ⁇ + 1 80 , where ⁇ is an integer ,
  • the central processing unit of the CT to be tested acquires each projection value of each detector pixel at the different rotation angles, and maps all projection values corresponding to each detector pixel to 180°.
  • the boundary is divided into a first sequence and a second sequence.
  • the first sequence is the same as the number of data included in the second sequence.
  • the central processor acquires each column in the matrix ,, and further divides each column of data of the ⁇ into a first sequence and a second sequence according to different rotation angles.
  • Step 130 Calculate a correlation coefficient between the first sequence ⁇ ( ⁇ ) and the second sequence ⁇ ( ⁇ +180") in each detector pixel to form a cross-correlation sequence ?().
  • the first sequence ( ⁇ ) and the second sequence of each detector pixel point are calculated.
  • the number of correlations between ( ⁇ +180") specifically includes:
  • Step 140 Traverse the sequence of cross-correlation coefficients ((), and find a detector pixel corresponding to the maximum value in the sequence as a central pixel point.
  • the cross-correlation sequence ?() is a function of the pixel of the detector, and by traversing ?(), the detector pixel corresponding to the maximum value is found, that is, the central pixel point.
  • Step 150 Determine a projection rotation center of the CT to be tested according to a position of the central pixel point.
  • the position of the center pixel is used as the projection rotation center of the CT to be tested.
  • the object to be tested is rotated and scanned, and the projection values of the detecting units in the detector array under different rotation angles are obtained, and the projection values of the same detecting unit at different rotation angles are 180°.
  • Dividing into two sets of sequences respectively calculating the correlation coefficient of the two sets of sequences in each detecting unit, obtaining a sequence of cross-correlation numbers, traversing the sequence of the cross-correlation numbers, and searching for a detecting unit corresponding to the maximum value in the sequence,
  • the technical means of determining the rotation center of the projection is solved, and the measurement correction method needs to use a dedicated model.
  • the calculation is cumbersome and the symmetry relationship method has many restrictions and cannot be accurately measured. It is not necessary to use a dedicated model, but only through direct use. By scanning the measured object and performing simple calculations on the collected raw data, the projected rotation center of the CT to be tested can be quickly and accurately determined.
  • the method further includes: performing calibration on the original projection value of each detector pixel of the CT to be tested according to the determination result of the projection rotation center of the CT to be tested, according to the calibration
  • the projection data is generated to generate a corresponding reconstructed image.
  • the CT to be tested is a first generation CT machine that scans using a pan-and-rotation mode and includes only one detector unit. At each rotation angle, the ray source and the corresponding detector unit scan the object to be measured in a translational manner. The translation position corresponding to each rotation angle is M. Then, the projection values of the corresponding M detector pixel points at each rotation angle are M projection values of the detector unit at M translation positions.
  • the CT to be tested is a second generation CT machine that scans using a translation-rotation mode and includes only one detector unit, where i is an integer greater than one.
  • the source and the corresponding detector unit scan the object to be measured in a translational manner.
  • the translation position corresponding to each rotation angle is one, where J- is an integer greater than one.
  • i *j M.
  • the projection values of the corresponding M detector pixel points at each rotation angle are M projection values of the detector unit at j translation positions.
  • the CT to be tested is a third-to-fifth generation CT machine including M detector units scanned using a rotational-rotation method. Then, the projection values of the corresponding M detector pixels at each rotation angle are M projection values of the M detector units at corresponding rotation angles.
  • the high-voltage generator using the third-generation equidistant fan beam CT as the CTo ray source system to be tested is produced by German Yxlon, model MGG40, X.
  • the light pipe is produced by Philips (Philips), model Y-TU/100-D01;
  • the detector is an equidistant array produced by Anbangshi (Beijing) Technology Co., Ltd., model LSC-412, and the scintillator material is GOS (Gadolinium) Oxysulfide, sulfur oxidation ceremony
  • the detector array has a total of 1536 detector units, also That is, 1536 detector pixels, the pixel size is 0.3*0.6mm, and the pixel pitch is 0.4mm.
  • the X-ray operating voltage used in this experiment is 100kVp, and the current is 2.2mA.
  • the stepping motor drives the stage and rotates 360° for sampling.
  • the incremental step of the sampling angle is 0. ⁇ .
  • the object to be measured adopts a porous plexiglass model.
  • the reconstruction algorithm uses a standard FBP (Filtered Back-Projection) algorithm using an mmp (ramp) filter. All reconstruction results are not post-image processed, so they can effectively reflect the real situation determined by the projection rotation center.
  • the measured object is rotated and scanned by using the CT to be tested to obtain the projection value p( ⁇ , ) of each detector pixel under different rotation angles; wherein, n, and N are Integer, "e [0,3600], me [1,1536] , the sinogram of the original projection data of the acquired full-circumference scan without logarithm operation, as shown in Figure 3.
  • the traversal maximal value is obtained.
  • a full-week scan of the CT to be tested produces excessive redundant information. Therefore, it is possible to perform an incomplete full-week scan of the CT to be tested based on actual conditions or empirical values.
  • the present embodiment Based on the CT to be tested selected in the second embodiment, the present embodiment performs 210° rotation sampling, and the incremental step size of the sampling is 0.
  • the object to be measured adopts a porous plexiglass model.
  • the reconstruction algorithm uses the standard FBP algorithm, using a ramp filter. All reconstruction results are not post-image processed, so it can effectively reflect the real situation determined by the projection rotation center.
  • the object to be tested is rotated and scanned by the CT to be tested, and the projection values of the pixels of each detector are obtained under different rotation angles, where n, and N are integers, “e [0] , 2100], me [1,1536] , the acquired sinogram of the original projection data of the incomplete full-circumference scan without logarithm operation, as shown in Fig. 7.
  • the traversal maximal value is obtained.
  • the incomplete full-circle scan can obtain the same projection rotation center as the full-circle scan. It can be understood that the larger the maximum rotation angle is selected, the calculation is calculated. The projection of the rotation center of the more accurate, but the measurement of the amount of redundancy greater t
  • FIG. 9 is a structural diagram of a CT rotation center determining device according to a fourth embodiment of the present invention. As shown in Figure 9, the device includes:
  • the projection value obtaining unit 910 is configured to perform a rotation scan on the measured object according to a predetermined incremental step size using the CT to be tested, and obtain a projection value p( ⁇ , ) of each detector pixel point under different rotation angles ; wherein, n , and N are integers, "e[0,N], N> ⁇ , me[l, ], M is the pixel of the detector
  • M is an integer greater than 1;
  • a classification sequence obtaining unit 920 configured to acquire respective projection values of the detector pixel points at the different rotation angles, and divide the projection values into a first sequence ⁇ ( ⁇ ) and a second sequence ⁇ +180), Among them, win 180' ⁇
  • ⁇ is an integer
  • the ⁇ ⁇ cross-correlation sequence calculation unit 930 is configured to calculate the correlation coefficient between the first sequence ? m (W) and the second sequence ⁇ ( ⁇ + 180) in each detector pixel, respectively a sequence of relationship numbers? ( ); a central pixel point searching unit 940, configured to traverse the sequence of cross-correlation numbers ( ), and find a detector pixel corresponding to the maximum value in the sequence as a central pixel point;
  • the rotation center determining unit 950 is configured to determine a projection rotation center of the CT to be tested according to the position of the central pixel point.
  • the object to be tested is rotated and scanned, and the projection values of the detecting units in the detector array under different rotation angles are obtained, and the projection values of the same detecting unit at different rotation angles are 180°.
  • the cross-correlation coefficient obtains a sequence of cross-correlation coefficients, traverses the sequence of the cross-correlation coefficients, finds a detecting unit corresponding to the maximum value in the sequence, and serves as a central detecting unit, thereby determining a technical means of projecting the rotating center, and solving the measurement correcting method It is necessary to use a dedicated model, and the technical problems of cumbersome and symmetrical relationship methods are too many to be accurately measured. It is not necessary to use a dedicated model. It is only by directly scanning the measured object and performing simple calculation on the collected raw data. And accurately determine the projection rotation center of the CT to be tested.
  • the projection value acquiring unit is specifically configured to:
  • the ray source of the CT to be tested rotates according to a predetermined incremental step size, and the rotating table of the object to be measured is fixed, and the object to be tested is rotated and scanned by using the CT to be tested, and the projection of the pixel points of each detector is obtained under different rotation angles.
  • the rotating table on which the object to be measured is placed is rotated according to a predetermined incremental step, and the radiation source of the CT to be tested is fixed, and the object to be tested is rotated and scanned by using the CT to be tested, and the projection of the pixel points of each detector is obtained under different rotation angles.
  • the cross-correlation sequence calculation unit is specifically configured to: sequentially acquire a first sequence and a second sequence of a detector pixel ⁇ +180":); a sequence pl ⁇ and a sequence ⁇ ( ⁇ + 18 (0 is a cross-correlation operation to obtain a corresponding cross-correlation sequence; the maximum value in the cross-correlation sequence is obtained as the cross-correlation coefficient of the pixel of the detector.
  • the apparatus further includes: a calibration reconstruction unit:
  • the determining device of the CT rotating center provided by the embodiment of the present invention can be used to perform the determining method of the CT rotating center provided by any embodiment of the present invention, and has a corresponding functional module to achieve the same technology. Effect.
  • the embodiments of the present invention may be implemented by a program executable by a computer device, so that they may be stored in a storage device and executed by a processor, and the program may be stored in a computer.
  • the above mentioned storage medium may be a read only memory, a magnetic disk or an optical disk, etc.; or they may be separately fabricated into individual integrated circuit modules, or a plurality of modules or steps thereof may be fabricated into a single integrated circuit module.
  • the invention is not limited to any specific combination of hardware and software.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

提供一种CT投影旋转中心的确定方法和装置。所述方法包括：按照预定的递增步长Θ，使用待测CT对被测物体进行旋转扫描，获取不同旋转角度下，各探测器像素点的投影值，获取探测器像素点m在不同旋转角度下的各投影值，将各投影值分为第一序列和第二序列，分别计算各探测器像素点中的第一序列和第二序列之间的互相关系数，组成互相关系数序列R（m）；遍历互相关系数序列R（m），査找与R（m）序列中的最大值相对应的探测器像素点，作为中心像素点；根据中心像素点的位置，确定待测CT的投影旋转中心。所述方法和装置无需使用专用的校正模体，仅通过直接扫描被测物体，并对采集的原始数据进行简单的运算，即可快速而精确的确定待测CT的投影旋转中心。

Description

说 明 书
CT投影旋转中心的确定方法和装置
技术领域
本发明涉及医学数字成像和无损检测技术领域，尤其涉及一种 CT投影旋转 中心的确定方法和装置。 背景技术
CT (Computed Tomography, 计算机断层扫描）最早应用于医学影像， 随后 被引入工业无损检测领域， 由于其非插入、 无干扰的检测特性， CT在农林业、 地球物理、 化工等领域也得到了很好的应用。 根据射线源-探测器的运动采集数 据方式的不同， 可以将 CT分为五代， 目前较为常用的为平移-旋转的第一代笔 束扫描模式和旋转-旋转的第三代扇束扫描模式。 精确重建出被测截面的二维信 息的前提，需要获得准确的 CT投影几何参数，其中 CT投影旋转中心（Center Of Rotation, COR) 是非常重要的一项参数。 理论上， 射线源的焦点与投影旋转中 心的连线应该垂直于探测器， 并且对应于探测器的中心像素点， 但 CT的安装过 程中此三者难以完全精确调节一致， 若不对投影旋转中心进行校正， 所得到的 重建图像将存在伪影， 影响测量者真实地了解截面情况， 致使测量者产生误判。
现有的投影旋转中心的测量校正方法包括直接测量法、 模型校正法以及对 称关系法。 直接测量法是通过直接测量 CT的源-探测器及旋转台的旋转中心来 确定探测器的投影旋转中心， 该方法受限于测量误差， 实际应用中较难实现。 模型校正法则采用专用的薄壁圆筒模型或针状模型， 先对校正模型进行扫描， 通过迭代计算确定投影旋转中心， 然后对被测物进行正式测量， 使用专用模型 测量所得的投影中心数据来校正正式测量数据， 较为繁琐。 对称关系法等是利 用原始数据存在的对称性， 通过加减求合或求差的运算， 确定投影旋转中心。 但此类方法受被测物、 投影模糊及投影角度的限制， 也无法做到精确测量。 发明内容
有鉴于此， 本发明提供一种 CT投影旋转中心的确定方法和装置， 无需使用 专用的校正模体， 仅通过直接扫描被测物体， 并对采集的原始数据进行简单的 运算， 即可快速而精确的确定待测 CT的投影旋转中心。
在第一方面，本发明实施例提供了一种 CT投影旋转中心的确定方法，包括: 按照预定的递增步长^ 使用待测 CT对被测物体进行旋转扫描， 获取不同 旋转角度下， 各探测器像素点的投影值 p ^， ); 其中， n， 和 N均为整数， n [0,N], N>^， me[\,M], M为探测器像素点的总个数， M为大于 1的整数; 获取探测器像素点 在所述不同旋转角度下的各投影值， 将所述各投影值 分为第一序列 pm(W)和第二序列 ^ + 180'·)， 其中， έ为整数， e[0,^^];
Θ 分别计算各探测器像素点中的第一序列 Pmm和第二序列 ^„^+ΐ8(π之间 的互相关系数， 组成互相关系数序列 ?( );
遍历所述互相关系数序列 ， 査找与 ?( )序列中的最大值相对应的探测 器像素点， 作为中心像素点；
根据所述中心像素点的位置， 确定所述待测 CT的投影旋转中心。
在第一种可能的实现方式中， 所述按照预定的递增步长^ 使用待测 CT对 被测物体进行旋转扫描包括：
待测 CT的射线源按照预定的递增步长 旋转， 放置被测物体的旋转台固定 不动， 使用待测 CT对被测物体进行旋转扫描； 或者
放置被测物体的旋转台按照预定的递增步长 旋转， 待测 CT的射线源固定 不动， 使用待测 CT对被测物体进行旋转扫描。
进一步的， 所述计算各探测器像素点的第一序列 pm(w)和第二序列
^(^+180之间的互相关系数具体包括：
依次获取一个探测器像素点的第一序列 和第二序列 ^ ^+180"：)； 对序列 p l^和序列 ^(^ + 180'·)做互相关运算， 得到相应的互相关序列； 获取所述互相关序列中的最大值作为该探测器像素点的互相关系数。
进一步的，所述方法还包括：根据所述待测 CT的投影旋转中心的确定结果， 对所述待测 CT的各探测器像素点的原始投影值进行校准，根据校准后的投影数 据， 生成相应的重建图像。
在第二方面，本发明实施例提供了一种 CT投影旋转中心的确定装置，包括: 投影值获取单元， 用于按照预定的递增步长^ 使用待测 CT对被测物体进 行旋转扫描，获取不同旋转角度下，各探测器像素点的投影值 p( ， );其中， ^ 和 N均为整数， 《e[0,N]， N>-1^-, me[l, ], M为探测器像素点的总个数，
Θ
M为大于 1的整数;
分类序列获取单元， 用于获取探测器像素点 在所述不同旋转角度下的各 投影值， 将所述各投影值分为第一序列 ^(^)和第二序列 ^ ^ + 180")， 其中， k
Figure imgf000005_0001
互相关系数序列计算单元， 用于分别计算各探测器像素点中的第一序列 ^(^)和第二序列 ^(^ + 180")之间的互相关系数， 组成互相关系数序列 ?( )； 中心像素点査找单元， 用于遍历所述互相关系数序列 ?( )， 査找与 ?( )序 列中的最大值相对应的探测器像素点， 作为中心像素点；
旋转中心确定单元， 用于根据所述中心像素点的位置，确定所述待测 CT的 投影旋转中心。
在第一种可能的实现方式中， 所述投影值获取单元具体用于：
待测 CT的射线源按照预定的递增步长 旋转， 放置被测物体的旋转台固定 不动， 使用待测 CT对被测物体进行旋转扫描， 获取不同旋转角度下， 各探测器 像素点的投影值 p ^， ) ; 或者
放置被测物体的旋转台按照预定的递增步长 旋转， 待测 CT的射线源固定 不动， 使用待测 CT对被测物体进行旋转扫描， 获取不同旋转角度下， 各探测器 像素点的投影值 W , )。
进一步的， 所述互相关系数序列计算单元具体用于：
依次获取一个探测器像素点的第一序列 和第二序列 ^ ^+180 ； 对序列 p i^和序列 ^ ^+iscr)做互相关运算， 得到相应的互相关序列； 获取所述互相关序列中的最大值作为该探测器像素点的互相关系数。
进一步的， 所述装置还包括， 校准重建单元：
用于根据所述待测 CT的投影旋转中心的确定结果， 对所述待测 CT的各探 测器像素点的原始投影值进行校准， 根据校准后的投影数据， 生成相应的重建 图像。
本发明实施例通过使用待测 CT对被测物体进行旋转扫描，获取不同旋转角 度下， 探测器阵列中各探测单元的投影值， 对同一探测单元在不同旋转角度下 的投影值， 以 180° 为界分为两组序列， 分别计算各探测单元中的上述两组序列 的互相关系数， 得到互相关系数序列， 遍历上述互相关系数序列， 査找与序列 中的最大值相对应的探测单元， 作为中心探测单元， 进而确定投影旋转中心的 技术手段， 解决了测量校正方法需要使用专用模型， 计算繁琐以及对称关系法 限制条件较多， 无法精准测量的技术问题， 无需使用专用模型， 仅通过直接扫 描被测物体， 并对采集的原始数据进行简单的运算， 即可快速而精确的确定待 测 CT的投影旋转中心。
图 1是现有技术中一种 CT扫描投影的特性示意图；
图 2是本发明第一实施例的一种 CT投影旋转中心确定方法的流程图； 图 3 是本发明第
弦图；
图 4是本发明第二实施例的计算得到的互相关系数序列图；
图 5 是本发明第二实施例的全周采样时未经投影旋转中心校正的二维重建 图；
图 6是本发明第二实施例的全周采样时经过投影旋转中心校正后的二维重 图 7是本发明第三实施例的以多孔圆盘为测量对象的不完全全周扫描数据 投影正弦图；
图 8是本发明第三实施例的计算得到的互相关系数序列图；
图 9是本发明第四实施例的一种 CT投影旋转中心确定装置的结构图, 具体实 J ^r式
为了使本发明的目的、 技术方案和优点更加清楚， 下面结合附图对本发明 具体实施例作进一步的详细描述。 可以理解的是， 此处所描述的具体实施例仅 仅用于解释本发明， 而非对本发明的限定。 另外还需要说明的是， 为了便于描 述， 附图中仅示出了与本发明相关的部分而非全部内容。 在图 1中示出了一种 CT扫描投影的特性示意图。其中， 图 1中示出了第三 代扇束扫描模式下的 CT机。 该三代 CT包括射线源 110和探测器阵列 120， 其 中，射线源 110的发射信号具有一定的张角， 典型的张角角度范围为， 30〜45 ° ， 探测器阵列 120设置于射线源 110发射信号覆盖的扇形角之内，探测器阵列 120 中包括多个探测器单元， 每一个探测器单元为一探测器像素点。 在扫描时， 射 线源 110和探测器阵列 120同步的围绕着被测物体 130的轴心进行旋转。
理论上， 射线源 110 的焦点与投影旋转中心的连线 （中心射线） 应该垂直 于探测器阵列 120， 也即图 1 中的 Y轴位置， 但是实际上此三者难以完全精确 调节一致， 因此， 中心射线往往与 Y轴位置存在一定的偏移， 例如在图 1中， 射线源 110发生了一定的偏移， 因此， 中心射线 140与 Y轴偏移了 角的。 此 时， 如果以 Y轴为中心射线重建图像， 会存在伪影。
通过对 CT的投影过程进行分析可知， 当通过 CT对待测物体进行扫描时， 在进行投影数据全周采集（射线源 110旋转 360° 进行扫描）过程中， 投影信息 是冗余的， 每一条投影路径， 均有两条射线从正反方向通过， 因此， 全周扫描 原始数据中， 在不考虑噪声的前提下， 总存在两个完全相同的投影数据。 其中， 只有中心射线 140对应的投影路径， 其正反两组数据均出现在同一个探测器像 素点 150上， 其他投影路径对应的正反两组数据点， 需根据几何关系确定其出 现的探测器像素点以及相应的旋转角度差。 而且， 对于中心射线 140， 此两组相 同的数据是以 180度为界， 重复出现的。 因此， 通过对不同旋转角度下， 同一 探测器像素点的投影数据进行数据处理后， 获取其中相关度最高的探测器像素 点 150， 即可确定待测 CT的投影旋转中心。 第一实施例
图 2是本发明第一实施例的一种 CT投影旋转中心确定方法的流程图，本实 施例的方法可以由 CT投影旋转中心确定装置来执行， 该装置可通过硬件和 /或 软件的方式实现，一般可集成于 CT机内部。本实施例的方法具体包括如下步骤： 步骤 110、按照预定的递增步长^ 使用待测 CT对被测物体进行旋转扫描， 获取不同旋转角度下， 各探测器像素点的投影值 p( ， ); 其中， n， 和 N均 为整数， we[0,N]， N>^， me[\,M], M为探测器像素点的总个数， M为大于
Θ
1的整数。
在本实施例中，待测 CT的中央处理器控制待测 CT，按照预定的递增步长^ 对被测物体进行旋转扫描。
其中，按照预定的递增步长^使用待测 CT对被测物体进行旋转扫描包括: 待测 CT的射线源按照预定的递增步长 旋转， 放置被测物体的旋转台固定 不动， 使用待测 CT对被测物体进行旋转扫描； 或者
放置被测物体的旋转台按照预定的递增步长 旋转， 待测 CT的射线源固定 不动， 使用待测 CT对被测物体进行旋转扫描。
其中， 递增步长 可以根据实际情况进行设定， 作为示例而非限定， 可以将 设定为 0.05° ,0.1° 或者 0.2° 等。 在本实施例中， 可以使用待测 CT对被测物体进行全周扫描（N = ^)， 或
Θ 者不完全全周扫描 （^<N<^)。 当然， 本领域技术人员可以理解， N也可 以为大于 ^的整数， 对此并不限定。
Θ 在本实施例中，不同旋转角度下，各探测器像素点的投影值构成一个矩阵 P，
Figure imgf000010_0001
p
Figure imgf000010_0002
其中， 矩阵 Ρ的行向量代表同一旋转角度下， Μ个探测器像素点对应的 Μ 水投影值； 矩阵 Ρ的列向量代表 N 个不同的旋转角度下， 同一个探测器像素点
Figure imgf000010_0003
步骤 120、 获取探测器像素点 在所述不同旋转角度下的各投影值， 将所述 各投影值分为第一序列 ^(^)和第二序列 ^ ^+ 1 80 ， 其中， έ为整数，
Figure imgf000010_0004
在本实施例中，待测 CT的中央处理器获取各探测器像素点在所述不同旋转 角度下的各投影值， 并将对应于每一个探测器像素点的所有投影值， 以 180° 为 界， 分为第一序列和第二序列。
其中， 第一序列与第二序列中包括的数据个数相同。
在本实施例中， 中央处理器获取矩阵 Ρ中的每一列， 根据不同的旋转角度， 将 Ρ的每一列数据进一步分为第一序列和第二序列。
步骤 130、 分别计算各探测器像素点中的第一序列 ^(^)和第二序列 ^(^+180")之间的互相关系数， 组成互相关系数序列 ?( )。
在本实施例中， 计算各探测器像素点的第一序列 (^)和第二序列 (^+180")之间的互相关系数具体包括：
依次获取一个探测器像素点中的第一序列 pm,和第二序列 ^ ^ + 18(0； 对序列 p l^和序列 ^(^ + 180'·)做互相关运算， 得到相应的互相关序列； 获取所述互相关序列中的最大值作为该探测器像素点的互相关系数。
步骤 140、遍历所述互相关系数序列 ?( )，査找与 序列中的最大值相对 应的探测器像素点， 作为中心像素点。
在本实施例中， 互相关序列 ?( )为探测器像素点 的函数， 通过遍历 ?( )， 査找其中的最大值对应的探测器像素点， 即为中心像素点。
步骤 150、根据所述中心像素点的位置，确定所述待测 CT的投影旋转中心。 在本实施例中， 将待测 CT在旋转角度为 0° 时， 中心像素点所在的位置， 作为该待测 CT的投影旋转中心。
本发明实施例通过使用待测 CT对被测物体进行旋转扫描，获取不同旋转角 度下， 探测器阵列中各探测单元的投影值， 对同一探测单元在不同旋转角度下 的投影值， 以 180° 为界分为两组序列， 分别计算各探测单元中的上述两组序列 的互相关系数， 得到互相关系数序列， 遍历上述互相关系数序列， 査找与序列 中的最大值相对应的探测单元， 作为中心探测单元， 进而确定投影旋转中心的 技术手段， 解决了测量校正方法需要使用专用模型， 计算繁琐以及对称关系法 限制条件较多， 无法精准测量的技术问题， 无需使用专用模型， 仅通过直接扫 描被测物体， 并对采集的原始数据进行简单的运算， 即可快速而精确的确定待 测 CT的投影旋转中心。
在上述技术方案的基础上， 所述方法还包括：根据所述待测 CT的投影旋转 中心的确定结果，对所述待测 CT的各探测器像素点的原始投影值进行校准， 根 据校准后的投影数据， 生成相应的重建图像。 在本实施例的一个优选的实施方式中， 待测 CT为使用平移-旋转方式扫描 并且仅包括一个探测器单元的第一代 CT机。在每一个旋转角度下， 射线源与对 应的探测器单元按照平移的方式， 对待测物体进行扫描。 每一个旋转角度对应 的平移位置为 M个。则， 每个旋转角度下对应的 M个探测器像素点的投影值为 该探测器单元在 M个平移位置下的 M个投影值。
在本实施例的一个优选的实施方式中， 待测 CT为使用平移-旋转方式扫描 并且仅包括个探测器单元的第二代 CT机， 其中， i为大于 1的整数。在每一个 旋转角度下， 射线源与对应的 探测器单元按照平移的方式， 对待测物体进行扫 描。每一个旋转角度对应的平移位置为 ·个，其中， J-为大于 1的整数。 i *j = M。 则， 每个旋转角度下对应的 M个探测器像素点的投影值为该 /个探测器单元在 j 个平移位置下的 M个投影值。
在本实施例的一个优选的实施方式中， 待测 CT为使用旋转-旋转方式扫描 的包括 M个探测器单元的第三代-第五代 CT机。 则， 每个旋转角度下对应的 M 个探测器像素点的投影值为该 M个探测器单元在相应的旋转角度下的 M个投影 值。 第二实施例
在上述实施例的基础上， 在本实施例中， 采用第三代等距扇束 CT作为待测 CTo射线源系统的高压发生器为德国 Yxlon (依科视朗）出品， 型号为 MGG40, X光管为 Philips (飞利浦） 公司生产， 型号为 Y-TU/100-D01 ; 探测器为安邦世 (北京） 科技公司出品的等距阵列， 型号为 LSC-412 , 闪烁体材料为 GOS (Gadolinium Oxysulfide, 硫氧化礼）， 探测器阵列共有 1536个探测器单元， 也 即 1536个探测器像素点， 像素点尺寸为 0.3*0.6mm， 像素间距为 0.4mm。 本实 验采用的 X射线工作电压为 100kVp， 电流为 2.2mA， 步进电机驱动载物台， 旋 转 360° 进行采样， 采样角度的递增步长为 0. Γ 。 被测物体采用多孔有机玻璃 模型。重建算法采用标准的 FBP (Filtered Back-Projection,滤波反向投影）算法， 使用 mmp (斜坡） 滤波器。 所有重建结果均未经图像后处理， 因此可有效反应 投影旋转中心确定的真实情况。
按照 0. Γ 的递增步长， 使用上述待测 CT对被测物体进行旋转扫描， 获取 不同旋转角度下， 各探测器像素点的投影值 p(^， ) ; 其中， n， 和 N均为整 数， 《e [0,3600]， m e [1,1536] , 采集得到的未经过对数运算的全周扫描的原始投影 数据正弦图， 如图 3所示。
获取探测器像素点 在所述不同旋转角度下的各投影值， 将所述各投影值 分为第一序列 ?m(W)和第二序列 m(^ + 180")， 其中， έ为整数， A e [0，1800]。
分别计算各探测器像素点中的第一序列 Pmm和第二序列 ^„ ^+18(Π之间 的互相关系数， 组成互相关系数序列 ?( )， 图 4 中示出了计算得到的互相关系 数序列图。 在本实施例中， 对互相关系数序列 做了相应的归一化处理。
对互相关系数序列 ?( )进行遍历求极大值， 在 m=801号探测器像素点上， 产 生互相关系数的极大值， ?(801) = 0.9787。 由图 4可以看出， 中心射线的互相关系 数显著高于周围像素点， 因此， 可以确定投影旋转中心对应的中心像素为 801 号像素， 该像素与探测器的物理中心像素 1536/2=768偏移了 33个像素点。
未经投影旋转中心校正的重建图像如图 5 所示， 可以看出， 重建图像中出 现了明显的伪影。 使用 801 号像素进行投影旋转中心校正后的重建图像如图 6 所示， 可以看出： 由于投影旋转中心得到了精确的校正， 伪影得到了消除， 重 建图像的质量得到极大提升。 第三实施例
一般来说， 对待测 CT做全周扫描会产生过多的冗余信息， 因此， 可以根据 实际情况或者经验值对待测 CT做不完全全周扫描。
在与实施例二选择完全相同的待测 CT的基础上， 本实施例进行 210° 旋转 采样， 采样的递增步长为 0.Γ 。 被测物体采用多孔有机玻璃模型。 重建算法采 用标准的 FBP算法， 使用 ramp滤波器。所有重建结果均未经图像后处理， 因此 可有效反映投影旋转中心确定的真实情况。
按照 0.Γ 的递增步长， 使用上述待测 CT对被测物体进行旋转扫描， 获取 不同旋转角度下， 各探测器像素点的投影值 其中， n， 和 N均为整 数， 《e [0, 2100]， m e [1,1536] , 采集得到的未经过对数运算的不完全全周扫描的原 始投影数据正弦图， 如图 7所示。
获取探测器像素点 在所述不同旋转角度下的各投影值， 将所述各投影值 分为第一序列 ?m(W)和第二序列 m(^ + 18(T)， 其中， έ为整数， e [0，300]。
分别计算各探测器像素点中的第一序列 pmm和第二序列 ^(^+18(0之间 的互相关系数， 组成互相关系数序列 ?( )， 图 8 中示出了计算得到的互相关系 数序列图。 在本实施例中， 对互相关系数序列 做了相应的归一化处理。
对互相关系数序列 ?( )进行遍历求极大值， 在 m=801号探测器像素点上， 产 生互相关系数的极大值， ?(801) = 0.9787。 由图 8可以看出， 中心射线的互相关系 数显著高于周围像素点， 因此， 可以确定投影旋转中心对应的中心像素为 801 号像素， 该像素与探测器的物理中心像素 1536/2=768偏移了 33个像素点。
可以看出， 设定合适的最大旋转角度后， 不完全全周扫描可以与全周扫描 得到相同的投影旋转中心， 可以理解的是， 最大旋转角度选择的越大， 计算得 到的投影旋转中心也就越准， 但是测量冗余量也就越大 t
应当说明的是， 本发明并不限于对圆盘模型的¾ ：有效， 本发明还可以用 于其它 CT无损检测领域中， 这是本领域普通技术人员容易理解的。 图 9为本发明第四实施例的一种 CT旋转中心确定装置的结构图。如图 9所 示， 所述装置包括：
投影值获取单元 910， 用于按照预定的递增步长^ 使用待测 CT对被测物 体进行旋转扫描， 获取不同旋转角度下， 各探测器像素点的投影值 p(^， ); 其 中， n， 和 N均为整数， "e[0，N]， N>^, me[l, ], M为探测器像素点的
Θ 总个数， M为大于 1的整数；
分类序列获取单元 920，用于获取探测器像素点 在所述不同旋转角度下的 各投影值， 将所述各投影值分为第一序列 ^(^)和第二序列 ^ ^+180 )， 其中， 勝 180'··
έ为整数，
~~ θ ~ 互相关系数序列计算单元 930，用于分别计算各探测器像素点中的第一序列 ?m(W)和第二序列 ^(^ + 180)之间的互相关系数， 组成互相关系数序列 ?( )； 中心像素点査找单元 940，用于遍历所述互相关系数序列 ?( )，査找与 序列中的最大值相对应的探测器像素点， 作为中心像素点；
旋转中心确定单元 950， 用于根据所述中心像素点的位置， 确定所述待测 CT的投影旋转中心。
本发明实施例通过使用待测 CT对被测物体进行旋转扫描，获取不同旋转角 度下， 探测器阵列中各探测单元的投影值， 对同一探测单元在不同旋转角度下 的投影值， 以 180° 为界分为两组序列， 分别计算各探测单元中的上述两组序列 的互相关系数， 得到互相关系数序列， 遍历上述互相关系数序列， 査找与序列 中的最大值相对应的探测单元， 作为中心探测单元， 进而确定投影旋转中心的 技术手段， 解决了测量校正方法需要使用专用模型， 计算繁琐以及对称关系法 限制条件较多， 无法精准测量的技术问题， 无需使用专用模型， 仅通过直接扫 描被测物体， 并对采集的原始数据进行简单的运算， 即可快速而精确的确定待 测 CT的投影旋转中心。
在上述各实施例的基础上， 所述投影值获取单元具体用于：
待测 CT的射线源按照预定的递增步长 旋转， 放置被测物体的旋转台固定 不动， 使用待测 CT对被测物体进行旋转扫描， 获取不同旋转角度下， 各探测器 像素点的投影值 p ^， ) ; 或者
放置被测物体的旋转台按照预定的递增步长 旋转， 待测 CT的射线源固定 不动， 使用待测 CT对被测物体进行旋转扫描， 获取不同旋转角度下， 各探测器 像素点的投影值 p( , )。
在上述各实施例的基础上， 所述互相关系数序列计算单元具体用于： 依次获取一个探测器像素点的第一序列 和第二序列 ^ ^+180"：)； 对序列 p l^和序列 ^(^ + 18(0做互相关运算， 得到相应的互相关序列； 获取所述互相关序列中的最大值作为该探测器像素点的互相关系数。
在上述各实施例的基础上，， 所述装置还包括， 校准重建单元：
用于根据所述待测 CT的投影旋转中心的确定结果， 对所述待测 CT的各探 测器像素点的原始投影值进行校准， 根据校准后的投影数据， 生成相应的重建 图像。
本发明实施例所提供的 CT旋转中心的确定装置可以用于执行本发明任意 实施例提供的 CT旋转中心的确定方法， 具备相应的功能模块， 达到同样的技术 效果。 显然， 本领域技术人员应该明白， 本发明实施例可以用计算机装置可执行 的程序来实现， 从而可以将它们存储在存储装置中由处理器来执行， 所述的程 序可以存储于一种计算机可读存储介质中， 上述提到的存储介质可以是只读存 储器， 磁盘或光盘等； 或者将它们分别制作成各个集成电路模块， 或者将它们 中的多个模块或步骤制作成单个集成电路模块来实现。 这样， 本发明不限制于 任何特定的硬件和软件的结合。
以上所述仅为本发明的优选实施例， 并不用于限制本发明， 对于本领域技 术人员而言， 本发明可以有各种改动和变化。 凡在本发明的精神和原理之内所 作的任何修改、 等同替换、 改进等， 均应包含在本发明的保护范围之内。

Claims

权 利 要 求 书
1、 一种 CT投影旋转中心的确定方法， 其特征在于， 包括：
按照预定的递增步长^ 使用待测 CT对被测物体进行旋转扫描， 获取不同 旋转角度下， 各探测器像素点的投影值 p(^， ); 其中， n， 和 N均为整数， n^[0,N], Ν>-1^, ,η ,Μ], Μ为探测器像素点的总个数， Μ为大于 1 的整 数；
获取探测器像素点 在所述不同旋转角度下的各投影值， 将所述各投影值 分为第一序列 ^)和第二序列 ^^ + 18(Τ)， 其中， 为整数， [0，^ ;
Θ 分别计算各探测器像素点中的第一序列 Α ^)和第二序列 ^+ΐ8(Τ)之间 的互相关系数， 组成互相关系数序列
遍历所述互相关系数序列 R( )， 查找与 R( )序列中的最大值相对应的探测 器像素点， 作为中心像素点；
根据所述中心像素点的位置， 确定所述待测 CT的投影旋转中心。
2、 根据权利要求 1所述的 CT投影旋转中心的确定方法， 其特征在于， 所 述按照预定的递增歩长^ 使用待测 CT对被测物体进行旋转扫描包括：
待测 CT的射线源按照预定的递增步长 e旋转，放置被测物体的旋转台固定 不动， 使用待测 CT对被测物体进行旋转扫描； 或者
放置被测物体的旋转台按照预定的递增步长 ^旋转，待测 CT的射线源固定 不动， 使用待测 CT对被测物体进行旋转扫描。
3、 根据权利要求 2所述的 CT投影旋转中心的确定方法， 其特征在于， 所 述计算各探测器像素点的第一序列 Pm 和第二序列 (^+18(Τ)之间的互相关 系数具体包括：
依次获取一个探测器像素点的第一序列 (^)和第二序列 ^ ^+18(Π； 对序列 Pm 和序列 ^^+iscr)做互相关运算， 得到相应的互相关序列； 获取所述互相关序列中的最大值作为该探测器像素点的互相关系数。
4、 根据权利要求 2所述的 CT投影旋转中心的确定方法， 其特征在于， 所 述方法还包括： 根据所述待测 CT的投影旋转中心的确定结果， 对所述待测 CT 的各探测器像素点的原始投影值进行校准， 根据校准后的投影数据， 生成相应 的重建图像。
5、 一种 CT投影旋转中心的确定装置， 其特征在于， 包括：
投影值获取单元， 用于按照预定的递增步长^ 使用待测 CT对被测物体进 行旋转扫描， 获取不同旋转角度下， 各探测器像素点的投影值 p( ， ); 其中， n, 和 N均为整数， 《e[0，N]， N>-1^-, me[l, ], M为探测器像素点的总个 数， M为大于 1的整数；
分类序列获取单元， 用于获取探测器像素点 在所述不同旋转角度下的各 投影值， 将所述各投影值分为第一序列 ^(^)和第二序列 ^+18(Π， 其中， k
Figure imgf000019_0001
互相关系数序列计算单元， 用于分别计算各探测器像素点中的第一序列 ?m(W)和第二序列 (^ + 18(Π之间的互相关系数， 组成互相关系数序列 ?( )； 中心像素点査找单元， 用于遍历所述互相关系数序列 ?( )， 査找与 序 列中的最大值相对应的探测器像素点， 作为中心像素点；
旋转中心确定单元， 用于根据所述中心像素点的位置， 确定所述待测 CT 的投影旋转中心。
6、 根据权利要求 5所述的 CT投影旋转中心的确定装置， 其特征在于， 所 述投影值获取单元具体用于： 待测 CT的射线源按照预定的递增步长 旋转，放置被测物体的旋转台固 不动， 使用待测 CT对被测物体进行旋转扫描， 获取不同旋转角度下， 各探 器像素点的投影值 P( ， ) ; 或者
放置被测物体的旋转台按照预定的递增步长 旋转，待测 CT的射线源固 不动， 使用待测 CT对被测物体进行旋转扫描， 获取不同旋转角度下， 各探 器像素点的投影值 p( ， )。
7、 根据权利要求 6所述的 CT投影旋转中心的确定装置， 其特征在于， 述互相关系数序列计算单元具体用于：
依次获取一个探测器像素点的第一序列 和第二序列 ^ ^+180"：)； 对序列 p l^和序列 ^( + ISO")做互相关运算， 得到相应的互相关序列; 获取所述互相关序列中的最大值作为该探测器像素点的互相关系数。
8、 根据权利要求 6所述的 CT投影旋转中心的确定装置， 其特征在于， 述装置还包括， 校准重建单元：
用于根据所述待测 CT的投影旋转中心的确定结果，对所述待测 CT的各 测器像素点的原始投影值进行校准， 根据校准后的投影数据， 生成相应的重 图像。
PCT/CN2013/087534 2013-10-17 2013-11-20 Ct投影旋转中心的确定方法和装置 WO2015054946A1 (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201310488193.5 2013-10-17
CN201310488193.5A CN103606144A (zh) 2013-10-17 2013-10-17 Ct投影旋转中心的确定方法和装置

Publications (1)

Publication Number Publication Date
WO2015054946A1 true WO2015054946A1 (zh) 2015-04-23

Family

ID=50124364

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2013/087534 WO2015054946A1 (zh) 2013-10-17 2013-11-20 Ct投影旋转中心的确定方法和装置

Country Status (2)

Country Link
CN (1) CN103606144A (zh)
WO (1) WO2015054946A1 (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104007131A (zh) * 2014-05-23 2014-08-27 中国科学院过程工程研究所 便携式ct扫描设备、便携式ct系统和ct检测方法
CN105510362B (zh) * 2015-12-07 2019-01-08 华中农业大学 基于微型ct的水稻分蘖性状无损测量装置及其测量方法
CN106651977B (zh) * 2016-09-30 2020-03-31 重庆大学 基于重建图像梯度的l0范数最小化的锥束ct旋转中心标定方法
CN111242913B (zh) * 2020-01-08 2023-04-11 浙江大学 一种获取肋骨展开图像的方法、系统、装置及存储介质

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101584587A (zh) * 2008-05-21 2009-11-25 中国科学院过程工程研究所 一种ct投影中心的自动校正方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3846576B2 (ja) * 2002-07-03 2006-11-15 株式会社島津製作所 コンピュータ断層撮影装置
CN102539460B (zh) * 2012-01-06 2014-07-02 公安部第一研究所 一种ct系统投影旋转中心定位方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101584587A (zh) * 2008-05-21 2009-11-25 中国科学院过程工程研究所 一种ct投影中心的自动校正方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LI, BAOLEI ET AL.: "Automatic determination of the projection center of rotation of the 2D-CT system based on symmetry projection", SEMINAR, PROCEEDINGS, NATIONAL RAY DIGITAL IMAGING AND CT NEW TECHNOLOGY, 30 June 2009 (2009-06-30), CHINA , SHANGHAI, pages 64 - 72 *
MENG, FANYONG ET AL.: "Determination of the center of rotation for computed tomography based on projection data", THE 13TH CHINA STEREOLOGY AND IMAGE ANALYSIS CONFERENCE PROCEEDINGS, 23 September 2013 (2013-09-23), pages 390 - 395 *

Also Published As

Publication number Publication date
CN103606144A (zh) 2014-02-26

Similar Documents

Publication Publication Date Title
RU2606561C2 (ru) Компенсация усечения для итерационной реконструкции в компьютерной томографии (кт) с коническим пучком в комбинированных системах офэкт/кт
JP6482815B2 (ja) Ｘ線コンピュータ断層撮影装置及びリファレンス補正プログラム
RU2619839C1 (ru) Способ исследования с помощью кт и устройство для кт
JP2007117736A5 (zh)
WO2015074254A1 (zh) Ct探测器偏转角的确定方法和装置
WO2007148263A1 (en) Method and system for error compensation
US20180174335A1 (en) X- Ray CT Device and Sequential Correction Parameter Determination Method
JP2005021675A (ja) 断層撮影装置
JP7219888B2 (ja) Ｘ線画像システムの較正
WO2010135901A1 (zh) 伪双能欠采样物质识别系统和方法
AU2016200833B2 (en) A computed tomography imaging process and system
CN111728632B (zh) 射线探测装置、射线探测方法和ct图像重建方法
JP2009534631A (ja) ダーティアイソトープのｐｅｔ再構成
US20130235970A1 (en) Image Processing Device of a Computer Tomography System
WO2015054946A1 (zh) Ct投影旋转中心的确定方法和装置
JP2020531083A (ja) 撮像システム拡張視野
JP4977007B2 (ja) 複数の焦点取得方法及び装置
US7769217B2 (en) Fast iterative 3D PET image reconstruction using a set of 2D linogram transformations
Yang et al. A review of geometric calibration for different 3-D X-ray imaging systems
Fin et al. A practical way to improve contrast‐to‐noise ratio and quantitation for statistical‐based iterative reconstruction in whole‐body PET imaging
WO2017188559A1 (en) Method of image reconstruction in computed tomography
JP4792187B2 (ja) オブジェクトデータセットからの横断面分布の導出
WO2016190218A1 (ja) 測定方法および放射線撮影装置
Ghani et al. An empirical method for geometric calibration of a photon counting detector-based cone beam CT system
CN112748135B (zh) 一种外部ct扫描模式下转台旋转中心标定方法及装置

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 13895737

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 13895737

Country of ref document: EP

Kind code of ref document: A1