WO2017148115A1 - Phantom projection simulation method - Google Patents

Phantom projection simulation method Download PDF

Info

Publication number
WO2017148115A1
WO2017148115A1 PCT/CN2016/096857 CN2016096857W WO2017148115A1 WO 2017148115 A1 WO2017148115 A1 WO 2017148115A1 CN 2016096857 W CN2016096857 W CN 2016096857W WO 2017148115 A1 WO2017148115 A1 WO 2017148115A1
Authority
WO
WIPO (PCT)
Prior art keywords
boundary
projection
phantom
simulation algorithm
projection simulation
Prior art date
Application number
PCT/CN2016/096857
Other languages
French (fr)
Chinese (zh)
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 WO2017148115A1 publication Critical patent/WO2017148115A1/en

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects

Definitions

  • the invention relates to a phantom projection simulation algorithm and belongs to the technical field of medical image processing.
  • X-ray cone beam 3D reconstruction is a hotspot in the field of medical image processing. This current method mainly needs to be improved in accuracy and computational efficiency.
  • the algorithm involves high difficulty and long calculation time is a technology that is currently plaguing many experts. point.
  • the object of the present invention is to provide a phantom projection simulation algorithm, which adopts the Feldkamp algorithm to further improve the computational efficiency of projection simulation by examining the ray path.
  • a phantom projection simulation algorithm comprising the following steps:
  • the above-described phantom projection simulation algorithm is characterized in that, in the step 1), the detector adopts a virtual XCT imaging system.
  • the aforementioned phantom projection simulation algorithm is characterized in that the tissue boundary appears in the step 2) It is based on the fact that the CT value of the current point is the same as the CT value of the previous point.
  • the aforementioned phantom projection simulation algorithm is characterized in that if a tissue boundary occurs, a boundary approximation strategy is adopted to reduce the discrete error.
  • the above-described phantom projection simulation algorithm is characterized in that the boundary approximation strategy continuously divides the initial step size across the boundary and decides whether to convert the moving direction until reaching a preset maximum dichotomy number, and finally at the current The position is restored to the initial step size.
  • the above-mentioned phantom projection simulation algorithm is characterized in that the calculation formula of the projection result in the step 3) is adopted.
  • l N represents the last boundary position of the penetration path L
  • E min and E max are the minimum and maximum values in the energy spectrum distribution, respectively.
  • the method determines the correct tissue boundary by examining the CT value on the path, and uses the Feldkamp algorithm to reconstruct the image, thereby improving the computational efficiency of the projection simulation.
  • the invention relates to a phantom projection simulation algorithm, comprising the following steps:
  • the detector uses a virtual XCT imaging system.
  • the boundary approximation strategy is adopted to reduce the discrete error.
  • the boundary approximation strategy is to continuously divide the initial step size across the boundary and decide whether to change the moving direction until reaching The preset maximum number of dichotomies, and finally the initial step size is restored at the current position.
  • the method determines the correct tissue boundary by examining the CT value on the path, and uses the Feldkamp algorithm to reconstruct the image, thereby improving the computational efficiency of the projection simulation.

Abstract

A phantom projection simulation method is characterized by comprising the following steps: 1) setting a parameter of a projection procedure detector; 2) determining a location of a source point of an X-ray and a penetration path L of the X-ray, examining CT values of all points on the path L, and determining a tissue boundary; 3) computing a power spectral distribution function and a projection result; and 4) adopting a Feldkamp algorithm to reconstruct an image. The embodiment determines the accurate tissue boundary by examining the CT values on the examination path, and increases efficiency of computing projection simulation by adopting the Feldkamp algorithm to reconstruct the image.

Description

一种模体投影仿真算法A phantom projection simulation algorithm 技术领域Technical field
本发明涉及一种模体投影仿真算法,属于医学图像处理技术领域。The invention relates to a phantom projection simulation algorithm and belongs to the technical field of medical image processing.
背景技术Background technique
X射线锥形束三维重建是目前医学图像处理领域研究的热点,这一块目前的方法主要在精确度以及计算效率上还有待提高,算法涉及难度高,计算时间较长是目前困扰很多专家的技术点。X-ray cone beam 3D reconstruction is a hotspot in the field of medical image processing. This current method mainly needs to be improved in accuracy and computational efficiency. The algorithm involves high difficulty and long calculation time is a technology that is currently plaguing many experts. point.
发明内容Summary of the invention
为解决现有技术的不足,本发明的目的在于提供一种模体投影仿真算法,通过对射线路径的考察,采用Feldkamp算法,从而进一步提高投影仿真的计算效率。In order to solve the deficiencies of the prior art, the object of the present invention is to provide a phantom projection simulation algorithm, which adopts the Feldkamp algorithm to further improve the computational efficiency of projection simulation by examining the ray path.
为了实现上述目标,本发明采用如下的技术方案:In order to achieve the above objectives, the present invention adopts the following technical solutions:
一种模体投影仿真算法,其特征是,包括如下步骤:A phantom projection simulation algorithm, comprising the following steps:
1)设置投影过程检测器的参数;1) setting the parameters of the projection process detector;
2)确定X射线源点位置以及X射线的穿透路径L,考察路径L上的所有点的CT值,确定组织边界;2) determining the X-ray source point position and the X-ray penetration path L, examining the CT values of all points on the path L, and determining the tissue boundary;
3)计算X射线能谱分布函数以及投影结果;3) calculating an X-ray energy spectrum distribution function and a projection result;
4)采用Feldkamp算法重建图像。4) Reconstruct the image using the Feldkamp algorithm.
前述的一种模体投影仿真算法,其特征是,所述步骤1)中,检测器采用虚拟XCT成像系统。The above-described phantom projection simulation algorithm is characterized in that, in the step 1), the detector adopts a virtual XCT imaging system.
前述的一种模体投影仿真算法,其特征是,所述步骤2)中组织边界出现的 依据是当前点的CT值与前一点的CT值相同。The aforementioned phantom projection simulation algorithm is characterized in that the tissue boundary appears in the step 2) It is based on the fact that the CT value of the current point is the same as the CT value of the previous point.
前述的一种模体投影仿真算法,其特征是,若出现组织边界,采用边界逼近策略减少离散误差。The aforementioned phantom projection simulation algorithm is characterized in that if a tissue boundary occurs, a boundary approximation strategy is adopted to reduce the discrete error.
前述的一种模体投影仿真算法,其特征是,所述边界逼近策略为对跨越边界的初始步长不断做二分并决定是否转换移动方向,直至达到预先设定的最大二分次数,最后在当前位置恢复初始步长。The above-described phantom projection simulation algorithm is characterized in that the boundary approximation strategy continuously divides the initial step size across the boundary and decides whether to convert the moving direction until reaching a preset maximum dichotomy number, and finally at the current The position is restored to the initial step size.
前述的一种模体投影仿真算法,其特征是,所述步骤3)中投影结果的计算公式采用
Figure PCTCN2016096857-appb-000001
其中,lN表示穿透路径L的最后一个边界位置,Emin和Emax分别为能谱分布中的最小值和最大值。
The above-mentioned phantom projection simulation algorithm is characterized in that the calculation formula of the projection result in the step 3) is adopted.
Figure PCTCN2016096857-appb-000001
Wherein, l N represents the last boundary position of the penetration path L, and E min and E max are the minimum and maximum values in the energy spectrum distribution, respectively.
本发明所达到的有益效果:本方法通过考察考察路径上的CT值,确定正确的组织边界,采用Feldkamp算法进行重建图像,提高了投影仿真的计算效率。The beneficial effects achieved by the invention: the method determines the correct tissue boundary by examining the CT value on the path, and uses the Feldkamp algorithm to reconstruct the image, thereby improving the computational efficiency of the projection simulation.
具体实施方式detailed description
下面对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。The invention is further described below. The following examples are only intended to more clearly illustrate the technical solutions of the present invention, and are not intended to limit the scope of the present invention.
本发明涉及一种模体投影仿真算法,包括如下步骤:The invention relates to a phantom projection simulation algorithm, comprising the following steps:
1)设置投影过程检测器的参数,本实施例中检测器采用虚拟XCT成像系统。1) Set the parameters of the projection process detector. In this embodiment, the detector uses a virtual XCT imaging system.
2)确定X射线源点位置以及X射线的穿透路径L,考察路径L上的所有点的CT值,确定组织边界,组织边界出现的依据是当前点的CT值与前一点的CT值相同。2) Determine the X-ray source point position and the X-ray penetration path L, examine the CT value of all points on the path L, and determine the tissue boundary. The basis of the tissue boundary is that the CT value of the current point is the same as the CT value of the previous point. .
若出现组织边界,采用边界逼近策略减少离散误差,本实施例中边界逼近策略为对跨越边界的初始步长不断做二分并决定是否转换移动方向,直至达到 预先设定的最大二分次数,最后在当前位置恢复初始步长。If the boundary of the tissue appears, the boundary approximation strategy is adopted to reduce the discrete error. In this embodiment, the boundary approximation strategy is to continuously divide the initial step size across the boundary and decide whether to change the moving direction until reaching The preset maximum number of dichotomies, and finally the initial step size is restored at the current position.
3)计算X射线能谱分布函数以及投影结果,本实施例中投影结果的计算公式采用
Figure PCTCN2016096857-appb-000002
其中,lN表示穿透路径L的最后一个边界位置,Emin和Emax分别为能谱分布中的最小值和最大值。
3) Calculating the X-ray energy spectrum distribution function and the projection result, the calculation formula of the projection result in this embodiment is adopted
Figure PCTCN2016096857-appb-000002
Wherein, l N represents the last boundary position of the penetration path L, and E min and E max are the minimum and maximum values in the energy spectrum distribution, respectively.
4)采用Feldkamp算法重建图像。4) Reconstruct the image using the Feldkamp algorithm.
本发明所达到的有益效果:本方法通过考察考察路径上的CT值,确定正确的组织边界,采用Feldkamp算法进行重建图像,提高了投影仿真的计算效率。The beneficial effects achieved by the invention: the method determines the correct tissue boundary by examining the CT value on the path, and uses the Feldkamp algorithm to reconstruct the image, thereby improving the computational efficiency of the projection simulation.
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。 The above is only a preferred embodiment of the present invention, and it should be noted that those skilled in the art can make several improvements and modifications without departing from the technical principles of the present invention. It should also be considered as the scope of protection of the present invention.

Claims (6)

  1. 一种模体投影仿真算法,其特征是,包括如下步骤:A phantom projection simulation algorithm, comprising the following steps:
    1)设置投影过程检测器的参数;1) setting the parameters of the projection process detector;
    2)确定X射线源点位置以及X射线的穿透路径L,考察路径L上的所有点的CT值,确定组织边界;2) determining the X-ray source point position and the X-ray penetration path L, examining the CT values of all points on the path L, and determining the tissue boundary;
    3)计算X射线能谱分布函数以及投影结果;3) calculating an X-ray energy spectrum distribution function and a projection result;
    4)采用Feldkamp算法重建图像。4) Reconstruct the image using the Feldkamp algorithm.
  2. 根据权利要求1所述的一种模体投影仿真算法,其特征是,所述步骤1)中,检测器采用虚拟XCT成像系统。A phantom projection simulation algorithm according to claim 1, wherein in the step 1), the detector employs a virtual XCT imaging system.
  3. 根据权利要求1所述的一种模体投影仿真算法,其特征是,所述步骤2)中组织边界出现的依据是当前点的CT值与前一点的CT值相同。A phantom projection simulation algorithm according to claim 1, wherein the tissue boundary appears in the step 2) based on the fact that the CT value of the current point is the same as the CT value of the previous point.
  4. 根据权利要求3所述的一种模体投影仿真算法,其特征是,若出现组织边界,采用边界逼近策略减少离散误差。A phantom projection simulation algorithm according to claim 3, wherein if a tissue boundary occurs, a boundary approximation strategy is employed to reduce the discrete error.
  5. 根据权利要求4所述的一种模体投影仿真算法,其特征是,所述边界逼近策略为对跨越边界的初始步长不断做二分并决定是否转换移动方向,直至达到预先设定的最大二分次数,最后在当前位置恢复初始步长。A phantom projection simulation algorithm according to claim 4, wherein the boundary approximation strategy is to continuously divide the initial step size across the boundary and decide whether to switch the moving direction until the preset maximum dichotomy is reached. The number of times, and finally the initial step size is restored at the current position.
  6. 根据权利要求1所述的一种模体投影仿真算法,其特征是,所述步骤3)中投影结果的计算公式采用
    Figure PCTCN2016096857-appb-100001
    其中,lN表示穿透路径L的最后一个边界位置,Emin和Emax分别为能谱分布中的最小值和最大值。
    A phantom projection simulation algorithm according to claim 1, wherein the calculation formula of the projection result in the step 3) is
    Figure PCTCN2016096857-appb-100001
    Wherein, l N represents the last boundary position of the penetration path L, and E min and E max are the minimum and maximum values in the energy spectrum distribution, respectively.
PCT/CN2016/096857 2016-02-29 2016-08-26 Phantom projection simulation method WO2017148115A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201610108658.3A CN105719338A (en) 2016-02-29 2016-02-29 Motif projection simulation algorithm
CN2016101086583 2016-02-29

Publications (1)

Publication Number Publication Date
WO2017148115A1 true WO2017148115A1 (en) 2017-09-08

Family

ID=56156186

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2016/096857 WO2017148115A1 (en) 2016-02-29 2016-08-26 Phantom projection simulation method

Country Status (2)

Country Link
CN (1) CN105719338A (en)
WO (1) WO2017148115A1 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105719338A (en) * 2016-02-29 2016-06-29 江苏美伦影像系统有限公司 Motif projection simulation algorithm

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102609908A (en) * 2012-01-13 2012-07-25 中国人民解放军信息工程大学 Base image TV model based CT (Computed Tomography) beam hardening correcting method
US20130267841A1 (en) * 2009-02-11 2013-10-10 Siemens Medical Solutions Usa, Inc. Extracting Application Dependent Extra Modal Information from an Anatomical Imaging Modality for use in Reconstruction of Functional Imaging Data
CN105719338A (en) * 2016-02-29 2016-06-29 江苏美伦影像系统有限公司 Motif projection simulation algorithm

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6504892B1 (en) * 2000-10-13 2003-01-07 University Of Rochester System and method for cone beam volume computed tomography using circle-plus-multiple-arc orbit
US6560308B1 (en) * 2001-10-26 2003-05-06 Kabushiki Kaisha Toshiba Method and system for approximating missing data in cone beam x-ray CT reconstruction
CA2605836A1 (en) * 2005-04-25 2006-11-02 University Of Rochester Method and apparatus of global de-noising for ct imaging
CN101126722B (en) * 2007-09-30 2011-03-16 西北工业大学 Cone-beam CT beam hardening calibration method based on registration model emulation
CN101472381B (en) * 2007-12-29 2013-03-27 Ge医疗系统环球技术有限公司 Method for controlling X ray exposure of X ray CT system
CN102521853A (en) * 2011-12-08 2012-06-27 苏州生物医学工程技术研究所 Reconstruction algorithm for back projection weight cone-beam CT (Computed Tomography)

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130267841A1 (en) * 2009-02-11 2013-10-10 Siemens Medical Solutions Usa, Inc. Extracting Application Dependent Extra Modal Information from an Anatomical Imaging Modality for use in Reconstruction of Functional Imaging Data
CN102609908A (en) * 2012-01-13 2012-07-25 中国人民解放军信息工程大学 Base image TV model based CT (Computed Tomography) beam hardening correcting method
CN105719338A (en) * 2016-02-29 2016-06-29 江苏美伦影像系统有限公司 Motif projection simulation algorithm

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
TANG, SHAOJIE ET AL.,: "Analytic Calculation of CT Phantom Projections and Its Application in Beam Hardening Correction", JOURNAL OF XI`AN JIAOTONG UNIVERSITY, vol. 42, no. 6, 1 July 2008 (2008-07-01) *
TANG, SHAOJIE ET AL.,: "Simulation Calculation of Phantom Projections Based on Physics Model of Medical X-Ray Imaging", JOURNAL OF XI`AN JIAOTONG UNIVERSITY, vol. 40, no. 8, 31 August 2006 (2006-08-31), XP055602273 *

Also Published As

Publication number Publication date
CN105719338A (en) 2016-06-29

Similar Documents

Publication Publication Date Title
US8958527B2 (en) Method for acquiring morphology of a breast
US9947101B2 (en) Radiographic image analysis device and method, and recording medium having program recorded therein
US9554766B2 (en) Reducing artifacts in an image data set and X-ray device
EP3234919A1 (en) System and method for image reconstruction
CN105469366B (en) A kind of analytic method of abatement CT image metal artifacts
JP6313024B2 (en) Method and system for automatically determining a localizer within a scout image
JP6156847B2 (en) Radiation image processing apparatus and method, and program
CN111757701B (en) System and method for detecting anatomical regions
WO2017148115A1 (en) Phantom projection simulation method
US20220054862A1 (en) Medical image processing device, storage medium, medical device, and treatment system
RU2016132366A (en) Auto-calibration of a computer tomograph without calibration targets
JP2015514448A (en) Suppression of extracorporeal high-contrast objects based on stereo X-ray tubes
JP2018143699A5 (en)
JP2006167161A5 (en)
WO2020118894A1 (en) Image reconstruction method, computer device and computer readable storage medium
CN111179373A (en) Medical image bone removal model construction method and bone information removal method
CN110288540A (en) A kind of online imaging standards method of carbon-fibre wire radioscopic image
US20150063541A1 (en) X-ray image processing apparatus, x-ray diagnostic apparatus, and x-ray image processing method
WO2021044757A1 (en) Image processing device, image processing method, and program
TWI645836B (en) Particle beam therapy apparatus and digital reconstructed radiography image creation method
McSkimming et al. Adaptive kernel-based scatter correction for multi-source stationary CT with non-circular geometry
JP2013116213A (en) X-ray ct apparatus
KR102003197B1 (en) Scatter-correction system and method for a single radiographic projection
CN112446931A (en) Reconstruction data processing method and device, medical imaging system and storage medium
JP6299045B2 (en) Image processing apparatus, image processing method, and X-ray imaging apparatus

Legal Events

Date Code Title Description
NENP Non-entry into the national phase

Ref country code: DE

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

Ref document number: 16892300

Country of ref document: EP

Kind code of ref document: A1

122 Ep: pct application non-entry in european phase

Ref document number: 16892300

Country of ref document: EP

Kind code of ref document: A1