KR20130124895A - 단층 화상의 윤곽 추출 방법, 윤곽 화상 형성 프로그램을 기록한 컴퓨터 판독 가능한 기록 매체, 및 윤곽 화상 형성 장치 - Google Patents

단층 화상의 윤곽 추출 방법, 윤곽 화상 형성 프로그램을 기록한 컴퓨터 판독 가능한 기록 매체, 및 윤곽 화상 형성 장치 Download PDF

Info

Publication number
KR20130124895A
KR20130124895A KR1020130047979A KR20130047979A KR20130124895A KR 20130124895 A KR20130124895 A KR 20130124895A KR 1020130047979 A KR1020130047979 A KR 1020130047979A KR 20130047979 A KR20130047979 A KR 20130047979A KR 20130124895 A KR20130124895 A KR 20130124895A
Authority
KR
South Korea
Prior art keywords
projection data
sinogram
image
binarization
image forming
Prior art date
Application number
KR1020130047979A
Other languages
English (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
Priority claimed from EP12167018.6A external-priority patent/EP2662022A1/en
Application filed by 도쿠리츠교세이호진 고쿠리츠쥰칸키뵤 겐큐센터, 니혼 메디피직스 가부시키가이샤 filed Critical 도쿠리츠교세이호진 고쿠리츠쥰칸키뵤 겐큐센터
Publication of KR20130124895A publication Critical patent/KR20130124895A/ko

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00—Geometric image transformations in the plane of the image
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • G06T7/12—Edge-based segmentation
    • 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/037—Emission tomography
    • 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/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/507—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for determination of haemodynamic parameters, e.g. perfusion 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/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10081—Computed x-ray tomography [CT]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Pathology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Biophysics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Nuclear Medicine (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

본 발명은 피험자의 핵의학 화상의 프로젝션 데이터를 이용하여 윤곽 화상을 양호한 정밀도로 작성하는 방법을 제공한다.
산란 보정 및 흡수 보정을 실시하지 않은 피험자의 핵의학 프로젝션 데이터를 이용하여 역치법에 의해 화소값을 2치화 처리하여, 윤곽 정보를 추출한다. 필요에 따라 보간 처리를 실시한 후, 화상 재구성 처리를 실시하여 재구성 화상을 얻는다. 얻어진 재구성 화상에 대해 제2 2치화 처리를 실시하여, 윤곽 화상을 얻는다.

Description

단층 화상의 윤곽 추출 방법, 윤곽 화상 형성 프로그램을 기록한 컴퓨터 판독 가능한 기록 매체, 및 윤곽 화상 형성 장치{METHOD OF EXTRACTING CONTOUR OF TOMOGRAM, COMPUTER-READABLE RECORDING MEDIUM HAVING CONTOUR IMAGE FORMING PROGRAM, AND CONTOUR IMAGE FORMING APPARATUS}
본 발명은, 양전자 방출형 단층 촬상(이하,「PET」로 칭한다) 화상 및 단광자 방출형 단층 촬상(이하,「SPECT」로 칭한다) 화상으로 대표되는 핵의학 화상의 화상 처리 방법, 당해 방법을 실시하기 위한 프로그램, 및 화상 처리 장치에 관한 것이다. 보다 자세하게는, PET 화상 및 SPECT 화상으로 대표되는 핵의학 화상을 이용한, 단층 화상의 윤곽 추출 방법, 당해 방법을 실시하기 위한 프로그램, 및 화상 처리 장치에 관한 것이다.
PET 화상 및 SPECT 화상으로 대표되는 핵의학 화상은, 심장 질환이나 암을 비롯한 여러 질환의 진단에 유효하다. 이들 화상은, 특정 방사성 동위 원소로 라벨 된 약제(이하,「방사성 의약품」으로 칭한다)를 투여하고, 당해 약제로부터 방출된 γ선을 전용 카메라에 의해 검출하여, 재구축함으로써 얻어진다. 핵의학 화상은, 질환에 대한 특이도나 감도가 높은 우수한 성질을 갖고 있을 뿐만 아니라, 병변부의 기능에 관한 정보를 얻을 수 있다는, 다른 진단 화상에는 없는 특징을 갖고 있다.
상술한 바와 같이, 핵의학 화상은 피험자에게 투여된 방사성 의약품으로부터 방출된 γ선을 검출함으로써 얻어지지만, γ선은 생체 내에 있어서 감쇠 및 산란의 영향을 받기 때문에, 검출된 γ선의 강도는 체표면으로부터의 깊이에 따라 증감된 것이 된다. 따라서, 검출된 γ선 강도를 이용하여 핵의학 화상을 재구축할 때에는, 생체에서의 γ선의 감쇠 및 산란을 고려한 보정을 실시할 필요가 있다. 이것은 핵의학 화상을 생체의 정량 평가에 이용할 때 특히 중요하다.
이와 같은 보정에는, 많은 경우, 대상 부위에 있어서의 흡수 계수(μ값)의 분포를 나타내는 흡수맵(μ맵)이 이용된다. 즉, μ맵으로부터 산란선 분포를 추정한 후에, 프로젝션 데이터에 당해 산란선을 고려한 보정을 실시하고, 그 후에, μ맵을 이용한 화상 재구성을 실시하여 감쇠 보정이 실시된다(비특허문헌 1).
상술한 보정에 이용되는 μ맵의 작성시에는, 피험자의 윤곽을 추출하는 것이 필요로 된다. 이 윤곽 추출은, 많은 경우, 피험자의 핵의학 화상 데이터를 이용하여 역치법에 의해 실시된다. 이 경우, 윤곽 추출에 이용하는 핵의학 화상 데이터는, 양호한 화상을 이용하여 윤곽 추출을 실시하고자 하는 요청에 의해, 감쇠 및 산란선 보정을 실시한 화상을 이용하는 것이 통상적이다. 그리고, 이 윤곽 추출용 화상의 화상 재구성시에는, 피험자의 형상을 타원 근사하여 얻어진 μ맵을 이용한 산란 및 흡수 보정 방법이나, μ맵을 이용하지 않는 산란 및 흡수 보정 방법이 실시되고 있었다.
피험자의 재구성 화상으로부터 역치법에 의해 윤곽을 추출하는 방법은, 가장 간편한 방법이기 때문에, 현재 널리 이용되고 있다. 그러나, 감쇠 및 산란 보정을 실시한 화상을 이용하는, 역치를 이용한 종래의 윤곽 추출 방법에서는, 윤곽을 과대 또는 과소 평가해 버리는 경우가 있었다. 그리고, 추출된 체윤곽의 형상이 부정확해지면, 산란 보정이나 감쇠 보정이 부정확해지기 때문에 바람직하지 않다.
본 발명은 이와 같은 사정을 감안하여 이루어진 것으로, 피험자의 핵의학 화상을 이용하여 정확한 윤곽 추출을 실시하는 방법, 당해 방법을 실시하기 위한 프로그램, 및 화상 처리 장치를 제공하는 것을 목적으로 하였다.
발명자는 검토한 결과, 산란선 보정을 실시하지 않은 투영 데이터 세트를 이용하여 역치법에 의해 2치화 처리를 실시하여, 당해 데이터로부터 재구성 화상을 얻음으로써, 핵의학 화상 데이터로부터 양호한 정밀도로 피험자의 윤곽 화상을 얻는 것이 가능한 것을 알아내어, 본 발명을 완성시켰다.
또한, 본 명세서에 있어서, 투영 데이터 세트란, 한 번의 촬상 처리에 의해 얻어진 각 투영 각도에서의 투영 화상 세트를 말한다. 예를 들면, 환자에 대해, 2°단계에서 360°의 투영을 실시한 경우, 투영 각도 2°, 4°, 6°, …, 358°, 360°에서 얻어진 180장의 투영 화상 세트가, 본 명세서에 있어서의 투영 데이터 세트가 된다.
또한, 사이노그램이란, 투영 데이터 세트를 구성하는 각 투영 화상을 투영 각도의 순서로 배열하여 표시한 데이터를 가리킨다.
또한, 수집 각도 범위는, 데이터 수집을 실시하는 투영 각도의 범위를 가리킨다. 예를 들면, 수집 각도 0°내지 2°단계에서 180°까지의 투영을 실시한 경우는, 수집 각도 범위 0∼180°가 된다. 또한, 예를 들면 수집 각도 범위의 하한값을 0°로 하고, 상한값을 180°이상으로 한 경우는, 수집 각도 범위를 0∼180°이상으로 표현하는 것으로 한다.
본 발명의 일측면에 따른 윤곽 화상 형성 방법은, 피험자의 투영 데이터로부터 산란 및 흡수 보정을 실시하지 않고 윤곽 화상을 생성하는 방법으로서, 컴퓨터에 의해, 수집 각도 범위 0∼180°이상의 핵의학 화상의 투영 데이터 세트를 취득하는 데이터 취득 단계와, 취득한 상기 투영 데이터 세트로부터 사이노그램을 작성하는 사이노그램 작성 단계와, 작성한 사이노그램에 대해 제1 2치화 처리를 실시하는 제1 2치화 단계와, 상기 제1 2치화 처리 후의 상기 사이노그램을 역투영하여 화상을 재구성하는 화상 재구성 단계와, 재구성된 화상에 대해 제2 2치화 처리를 실시하는 제2 2치화 단계를 순차 실행한다.
본 발명의 다른 일측면에 따른 윤곽 화상 형성 프로그램은, 컴퓨터에 판독 입력되어 당해 컴퓨터를 피험자의 투영 데이터로부터 산란 및 흡수 보정을 실시하지 않고 윤곽 화상을 생성하는 윤곽 화상 생성 장치로서 동작시키기 위한 프로그램으로서, 컴퓨터에 수집 각도 범위 0∼180°이상의 핵의학 화상의 투영 데이터 세트를 취득하는 데이터 취득 단계와, 취득한 상기 투영 데이터 세트로부터 사이노그램을 작성하는 사이노그램 작성 단계와, 상기 사이노그램에 대해 제1 2치화 처리를 실시하는 제1 2치화 단계와, 제1 2치화 후의 상기 사이노그램을 역투영하여 화상을 재구성하는 화상 재구성 단계와, 재구성된 화상에 대해 제2 2치화 처리를 실시하는 제2 2치화 단계를 순차 실행시킨다.
본 발명의 또 다른 일측면에 따른 윤곽 화상 형성 장치는, 피험자의 투영 데이터로부터 산란 및 흡수 보정을 실시하지 않고 윤곽 화상을 생성하기 위한 윤곽 화상 형성 장치로서, 핵의학 화상의 투영 데이터 세트를 취득하는 기능을 갖는 데이터 취득부와, 투영 데이터 세트로부터 사이노그램을 작성하는 사이노그램 작성부와, 사이노그램에 대해 데이터의 2치화 처리를 실시하는 기능을 갖는 제1 2치화 처리부와, 사이노그램을 역투영하여 화상을 재구성하는 기능을 갖는 화상 재구성 처리부와, 재구성된 화상에 대해 화상의 2치화 처리를 실시하는 기능을 갖는 제2 2치화 처리부를 구비한다.
본 발명은 산란 및 흡수 보정을 실시하지 않은 투영 데이터를 이용하여 간편한 역치법에 의해 양호한 정밀도로 윤곽을 추출할 수 있는 것을 알아낸 점에 특징을 갖는 것이다. 상기에도 서술한 바와 같이, 양호한 화상을 이용하여 화상의 윤곽을 구하고자 한 요청으로부터, 핵의학 화상에서의 윤곽 추출은 종래는 산란 보정이나 흡수 보정을 실시한 단층 화상 상에서 실시되는 것이 기술 상식이었다. 그러나, 발명자가 검토한 결과, 산란선 보정을 실시함으로써 화상의 윤곽이 변화되어 버리는 경우가 있어, 이것에 기인하여, 결정된 윤곽이 피험자의 화상을 양호한 정밀도로 재현할 수 없는 경우가 있는 것을 알아내었다. 또한, 발명자가 검토를 실시한 결과, 산란선 보정 및 흡수 보정을 실시하지 않은 투영 데이터에 대해 역치법에 의해 윤곽 추출을 실시한 경우여도, 양호한 정밀도로 윤곽 추출을 실시하는 것이 가능한 것을 알아내었다. 그리고, 발명자는 이 지견에 착안하여, 본 발명을 완성하기에 이른 것이다.
본 발명에 따른 윤곽 화상 형성 방법, 윤곽 화상 형성 프로그램, 및, 윤곽 화상 형성 장치에 있어서, 제1 2치화 처리는, 일정한 역치 이하의 값을 0으로 하고, 당해 역치보다 큰 값을 어느 일정값(예를 들면, 1)으로 하는, 일반적으로 널리 실시되고 있는 2치화 처리와 동일한 처리에 의해 행해진다. 예를 들면, 역치를 최대값의 25％로 한 경우, 최대값의 25％이하의 화소값을 갖는 화소에 대해서는 화소값을 0으로 하고, 최대값의 25％보다 큰 화소값을 갖는 화소에 대해서는 화소를 1로 한다.
또한, 제1 2치화 처리에 있어서 사용되는 역치는, 2치화 처리 후의 투영 데이터 세트를 이용하여 재구성된 화상과, 2치화 처리를 실시하지 않고 작성한 재구성 화상의 비교에 기초하여 결정된 값을 사용할 수 있다. 당해 역치의 결정 방법에 대해서는 후술한다.
또한, 화상 재구성 단계 및 화상 재구성부에 있어서 실시되는 화상 재구성은, 핵의학 분야에서 널리 일반적으로 실시되고 있는 화상 재구성 방법을 이용할 수 있다. 본 발명에 있어서는, 2치화 처리 후의 투영 데이터 세트를 이용하여 화상 재구성을 실시하므로, 얻어진 화상은 중심 방향의 화소값이 커진 것 같은 화상이 된다. 여기서 본 발명에서는, 작성된 재구성 화상에 대해 제2 2치화 처리를 실시하여, 윤곽 화상 내부의 화소값을 균일화한다. 이 제2 2치화 처리에 사용하는 역치는, 윤곽 내부의 데이터를 컷하지 않고, 윤곽 내부의 화소값을 균일화할 수 있는 값이면, 특별히 한정할 필요는 없다. 통상적으로는, 최대값의 25∼50％ 정도의 값을 사용하면 된다.
본 발명에 따른 윤곽 화상 형성 방법은, 컴퓨터에 의해 제1 2치화 처리 후의 사이노그램에 대해 데이터 보간을 실시하는 보간 단계를 추가로 포함하고, 화상 재구성 단계는, 당해 보간 단계 완료 후의 사이노그램을 이용하여 역투영을 실시하는 구성으로 해도 된다. 동일하게, 본 발명에 따른 윤곽 화상 형성 프로그램은, 컴퓨터에 제1 2치화 단계 후의 사이노그램에 대해 데이터 보간을 실시하는 보간 단계를 추가로 실행시키고, 화상 재구성 단계는, 보간 단계 완료 후의 사이노그램을 이용하여 역투영을 실시하는 구성으로 해도 된다. 또한, 본 발명에 따른 윤곽 화상 형성 장치는, 사이노그램의 데이터 보간을 실시하는 기능을 갖는 보간 처리부를 추가로 포함한 구성으로 해도 된다.
이와 같은 구성으로 함으로써, 추출된 윤곽이 매끄러운 형상이 되어, 보다 양호한 윤곽 화상을 얻는 것이 가능해진다.
또한, 보간 단계에 있어서의 보간 처리는, 화상 처리 분야에 있어서 통상적으로 널리 이용되고 있는 데이터 보간 방법을 이용하여 실시할 수 있다. 예를 들면, 윤곽을 푸리에 급수 근사와 같이 주기 함수로 근사하는 방법이나, 2차 또는 3차의 Spline 보간 등의 방법을 이용하는 것이 가능하다.
또한, 윤곽 화상 형성 장치에서의 보간 처리부는, 상술한 바와 같은 보간 처리를 실행할 수 있는 기능을 구비하고 있는 것이면 된다.
본 발명에 따른 윤곽 화상 형성 방법 및 윤곽 화상 형성 프로그램에 있어서, 투영 데이터 세트의 수집 각도 범위의 상한값은, 180°이상이면 특별히 제한할 필요는 없지만, 360°이상으로 하면, 윤곽 추출의 정밀도가 보다 향상되기 때문에 바람직하다. 또한, 데이터의 낭비를 줄인다는 의미에서는, 수집 각도 범위는 0∼360°로 하는 것이 가장 바람직하다.
본 발명에 따른 윤곽 화상 형성 방법에 있어서, 데이터 취득 단계에서 취득하는 투영 데이터 세트의 수집 각도 범위를 0∼360°이상(바람직하게는 0∼360°)으로 하고, 수집 각도 범위 0∼180°의 범위에 있어서의 투영 데이터 세트를 구성하는 각각의 투영 데이터의 각 화소에 대해, 당해 투영 데이터와 180°대향하는 투영 데이터에 있어서의 대응하는 화소 사이에서 화소값의 합, 기하 평균, 또는 산술 평균을 산출하고, 얻어진 값을 당해 투영 데이터에 있어서의 화소값과 치환하는 평균화 단계를 추가로 실행하고, 사이노그램 작성 단계는 평균화 단계 완료 후의 수집 각도 범위 0∼180°의 범위에 있어서의 투영 데이터 세트를 이용하여 실시한다는 구성으로 해도 된다. 이 경우, 평균화 단계에서 실시하는 연산은, 투영 데이터 세트에 포함되는 모든 화소에 대해 전부 동일한 것을 사용한다. 즉, 예를 들면 어느 화소에 대해 기하 평균에 의해 화소값의 평균화 처리를 실행한 경우는, 나머지의 모든 화소에 대해서도 동일하게 기하 평균을 실시한다.
동일하게, 본 발명에 따른 윤곽 화상 형성 프로그램은, 데이터 취득 단계에 있어서 취득하는 투영 데이터 세트의 수집 각도 범위를 0∼360°이상(바람직하게는 0∼360°)으로 하고, 수집 각도 범위 0∼180°의 범위에 있어서의 투영 데이터 세트를 구성하는 각각의 투영 데이터의 각 화소에 대해, 당해 투영 데이터와 180°대향하는 투영 데이터에 있어서의 대응하는 화소 사이에서 화소값의 합, 기하 평균, 또는 산술 평균을 산출하고, 얻어진 값을 당해 투영 데이터에 있어서의 화소값과 치환하는 평균화 단계를 컴퓨터에 추가로 실행시키고, 사이노그램 작성 단계는 평균화 단계 완료 후의 수집 각도 범위 0∼180°의 범위에 있어서의 투영 데이터 세트를 이용하여 실시한다는 구성으로 해도 된다. 이 경우에 있어서도, 평균화 단계에서 실시하는 연산은, 투영 데이터 세트에 포함되는 모든 화소에 대해 전부 동일한 것을 사용하는 것은 당연하다.
동일하게, 본 발명에 따른 윤곽 화상 형성 장치는, 투영 데이터 세트를 구성하는 각각의 투영 데이터의 각 화소에 대해, 당해 투영 데이터와 180°대향하는 투영 데이터에 있어서의 대응하는 화소 사이에서 화소값의 합, 기하 평균, 또는 산술 평균을 산출하고, 얻어진 값을 당해 투영 데이터에 있어서의 화소값과 치환하는 기능을 갖는 화소값 평균화부를 추가로 구비한 구성으로 해도 된다. 이 경우도 상기와 같이, 평균화부에 있어서 사용하는 연산은, 투영 데이터 세트에 포함되는 모든 화소에 대해 전부 동일한 것을 사용한다.
이와 같은 구성으로 함으로써, 윤곽이 보다 선명해지기 때문에, 작성된 윤곽 화상의 정밀도를 보다 향상시키는 것이 가능해진다.
또한, 상기에 설명한 평균화 단계는, 투영 데이터 세트를 각도 방향으로 재배열하여 사이노그램으로부터, 각 투영 데이터에 대해 실시하는 구성으로 해도 된다.
즉, 본 발명에 따른 윤곽 화상 형성 방법에 있어서, 데이터 취득 단계에 있어서 취득하는 투영 데이터 세트의 수집 각도 범위를 0∼360°로 하고, 수집 각도 범위 0∼180°의 범위에 있어서의 사이노그램을 구성하는 각각의 투영 데이터의 각 화소에 대해, 당해 투영 데이터와 180°대향하는 투영 데이터에 있어서의 대응하는 화소 사이에서 화소값의 합, 기하 평균, 또는 산술 평균을 산출하고, 얻어진 값을 당해 투영 데이터에 있어서의 화소값과 치환하는 평균화 단계를 추가로 실행하여, 투영 각도 180°이상의 투영 데이터를 사이노그램으로부터 삭제하고, 제1 2치화 단계는, 평균화 단계 완료 후의 사이노그램을 이용하여 실시한다는 구성으로 해도 된다. 이 경우, 평균화 단계에서 실시하는 연산은, 투영 데이터 세트에 포함되는 모든 화소에 대해 전부 동일한 것을 사용한다. 즉, 예를 들면 어느 화소에 대해 기하 평균에 의해 화소값의 평균화 처리를 실행한 경우는, 나머지의 모든 화소에 대해서도, 동일하게 기하 평균을 실시한다.
동일하게, 본 발명에 따른 윤곽 화상 형성 프로그램은, 데이터 취득 단계에 있어서 취득하는 투영 데이터 세트의 수집 각도 범위를 0∼360°로 하고, 수집 각도 범위 0∼180°의 범위에 있어서의 사이노그램을 구성하는 각각의 투영 데이터의 각 화소에 대해, 당해 투영 데이터와 180°대향하는 투영 데이터에 있어서의 대응하는 화소 사이에서 화소값의 합, 기하 평균, 또는 산술 평균을 산출하고, 얻어진 값을 당해 투영 데이터에 있어서의 화소값과 치환하는 평균화 단계를 컴퓨터에 추가로 실행시켜, 투영 각도 180°이상의 투영 데이터를 사이노그램으로부터 삭제하고, 제1 2치화 단계는, 평균화 단계 완료 후의 사이노그램을 이용하여 실시하는 구성으로 해도 된다. 이 경우에 있어서도, 평균화 단계에서 실시하는 연산은 투영 데이터 세트에 포함되는 모든 화소에 대해 전부 동일한 것을 사용하는 것은 당연하다.
동일하게, 본 발명에 따른 윤곽 화상 형성 장치는, 사이노그램을 구성하는 각각의 투영 데이터의 각 화소에 대해, 당해 투영 데이터와 180°대향하는 투영 데이터에 있어서의 대응하는 화소 사이에서 화소값의 합, 기하 평균, 또는 산술 평균을 산출하고, 얻어진 값을 당해 투영 데이터에 있어서의 화소값과 치환하는 기능을 갖는 화소값 평균화부를 추가로 구비한 구성으로 해도 된다. 이 경우도 상기와 같이, 평균화부에 있어서 사용하는 연산은, 투영 데이터 세트에 포함되는 모든 화소에 대해 전부 동일한 것을 사용한다.
본 발명의 또 다른 일측면에 따른 윤곽 화상 형성 프로그램을 기록한 컴퓨터 판독 가능한 기록 매체는, 본 발명에 따른 윤곽 화상 형성 프로그램을 컴퓨터로 판독 가능한 기록 매체에 저장한 것이다.
본 발명에 의해, 피험자의 투영 데이터로부터, 양호한 윤곽 화상을 작성하는 것이 가능해졌다.
도 1은 본 발명에 따른 윤곽 화상 형성 장치의 바람직한 양태에 있어서의 처리의 개요를 나타내는 플로우 차트이다.
도 2는 본 발명에 따른 윤곽 화상 형성 장치의 바람직한 양태에 있어서의 기능 블록도이다.
도 3은 본 발명에 따른 윤곽 화상 형성 장치의 바람직한 양태에 있어서의 시스템 구성을 나타내는 도면이다.
도 4는 제1 2치화 처리에 이용하는 역치 결정 방법의 바람직한 일례에 있어서의 처리의 흐름을 나타내는 플로우 차트이다.
도 5는 두부 핵의학 화상 데이터에 있어서의 화소값 평균화 처리의 모식도이다.
도 6은 본 발명에 따른 윤곽 화상 형성 프로그램의 바람직한 양태에 있어서의 구성을 나타내는 도면이다.
도 7의 a는 본 발명에 따른 방법에 의해 얻어진 두부 윤곽 화상의 일례이고, b는 윤곽 화상 작성에 이용한 투영 데이터 세트에 의해 얻어진 재구성 화상(Filtered Back Projection법 사용)이다.
이하, 본 발명에 대해, 도면을 참조하여 설명한다. 또한, 이하에 나타내는 예는, 어디까지나 바람직한 형태에 대해 설명하는 것으로, 본 발명의 내용은 이들 기재에 의해 한정되는 것은 결코 아니다.
도 1은 본 발명에 따른 윤곽 화상 형성 장치의 바람직한 양태에 있어서의 처리의 개요를 나타내는 플로우 차트이고, 도 2는 본 발명에 따른 윤곽 화상 형성 장치의 바람직한 양태에 있어서의 기능 블록도이다. 본 발명에 따른 윤곽 화상 형성 장치(10)는, 윤곽 화상 형성 프로그램(300)을 판독 입력한 컴퓨터에 의해 구성할 수 있고, 윤곽 화상 형성 장치(10)를 동작시킴으로써, 본 발명에 따른 윤곽 화상 형성 방법을 실현할 수 있다.
바람직한 양태에 있어서, 윤곽 화상 형성 장치(10)는, 데이터 취득부(20), 사이노그램 작성부(30), 제1 2치화 처리부(40), 윤곽 추출부(50), 데이터 보간 처리부(60), 화상 재구성부(70), 제2 2치화 처리부(80), 및 출력부(90)에 의해 구성된다. 그리고, 바람직한 양태에 있어서, 윤곽 화상 형성 장치(10)는 SPECT 장치 등의 핵의학 화상 촬상 장치(100)와 전기 통신 회선을 통해서 접속된다.
도 3은, 본 발명에 따른 윤곽 화상 형성 장치(10)의 가장 바람직한 양태에 있어서의 시스템 구성이다. 바람직한 양태에 있어서, 윤곽 화상 형성 장치(10)는, CPU(230)와, 메모리(240)와, 모니터 등의 출력 기기(250)와, 통신 인터페이스(260)와, 키보드 등의 입력 장치(270)가, 버스(280)를 개재하여 접속되어 있다. 윤곽 화상 형성 장치(10)는, 그 밖에도 CD-ROM 드라이브나 USB 인터페이스 등을 구비하고 있어도 된다. 통신 인터페이스(260)는 핵의학 화상 촬상 장치(100)와 접속하기 위해서 사용된다. 또한, 메모리(240)에는 본 발명에 따른 윤곽 화상 형성 프로그램(300)이 기억되어 있다.
본 발명에 따른 윤곽 화상 형성 장치(10)는, 우선, 데이터 취득부(20)에 데이터 취득 단계를 실행시켜, 처리에 제공하는 핵의학 화상의 투영 데이터 세트를 컴퓨터 시스템에 입력한다(단계 S01). 여기서, 윤곽을 추출하기 위해서 필요한 정보를 얻는다는 목적에서, 투영 데이터 세트의 수집 각도 범위는 0∼180°이상으로 하고, 바람직하게는 0∼360°이상, 특히 바람직하게는 0∼360°로 한다.
투영 데이터 세트는, SPECT 장치나 PET 장치라는 통상적인 핵의학 화상 촬상 장치에 의해 촬상된 것을 이용할 수 있다. 투영 데이터 세트는 DICOM 형식 등의, 컴퓨터로 판독 가능한 형식으로 저장된 것을 핵의학 화상 촬상 장치(100)로부터, 네트워크를 통해서 직접 입력함으로써, 윤곽 화상 형성 장치(10)에 입력된다. 또한, 당해 투영 데이터 세트는, 하드 디스크, CD-ROM, DVD 등의, 컴퓨터로 판독 가능한 기록 매체에 저장된 상태로 취득한 것을, 컴퓨터 시스템에 설치된 판독 장치에 의해 판독 입력하는 방법으로 입력되는 것이어도 된다.
다음으로, 윤곽 화상 형성 장치(10)는, 사이노그램 작성부(30)에 사이노그램 작성 단계를 실행시켜, 투영 데이터 세트로부터 사이노그램을 작성한다(단계 S02).
사이노그램은 어느 1단면으로부터의 투영 데이터를 투영 각도의 순서로 배열하여 얻어지는 2차원 데이터를 말하고, 공지된 방법에 의해 작성할 수 있다. 예를 들면, 4°마다의 각도 간격으로 0∼360°(즉, 투영 각도 범위 0∼360°)에 걸쳐서 얻어진 투영 데이터(x, y)를 고려한 경우, 어느 y값(체축 방향의 어느 위치)에 있어서의 투영 데이터(x)를 0°에서 360°에 걸쳐서 2차원 좌표 위에서 아래로 순차 배치함으로써 임의의 단층면(y)에 있어서의 사이노그램을 작성할 수 있다.
사이노그램 작성 단계가 완료되면, 윤곽 화상 형성 장치(10)는 제1 2치화 처리부(40)에 제1 2치화 처리 단계를 실행시킨다(단계 S03). 이 경우에 있어서의 2치화 처리는, 예를 들면, 설정된 역치 이하의 화소값을 0으로, 역치보다 큰 화소값을 1로 치환한다는, 통상적으로 행해지는 2치화 처리 방법에 의해 실시할 수 있다.
제1 2치화 처리 단계에서 사용되는 역치의 값은, 재구성 화상과의 비교에 기초하여 별도 결정된 값을 사용할 수 있다. 도 4에, 제1 2치화 처리 단계에서 사용되는 역치 결정 방법의 바람직한 양태에 있어서의 플로우 차트를 나타낸다. 역치의 결정시에는, 우선은 일정한 초기값에 의해 2치화 처리를 실시하고(단계 S11), 역투영 처리에 의해 화상을 재구성한다(단계 S12). 여기서, 역치의 초기값은 사용자가 경험에 기초하여 임의로 결정한 값을 사용할 수 있다. 통상적으로는, 화소값의 최대값의 25％ 정도의 값을 사용하면 된다.
화상의 재구성 처리가 완료되면, 별도로 작성한 피험자의 재구성 화상과의 비교를 실시하여(단계 S13), 양자의 화상간에 있어서의 형상의 오차가 허용 범위 이내인지 여부를 판정한다(단계 S14). 여기서 사용하는 재구성 화상은, 처리에 사용하고 있는 투영 데이터 세트로부터, 산란, 흡수 보정을 수반하지 않는 공지된 화상 재구성법에 의해 생성된다. 이 경우의 화상 재구성법의 일례로는, Filtered back projection법 등을 들 수 있다. 또한, 오차의 판정은 양자의 중첩 화상을 사용하여, 사용자가 육안으로 판단할 수 있다. 물론, 잔차 제곱의 합 등의 오차의 지표가 되는 값을 계산하여 그 허용값과의 비교를 실시한다는 방법을 이용하여도 된다. 그리고, 오차가 허용 범위보다 큰 것으로 판단되었을 경우(단계 S14에서 no)는, 역치의 값을 변화시켜(단계 S15), 단계 S11∼S14의 처리를 반복한다. 오차가 허용 범위 내인 것으로 판단된 경우(단계 S14에서 yes)는, 그 때의 역치를 처리에 사용하는 역치로서 채용한다.
도 1로 돌아와, 윤곽 화상 형성 장치(10)에 있어서의 처리에 대해 설명한다.
제1 2치화 처리 단계가 완료되면, 윤곽 추출부(50)에 의해, 윤곽 추출 단계를 실시하여, 사이노그램 상에서의 윤곽 데이터의 추출을 실시한다(단계 S04). 바람직한 양태에 있어서, 윤곽 추출 단계는 2치화 처리 후의 투영 데이터에 대해 미분화 처리를 실시하고, 각 좌표축 상에 있어서의 좌우 양단으로부터의 최초의 피크 위치를 추출한다는 방법을 이용할 수 있다. 이 미분화 처리는, 투영 데이터 세트로부터 작성한 사이노그램에 소벨 필터나 프리윗 필터 등의 공지된 미분화 필터를 적용함으로써 실시할 수 있다. 이 처리에 의해, 사이노그램의 좌우에 2개의 윤곽 데이터가 추출된다.
윤곽 추출 단계가 완료되면, 데이터 보간 처리부(60)에 의해, 좌우 각각의 윤곽 데이터에 대해 보간 단계를 실행하여 보간 처리를 실시한다(단계 S05). 여기서의 보간 처리는, 윤곽을 푸리에 급수 근사와 같이 주기 함수로 근사하는 방법이나, 2차 또는 3차의 Spline 보간 등의 공지된 방법을 이용할 수 있다. 이 처리에 의해, 추출된 윤곽 데이터를 매끄러운 데이터로 변환할 수 있다.
보간 단계가 완료되면, 좌우 2개의 윤곽 데이터에 의해 사이에 낀 영역의 화소에 일정한 화소값(예를 들면, 1)을 배치한 다음, 화상 재구성부(70)에 의해 화상 재구성 단계를 실행하여, 사이노그램을 역투영함으로써, 재구성 화상을 생성한다(단계 S06). 여기서 형성된 재구성 화상은, 화상의 중심을 향해 화소값이 커져 있다. 여기서, 제2 2치화 처리부(80)를 동작시켜 제2 2치화 단계를 실행하여, 윤곽 화상 내의 화소값을 일정한 값으로 한다(단계 S07). 제2 2치화 단계에 있어서의 역치는, 윤곽 화상 데이터를 결손시키지 않는 정도의 크기이면 되고, 2치화 처리 전후의 화상을 비교함으로써 결정할 수 있다.
제2 2치화 처리 단계가 완료되면, 출력부(90)에 의해, 윤곽 화상을 디스플레이라는 출력 기기에 출력한다(단계 S08). 이들 일련의 처리를 실행함으로써, 윤곽 화상을 형성할 수 있다.
또한, 본 발명에 따른 윤곽 화상 형성 장치(10)는, 360°투영 데이터 세트를 구성하는 각각의 투영 데이터의 각 화소에 대해, 당해 투영 데이터와 180°대향하는 투영 데이터에 있어서의 대응하는 화소 사이에서 화소값의 합, 기하 평균, 또는 산술 평균을 산출하고, 얻어진 값을 당해 투영 데이터에 있어서의 화소값과 치환하는 기능을 갖는 화소값 평균화부(도시하지 않음)를 추가로 구비한 것이어도 된다. 이 경우, 데이터 취득 단계(단계 S01)에서 취득된 투영 데이터 세트에 대해서 화소값 평균화 단계(도시하지 않음)를 실행시키고, 사이노그램 작성 단계(단계 S02)에서는, 화소값 평균화 단계 실행 후의 수집 각도 범위 0∼180°의 범위에 있어서의 투영 데이터 세트를 이용하여, 사이노그램을 작성한다. 이러한 구성으로 함으로써, 평균화 처리 후의 투영 데이터 세트에서는 윤곽의 편차가 평활화되어 윤곽이 보다 선명해지기 때문에, 작성된 윤곽 화상의 정밀도를 보다 향상시키는 것이 가능해진다. 또한, 상술한 화소값 평균화 단계는, 사이노그램의 형상으로 한 투영 데이터 세트에 대해서, 개개의 투영 데이터를 이용하여 실시할 수도 있다.
화소값 평균화부는, 후술하는 화소값 평균화 모듈(도시하지 않음)을 장착한 윤곽 화상 형성 프로그램을 컴퓨터에 판독 입력함으로써 구성된다.
도 5에, 두부 핵의학 화상에 대해서 화소값 평균화 단계를 실시하는 경우의 예를 모식적으로 나타낸다. 도 5에 있어서, 투영상 A와 투영상 B는 서로 180°대향하고 있다. 이 경우에 있어서, 투영상 A와 투영상 B 사이에서, a1과 b1, a2와 b2, a3과 b3, …, an-1과 bn-1, an과 bn …과 같이 대응하는 화소끼리 화소값의 합, 기하 평균, 또는 산술 평균을 계산한다. 그리고, 각각의 화소에 대해 얻어진 계산값을, 투영상 A에 있어서의 a1, a2, a3, …, an-1, an…에 있어서의 화소값과 각각 치환한다. 이 처리를 투영 각도 0°∼180°에 걸쳐 순차 실시함으로써, 평균화 단계가 완료된다. 또한, 당연한 것이지만, 평균화 단계에 있어서 사용하는 연산은, 투영 데이터 세트에 포함되는 모든 화소에 대해 전부 동일한 것을 사용한다. 따라서, 화소 a1과 b1 사이에서 기하 평균에 의해 화소값의 평균화 처리를 실시한 경우는, a2와 b2 사이 등, 다른 대응 화소 사이에 대해서도, 동일하게 기하 평균에 의해 평균화 처리를 실시한다.
이상, 본 발명에 따른 윤곽 화상 형성 장치(10)의 바람직한 양태에 있어서의 구성과 동작에 대해 설명하였다.
다음으로, 본 발명에 따른 윤곽 화상 형성 프로그램(300)에 대해 설명한다. 상술한 바와 같이, 본 발명에 따른 윤곽 화상 형성 장치(10)는, 본 발명에 따른 윤곽 화상 형성 프로그램(300)을 판독 입력한 컴퓨터로서 구성할 수 있다.
도 6에, 본 발명에 따른 윤곽 화상 형성 프로그램(300)의 바람직한 양태에 있어서의 구성을 나타낸다. 바람직한 양태에 있어서, 윤곽 화상 형성 프로그램(300)은, 처리를 통괄하는 메인 모듈(310)과, 데이터 취득 모듈(320)과, 사이노그램 작성 모듈(330)과, 제1 2치화 처리 모듈(340)과, 윤곽 추출 모듈(350)과, 데이터 보간 처리 모듈(360)과, 화상 재구성 모듈(370)과, 제2 2치화 처리 모듈(380)과, 출력 모듈(390)에 의해 구성된다.
데이터 취득 모듈(320)은, 단계 S01에 따른 처리를 컴퓨터에 실행시킨다.
사이노그램 작성 모듈(330)은, 단계 S02와 관련된 처리를 컴퓨터에 실행시킨다.
제1 2치화 처리 모듈(340)은, 단계 S03에 따른 처리를 컴퓨터에 실행시킨다.
윤곽 추출 모듈(350)은, 단계 S04와 관련된 처리를 컴퓨터에 실행시킨다.
데이터 보간 처리 모듈(360)은, 단계 S05와 관련된 처리를 컴퓨터에 실행시킨다.
화상 재구성 모듈(370)은, 단계 S06에 따른 처리를 컴퓨터에 실행시킨다.
제2 2치화 처리 모듈(380)은, 단계 S07에 따른 처리를 컴퓨터에 실행시킨다.
출력 모듈(390)은, 단계 S08에 따른 처리를 컴퓨터에 실행시킨다.
윤곽 화상 형성 프로그램(300)은, 화소값 평균화 모듈(도시하지 않음)을 추가로 구비한 것이어도 된다. 화소값 평균화 모듈은, 컴퓨터에 의해 판독 입력됨으로써, 컴퓨터를 상술한 화소값 평균화부로서 기능시켜, 화소값 평균화 단계를 실행시킨다.
본 발명의 실시에 의해 얻어진 윤곽 화상은, μ맵을 형성하기 위해서 바람직하게 이용할 수 있다. 예를 들면 본 발명의 실시에 의해 얻어진 윤곽 화상에 있어서 0 이외의 화소값을 갖는 부분을 임의의 흡수 계수로 치환함으로써, 간단하게 μ맵을 형성하는 것이 가능해진다. 본 발명의 실시에 의해 얻어진 윤곽 화상은 정확하게 피험자 화상의 윤곽을 추적하고 있으므로, 본 윤곽 화상을 이용하여 작성된 μ맵에 의해, 보다 양호한 정밀도로 산란 보정 및 흡수 보정을 실시하는 것이 가능하다.
실시예
I-123 IMP(니혼 메디피직스 주식회사 제조) 투여에 의한 두부 SPECT의 투영 데이터 세트(투영 각도 피치：4°, 투영 데이터 장수：90장, 투영 각도 범위：0∼360°)를 이용하였다. 투영 데이터 세트는, 좌우 방향에 표시된 각 투영 각도에 있어서의 투영상을 투영 각도순으로 세로 방향으로 배열하여 작성한, 이른바 사이노그램의 형태로 이용하였다(즉, 사이노그램 작성 단계 처리가 끝난 데이터).
상기 투영 데이터 세트에 대해, 화소값의 최고값의 25％ 역치로 한 제1 2치화 처리를 실시하였다. 2치화 후의 투영 데이터 세트를 이용하여 360°역투영을 실시해, 화상의 재구성을 실시하였다. 얻어진 재구성 화상에 대해, 화소값의 최대값의 50％를 역치로 하여 제2 2치화 처리를 실시하여, 윤곽 화상을 작성하였다(도 7a). 얻어진 윤곽 화상은, 이용한 투영 데이터 세트에 의해 얻어진 재구성 화상(도 7b)의 윤곽을 양호하게 재현할 수 있었다. 이상의 결과로부터, 본 발명에 따른 방법에 의해, 윤곽 화상을 양호한 정밀도로 재현할 수 있는 것이 확인되었다.
본 발명은 화상 처리 소프트웨어 및 화상 진단 기기의 분야에서 이용할 수 있다.
10 윤곽 화상 형성 장치
20 데이터 취득부
30 사이노그램 작성부
40 제1 2치화 처리부
50 윤곽 추출부
60 데이터 보간 처리부
70 화상 재구성부
80 제2 2치화 처리부
90 출력부
100 핵의학 화상 촬상 장치
120 투영상 A
130 투영상 B
150 피험자
230 CPU
240 메모리
250 모니터
260 통신 인터페이스
270 키보드
280 버스
300 윤곽 화상 형성 프로그램
310 메인 모듈
320 데이터 취득 모듈
330 사이노그램 작성 모듈
340 제1 2치화 처리 모듈
350 윤곽 추출 모듈
360 데이터 보간 처리 모듈
370 화상 재구성 모듈
380 제2 2치화 처리 모듈
390 출력 모듈

Claims (13)

  1. 피험자의 투영 데이터로부터 산란 및 흡수 보정을 실시하지 않고 윤곽 화상을 생성하는 방법으로서,
    수집 각도 범위 0∼180°이상의 핵의학 화상의 투영 데이터 세트를 취득하는 데이터 취득 단계와,
    취득한 상기 투영 데이터 세트로부터 사이노그램을 작성하는 사이노그램 작성 단계와,
    작성한 사이노그램에 대해 제1 2치화 처리를 실시하는 제1 2치화 단계와,
    상기 제1 2치화 처리 후의 상기 사이노그램을 역투영하여 화상을 재구성하는 화상 재구성 단계와,
    재구성된 화상에 대해 제2 2치화 처리를 실시하는 제2 2치화 단계를 컴퓨터에 의해 순차 실행하는 것을 특징으로 하는 윤곽 화상 형성 방법.
  2. 제 1 항에 있어서,
    제1 2치화 단계 후의 사이노그램에 대해 데이터 보간을 실시하는 보간 단계를 추가로 포함하고,
    화상 재구성 단계는, 보간 단계 완료 후의 사이노그램을 이용하여 역투영을 실시하는 것인 윤곽 화상 형성 방법.
  3. 제 1 항 또는 제 2 항에 있어서,
    투영 데이터 세트에 있어서의 수집 각도 범위의 상한값이 360°이상인 윤곽 화상 형성 방법.
  4. 제 3 항에 있어서,
    데이터 취득 단계 후, 사이노그램 작성 단계 전에 있어서, 수집 각도 범위 0∼180°까지의 범위에 있어서의 투영 데이터 세트를 구성하는 각각의 투영 데이터의 각 화소에 대해, 당해 투영 데이터와 180°대향하는 투영 데이터에 있어서의 대응하는 화소 사이에서 화소값의 합, 기하 평균, 또는 산술 평균을 산출하고, 얻어진 값을 당해 투영 데이터에 있어서의 화소값과 치환하는 평균화 단계를 추가로 실행하고, 여기서, 평균화 단계에서 실시하는 연산은 투영 데이터 세트에 포함되는 모든 화소에 대해 전부 동일한 것이고,
    사이노그램 작성 단계는 평균화 단계 완료 후의 수집 각도 0°∼180°까지의 투영 데이터 세트를 이용하여 실시하는 것인 윤곽 화상 형성 방법.
  5. 제 3 항에 있어서,
    사이노그램 작성 단계 후, 제1 2치화 단계 전에 있어서, 수집 각도 범위 0∼180°까지의 범위에 있어서의 사이노그램을 구성하는 각각의 투영 데이터의 각 화소에 대해, 당해 투영 데이터와 180°대향하는 투영 데이터에 있어서의 대응하는 화소 사이에서 화소값의 합, 기하 평균, 또는 산술 평균을 산출하고, 얻어진 값을 당해 투영 데이터에 있어서의 화소값과 치환하는 평균화 단계를 추가로 실행하고, 또한, 수집 각도 180°이상의 투영 데이터를 당해 사이노그램으로부터 삭제하고, 여기서, 평균화 단계에서 실시하는 연산은 투영 데이터 세트에 포함되는 모든 화소에 대해 전부 동일한 것이고,
    제1 2치화 단계는 평균화 단계 완료 후의 사이노그램을 이용하여 실시하는 것인 윤곽 화상 형성 방법.
  6. 컴퓨터에 판독 입력되어 당해 컴퓨터를 피험자의 투영 데이터로부터 산란 및 흡수 보정을 실시하지 않고 윤곽 화상을 생성하는 윤곽 화상 생성 장치로서 동작시키기 위한 프로그램을 기록한 컴퓨터 판독 가능한 기록 매체로서,
    컴퓨터에,
    수집 각도 범위 0∼180°이상의 핵의학 화상의 투영 데이터 세트를 취득하는 데이터 취득 단계와,
    취득한 상기 투영 데이터 세트로부터 사이노그램을 작성하는 사이노그램 작성 단계와,
    상기 사이노그램에 대해 제1 2치화 처리를 실시하는 제1 2치화 단계와,
    제1 2치화 후의 상기 사이노그램을 역투영하여 화상을 재구성하는 화상 재구성 단계와,
    재구성된 화상에 대해 제2 2치화 처리를 실시하는 제2 2치화 단계를 순차 실행시키는 것을 특징으로 하는 윤곽 화상 형성 프로그램을 기록한 컴퓨터 판독 가능한 기록 매체.
  7. 제 6 항에 있어서,
    컴퓨터에,
    제1 2치화 단계 후의 사이노그램에 대해 데이터 보간을 실시하는 보간 단계를 추가로 실행시키고,
    화상 재구성 단계는, 보간 단계 완료 후의 사이노그램을 이용하여 역투영을 실시하는 것인 윤곽 화상 형성 프로그램을 기록한 컴퓨터 판독 가능한 기록 매체.
  8. 제 6 항 또는 제 7 항에 있어서,
    투영 데이터 세트에 있어서의 수집 각도 범위의 상한값이 360°이상인 윤곽 화상 형성 프로그램을 기록한 컴퓨터 판독 가능한 기록 매체.
  9. 제 8 항에 있어서,
    데이터 취득 단계 후, 사이노그램 작성 단계 전에 있어서, 컴퓨터에 수집 각도 범위 0∼180°까지의 범위에 있어서의 투영 데이터 세트를 구성하는 각각의 투영 데이터의 각 화소에 대해, 당해 투영 데이터와 180°대향하는 투영 데이터에 있어서의 대응하는 화소 사이에서 화소값의 합, 기하 평균, 또는 산술 평균을 산출하고, 얻어진 값을 당해 투영 데이터에 있어서의 화소값과 치환하는 평균화 단계를 추가로 실행시키고, 여기서, 평균화 단계에서 실시하는 연산은 투영 데이터 세트에 포함되는 모든 화소에 대해 전부 동일한 것이고,
    컴퓨터에 평균화 단계 완료 후의 수집 각도 0°∼180°까지의 투영 데이터 세트를 이용하여 사이노그램 작성 단계를 실행시키는 것인 윤곽 화상 형성 프로그램을 기록한 컴퓨터 판독 가능한 기록 매체.
  10. 제 8 항에 있어서,
    사이노그램 작성 단계 후, 제1 2치화 단계 전에 있어서, 컴퓨터에 수집 각도 범위 0∼180°까지의 범위에 있어서의 사이노그램을 구성하는 각각의 투영 데이터의 각 화소에 대해, 당해 투영 데이터와 180°대향하는 투영 데이터에 있어서의 대응하는 화소 사이에서 화소값의 합, 기하 평균, 또는 산술 평균을 산출하고, 얻어진 값을 당해 투영 데이터에 있어서의 화소값과 치환하는 평균화 단계를 추가로 실행시키고, 또한, 수집 각도 180°이상의 투영 데이터를 당해 사이노그램으로부터 삭제하고, 여기서, 평균화 단계에서 실시하는 연산은 투영 데이터 세트에 포함되는 모든 화소에 대해 전부 동일한 것이고,
    컴퓨터에 평균화 단계 완료 후의 사이노그램을 이용하여 제1 2치화 단계를 실행시키는 것인 윤곽 화상 형성 프로그램을 기록한 컴퓨터 판독 가능한 기록 매체.
  11. 피험자의 투영 데이터로부터 산란 및 흡수 보정을 실시하지 않고 윤곽 화상을 생성하기 위한 윤곽 화상 형성 장치로서,
    핵의학 화상의 투영 데이터 세트를 취득하는 기능을 갖는 데이터 취득부와,
    투영 데이터 세트로부터 사이노그램을 작성하는 사이노그램 작성부와,
    사이노그램에 대해 데이터의 2치화 처리를 실시하는 기능을 갖는 제1 2치화 처리부와,
    사이노그램을 역투영하여 화상을 재구성하는 기능을 갖는 화상 재구성 처리부와,
    재구성된 화상에 대해 화상의 2치화 처리를 실시하는 기능을 갖는 제2 2치화 처리부를 구비한 것을 특징으로 하는 윤곽 화상 형성 장치.
  12. 제 11 항에 있어서,
    사이노그램의 데이터 보간을 실시하는 기능을 갖는 보간 처리부를 추가로 포함한 것인 윤곽 화상 형성 장치.
  13. 제 11 항 또는 제 12 항에 있어서,
    투영 데이터 세트 또는 사이노그램을 구성하는 각각의 투영 데이터의 각 화소에 대해, 당해 투영 데이터와 180°대향하는 투영 데이터에 있어서의 대응하는 화소 사이에서 화소값의 합, 기하 평균, 또는 산술 평균을 산출하고, 얻어진 값을 당해 투영 데이터에 있어서의 화소값과 치환하는 기능을 갖는 화소값 평균화부를 추가로 구비한 것이고, 여기서, 평균화부에 있어서 사용하는 연산은 투영 데이터 세트에 포함되는 모든 화소에 대해 전부 동일한 것을 사용하는 것인 윤곽 화상 형성 장치.
KR1020130047979A 2012-05-07 2013-04-30 단층 화상의 윤곽 추출 방법, 윤곽 화상 형성 프로그램을 기록한 컴퓨터 판독 가능한 기록 매체, 및 윤곽 화상 형성 장치 KR20130124895A (ko)

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
EP12167018.6 2012-05-07
EP12167018.6A EP2662022A1 (en) 2012-05-07 2012-05-07 Method of extracting contour of tomogram
JP2012134231 2012-05-29
JPJP-P-2012-134231 2012-05-29

Publications (1)

Publication Number Publication Date
KR20130124895A true KR20130124895A (ko) 2013-11-15

Family

ID=49512555

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020130047979A KR20130124895A (ko) 2012-05-07 2013-04-30 단층 화상의 윤곽 추출 방법, 윤곽 화상 형성 프로그램을 기록한 컴퓨터 판독 가능한 기록 매체, 및 윤곽 화상 형성 장치

Country Status (4)

Country Link
US (1) US9202293B2 (ko)
JP (1) JP2014006246A (ko)
KR (1) KR20130124895A (ko)
CA (1) CA2812998A1 (ko)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103908278B (zh) * 2013-01-06 2017-07-28 上海联影医疗科技有限公司 图像重建方法及装置、医疗成像系统
US10188358B2 (en) * 2014-05-15 2019-01-29 General Electric Company System and method for subject shape estimation
JP2016180649A (ja) * 2015-03-24 2016-10-13 日本メジフィジックス株式会社 画像処理装置、画像処理方法、及びプログラム
US10078889B2 (en) 2015-08-25 2018-09-18 Shanghai United Imaging Healthcare Co., Ltd. System and method for image calibration
WO2020006032A1 (en) * 2018-06-26 2020-01-02 The Medical College Of Wisconsin, Inc. Systems and methods for accelerated online adaptive radiation therapy

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2169180B (en) * 1984-12-28 1988-06-15 Toshiba Kk Ct apparatus and operating method therefor
US4803639A (en) * 1986-02-25 1989-02-07 General Electric Company X-ray inspection system
JP3925278B2 (ja) * 2002-04-04 2007-06-06 株式会社島津製作所 核医学イメージング装置
GB0324374D0 (en) 2003-10-17 2003-11-19 Hammersmith Imanet Ltd Method of, and software for, conducting motion correction for a tomographic scanner
US7254209B2 (en) * 2003-11-17 2007-08-07 General Electric Company Iterative CT reconstruction method using multi-modal edge information
US7203267B2 (en) * 2004-06-30 2007-04-10 General Electric Company System and method for boundary estimation using CT metrology
JP2006068383A (ja) * 2004-09-03 2006-03-16 Hitachi Medical Corp 画像処理装置
US7176916B2 (en) 2005-04-15 2007-02-13 T.I.E.S., Inc. Object identifying system for segmenting unreconstructed data in image tomography
JP2007248121A (ja) * 2006-03-14 2007-09-27 Fujita Gakuen 断層画像の輪郭抽出方法、プログラム、および装置
US7778386B2 (en) * 2006-08-28 2010-08-17 General Electric Company Methods for analytic reconstruction for mult-source inverse geometry CT
US8284892B2 (en) * 2008-12-22 2012-10-09 General Electric Company System and method for image reconstruction
US8620053B2 (en) * 2009-11-04 2013-12-31 Siemens Medical Solutions Usa, Inc. Completion of truncated attenuation maps using maximum likelihood estimation of attenuation and activity (MLAA)
US8989462B2 (en) * 2010-02-11 2015-03-24 Emory University Systems, methods and computer readable storage mediums storing instructions for applying multiscale bilateral filtering to magnetic resonance (RI) images
JP4937397B2 (ja) * 2010-10-25 2012-05-23 富士フイルム株式会社 医用画像診断支援装置および方法、並びにプログラム

Also Published As

Publication number Publication date
US20130294671A1 (en) 2013-11-07
US9202293B2 (en) 2015-12-01
JP2014006246A (ja) 2014-01-16
CA2812998A1 (en) 2013-11-07

Similar Documents

Publication Publication Date Title
Kapucu et al. EANM procedure guideline for brain perfusion SPECT using 99m Tc-labelled radiopharmaceuticals, version 2
Harms et al. Automatic generation of absolute myocardial blood flow images using [15 O] H 2 O and a clinical PET/CT scanner
CN105120738B (zh) 狭窄治疗规划
CN110599472B (zh) 计算spect定量断层图像中suv归一化系数的方法及系统
KR20080081005A (ko) 두부 변성 질환의 검출 방법, 검출 프로그램, 및 검출 장치
KR20130124895A (ko) 단층 화상의 윤곽 추출 방법, 윤곽 화상 형성 프로그램을 기록한 컴퓨터 판독 가능한 기록 매체, 및 윤곽 화상 형성 장치
EP2312337A1 (en) Technique for detecting cranial nerve disease
JP2019537035A5 (ko)
Germano et al. Attenuation correction in cardiac SPECT: the boy who cried wolf?
CN110996800A (zh) 用于确定pet成像动力学参数的系统、方法
JP4997575B2 (ja) 減弱係数マップ作成装置、減弱係数マップ作成方法およびプログラム
Abdelmoneim et al. Assessment of Myocardial Perfusion during Adenosine Stress Using Real Time Three‐Dimensional and Two‐Dimensional Myocardial Contrast Echocardiography: Comparison with Single‐Photon Emission Computed Tomography
JP5930193B2 (ja) Ｓｐｅｃｔ画像の画像変換装置、画像変換プログラム、および画像変換方法
US9786069B2 (en) Refined reconstruction of time-varying data
JP2007248121A (ja) 断層画像の輪郭抽出方法、プログラム、および装置
CN106662657A (zh) 头部核医学图像的评价方法
Piccinelli et al. Multimodality Image Fusion for Coronary Artery Disease Detection Concepts and Latest Developments
JP2013063324A5 (ko)
EP2662022A1 (en) Method of extracting contour of tomogram
JP6995642B2 (ja) 画像処理装置、画像処理方法およびプログラム
Okwuosa et al. Left ventricular mass from gated SPECT myocardial perfusion imaging: comparison with cardiac computed tomography
JPWO2009050837A1 (ja) 画像処理プログラム、画像処理方法および画像処理装置
JP4494510B1 (ja) 核医学画像の画像処理方法、画像処理装置、及び画像処理プログラム
TW201923381A (zh) 影像處理程式製品、影像處理裝置及影像處理方法
KR102698910B1 (ko) 의료 영상 처리 방법

Legal Events

Date Code Title Description
WITN Application deemed withdrawn, e.g. because no request for examination was filed or no examination fee was paid