KR101763323B1 - A reducing method of a missing wedge effect in 3D internal image by an oblique type CT system. - Google Patents

A reducing method of a missing wedge effect in 3D internal image by an oblique type CT system. Download PDF

Info

Publication number
KR101763323B1
KR101763323B1 KR1020160081979A KR20160081979A KR101763323B1 KR 101763323 B1 KR101763323 B1 KR 101763323B1 KR 1020160081979 A KR1020160081979 A KR 1020160081979A KR 20160081979 A KR20160081979 A KR 20160081979A KR 101763323 B1 KR101763323 B1 KR 101763323B1
Authority
KR
South Korea
Prior art keywords
wedge
missing wedge
angular range
missing
function
Prior art date
Application number
KR1020160081979A
Other languages
Korean (ko)
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 한국생산기술연구원
Priority to KR1020160081979A priority Critical patent/KR101763323B1/en
Application granted granted Critical
Publication of KR101763323B1 publication Critical patent/KR101763323B1/en

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/40—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis
    • A61B6/4021—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis involving movement of the focal spot
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computerised tomographs
    • A61B6/032—Transmission computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/00—3D [Three Dimensional] image rendering
    • 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
    • G06T17/10—Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10081—Computed x-ray tomography [CT]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20076—Probabilistic image processing

Abstract

The present invention can reduce line effects in a missing wedge angle region by solving the degradation of quality after a 3D internal shape image restoration caused by a missing wedge in an oblique type CT apparatus in which radiation and an object to be inspected are arranged at a specific angle. The present invention relates to a method for reducing a missing wedge influence in the 3D image processing of the oblique type CT apparatus, which includes the steps of: obtaining projection data for each rotation angle; restoring a 3D CT image; acquiring Fourier image data by performing Fourier transform; and minimizing total variation with regard to the Fourier image data, wherein the steps are repeatedly performed until a preset condition is satisfied.

Description

오블리크타입 CT장치의 3차원내부형상이미지처리 중의 미싱웨지영향 저감방법{A reducing method of a missing wedge effect in 3D internal image by an oblique type CT system.}Technical Field [0001] The present invention relates to a method of reducing the influence of a missing wedge during image processing of an internal shape image of an oblique type CT apparatus.

본 발명은, 오블리크타입 CT장치의 3차원내부형상이미지처리 중의 미싱웨지영향 저감방법에 대한 것으로, 더욱 상세하게는, 특히 방사선과 피검사체가 특정한 각도를 이루며 배치되는 오블리크타입의 CT장치를 대상으로 하여, 방사선이 통과하지 못하는 피검사체의 두께 등의 문제로부터 야기되는 특정 영역 뷰에 대한 데이터의 상실(미싱웨지)로부터 기인하는 3차원내부형상이미지복원 후의 단층영상의 품질저하 문제를 해결하여 미싱웨지 각도영역에서 라인효과를 개선할 수 있는 오블리크타입 CT장치의 3차원내부형상이미지처리 중의 미싱웨지영향 저감방법에 관한 것이다.The present invention relates to a method of reducing the influence of a missing wedge during three-dimensional internal shape image processing of an oblique type CT apparatus, and more particularly, to an oblique type CT apparatus in which radiation and an object to be inspected are arranged at a specific angle, Dimensional image after restoration of a three-dimensional internal shape image resulting from loss of data (missing wedge) for a specific region view caused by a problem such as thickness of an object that does not pass through the radiation, The present invention relates to a method of reducing the influence of a missing wedge during three-dimensional internal shape image processing of an oblique type CT apparatus capable of improving a line effect in a wedge angle region.

CT(Computed Tomography) 장치는, X선 등의 방사선의 출력원과 피검사체를 투과한 방사선의 수용부를 피검사체의 주위로 회전시키며, 수집된 데이터를 처리하여 피검사체에 대한 2, 3차원 이미지를 획득하는 장치이다. A CT (Computed Tomography) apparatus rotates a receiving portion of a radiation transmitted through an output source of a radiation such as an X-ray and an object to be examined to the periphery of the subject and processes the collected data to generate a two- or three- .

피검사체에 투과된 방사선의 양은 방사선의 진행 경로 상에 놓인 피검사체의 재료 및 진행거리에 따라 달라지며, 이러한 방사선의 양을 회전각도에 따라 그래프로 나타내면, 사이노그램(sinogram)을 얻을 수 있고, 이 사이노그램을 특정한 연산처리하여 피검사체에 대한 3차원영상을 얻을 수 있다.The amount of the radiation transmitted to the subject differs depending on the material and the traveling distance of the subject placed on the path of the radiation, and when the amount of such radiation is plotted according to the rotation angle, a sinogram can be obtained , And a three-dimensional image of the subject can be obtained by performing specific arithmetic processing on the sinogram.

이러한 CT장치는, 의료 분야에서도 많이 활용되지만, 정밀한 전자부품의 품질검사 등의 분야에서도 활용될 수 있다. 특히 최근에는, 1개의 패키치 내부에 복수 개의 칩을 수직으로 올려서 구성하는 3차원 패키징 기술이 널리 사용되고 있으며, CT장치를 이용하여, 복수 개의 칩간 전기 접속여부를 검사할 수 있다. Such a CT apparatus is widely used in the medical field, but can also be utilized in fields such as quality inspection of precision electronic parts. Particularly, in recent years, a three-dimensional packaging technique in which a plurality of chips are vertically raised inside a package is widely used, and it is possible to check whether a plurality of chips are electrically connected using a CT device.

대한민국 등록특허 10-1531654호(이하 인용문헌1이라 함)에서는, 투시 데이터로부터 압축센싱 기법을 이용하여 3차원적으로 재구성되는 CT(computed tomography) 볼륨 영상의 값 (i,j,k)를 구하기 위하여, 투시 데이터로부터 정확도를 강화시키기 위한 정확도 항(fidelity term)과, 희소성(sparsity)을 촉진시키는 규칙항(regularization term)을 포함하는 목적 함수를 구하는 단계, 경사 투시 알고리즘에 따라 상기 목적 함수를 에 대해 미분한 경사 함수를 구하는 단계 및 상기 경사 함수로부터 컨벡스 최소화를 위한 반복적 연산을 통하여 수렴되는 최종 결과값 (i,j,k)을 구하여 3차원 영상을 재구성하는 단계를 포함하는 3차원 반복적 영상 재구성 과정에서, 상기 목적 함수는 규칙항에 전체변이(total variation)와, 전체변이를 에 대해 미분한 경사 전체변이(gradient total variation)를 함께 적용하여 구해지는 것을 특징으로 하는 3차원 반복적 영상 재구성의 수렴속도 개선방법이 개시되어 있다.In Korean Patent No. 10-1531654 (hereinafter referred to as Reference 1), a value (i, j, k) of a reconstructed CT (volume) computed tomography volume image is obtained from perspective data using a compression sensing technique Obtaining an objective function including a fidelity term for enhancing the accuracy from the perspective data and a regularization term for promoting sparsity, calculating the objective function according to the oblique perspective algorithm, And reconstructing the three-dimensional image by obtaining a final result value (i, j, k) converged through an iterative operation for minimizing convex from the slope function, In the process, the objective function has a total variation in the rule term and a gradient total variation with respect to the total variation A method of improving convergence speed of a three-dimensional repetitive image reconstruction is disclosed.

구체적으로, 통상의 CT 시스템에서는 입사방사선 조사방향과 피검사체의 회전축이 상호 수직으로 구성되어 있기 때문에, 투영데이터는 피검사체 회전축에 대해 전체 360도 각도 모두에 대해 투영데이터가 획득되지 되지 못하고, 미싱웨지가 발생한다. 이러한 현상에 대응하기 위한 종래의 미싱웨지 처리 알고리즘은, 미싱웨지영역을 개선하는 것이 아니라, 복원된 단층영상을 전반적으로 개선하고자 하는 목적을 갖는다.Specifically, in the conventional CT system, since the incident radiation direction and the rotation axis of the subject are perpendicular to each other, the projection data can not obtain projection data for all the 360 degree angles with respect to the subject's rotation axis, A wedge occurs. Conventional sewing wedge processing algorithms for coping with such a phenomenon have an object to improve the restored tomographic image, not to improve the sewing wedge area.

KR 10-1531654 B1KR 10-1531654 B1

인용문헌1의 기술은, 피검사체에 대한 3차원단층영상의 처리 및 이러한 처리를 가속화할 수 있는 구성에 대해 개시하고 있을 뿐, CT 기술분야의 샘플링 단계에서 필연적으로 발생하는 미싱웨지의 존재로 인한 복원된 입체 단층 영상의 품질 저하를 해결하기 위한 기술구성을 개시하고 있지 않다는 제1문제점, 나아가 입사방사선의 방향과 피검사체의 회전축이 직각이 아닌 각도로 되는 오블리크타입의 CT장치에 있어 발생하는 기울어진 각도에 의한 미싱웨지에 대한 처리방안이 개시되고 있지 않는다는 제2문제점을 갖는다.The technology of Patent Document 1 discloses a process for processing a three-dimensional tomographic image of an object and a structure capable of accelerating the process, and it is possible to prevent the occurrence of a defect due to the presence of a missing wedge inevitably occurring in a sampling step in the field of CT technology The first problem is that it does not disclose a technique for solving the degradation of the reconstructed stereoscopic tomographic image, and further, there is a problem in the oblique type CT apparatus in which the direction of the incident radiation and the rotation axis of the test subject are not perpendicular to each other There is a second problem that a processing scheme for a missing wedge by an inclined angle is not disclosed.

또한, 종래 알려진 미싱웨지 처리 알고리즘들은, 빔과 물체회전축이 수직인 구성의 CT장치와 관련한 구성에 대하여 직접 적용이 가능할 뿐, 빔 방향과 피검사체의 회전축이 기울어진 구성의 CT장치에 대하여는 직접 적용하기 난해하다는 문제를 갖는다. In addition, the known sewing wedge processing algorithms can be directly applied to the configuration related to the CT apparatus having a configuration in which the beam and the object rotation axis are perpendicular to each other, and directly applied to the CT apparatus in which the beam direction and the rotation axis of the object are inclined There is a problem that it is difficult to do the following.

상기와 같은 문제를 해결하기 위해 안출된 본 발명은, 피검사체 상에 조사될 입사방사선을 형성하는 입사방사선형성부, 상기 피검사체에 상기 입사방사선이 투과되어 생성되는 투영데이터를 획득하는 투영데이터검출부 및 상기 투영데이터를 처리하여, 상기 피검사체의 3차원단층영상을 복원하는 3차원영상처리부를 포함하고, 상기 입사방사선의 조사방향과 상기 피검사체의 회전축이 소정의 각도를 이루도록 배치되는 오블리크타입 CT장치의 3차원영상처리 중의 미싱웨지(missing wedge)영향 저감방법으로서, (1) 상기 투영데이터검출부가, 상기 피검사체가 상기 회전축을 중심으로 회전함에 따라 생성되는 회전각도별 투영데이터를 획득하는 단계, (2) 상기 3차원영상처리부가 상기 소정의 각도 및 상기 회전각도별 투영데이터를 처리하여 상기 피검사체에 대한 3차원단층영상을 복원하는 단계, (3) 상기 3차원영상처리부가 상기 복원된 3차원단층영상에 대한 푸리에변환을 수행하여 푸리에영상데이터를 획득하는 단계, (4) 상기 3차원영상처리부가 상기 푸리에영상데이터에 대한 전체분산(total variation)을 최소화하는 단계, (5) 상기 (4)단계 이후의 상기 푸리에영상데이터를 푸리에역변환하여 출력영상으로 하는 단계를 포함하고, 상기 (2)단계 내지 상기 (4)단계는, 소정의 조건을 만족할 때까지 반복수행되는 구성으로 되는 것을 제공한다.According to an aspect of the present invention, there is provided an image processing apparatus including an incident radiation forming unit for forming incident radiation to be irradiated onto a subject, a projection data detecting unit for obtaining projection data generated by transmitting the incident radiation to the subject, And a three-dimensional image processing unit for processing the projection data and reconstructing a three-dimensional tomographic image of the subject, wherein the oblique type imaging unit is arranged so that the irradiation direction of the incident radiation and the rotation axis of the subject coincide at a predetermined angle A method of reducing a missing wedge effect in a three-dimensional image processing of a CT apparatus, the method comprising the steps of: (1) obtaining projection data for each rotation angle generated as the subject rotates about the rotation axis (2) The three-dimensional image processing section processes the projection data by the predetermined angle and the rotation angle, (3) performing a Fourier transform on the reconstructed 3-D tomographic image to obtain Fourier image data, (4) reconstructing the 3-dimensional tomographic image of the 3- Minimizing total variation of the Fourier image data; (5) performing Fourier inverse transform of the Fourier image data after the step (4) to produce an output image; The step (4) may be repeatedly performed until a predetermined condition is satisfied.

또한, 상기 (4)단계는, (4-1) 상기 전체분산과 소정의 기준값을 크기 비교하는 단계, (4-2) 상기 전체분산이 상기 소정의 기준값보다 작은 경우, 상기 3차원영상처리부가 상기 푸리에영상데이터를 출력영상의 기초로 제공하는 단계, (4-3) 상기 전체분산이 상기 소정의 기준값보다 큰 경우, 상기 3차원영상처리부가 미싱웨지보상을 수행하는 단계를 포함할 수 있다.The step (4) includes the steps of (4-1) comparing the total variance with a predetermined reference value, (4-2) if the total variance is smaller than the predetermined reference value, Providing the Fourier image data as a basis of an output image; and (4-3) performing the missing wedge compensation if the total variance is greater than the predetermined reference value.

또한, 상기 (4-3)단계는, (4-3-1) 상기 미싱웨지의 각도 범위를 입력받는 단계, (4-3-2) 상기 미싱웨지의 각도 범위 내의 소정의 각도값들에 대하여, 미싱웨지의 각도 범위 밖의 값들을 이용한 보간(interpolation)을 통해, 상기 미싱웨지의 각도 범위 내의 각도별 부분푸리에값을 추정하여, 상기 푸리에영상데이터를 업데이트하는 단계를 포함할 수 있다.The step (4-3) includes the steps of (4-3-1) receiving the angular range of the missing wedge, (4-3-2) selecting the angular range of the missing wedge And estimating a partial Fourier value per angle within the angular range of the missing wedge, through interpolation using values outside the angular range of the missing wedge, to update the Fourier image data.

또한, 상기 (4-3)단계는, (4-3-a) 상기 미싱웨지의 각도 범위를 입력받는 단계, (4-3-b) 상기 미싱웨지의 각도 범위 내의 각도별로, 푸리에변환값의 가중치를 계산하는 단계, (4-3-c) 상기 미싱웨지의 각도 범위 내의 각도별 상기 푸리에변환값의 가중치를 반영하여, 상기 미싱웨지의 각도 범위 내의 각도별 부분푸리에값을 추정하여 상기 푸리에영상데이터를 업데이트하는 단계를 포함할 수 있다.The step (4-3) includes the steps of (4-3-a) receiving the angular range of the missing wedge, (4-3-b) (4-3-c) estimating a partial Fourier value per angle in the angular range of the missing wedge, by reflecting the weight of the Fourier transform value for each angle within the angular range of the missing wedge, And updating the data.

또한, 본 발명은, 피검사체 상에 조사될 입사방사선을 형성하는 입사방사선형성부, 상기 피검사체에 상기 입사방사선이 투과되어 생성되는 투영데이터를 획득하는 투영데이터검출부 및 상기 투영데이터를 처리하여, 상기 피검사체의 3차원단층영상을 복원하는 3차원영상처리부를 포함하고, 상기 입사방사선의 조사방향과 상기 피검사체의 회전축이 소정의 각도를 이루도록 배치되는 오블리크타입 CT장치의 3차원영상처리 중의 미싱웨지(missing wedge)영향 저감방법으로서, (i) 상기 투영데이터검출부가, 상기 피검사체가 상기 회전축을 중심으로 회전함에 따라 생성되는 회전각도별 투영데이터를 획득하는 단계, (ii) 상기 3차원영상처리부가 상기 소정의 각도 및 상기 회전각도별 투영데이터를 처리하여 상기 피검사체에 대한 3차원단층영상을 복원하는 단계, (iii) 상기 3차원영상처리부가 상기 복원된 3차원단층영상을 전방투영(forward projection)하여, 연산투영데이터를 획득하는 단계, (iv) 상기 3차원영상처리부가 상기 연산투영데이터에 대하여, 역투영(back projection)하여, 수정3차원단층영상을 획득하고, 상기 3차원단층영상을 상기 수정3차원단층영상으로 업데이트하는 단계, (v) 상기 3차원영상처리부가 상기 수정3차원단층영상에 대한 전체분산(total variation)을 최소화하는 단계를 포함하고, 상기 (ii)단계 내지 상기 (v)단계는, 소정의 조건을 만족할 때까지 반복수행되는 것을 특징으로 하는 것을 제공한다.According to another aspect of the present invention, there is provided a projection exposure apparatus comprising: an incident radiation forming unit that forms incident radiation to be irradiated onto a subject; a projection data detecting unit that obtains projection data generated by transmitting the incident radiation to the subject; And a three-dimensional image processing unit for reconstructing a three-dimensional tomographic image of the subject, wherein the three-dimensional image processing unit of the oblique-type CT apparatus in which the irradiation direction of the incident radiation and the rotation axis of the subject are arranged at a predetermined angle A method of reducing a missing wedge effect, the method comprising: (i) obtaining projection data for each rotation angle generated as the subject rotates about the rotation axis; (ii) The image processing unit processes the projection data by the predetermined angle and the rotation angle to reconstruct a three-dimensional tomographic image of the subject (Iii) the three-dimensional image processing unit forward projection the reconstructed three-dimensional tomographic image to obtain operation projection data, (iv) the three-dimensional image processing unit Dimensional tomographic image to obtain a corrected three-dimensional tomographic image, and updating the three-dimensional tomographic image to the corrected three-dimensional tomographic image, (v) (Ii) to (v) are repeatedly performed until a predetermined condition is satisfied. The method according to claim 1,

또한, 본 발명은, 피검사체에 대한 조사방향이 상기 피검사체의 회전축과 소정의 각도를 이루도록 되는 입사방사선을 형성하는 입사방사선형성부, 상기 피검사체에 상기 입사방사선이 투과되어 생성되는 투영데이터를 획득하는 투영데이터검출부, 상기 투영데이터를 처리하여, 상기 피검사체의 3차원단층영상을 복원하는 3차원영상처리부를 포함하고, 상기 3차원영상처리부는, 3차원영상처리 중의 미싱웨지(missing wedge)영향 을 저감하기 위해, 상기 소정의 각도 및 회전각도별 투영데이터를 처리하여 상기 피검사체에 대한 3차원단층영상을 복원하는 제1기능, 상기 복원된 3차원단층영상에 대한 푸리에변환을 수행하여 푸리에영상데이터를 획득하는 제2기능, 상기 푸리에영상데이터에 대한 전체분산(total variation)을 최소화하는 제3기능을 순차적으로 수행하고, 상기 제1기능 내지 제3기능은, 소정의 조건을 만족할 때까지 반복수행되는 것을 특징으로 하는 미싱웨지 영향 저감이 가능한 오블리크타입 CT장치를 제공한다.According to another aspect of the present invention, there is provided a projection exposure apparatus, comprising: an incident radiation forming unit that forms incident radiation whose irradiation direction is such that the irradiation direction of the object is a predetermined angle with the rotation axis of the object to be inspected; Dimensional image processing unit for processing the projection data and reconstructing a three-dimensional tomographic image of the subject, wherein the three-dimensional image processing unit includes a missing wedge during a three-dimensional image processing, A first function for processing the projection data for each of the predetermined angles and rotation angles to reduce the influence to reconstruct a three-dimensional tomographic image of the subject, a second function for performing a Fourier transform on the reconstructed three- A second function of acquiring image data, and a third function of minimizing total variation of the Fourier image data, And, the first feature to the third feature, provides an oblique type CT apparatus capable of reducing the missing wedge effect, characterized in that is carried out repeatedly until they are satisfied a predetermined condition.

또한, 본 발명은, 피검사체에 대한 조사방향이 상기 피검사체의 회전축과 소정의 각도를 이루도록 되 입사방사선을 형성하는 입사방사선형성부, 상기 피검사체에 상기 입사방사선이 투과되어 생성되는 투영데이터를 획득하는 투영데이터검출부, 상기 투영데이터를 처리하여, 상기 피검사체의 3차원단층영상을 복원하는 3차원영상처리부를 포함하고, 상기 3차원영상처리부는, 3차원영상처리 중의 미싱웨지(missing wedge)영향 을 저감하기 위해, 상기 소정의 각도 및 회전각도별 투영데이터를 처리하여 상기 피검사체에 대한 3차원단층영상을 복원하는 제a기능, 상기 복원된 3차원단층영상을 전방투영(forward projection)하여, 연산투영데이터를 획득하는 제b기능, 상기 연산투영데이터에 대하여, 역투영(back projection)하여, 수정3차원단층영상을 획득하고, 상기 3차원단층영상을 상기 수정3차원단층영상으로 업데이트하는 제c기능, 및 상기 수정3차원단층영상에 대한 전체분산(total variation)을 최소화하는 제d기능을 순차적으로 수행하고, 상기 제a기능 내지 제d기능은, 소정의 조건을 만족할 때까지 반복수행할 수 있는 미싱웨지 영향 저감이 가능한 오블리크타입 CT장치를 제공한다.According to another aspect of the present invention, there is provided a projection exposure apparatus, comprising: an incident radiation forming unit for forming an incident radiation such that the irradiation direction of the object is a predetermined angle with the rotation axis of the object; Dimensional image processing unit for processing the projection data and reconstructing a three-dimensional tomographic image of the subject, wherein the three-dimensional image processing unit includes a missing wedge during a three-dimensional image processing, An a-function to reconstruct a three-dimensional tomographic image of the subject by processing the projection data for each of the predetermined angles and rotation angles in order to reduce the influence, a forward projection of the reconstructed 3-dimensional tomographic image A b-th function for obtaining operation projection data, back projection for the operation projection data to obtain a corrected three-dimensional tomographic image, A c-th function for updating the original tomographic image to the corrected 3-dimensional tomographic image, and a d-function for minimizing the total variation of the corrected 3-dimensional tomographic image, d function provides an oblique type CT apparatus capable of reducing the influence of a missing wedge that can be repeatedly performed until a predetermined condition is satisfied.

본 발명은, 방사선이 통과하지 못하는 피검사체의 두께 등의 문제로부터 야기되는 특정 영역 뷰에 대한 데이터의 상실(미싱웨지) 문제에 기인하는 3D복원 후의 단층영상의 품질저하 문제를 해결할 수 있고, 나아가, 방사선소스와 피검사체를 근접시킬 수 있어, 높은 해상도의 단층영상을 얻을 수 있는 장점을 가진 오블리크타입의 CT장치에 특화된 미싱웨지의 영향 최소화 방법을 개시하고 있다는 효과를 갖는다. 또한, 종래의 알고리즘이 미싱웨지영역을 특정적으로 개선하는 것이 아니라, 영상을 전반적으로 개선하는 특징(미싱웨지 해당되지 않는 나머지부분에 대해서도, 밝기의 저하가 함께 발생한다)이 있는 반면에, 본 발명은 미싱웨지영역에 대하여 집중적으로 처리(미싱웨지 해당되지 않는 나머지부분에 대해서는, 밝기의 저하가 거의 없다)할 수 있어, 미싱 각도에서 라인효과가 개선되므로, 최종적인 영상품질이 더 우수하다는 장점을 갖는다.The present invention solves the problem of degradation of the tomographic image after 3D reconstruction due to the loss of data (missing wedge) for a specific area view, which is caused by the problem of the thickness of the subject that does not pass through the radiation, , A method of minimizing the influence of a missing wedge that is specialized in an oblique type CT apparatus having the advantage that a radiation source and an object can be brought close to each other and a high resolution tomographic image can be obtained. In addition, the conventional algorithm does not specifically improve the missing wedge region, but has a feature of improving the image overall (a decrease in brightness is also caused in the remaining portion that does not correspond to the missing wedge) The invention has the advantage that the final image quality is better because the line effect can be improved intensively on the missing wedge area (there is little decrease in brightness for the remaining part that does not correspond to the missing wedge) Respectively.

도 1은, 본 발명의 오블리크타입의 CT장치의 일실시를 나타내는 모식도이다.
도 2는, 본 발명의 미싱웨지 영향 저감방법을 실시한 결과를 나타내는 그림이다.1 is a schematic diagram showing an embodiment of an oblique type CT apparatus of the present invention.
Fig. 2 is a diagram showing the result of performing the method of reducing the influence of the missing wedge of the present invention. Fig.

이하에서는 첨부한 도면을 참조하여 본 발명을 설명하기로 한다. 그러나 본 발명은 여러 가지 상이한 형태로 구현될 수 있으며, 따라서 여기에서 설명하는 실시예로 한정되는 것은 아니다. 그리고 도면에서 본 발명을 명확하게 설명하기 위해서 설명과 관계없는 부분은 생략하였으며, 명세서 전체를 통하여 유사한 부분에 대해서는 유사한 도면 부호를 붙였다.DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS Hereinafter, the present invention will be described with reference to the accompanying drawings. The present invention may, however, be embodied in many different forms and should not be construed as limited to the embodiments set forth herein. In order to clearly illustrate the present invention, parts not related to the description are omitted, and similar parts are denoted by like reference characters throughout the specification.

명세서 전체에서, 어떤 부분이 다른 부분과 "연결(접속, 접촉, 결합)"되어 있다고 할 때, 이는 "직접적으로 연결"되어 있는 경우뿐 아니라, 그 중간에 다른 부재를 사이에 두고 "간접적으로 연결"되어 있는 경우도 포함한다. 또한, 어떤 부분이 어떤 구성요소를 "포함"한다고 할 때, 이는 특별히 반대되는 기재가 없는 한다른 구성요소를 제외하는 것이 아니라 다른 구성요소를 더 구비할 수 있다는 것을 의미한다.Throughout the specification, when a part is referred to as being "connected" (connected, connected, coupled) with another part, it is not only the case where it is "directly connected" "Is included. Also, when a component is referred to as "comprising ", it means that it can include other components, not excluding components that are not specially contradictory.

본 명세서에서 사용한 용어는 단지 특정한 실시예를 설명하기 위해 사용된 것으로, 본 발명을 한정하려는 의도가 아니다. 단수의 표현은 문맥상 명백하게 다르게 뜻하지 않는 한, 복수개의 표현을 포함한다. 본 명세서에서, '포함하다' 또는 가지다' 등의 용어는 명세서상에 기재된 특징, 숫자, 단계, 동작, 구성요소, 부품 또는 이들을 조합한 것이 존재함을 지정하려는 것이지, 하나 또는 그 이상의 다른 특징들이나 숫자, 단계, 동작, 구성요소, 부품 또는 이들을 조합한 것들의 존재 또는 부가 가능성을 미리 배제하지 않는 것으로 이해되어야 한다. The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. The singular expressions include plural expressions unless the context clearly indicates otherwise. In this specification, the terms "comprises" or "having", etc. are intended to specify the presence of stated features, integers, steps, operations, elements, parts or combinations thereof, But do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, parts, or combinations thereof.

이하 첨부된 도면을 참고하여 본 발명의 실시예를 상세히 설명하기로 한다.Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings.

본 발명은, 오블리크타입의 CT장치에 직접 적용할 수 있는 미싱웨지 구역의 데이터 부족 문제를 해결하기 위한 기술이다. The present invention is a technique for solving the data shortage problem of a missing wedge zone that can be directly applied to an oblique type CT apparatus.

이에, 본 발명이 적용되는 CT장치는, 피검사체 상에 조사될 입사방사선을 형성하는 입사방사선형성부, 상기 피검사체에 상기 입사방사선이 투과되어 생성되는 투영데이터를 획득하는 투영데이터검출부 및 상기 투영데이터를 처리하여, 상기 피검사체의 3차원단층영상을 복원하는 3차원영상처리부의 필수구성요소를 구비하게 되며, 각 구성요소에 대해서는, 공지의 구성을 채택할 수 있음은 물론이다. A CT apparatus to which the present invention is applied includes an incident radiation forming unit for forming incident radiation to be irradiated on a subject, a projection data detecting unit for obtaining projection data generated by transmitting the incident radiation to the subject, Dimensional image processing unit for processing the data and reconstructing the three-dimensional tomographic image of the subject, and it is needless to say that a known configuration can be adopted for each component.

특히, 상기 입사방사선의 조사방향과 상기 피검사체의 회전축이 소정의 각도를 이루도록 배치될 수 있을 것이다. 이러한 구성에서는, 입사방사선형성부와 피검사체는 상호 근접한 위치에 배치될 수 있어, 나노 또는 마이크로 스케일의 피검사체 피검사영역에 대한 단층영상을 복원할 수 있다는 장점이 있게 된다. In particular, the irradiation direction of the incident radiation and the rotation axis of the subject may be arranged to form a predetermined angle. In such a configuration, the incident radiation forming unit and the subject can be disposed at positions close to each other, and thus it is possible to restore a tomographic image of the nano- or microscale inspection target area to be inspected.

오블리크 CT 시스템에서는 입사방사선의 조사방향과 피검사체의 회전축이 특정한 각도를 이루도록 구성되어 있고, 투영데이터는 피검사체 회전으로부터 획득되므로, 미싱웨지가 반드시 발생하게 된다. 즉, 입사방사선의 조사방향과 피검사체의 회전축이 이루는 특정한 각도(틸팅각)로 인해, 투영데이터검출부는, 투영데이터검출부에 의해 검출된 투영데이터들은 동일한 회전축 방향정보가 아닌 기울어진 축 정보임을 감안하여, 재구성되어야 하며, 이는 입사방사선의 조사방향과 피검사체의 회전축이 이루는 특정한 각도(틸팅각)에 대한 정보가 부족하다는 것을 의미한다.In the oblique CT system, the irradiation direction of the incident radiation and the rotation axis of the object to be inspected are configured to form a certain angle, and the projection data is obtained from the rotation of the object, so that the missing wedge is always generated. That is, due to the specific angle (tilting angle) between the irradiation direction of the incident radiation and the rotation axis of the subject, the projection data detection unit determines that the projection data detected by the projection data detection unit is not the same rotation axis direction information but the inclined axis information This means that there is a lack of information on a specific angle (tilting angle) between the irradiation direction of incident radiation and the rotation axis of the subject.

오블리크타입 CT장치에서는, 투영데이터검출부가 검출한 각도별 투영데이터마다 틸트된 회전축 방향 데이터로 재구성을 수행해야 하며, 틸팅각에 인한 데이터 손실로 미싱웨지가 발생하는 것이다. 이에, 단층영상의 복원(재구성)이후, 후처리 단계에서 미싱웨지를 개선하여야 한다. In the oblique type CT apparatus, reconstruction must be performed with tilted rotation axis direction data for each projection data of each angle detected by the projection data detection unit, and a missing wedge is generated due to data loss due to the tilting angle. Therefore, after restoration (reconstruction) of the tomographic image, the sewing wedge should be improved in the post-processing step.

일반 CT구성에서는, 복원이후 후처리과정에서, 전방투영(forward projection)과 역투영(back projection)의 작업을 통해, 전체분산(Total Variation)을 최소화 알고리즘을 적용하고, 이 때, 미싱웨지 영역의 데이터가 없다는 것을 전제로 알고리즘이 구성되는 것이지만, 오블리크타입 CT장치의 구성에서는 미싱웨지 영역 데이터가 중첩되어 투영데이터가 만들어지고, 결국, 단층영상이 되는 것이기 때문에, 미싱웨지 각도의 영향값을 변화(보상)하면서 전체분산(Total Variation)연산을 수행해야 한다.In the general CT configuration, a total variation minimization algorithm is applied through forward projection and back projection operations in post-restoration post-processing. At this time, However, in the configuration of the oblique type CT apparatus, since the missing data of the wedge region is superimposed to generate the projection data and eventually to the tomographic image, the influence value of the missing wedge angle is changed (Compensation), and perform a total variance calculation.

이하에서는, 미싱웨지에 대한 처리 관련하여, 총 두 가지의 실시예를 제안하는데, 하나는, 3D푸리에영역에서 알고리즘을 수행하는 것(제1실시예)이고, 다른 하나는, 2D푸리에영역과 3D푸리에영역을 전방투영과 역투영의 방법으로 전환하면서, 수행(제2실시예)하는 것이다.Hereinafter, with respect to the processing for the missing wedge, two total embodiments are proposed, one for performing the algorithm in the 3D Fourier domain (the first embodiment) and the other for the 2D Fourier domain and 3D (The second embodiment) while converting the Fourier domain into the forward projection and the reverse projection.

<< 제1실시예First Embodiment >>

이하, 오블리크타입 CT장치의 3차원영상처리 중의 미싱웨지(missing wedge)영향 저감방법 중 제1실시예에 대해 단계별로 설명하기로 한다.Hereinafter, a method of reducing the effect of missing wedge during 3D image processing of an oblique type CT apparatus will be described step by step.

(1) 상기 투영데이터검출부가, 상기 피검사체가 상기 회전축을 중심으로 회전함에 따라 생성되는 회전각도별 투영데이터를 획득한다. (1) The projection data detecting unit obtains projection data for each rotation angle generated as the subject rotates about the rotation axis.

(2) 상기 3차원영상처리부가 상기 소정의 각도(틸팅각) 및 상기 회전각도별 투영데이터를 처리하여 상기 피검사체에 대한 3차원단층영상을 복원한다. 오블리크타입의 CT 시스템에서는 미싱웨지현상으로 인한 단층영상의 품질저하를 복원 과정 이후에 확인할 수 있으므로, 복원된 3차원단층영상의 측면영역을 추출하여 각 추출단면에서 미싱웨지 개선 알고리즘을 적용하는 것이다.(2) The three-dimensional image processing unit processes the predetermined angle (tilting angle) and projection data for each rotation angle to reconstruct a three-dimensional tomographic image of the subject. In the oblique type CT system, since the degradation of the tomographic image due to the missing wedge phenomenon can be confirmed after the reconstruction process, the lateral area of the reconstructed 3-dimensional tomographic image is extracted and the missing wedge improvement algorithm is applied to each extracted section .

(3) 상기 3차원영상처리부가 상기 복원된 3차원단층영상에 대한 푸리에변환(3D푸리에변환)을 수행하여 푸리에영상데이터(3D푸리에변환된 3차원단층영상)를 획득한다. (3) The 3D image processing unit performs Fourier transformation (3D Fourier transform) on the reconstructed 3D tomographic image to obtain Fourier image data (3D Fourier transformed 3D tomographic image).

(4) 상기 3차원영상처리부가 상기 푸리에영상데이터에 대한 전체분산(total variation)을 최소화한다. 이런 과정을 통해, 미싱웨지에 대한 각도 영상효과를 감쇠하게 됨으로써, 미싱 각도에서의 라인효과가 개선될 수 있게 되는 것이다. (4) The 3D image processing unit minimizes the total variation of the Fourier image data. Through this process, the angle effect on the missing wedge is attenuated, so that the line effect at the missing angle can be improved.

(5) 이터레이션이 종료된 경우, 최종의 푸리에영상데이터를 푸리에역변환(3D푸리에역변환)하여 출력영상으로서, 최종 제공한다.(5) When the iteration ends, the final Fourier image data is Fourier inverse transformed (3D Fourier inverse transform) and finally provided as an output image.

이 중에서 전체분산의 최소화를 위한 (2) 내지 (4)단계는, 이터레이션의 방식으로 수행할 수 있으며, 이터레이션의 종료 조건의 일실시예로서는, 다음과 같은 것을 고려할 수 있지만, 이러한 것에 한정할 것은 아니다. Among these, steps (2) to (4) for minimizing the total dispersion can be performed by the iteration method. As one embodiment of the iteration end condition, the following can be considered, It is not.

구체적으로는, (4)단계를, (4-1) 상기 전체분산과 소정의 기준값을 크기 비교하는 단계, (4-2) 상기 전체분산이 상기 소정의 기준값보다 작은 경우, 상기 3차원영상처리부가 상기 푸리에영상데이터를 출력영상의 기초로 제공하는 단계(이터레이션의 종료), (4-3) 상기 전체분산이 상기 소정의 기준값보다 큰 경우, 상기 3차원영상처리부가 미싱웨지보상을 수행(이터레이션의 계속)하는 단계로 수행할 수 있다. Specifically, the step (4) may include: (4-1) comparing the total dispersion with a predetermined reference value; (4-2) when the total dispersion is smaller than the predetermined reference value, (4) if the total variance is greater than the predetermined reference value, the 3D image processing unit performs a missing wedge compensation (step (4)) of providing the Fourier image data on the basis of the output image The continuation of the iteration).

상기 (4)단계를 포함하는 방법은, 소위 TV(Total Variation)법으로서, 콘빔 CT 시스템에서 각도에 대한 보상방법으로 사용되는 것이다. 본 발명은 이러한 TV법을 오블리크타입의 CT 시스템에 적용하는 방법을 제안하는 것이다. The method including the step (4) is a so-called TV (Total Variation) method, which is used as a compensation method for an angle in a cone beam CT system. The present invention proposes a method of applying such a TV method to an oblique type CT system.

다만, 본 발명에서 적용되는 TV법이라 하여, 전체분산(total variation)만을 목적함수로 하여, 이를 최소화하는 것으로 한정할 필요는 없다. 목적함수는, 희소성과 관련되는 항 및 정확도와 관련되는 항들을 더 포함할 수 있으며, 특정한 형식에 한정되지 않는다는 것이다. 후술하는 실험예1에서 적용한 목적함수는 아래 수식1에 기재된 것과 같다. 여기서 u는 복원대상이며, 첫번째 항인 TV(u)는 전체분산을 정규화하는 항이고, 두번째 항은, 희소성관련 항이다. 또한, 세번째 항 중 Fp 는 부분푸리에행렬을, fp 는 부분 푸리에계수를 나타낸다.However, the TV method applied in the present invention does not need to limit the total variance to an objective function and to minimize it. The objective function may further include terms related to terms and accuracy associated with scarcity, and is not limited to any particular format. The objective function applied in Experimental Example 1 to be described later is as shown in Equation 1 below. Where u is the restoration object, TV (u), the first term, is the term that normalizes the total variance, and the second term is the term related to the scarcity. In the third term, F p denotes a partial Fourier matrix, f p Represents a partial Fourier coefficient.

<수식 1>&Lt; Formula 1 >

Figure 112016063178590-pat00001
Figure 112016063178590-pat00001

미싱웨지보상과 관련되는 (4-3)단계의 수행에 있어서, 구체적으로는, 복원된 3차원단층영상을 푸리에도메인으로 변환한 이후, 미싱웨지 각도에 대한 부분 푸리에값을 추정한 이후에, 원래 영상도메인으로 복원을 수행하는 것인데, 구체적인 실시예로 다음의 두 가지를 고려할 수 있다.In the performance of the step (4-3) related to the missing wedge compensation, specifically, after converting the reconstructed three-dimensional tomographic image into the Fourier domain, after estimating the partial Fourier value with respect to the missing wedge angle, The restoration is performed to the image domain, and the following two methods can be considered as concrete embodiments.

하나는, (4-3-1) 상기 미싱웨지의 각도 범위를 입력받고, (4-3-2) 상기 미싱웨지의 각도 범위 내의 소정의 각도값들에 대하여, 미싱웨지의 각도 범위 밖의 값들을 이용한 보간(interpolation)을 통해, 상기 미싱웨지의 각도 범위 내의 각도별 부분푸리에값을 추정하여, 상기 푸리에영상데이터(3D푸리에변환된 3차원단층영상)를 업데이트하는 것이다. 보간의 알고리즘은 공지의 것, 무엇이든 적용할 수 있다. (4-3-1) receives the angular range of the missing wedge, (4-3-2) obtains values outside the angular range of the missing wedge with respect to predetermined angular values within the angular range of the missing wedge (3D Fourier transformed three-dimensional tomographic image) by estimating a partial Fourier value per angle within the angular range of the missing wedge through interpolation using the Fourier transform. The algorithm of interpolation can be any known one.

둘은, (4-3)단계를, (4-3-a) 상기 미싱웨지의 각도 범위를 입력받는 단계 및(4-3), (4-3-a) receiving the angular range of the missing wedge, and

(4-3-b) 상기 미싱웨지의 각도 범위 내의 각도별로, 푸리에변환값의 가중치를 계산하는 단계, (4-3-c) 상기 미싱웨지의 각도 범위 내의 각도별 상기 푸리에변환값의 가중치를 반영하여, 상기 미싱웨지의 각도 범위 내의 각도별 부분푸리에값을 추정하여 상기 푸리에영상데이터(3D푸리에변환된 3차원단층영상)를 업데이트하는 단계로 구성하는 것이다.(4-3-b) calculating a weight of the Fourier transform value for each angle within the angular range of the missing wedge, (4-3-c) calculating a weight of the Fourier transform value for each angle within the angular range of the missing wedge And updating the Fourier image data (3D Fourier transformed three-dimensional tomographic image) by estimating a partial Fourier value per angle within the angular range of the missing wedge.

미싱웨지의 각도 범위는, 푸리에도메인에서의 샘플링 최대 기울임각도로 표시될 수 있다. 다만, 영상도메인(오브젝트도메인)에서의 각종 파라미터(일례로, 틸팅각 및 토모그래프를 위한 장치구성조건)를 입력받고, 푸리에도메인에서의 샘플링 최대 틸트각도를 자동 계산하는 방식을 배제하지는 않는다.The angular extent of the missing wedge may be expressed as the maximum sampling angle of inclination in the Fourier domain. However, it does not exclude a method of receiving various parameters (for example, device configuration conditions for tilting angle and tomograph) in the image domain (object domain) and automatically calculating the sampling maximum tilt angle in the Fourier domain.

미싱웨지의 각도 범위 내의 각도(정보가 보충되어야 하는 각도)는, 사전에 직접 개별적으로 결정되어 입력될 수도 있으나, 인접한 각도들 간의 차이(증분, increment)을 입력받아 결정될 수 있음은 물론이다. 이때의 증분은, 상호 상이하거나 동일하게 설정할 수 있다. 미묘한 각도 변경에 의해 투영데이터가 바뀌는 각도 구역에 있어서는, 증분을 작게 설정하고, 다소 각도가 많이 변경되더라도 투영데이터가 크게 바뀌지 않는 각도 구역에 있어서는, 증분을 크게 설정하는 등의 실시예를 고려할 수 있다.It is needless to say that the angles within the angular range of the missing wedge (the angles to which the information should be supplemented) may be directly determined and input directly in advance, but may be determined by inputting the difference between adjacent angles. The increments at this time can be set to be different or the same. It is possible to consider an embodiment in which the increment is set small in an angle area where projection data is changed by a subtle angle change and an increment is set large in an angle area in which the projection data does not largely change even if the angle is changed a little .

<실험예><Experimental Example>

제1실시예의 구성에 대한 미싱웨지영향 저감성능을 확인하기 위해 도2(a)에 도시된 것과 같은 영상을 준비하고, ± 60° 근처에서 미싱웨지 효과를 추가하여 준비하였다.An image as shown in Fig. 2 (a) was prepared and a missing wedge effect was added in the vicinity of 占 0 to prepare a sewing wedge effect reducing performance for the configuration of the first embodiment.

도2(b)는 본 발명의 기술을 적용하지 않은 비교예1에 대한 실험결과를 나타낸 것이다. 미싱웨지 효과에 의해서 바람직하지 않은 무늬 패턴이 포함되도록 복원된 결과를 보여준다.FIG. 2 (b) shows the results of an experiment for Comparative Example 1 to which the technique of the present invention is not applied. And the result is restored to include an undesired pattern due to the missing wedge effect.

도2(c)는, 고속푸리에변환영역에서 미싱된 각도 영역의 데이터를 주변 데이터로부터 정규화하여 전체분산(TV)값을 최소화 하도록 수행하는 것을 주요 특징으로 하는 본 발명의 실시예 1이 적용된 결과를 나타낸다. 제안한 알고리즘으로 미싱웨지 영향이 줄어든 결과를 보여준다. 2 (c) shows a result obtained by applying the first embodiment of the present invention, which is characterized in that the data of the angular region missing in the fast Fourier transform domain is normalized from the neighboring data to minimize the total dispersion (TV) value . The proposed algorithm shows the effect of missing wedge effect.

표1에는, 비교예1과 제1실시예에 대해 정량적인 성능 비교결과가 표시되어 있다. 특히, 제안한 알고리즘을 적용한 제1실시예에서는, 복원영상의 미싱웨지 효과를 줄이면서 원 샘플영상과의 correlation이 0.9869 로서 매우 높음을 확인할 수 있었다.Table 1 shows quantitative performance comparison results for Comparative Example 1 and the first embodiment. Particularly, in the first embodiment using the proposed algorithm, it is confirmed that the correlation with the original sample image is very high as 0.9869 while reducing the missing wedge effect of the reconstructed image.

비교
구분compare
division CorrelationCorrelation Relative Errors (%)Relative Errors (%) SNRSNR 비교예1Comparative Example 1 0.97940.9794 38.7338.73 7.27.2 제1실시예First Embodiment 0.98690.9869 21.1221.12 12.412.4

<< 제2실시예Second Embodiment >>

이하, 오블리크타입 CT장치의 3차원영상처리 중의 미싱웨지(missing wedge)영향 저감방법 중 제2실시예에 대해 단계별로 설명하기로 한다.Hereinafter, the second embodiment of the method for reducing the missing wedge effect during 3D image processing of the oblique type CT apparatus will be described step by step.

제1실시예의 그것과 장치 구성은 동일하게 할 수 있다.The configuration of the apparatus and the apparatus of the first embodiment can be the same.

(i) 상기 투영데이터검출부가, 상기 피검사체가 상기 회전축을 중심으로 회전함에 따라 생성되는 회전각도별 투영데이터를 획득한다. (i) The projection data detecting unit obtains projection data for each rotation angle generated as the subject rotates about the rotation axis.

(ii) 상기 3차원영상처리부가 상기 소정의 각도(틸팅각) 및 상기 회전각도별 투영데이터를 처리하여 상기 피검사체에 대한 3차원단층영상을 복원한다. 오블리크타입의 CT 시스템에서는 미싱웨지현상으로 인한 단층영상의 품질저하를 복원 과정 이후에 확인할 수 있으므로, 복원된 3차원단층영상의 측면영역을 추출하여 각 추출단면에서 미싱웨지 개선 알고리즘을 적용하는 것이다.(ii) The three-dimensional image processing unit processes the predetermined angle (tilting angle) and projection data for each rotation angle to reconstruct a three-dimensional tomographic image of the subject. In the oblique type CT system, since the degradation of the tomographic image due to the missing wedge phenomenon can be confirmed after the reconstruction process, the lateral area of the reconstructed 3-dimensional tomographic image is extracted and the missing wedge improvement algorithm is applied to each extracted section .

(iii) 상기 3차원영상처리부가 상기 복원된 3차원단층영상을 전방투영(forward projection)하여, 연산투영데이터를 획득한다. (iii) The 3D image processing unit forward-projects the reconstructed 3D tomographic image to obtain operation projection data.

(iv) 상기 3차원영상처리부가 상기 연산투영데이터에 대하여, 역투영(back projection)하여, 수정3차원단층영상을 획득하고, 상기 3차원단층영상을 상기 수정3차원단층영상으로 업데이트한다. (iv)단계 푸리에 도메인에서 수행될 수 있다.(iv) The 3-dimensional image processing unit back-projects the computed projection data to obtain a corrected 3-dimensional tomographic image, and updates the 3-dimensional tomographic image to the corrected 3-dimensional tomographic image. (iv) step Fourier domain.

(v) 상기 3차원영상처리부가 상기 수정3차원단층영상에 대한 전체분산(total variation)을 최소화한다. 이를 위해, (ii)단계 내지 (v)단계는, 소정의 조건을 만족할 때까지 반복수행되는 이터레이션의 방법으로 수행할 수 있다는 것은 제1실시예의 그것과 동일하다. (v)단계는 푸리에 도메인에서 수행될 수 있다.(v) the 3-dimensional image processing unit minimizes the total variation of the corrected 3-dimensional tomographic image. For this purpose, it is the same as that of the first embodiment that the steps (ii) and (v) can be performed by a repeated iteration method until the predetermined condition is satisfied. Step (v) may be performed in the Fourier domain.

이터레이션의 종료 조건을 포함하는 구체적인 방식에 관하여, (v)단계를 다음과 같이 진행한다.With respect to the concrete method including the termination condition of the iteration, the step (v) proceeds as follows.

(v-1) 상기 전체분산과 소정의 기준값을 크기 비교하고, (v-2) 상기 전체분산이 상기 소정의 기준값보다 작은 경우, 상기 3차원영상처리부가 상기 수정3차원단층영상을 출력영상으로 활용하며, (v-3) 상기 전체분산이 상기 소정의 기준값보다 큰 경우, 상기 3차원영상처리부가 미싱웨지보상을 수행한다. (v-1) comparing the total variance with a predetermined reference value, (v-2) if the total variance is smaller than the predetermined reference value, the 3D image processing unit converts the corrected 3D tomographic image into an output image (V-3) when the total variance is larger than the predetermined reference value, the 3D image processing unit performs missing wedge compensation.

상기 (v)단계를 포함하는 방법은, 소위 TV(Total Variation)법으로서, 콘빔 CT 시스템에서 각도에 대한 보상방법으로 사용되는 것이다. 본 발명은 이러한 TV법을 오블리크타입의 CT 시스템에 적용하는 방법을 제안하는 것이다. The method including the step (v) is a so-called TV (Total Variation) method, which is used as a compensation method for an angle in a cone beam CT system. The present invention proposes a method of applying such a TV method to an oblique type CT system.

다만, 본 발명에서 적용되는 TV법이라 하여, 전체분산(total variation)만을 목적함수로 하여, 이를 최소화하는 것으로 한정할 필요는 없다. 목적함수는, 희소성과 관련되는 항 및 정확도와 관련되는 항들을 더 포함할 수 있으며, 특정한 형식에 한정되지 않는다는 것이다. 전술한 실험예1에서 적용한 목적함수는 수식1에 기재된 것과 같다. 여기서 u는 복원대상이며, 첫번째 항인 TV(u)는 전체분산을 정규화하는 항이고, 두번째 항은, 희소성관련 항이다. 또한, 세번째 항 중 Fp 는 부분푸리에행렬을, fp 는 부분 푸리에계수를 나타낸다.However, the TV method applied in the present invention does not need to limit the total variance to an objective function and to minimize it. The objective function may further include terms related to terms and accuracy associated with scarcity, and is not limited to any particular format. The objective function applied in Experimental Example 1 described above is the same as described in Equation (1). Where u is the restoration object, TV (u), the first term, is the term that normalizes the total variance, and the second term is the term related to the scarcity. In the third term, F p represents a partial Fourier matrix, and f p represents a partial Fourier coefficient.

미싱웨지보상과 관련되는 (v-3)단계의 수행에 있어, 구체적인 실시예로 다음의 두 가지를 고려할 수 있다.In carrying out the step (v-3) related to the missing wedge compensation, the following two concrete examples can be considered.

하나는, (v-3-1) 상기 3차원영상처리부가 상기 수정3차원단층영상을 전방투영하여 수정연산투영데이터를 획득한다. 이 후, (v-3-2) 상기 미싱웨지의 각도 범위를 입력받고, (v-3-3) 상기 미싱웨지의 각도 범위 내의 소정의 각도값들에 대하여, 미싱웨지의 각도 범위 밖의 값들을 이용한 보간(interpolation)을 통해, 상기 (v-3-1)단계에서의 수정연산투영데이터를 업데이트한다. 보간의 알고리즘은 공지의 것, 무엇이든 적용할 수 있다. In (v-3-1), the 3D image processing unit front-projects the corrected 3D tomographic image to acquire correction operation projection data. Then, (v-3-2) receiving the angular range of the missing wedge, (v-3-3) calculating, for predetermined angular values within the angular range of the missing wedge, values outside the angular range of the missing wedge The correction operation projection data in the (v-3-1) step is updated through interpolation using the image data. The algorithm of interpolation can be any known one.

다른 하나는, 상기 (v-3)단계를 다음과 같이 구성하는 것이다.The other is to configure the step (v-3) as follows.

(v-3-a) 상기 3차원영상처리부가 상기 수정3차원단층영상을 전방투영하여 수정연산투영데이터를 획득한다. 이후, (v-3-b) 상기 미싱웨지의 각도 범위를 입력받고, (v-3-c) 상기 미싱웨지의 각도 범위 내의 각도별로, 푸리에변환값의 가중치를 계산하며, (v-3-d) 상기 미싱웨지의 각도 범위 내의 각도별 상기 푸리에변환값의 가중치를 반영하여, 상기 수정연산투영데이터를 업데이트한다.(v-3-a) The three-dimensional image processing unit forward-projects the corrected three-dimensional tomographic image to obtain corrected operation projection data. (V-3-c) calculating a weight of the Fourier transform value for each angle within the angular range of the missing wedge, (v-3-b) d) update the modified operation projection data by reflecting the weight of the Fourier transform value by angle within the angular range of the missing wedge.

두 경우 모두, 미싱웨지의 각도 범위는, 푸리에도메인에서의 샘플링 최대 틸트각도로 표시될 수 있다. 다만, 영상도메인(오브젝트도메인)에서의 각종 파라미터(일례로, 틸팅각 및 토모그래프를 위한 장치구성조건)를 입력받고, 푸리에도메인에서의 샘플링 최대 틸트각도를 자동 계산하는 방식을 배제하지는 않는다.In both cases, the angular extent of the missing wedge can be represented by the sampling maximum tilt angle in the Fourier domain. However, it does not exclude a method of receiving various parameters (for example, device configuration conditions for tilting angle and tomograph) in the image domain (object domain) and automatically calculating the sampling maximum tilt angle in the Fourier domain.

미싱웨지의 각도 범위 내의 각도(정보가 보충되어야 하는 각도)는, 사전에 직접 개별적으로 결정되어 입력될 수도 있으나, 인접한 각도들 간의 차이(증분, increment)을 입력받아 결정될 수 있음은 물론이다. 이때의 증분은, 상호 상이하거나 동일하게 설정할 수 있다. 미묘한 각도 변경에 의해 투영데이터가 바뀌는 각도 구역에 있어서는, 증분을 작게 설정하고, 다소 각도가 많이 변경되더라도 투영데이터가 크게 바뀌지 않는 각도 구역에 있어서는, 증분을 크게 설정하는 등의 실시예를 고려할 수 있다.It is needless to say that the angles within the angular range of the missing wedge (the angles to which the information should be supplemented) may be directly determined and input directly in advance, but may be determined by inputting the difference between adjacent angles. The increments at this time can be set to be different or the same. It is possible to consider an embodiment in which the increment is set small in an angle area where projection data is changed by a subtle angle change and an increment is set large in an angle area in which the projection data does not largely change even if the angle is changed a little .

또한, (iv)단계에서의 역투영(back projection)은, 가중역투영법(weighted back projection), 필터역투영법(filtered back projection), 가중필터역투영법(weighted filtered back projection) 중 하나를 적용할 수 있다. 바람직하게는, 고속 푸리에 변환(Fast Fourier Transform, FFT)을 이용한 WBP(Weighted Back Projection)방법을 사용하는데, 이 방법은 현재 가장 널리 쓰이는 방법으로 으로서, 라돈변환(radon transform) 형을 기본 수식으로 사용하고 있다. In addition, the back projection in the step (iv) can be applied to one of weighted back projection, filtered back projection, and weighted filtered back projection have. Preferably, a weighted back projection (WBP) method using Fast Fourier Transform (FFT) is used. This method is currently the most widely used method, and a radon transform type is used as a basic formula .

이하, 본 발명의 미싱웨지 영향 저감이 가능한 오블리크타입 CT장치에 대해 설명하고자 한다. 다만, 전술한 미싱웨지영향 저감방법과 실질적으로 동일한 것으로 인정되는 부분에 대해서는 상술을 생략하기로 한다.Hereinafter, an oblique type CT apparatus capable of reducing the influence of the missing wedge of the present invention will be described. However, the parts that are substantially equivalent to the above-described method of reducing the influence of the missing wedge will not be described in detail.

본 발명의 미싱웨지 영향 저감이 가능한 오블리크타입 CT장치에 대해서도 크게 두 가지 실시예가 존재한다. 하나는, 3D푸리에영역에서 알고리즘을 수행하는 것(실시예A)이고, 다른 하나는, 2D푸리에영역과 3D푸리에영역을 전방투영과 역투영의 방법으로 전환하면서, 수행(실시예B)하는 것이다.There are largely two embodiments of the oblique type CT apparatus capable of reducing the influence of the missing wedge of the present invention. One is to perform the algorithm in the 3D Fourier domain (embodiment A), and the other is to perform the 2D Fourier domain and the 3D Fourier domain with the forward projection and the reverse projection method (example B) .

<< 미싱웨지Sewing machine 영향  effect 저감이Reduction 가능한  possible 오블리크타입Oblique type CT장치CT device - - 실시예AExample A >>

미싱웨지 영향 저감이 가능한 오블리크타입 CT장치는, 피검사체에 대한 조사방향이 상기 피검사체의 회전축과 소정의 각도를 이루도록 되는 입사방사선을 형성하는 입사방사선형성부, 상기 피검사체에 상기 입사방사선이 투과되어 생성되는 투영데이터를 획득하는 투영데이터검출부, 상기 투영데이터를 처리하여, 상기 피검사체의 3차원단층영상을 복원하는 3차원영상처리부를 주요 구성요소로 가지면, 이와 관련하여서는, 공지의 구성을 채택하는 것이 가능하다.An oblique type CT apparatus capable of reducing the influence of a missing wedge includes an incident radiation forming unit that forms incident radiation whose irradiation direction with respect to the object to be inspected is at a predetermined angle with the rotation axis of the object to be inspected, Dimensional image processing unit for processing the projection data and reconstructing a three-dimensional tomographic image of the subject, as a main component, and in the related art, a known configuration It is possible to adopt.

특히 3차원영상처리부는, 본 발명의 핵심 구성으로서, 3차원영상처리 중의 미싱웨지(missing wedge)영향 을 저감하기 위해, 상기 소정의 각도 및 회전각도별 투영데이터를 처리하여 상기 피검사체에 대한 3차원단층영상을 복원하는 제1기능, 상기 복원된 3차원단층영상에 대한 푸리에변환을 수행하여 푸리에영상데이터를 획득하는 제2기능, 상기 푸리에영상데이터에 대한 전체분산(total variation)을 최소화하는 제3기능을 순차적으로 수행하고, 상기 제1기능 내지 제3기능은, 소정의 조건을 만족할 때까지 반복수행한다. Particularly, as a core configuration of the present invention, in order to reduce the influence of missing wedge during three-dimensional image processing, the three-dimensional image processing unit processes the projection data for the predetermined angle and rotation angle, A first function of reconstructing a 3D tomographic image, a second function of acquiring Fourier image data by performing a Fourier transform on the reconstructed 3D tomographic image, a second function of minimizing total variation of the Fourier image data, 3 functions sequentially, and the first function to the third function are repeatedly performed until a predetermined condition is satisfied.

상기 제3기능은, 상기 전체분산과 소정의 기준값을 크기 비교하는 제3-1기능,The third function may further include a third function for comparing magnitudes of the total variance and a predetermined reference value,

상기 전체분산이 상기 소정의 기준값보다 작은 경우, 상기 푸리에영상데이터를 출력영상의 기초로 제공하는 제3-2기능, 및 상기 전체분산이 상기 소정의 기준값보다 큰 경우, 미싱웨지보상을 수행하는 제3-3기능을 포함하도록 한다.A third function for providing the Fourier image data as a basis of an output image when the total variance is smaller than the predetermined reference value, and a third function for providing the Fourier image data as a basis of an output image, 3-3 Include functions.

나아가, 미싱웨지보상과 관련되는 (v-3)단계의 수행에 있어, 구체적인 실시예로 다음의 두 가지를 고려할 수 있다.Further, in carrying out the step (v-3) related to the missing wedge compensation, the following two concrete examples can be considered.

하나는, 상기 수정3차원단층영상을 전방투영하여 수정연산투영데이터를 획득하는 제3-3-1기능, 상기 미싱웨지의 각도 범위를 입력받는 제3-3-2기능, 및 상기 미싱웨지의 각도 범위 내의 소정의 각도값들에 대하여, 미싱웨지의 각도 범위 밖의 값들을 이용한 보간(interpolation)을 통해, 상기 푸리에영상데이터를 업데이트하는 제3-3-3기능을 포함하도록 할 수 있다.One is a 3-3-1 function for projecting the corrected 3D tomographic image forward to acquire correction operation projection data, a 3-3-2 function receiving an angular range of the missing wedge, and a 3- And a third 3-3 function for updating the Fourier image data by interpolation using values outside the angular range of the missing wedge, for predetermined angular values within the angular range.

다른 하나는, 상기 미싱웨지의 각도 범위를 입력받는 제3-3-a기능, 상기 미싱웨지의 각도 범위 내의 각도별로, 푸리에변환값의 가중치를 계산하는 제3-3-b기능, 및 상기 미싱웨지의 각도 범위 내의 각도별 상기 푸리에변환값의 가중치를 반영하여, 상기 푸리에영상데이터를 업데이트하는 제3-3-c기능을 포함하도록 할 수 있다.The other is a 3-3-a function receiving the angular range of the missing wedge, a 3-3-b function calculating the weight of the Fourier transform value for each angle within the angular range of the missing wedge, And a 3-3-c function for updating the Fourier image data by reflecting a weight of the Fourier transform value for each angle within the angular range of the wedge.

미싱웨지의 각도 범위는, 푸리에도메인에서의 샘플링 최대 기울임각도로 표시될 수 있다. 다만, 영상도메인(오브젝트도메인)에서의 각종 파라미터(일례로, 틸팅각 및 토모그래프를 위한 장치구성조건)를 입력받고, 푸리에도메인에서의 샘플링 최대 틸트각도를 자동 계산하는 방식을 배제하지는 않는다.The angular extent of the missing wedge may be expressed as the maximum sampling angle of inclination in the Fourier domain. However, it does not exclude a method of receiving various parameters (for example, device configuration conditions for tilting angle and tomograph) in the image domain (object domain) and automatically calculating the sampling maximum tilt angle in the Fourier domain.

미싱웨지의 각도 범위 내의 각도(정보가 보충되어야 하는 각도)는, 사전에 직접 개별적으로 결정되어 입력될 수도 있으나, 인접한 각도들 간의 차이(증분, increment)을 입력받아 결정될 수 있음은 물론이다. 이때의 증분은, 상호 상이하거나 동일하게 설정할 수 있다. 미묘한 각도 변경에 의해 투영데이터가 바뀌는 각도 구역에 있어서는, 증분을 작게 설정하고, 다소 각도가 많이 변경되더라도 투영데이터가 크게 바뀌지 않는 각도 구역에 있어서는, 증분을 크게 설정하는 등의 실시예를 고려할 수 있다.It is needless to say that the angles within the angular range of the missing wedge (the angles to which the information should be supplemented) may be directly determined and input directly in advance, but may be determined by inputting the difference between adjacent angles. The increments at this time can be set to be different or the same. It is possible to consider an embodiment in which the increment is set small in an angle area where projection data is changed by a subtle angle change and an increment is set large in an angle area in which the projection data does not largely change even if the angle is changed a little .

<< 미싱웨지Sewing machine 영향  effect 저감이Reduction 가능한  possible 오블리크타입Oblique type CT장치CT device - - 실시예BExample B >>

본 발명의 미싱웨지 영향 저감이 가능한 오블리크타입 CT장치는, 입사방사선형성부, 투영데이터검출부, 3차원영상처리부를 주요구성요소로 하고, 특히 3차원영상처리부는, 3차원영상처리 중의 미싱웨지(missing wedge)영향 을 저감하기 위해, 상기 소정의 각도 및 회전각도별 투영데이터를 처리하여 상기 피검사체에 대한 3차원단층영상을 복원하는 제a기능, 상기 복원된 3차원단층영상을 전방투영(forward projection)하여, 연산투영데이터를 획득하는 제b기능, 상기 연산투영데이터에 대하여, 역투영(back projection)하여, 수정3차원단층영상을 획득하고, 상기 3차원단층영상을 상기 수정3차원단층영상으로 업데이트하는 제c기능, 및 상기 수정3차원단층영상에 대한 전체분산(total variation)을 최소화하는 제d기능을 순차적으로 수행하며, 상기 제a기능 내지 제d기능은, 소정의 조건을 만족할 때까지 반복수행할 수 있다. The oblique type CT apparatus capable of reducing the influence of the missing wedge according to the present invention includes an incident radiation forming section, a projection data detecting section, and a three-dimensional image processing section as main components. Particularly, the three- (a) a function of processing the projection data by the predetermined angle and rotation angle to reduce a missing wedge effect to reconstruct a three-dimensional tomographic image of the subject, forward projection of the corrected three-dimensional tomographic image to acquire operation projection data; back projection of the computed projection data to obtain a corrected three-dimensional tomographic image; And a d-th function for minimizing a total variation of the corrected three-dimensional tomographic image, wherein the a-th function to the d-th function are sequentially performed, It can be repeated until you are satisfied with the defined conditions.

상기 제c기능 및 상기 제d기능은, 푸리에 도메인에서 수행할 수 있음은 물론이다.It is a matter of course that the c-th function and the d-th function can be performed in the Fourier domain.

또한, 제d기능은, 상기 전체분산과 소정의 기준값을 크기 비교하는 제d-1기능, 상기 전체분산이 상기 소정의 기준값보다 작은 경우, 상기 3차원영상처리부가 상기 수정3차원단층영상을 출력영상으로 활용하는 제d-2기능, 및 상기 전체분산이 상기 소정의 기준값보다 큰 경우, 상기 3차원영상처리부가 미싱웨지보상을 수행하는 제d-3기능을 포함할 수 있다.The d-th function may further include a d-1 function for comparing magnitudes of the total variance with a predetermined reference value, and a function for outputting the corrected three-dimensional tomographic image when the total variance is smaller than the predetermined reference value A d-2 function used as an image, and a d-3 function performing the missing wedge compensation if the total variance is greater than the predetermined reference value.

미싱웨지보상과 관련되는 (d-3)단계의 수행에 있어, 구체적인 실시예로 다음의 두 가지를 고려할 수 있다.In carrying out the step (d-3) related to the missing wedge compensation, the following two concrete examples can be considered.

하나는, 상기 수정3차원단층영상을 전방투영하여 수정연산투영데이터를 획득하는 제d-3-1기능, 상기 미싱웨지의 각도 범위를 입력받는 제d-3-2기능, 및 상기 미싱웨지의 각도 범위 내의 소정의 각도값들에 대하여, 미싱웨지의 각도 범위 밖의 값들을 이용한 보간(interpolation)을 통해, 상기 수정연산투영데이터를 업데이트하는 제d-3-3기능을 포함하도록 할 수 있다. 상기 미싱웨지의 각도 범위는, 푸리에도메인에서의 샘플링 최대 틸트각도로 표시될 수 있다.A d-3-1 function for obtaining the corrected operation projection data by projecting the corrected 3D tomographic image forward, a d-3-2 function receiving the angular range of the missing wedge, and a d- 3-3 function for updating the modified operation projection data by interpolation using values outside the angular range of the missing wedge, for predetermined angular values within the angular range. The angular extent of the missing wedge may be denoted by the sampling maximum tilt angle in the Fourier domain.

다른 하나는, 제d-3기능을, 상기 수정3차원단층영상을 전방투영하여 수정연산투영데이터를 획득하는 제d-3-a기능, 상기 미싱웨지의 각도 범위를 입력받는 제d-3-b기능, 상기 미싱웨지의 각도 범위 내의 각도별로, 푸리에변환값의 가중치를 계산하는 제d-3-c기능, 및 상기 미싱웨지의 각도 범위 내의 각도별 상기 푸리에변환값의 가중치를 반영하여, 상기 수정연산투영데이터를 업데이트하는 제d-3-d기능을 포함하도록 하는 것이다. 상기 미싱웨지의 각도 범위는, 푸리에도메인에서의 샘플링 최대 틸트각도로 표시될 수 있다.The other is a d-3-a function for projecting the corrected three-dimensional tomographic image forward by obtaining the corrected operation projection data, a d-3-a function for receiving the angular range of the missing wedge, b-function, a d-3-c function for calculating a weight of the Fourier transform value for each angle within the angular range of the missing wedge, and a weighting value of the Fourier transform value for each angle in the angular range of the missing wedge, And a d-3-d function for updating the correction operation projection data. The angular extent of the missing wedge may be denoted by the sampling maximum tilt angle in the Fourier domain.

미싱웨지의 각도 범위 내의 각도(정보가 보충되어야 하는 각도)는, 사전에 직접 개별적으로 결정되어 입력될 수도 있으나, 인접한 각도들 간의 차이(증분, increment)을 입력받아 결정될 수 있음은 물론이다. 이때의 증분은, 상호 상이하거나 동일하게 설정할 수 있다. 미묘한 각도 변경에 의해 투영데이터가 바뀌는 각도 구역에 있어서는, 증분을 작게 설정하고, 다소 각도가 많이 변경되더라도 투영데이터가 크게 바뀌지 않는 각도 구역에 있어서는, 증분을 크게 설정하는 등의 실시예를 고려할 수 있다.It is needless to say that the angles within the angular range of the missing wedge (the angles to which the information should be supplemented) may be directly determined and input directly in advance, but may be determined by inputting the difference between adjacent angles. The increments at this time can be set to be different or the same. It is possible to consider an embodiment in which the increment is set small in an angle area where projection data is changed by a subtle angle change and an increment is set large in an angle area in which the projection data does not largely change even if the angle is changed a little .

역투영(back projection)의 방식으로는, 가중역투영법(weighted back projection), 필터역투영법(filtered back projection), 가중필터역투영법(weighted filtered back projection) 중 하나를 적용할 수 있다.As a method of back projection, one of a weighted back projection, a filtered back projection, and a weighted filtered back projection can be applied.

본 발명을 첨부된 도면과 함께 설명하였으나, 이는 본 발명의 요지를 포함하는 다양한 실시 형태 중의 하나의 실시예에 불과하며, 당업계에서 통상의 지식을 가진 자가 용이하게 실시할 수 있도록 하는 데에 그 목적이 있는 것으로, 본 발명은 상기 설명된 실시예에만 국한되는 것이 아님은 명확하다. 따라서, 본 발명의 보호범위는 하기의 청구범위에 의해 해석되어야 하며, 본 발명의 요지를 벗어나지 않는 범위 내에서의 변경, 치환, 대체 등에 의해 그와 동등한 범위 내에 있는 모든 기술 사상은 본 발명의 권리범위에 포함될 것이다. 또한, 도면의 일부 구성은 구성을 보다 명확하게 설명하기 위한 것으로 실제보다 과장되거나 축소되어 제공된 것임을 명확히 한다.While the present invention has been particularly shown and described with reference to exemplary embodiments thereof, it should be understood that various changes and modifications will be apparent to those skilled in the art. Obviously, the invention is not limited to the embodiments described above. Accordingly, the scope of protection of the present invention should be construed according to the following claims, and all technical ideas which fall within the scope of equivalence by alteration, substitution, substitution, and the like within the scope of the present invention, Range. In addition, it should be clarified that some configurations of the drawings are intended to explain the configuration more clearly and are provided in an exaggerated or reduced size than the actual configuration.

100 : 입사방사선형성부
101 : 피검사체
102 : 투영데이터검출부
103 : 스테이지100: incident radiation forming part
101: Subject
102: projection data detector
103: stage

Claims (32)

피검사체 상에 조사될 입사방사선을 형성하는 입사방사선형성부, 상기 피검사체에 상기 입사방사선이 투과되어 생성되는 투영데이터를 획득하는 투영데이터검출부 및 상기 투영데이터를 처리하여, 상기 피검사체의 3차원단층영상을 복원하는 3차원영상처리부를 포함하고, 상기 입사방사선의 조사방향과 상기 피검사체의 회전축이 소정의 각도를 이루도록 배치되는 오블리크타입 CT장치의 3차원영상처리 중의 미싱웨지(missing wedge)영향 저감방법에 있어서,
(1) 상기 투영데이터검출부가, 상기 피검사체가 상기 회전축을 중심으로 회전함에 따라 생성되는 회전각도별 투영데이터를 획득하는 단계;
(2) 상기 3차원영상처리부가 상기 소정의 각도 및 상기 회전각도별 투영데이터를 처리하여 상기 피검사체에 대한 3차원단층영상을 복원하는 단계;
(3) 상기 3차원영상처리부가 상기 복원된 3차원단층영상에 대한 푸리에변환을 수행하여 푸리에영상데이터를 획득하는 단계;
(4) 상기 3차원영상처리부가 상기 푸리에영상데이터에 대한 전체분산(total variation)을 최소화하는 단계;
(5) 상기 (4)단계 이후의 상기 푸리에영상데이터를 푸리에역변환하여 출력영상으로 하는 단계;
를 포함하고, 상기 (2)단계 내지 상기 (4)단계는, 소정의 조건을 만족할 때까지 반복수행되는 것을 특징으로 하는 오블리크타입 CT장치의 3차원영상처리 중의 미싱웨지(missing wedge)영향 저감방법
A projection data detecting unit that obtains projection data generated by transmitting the incident radiation to the object to be inspected, and a projection data detecting unit that processes the projection data to detect a three-dimensional Dimensional image processing of an oblique-type CT apparatus having a three-dimensional image processing unit for reconstructing a tomographic image and arranged such that an irradiation direction of the incident radiation and a rotation axis of the object to be examined form a predetermined angle, a missing wedge, In an impact mitigation method,
(1) the projection data detecting unit obtains projection data for each rotation angle generated as the subject rotates around the rotation axis;
(2) reconstructing a three-dimensional tomographic image of the subject by processing the projection data according to the predetermined angle and the rotation angle of the three-dimensional image processing unit;
(3) obtaining the Fourier image data by performing the Fourier transform on the reconstructed three-dimensional tomographic image;
(4) minimizing total variation of the three-dimensional image processing unit with respect to the Fourier image data;
(5) transforming the Fourier image data after the step (4) into an output image by Fourier inversely transforming the Fourier image data;
(2) to (4) are repeatedly performed until a predetermined condition is satisfied. The method according to claim 1, wherein the oblique type CT apparatus is a three-dimensional image processing apparatus, Way
청구항 1에 있어서,
상기 (4)단계는,
(4-1) 상기 전체분산과 소정의 기준값을 크기 비교하는 단계;
(4-2) 상기 전체분산이 상기 소정의 기준값보다 작은 경우, 상기 3차원영상처리부가 상기 푸리에영상데이터를 출력영상의 기초로 제공하는 단계;
(4-3) 상기 전체분산이 상기 소정의 기준값보다 큰 경우, 상기 3차원영상처리부가 미싱웨지보상을 수행하는 단계;
를 포함하는 오블리크타입 CT장치의 3차원영상처리 중의 미싱웨지영향 저감방법.
The method according to claim 1,
The step (4)
(4-1) comparing the total variance with a predetermined reference value;
(4-2) if the total variance is smaller than the predetermined reference value, the 3D image processing unit provides the Fourier image data as a basis of the output image;
(4-3) if the total variance is greater than the predetermined reference value, performing the missing wedge compensation by the 3D image processing unit;
A method of reducing the impact of a missing wedge during a three-dimensional image processing of an oblique type CT apparatus including the method.
청구항 2에 있어서,
상기 (4-3)단계는,
(4-3-1) 상기 미싱웨지의 각도 범위를 입력받는 단계;
(4-3-2) 상기 미싱웨지의 각도 범위 내의 소정의 각도값들에 대하여, 미싱웨지의 각도 범위 밖의 값들을 이용한 보간(interpolation)을 통해, 상기 미싱웨지의 각도 범위 내의 각도별 부분푸리에값을 추정하여, 상기 푸리에영상데이터를 업데이트하는 단계;
를 포함하는 것을 특징으로 하는 오블리크타입 CT장치의 3차원영상처리 중의 미싱웨지영향 저감방법.
The method of claim 2,
The step (4-3)
(4-3-1) receiving an angular range of the missing wedge;
(4-3-2) For certain angular values within the angular range of the missing wedge, by interpolation using values outside the angular range of the missing wedge, a partial Fourier value per angle within the angular range of the missing wedge Updating the Fourier image data;
Wherein the oblique-type CT apparatus is a three-dimensional image processing apparatus.
청구항 2에 있어서,
상기 (4-3)단계는,
(4-3-a) 상기 미싱웨지의 각도 범위를 입력받는 단계;
(4-3-b) 상기 미싱웨지의 각도 범위 내의 각도별로, 푸리에변환값의 가중치를 계산하는 단계;
(4-3-c) 상기 미싱웨지의 각도 범위 내의 각도별 상기 푸리에변환값의 가중치를 반영하여, 상기 미싱웨지의 각도 범위 내의 각도별 부분푸리에값을 추정하여 상기 푸리에영상데이터를 업데이트하는 단계;
를 포함하는 것을 특징으로 하는 오블리크타입 CT장치의 3차원영상처리 중의 미싱웨지영향 저감방법.
The method of claim 2,
The step (4-3)
(4-3-a) receiving an angular range of the missing wedge;
(4-3-b) calculating weights of Fourier transform values for angles within an angular range of the missing wedge;
(4-3-c) updating the Fourier image data by estimating a partial Fourier value per angle within an angular range of the missing wedge, by reflecting a weight of the Fourier transform value for each angle within an angular range of the missing wedge;
Wherein the oblique-type CT apparatus is a three-dimensional image processing apparatus.
청구항 3 또는 청구항 4에 있어서,
상기 미싱웨지의 각도 범위 내의 각도는, 각도간 증분을 입력받아 결정되는 것을 특징으로 하는 오블리크타입 CT장치의 3차원영상처리 중의 미싱웨지영향 저감방법.

The method according to claim 3 or 4,
Wherein an angle within the angular range of the missing wedge is determined by inputting an increment between angles.

청구항 5에 있어서,
상기 각도간 증분은, 상호 상이하거나 동일한 것을 특징으로 하는 오블리크타입 CT장치의 3차원영상처리 중의 미싱웨지영향 저감방법.
The method of claim 5,
Wherein the increment between the angles is different or the same between the angles. The method of reducing the impact of a missing wedge during three-dimensional image processing in an oblique type CT apparatus.
청구항 3 또는 청구항 4에 있어서,
상기 미싱웨지의 각도 범위는,
푸리에도메인에서의 샘플링 최대 틸트각도로 표시되는 것을 특징으로 하는 오블리크타입 CT장치의 3차원영상처리 중의 미싱웨지영향 저감방법.
The method according to claim 3 or 4,
The angular range of the sewing wedge,
Wherein the sampling tilt angle in the Fourier domain is represented by a sampling maximum tilt angle.
피검사체 상에 조사될 입사방사선을 형성하는 입사방사선형성부, 상기 피검사체에 상기 입사방사선이 투과되어 생성되는 투영데이터를 획득하는 투영데이터검출부 및 상기 투영데이터를 처리하여, 상기 피검사체의 3차원단층영상을 복원하는 3차원영상처리부를 포함하고, 상기 입사방사선의 조사방향과 상기 피검사체의 회전축이 소정의 각도를 이루도록 배치되는 오블리크타입 CT장치의 3차원영상처리 중의 미싱웨지(missing wedge)영향 저감방법에 있어서,
(i) 상기 투영데이터검출부가, 상기 피검사체가 상기 회전축을 중심으로 회전함에 따라 생성되는 회전각도별 투영데이터를 획득하는 단계;
(ii) 상기 3차원영상처리부가 상기 소정의 각도 및 상기 회전각도별 투영데이터를 처리하여 상기 피검사체에 대한 3차원단층영상을 복원하는 단계;
(iii) 상기 3차원영상처리부가 상기 복원된 3차원단층영상을 전방투영(forward projection)하여, 연산투영데이터를 획득하는 단계;
(iv) 상기 3차원영상처리부가 상기 연산투영데이터에 대하여, 역투영(back projection)하여, 수정3차원단층영상을 획득하고, 상기 3차원단층영상을 상기 수정3차원단층영상으로 업데이트하는 단계;
(v) 상기 3차원영상처리부가 상기 수정3차원단층영상에 대한 전체분산(total variation)을 최소화하는 단계;
를 포함하고,
상기 (ii)단계 내지 상기 (v)단계는, 소정의 조건을 만족할 때까지 반복수행되는 것을 특징으로 하는 오블리크타입 CT장치의 3차원영상처리 중의 미싱웨지(missing wedge)영향 저감방법
A projection data detecting unit that obtains projection data generated by transmitting the incident radiation to the object to be inspected, and a projection data detecting unit that processes the projection data to detect a three-dimensional Dimensional image processing of an oblique-type CT apparatus having a three-dimensional image processing unit for reconstructing a tomographic image and arranged such that an irradiation direction of the incident radiation and a rotation axis of the object to be examined form a predetermined angle, a missing wedge, In an impact mitigation method,
(i) the projection data detecting unit obtains projection data for each rotation angle generated as the subject rotates about the rotation axis;
(ii) reconstructing a three-dimensional tomographic image of the subject by processing the projection data by the predetermined angle and the rotation angle of the three-dimensional image processing unit;
(iii) the three-dimensional image processing unit forward-projects the restored three-dimensional tomographic image to obtain operation projection data;
(iv) back-projecting the computed projection data to obtain a corrected three-dimensional tomographic image, and updating the three-dimensional tomographic image with the corrected three-dimensional tomographic image;
(v) minimizing total variation of the 3-dimensional tomographic image by the 3-dimensional image processing unit;
Lt; / RTI &gt;
Wherein the step (ii) and the step (v) are repeatedly performed until a predetermined condition is satisfied. The method for reducing the missing wedge effect in a three-dimensional image processing of an oblique type CT apparatus
청구항 8에 있어서,
상기 (v)단계는,
(v-1) 상기 전체분산과 소정의 기준값을 크기 비교하는 단계;
(v-2) 상기 전체분산이 상기 소정의 기준값보다 작은 경우, 상기 3차원영상처리부가 상기 수정3차원단층영상을 출력영상의 기초로 제공하는 단계;
(v-3) 상기 전체분산이 상기 소정의 기준값보다 큰 경우, 상기 3차원영상처리부가 미싱웨지보상을 수행하는 단계;
를 포함하는 오블리크타입 CT장치의 3차원영상처리 중의 미싱웨지영향 저감방법.

The method of claim 8,
The step (v)
(v-1) comparing the total variance with a predetermined reference value;
(v-2) if the total variance is smaller than the predetermined reference value, the 3D image processing unit providing the corrected 3D tomographic image as a basis of the output image;
(v-3) performing the missing wedge compensation if the total variance is greater than the predetermined reference value;
A method of reducing the impact of a missing wedge during a three-dimensional image processing of an oblique type CT apparatus including the method.

청구항 9에 있어서,
상기 (v-3)단계는,
(v-3-1) 상기 3차원영상처리부가 상기 수정3차원단층영상을 전방투영하여 수정연산투영데이터를 획득하는 단계;
(v-3-2) 상기 미싱웨지의 각도 범위를 입력받는 단계;
(v-3-3) 상기 미싱웨지의 각도 범위 내의 소정의 각도값들에 대하여, 미싱웨지의 각도 범위 밖의 값들을 이용한 보간(interpolation)을 통해, 상기 (v-3-1)단계에서의 수정연산투영데이터를 업데이트하는 단계;
를 포함하는 것을 특징으로 하는 오블리크타입 CT장치의 3차원영상처리 중의 미싱웨지영향 저감방법.
The method of claim 9,
The step (v-3)
(v-3-1) obtaining the corrected operation projection data by projecting the corrected three-dimensional tomographic image forward;
(v-3-2) receiving the angular range of the missing wedge;
(v-3-3) correcting for the predetermined angular values within the angular range of the missing wedge, by interpolation using values outside the angular range of the missing wedge, Updating the operation projection data;
Wherein the oblique-type CT apparatus is a three-dimensional image processing apparatus.
청구항 9에 있어서,
상기 (v-3)단계는,
(v-3-a) 상기 3차원영상처리부가 상기 수정3차원단층영상을 전방투영하여 수정연산투영데이터를 획득하는 단계;
(v-3-b) 상기 미싱웨지의 각도 범위를 입력받는 단계;
(v-3-c) 상기 미싱웨지의 각도 범위 내의 각도별로, 푸리에변환값의 가중치를 계산하는 단계;
(v-3-d) 상기 미싱웨지의 각도 범위 내의 각도별 상기 푸리에변환값의 가중치를 반영하여, 상기 수정연산투영데이터를 업데이트하는 단계;
를 포함하는 것을 특징으로 하는 오블리크타입 CT장치의 3차원영상처리 중의 미싱웨지영향 저감방법.
The method of claim 9,
The step (v-3)
(v-3-a) obtaining the corrected operation projection data by projecting the corrected three-dimensional tomographic image forward by the three-dimensional image processing unit;
(v-3-b) receiving the angular range of the missing wedge;
(v-3-c) calculating weights of Fourier transform values for angles within an angular range of the missing wedge;
(v-3-d) updating the corrected operation projection data by reflecting a weight of the Fourier transform value by angle within an angular range of the missing wedge;
Wherein the oblique-type CT apparatus is a three-dimensional image processing apparatus.
청구항 10 또는 청구항 11에 있어서,
상기 미싱웨지의 각도 범위 내의 각도는, 각도간 증분을 입력받아 결정되는 것을 특징으로 하는 오블리크타입 CT장치의 3차원영상처리 중의 미싱웨지영향 저감방법.
The method according to claim 10 or 11,
Wherein an angle within the angular range of the missing wedge is determined by inputting an increment between angles.
청구항 12에 있어서,
상기 각도간 증분은, 상호 상이하거나 동일한 것을 특징으로 하는 오블리크타입 CT장치의 3차원영상처리 중의 미싱웨지영향 저감방법.

The method of claim 12,
Wherein the increment between the angles is different or the same between the angles. The method of reducing the impact of a missing wedge during three-dimensional image processing in an oblique type CT apparatus.

청구항 8에 있어서,
상기 (iv)단계에서의 역투영(back projection)은, 가중역투영법(weighted back projection), 필터역투영법(filtered back projection), 가중필터역투영법(weighted filtered back projection) 중 하나를 적용하는 것을 특징으로 하는 오블리크타입 CT장치의 3차원영상처리 중의 미싱웨지영향 저감방법.
The method of claim 8,
The back projection in the step (iv) is performed by applying one of weighted back projection, filtered back projection, and weighted filtered back projection. A method of reducing the influence of a missing wedge during a three-dimensional image processing of an oblique type CT apparatus using the same.
청구항 8에 있어서,
상기 (iv)단계 및 상기 (v)단계는 푸리에 도메인에서 수행되는 것을 특징으로 하는 오블리크타입 CT장치의 3차원영상처리 중의 미싱웨지영향 저감방법.
The method of claim 8,
Wherein the step (iv) and the step (v) are performed in a Fourier domain.
청구항 10 또는 청구항 11에 있어서,
상기 미싱웨지의 각도 범위는,
푸리에도메인에서의 샘플링 최대 틸트각도로 표시되는 것을 특징으로 하는 오블리크타입 CT장치의 3차원영상처리 중의 미싱웨지영향 저감방법.
The method according to claim 10 or 11,
The angular range of the sewing wedge,
Wherein the sampling tilt angle in the Fourier domain is represented by a sampling maximum tilt angle.
피검사체에 대한 조사방향이 상기 피검사체의 회전축과 소정의 각도를 이루도록 되는 입사방사선을 형성하는 입사방사선형성부;
상기 피검사체에 상기 입사방사선이 투과되어 생성되는 투영데이터를 획득하는 투영데이터검출부;
상기 투영데이터를 처리하여, 상기 피검사체의 3차원단층영상을 복원하는 3차원영상처리부;
를 포함하고, 상기 3차원영상처리부는,
3차원영상처리 중의 미싱웨지(missing wedge)영향 을 저감하기 위해,
상기 소정의 각도 및 회전각도별 투영데이터를 처리하여 상기 피검사체에 대한 3차원단층영상을 복원하는 제1기능,
상기 복원된 3차원단층영상에 대한 푸리에변환을 수행하여 푸리에영상데이터를 획득하는 제2기능,
상기 푸리에영상데이터에 대한 전체분산(total variation)을 최소화하는 제3기능을 순차적으로 수행하고,
상기 제1기능 내지 제3기능은, 소정의 조건을 만족할 때까지 반복수행되는 것을 특징으로 하는 미싱웨지 영향 저감이 가능한 오블리크타입 CT장치.

An incident radiation forming unit for forming incident radiation whose irradiation direction with respect to the object to be inspected has a predetermined angle with the rotation axis of the object to be inspected;
A projection data detector for obtaining projection data generated by transmitting the incident radiation to the object;
A three-dimensional image processing unit for processing the projection data to reconstruct a three-dimensional tomographic image of the subject;
Wherein the three-dimensional image processing unit includes:
In order to reduce the effect of missing wedge during 3D image processing,
A first function of processing projection data for each of the predetermined angles and rotation angles to restore a three-dimensional tomographic image of the subject,
A second function of performing Fourier transform on the restored 3-D tomographic image to acquire Fourier image data,
Sequentially performing a third function for minimizing a total variation of the Fourier image data,
Wherein the first function to the third function are repeatedly performed until a predetermined condition is satisfied, wherein the oblique type CT apparatus is capable of reducing the influence of the missing wedge.

청구항 17에 있어서,
상기 제3기능은,
상기 전체분산과 소정의 기준값을 크기 비교하는 제3-1기능,
상기 전체분산이 상기 소정의 기준값보다 작은 경우, 상기 푸리에영상데이터를 출력영상의 기초로 제공하는 제3-2기능, 및
상기 전체분산이 상기 소정의 기준값보다 큰 경우, 미싱웨지보상을 수행하는 제3-3기능,
을 포함하는 미싱웨지 영향 저감이 가능한 오블리크타입 CT장치.
18. The method of claim 17,
The third function may include:
A 3-1 function of comparing the total variance with a predetermined reference value,
A 3-2 function of providing the Fourier image data as a basis of an output image when the total variance is smaller than the predetermined reference value,
A third 3-3 function for performing missing wedge compensation when the total variance is greater than the predetermined reference value,
An oblique type CT device capable of reducing the impact of the sewing wedge.
청구항 18에 있어서,
상기 제3-3기능은,
상기 미싱웨지의 각도 범위를 입력받는 제3-3-1기능, 및
상기 미싱웨지의 각도 범위 내의 소정의 각도값들에 대하여, 미싱웨지의 각도 범위 밖의 값들을 이용한 보간(interpolation)을 통해, 상기 푸리에영상데이터를 업데이트하는 제3-3-2기능,
을 포함하는 것을 특징으로 하는 미싱웨지 영향 저감이 가능한 오블리크타입 CT장치.
19. The method of claim 18,
The third 3-
A 3-3-1 function receiving the angular range of the missing wedge, and
A third 3-3 function for updating the Fourier image data by interpolation using values outside the angular range of the missing wedge, for predetermined angular values within the angular range of the missing wedge,
And an oblique type CT apparatus capable of reducing the influence of the missing wedge.
청구항 18에 있어서,
상기 제3-3기능은,
상기 미싱웨지의 각도 범위를 입력받는 제3-3-a기능,
상기 미싱웨지의 각도 범위 내의 각도별로, 푸리에변환값의 가중치를 계산하는 제3-3-b기능, 및
상기 미싱웨지의 각도 범위 내의 각도별 상기 푸리에변환값의 가중치를 반영하여, 상기 푸리에영상데이터를 업데이트하는 제3-3-c기능,
을 포함하는 것을 특징으로 하는 미싱웨지 영향 저감이 가능한 오블리크타입 CT장치.

19. The method of claim 18,
The third 3-
A 3-3-a function receiving the angular range of the missing wedge,
A 3-3-b function for calculating a weight of the Fourier transform value for each angle within the angular range of the missing wedge, and
A 3-3-c function for updating the Fourier image data by reflecting a weight of the Fourier transform value for each angle within an angular range of the missing wedge,
And an oblique type CT apparatus capable of reducing the influence of the missing wedge.

청구항 19 또는 청구항 20에 있어서,
상기 미싱웨지의 각도 범위 내의 각도는, 각도간 증분을 입력받아 결정되는 것을 특징으로 하는 미싱웨지 영향 저감이 가능한 오블리크타입 CT장치.
The method according to claim 19 or 20,
Wherein an angle within the angular range of the missing wedge is determined by inputting an increment between angles.
청구항 21에 있어서,
상기 각도간 증분은, 상호 상이하거나 동일한 것을 특징으로 하는 미싱웨지 영향 저감이 가능한 오블리크타입 CT장치.

23. The method of claim 21,
The oblique type CT apparatus capable of reducing the influence of the missing wedge is characterized in that the increment between the angles is different or the same.

청구항 19 또는 청구항 20에 있어서,
상기 미싱웨지의 각도 범위는,
푸리에도메인에서의 샘플링 최대 틸트각도로 표시되는 것을 특징으로 하는 미싱웨지 영향 저감이 가능한 오블리크타입 CT장치.
The method according to claim 19 or 20,
The angular range of the sewing wedge,
And a sampling tilt angle in a Fourier domain. The oblique type CT apparatus is capable of reducing the influence of a missing wedge.
피검사체에 대한 조사방향이 상기 피검사체의 회전축과 소정의 각도를 이루도록 되는 입사방사선을 형성하는 입사방사선형성부;
상기 피검사체에 상기 입사방사선이 투과되어 생성되는 투영데이터를 획득하는 투영데이터검출부;
상기 투영데이터를 처리하여, 상기 피검사체의 3차원단층영상을 복원하는 3차원영상처리부;
를 포함하고, 상기 3차원영상처리부는,
3차원영상처리 중의 미싱웨지(missing wedge)영향 을 저감하기 위해,
상기 소정의 각도 및 회전각도별 투영데이터를 처리하여 상기 피검사체에 대한 3차원단층영상을 복원하는 제a기능,
상기 복원된 3차원단층영상을 전방투영(forward projection)하여, 연산투영데이터를 획득하는 제b기능,
상기 연산투영데이터에 대하여, 역투영(back projection)하여, 수정3차원단층영상을 획득하고, 상기 3차원단층영상을 상기 수정3차원단층영상으로 업데이트하는 제c기능, 및
상기 수정3차원단층영상에 대한 전체분산(total variation)을 최소화하는 제d기능을 순차적으로 수행하고,
상기 제a기능 내지 제d기능은, 소정의 조건을 만족할 때까지 반복수행되는 것을 특징으로 하는 미싱웨지 영향 저감이 가능한 오블리크타입 CT장치.

An incident radiation forming unit for forming incident radiation whose irradiation direction with respect to the object to be inspected has a predetermined angle with the rotation axis of the object to be inspected;
A projection data detector for obtaining projection data generated by transmitting the incident radiation to the object;
A three-dimensional image processing unit for processing the projection data to reconstruct a three-dimensional tomographic image of the subject;
Wherein the three-dimensional image processing unit includes:
In order to reduce the effect of missing wedge during 3D image processing,
An a-function for processing the projection data for each of the predetermined angles and rotation angles to reconstruct a three-dimensional tomographic image of the subject,
A b-th function for forward-projecting the reconstructed three-dimensional tomographic image to obtain operation projection data,
A c-th function for performing back projection on the operation projection data to obtain a corrected three-dimensional tomographic image, and updating the three-dimensional tomographic image to the corrected three-dimensional tomographic image, and
Sequentially performing a d-function for minimizing a total variation of the corrected 3D tomographic image,
Wherein the a-th function and the d-th function are repeatedly performed until a predetermined condition is satisfied.

청구항 24에 있어서,
상기 제d기능은,
상기 전체분산과 소정의 기준값을 크기 비교하는 제d-1기능,
상기 전체분산이 상기 소정의 기준값보다 작은 경우, 상기 3차원영상처리부가 상기 수정3차원단층영상을 출력영상으로 활용하는 제d-2기능, 및
상기 전체분산이 상기 소정의 기준값보다 큰 경우, 상기 3차원영상처리부가 미싱웨지보상을 수행하는 제d-3기능,
을 포함하는 미싱웨지 영향 저감이 가능한 오블리크타입 CT장치.27. The method of claim 24,
The d-
A d-1 function for comparing the total variance with a predetermined reference value,
A d-2 function in which the 3D image processing unit uses the corrected 3D tomographic image as an output image when the total variance is smaller than the predetermined reference value, and
A d-3 function for performing the missing wedge compensation by the 3D image processing unit if the total variance is greater than the predetermined reference value,
An oblique type CT device capable of reducing the impact of the sewing wedge. 청구항 25에 있어서,
상기 제d-3기능은,
상기 수정3차원단층영상을 전방투영하여 수정연산투영데이터를 획득하는 제d-3-1기능,
상기 미싱웨지의 각도 범위를 입력받는 제d-3-2기능, 및
상기 미싱웨지의 각도 범위 내의 소정의 각도값들에 대하여, 미싱웨지의 각도 범위 밖의 값들을 이용한 보간(interpolation)을 통해, 상기 수정연산투영데이터를 업데이트하는 제d-3-3기능,
을 포함하는 것을 특징으로 하는 미싱웨지 영향 저감이 가능한 오블리크타입 CT장치.
26. The method of claim 25,
The (d-3)
A d-3-1 function of forwardly projecting the corrected three-dimensional tomographic image to acquire correction operation projection data,
A d-3-2 function receiving the angular range of the missing wedge, and
A d-3-3 function for updating the modified operation projection data through interpolation using values outside the angular range of the missing wedge, for predetermined angular values within the angular range of the missing wedge,
And an oblique type CT apparatus capable of reducing the influence of the missing wedge.
청구항 25에 있어서,
상기 제d-3기능은,
상기 수정3차원단층영상을 전방투영하여 수정연산투영데이터를 획득하는 제d-3-a기능,
상기 미싱웨지의 각도 범위를 입력받는 제d-3-b기능,
상기 미싱웨지의 각도 범위 내의 각도별로, 푸리에변환값의 가중치를 계산하는 제d-3-c기능, 및
상기 미싱웨지의 각도 범위 내의 각도별 상기 푸리에변환값의 가중치를 반영하여, 상기 수정연산투영데이터를 업데이트하는 제d-3-d기능,
을 포함하는 것을 특징으로 하는 미싱웨지 영향 저감이 가능한 오블리크타입 CT장치.

26. The method of claim 25,
The (d-3)
A d-3-a function for obtaining the corrected operation projection data by projecting the corrected three-dimensional tomographic image forward,
A d-3-b function receiving the angular range of the missing wedge,
A d-3-c function for calculating weights of Fourier transform values for angles within the angular range of the missing wedge, and
A d-3-d function for updating the corrected operation projection data by reflecting a weight of the Fourier transform value by angle within an angular range of the missing wedge,
And an oblique type CT apparatus capable of reducing the influence of the missing wedge.

청구항 26 또는 청구항 27에 있어서,
상기 미싱웨지의 각도 범위 내의 각도는, 각도간 증분을 입력받아 결정되는 것을 특징으로 하는 미싱웨지 영향 저감이 가능한 오블리크타입 CT장치.
26. The method of claim 26 or 27,
Wherein an angle within the angular range of the missing wedge is determined by inputting an increment between angles.
청구항 28에 있어서,
상기 각도간 증분은, 상호 상이하거나 동일한 것을 특징으로 하는 미싱웨지 영향 저감이 가능한 오블리크타입 CT장치.
29. The method of claim 28,
The oblique type CT apparatus capable of reducing the influence of the missing wedge is characterized in that the increment between the angles is different or the same.
청구항 24에 있어서,
상기 역투영(back projection)은, 가중역투영법(weighted back projection), 필터역투영법(filtered back projection), 가중필터역투영법(weighted filtered back projection) 중 하나를 적용하는 것을 특징으로 하는 미싱웨지 영향 저감이 가능한 오블리크타입 CT장치.
27. The method of claim 24,
Wherein the back projection applies one of a weighted back projection, a filtered back projection, and a weighted filtered back projection. This possible oblique type CT device.
청구항 24에 있어서,
상기 제c기능 및 상기 제d기능은, 푸리에 도메인에서 수행하는 것을 특징으로 하는 미싱웨지 영향 저감이 가능한 오블리크타입 CT장치.
27. The method of claim 24,
Wherein the c-th function and the d-th function are performed in a Fourier domain.
청구항 26 또는 청구항 27에 있어서,
상기 미싱웨지의 각도 범위는,
푸리에도메인에서의 샘플링 최대 틸트각도로 표시되는 것을 특징으로 하는 미싱웨지 영향 저감이 가능한 오블리크타입 CT장치.


26. The method of claim 26 or 27,
The angular range of the sewing wedge,
And a sampling tilt angle in a Fourier domain. The oblique type CT apparatus is capable of reducing the influence of a missing wedge.


KR1020160081979A 2016-06-29 2016-06-29 A reducing method of a missing wedge effect in 3D internal image by an oblique type CT system. KR101763323B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020160081979A KR101763323B1 (en) 2016-06-29 2016-06-29 A reducing method of a missing wedge effect in 3D internal image by an oblique type CT system.

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020160081979A KR101763323B1 (en) 2016-06-29 2016-06-29 A reducing method of a missing wedge effect in 3D internal image by an oblique type CT system.

Publications (1)

Publication Number Publication Date
KR101763323B1 true KR101763323B1 (en) 2017-08-01

Family

ID=59650140

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020160081979A KR101763323B1 (en) 2016-06-29 2016-06-29 A reducing method of a missing wedge effect in 3D internal image by an oblique type CT system.

Country Status (1)

Country Link
KR (1) KR101763323B1 (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000000236A (en) 1998-05-14 2000-01-07 General Electric Co <Ge> Artifact correction for high attenuation object
JP2005127886A (en) 2003-10-24 2005-05-19 Shimadzu Corp Inclined x-ray ct equipment and method for calibrating rotating center axis in the same
JP2011227028A (en) 2010-04-23 2011-11-10 Shimadzu Corp Tomogram image reconstruction method and x-ray ct apparatus

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000000236A (en) 1998-05-14 2000-01-07 General Electric Co <Ge> Artifact correction for high attenuation object
JP2005127886A (en) 2003-10-24 2005-05-19 Shimadzu Corp Inclined x-ray ct equipment and method for calibrating rotating center axis in the same
JP2011227028A (en) 2010-04-23 2011-11-10 Shimadzu Corp Tomogram image reconstruction method and x-ray ct apparatus

Similar Documents

Publication Publication Date Title
JP6301081B2 (en) X-ray computed tomography apparatus, image reconstruction method, and structure of reconstruction filter
US8731266B2 (en) Method and system for correcting artifacts in image reconstruction
EP2727083B1 (en) Iterative image reconstruction
US8774355B2 (en) Method and apparatus for direct reconstruction in tomosynthesis imaging
US7747057B2 (en) Methods and apparatus for BIS correction
CN109920020B (en) Cone beam CT (computed tomography) pathologic projection reconstruction artifact suppression method
JP6482934B2 (en) Image processing apparatus, radiation detection apparatus, and image processing method
WO2015010632A1 (en) X-ray dual energy ct reconstruction method
Abu Anas et al. Comparison of ring artifact removal methods using flat panel detector based CT images
US7656990B2 (en) Adaptive anisotropic filtering of projection data for computed tomography
Rashid et al. An improved method for the removal of ring artifacts in high resolution CT imaging
KR101591381B1 (en) Method for reducing metal artifact in computed tomography
Ametova et al. Software-based compensation of instrument misalignments for X-ray computed tomography dimensional metrology
CN110916708A (en) CT scanning projection data artifact correction method and CT image reconstruction method
US8098920B2 (en) Method of extrapolating a grey-level bidimensional image in a radiography device and a radiography device implementing such a method
CN111553960B (en) Ring artifact rapid correction method based on projection mean image
KR101763323B1 (en) A reducing method of a missing wedge effect in 3D internal image by an oblique type CT system.
US9613437B2 (en) Image producing apparatus, radiation tomography apparatus, and program
Yang et al. Fast iterative beam hardening correction based on frequency splitting in computed tomography
US11080900B2 (en) Method and apparatus for metal artifact reduction in industrial 3-dimensional cone beam computed tomography
KR20160056979A (en) Method and system of iterative image reconstruction for interior tomography using exterior projection data
Dennerlein et al. Geometric jitter compensation in cone-beam CT through registration of directly and indirectly filtered projections
Recur et al. Mojette reconstruction from noisy projections
KR102225043B1 (en) Method and apparatus for metal artifact reduction in industrial 3-dimensional cone beam computed tomography
CN112884862B (en) Cone beam CT temperature drift correction method and system based on centroid projection trajectory fitting

Legal Events

Date Code Title Description
E701 Decision to grant or registration of patent right
GRNT Written decision to grant