WO2017148115A1 - Phantom projection simulation method - Google Patents
Phantom projection simulation method Download PDFInfo
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three 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
本发明涉及一种模体投影仿真算法,属于医学图像处理技术领域。The invention relates to a phantom projection simulation algorithm and belongs to the technical field of medical image processing.
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)中投影结果的计算公式采用其中,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. 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.
下面对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。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射线能谱分布函数以及投影结果,本实施例中投影结果的计算公式采用其中,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 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)
- 一种模体投影仿真算法,其特征是,包括如下步骤: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所述的一种模体投影仿真算法,其特征是,所述步骤1)中,检测器采用虚拟XCT成像系统。A phantom projection simulation algorithm according to claim 1, wherein in the step 1), the detector employs a virtual XCT imaging system.
- 根据权利要求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.
- 根据权利要求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.
- 根据权利要求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.
- 根据权利要求1所述的一种模体投影仿真算法,其特征是,所述步骤3)中投影结果的计算公式采用其中,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 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.
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)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105719338A (en) * | 2016-02-29 | 2016-06-29 | 江苏美伦影像系统有限公司 | Motif projection simulation algorithm |
Citations (3)
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)
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) |
-
2016
- 2016-02-29 CN CN201610108658.3A patent/CN105719338A/en active Pending
- 2016-08-26 WO PCT/CN2016/096857 patent/WO2017148115A1/en active Application Filing
Patent Citations (3)
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)
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 |