KR101689866B1 - 영상 처리 방법 및 장치와 이를 채용한 의료영상시스템 - Google Patents

영상 처리 방법 및 장치와 이를 채용한 의료영상시스템 Download PDF

Info

Publication number
KR101689866B1
KR101689866B1 KR1020100073697A KR20100073697A KR101689866B1 KR 101689866 B1 KR101689866 B1 KR 101689866B1 KR 1020100073697 A KR1020100073697 A KR 1020100073697A KR 20100073697 A KR20100073697 A KR 20100073697A KR 101689866 B1 KR101689866 B1 KR 101689866B1
Authority
KR
South Korea
Prior art keywords
image
materials
thickness
total thickness
subject
Prior art date
Application number
KR1020100073697A
Other languages
English (en)
Other versions
KR20120011679A (ko
Inventor
강동구
한석민
이성덕
성영훈
김성수
오현화
권재현
Original Assignee
삼성전자주식회사
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 삼성전자주식회사 filed Critical 삼성전자주식회사
Priority to KR1020100073697A priority Critical patent/KR101689866B1/ko
Priority to EP11812809.9A priority patent/EP2598031A4/en
Priority to CN2011800372808A priority patent/CN103037768A/zh
Priority to PCT/KR2011/005631 priority patent/WO2012015280A2/en
Priority to US13/813,007 priority patent/US9008351B2/en
Priority to JP2013521719A priority patent/JP5536284B2/ja
Publication of KR20120011679A publication Critical patent/KR20120011679A/ko
Application granted granted Critical
Publication of KR101689866B1 publication Critical patent/KR101689866B1/ko

Links

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/48—Diagnostic techniques
    • A61B6/482—Diagnostic techniques involving multiple energy imaging
    • 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/502—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 diagnosis of breast, i.e. mammography
    • 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/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/50—Depth or shape recovery
    • G06T7/529—Depth or shape recovery from texture
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • 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/10116—X-ray image
    • 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/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30068—Mammography; Breast

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Surgery (AREA)
  • Molecular Biology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Biophysics (AREA)
  • Veterinary Medicine (AREA)
  • Optics & Photonics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Quality & Reliability (AREA)
  • Physiology (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Primary Health Care (AREA)
  • Epidemiology (AREA)
  • Databases & Information Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)

Abstract

영상 처리 방법 및 장치와 이를 채용한 의료영상시스템이 개시되어 있다. 영상 처리 장치에서 전체 두께 추정부는 적어도 두가지 물질로 이루어지는 소정 피사체로부터 얻어지는 적어도 두가지 에너지 대역의 방사선 영상으로부터 상기 피사체의 전체 두께를 추정한다. 진단영상 생성부는 추정된 전체 두께와 방사선 영상에서 국부 영역의 두께를 비교하고, 비교결과에 따라서 관심 영역을 추출하여 진단영상을 생성한다.

Description

영상 처리 방법 및 장치와 이를 채용한 의료영상시스템{Method and apparatus of processing image and medical image system employing the same}
영상 처리방법 및 장치와 이를 채용한 의료영상시스템에 관한 것이다.
일반적으로 소정 피사체에 대한 방사선 영상 예를 들면 엑스-레이(X-ray) 영상은 엑스-레이의 투과도에 의하여 결정된다. 엑스-레이의 투과도는 피사체의 종류, 밀도, 엑스-레이의 에너지 대역에 따라서 달라진다. 따라서, 한장의 엑스-레이 영상만으로는 피사체의 종류를 구분하는 것이 어렵고, 특히 엑스-레이 영상에 포함된 정상 물질과 비정상 물질을 구분하는 것이 어렵다.
복수개의 물질이 혼재된 다중 에너지 방사선 영상으로부터 비정상 조직을 추출하고, 추출된 비정상 조직으로 이루어진 진단영상을 얻는 영상 처리 방법 및 장치를 제공한다.
상기 영상 처리 장치를 채용한 의료영상시스템을 제공한다.
본 발명의 한 측면(aspect)에 따르면, 적어도 두가지 물질로 이루어지는 소정 피사체로부터 얻어지는 적어도 두가지 에너지 대역의 방사선 영상으로부터 상기 피사체의 전체 두께를 추정하는 단계; 및 상기 추정된 전체 두께와 상기 방사선 영상에서 국부 영역의 두께를 비교하고, 비교결과에 따라서 관심 영역을 추출하여 진단영상을 생성하는 단계를 포함하는 영상 처리 방법이 제공된다.
본 발명의 다른 측면에 따르면, 적어도 두가지 물질로 이루어지는 소정 피사체로부터 얻어지는 적어도 두가지 에너지 대역의 방사선 영상으로부터 상기 피사체의 전체 두께를 추정하는 전체 두께 추정부; 및 상기 추정된 전체 두께와 상기 방사선 영상에서 국부 영역의 두께를 비교하고, 비교결과에 따라서 관심 영역을 추출하여 진단영상을 생성하는 진단영상 생성부를 포함하는 영상 처리 장치가 제공된다.
여기서, 상기 관심 영역은 비정상 조직이고, 상기 두가지 물질은 정상 조직의 구성물질일 수 있다.
여기서, 상기 전체 두께 추정시 두가지 물질 조합을 갖는 연속두께 팬텀을 이용하여 얻어는 전체 두께 모델을 적용하여 수행될 수 있다.
여기서, 상기 전체 두께 모델은 상기 두가지 물질 중 어느 하나의 물질의 밀도를 파라미터로 포함할 수 있다.
여기서, 상기 전체 두께 추정시 상기 두가지 물질의 밀도를 추정하고, 상기 추정된 두가지 물질의 밀도와 각 물질의 감쇠 베이시스를 이용하여 상기 다중 에너지 방사선 영상에서 상기 두가지 물질을 분리하고, 상기 분리된 두가지 물질의 두께를 더하여 상기 피사체의 전체 두께를 추정할 수 있다.
본 발명의 또 다른 측면에 따르면, 상기한 영상 처리장치를 채용하는 의료영상시스템이 제공된다.
여기서, 상기 의료영상시스템은 상기 소정 피사체에 적어도 두가지 에너지 대역의 방사선을 여기서 조사하여 상기 방사선 영상을 획득하는 방사선 영상 촬영부를 더 포함할 수 있다.
여기서, 상기 의료영상시스템은 상기 생성된 진단영상을 저장하거나, 상기 생성된 진단영상으로부터 얻어지는 진단정보를 상기 진단영상에 대응시켜 저장하는 저장부를 더 포함할 수 있다.
여기서, 상기 의료영상시스템은 상기 생성된 진단영상을 전송하거나, 상기 생성된 진단영상으로부터 얻어지는 진단정보를 상기 진단영상에 대응시켜 전송하는 통신부를 더 포함할 수 있다.
복수개의 물질이 혼재된 다중 에너지 방사선 영상에 대하여 관심영역 예를 들면 비정상 조직을 추출하고, 추출된 비정상 조직을 강조 처리하여 선명도가 높은 진단영상을 얻음으로써, 한장의 방사선 영상을 사용하더라도 진단의 정확도를 높일 수 있다.
도 1은 정상 물질과 비정상 물질의 에너지대역별 엑스-레이 감쇠 계수를 나타낸 그래프이다.
도 2는 인체 구성 물질의 감쇠 베이시스 모델을 나타낸다.
도 3은 3가지 물질을 분리할 때 피사체의 전체 두께에 대한 정보를 이용하는 방법의 개념을 설명하는 도면이다.
도 4는 본 발명의 실시예에 따른 영상 처리 방법을 설명하는 흐름도이다.
도 5는 총분산을 이용하여 유방의 전체 두께의 균일성을 최대화하는 예를 보여주는 도면이다.
도 6은 두가지 물질 조합을 갖는 연속두께 팬텀을 이용하여 하나의 파라미터를 갖는 전체 두께 모델을 설정하는 방법의 개념을 보여주는 도면이다.
도 7은 본 발명의 실시예에 따른 의료영상시스템의 구성을 보여주는 블럭도이다.
도 8은 본 발명의 실시예에 따른 영상 처리 장치의 구성을 보여주는 블럭도이다.
이하, 본 발명의 실시예에 따른 영상 처리 방법 및 장치와 이를 채용한 의료영상시스템을 첨부된 도면을 참조하여 상세하게 설명한다.
먼저, 방사선 영상 예를 들면 엑스-레이 영상과 관련하여 입사 세기와 투과 세기와의 관계를 설명하기로 한다.
엑스-레이의 투과도는 피사체의 종류, 밀도, 엑스-레이의 에너지 대역에 따라서 달라지며, 입사 세기(incident intensity)와 투과 세기(transmitted intensity)는 다음 수학식 1과 같이 나타낼 수 있다.
Figure 112010049307952-pat00001
여기서, I(E)는 투과 세기, I0(E)는 입사 세기를 각각 나타낸다.
Figure 112010049307952-pat00002
는 에너지대역(E)에서의 질량 감쇠 계수(mass attenuation coefficient)를 나타내고,
Figure 112010049307952-pat00003
는 피사체의 밀도를 나타내고, x는 피사체의 두께를 나타낸다. 상기 수학식 1은 Beer-Lambert Law라고도 일컫는다.
도 1은 정상 물질과 비정상 물질의 에너지대역별 엑스-레이 감쇠 계수를 나타낸 그래프이다. 도 1을 참조하면, 에너지대역별로 지방조직(Fat), 선상조직 (Glandular), 침윤성 유관상피암 조직(IDC; Infiltrating Ductal Carcinoma)의 엑스-레이 감쇠 계수가 서로 상이함을 알 수 있다.
한편, 적어도 두가지 이상의 에너지 대역의 방사선, 예를 들면 엑스-레이를 사용하는 경우, 각 에너지 대역별 영상(I1, I2,..., In, 여기서 n은 에너지대역)은 다음 수학식 2와 같이 나타낼 수 있다.
Figure 112010049307952-pat00004
여기서, wn은 에너지대역(n)의 엑스-레이 입사 세기,
Figure 112010049307952-pat00005
로서 각 위치벡터(r)에서의 물질 밀도 투사영상(material density projections)를 각각 나타낸다.
상기한 수학식 2에서 각 에너지 대역별 영상(I1, I2,..., In, 여기서 n은 에너지대역)가 얻어질 경우, 각 에너지 대역별 영상으로부터 물질 밀도 투사영상(Ci)를 구함으로써, 정상 물질과 비정상 물질을 분리할 수 있다.
한편, 인체의 구성 물질을 가정하면, 인체 구성 물질의 감쇠 베이시스의 수가 두개 즉, 광전흡수(photoelectric absorption) 및 컴프턴 산란(Compton scattering)이기 때문에, 최대 두가지의 물질이 섞여있는 경우 각 물질을 분리할 수 있다.
도 2는 인체 구성 물질 예를 들면 물(water))의 감쇠 베이시스 모델을 나타낸다. 도 2를 참조하면, 물리적 모델은 다음 수학식 3
Figure 112010049307952-pat00006
과 같이 나타낼 수 있다.
한편, 실험적 모델은 다음 수학식 4
Figure 112010049307952-pat00007
과 같이 나타낼 수 있다.
상기 수학식 3 및 수학식 4에 있어서,
Figure 112010049307952-pat00008
는 에너지대역 E에서 위치벡터 r의 감쇠 계수,
Figure 112010049307952-pat00009
,
Figure 112010049307952-pat00010
는 각 조직의 위치벡터 r의 길이(두께),
Figure 112010049307952-pat00011
,
Figure 112010049307952-pat00012
는 에너지대역(E)의 광전흡수 및 컴프턴 산란에 의한 베이시스 함수 즉 감쇠 계수,
Figure 112010049307952-pat00013
,
Figure 112010049307952-pat00014
는 에너지대역(E)의 물과 뼈의 감쇠 계수를 각각 나타낸다.
한편, 인체의 장기 중 예를 들어 유방(breast)을 구성하는 정상 조직은 크게 지방조직(adipose tissue)과 선상조직(glandular tissue)으로 구성되고, 종괴(mass)와 같은 비정상 조직이 함께 섞여 있을 경우, 다중 에너지 방사선 영상에는 총 3가지 물질이 존재하게 된다. 이러한 경우, 인체 구성 물질의 감쇠 베이시스인 광전 흡수와 컴프턴 산란에 피사체의 전체 두께 정보를 포함시킴으로써, 3가지 물질을 분리하는 것이 가능하다.
구체적으로, 인체 구성 물질의 감쇠 베이시스인 광전 흡수와 컴프턴 산란을 이용하는 경우 에너지대역 E에서 위치벡터 r의 감쇠 계수(
Figure 112010049307952-pat00015
)는 다음 수학식 5
Figure 112010049307952-pat00016
Figure 112010049307952-pat00017
와 같이 나타낼 수 있으므로, 지방조직(A), 선상조직(G), 비정상조직(C)과 같은 3가지 물질을 분리하는 것이 어렵다. 여기서,
Figure 112010049307952-pat00018
,
Figure 112010049307952-pat00019
,
Figure 112010049307952-pat00020
는 각 조직의 두께를 나타내는 비례상수이고,
Figure 112010049307952-pat00021
,
Figure 112010049307952-pat00022
,
Figure 112010049307952-pat00023
는 각각 지방조직(A), 선상조직(G), 비정상조직(C)의 감쇠 계수를 나타낸다. α,β는 각각 비례상수를 나타내고, fph(E)는 에너지대역(E)에서의 광전흡수에 의한 베이시스 함수 즉, 감쇠계수를 나타내고, fCo(E)는 에너지대역(E)에서의 컴프턴 산란에 의한 베이시스 함수 즉, 감쇠계수를 나타낸다.
한편, 피사체의 전체 두께(T)는 다음 수학식 6
Figure 112010049307952-pat00024
과 같이 나타낼 수 있다.
상기한 바와 같은 피사체의 전체 두께(T)를 이용하게 되면, 에너지대역 E에서 위치벡터 r의 감쇠 계수(
Figure 112010049307952-pat00025
)는 다음 수학식 7
Figure 112010049307952-pat00026
과 같이 나타낼 수 있으므로, 3가지 물질의 분리가 가능해진다.
도 3은 3가지 물질을 분리할 때 피사체의 전체 두께에 대한 정보를 이용하는 방법의 개념을 설명하는 도면이다. 이는 피사체 예를 들면 압착 유방의 전체 두께는 거의 일정하다는 사실을 이용한 것이다. 구체적으로, 선상 조직(G) 및 지방 조직(A)과 같은 정상 물질의 감쇠 베이시스를 이용하여, 방사선 영상에서 선상 조직과 지방 조직을 분리한다. 분리된 선상 조직과 지방 조직을 더하면 전체 두께 영상을 얻을 수 있다. 즉, 전체 두께 영상에서 정상 물질의 영역에서는 피사체 즉, 유방의 전체 두께와 동일한 값으로 표현되는 반면, 비정상 물질의 영역에서는 정상 물질의 두께와는 다른 값으로 표현된다.
도 3을 참조하면, 지방조직(A)과 선상조직(G)으로 이루어지는 영역(310)과 지방조직(A), 선상조직(C)과 비정상조직인 암 조직(C)으로 이루어지는 영역(320)의 그레이 스케일이 각각 세기(330) 및 세기(340)로 나타날 경우, 지방조직(A)과 선상조직(G)에 대한 감쇠 베이시스를 이용하여, 그레이 스케일 세기(330) 및 그레이 스케일 세기(340)로부터 각각 지방조직(A)과 선상조직(G)을 분리한다. 두께(350)는 그레이 스케일 세기(330)으로부터 분리된 지방조직(A)과 선상조직(G)을, 두께(360)은 그레이 스케일 세기(330)으로부터 분리된 지방조직(A)과 선상조직(G)을 나타낸다. 이와 같이 분리된 선상 조직과 지방 조직을 더하면 정상 물질의 영역(310)에 대해서는 유방의 전체 두께와 동일한 두께(370)가 되고, 비정상 물질의 영역(320)에 대해서는 유방의 전체 두께와는 다른 두께(380)가 된다.
도 4는 본 발명의 실시예에 따른 영상 처리 방법을 설명하는 흐름도이다.
도 4를 참조하면, 410 단계에서는 소정의 피사체에 대한 적어도 두가지 에너지 대역의 다중 에너지 방사선 영상을 입력한다.
430 단계에서는 다중 에너지 방사선 영상을 이루는 적어도 두가지 물질의 밀도를 추정한다.
450 단계에서는 추정된 각 물질의 밀도와 해당 물질의 감쇠 베이시스를 이용하여 다중 에너지 방사선 영상으로부터 두가지 물질 예를 들면, 선상조직과 지방조직을 분리한다.
470 단계에서는 분리된 선상 조직과 지방 조직을 더하여 피사체의 전체 두께를 추정한다. 한편, 도시되지는 않았으나, 470 단계에서 추정된 전체 두께에 대한 정보를 이용하여 피사체의 균일도를 측정하고, 측정된 균일도에 근거하여 물질의 밀도와 같은 파라미터를 업데이트시켜 430 단계로 피드백한다.
490 단계에서는 추정된 전체 두께와 다중 에너지 방사선 영상의 국부영역에서 산출된 두께를 비교하여 종괴(mass) 혹은 미세석회화(microcalcification) 조직과 같은 비정상 영역을 추출한다. 추출된 비정상 영역을 이용하여 진단영상을 생성한다.
한편, 상기 430 단계, 450 단계 및 470 단계를 전체 두께 모델로 치환하여 일정한 개수의 파라미터를 갖는 다항식으로 표현하고, 각 파라미터를 전체 두께의 균일성을 최대화할 수 있도록 최적화시킴으로써 전체 두께 영상을 얻을 수도 있다. 여기서, 1차 다항식으로 표현된 전체 두께 모델은 다음 수학식 8
Figure 112010049307952-pat00027
과 같이 나타낼 수 있고, 2차 다항식으로 표현된 전체 두께 모델은 다음 수학식 9
Figure 112010049307952-pat00028
와 같이 나타낼 수 있다.
상기 수학식 8 및 수학식 9에 있어서, IL 및 IH는 각각 낮은 에너지대역과 높은 에너지대역의 방사선 영상을 나타내고, a1, a2, a3, a4, a5, a6는 각각 최적화된 비례상수를 나타낸다.
도 5는 총분산(total variation)을 이용하여 유방의 전체 두께의 균일성을 최대화하는 예를 보여주는 도면이다. 도 5를 참조하면, 유방에 대한 방사선 영상(510)에서 총분산(L1-norm)을 최소화시켜 1차 다항식 영상(530)으로 표현함으로써 유방의 전체 두께의 균일성을 최대화시킬 수 있다.
한편, 피사체 즉, 유방의 전체 두께는 상기한 수학식 8 및 수학식 9에서와 같이 다항식으로 근사화가 가능하다. 이때 다항식의 차수가 높을수록 근사 정확도가 증가하는 반면, 다수의 파라미터를 최적화하여야 하므로 수렴 능력이 저하될 수 있다. 따라서, 근사 정확도가 높으면서도 파라미터의 수가 적은 전체 두께 모델을 설정하는 것이 바람직하다. 이를 위하여, 두가지 물질 조합을 갖는 연속두께 팬텀을 이용하여 하나의 파라미터를 갖는 전체 두께 모델을 설정할 수 있다.
도 6은 두가지 물질 조합을 갖는 연속두께 팬텀을 이용하여 하나의 파라미터를 갖는 전체 두께 모델을 설정하는 방법의 개념을 보여주는 도면이다.
도 6을 참조하면, 피사체 예를 들면 유방의 주요 구성 물질은 지방조직(adipose tissue)와 선상조직(glandular tissue)이며, 각 조직의 밀도는 개인마다 다를 수 있다. 이와 같은 밀도 변화를 반영할 수 있는 전체 두께 모델을 구하기 위해서는, 먼저 어느 하나의 물질, 예를 들면 선상조직의 평균 기준밀도를 갖는 연속두께 팬텀(600)을 촬영하여 각 에너지대역별 팬텀 영상 I(n)과 이에 대응하는 두께 영상 T(m)을 획득한다. 다음, 팬텀 영상 I(n)이 선상조직의 밀도(d)가 특정값 즉, 평균 기준밀도(r)에서 얻은 것이라고 가정할 때 이에 대응하는 두께 영상 T'(m)을 획득한다. 이때, T'(m)은 T(m)/r로 계산된다. 만약, 팬텀 영상 I(n)과 두께 영상 T'(m) 사이의 근사 다항식 계수를 P'이라고 하면, 이는 선상조직의 밀도(d)가 r인 부분에서 실제로 촬영한 팬텀 영상과 이에 대응하는 두께 영상을 이용하여 얻어낸 다항식 계수와 거의 일치할 것이다. 따라서, 선상조직의 모든 밀도에 대하여 팬텀 영상을 촬영할 필요없이 평균 기준밀도의 팬텀 영상 한장으로 선상조직의 임의의 밀도에서의 다항식 계수를 계산할 수 있다.
결과적으로, 밀도(d)와 같은 하나의 파라미터를 이용하여 피사체의 전체 두께를 표현할 수 있고, 피사체의 전체 두께는 수학식 8 및 수학식 9의 다항식 모델에 적용할 수 있다. 즉, 피사체의 전체 두께를 f(x)라고 했을 때, f(x)=수학식 8, 또는 f(x)=수학식 9로 근사화시킬 수 있다.
도 7은 본 발명의 실시예에 따른 의료영상시스템의 구성을 보여주는 블럭도이다. 도 7에 도시된 의료영상시스템(700)은 방사선 영상 촬영부(710), 영상처리부(730), 디스플레이부(750), 저장부(770) 및 통신부(790)를 포함할 수 있다. 여기서, 의료영상시스템(700)은 영상처리부(730)만으로 구현될 수 있다. 즉, 방사선 영상 촬영부(710), 디스플레이부(750), 저장부(770), 및 통신부(790)는 옵션으로 구비될 수 있다. 한편, 영상처리부(730)는 적어도 하나의 프로세서로 구현될 수 있다.
도 7을 참조하면, 방사선 영상 촬영부(710)는 적어도 두가지 이상의 서로 다른 에너지 대역의 방사선을 피사체에 조사하여 피사체에 대한 다중 에너지 방사선 영상을 촬영한다. 피사체에서 동일한 조직에 대하여 조사되는 방사선 예를 들면, 엑스-레이의 에너지 대역이 다르면, 해당 조직에 흡수되는 정도가 다르다. 이러한 특성을 이용하여 두가지 이상의 에너지 대역의 엑스-레이를 각각 해당 조직에 조사하여 각 에너지대역별 흡수특성이 반영된 다중 에너지 방사선 영상을 획득한다.
한편, 방사선 영상 촬영부(710)가 의료영상시스템(700)에 포함되지 않을 경우에는, 외부로부터 제공되는 다중 에너지 방사선 영상이 영상처리부(730)로 입력된다.
영상처리부(730)는 방사선 영상 촬영부(710) 또는 외부로부터 제공되는 다중 에너지 방사선 영상에서 피사체의 전체 두께와 다중 에너지 방사선 영상의 국부영역에서 산출된 두께를 비교하여 비정상 조직을 추출하고, 추출된 비정상 조직을 포함하는 진단영상을 생성한다. 즉, 영상처리부(730)에서 생성된 진단영상에는 정상조직이 제거되어 있다. 한편, 영상처리부(730)는 다중 에너지 방사선 영상에 대하여 잡음 저감 처리를 수행하거나, 진단영상에 대하여 잡음 저감 처리 및/또는 대조도 향상 처리를 수행할 수 있다. 한편, 영상처리부(730)는 영상 판독 기능을 더 구비하여, 진단영상으로부터 필요로 하는 진단정보를 생성할 수 있다.
디스플레이부(750)는 모니터 등으로 구현되며, 영상처리부(730)에서 생성된 방사선 진단영상을 디스플레이하거나, 진단정보를 진단영상과 함께 디스플레이할 수 있다.
저장부(770)는 메모리 등으로 구현되며, 영상처리부(730)에서 생성된 진단영상을 저장하거나, 영상처리부(730)에서 얻어진 진단정보를 진단영상에 대응시켜 저장할 수 있다.
통신부(790)는 원격지에 위치한 다른 의료영상시스템 혹은 병원의 의사와 같은 전문가에게 영상처리부(730)에서 생성된 방사선 진단영상 혹은 진단정보와 결합된 진단영상을 유선 혹은 무선으로 전송하거나, 외부로부터 제공되는 다중 에너지 방사선 영상을 수신하여 영상처리부(730)로 입력할 수 있다. 특히, 다중 에너지 방사선 영상을 송신해 준 다른 의료영상시스템 혹은 전문가에게 방사선 진단영상 혹은 진단정보와 결합된 진단영상을 유선 혹은 무선으로 전송할 수 있다.
한편, 저장부(770)와 통신부(790)는 영상 판독 및 검색 기능을 더 포함시켜 PACS(Picture Archiving Communication System)로 일체화될 수 있다.
또는, 영상처리부(730), 저장부(770) 및 통신부(790)를 PACS로 일체화시킬 수 있다.
한편, 도 7에 도시된 의료영상시스템은 방사선 예를 들면 엑스-레이를 사용하는 모든 영상진단시스템일 수 있다. 예를 들면, 인체 중 뼈 이외의 연조직만으로 구성된 가슴 근육의 병변을 감별하는데 이용되는 마모그래피(mammography) 영상진단시스템일 수 있다.
도 8은 본 발명의 실시예에 따른 영상 처리 장치의 구성을 보여주는 블럭도이다. 영상 처리 장치(730)는 전체 두께 추정부(810)와 진단영상 생성부(830)를 포함할 수 있다. 여기서, 전체 두께 추정부(810)와 진단영상 생성부(830)는 적어도 하나의 프로세서로 구현될 수 있다.
도 8을 참조하면, 전체 두께 추정부(810)는 입력되는 다중 에너지 방사선 영상으로부터 피사체의 전체 두께를 추정한다. 이때, 두가지 물질 조합을 갖는 연속 두께 팬텀을 이용하여 전체 두께가 얻어지며, 두가지 물질 중 어느 하나의 물질에 대한 밀도를 파라미터로 하는 전체 두께 모델을 다중 에너지 방사선 영상에 적용할 수 있다.
진단영상 생성부(830)는 전체 두께 추정부(810)에서 얻어지는 피사체의 전체 두께와 다중 에너지 방사선 영상의 국부영역에서 산출된 두께를 비교하고, 비교결과에 따라서 관심 영역 예를 들면, 비정상 조직을 추출하고, 추출된 비정상 조직을 포함하는 진단영상을 생성한다.
본 발명은 또한 컴퓨터로 읽을 수 있는 기록매체에 컴퓨터가 읽을 수 있는 코드로서 구현하는 것이 가능하다. 컴퓨터가 읽을 수 있는 기록매체는 컴퓨터 시스템에 의하여 읽혀질 수 있는 데이터가 저장되는 모든 종류의 기록장치를 포함한다. 컴퓨터가 읽을 수 있는 기록 매체의 예로는 ROM, RAM, CD-ROM, 자기 테이프, 플로피 디스크, 광 데이터 저장 장치 등이 있으며, 또한 캐리어 웨이브 (예를 들어 인터넷을 통한 전송)의 형태로 구현되는 것도 포함한다. 또한 컴퓨터가 읽을 수 있는 기록 매체는 네트워크로 연결된 컴퓨터 시스템에 분산되어, 분산 방식으로 컴퓨터가 읽을 수 있는 코드가 저장되고 실행될 수 있다.
상기한 설명에서 많은 사항이 구체적으로 기재되어 있으나, 그들은 발명의 범위를 한정하는 것이라기보다, 실시예의 예시로서 해석되어야 한다. 예들 들어, 본 발명이 속하는 기술 분야에서 통상의 지식을 가진 자라면, 본 발명의 실시예에서 다양하게 변형할 수 있음을 알 수 있을 것이다. 따라서, 본 발명의 범위는 설명된 실시예에 의하여 정하여 질 것이 아니고 특허 청구범위에 기재된 기술적 사상에 의해 정하여져야 한다.
600 : 연속두께 팬텀 700 : 의료영상시스템
710 : 방사선영상 촬영부 730 : 영상처리부
750 : 디스플레이부 770 : 저장부
790 : 통신부 810 : 전체 두께 추정부
830 : 진단영상 생성부

Claims (16)

  1. 영상 처리 장치의 영상 처리 방법에 있어서,
    적어도 두가지 물질로 이루어지는 소정 피사체로부터 얻어지는 적어도 두가지 에너지 대역의 방사선 영상으로부터 상기 피사체의 전체 두께를 추정하는 단계; 및
    상기 추정된 전체 두께와 상기 방사선 영상에서 국부 영역의 두께를 비교하고, 비교결과에 따라서 관심 영역을 추출하여 진단영상을 생성하는 단계를 포함하는 영상 처리 방법.
  2. 제1 항에 있어서, 상기 관심 영역은 비정상 조직인 영상 처리 방법.
  3. 제2 항에 있어서, 상기 두가지 물질은 정상 조직의 구성물질인 영상 처리 방법.
  4. 제1 항에 있어서, 상기 전체 두께 추정단계는 두가지 물질 조합을 갖는 연속두께 팬텀을 이용하여 얻어는 전체 두께 모델을 적용하여 수행되는 영상 처리 방법.
  5. 제4 항에 있어서, 상기 전체 두께 모델은 상기 두가지 물질 중 어느 하나의 물질의 밀도를 파라미터로 포함하는 영상 처리 방법.
  6. 제1 항에 있어서, 상기 전체 두께 추정단계는
    상기 두가지 물질의 밀도를 추정하는 단계;
    상기 추정된 두가지 물질의 밀도와 각 물질의 감쇠 베이시스를 이용하여 상기 방사선 영상에서 상기 두가지 물질을 분리하는 단계; 및
    상기 분리된 두가지 물질의 두께를 더하여 상기 피사체의 전체 두께를 추정하는 단계를 포함하는 영상 처리 방법.
  7. 적어도 두가지 물질로 이루어지는 소정 피사체로부터 얻어지는 적어도 두가지 에너지 대역의 방사선 영상 데이터를 획득하는 방사선 영상 촬영부; 및
    상기 방사선 영상 데이터로부터 상기 피사체의 전체 두께를 추정하고, 상기 추정된 전체 두께와 상기 방사선 영상 데이터에서 국부 영역의 두께를 비교하고, 비교결과에 따라서 관심 영역을 추출하여 진단영상을 생성하는 적어도 하나의 프로세서를 포함하는 영상 처리 장치.
  8. 제7 항에 있어서, 상기 관심 영역은 비정상 조직인 영상 처리 장치.
  9. 제8 항에 있어서, 상기 두가지 물질은 정상 조직의 구성물질인 영상 처리 장치.
  10. 제7 항에 있어서, 상기 적어도 하나의 프로세서는 두가지 물질 조합을 갖는 연속두께 팬텀을 이용하여 얻어는 전체 두께 모델을 적용하여 수행되는, 영상 처리 장치.
  11. 제10 항에 있어서, 상기 전체 두께 모델은 상기 두가지 물질 중 어느 하나의 물질의 밀도를 파라미터로 포함하는 영상 처리 장치.
  12. 제7 항에 있어서, 상기 적어도 하나의 프로세서는 상기 두가지 물질의 밀도를 추정하고, 상기 추정된 두가지 물질의 밀도와 각 물질의 감쇠 베이시스를 이용하여 상기 방사선 영상 데이터에서 상기 두가지 물질을 분리하고, 상기 분리된 두가지 물질의 두께를 더하여 상기 피사체의 전체 두께를 추정하는, 영상 처리 장치.
  13. 제7 항 내지 제12 항 중 어느 한 항에 기재된 영상 처리 장치를 채용하는 의료영상시스템.
  14. 제13 항에 있어서, 상기 방사선 영상 촬영부는 상기 소정 피사체에 적어도 두가지 에너지 대역의 방사선을 조사하여 상기 방사선 영상 데이터를 획득하는, 의료영상시스템.
  15. 제13 항에 있어서, 상기 생성된 진단영상을 저장하거나, 상기 생성된 진단영상으로부터 얻어지는 진단정보를 상기 진단영상에 대응시켜 저장하는 저장부를 더 포함하는 의료영상시스템.
  16. 제13 항에 있어서, 상기 생성된 진단영상을 전송하거나, 상기 생성된 진단영상으로부터 얻어지는 진단정보를 상기 진단영상에 대응시켜 전송하는 통신부를 더 포함하는 의료영상시스템.
KR1020100073697A 2010-07-29 2010-07-29 영상 처리 방법 및 장치와 이를 채용한 의료영상시스템 KR101689866B1 (ko)

Priority Applications (6)

Application Number Priority Date Filing Date Title
KR1020100073697A KR101689866B1 (ko) 2010-07-29 2010-07-29 영상 처리 방법 및 장치와 이를 채용한 의료영상시스템
EP11812809.9A EP2598031A4 (en) 2010-07-29 2011-07-29 Method and apparatus for processing image, and medical image system employing the apparatus
CN2011800372808A CN103037768A (zh) 2010-07-29 2011-07-29 用于处理图像的方法和设备以及采用该设备的医学图像***
PCT/KR2011/005631 WO2012015280A2 (en) 2010-07-29 2011-07-29 Method and apparatus for processing image, and medical image system employing the apparatus
US13/813,007 US9008351B2 (en) 2010-07-29 2011-07-29 Method and apparatus for processing image, and medical imaging system employing the apparatus
JP2013521719A JP5536284B2 (ja) 2010-07-29 2011-07-29 映像処理方法及びその装置、並びにそれを用いた医療映像システム

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020100073697A KR101689866B1 (ko) 2010-07-29 2010-07-29 영상 처리 방법 및 장치와 이를 채용한 의료영상시스템

Publications (2)

Publication Number Publication Date
KR20120011679A KR20120011679A (ko) 2012-02-08
KR101689866B1 true KR101689866B1 (ko) 2016-12-27

Family

ID=45530637

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020100073697A KR101689866B1 (ko) 2010-07-29 2010-07-29 영상 처리 방법 및 장치와 이를 채용한 의료영상시스템

Country Status (6)

Country Link
US (1) US9008351B2 (ko)
EP (1) EP2598031A4 (ko)
JP (1) JP5536284B2 (ko)
KR (1) KR101689866B1 (ko)
CN (1) CN103037768A (ko)
WO (1) WO2012015280A2 (ko)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102012215515A1 (de) * 2012-08-31 2014-03-06 Klinikum Der Universität München Verfahren zur Detektion von Schäden an Silikonimplantaten und Computertomographiegerät
KR101479212B1 (ko) 2012-09-05 2015-01-06 삼성전자 주식회사 엑스선 영상 장치 및 엑스선 영상 생성 방법
KR102025753B1 (ko) 2012-09-21 2019-09-26 삼성전자주식회사 객체 정보 추정 장치 및 동작 방법
WO2014097026A1 (en) * 2012-12-21 2014-06-26 Koninklijke Philips N.V. Breast thickness measurement in mammography
CN103473536A (zh) * 2013-09-13 2013-12-25 天津大学 最大粘聚超像素网格化方法与超像素级图像处理方法
KR102393294B1 (ko) * 2014-09-26 2022-05-03 삼성전자주식회사 의료 영상 장치 및 의료 영상 장치의 제어 방법
JP6745756B2 (ja) * 2017-05-18 2020-08-26 富士フイルム株式会社 放射線画像撮影システム、放射線画像撮影方法、放射線画像撮影プログラム、及び導出装置
JP6745755B2 (ja) * 2017-05-18 2020-08-26 富士フイルム株式会社 放射線画像撮影システム、放射線画像撮影方法、放射線画像撮影プログラム、及び体厚推定装置
KR20210126647A (ko) * 2019-05-27 2021-10-20 다이아트렌드 코포레이션 데이터 처리 장치 및 데이터 처리 방법, 및, 그 장치 또는 방법을 탑재 또는 실시하는 x선 시스템
JP7373323B2 (ja) * 2019-09-02 2023-11-02 キヤノン株式会社 画像処理装置、放射線撮像システム、画像処理方法及びプログラム
KR102336229B1 (ko) * 2020-01-20 2021-12-06 연세대학교 원주산학협력단 3차원 레이저 스캐너를 이용한 단일 에너지 물질 분리 시스템, 장치 및 방법
KR102357455B1 (ko) * 2020-03-05 2022-02-04 한국과학기술원 X선 필터, x선 영상 촬영 시스템 및 x선 영상 촬영 방법

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050187478A1 (en) 2001-07-16 2005-08-25 Art, Advanced Research Technologies Inc. Multi-wavelength imaging of highly turbid media
JP2006519625A (ja) 2002-11-27 2006-08-31 ホロジック， インコーポレイテッド 組織露光制御、トモシンセシス及び動的視野処理を具備するフルフィールドマンモグラフィー
JP2007111525A (ja) 2005-10-17 2007-05-10 Siemens Ag Ｘ線コンピュータ断層撮影におけるコンピュータ断層撮影画像作成方法

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS63247870A (ja) 1987-04-02 1988-10-14 Konica Corp Ｘ線画像処理方法
US5247559A (en) 1991-10-04 1993-09-21 Matsushita Electric Industrial Co., Ltd. Substance quantitative analysis method
JP2973643B2 (ja) * 1991-10-04 1999-11-08 松下電器産業株式会社 物質の定量測定法
JPH05161631A (ja) 1991-12-18 1993-06-29 Toshiba Corp 放射線診断装置
CA2350648A1 (en) 1998-11-20 2000-06-02 Direct Radiography Corporation Interactive digital radiographic system
FR2801776B1 (fr) 1999-12-03 2002-04-26 Commissariat Energie Atomique Procede d'utilisation d'un systeme d'osteodensitometrie, par rayonnement x bi-energie, a faisceau conique
US6674835B2 (en) * 2001-10-12 2004-01-06 General Electric Co. Methods and apparatus for estimating a material composition of an imaged object
US6950492B2 (en) * 2003-06-25 2005-09-27 Besson Guy M Dynamic multi-spectral X-ray projection imaging
WO2006055830A2 (en) 2004-11-15 2006-05-26 Hologic, Inc. Matching geometry generation and display of mammograms and tomosynthesis images
WO2006138717A2 (en) * 2005-06-16 2006-12-28 California Pacific Medical Center Research Institute Methods and apparatus for determining proportions of body materials
US20080167552A1 (en) 2007-01-04 2008-07-10 General Electric Company System and method of generating an image of a contrast agent injected into an imaged subject
US8031835B2 (en) * 2008-08-07 2011-10-04 Xcision Medical Systems Llc Method and system for translational digital tomosynthesis mammography

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050187478A1 (en) 2001-07-16 2005-08-25 Art, Advanced Research Technologies Inc. Multi-wavelength imaging of highly turbid media
JP2006519625A (ja) 2002-11-27 2006-08-31 ホロジック， インコーポレイテッド 組織露光制御、トモシンセシス及び動的視野処理を具備するフルフィールドマンモグラフィー
JP2007111525A (ja) 2005-10-17 2007-05-10 Siemens Ag Ｘ線コンピュータ断層撮影におけるコンピュータ断層撮影画像作成方法

Also Published As

Publication number Publication date
EP2598031A4 (en) 2018-03-21
WO2012015280A9 (en) 2012-07-05
CN103037768A (zh) 2013-04-10
JP2013541349A (ja) 2013-11-14
KR20120011679A (ko) 2012-02-08
WO2012015280A3 (en) 2012-05-31
US20130129180A1 (en) 2013-05-23
EP2598031A2 (en) 2013-06-05
WO2012015280A2 (en) 2012-02-02
JP5536284B2 (ja) 2014-07-02
US9008351B2 (en) 2015-04-14

Similar Documents

Publication Publication Date Title
KR101689866B1 (ko) 영상 처리 방법 및 장치와 이를 채용한 의료영상시스템
EP2598032B1 (en) Method and apparatus for processing image and medical image system employing the apparatus
KR101747300B1 (ko) 영상을 처리하는 방법, 이를 수행하는 영상처리장치 및 의료영상시스템
US8363915B2 (en) Device, method and computer readable recording medium containing program for separating image components
US8611633B2 (en) Method for processing image, image processing apparatus and medical image system for performing the same
US10064591B2 (en) System, method and computer readable medium for preview of low-dose x-ray projection and tomographic images
Stolzmann et al. Monoenergetic computed tomography reconstructions reduce beam hardening artifacts from dental restorations
JP2001238870A (ja) 器官の厚みを補償する方法および装置
KR101820844B1 (ko) 진단영상 생성장치, 의료영상시스템 및 이미지를 처리하는 방법
Wang et al. Low‐dose preview for patient‐specific, task‐specific technique selection in cone‐beam CT
US20180253838A1 (en) Systems and methods for medical imaging of patients with medical implants for use in revision surgery planning
JP2015089428A (ja) 放射線画像処理装置および方法並びにプログラム
US20150172681A1 (en) Medical image compression system and method using visually lossless compression
KR101853814B1 (ko) 방사선 영상 촬영장치, 의료영상시스템 및 방사선 영상을 촬영하는 방법
KR101364904B1 (ko) 4차원 영상합성법의 재구성 장치 및 방법
Boos et al. Metal artifact reduction (MAR) based on two-compartment physical modeling: evaluation in patients with hip implants
CN111369635A (zh) 预测截断图像的***和方法、准备数据的方法及其介质
Putera et al. Simulation of Three-Dimensional Images and Estimation of Lung Volumes from Two-Dimensional MRI and CT Images
Ott Iterative reconstruction in CT: using mathematical model observers to determine low dose images' trustworthiness
KR20180014932A (ko) 의료방사선 장치의 촬영조건 추천 방법 및 장치

Legal Events

Date Code Title Description
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20191128

Year of fee payment: 4