KR100438903B1 - 초음파 영상으로부터 대상 장기의 지방량을 정량적으로 측정하는 장비 - Google Patents

초음파 영상으로부터 대상 장기의 지방량을 정량적으로 측정하는 장비 Download PDF

Info

Publication number
KR100438903B1
KR100438903B1 KR10-2002-0005642A KR20020005642A KR100438903B1 KR 100438903 B1 KR100438903 B1 KR 100438903B1 KR 20020005642 A KR20020005642 A KR 20020005642A KR 100438903 B1 KR100438903 B1 KR 100438903B1
Authority
KR
South Korea
Prior art keywords
brightness value
fat
target
target area
value
Prior art date
Application number
KR10-2002-0005642A
Other languages
English (en)
Other versions
KR20030065155A (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 KR10-2002-0005642A priority Critical patent/KR100438903B1/ko
Priority to US10/205,104 priority patent/US6656121B2/en
Publication of KR20030065155A publication Critical patent/KR20030065155A/ko
Application granted granted Critical
Publication of KR100438903B1 publication Critical patent/KR100438903B1/ko

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/48Other medical applications
    • A61B5/4869Determining body composition
    • A61B5/4872Body fat
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Surgery (AREA)
  • Physics & Mathematics (AREA)
  • Veterinary Medicine (AREA)
  • Pathology (AREA)
  • Public Health (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

초음파 영상으로부터 대상 장기의 지방량을 정량적으로 측정하는 장비를 제공한다. 본 발명의 일 관점에 의한 장비는, 복부 단면 초음파 영상을 얻는 초음파 영상 촬영부, 초음파 영상의 데이터를 저장하는 기억 저장부, 및 연산부를 포함하여 구성되되, 연산부는, 초음파 영상 데이터를 이용하여 간을 포함하는 영역으로 설정된 대상 영역의 화소 분포 패턴을 얻어 각각의 화소의 명도값들에 대한 화소의 수의 히스토그램 분포를 얻은 후 히스토그램 분포의 변곡점에 위치하는 명도값을 대상 영역의 대표 명도값으로 추출하고, 영상 데이터에서 복부 피하 지방 영역 및 신장을 포함하는 영역으로 설정된 비교 영역을 가로지는 경로 상의 영상의 화소의 수의 명도값들의 분포 곡선을 얻고 명도값들의 평균 또는 명도값들의 분포 곡선의 최대점들에서의 명도값들의 평균 또는 최소점들에서의 명도값들의 평균을 비교 영역의 대표 명도값으로 추출하고, 대상 영역의 대표 명도값을 비교 영역의 대표 명도값에 대한 상대적인 수치로 정량화하여 대상 영역의 정량화된 대표 명도값을 얻고, 대상 영역의 정량화된 대표 명도값을 대상 영역의 지방량으로 대응시키는 연산을 수행한다.

Description

초음파 영상으로부터 대상 장기의 지방량을 정량적으로 측정하는 장비{Apparatus of quantitatively measuring fat content in target organ from ultrasound visual image}
본 발명은 인체에 대한 초음파 단면 영상 장비에 관한 것으로, 특히, 초음파 단면 영상의 화소들의 명도 분포(gray level distribution)로부터 대상 장기의 대표 명도를 정량화하여 대상 장기의 지방량을 측정하는 장비에 관한 것이다.
인체의 장기 조직의 지방화는 여러 가지 합병증을 유발할 수 있고, 지방화 그 자체로도 장기의 기능 장애의 위험성을 내포하고 있다. 예를 들어, 지방간은 인체 내의 간 세포가 비정상적으로 과다하게 지방을 포함하고 있는 상태로, 그 자체는 통증을 유발하지 않고 별다른 자각 증상도 없는 특징이 있다. 미미한 정도의 지방간은 건강한 사람에게도 발견될 수 있으나, 지방화 정도가 심할 경우 합병증이나 간경변 등으로 발전하는 경우가 빈번하다. 최근에는 현대인의 식생활 및 비만 문제와 관련하여 이러한 장기의 지방화 정도를 측정하는 데 관심이 집중되고 있다. 이에 따라, 이러한 장기의 지방화 정도를 간편하고 저렴하게 측정하는 방법의 개발에 관심이 집중되고 있다.
지방간은 간 실질 내의 지방량과 매우 밀접한 관련을 가진다고 알려져 있으며, 정량적인 지방량 측정을 위한 여러 가지 방법이 제안되고 있거나 사용되고 있다. 그럼에도 불구하고, 간편하고 저렴하게 장기의 지방량 정도를 측정하는 방법의 개발을 계속적으로 요구되고 있다.
예를 들어, 컴퓨터 단층 촬영(CT:Computed tomography)은 이러한 인체 장기의 단면 영상을 얻는 데 많이 이용되고 있다. CT는 촬영한 영상 해상도가 높고 영상 내에서 지방 및 그 외의 성분의 비율을 촬영 기술 및 촬영 조건에 비교적 영향 받지 않고 정확히 분리해 낼 수 있는 방법으로 최근 많은 연구가 집중되고 있다. 그럼에도 불구하고, CT 촬영은 매우 많은 비용을 요구하고 있으며, 방사선 피폭을 감수하여야 하므로 안정성 등의 문제가 제기될 수 있다.
또한, 이러한 장기의 지방량을 측정하기 위해서 간 생체 검사가 이루어지고 있다. 간 생체 검사는 장기, 예를 들어, 간 조직을 떼어내어 조직 및 성분을 분석하는 방법으로, 모든 간질환 진단 및 검사에서 가장 근본적이고 정확한 결과를 제공한다. 그럼에도 불구하고, 절차가 복잡하고 침습적(invasive)이다. 또한, 질병의 징후가 나타나지 않는 경우에는 간 생체 검사가 행해지지 않으므로 예방 조치를 위해서는 적용되기 어렵다.
초음파 촬영법은 장기의 단면 영상, 예를 들어, 복부 단면을 비롯한 생체 각 부위에 대한 영상을 얻어 이에 대한 진단에 광범위하게 이용되고 있다. 이러한 초음파 촬영법은 상대적으로 저렴하고 촬영이 간편하고 근본적으로 방사선을 이용하지 않아 안전한 방법이면서 비침습적이다. 그럼에도 불구하고, 초음파 영상의 품질에 대한 촬영 조건 및 촬영 기술의 영향이 절대적이며, 초음파가 통과하는 부분별 반사 특성 및 해상도가 일정하지 않는 단점이 있다. 이에 따라, 촬영된 영상을 이용하여 각종 진단이 행해질 때 진단자의 자의적 해석이 개입될 여지가 많으며, 이에 따라, 진단의 신뢰성이 완벽하게 보장되지 못하고 있다.
현재, 초음파 촬영법에 의한 영상을 이용하여 지방간 지수를 임상적으로 판단하는 데 있어서 정상, 경도, 중등도, 고도의 대략 4단계로 구별하고 있는 데, 진단 경험이 있는 전문의들 간의 오차가 1 단계 정도 틀리는 경우가 다분하며 이러한 오차가 발생할 확률은 최대 20% 정도에 이르고 있다.
따라서, 간편하고 저렴한 초음파 촬영 영상을 이용하면서도 촬영 조건 등에 의존하지 않고 정량화된 지방량 값이 제공될 수 있는 방법의 개발이 크게 요구되고 있다. 초음파 촬영 영상만을 이용하여 지방량을 측정하기 위한 많은 시도가 보고되고 있으나, 초음파 촬영 방법이 가지는 근원적인 문제, 예를 들어, 명도 분포(gray level distribution) 및 해상도의 불규칙성에 의한 영상 특성 왜곡에 의해서 지방량을 정량적으로 추출하기 어려운 점 및 신뢰성을 확보하기 어려운 점 등을 극복하기는 현재로서는 어렵다.
본 발명이 이루고자 하는 기술적 과제는, 초음파 촬영에 의해서 얻어지는 영상으로부터 대상 장기의 명도 분포 특성을 정량화하여 대표 명도값을 추출하고 이러한 대표 명도값으로부터 지방량을 정량적으로 측정하는 장비를 제공하는 데 있다.
도 1은 본 발명의 실시예에 의한 초음파 영상으로부터 대상 장기의 지방량을 정량적으로 측정하는 방법을 개략적으로 나타내는 흐름도이다.
도 2는 본 발명의 실시예에 의한 정량화된 대표 명도값을 얻는 방법을 개략적으로 나타내는 흐름도이다.
도 3은 본 발명의 실시예에 의한 정량화된 대표 명도값으로부터 실제 지방량을 결정하는 방법을 개략적으로 나타내는 흐름도이다.
도 4는 본 발명의 실시예에 의한 초음파 영상으로부터 대상 장기의 지방량을 정량적으로 측정하는 방법에 이용되는 장비를 설명하기 위해서 개략적으로 도시한 도면이다.
도 5는 초음파 영상 촬영으로 얻어진 영상의 일례를 모식적으로 도시한 도면이다.
도 6은 복부 피하 지방 영역인 제1비교 영역을 모식적으로 도시한 도면이다.
도 7은 도 6의 종단 경로 상의 명도 분포를 곡선으로 개략적으로 도시한 도면이다.
도 8은 도 6의 종단 경로 상의 명도값에 대한 화소 수의 분포를 개략적으로 도시한 히스토그램(histogram) 분포도이다.
도 9는 신장 영역인 제2비교 영역을 모식적으로 도시한 도면이다.
도 10은 도 9의 종단 경로 상의 명도 분포를 곡선으로 개략적으로 도시한 도면이다.
도 11은 도 9의 종단 경로 상의 명도값에 대한 화소 수의 분포를 개략적으로 도시한 히스토그램 분포도이다.
도 12는 간 실질인 대상 영역을 모식적으로 도시한 도면이다.
도 13은 도 12의 설정 영역 내의 명도값에 대한 화소 수의 분포를 개략적으로 도시한 히스토그램 분포도이다.
도 14는 도 13의 히스토그램 분포로부터 필터링되어 분리된 히스토그램 분포 곡선 그래프(graph)이다.
도 15는 지방간 지수와 정량화된 대표 명도값의 상관 관계를 나타내는 그래프이다.
<도면의 주요 부호에 대한 간략한 설명>
510: 간 실질인 대상 영역,
530: 복부 피하 지방 영역인 제1비교 영역,
550: 신장인 제2비교 영역.
상기의 기술적 과제들을 달성하기 위한 본 발명의 일 관점은, 복부 단면 초음파 영상을 얻는 초음파 영상 촬영부, 상기 초음파 영상의 데이터를 저장하는 기억 저장부, 및 연산부를 포함하여 구성되는 초음파 단면 영상으로부터 대상 장기의 지방량을 정량적으로 측정하는 장비를 제시한다.상기 장비의 연산부는 상기 초음파 영상 데이터를 이용하여 상기 영상 데이터에서 간을 포함하는 영역으로 설정된 대상 영역의 화소 분포 패턴을 얻어 각각의 화소의 명도값들에 대한 상기 화소의 수의 히스토그램 분포를 얻은 후 상기 히스토그램 분포의 변곡점에 위치하는 명도값을 상기 대상 영역의 대표 명도값으로 추출하고, 상기 영상 데이터에서 복부 피하 지방 영역 및 신장을 포함하는 영역으로 설정된 비교 영역을 가로지는 경로 상의 상기 영상의 화소의 수의 명도값들의 분포 곡선을 얻고 상기 명도값들의 평균 또는 상기 명도값들의 분포 곡선의 최대점들에서의 명도값들의 평균 또는 최소점들에서의 명도값들의 평균을 상기 비교 영역의 대표 명도값으로 추출하고, 상기 대상 영역의 대표 명도값을 상기 비교 영역의 대표 명도값에 대한 상대적인 수치로 정량화하여 상기 대상 영역의 정량화된 대표 명도값을 얻고, 상기 대상 영역의 정량화된 대표 명도값을 상기 대상 영역의 지방량으로 대응시키는 연산을 수행한다.상기 연산부는 상기 대상 영역의 히스토그램 분포가 다수의 변곡점들을 나타낼 때 상기 대상 영역의 대표 명도값을 상기 변곡점들 중 상기 히스토그램 분포에서 화소의 명도값의 최대값 보다 더 밝은 쪽의 첫 번째 변곡점에 위치하는 화소의 명도값으로 추출하는 연산을 수행할 수 있다.상기 연산부는 상기 대상 영역의 대표 명도값을 추출하기 이전에 상기 대상 영역의 히스토그램 분포를 가우시안 브로드닝(Gaussian broadening)하여 밝은 쪽의 최대값을 포함하는 분포 곡선을 필터링하여 상기 히스토그램 분포로부터 노이즈를 제거하는 연산을 수행할 수 있다.
본 발명에 따르면, 초음파 단면 영상으로부터 지방량 정도를 효과적으로 그리고 신뢰성있게 측정할 수 있어, 저렴한 비용으로 안전한 지방량 측정 방법을 제공할 수 있다.
이하, 첨부 도면을 참조하여 본 발명의 실시예를 상세히 설명한다. 그러나, 본 발명의 실시예들은 여러 가지 다른 형태로 변형될 수 있으며, 본 발명의 범위가 아래에서 상술하는 실시예들로 인해 한정되어지는 것으로 해석되어져서는 안된다. 본 발명의 실시예들은 당업계에서 평균적인 지식을 가진 자에게 본 발명을 보다 완전하게 설명하기 위해서 제공되어지는 것이다. 따라서, 도면에서의 요소의 형상 등은 보다 명확한 설명을 강조하기 위해서 과장되어진 것이며, 도면 상에서 동일한 부호로 표시된 요소는 동일한 요소를 의미한다.
본 발명의 실시예에서는, 인체의 관심 장기, 즉, 대상 장기 내에 지방이 많아지면서 나타나는 관심 장기의 초음파에 대한 반사 굴절 특성의 변화를 초음파 촬영 영상을 이용하여 지표화하고, 이러한 변화로부터 지방량 정도를 추산하는 방법을 제시한다.
초음파 촬영 영상, 예를 들어, 인체의 초음파 촬영 복부 단면 영상은 각각의 생체 조직의 초음파에 대한 반사 굴절 특성이 명도(gray level)로 변환되어 다시 2차원 영상으로 재구성된 것이다. 생체 조직의 주요 구성 성분으로는 지방과 수분을 들 수 있는 데, 초음파는 지방과 수분에 대해서 상이한 반사 특성을 나타내게 된다. 이에 따라, 특히 지방 영역은 촬영된 영상에서 밝게 나타나게 된다. 따라서, 촬영된 복부 단면 영상의 각각의 화소(pixel)들의 명도 또는 밝기 정도는 그 화소 영역에 대응하는 생체 조직 영역 내의 지방 함량에 비례하게 된다.
그러나, 화소의 이러한 명도 또는 밝기 정도 자체는 관심 장기의 지방량을 직접적으로 대변한다고 간주되기 어렵다. 이는 명도 자체는 초음파 복부 단면 촬영 조건 및 촬영 기술에 따라 상대적으로 달라질 수 있으므로, 이를 비교할 수 있는 지방량 지표로 이용하기는 어렵다. 즉, 동일한 대상에 대해 촬영한 초음파 복부 단면 영상의 동일한 위치에 대한 서로 다른 촬영 조건의 두 영상의 명도 또는 동일한 위치의 화소의 명도가 서로 다르게 나타날 수 있다. 따라서, 초음파 영상의 명도를 지방량으로 직접 사용하기는 어렵다.
이를 극복하기 위해서 본 발명의 실시예에서는 관심 장기 영역, 즉, 대상 영역의 대표 명도값을 추출하고 이를 정량화하여 지방량을 추산하는 방법을 제시한다. 즉, 정량화된 대표 명도값에 대한 초음파 단면 영상이 획득되는 조건, 예를 들어, 촬영 조건 등에 의한 영향을 배제하여 대표 명도값을 지방량에 대한 지표값으로 이용할 수 있도록 한다.
도 1을 참조하면, 본 발명의 실시예에 의한 대상 장기의 지방량을 측정하는 방법은, 먼저, 복부 단면 초음파 영상을 얻은(110) 후. 대상 영역, 예를 들어, 간 실질에 대한 정량화된 대표 명도값(representative gray level)을 얻은(120) 후, 간 실질에 대한 정량화된 대표 명도값으로부터 간 실질의 실제 지방량 정도를 계산하여 결정한다(130). 도 1은 본 발명의 실시예에 의한 초음파 영상으로부터 대상 장기의 지방량을 정량적으로 측정하는 방법을 개략적으로 나타내는 흐름도이다.
도 5에 도시된 바와 같은 복부 단면 초음파 영상은 인체의 복부에 대한 초음파 영상 촬영으로 얻어진다. 이러한 초음파 복부 단면 촬영 시에 촬영 각도와 밀착정도를 조정하여 영상 내의 명도 분포가 가능한 한 일정하도록 하는 것이 바람직하다. 또한, 이러한 영상 내에 촬영하고자 하는 장기, 예를 들어, 간 실질, 신장, 복부 피하 지방 영역 등의 단면이 빠짐없이 포함될 수 있도록 하는 것이 바람직하다.
이러한 초음파 영상 촬영은 도 4에 개략적으로 도시한 바와 같이 구성되는 대상 장기의 지방량을 정량적으로 측정하는 방법에 이용되는 장비의 초음파 영상 촬영부(440)에 의해서 이루어진다. 도 4는 본 발명의 실시예에 의한 초음파 영상으로부터 대상 장기의 지방량을 정량적으로 측정하는 방법에 이용되는 장비를 개략적으로 나타낸다. 이러한 초음파 영상 촬영부(440)는 일반적인 초음파 영상 촬영 장비로 이루어질 수 있다. 이러한 초음파 영상 촬영부(440)에 의해서 얻어지는 영상 데이터는 기억 저장부(430)에 데이터로 저장된다. 이러한 기억 저장부(430)에 저장된 데이터들은 입 출력부(410)의 제어에 의해서 연산부(420)에서 대표 명도값을 계산하고 결정하고 정량화하고 이로부터 지방량을 추산하는 데 이용된다.
이와 같이 얻어진 초음파 복부 단면 영상은 0에서 255 사이의 명도값들을 가지는 화소들로 이루어지며, 화소 수는 대략 500 ×700으로 설정될 수 있다. 이러한 데이터의 처리는 상기한 장비의 연산부(420)에 의해서 이루어질 수 있으며, 이러한 화상의 출력은 입 출력부(410)의 모니터(monitor) 등에 이루어질 수 있다.
초음파 복부 단면 영상으로부터 대상 영역, 예컨대, 간 실질의 대표 명도값을 추출하기 위해서, 먼저, 초음파 단면 영상 내에 대상 영역을 설정한다. 도 2는 본 발명의 실시예에 의한 정량화된 대표 명도값을 얻는 방법을 개략적으로 나타내는 흐름도이다. 구체적으로, 도 5에 도시된 바와 같은 초음파 복부 단면 영상 내의간 실질을 대상 영역(510)으로 설정한다. 그리고, 대상 영역(510)에 대해 추출되는 대표 명도값을 정량화하기 위한 비교 영역(530), 550)을 함께 설정한다(도 2의 111).
대상 영역(510)인 간 실질을 정량화하기 위한 비교 대상으로는 인체 내에 지방량이 최대인 장기를 선택할 수 있고, 또는, 인체 내의 지방량이 최소인 장기를 선택할 수 있다. 본 발명의 실시예에서는 지방량이 최대인 장기인 복부 피하 지방 영역(530)을 제1비교 영역(530)으로 설정하고, 지방량이 최소인 장기인 신장을 제2비교 영역(550)으로 설정하는 경우를 예로 들어 설명한다.
먼저, 설정된 영역이 간 실질이 아닌 비교 영역(530, 550)일 경우에, 간 실질의 대상 영역(510)으로부터 추출될 대상 영역의 대표 명도값을 정량화하는 기준으로 이용하기 위해서 비교 영역(530, 550)의 대표 명도값들을 추출한다(115). 이러한 비교 영역의 대표 명도값은 비교 영역(530, 550) 내의 화소들의 명도 분포들을 얻은 후, 이로부터 비교 영역의 대표 명도값을 결정한다.
구체적으로 설명하면, 복부 피하 지방 영역인 제1비교 영역(530)은 도 6에 도시된 바와 같이 모식적으로 도시될 수 있다. 이러한 복부 피하 지방 영역의 종단 경로를 설정한 후, 도 7에 도시된 바와 같이 이러한 종단 경로 상의 명도 분포를 곡선으로 도시할 수 있다. 도 7에 도시된 바와 같은 종단 경로 상의 명도값의 분포곡선의 최대점들(61, 62, 63)에 해당되는 명도값들을 추출한 후 평균하여 복부 피하 지방 영역인 제1비교 영역(530)의 대표 명도값으로 설정한다.
복부 피하 지방 영역은 지방량이 많은 영역이므로 이를 반영하도록 명도값이높은 부분을 추출하는 것이 대상 영역의 대표 명도값을 정량화하기 위한 비교 대상으로 적절하다. 지방 부분은 초음파 영상에서 밝게, 즉, 명도값이 높게 나타난다. 이러한 명도값의 추출을 위한 영역의 설정 및 연산은 도 4에 도시된 바와 같은 장비의 연산부(420), 예를 들어, 마이크로 프로세서(micro processor)를 통한 연산으로 이루어질 수 있다.
이러한 복부 피하 지방 영역인 제1비교 영역(530)에 대한 대표 명도값은 상기한 바와 같은 추출 방법에 의해서 대략 228.9로 결정될 수 있다.
한편, 이러한 복부 피하 지방 영역인 제1비교 영역(530)에 대한 대표 명도값은 도 6의 종단 경로 상의 명도값의 히스토그램 분포를 얻어 추출될 수 있다. 도 8은 종단 경로 상의 명도값에 대한 화소의 수의 히스토그램 분포를 나타낸다. 이러한 복부 피하 지방 영역은 지방량이 매우 높은 영역이므로 히스토그램 분포에서 명도값이 높은 영역을 설정하고 이들의 명도값을 평균하여 제1비교 영역(530)에 대한 대표 명도값으로 결정한다.
이와 같이 지방량이 많은 부분인 복부 피하 지방 영역을 제1비교 영역(530)으로 설정하여, 이러한 제1비교 영역(530)에 대한 대표 명도값을 얻고, 또한, 지방량이 적은 부분인 신장을 제2비교 영역(550)으로 설정한다. 도 9는 신장으로 설정된 제2비교 영역을 나타낸다. 이러한 제2비교 영역(도 5의 550)에 대한 종단 경로를 설정한 후, 도 10에 도시된 바와 같이 이러한 종단 경로 상의 명도 분포를 곡선으로 도시한다. 도 10에 도시된 바와 같은 종단 경로 상의 명도값의 분포 곡선의 최소점들(91, 92, 93, 94, 95, 96)에 해당되는 명도값들을 추출한 후 평균하여 신장인 제2비교 영역(550)의 대표 명도값으로 설정한다.
신장은 지방량이 적은 영역이므로 이를 반영하도록 명도값이 낮은 부분을 추출하는 것이 대상 영역의 대표 명도값을 정량화하기 위한 비교 대상으로 적절하다. 이러한 명도값의 추출을 위한 영역의 설정 및 연산은 도 4에 도시된 바와 같은 장비의 연산부(420), 예를 들어, 마이크로 프로세서를 통한 연산으로 이루어질 수 있다.
이러한 신장 영역인 제2비교 영역(550)에 대한 대표 명도값은 상기한 바와 같은 추출 방법에 의해서 대략 47.7로 결정될 수 있다.
한편, 이러한 신장인 제2비교 영역(550)에 대한 대표 명도값은 도 9의 종단 경로 상의 명도값의 히스토그램 분포를 얻어 추출될 수 있다. 도 11은 종단 경로 상의 명도값에 대한 화소의 수의 히스토그램 분포를 나타낸다. 이러한 신장 영역은 지방량이 매우 낮은 영역이므로 히스토그램 분포에서 명도값이 낮은 영역을 설정하고 이들의 명도값을 평균하여 제2비교 영역(550)에 대한 대표 명도값으로 결정할 수 있다.
이와 같이 대상 영역(510)에 대해 추출될 대표 명도값의 비교 대상으로 이용될 비교 영역(530, 550)의 대표 명도값을 구한 후, 대상 영역(510)의 대표 명도값을 구한다(도 2의 113).
구체적으로, 도 12에 도시된 바와 같은 대상 영역(510)인 간 실질은 비교적 지방의 분포가 전체에 걸쳐 균일한 것으로 알려져 있으므로, 대상 영역(510)의 일부를 설정 영역(121)으로 설정하여 이러한 설정 영역(121) 내의 화소들의 명도값을분석하여 대표 명도값을 추출할 수 있다. 도 12는 대상 영역(510)인 간 실질을 모식적으로 나타낸다.
이러한 대상 영역(510) 내의 설정 영역(121) 내의 화소들의 명도값을 분석하여 각각의 화소의 명도값에 대한 화소의 수의 히스토그램 분포를 얻은 후 이러한 히스토그램 분포로부터 대표 명도값을 추출한다(도 2의 113). 구체적으로 예를 들면, 도 13에 도시된 바와 같은 명도값에 대한 화소 수의 히스토그램 분포를 얻은 후, 노이즈(noise)를 제거한다. 이러한 노이즈를 제거하는 데에는 가우시안 브로드닝(Gaussian broadening)을 이용한다. 도 13은 대상 영역의 히스토그램 분포를 개략적으로 나타낸다.
도 13에 도시된 바와 같은 히스토그램 분포는 두 개의 분포 곡선들(201, 203)로 필터링(filtering)될 수 있다. 지방량에 관심을 맞추고 있으므로, 지방량에 연관되는 밝은 쪽의 명도값 영역의 분포 곡선(203)을 따로 분리한다. 도 14는 필터링되어 분리된 히스토그램 분포 곡선을 개략적으로 나타낸다. 이와 같이 분리된 히스토그램 분포 곡선(203)의 밝은 쪽의 첫 번째 변곡점(205)을 선택하여 이러한 지점의 명도값을 대표 명도값으로 결정한다. 이 경우, 간 실질인 대상 영역(510)에 대한 대표 명도값은 예를 들어, 198.6으로 추출될 수 있다.
도 2를 다시 참조하면, 상기한 바와 같이 대상 영역(510)에 대한 대표 명도값을 추출한 후, 이러한 대상 영역(510)의 대표 명도값을 정량화한다(117). 이러한 정량화는 대상 영역(510)이 포함된 초음파 복부 단면 영역의 다른 영역들에 대해서, 상기한 대상 영역(510)의 대표 명도값을 상대적인 수치로 정량화함으로써 이루어질 수 있다. 예를 들어, 상기한 바와 같은 지방량이 최대인 복부 피하 지방 영역인 제1비교 영역(530)의 대표 명도값과 지방량이 최소인 신장인 제2비교 영역(550)의 대표 명도값에 대해서 상대적인 수치로 대상 영역(510)의 대표 명도값을 정량화 또는 표준화한다.
예를 들어, 정량화된 대상 영역(510)의 대표 명도값을 C'으로 두고, 대상 영역(510)의 대표 명도값을 C로 두고, 제1비교 영역(530)의 대표 명도값을 A로 두고, 제2비교 영역(550)의 대표 명도값을 B로 두면, C'는 f(C;A,B)인 함수로 계산되어 정량화될 수 있고 f(C;A,B)=(C-B)/(A-B)로 설정하여 정량화된 대상 영역(510)의 대표 명도값인 C'를 계산하여 정량화된 수치로 얻을 수 있다. 이러한 함수 f(C;A,B)는 상기한 바와 같이 주어지는 함수뿐만 아니라 단조 증가(monotonically increase)하는 함수이면 어떤 함수라도 가능하다.
상기한 바와 같은 함수에 의해서 계산되어 정량화되는 대상 영역(510)의 대표 명도값은 (198.6-47.7)/(228.9-47.7)=0.833으로 얻어진다. 이와 같이 정량화된 대상 영역(510)의 대표 명도값은 0과 1사이의 값으로 비교 영역에 대한 상대적인 값으로 정량화되므로, 다른 초음파 단면 영상으로부터 얻어지는 수치와도 비교되는 것이 가능해진다.
도 1을 참조하면, 이와 같이 얻어진 정량화된 대상 영역(510)의 대표 명도값으로부터 실제 지방량을 얻는다(130). 도 3은 본 발명의 실시예에 의한 정량화된 대표 명도값으로부터 실제 지방량을 결정하는 방법을 개략적으로 나타내는 흐름도이다. 구체적으로, 대상 영역(510)의 정량화된 대표 명도값과 실제 대상 영역(510)의 지방량에 대한 상관 함수 관계를 먼저 얻는다. 이를 위해서, 먼저, 다수의 초음파 단면 영상 표본들에 대한 데이터를 수집한다(131). 이후에, 표본들 각각에 대한 대상 영역, 즉, 간 실질에 대한 지방간 지수들과 정량화된 대표 명도값을 얻는다(133). 이때, 정량화된 대표 명도값들은 앞서 설명한 바와 같이 얻을 수 있다.
이후에, 표본들의 정량화된 대표 명도값들과 지방간 지수들 간의 상관 관계를 얻는다(135). 예를 들어, 도 15에 도시된 바와 같이 지방간 지수에 대한 정량화된 대표 명도값, 즉, 정량화된 값이므로 지방화 지표로 간주할 수 있는 값의 분포를 얻는다. 도 15는 지방간 지수와 정량화된 대표 명도값의 상관 관계를 나타내는 그래프(graph)이다. 이때 이용되는 지방간 지수는 실제 임상 전문의의 판정 결과를 이용한다. 이후에, 지방간 지수와 정량화된 대표 명도값들 간의 비례 함수를 근사 지수법(least square)법으로 구한다. 이와 같은 비례 함수는 도 15의 비례 증가 직선(150)으로 나타낼 수 있다.
이와 같은 상관 관계 함수에 대상 영역(510)의 정량화된 대표 명도값을 대입하면 해당되는 지방간 지수가 얻어질 수 있다(137). 따라서, 본 발명의 실시예에서 얻어지는 정량화된 대표 명도값은 바로 지방화 지표로 추정 또는 간주될 수 있다.
이상, 본 발명을 구체적인 실시예를 통하여 상세히 설명하였으나, 본 발명은 이에 한정되지 않고, 본 발명의 기술적 사상 내에서 당 분야의 통상의 지식을 가진 자에 의해 그 변형이나 개량이 가능함이 명백하다.
상술한 본 발명에 따르면, 초음파 단면 영상을 이루는 화소들의 명도 분포를 선형적으로 재조정함으로써, 초음파 영상이 촬영 조건 및 촬영 기술에 의하여 받는 영향을 극복할 수 있고, 관심 장기의 대상 영역의 대표 명도값을 얻을 수 있다. 이러한 대표 명도값을 정량화함으로써 초음파 영상들 간의 비교가 가능해지고, 이에 따라, 관심 장기의 지방량 정도를 정량적으로 측정할 수 있다.

Claims (16)

  1. 삭제
  2. 삭제
  3. 삭제
  4. 삭제
  5. 삭제
  6. 삭제
  7. 삭제
  8. 삭제
  9. 삭제
  10. 삭제
  11. 삭제
  12. 삭제
  13. 삭제
  14. 복부 단면 초음파 영상을 얻는 초음파 영상 촬영부;
    상기 초음파 영상의 데이터를 저장하는 기억 저장부; 및
    상기 초음파 영상 데이터를 이용하여
    상기 영상 데이터에서 간을 포함하는 영역으로 설정된 대상 영역의 화소 분포 패턴을 얻어 각각의 화소의 명도값들에 대한 상기 화소의 수의 히스토그램 분포를 얻은 후 상기 히스토그램 분포의 변곡점에 위치하는 명도값을 상기 대상 영역의 대표 명도값으로 추출하고,
    상기 영상 데이터에서 복부 피하 지방 영역 및 신장을 포함하는 영역으로 설정된 비교 영역을 가로지는 경로 상의 상기 영상의 화소의 수의 명도값들의 분포 곡선을 얻고 상기 명도값들의 평균 또는 상기 명도값들의 분포 곡선의 최대점들에서의 명도값들의 평균 또는 최소점들에서의 명도값들의 평균을 상기 비교 영역의 대표 명도값으로 추출하고,
    상기 대상 영역의 대표 명도값을 상기 비교 영역의 대표 명도값에 대한 상대적인 수치로 정량화하여 상기 대상 영역의 정량화된 대표 명도값을 얻고,
    상기 대상 영역의 정량화된 대표 명도값을 상기 대상 영역의 지방량으로 대응시키는 연산을 수행하는 연산부를 포함하는 것을 특징으로 하는 초음파 단면 영상으로부터 대상 장기의 지방량을 정량적으로 측정하는 장비.
  15. 제14항에 있어서, 상기 연산부는
    상기 대상 영역의 히스토그램 분포가 다수의 변곡점들을 나타낼 때 상기 대상 영역의 대표 명도값을 상기 변곡점들 중 상기 히스토그램 분포에서 화소의 명도값의 최대값 보다 더 밝은 쪽의 첫 번째 변곡점에 위치하는 화소의 명도값으로 추출하는 연산을 수행하는 것을 특징으로 하는 초음파 영상으로부터 대상 장기의 지방량을 정량적으로 측정하는 장비.
  16. 제14항에 있어서, 상기 연산부는
    상기 대상 영역의 대표 명도값을 추출하기 이전에 상기 대상 영역의 히스토그램 분포를 가우시안 브로드닝(Gaussian broadening)하여 밝은 쪽의 최대값을 포함하는 분포 곡선을 필터링하여 상기 히스토그램 분포로부터 노이즈를 제거하는 연산을 수행하는 것을 특징으로 하는 초음파 영상으로부터 대상 장기의 지방량을 정량적으로 측정하는 장비.
KR10-2002-0005642A 2002-01-31 2002-01-31 초음파 영상으로부터 대상 장기의 지방량을 정량적으로 측정하는 장비 KR100438903B1 (ko)

Priority Applications (2)

Application Number Priority Date Filing Date Title
KR10-2002-0005642A KR100438903B1 (ko) 2002-01-31 2002-01-31 초음파 영상으로부터 대상 장기의 지방량을 정량적으로 측정하는 장비
US10/205,104 US6656121B2 (en) 2002-01-31 2002-07-24 Method of quantitatively measuring fat content in target organ from ultrasound visual image

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR10-2002-0005642A KR100438903B1 (ko) 2002-01-31 2002-01-31 초음파 영상으로부터 대상 장기의 지방량을 정량적으로 측정하는 장비

Publications (2)

Publication Number Publication Date
KR20030065155A KR20030065155A (ko) 2003-08-06
KR100438903B1 true KR100438903B1 (ko) 2004-07-02

Family

ID=27607058

Family Applications (1)

Application Number Title Priority Date Filing Date
KR10-2002-0005642A KR100438903B1 (ko) 2002-01-31 2002-01-31 초음파 영상으로부터 대상 장기의 지방량을 정량적으로 측정하는 장비

Country Status (2)

Country Link
US (1) US6656121B2 (ko)
KR (1) KR100438903B1 (ko)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101745502B1 (ko) * 2015-09-15 2017-06-13 경상북도 (관련부서:경상북도축산기술연구소장) 근내지방도 측정 시스템 및 이를 이용한 근내지방도 측정 방법
KR20220063096A (ko) * 2020-11-09 2022-05-17 지멘스 메디컬 솔루션즈 유에스에이, 인크. 물리적 양에 대한 초음파-도출된 프록시

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8094893B2 (en) * 2002-12-02 2012-01-10 Koninklijke Philips Electronics N.V. Segmentation tool for identifying flow regions in an image system
JP4137691B2 (ja) * 2003-04-30 2008-08-20 株式会社リコー 光記録媒体
KR100884745B1 (ko) * 2006-11-07 2009-02-20 한국전자통신연구원 초음파 영상으로부터 대상 장기의 조잡도를 측정하는 방법및 장치
KR100922724B1 (ko) 2006-12-06 2009-10-22 한국전자통신연구원 초음파 촬영 장비의 규준화를 위한 팬텀 장치 및 상기 팬텀장치를 이용한 초음파 영상 획득 방법
WO2011125549A1 (ja) * 2010-03-31 2011-10-13 公立大学法人大阪府立大学 脂肪組織画像表示装置
US8512247B2 (en) 2010-06-25 2013-08-20 John C. Hill System for non-invasive determination of glycogen stores
JP5943373B2 (ja) * 2010-10-26 2016-07-05 東芝メディカルシステムズ株式会社 超音波診断装置、超音波画像処理装置、医用画像診断装置、医用画像処理装置及び医用画像処理プログラム
US9579079B2 (en) 2011-01-05 2017-02-28 Koninklijke Philips Electronics N.V. Device and method for determining actual tissue layer boundaries of a body
JP2015181795A (ja) * 2014-03-25 2015-10-22 株式会社アドバンテスト 測定装置、測定方法、および超音波測定装置
US9642593B2 (en) 2014-09-19 2017-05-09 MuscleSound, LLC System and method for non-invasive determination of human body fat
US11013490B2 (en) 2016-11-15 2021-05-25 Musclesound, Inc. Non-invasive determination of muscle tissue size
US11064971B2 (en) 2016-11-30 2021-07-20 Musclesound, Inc. Non-Invasive determination of muscle tissue quality and intramuscular fat
US11096658B2 (en) 2017-02-02 2021-08-24 Musclesound, Inc. Non-invasive determination of pennation angle and/or fascicle length
US11160493B2 (en) 2017-03-03 2021-11-02 Musclesound, Inc. System and method for determining a subject's muscle fuel level, muscle fuel rating, and muscle energy status
CN113450713B (zh) * 2020-03-25 2022-08-12 北京小米移动软件有限公司 屏幕显示方法及装置、灰阶映射信息生成方法及装置
TWI779708B (zh) * 2021-07-07 2022-10-01 瑞昱半導體股份有限公司 基於顯示內容優化畫面的方法與相關的顯示器控制晶片

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR910007487A (ko) * 1989-10-17 1991-05-30 국립항공우주국 반사된 초음파 펄스의 특성을 나타내는 방법 및 장치
KR20010099267A (ko) * 2001-09-17 2001-11-09 오세택 휴대용 체지방 측정기
KR20030033318A (ko) * 2001-10-20 2003-05-01 (주)에이치비메디컬스 지방 흡입 및 지방층 측정 시스템

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5339815A (en) * 1992-12-22 1994-08-23 Cornell Research Foundation, Inc. Methods and apparatus for analyzing an ultrasonic image of an animal or carcass
US5398290A (en) * 1993-05-03 1995-03-14 Kansas State University Research Foundation System for measurement of intramuscular fat in cattle
US5941825A (en) * 1996-10-21 1999-08-24 Philipp Lang Measurement of body fat using ultrasound methods and devices
KR100318699B1 (ko) 1999-06-30 2001-12-28 강중길 초음파 스캐닝을 통한 골다공증 진단기

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR910007487A (ko) * 1989-10-17 1991-05-30 국립항공우주국 반사된 초음파 펄스의 특성을 나타내는 방법 및 장치
KR20010099267A (ko) * 2001-09-17 2001-11-09 오세택 휴대용 체지방 측정기
KR20030033318A (ko) * 2001-10-20 2003-05-01 (주)에이치비메디컬스 지방 흡입 및 지방층 측정 시스템

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101745502B1 (ko) * 2015-09-15 2017-06-13 경상북도 (관련부서:경상북도축산기술연구소장) 근내지방도 측정 시스템 및 이를 이용한 근내지방도 측정 방법
KR20220063096A (ko) * 2020-11-09 2022-05-17 지멘스 메디컬 솔루션즈 유에스에이, 인크. 물리적 양에 대한 초음파-도출된 프록시
KR102661893B1 (ko) 2020-11-09 2024-04-26 지멘스 메디컬 솔루션즈 유에스에이, 인크. 물리적 양에 대한 초음파-도출된 프록시

Also Published As

Publication number Publication date
US6656121B2 (en) 2003-12-02
US20030144592A1 (en) 2003-07-31
KR20030065155A (ko) 2003-08-06

Similar Documents

Publication Publication Date Title
KR100438903B1 (ko) 초음파 영상으로부터 대상 장기의 지방량을 정량적으로 측정하는 장비
US8605979B2 (en) Automatic detection and quantification of plaque in the coronary arteries of subjects from CT scans
Loizou et al. Snakes based segmentation of the common carotid artery intima media
JP2016531709A (ja) 疾患を診断するための画像解析技術
EP2117428B1 (en) Identification and analysis of lesions in medical imaging
JP4216496B2 (ja) 脳組織内毛細血管の血流動態に関するインデックス演算方法、装置及びプログラムコード
EP1302163A2 (en) Method and apparatus for calculating an index of local blood flows
US20020147392A1 (en) Technique for manipulating medical images
US8649843B2 (en) Automated calcium scoring of the aorta
JP2013542046A (ja) 超音波画像処理のシステムおよび方法
US20160035091A1 (en) Methods and Apparatus for Quantifying Inflammation
Huang et al. Stability analysis of fractal dimension in retinal vasculature
KR100440255B1 (ko) 대상 장기의 지방량 측정 시스템 및 지방량 측정프로그램을 기록한 기록매체
CN110236544B (zh) 基于相关系数的中风灌注成像病变区域检测系统及方法
Molinari et al. Accurate and automatic carotid plaque characterization in contrast enhanced 2-D ultrasound images
JP4714228B2 (ja) 脳組織内毛細血管の血流動態に関するインデックス演算方法、装置及び記憶媒体
Giannetti et al. Strain elastography in Crohn’s disease: the role of visual observation and semiquantitative parameters
US20080107317A1 (en) Method and apparatus for calculating coarseness level in target organ from ultrasonic image
KR100884745B1 (ko) 초음파 영상으로부터 대상 장기의 조잡도를 측정하는 방법및 장치
KR102258070B1 (ko) 발의 유형 평가 방법 및 이를 이용한 발의 유형 평가용 디바이스
Trombini et al. A digital image processing approach for hepatic diseases staging based on the glisson’s capsule
Simion et al. A Non-invasive Diagnosis Tool Based on Hepatorenal Index For Hepatic Steatosis
JP2005046527A (ja) 超音波肝画像を用いた定量的肝組織性状診断方法
Wu et al. A fully automated multi-modal computer aided diagnosis approach to coronary calcium scoring of MSCT images
CN115409834B (zh) 一种用于跟腱病诊断的特征提取方法、系统和存储介质

Legal Events

Date Code Title Description
A201 Request for examination
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: 20110609

Year of fee payment: 8

LAPS Lapse due to unpaid annual fee