KR20170026984A - The beam-hardening correction method in X-ray CT - Google Patents
The beam-hardening correction method in X-ray CT Download PDFInfo
- 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
Links
- 238000012937 correction Methods 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 title claims abstract description 30
- 239000000463 material Substances 0.000 claims abstract description 43
- 238000002083 X-ray spectrum Methods 0.000 claims abstract description 8
- 238000002591 computed tomography Methods 0.000 claims description 27
- 238000004891 communication Methods 0.000 claims description 4
- 230000006870 function Effects 0.000 description 8
- 230000009977 dual effect Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 229910052704 radon Inorganic materials 0.000 description 2
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 2
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 239000007943 implant Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000007769 metal material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/56—Details of data transmission or power supply, e.g. use of slip rings
- A61B6/563—Details of data transmission or power supply, e.g. use of slip rings involving image data transmission via a network
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Networks & Wireless Communication (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
Description
본 발명은 엑스레이 컴퓨터 단층촬영 영상의 인공물 보정방법에 관한 것으로서, 더욱 상세하게는 고밀도 물질을 영역에 대한 인공물 보정 파라미터를 이용하여 영상의 인공물을 보정하는 방법 및 그 장치에 관한 것이다.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.
본 발명이 해결하고자 하는 첫 번째 과제는 고밀도 물질을 영역에 대한 인공물 보정 파라미터를 이용하여 영상의 인공물을 보정하는 영상의 인공물 보정방법을 제공하는 것이다.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
통신부(120)는 엑스레이 컴퓨터 단층촬영 영상정보를 수신한다. 또한, 영상의 인공물이 보정된 영상을 사용자에게 제공할 수 있다.The
처리부(110)는 엑스레이 컴퓨터 단층촬영 영상으로부터 고밀도 물질을 포함하는 영역을 분할하고, 상기 분할된 고밀도 물질을 포함하는 영역과 인공물 보정 파라미터를 이용하여 영상의 인공물을 보정하는 처리부를 포함한다.The
보다 구체적으로, 영상의 인공물을 보정함에 있어서, 분할된 고밀도 영역과 인공물 보정 파라미터에 대응하는 인공물 보정함수인 빔 경화 보정 함수를 도출한다. 빔 경화 보정 함수를 엑스레이 컴퓨터 단층촬영 영상에 적용하여 인공물에 의한 영상 왜곡이 없는 영상을 생성할 수 있다.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,λ)을 제거함으로써 영상의 인공물이 보정된 영상(fCT-φD,λ)을 얻을 수 있다.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.
여기서, X-ray 데이터 는 다음과 같이 정의된다.Here, the X-ray data Is defined as follows.
상기 수학식 1,2에 사용된 기호는 다음과 같다.The symbols used in the above equations (1) and (2) are as follows.
1. :위치, 에너지 E에서의 감쇄계수One. :location , The attenuation coefficient at energy E
2. :X-ray 에너지 스펙트럼2. : X-ray energy spectrum
3. , :라돈변환(Radon transform)3. , : Radon transform
4. :Backprojection 연산자(operator)4. : Backprojection operator (operator)
5. :Riesz potential 연산자5. : 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).
여기서, 사용되는 기호는 다음과 같다.Here, the symbols used are as follows.
1. :고밀도 물질을 포함하는 영역, 여기서, Di는 D의 connected subdomain이다.One. : A region containing a dense material, where Di is a connected subdomain of D;
2.: D에서의 characteristic 함수. 즉, on D, otherwise.2. The characteristic function at D. In other words, on D, otherwise.
3. : 인공물 보정 파라미터, N 개의 서로 다른 물질의 감쇄계수와 X-ray 에너지 스펙트럼에 의존하는 파라미터이다.3. : 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.
인공물 보정 파라미터(λ)는 보정 영상(fCT-φD,λ)에서 인공물을 최소화하도록 설정된다. 이는 다음 수학식 4와 같다.The artifact correction parameter (lambda) is set to minimize artifacts in the corrected image (f CT -φ D, lambda ). This is shown in Equation (4).
여기서, 는 D의 바깥영역이다. 인공물 보정 파라미터에 따른 인공물은 도 4와 같다. 도 4에서와 같이, 인공물이 최소가 되는 인공물 보정 파라미터(λ=71.65)를 설정할 수 있다. here, 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.
와 는 각각 Laplacian과 gradient 연산자를 나타낸다.
Wow 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.
상기 인공물 보정 파라미터는,
보정 영상에서 인공물을 최소화하도록 설정되는 것을 특징으로 하는 방법.The method according to claim 1,
Wherein the artificial-
And to minimize artifacts in the corrected image.
상기 인공물 보정 파라미터는,
인공물에 가중치를 부여하는 것을 특징으로 하는 방법.The method according to claim 1,
Wherein the artificial-
Characterized in that the artifact is weighted.
엑스레이 컴퓨터 단층촬영 영상으로부터 고밀도 물질을 포함하는 영역을 분할하고, 상기 분할된 고밀도 물질을 포함하는 영역과 인공물 보정 파라미터를 이용하여 영상의 인공물을 보정하는 처리부를 포함하고,
상기 인공물 보정 파라미터는 물질의 감쇄계수 및 엑스레이 스펙트럼에 의해 설정되는 것을 특징으로 하는 장치.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.
상기 인공물 보정 파라미터는,
보정 영상에서 인공물을 최소화하도록 설정되는 것을 특징으로 하는 장치.6. The method of claim 5,
Wherein the artificial-
And is set to minimize artifacts in the corrected image.
상기 인공물 보정 파라미터는,
인공물에 가중치를 부여하는 것을 특징으로 하는 장치.6. The method of claim 5,
Wherein the artificial-
And weighting the artifacts.
엑스레이 컴퓨터 단층촬영 영상 정보를 수신하는 통신부를 더 포함하는 장치.
6. The method of claim 5,
And a communication unit for receiving the X-ray computed tomographic image information.
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)
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)
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)
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 |
-
2015
- 2015-09-01 KR KR1020150123301A patent/KR101717433B1/en active IP Right Grant
Patent Citations (6)
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)
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 | |
CN102667852B (en) | Strengthen view data/dosage to reduce | |
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 | |
CN111815521B (en) | Cone beam CT metal artifact correction algorithm based on priori image | |
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 | |
Byl et al. | Photon‐counting normalized metal artifact reduction (NMAR) in diagnostic CT | |
WO2018126434A1 (en) | Ct image shadow correction method and apparatus, and electronic device | |
Chang et al. | Prior-guided metal artifact reduction for iterative X-ray computed tomography | |
CN111402150B (en) | CT image metal artifact removal method and device | |
CA2729607A1 (en) | Incorporation of mathematical constraints in methods for dose reduction and image enhancement in tomography | |
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 | |
Yang et al. | Fast iterative beam hardening correction based on frequency splitting in computed tomography | |
Anhaus et al. | A nonlinear scaling‐based normalized metal artifact reduction to reduce low‐frequency artifacts in energy‐integrating and photon‐counting CT | |
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 | |
KR102248191B1 (en) | Method and apparatus for correcting cone-beam artifacts in a cone-beam computed tomography image, cone-beam computed tomography apparatus including the same |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
E701 | Decision to grant or registration of patent right | ||
GRNT | Written decision to grant |