KR20170026984A - The beam-hardening correction method in X-ray CT - Google Patents

The beam-hardening correction method in X-ray CT Download PDF

Info

Publication number
KR20170026984A
KR20170026984A KR1020150123301A KR20150123301A KR20170026984A KR 20170026984 A KR20170026984 A KR 20170026984A KR 1020150123301 A KR1020150123301 A KR 1020150123301A KR 20150123301 A KR20150123301 A KR 20150123301A KR 20170026984 A KR20170026984 A KR 20170026984A
Authority
KR
South Korea
Prior art keywords
image
artifact
ray
artifact correction
computed tomography
Prior art date
Application number
KR1020150123301A
Other languages
Korean (ko)
Other versions
KR101717433B1 (en
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 KR1020150123301A priority Critical patent/KR101717433B1/en
Publication of KR20170026984A publication Critical patent/KR20170026984A/en
Application granted granted Critical
Publication of KR101717433B1 publication Critical patent/KR101717433B1/en

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5258Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computerised tomographs
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/56Details of data transmission or power supply, e.g. use of slip rings
    • A61B6/563Details of data transmission or power supply, e.g. use of slip rings involving image data transmission via a network

Abstract

The present invention relates to a beam-hardening artifact correction method in an X-ray CT, comprising: a step of dividing an area including high density material from an X-ray CT image; and a step of correcting an artifact of an image by using the area including the divided high density material and an artifact correction parameter. The artifact correction parameter can effectively reduce image distortion set by attenuation of the material and an X-ray spectrum.

Description

엑스레이 컴퓨터 단층촬영 환경에서 빔 경화현상에 의한 인공물 보정방법 {The beam-hardening correction method in X-ray CT}[0001] The present invention relates to a method for correcting an artifact caused by beam hardening in an X-ray computed tomography (X-ray CT)

본 발명은 엑스레이 컴퓨터 단층촬영 영상의 인공물 보정방법에 관한 것으로서, 더욱 상세하게는 고밀도 물질을 영역에 대한 인공물 보정 파라미터를 이용하여 영상의 인공물을 보정하는 방법 및 그 장치에 관한 것이다.The present invention relates to a method for correcting an artifact of an X-ray computed tomography image, and more particularly, to a method and apparatus for correcting an artifact of an image using an artifact correction parameter for a high-density material region.

현재 의료분야에서 필수적으로 사용되고 있는 CT의 가장 큰 문제점 중의 하나는 빔 경화(beam hardening) 현상에 의해 발생하는 흰색-검은색 줄무늬(streaking)로 나타나는 인공물이다. 빔 경화 현상은 낮은 에너지의 X-ray를 가지는 빔이 높은 에너지를 가지는 X-ray 빔에 비해 물체를 통과할 때 더 쉽게 흡수되는 현상을 말한다. 이로 인해 X-ray 데이터는 물체의 두께에 대해 선형적인 관계를 갖지 않는다. 빔 경화 현상은 저밀도 물질 (인체내 조직 등)에 비해 금속성 물질 (치아의 보철물, 임플란트, 인공관절 등) 과 같은 고밀도 물질에서 더욱 뚜렷하게 나타난다. 인체 내에 고밀도 물질이 존재할 경우, 선형관계를 가정으로 한 FBP 복원 알고리즘은 복원과정에서 심각한 인공물을 발생시켜 CT 영상을 왜곡시킨다. CT는 고품질의 의료영상을 제공하나, 빔 경화 현상에 의한 인공물은 임상적 유용성을 심각하게 떨어뜨린다.One of the biggest problems of CT, which is currently used in the medical field, is artifacts that appear as white-black streaking caused by beam hardening phenomenon. Beam curing is a phenomenon in which a beam with low energy X-rays is more easily absorbed when passing through an object than an X-ray beam with high energy. This causes the X-ray data to have no linear relationship to the thickness of the object. The beam hardening phenomenon is more pronounced in high density materials such as metallic materials (tooth prostheses, implants, artificial joints, etc.) compared to low density materials (such as tissues in the human body). When high density materials are present in the human body, the FBP reconstruction algorithm assuming a linear relationship generates serious artifacts during reconstruction, distorting CT images. CT provides high quality medical images, but artifacts due to beam hardening severely degrade clinical utility.

하드웨어 필터링(hardware filtering)은 종래 대표적인 빔 경화 보정방식중 하나로 X-ray 소스 앞에 알루미늄 얇은 판을 위치시켜 X-ray가 물체에 통과하기 전에 낮은 에너지의 빔을 사전에 차단하여 빔 경화 현상을 완화시키는 방식이다. 이 방식은 빔 경화 현상에 의한 인공물을 줄일 수 있지만, 필터링으로 인해 X-ray 선량이 감소하여 영상에 노이즈를 증가시켜 영상의 질을 떨어뜨리는 단점이 있다. Hardware filtering is one of the conventional beam curing correcting methods. The aluminum thin plate is placed in front of the X-ray source to prevent the beam hardening phenomenon by blocking the low energy beam before the X-ray passes through the object Method. Although this method can reduce artifacts due to beam hardening, the X-ray dose is reduced due to filtering, which results in an increase in noise in the image and a deterioration in image quality.

이중 에너지(dual energy)를 이용하여 빔 경화 현상을 보정하는 방법은 물질의 감쇄계수가 Compton scatter 와 Photoelectric interaction의 일차결합으로 표현된다는 가정 하에서 CT영상을 얻는다. 두 기저(basis)의 계수는 낮은 에너지와 높은 에너지를 가지는 두 개의 X-ray 빔을 이용하여 구할 수 있으며, 특정 에너지 레벨에서 복원된 물질 감쇄계수 (CT 영상)을 얻을 수 있다. 가상의 단일 에너지를 이용한 CT 영상을 얻었으므로 빔 경화 현상을 완화시킬 수 있다. 하지만, 이중에너지를 사용하여 X-ray 선량이 증가한다는 단점이 있다. The method of correcting beam hardening phenomenon using dual energy is obtained by assuming that the attenuation coefficient of material is represented by the primary combination of Compton scatter and Photoelectric interaction. Two basis coefficients can be obtained using two X-ray beams with low energy and high energy, and a restored material attenuation coefficient (CT image) at a specific energy level can be obtained. Since a CT image using a virtual single energy is obtained, the beam curing phenomenon can be mitigated. However, there is a disadvantage that the dose of X-ray is increased by using dual energy.

그 외 반복적 복원방법(iteration method), interpolation 또는 inpainting 등을 이용한 영상처리 방식이 있다. 이들 영상처리 방식은 모델의 복잡성에 의해 상당한 복원시간을 요구하거나, 인공물의 제거와 함께 다른 영상왜곡을 유발하기에 의사들의 부정적인 평가가 많다. Iteration method의 경우 일반적으로 인체 내에 사용되어지는 고물질의 정확한 구성성분을 사전정보로 이용하며 긴 복원시간 때문에 임상에 사용하기까지는 더 많은 연구가 필요하다.There are other image processing methods such as iterative method, interpolation or inpainting. These image processing methods require a considerable restoration time due to the complexity of the model, or cause other image distortions along with artifact removal. In the case of the Iteration method, it is generally necessary to use more precise components of the high-quality materials used in the human body as preliminary information, and further studies are required until clinical use because of the long restoration time.

한국공개특허공보 제10-2008-0010674호 "다항식모델 기반 엑스레이 영상의 밝기 및 형상왜곡보정방법"Korean Patent Application No. 10-2008-0010674 "Method of Correcting Brightness and Shape Distortion of X-ray Images Based on Polynomial Model"

본 발명이 해결하고자 하는 첫 번째 과제는 고밀도 물질을 영역에 대한 인공물 보정 파라미터를 이용하여 영상의 인공물을 보정하는 영상의 인공물 보정방법을 제공하는 것이다.The first problem to be solved by the present invention is to provide an artifact correction method of an image that corrects an artifact of an image using an artifact correction parameter for a high density material region.

본 발명이 해결하고자 하는 두 번째 과제는 고밀도 물질을 영역에 대한 인공물 보정 파라미터를 이용하여 영상의 인공물을 보정하는 영상의 인공물 보정장치를 제공하는 것이다.A second object of the present invention is to provide an image artifact correction apparatus for correcting an artifact of an image using an artifact correction parameter for a high density material region.

본 발명은 상기 첫 번째 과제를 달성하기 위하여, 엑스레이 컴퓨터 단층촬영 영상의 인공물 보정방법에 있어서, 엑스레이 컴퓨터 단층촬영 영상으로부터 고밀도 물질을 포함하는 영역을 분할하는 단계; 및 상기 분할된 고밀도 물질을 포함하는 영역과 인공물 보정 파라미터를 이용하여 영상의 인공물을 보정하는 단계를 포함하고, 상기 인공물 보정 파라미터는 물질의 감쇄계수 및 엑스레이 스펙트럼에 의해 설정되는 것을 특징으로 하는 방법을 제공한다.According to an aspect of the present invention, there is provided a method of correcting an artifact of an X-ray computed tomography image, the method comprising: dividing a region including a high density material from an X-ray computed tomography image; And correcting an artifact of an image using an artifact correction parameter and a region including the divided high density material, wherein the artifact correction parameter is set by an attenuation coefficient and an X-ray spectrum of the material. to provide.

본 발명의 실시예에 의하면, 상기 인공물 보정 파라미터는, 보정 영상에서 인공물을 최소화하도록 설정되는 것을 특징으로 하는 방법일 수 있다.According to an embodiment of the present invention, the artifact correction parameter may be set to minimize artefacts in the corrected image.

본 발명의 실시예에 의하면, 상기 인공물 보정 파라미터는, 인공물에 가중치를 부여하는 것을 특징으로 하는 방법일 수 있다.According to an embodiment of the present invention, the artifact correction parameter may be a method characterized by weighting an artifact.

본 발명은 상기 두 번째 과제를 달성하기 위하여, 엑스레이 컴퓨터 단층촬영 영상의 인공물 보정장치에 있어서, 엑스레이 컴퓨터 단층촬영 영상으로부터 고밀도 물질을 포함하는 영역을 분할하고, 상기 분할된 고밀도 물질을 포함하는 영역과 인공물 보정 파라미터를 이용하여 영상의 인공물을 보정하는 처리부를 포함하고, 상기 인공물 보정 파라미터는 물질의 감쇄계수 및 엑스레이 스펙트럼에 의해 설정되는 것을 특징으로 하는 장치를 제공한다.According to another aspect of the present invention, there is provided an apparatus for correcting an artifact of an X-ray computed tomography image, the apparatus comprising: an X-ray computed tomography apparatus for dividing a region including a high density material into a region including the divided high density material, And a processor for correcting the artifacts of the image using the artifact correction parameters, wherein the artifact correction parameters are set by the attenuation coefficient and the X-ray spectrum of the material.

본 발명의 실시예에 의하면, 엑스레이 컴퓨터 단층촬영 영상 정보를 수신하는 통신부를 더 포함하는 장치일 수 있다.According to an embodiment of the present invention, the apparatus may further include a communication unit for receiving the X-ray computed tomographic image information.

본 발명에 따르면, 효율적으로 영상왜곡을 줄일 수 있다. 또한, 컴퓨터 단층촬영 영상에 바로 적용이 가능하여 인공물에 의한 영상왜곡을 간편하게 제거할 수 있다.According to the present invention, image distortion can be efficiently reduced. In addition, it can be applied directly to a computed tomography image, so that image distortion due to artifacts can be easily removed.

도 1은 본 발명의 일 실시예에 따른 엑스레이 컴퓨터 단층촬영 영상의 인공물 보정장치의 블록도이다.
도 2는 엑스레이 컴퓨터 단층촬영 영상에서 발생하는 인공물이다.
도 3 내지 4는 본 발명의 실시예에 따른 엑스레이 컴퓨터 단층촬영 영상의 인공물 보정과정을 나타낸 것이다.
도 5는 본 발명의 실시예에 따른 엑스레이 컴퓨터 단층촬영 영상의 인공물 보정방법의 흐름도이다.
1 is a block diagram of an artifact correction apparatus for an X-ray computed tomography image according to an embodiment of the present invention.
2 is an artifact arising from an X-ray computed tomography image.
3 to 4 illustrate an artifact correction process of an X-ray computed tomography image according to an embodiment of the present invention.
5 is a flowchart illustrating an artifact correction method of an X-ray computed tomography image according to an embodiment of the present invention.

본 발명에 관한 구체적인 내용의 설명에 앞서 이해의 편의를 위해 본 발명이 해결하고자 하는 과제의 해결 방안의 개요 혹은 기술적 사상의 핵심을 우선 제시한다.Prior to the description of the concrete contents of the present invention, for the sake of understanding, the outline of the solution of the problem to be solved by the present invention or the core of the technical idea is first given.

본 발명의 일 실시예에 따른 엑스레이 컴퓨터 단층촬영 영상의 인공물 보정방법은 엑스레이 컴퓨터 단층촬영 영상으로부터 고밀도 물질을 포함하는 영역을 분할하는 단계, 및 상기 분할된 고밀도 물질을 포함하는 영역과 인공물 보정 파라미터를 이용하여 영상의 인공물을 보정하는 단계를 포함하고, 상기 인공물 보정 파라미터는 물질의 감쇄계수 및 엑스레이 스펙트럼에 의해 설정되는 것을 특징으로 한다.An artifact correction method for an X-ray computed tomography image according to an exemplary embodiment of the present invention includes dividing a region including a high density material from an X-ray computed tomography image, and dividing an area including the divided high density material and an artifact correction parameter Wherein the artifact correction parameter is set by an attenuation coefficient and an X-ray spectrum of the material.

이하 첨부된 도면을 참조하여 본 발명이 속하는 기술 분야에서 통상의 지식을 가진 자가 본 발명을 용이하게 실시할 수 있는 실시 예를 상세히 설명한다. 그러나 이들 실시예는 본 발명을 보다 구체적으로 설명하기 위한 것으로, 본 발명의 범위가 이에 의하여 제한되지 않는다는 것은 당업계의 통상의 지식을 가진 자에게 자명할 것이다.DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings. It will be apparent to those skilled in the art, however, that these examples are provided to further illustrate the present invention, and the scope of the present invention is not limited thereto.

본 발명이 해결하고자 하는 과제의 해결 방안을 명확하게 하기 위한 발명의 구성을 본 발명의 바람직한 실시예에 근거하여 첨부 도면을 참조하여 상세히 설명하되, 도면의 구성요소들에 참조번호를 부여함에 있어서 동일 구성요소에 대해서는 비록 다른 도면상에 있더라도 동일 참조번호를 부여하였으며 당해 도면에 대한 설명시 필요한 경우 다른 도면의 구성요소를 인용할 수 있음을 미리 밝혀둔다. 아울러 본 발명의 바람직한 실시 예에 대한 동작 원리를 상세하게 설명함에 있어 본 발명과 관련된 공지 기능 혹은 구성에 대한 구체적인 설명 그리고 그 이외의 제반 사항이 본 발명의 요지를 불필요하게 흐릴 수 있다고 판단되는 경우, 그 상세한 설명을 생략한다.
BRIEF DESCRIPTION OF THE DRAWINGS The above and other features and advantages of the present invention will become more apparent from the following detailed description of the present invention when taken in conjunction with the accompanying drawings, It is to be noted that components are denoted by the same reference numerals even though they are shown in different drawings, and components of different drawings can be cited when necessary in describing the drawings. In the following detailed description of the principles of operation of the preferred embodiments of the present invention, it is to be understood that the present invention is not limited to the details of the known functions and configurations, and other matters may be unnecessarily obscured, A detailed description thereof will be omitted.

도 1은 본 발명의 일 실시예에 따른 엑스레이 컴퓨터 단층촬영 영상의 인공물 보정장치의 블록도이다.1 is a block diagram of an artifact correction apparatus for an X-ray computed tomography image according to an embodiment of the present invention.

본 발명의 일 실시예에 따른 엑스레이 컴퓨터 단층촬영 영상의 인공물 보정장치(100)는 처리부(110) 및 통신부(120)로 구성된다. 저장부를 더 포함할 수 있다.An artifact correction apparatus 100 for an X-ray computed tomography image according to an embodiment of the present invention includes a processing unit 110 and a communication unit 120. And may further include a storage unit.

통신부(120)는 엑스레이 컴퓨터 단층촬영 영상정보를 수신한다. 또한, 영상의 인공물이 보정된 영상을 사용자에게 제공할 수 있다.The communication unit 120 receives the X-ray computed tomography image information. Further, it is possible to provide the user with an image in which the artifact of the image is corrected.

처리부(110)는 엑스레이 컴퓨터 단층촬영 영상으로부터 고밀도 물질을 포함하는 영역을 분할하고, 상기 분할된 고밀도 물질을 포함하는 영역과 인공물 보정 파라미터를 이용하여 영상의 인공물을 보정하는 처리부를 포함한다.The processing unit 110 includes a processing unit for dividing an area including a high density material from an X-ray computed tomography image, and correcting an artifact of the image using an area including the divided high density material and an artifact correction parameter.

보다 구체적으로, 영상의 인공물을 보정함에 있어서, 분할된 고밀도 영역과 인공물 보정 파라미터에 대응하는 인공물 보정함수인 빔 경화 보정 함수를 도출한다. 빔 경화 보정 함수를 엑스레이 컴퓨터 단층촬영 영상에 적용하여 인공물에 의한 영상 왜곡이 없는 영상을 생성할 수 있다.More specifically, in correcting an artifact of an image, a beam hardening correction function which is an artifact correction function corresponding to the divided high density area and the artifact correction parameter is derived. By applying the beam hardening correction function to the X - ray computed tomography images, images without image distortion due to artifacts can be generated.

엑스레이 컴퓨터 단층촬영을 하고자 하는 대상에 도 2(a)와 같이 고밀도 물질이 포함되어 있는 경우, 복원된 엑스레이 컴퓨터 단층촬영 영상에는 도 2(b)와 같이, 영상왜곡이 발생한다. 즉, 고밀도 물질의 기하학적인 구조에 의해 streaking artifact가 발생한다.When a subject who intends to perform an X-ray computed tomography includes a high-density material as shown in FIG. 2 (a), image distortion occurs in the reconstructed X-ray computed tomography image as shown in FIG. 2 (b). That is, streaking artifacts occur due to the geometric structure of high-density materials.

상기 영상의 인공물을 제거하기 위하여, 도 3과 같은 과정을 수행한다. 컴퓨터 단층촬영 영상(fCT)에서 인공물에 의한 영상(φD,λ)을 제거함으로써 영상의 인공물이 보정된 영상(fCTD,λ)을 얻을 수 있다.In order to remove artifacts of the image, a process as shown in FIG. 3 is performed. By removing the image (D φ, λ) due to artefacts in the CT image (CT f) it is possible to obtain the image of the artifact-corrected image (f CT -φ D, λ) .

본 발명에 일 실시예에 따른 영상의 인공물을 보정은, 이중 에너지, 물질의 종류, CT 시스템에 사용되는 에너지 스펙트럼 등의 사전 정보나 X-ray 데이터(sinogram) 등의 이용없이 CT 영상만을 이용하여 CT 영상복원방식인 FBP에 의해 제공되는 CT 영상에서 나타나는 인공물에 의한 영상 왜곡에 바로 적용이 가능하다.The correction of artifacts of an image according to an embodiment of the present invention can be performed by using only CT images without using prior information such as dual energy, kind of material, energy spectrum used in CT system, X-ray data (sinogram) It can be directly applied to image distortion caused by artifacts in the CT image provided by the FBP, which is a CT image restoration method.

fCT를 FBP로 복원한 CT 영상이라 할때, 다음과 같이 나타낼 수 있다.f CT is the reconstructed CT image, it can be expressed as follows.

Figure pat00001
Figure pat00001

여기서, X-ray 데이터

Figure pat00002
는 다음과 같이 정의된다.Here, the X-ray data
Figure pat00002
Is defined as follows.

Figure pat00003
Figure pat00003

상기 수학식 1,2에 사용된 기호는 다음과 같다.The symbols used in the above equations (1) and (2) are as follows.

1.

Figure pat00004
:위치
Figure pat00005
, 에너지 E에서의 감쇄계수One.
Figure pat00004
:location
Figure pat00005
, The attenuation coefficient at energy E

2.

Figure pat00006
:X-ray 에너지 스펙트럼2.
Figure pat00006
: X-ray energy spectrum

3.

Figure pat00007
,
Figure pat00008
:라돈변환(Radon transform)3.
Figure pat00007
,
Figure pat00008
: Radon transform

4.

Figure pat00009
:Backprojection 연산자(operator)4.
Figure pat00009
: Backprojection operator (operator)

5.

Figure pat00010
:Riesz potential 연산자5.
Figure pat00010
: Riesz potential operator

상기 CT 영상으로부터 인공물의 영상왜곡을 제거하거 위하여, 고밀도 물질을 포함하는 영역과 인공물 보정 파라미터를 이용한다. 인공물 보정 파라미터는 물질의 감쇄계수 및 엑스레이 스펙트럼에 의해 설정된다. To remove image distortion of the artifact from the CT image, an area containing high density material and an artifact correction parameter are used. The artifact correction parameters are set by the attenuation coefficient of the material and the X-ray spectrum.

인공물의 영상왜곡을 제거하기 위해 빔 경화 보정 함수를 이용하고, 이는 다음 수학식 3과 같이 나타낼 수 있다.The beam hardening correction function is used to remove the image distortion of the artifact, and this can be expressed by the following equation (3).

Figure pat00011
Figure pat00011

여기서, 사용되는 기호는 다음과 같다.Here, the symbols used are as follows.

1.

Figure pat00012
:고밀도 물질을 포함하는 영역, 여기서, Di는 D의 connected subdomain이다.One.
Figure pat00012
: A region containing a dense material, where Di is a connected subdomain of D;

2.

Figure pat00013
: D에서의 characteristic 함수. 즉,
Figure pat00014
on D,
Figure pat00015
otherwise.2.
Figure pat00013
The characteristic function at D. In other words,
Figure pat00014
on D,
Figure pat00015
otherwise.

3.

Figure pat00016
: 인공물 보정 파라미터, N 개의 서로 다른 물질의 감쇄계수와 X-ray 에너지 스펙트럼에 의존하는 파라미터이다.3.
Figure pat00016
: Parameters that depend on artifact correction parameters, the attenuation coefficients of the N different materials, and the X-ray energy spectrum.

인공물 보정 파라미터에 적용되는 물질의 감쇄계수는 CT 영상을 생성하고자 하는 대상물체에 따라 달라질 수 있다. 수학식 4를 풀어 설정될 수도 있고, 사용자의 설정에 따를 수도 있다. The attenuation coefficient of the material applied to the artifact correction parameter may vary depending on the object on which the CT image is to be generated. May be set by solving Equation (4), or may be set according to the setting of the user.

여기에서 라돈 변환 R과 backprojection operator R*는 상용화된 CT에서 (예를 들어, 2D fan beam CT, 3D cone beam CT) 각각 projection operator와 backprojection operator로 대체될 수 있다. Here, the radon transformer R and the backprojection operator R * can be replaced by a projection operator and a backprojection operator, respectively, in commercialized CT (for example, 2D fan beam CT and 3D cone beam CT).

빔 경화 보정 함수 φD,λ는 빔 경화 현상에 의한 검은색 흰색 줄무늬의 인공물을 포함하며, 고밀도 물질을 포함하는 영역 D와 파라미터 λ에 의해 결정된다. D와 λ는 unknown이며, D는 thresholding 기법을 이용하여 CT 영상 fCT로부터 영역분할(segmentation)할 수 있다. The beam hardening correction function φ D, λ includes artefacts of black white stripe due to beam hardening phenomenon, and is determined by the area D and the parameter λ containing the high density material. D and λ are unknown, and D can be segmented from the CT image f CT using the thresholding technique.

인공물 보정 파라미터(λ)는 보정 영상(fCTD,λ)에서 인공물을 최소화하도록 설정된다. 이는 다음 수학식 4와 같다.The artifact correction parameter (lambda) is set to minimize artifacts in the corrected image (f CTD, lambda ). This is shown in Equation (4).

Figure pat00017
Figure pat00017

여기서,

Figure pat00018
는 D의 바깥영역이다. 인공물 보정 파라미터에 따른 인공물은 도 4와 같다. 도 4에서와 같이, 인공물이 최소가 되는 인공물 보정 파라미터(λ=71.65)를 설정할 수 있다. here,
Figure pat00018
Is the outer region of D. The artifacts according to the artifact correction parameters are shown in Fig. As shown in FIG. 4, an artifact correction parameter (? = 71.65) in which artifacts are minimized can be set.

상기 수학식 4에서의 W(x)는 줄무늬 형태의 인공물에 더 많은 penalty를 부과하는 weight함수로써 다음과 같이 주어진다.W (x) in Equation (4) is a weight function that imposes a greater penalty on the striped artifacts as follows.

Figure pat00019
Figure pat00019

Figure pat00020
Figure pat00021
는 각각 Laplacian과 gradient 연산자를 나타낸다.
Figure pat00020
Wow
Figure pat00021
Represent the Laplacian and gradient operators, respectively.

도 5는 본 발명의 실시예에 따른 엑스레이 컴퓨터 단층촬영 영상의 인공물 보정방법의 흐름도이다.5 is a flowchart illustrating an artifact correction method of an X-ray computed tomography image according to an embodiment of the present invention.

도 5에 대한 상세한 설명은 도 1 내지 도 4의 엑스레이 컴퓨터 단층촬영 영상의 인공물 보정장치에 대한 상세한 설명에 대응하는바, 중복되는 설명은 이하 생략하도록 한다.5 corresponds to a detailed description of an artifact correction apparatus of the X-ray computed tomography image of FIGS. 1 to 4, and a duplicate description will be omitted below.

510 단계는 엑스레이 컴퓨터 단층촬영 영상으로부터 고밀도 물질을 포함하는 영역을 분할하는 단계이다.Step 510 is a step of dividing the region containing the high-density material from the X-ray computed tomography image.

보다 구체적으로, 고밀도 물질을 포함하는 영역이 어디인지를 알기 위하여, 엑스레이 컴퓨터 단층촬영 영상으로부터 경계화(thresholing) 방법을 이용하여 상기 고밀도 물질을 포함하는 영역을 분할한다.More specifically, in order to know the region containing the high-density material, the region containing the high-density material is divided from the X-ray computed tomography image using a threshing method.

520 단계는 상기 분할된 고밀도 물질을 포함하는 영역과 인공물 보정 파라미터를 이용하여 영상의 인공물을 보정하는 단계이다.Step 520 is a step of correcting the artifacts of the image using the area including the divided high density material and the artifact correction parameters.

보다 구체적으로, 상기 인공물 보정 파라미터는 물질의 감쇄계수 및 엑스레이 스펙트럼에 의해 설정된다. 상기 인공물 보정 파라미터는 줄무늬 모양의 인공물에 가중치를 부여하여 보정영상에서 인공물을 최소화하도록 설정될 수 있다.
More specifically, the artifact correction parameter is set by the attenuation coefficient of the material and the X-ray spectrum. The artifact correction parameters may be set to minimize artifacts in the corrected image by weighting stripe artifacts.

본 발명의 실시예들은 다양한 컴퓨터 수단을 통하여 수행될 수 있는 프로그램 명령 형태로 구현되어 컴퓨터 판독 가능 매체에 기록될 수 있다. 상기 컴퓨터 판독 가능 매체는 프로그램 명령, 데이터 파일, 데이터 구조 등을 단독으로 또는 조합하여 포함할 수 있다. 상기 매체에 기록되는 프로그램 명령은 본 발명을 위하여 특별히 설계되고 구성된 것들이거나 컴퓨터 소프트웨어 당업자에게 공지되어 사용 가능한 것일 수도 있다. 컴퓨터 판독 가능 기록 매체의 예에는 하드 디스크, 플로피 디스크 및 자기 테이프와 같은 자기 매체(magnetic media), CD-ROM, DVD와 같은 광기록 매체(optical media), 플롭티컬 디스크(floptical disk)와 같은 자기-광 매체(magneto-optical media), 및 롬(ROM), 램(RAM), 플래시 메모리 등과 같은 프로그램 명령을 저장하고 수행하도록 특별히 구성된 하드웨어 장치가 포함된다. 프로그램 명령의 예에는 컴파일러에 의해 만들어지는 것과 같은 기계어 코드뿐만 아니라 인터프리터 등을 사용해서 컴퓨터에 의해서 실행될 수 있는 고급 언어 코드를 포함한다. 상기된 하드웨어 장치는 본 발명의 동작을 수행하기 위해 하나 이상의 소프트웨어 모듈로서 작동하도록 구성될 수 있으며, 그 역도 마찬가지이다.Embodiments of the present invention may be implemented in the form of program instructions that can be executed on various computer means and recorded on a computer readable medium. The computer-readable medium may include program instructions, data files, data structures, and the like, alone or in combination. The program instructions recorded on the medium may be those specially designed and constructed for the present invention or may be available to those skilled in the art of computer software. Examples of computer-readable media include magnetic media such as hard disks, floppy disks and magnetic tape; optical media such as CD-ROMs and DVDs; magnetic media such as floppy disks; Magneto-optical media, and hardware devices specifically configured to store and execute program instructions such as ROM, RAM, flash memory, and the like. Examples of program instructions include machine language code such as those produced by a compiler, as well as high-level language code that can be executed by a computer using an interpreter or the like. The hardware devices described above may be configured to operate as one or more software modules to perform the operations of the present invention, and vice versa.

이상과 같이 본 발명에서는 구체적인 구성 요소 등과 같은 특정 사항들과 한정된 실시예 및 도면에 의해 설명되었으나 이는 본 발명의 보다 전반적인 이해를 돕기 위해서 제공된 것일 뿐, 본 발명은 상기의 실시예에 한정되는 것은 아니며, 본 발명이 속하는 분야에서 통상적인 지식을 가진 자라면 이러한 기재로부터 다양한 수정 및 변형이 가능하다. As described above, the present invention has been described with reference to particular embodiments, such as specific constituent elements, and limited embodiments and drawings. However, it should be understood that the present invention is not limited to the above- And various modifications and changes may be made thereto by those skilled in the art to which the present invention pertains.

따라서, 본 발명의 사상은 설명된 실시예에 국한되어 정해져서는 아니되며, 후술하는 특허청구범위뿐 아니라 이 특허청구범위와 균등하거나 등가적 변형이 있는 모든 것들은 본 발명 사상의 범주에 속한다고 할 것이다.Accordingly, the spirit of the present invention should not be construed as being limited to the embodiments described, and all of the equivalents or equivalents of the claims, as well as the following claims, belong to the scope of the present invention .

100: 영상의 인공물 보정장치
110: 처리부
120: 통신부
100: Image artifact correction device
110:
120:

Claims (8)

엑스레이 컴퓨터 단층촬영 영상의 인공물 보정방법에 있어서,
엑스레이 컴퓨터 단층촬영 영상으로부터 고밀도 물질을 포함하는 영역을 분할하는 단계; 및
상기 분할된 고밀도 물질을 포함하는 영역과 인공물 보정 파라미터를 이용하여 영상의 인공물을 보정하는 단계를 포함하고,
상기 인공물 보정 파라미터는 물질의 감쇄계수 및 엑스레이 스펙트럼에 의해 설정되는 것을 특징으로 하는 방법.
An artifact correction method for an X-ray computed tomography image,
Dividing a region containing high density material from an x-ray computed tomography image; And
And correcting an artifact of an image using an artifact correction parameter and a region including the divided high density material,
Wherein the artifact correction parameter is set by the attenuation coefficient of the material and the X-ray spectrum.
제 1 항에 있어서,
상기 인공물 보정 파라미터는,
보정 영상에서 인공물을 최소화하도록 설정되는 것을 특징으로 하는 방법.
The method according to claim 1,
Wherein the artificial-
And to minimize artifacts in the corrected image.
제 1 항에 있어서,
상기 인공물 보정 파라미터는,
인공물에 가중치를 부여하는 것을 특징으로 하는 방법.
The method according to claim 1,
Wherein the artificial-
Characterized in that the artifact is weighted.
제 1 항 내지 제 3 항 중에 어느 한 항의 방법을 컴퓨터에서 실행시키기 위한 프로그램을 기록한 컴퓨터로 읽을 수 있는 기록매체.A computer-readable recording medium storing a program for causing a computer to execute the method according to any one of claims 1 to 3. 엑스레이 컴퓨터 단층촬영 영상의 인공물 보정장치에 있어서,
엑스레이 컴퓨터 단층촬영 영상으로부터 고밀도 물질을 포함하는 영역을 분할하고, 상기 분할된 고밀도 물질을 포함하는 영역과 인공물 보정 파라미터를 이용하여 영상의 인공물을 보정하는 처리부를 포함하고,
상기 인공물 보정 파라미터는 물질의 감쇄계수 및 엑스레이 스펙트럼에 의해 설정되는 것을 특징으로 하는 장치.
An artifact correction apparatus for an X-ray computed tomography image,
A processor for dividing an area containing a high density material from an X-ray computed tomography image and correcting an artifact of the image using an area including the divided high density material and an artifact correction parameter,
Wherein the artifact correction parameter is set by the attenuation coefficient and x-ray spectrum of the material.
제 5 항에 있어서,
상기 인공물 보정 파라미터는,
보정 영상에서 인공물을 최소화하도록 설정되는 것을 특징으로 하는 장치.
6. The method of claim 5,
Wherein the artificial-
And is set to minimize artifacts in the corrected image.
제 5 항에 있어서,
상기 인공물 보정 파라미터는,
인공물에 가중치를 부여하는 것을 특징으로 하는 장치.
6. The method of claim 5,
Wherein the artificial-
And weighting the artifacts.
제 5 항에 있어서,
엑스레이 컴퓨터 단층촬영 영상 정보를 수신하는 통신부를 더 포함하는 장치.


6. The method of claim 5,
And a communication unit for receiving the X-ray computed tomographic image information.


KR1020150123301A 2015-09-01 2015-09-01 The beam-hardening correction method in X-ray CT KR101717433B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020150123301A KR101717433B1 (en) 2015-09-01 2015-09-01 The beam-hardening correction method in X-ray CT

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020150123301A KR101717433B1 (en) 2015-09-01 2015-09-01 The beam-hardening correction method in X-ray CT

Publications (2)

Publication Number Publication Date
KR20170026984A true KR20170026984A (en) 2017-03-09
KR101717433B1 KR101717433B1 (en) 2017-03-17

Family

ID=58402741

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020150123301A KR101717433B1 (en) 2015-09-01 2015-09-01 The beam-hardening correction method in X-ray CT

Country Status (1)

Country Link
KR (1) KR101717433B1 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111476856A (en) * 2020-04-08 2020-07-31 中北大学 Multispectral CT imaging method
CN114129912A (en) * 2021-11-25 2022-03-04 北京理工大学 Method for verifying radiation treatment dose

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101982941B1 (en) 2017-12-18 2019-08-28 연세대학교 원주산학협력단 Method and Apparatus for removing artifact in CT image using fuzzy neural network

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001137228A (en) * 1999-09-07 2001-05-22 General Electric Co <Ge> Method and apparatus for two-pass cone beam image reconstruction
JP2006105975A (en) * 2004-09-30 2006-04-20 General Electric Co <Ge> Ct reconstruction method and system equipped with preliminary correction
KR20080010674A (en) 2006-07-27 2008-01-31 동명대학교산학협력단 Image distortion compensation by using a polynomial model in an x-ray digital tomosynthesis system
US7889834B2 (en) * 2008-04-10 2011-02-15 Siemens Aktiengesellschaft Method for preparing reconstructed CT image data records and CT system
KR20130028057A (en) * 2010-02-02 2013-03-18 플란메카 오이 Dental computed tomography apparatus
KR20130058015A (en) * 2010-04-29 2013-06-03 플란메드 오와이 Medical computed tomography imaging apparatus

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001137228A (en) * 1999-09-07 2001-05-22 General Electric Co <Ge> Method and apparatus for two-pass cone beam image reconstruction
JP2006105975A (en) * 2004-09-30 2006-04-20 General Electric Co <Ge> Ct reconstruction method and system equipped with preliminary correction
KR20080010674A (en) 2006-07-27 2008-01-31 동명대학교산학협력단 Image distortion compensation by using a polynomial model in an x-ray digital tomosynthesis system
US7889834B2 (en) * 2008-04-10 2011-02-15 Siemens Aktiengesellschaft Method for preparing reconstructed CT image data records and CT system
KR20130028057A (en) * 2010-02-02 2013-03-18 플란메카 오이 Dental computed tomography apparatus
KR20130058015A (en) * 2010-04-29 2013-06-03 플란메드 오와이 Medical computed tomography imaging apparatus

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111476856A (en) * 2020-04-08 2020-07-31 中北大学 Multispectral CT imaging method
CN111476856B (en) * 2020-04-08 2023-06-06 中北大学 Multispectral CT imaging method
CN114129912A (en) * 2021-11-25 2022-03-04 北京理工大学 Method for verifying radiation treatment dose

Also Published As

Publication number Publication date
KR101717433B1 (en) 2017-03-17

Similar Documents

Publication Publication Date Title
US8233586B1 (en) Iterative reduction of artifacts in computed tomography images using forward projection and an edge-preserving blur filter
Meyer et al. Normalized metal artifact reduction (NMAR) in computed tomography
JP6106365B2 (en) How to get breast morphological information
Wang et al. Metal artifact reduction in CT using fusion based prior image
US9934597B2 (en) Metal artifacts reduction in cone beam reconstruction
Wei et al. X-ray CT high-density artefact suppression in the presence of bones
CN106683143B (en) Image metal artifact correction method
KR20120116940A (en) Enhanced image data/dose reduction
US8768037B2 (en) Method for correcting image artifacts occuring as a result of an anti scatter grid
JP2007014759A (en) Method for reducing image-based artifact in pet/ct imaging
CN111815521B (en) Cone beam CT metal artifact correction algorithm based on priori image
Byl et al. Photon‐counting normalized metal artifact reduction (NMAR) in diagnostic CT
CN111402150B (en) CT image metal artifact removal method and device
Chang et al. Prior-guided metal artifact reduction for iterative X-ray computed tomography
WO2018126434A1 (en) Ct image shadow correction method and apparatus, and electronic device
KR101717433B1 (en) The beam-hardening correction method in X-ray CT
Yazdi et al. An opposite view data replacement approach for reducing artifacts due to metallic dental objects
CA2729607A1 (en) Incorporation of mathematical constraints in methods for dose reduction and image enhancement in tomography
Yang et al. Fast iterative beam hardening correction based on frequency splitting in computed tomography
EP3404618B1 (en) Poly-energetic reconstruction method for metal artifacts reduction
CN111080740A (en) Image correction method, device, equipment and medium
Pua et al. An image-based reduction of metal artifacts in computed tomography
US10074196B2 (en) Reconstructing a three-dimensional image dataset from two-dimensional projection images, X-ray device and computer program
Sohn et al. Analytical low-dose CBCT reconstruction using non-local total variation regularization for image guided radiation therapy
Anhaus et al. A nonlinear scaling‐based normalized metal artifact reduction to reduce low‐frequency artifacts in energy‐integrating and photon‐counting CT

Legal Events

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