[go: up one dir, main page]

KR100908248B1 - 탄성 영상 디스플레이 방법 - Google Patents

탄성 영상 디스플레이 방법 Download PDF

Info

Publication number
KR100908248B1
KR100908248B1 KR1020060088547A KR20060088547A KR100908248B1 KR 100908248 B1 KR100908248 B1 KR 100908248B1 KR 1020060088547 A KR1020060088547 A KR 1020060088547A KR 20060088547 A KR20060088547 A KR 20060088547A KR 100908248 B1 KR100908248 B1 KR 100908248B1
Authority
KR
South Korea
Prior art keywords
strain
displacement
elastic image
graph
ultrasound data
Prior art date
Legal status (The legal status 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 status listed.)
Expired - Fee Related
Application number
KR1020060088547A
Other languages
English (en)
Other versions
KR20080024327A (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 KR1020060088547A priority Critical patent/KR100908248B1/ko
Priority to EP07017841A priority patent/EP1900328B1/en
Priority to DE602007004543T priority patent/DE602007004543D1/de
Priority to US11/854,202 priority patent/US8133179B2/en
Priority to JP2007238204A priority patent/JP5596899B2/ja
Publication of KR20080024327A publication Critical patent/KR20080024327A/ko
Application granted granted Critical
Publication of KR100908248B1 publication Critical patent/KR100908248B1/ko
Priority to JP2014130256A priority patent/JP5855708B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Clinical applications
    • A61B8/0858Clinical applications involving measuring tissue layers, e.g. skin, interfaces
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/461Displaying means of special interest
    • A61B8/463Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/485Diagnostic techniques involving measuring strain or elastic properties
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object

Landscapes

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

Abstract

본 발명에 따른 탄성 영상 디스플레이 방법은, a) 복수의 조직들을 포함하는 대상체에 스트레스(stress) 인가 없이 초음파 신호를 다수의 스캔라인을 따라 송수신하여 제 1 초음파 데이터를 얻는 단계; b) 대상체에 소정의 스트레스를 인가하고 초음파 신호를 상기 다수의 스캔라인을 따라 송수신하여 제 2 초음파 데이터를 얻는 단계; c) 상기 제 1 초음파 데이터 및 상기 제 2 초음파 데이터를 이용하여 상기 복수의 조직들의 이동 변위(displacement)를 계산하는 단계; d) 상기 계산된 변위를 이용하여 탄성 영상을 위한 스트레인을 계산하는 단계; e) 상기 스트레인을 이용하여 탄성 영상을 형성하고 상기 스캔라인 별로 상기 조직들의 변화를 표시하는 그래프를 형성하는 단계; 및 f) 상기 탄성 영상 및 상기 그래프를 디스플레이하는 단계를 포함한다.
탄성 영상, 변위(displacement), 스트레스(stress), 스트레인(strain), 스캔라인(scanline), IIR(Infinite Impulse Response) 필터

Description

탄성 영상 디스플레이 방법{METHOD FOR DISPLAYING AN ELASTIC IMAGE}
도 1은 대상체에 스트레스의 인가 전, 후의 초음파 신호의 모양을 보여주는 개략도.
도 2a 내지 도 2c는 깊이에 따른 변위를 보여주는 그래프.
도 3은 스트레인 계산을 설명하기 위한 스프링 모델을 보이는 개략도.
도 4는 본 발명의 일실시예에 따라 탄성영상과 스트레인 그래프를 디스플레이한 예를 보여주는 개략도.
도 5a 내지 도 5e는 본 발명의 일실시예에 따라서 탄성 영상의 한 프레임에서의 각 스캔라인 별로 스트레인 값을 보여주는 막대 그래프의 예들을 보여주는 예시도.
도 6a 내지 도 6d는 본 발명의 다른 실시예에 따른 스트레인 그래프를 보여주는 예시도.
도 7은 본 발명의 일실시예에 따라 탄성 영상을 디스플레이하는 방법을 보여주는 흐름도.
도 8은 스트레인 그래프를 시간 축으로 확장하여 3차원으로 보여주는 그래프.
도 9는 초음파 프로브에 2차원 어레이를 사용하여 모든 스캔라인에 대해서 스트레스를 3차원으로 보여주는 그래프.
본 발명은 초음파 시스템에 관한 것으로, 더욱 상세하게는 초음파 시스템에서 탄성 영상과 대상체 조직의 변화를 표시하는 그래프를 디스플레이하는 방법에 관한 것이다.
초음파 시스템은 대상체에 초음파 신호를 송신하고, 대상체의 불연속면에서 반사되어 되돌아오는 초음파 신호를 수신하고, 수신된 초음파 신호를 전기적 신호로 변환하여 소정의 영상 장치를 통하여 출력함으로써 대상체의 내부 상태를 진단한다.
초음파 영상은 조직 사이의 임피던스 차이에 의한 반사계수를 이용하는 B-모드(Brightness-mode)로써 주로 표현된다. 그러나, 종양이나 암 조직과 같이 주위의 조직과 비교하여 반사계수가 차이 나지 않는 부분은 B-모드 영상에서 구별하기 쉽지 않다. 이에 반하여, 초음파 탄성 영상 기법은 조직의 기계적인 성질을 영상화하므로 암 조직과 같은 병변의 진단에 큰 도움을 준다. 예를 들어, 종양이나 암조직 은 주위의 연조직에 비해서 조직이 단단하므로 외부에서 같은 힘을 가했을 때 주위 조직보다 변형 정도가 작다.
조직의 탄성은 대상체의 조직을 압축하기 전 얻은 초음파 데이터와 조직을 압축하여 얻은 초음파 데이터를 이용하여 측정한다. 일반적으로 초음파 프로브에 장착된 압축판을 이용하여 사용자가 압력을 가하여 조직을 압축한다. 이 경우, 조직의 변형 정도는 사용자에 의해 가해지는 압축의 크기에 의해서 결정된다. 탄성 영상의 화질은 대상체에 가하는 압축 크기와 압축 속도에 따라 달라진다. 예를 들어 압축 크기가 작을 경우 종양이나 암 조직의 변형과 주변 조직의 변형의 차이가 뚜렷하게 나타나지 않고, 압축 정도가 큰 경우 압축에 의한 종양이나 암 조직과 주변조직의 비상관성이 커져 탄성 영상의 화질이 저하된다. 따라서, 좋은 화질의 탄성 영상을 얻기 위해서는 적당한 크기의 압축이 요구된다. 보통 압축에 의해서 조직의 변형이 0.5~3%가 될 때 탄성 영상의 질이 가장 좋다고 알려져 있다. 압축의 정도는 사용자마다 또는 측정시마다 다르게 적용되므로 어느 정도의 압축을 가해야 하는 지를 알려주는 지표가 필요하다.
본 발명은 대상체에 스트레스 인가 전, 후에 얻은 초음파 데이터에서 계산된 각 스캔라인에 대한 대상체 조직의 최대 변위를 표시하는 변위 그래프 또는 최대 변위를 이용하여 계산된 각 스캔라인의 스트레인 표시하는 스트레인 그래프를 탄성 영상과 함께 디스플레이하는 방법을 제공한다.
상기 목적을 달성하기 위한 본 발명에 따른 탄성 영상 디스플레이 방법은, a) 복수의 조직들을 포함하는 대상체에 스트레스(stress) 인가 없이 초음파 신호를 다수의 스캔라인을 따라 송수신하여 제 1 초음파 데이터를 얻는 단계; b) 대상체에 소정의 스트레스를 인가하고 초음파 신호를 상기 다수의 스캔라인을 따라 송수신하 여 제 2 초음파 데이터를 얻는 단계; c) 상기 제 1 초음파 데이터 및 상기 제 2 초음파 데이터를 이용하여 상기 복수의 조직들의 이동 변위(displacement)를 계산하는 단계; d) 상기 계산된 변위를 이용하여 탄성 영상을 위한 스트레인을 계산하는 단계; e) 상기 스트레인을 이용하여 탄성 영상을 형성하고 상기 스캔라인 별로 상기 조직들의 변화를 표시하는 그래프를 형성하는 단계; 및 f) 상기 탄성 영상 및 상기 그래프를 디스플레이하는 단계를 포함한다.
이하, 본 발명에 따른 탄성 영상 및 스트레인 그래프를 표시하는 방법에 대해서 첨부된 도면을 참조하여 설명한다.
도 1은 대상체에 스트레스 인가 전, 후의 초음파 신호 모양을 보이고 있다.
먼저, 단위 면적당 가해지는 힘인 스트레스(stress)를 인가하기 전에 대상체에 초음파 신호를 인가하여 기준신호(102)를 얻고, 대상체의 표면에 스트레스를 가하여 대상체를 압축한 후 초음파 신호(103)를 얻는다. 본 발명에 따른 스트레스는 초음파 프로브를 이용하여 대상체에 인가할 수 있다.
대상체에 압축을 가하면 대상체 내의 반사체들이 압축 방향으로 이동한다. 이러한 반사체의 이동에 따라 압축되기 전과 비교할 때 초음파 수신 신호의 이동이 나타난다. 따라서, 두 신호 사이의 이동을 계산함으로써 매질의 변위를 구할 수 있다. 이러한 변위는 대상체의 단단함에 따라 달라지므로 매질의 특성 값을 반영한다.
한편, 도 1에 보이는 바와 같이, 트랜스듀서(101)에 가까운 곳에서는 신호의 이동이 작지만 먼 곳(깊이 방향)에서는 변위가 누적되어 신호 간 이동이 크게 나타 난다.
조직에 일정한 압력을 한 방향으로 가하면 조직의 단단한 정도에 따라 변형된 정도가 다르게 나타난다. 따라서 외부 또는 내부에서 인가된 힘에 대하여 매질의 이동 변위(displacement)를 계산하고, 변위의 함수를 미분하여 기울기를 구하면 변형률 즉, 스트레인(strain)을 구할 수 있다. 이 스트레인 값에 기초하여 탄성영상을 구성한다.
스트레인을 계산하기 위한 이동 변위는 RF 데이터 또는 기저대역(baseband)의 IQ 데이터를 이용하여 교차상관(cross correlation)법 또는 자기상관(auto correlation)법으로 계산할 수 있다. 일반적으로, 초음파 송수신을 통하여 얻은 RF 또는 IQ 데이터는 샘플링하여 얻은 샘플링 데이터이므로 정확한 이동 변위를 계산하는 것이 쉽지않다. 따라서, 본 발명의 일실예에서는 샘플링 데이터를 이용하여 샘플링 간격 단위로 제 1 이동 변위(coarse displacement)를 계산하고, 샘플링 데이터를 보간(interpolation)하여 얻은 데이터에서 샘플링 간격 단위 내에서 제 2 이동 변위(fine displacement)를 계산한다. 도 2a는 본 발명에 따른 일실시예에서 깊이에 따른 이동 변위를 보여 주는 그래프이고, 도 2b는 샘플링 간격 단위 내에서의 깊이에 따른 이동 변위를 보여 주는 그래프이다. 스트레스 인가에 따른 대상체의 전체 변위는 제 1 이동 변위와 제 2 이동 변위의 합으로 결정되며, 도 2c는 깊이에 따른 제 1 이동 변위와 제 2 이동 변위를 합한 깊이 따른 조직의 전체 변위를 보여 준다.
스트레스의 인가에 따른 대상체 내의 조직의 스트레인(ε)은 도 3에 도시된 바와 같이 스트레스 인가 전의 대상체의 표면에서 임의의 깊이까지의 길이(L)와 스트레스 인가 후의 길이 변화(△L)의 비율을 이용하여 계산한다.
ε=△L/L
앞에서 언급한 방법으로 탄성 영상의 한 프레임에 대하여 최대 변위를 이용하여 각 스캔라인의 스트레인을 계산하고, 각 스캔라인 별로 스트레인 값을 나타내는 그래프를 탄성 영상과 함께 디스플레이한다. 본 발명의 일실시예에따라 탄성 영상과 함께 스트레인 그래프를 디스플레이하였지만, 본 발명의 다른 실시예에서는 탄성 영상과 함께 각 스캔라인 별로 최대변위를 나타내는 변위 그래프를 디스플레이할 수 있다.
도 4는 탄성 영상과 함께 각 스캔라인 별로 스트레인 값을 나타내는 스트레인 그래프를 디스플레이한 예를 보여주는 개략도이다. 도 4를 참조하면, 대상체 내에서 상대적으로 단단한 제 1 조직(410)과 상대적으로 연한 제 2 조직(420)에 대한 탄성 영상(400)에서의 각 스캔라인 별로 스트레인 값을 보여는 스트레인 그래프로서 막대 그래프(430)가 디스플레이된다. 도 4의 막대 그래프에서 보이는 바와 같이, 상대적으로 단단한 제 1 조직(410)이 있는 부분에 대해서는 스트레인 값을 나타내는 막대의 크기 짧게 표시되며, 상대적으로 연한 조직인 제 1 조직(410)의 주변 조직에서의 스트레인 값을 나타내는 막대의 크기는 길게 표시된다.
한편, 최상의 화질을 갖는 탄성 영상을 얻기 위해서 평균 스트레인 값을 0.5 내지 3%로 조절하는 것이 바람직하다. 평균 스트레인은 본 발명에 따라 각 스캔라 인 별로 계산된 스트레인을 이용하여 계산되고, 이렇게 계산된 평균 스트레인 값을 탄성 영상이 제공되는 화면상에 표시하여 사용자가 적절한 스트레인 값이 되도록 스트레스를 조절할 수 있다.
도 5a 내지 도 5e는 본 발명의 일실시예에 따라서 탄성 영상의 한 프레임에서의 각 스캔라인 별로 스트레인 값을 보여주는 막대 그래프의 예들을 보여주는 예시도이다. 도 5a에 도시된 막대 그래프를 보면, 막대 그래프에서 가운데 부분에 있는 막대의 크기가 작게 양쪽 옆으로 갈수록 막대의 크기가 크게 도시된 것으로부터 대상체를 구성하는 조직의 가운데 부분에 상대적으로 단단한 조직이 존재하는 것을 알 수 있다. 또한, 도 5b에 도시된 막대 그래프로부터 대상체를 구성하는 조직의 가운데 부분에 상대적으로 연한 조직이 존재하는 것을 알 수 있다. 도 5c는 대상체에 스트레스가 오른쪽보다 왼쪽에 더 세게 인가되었을 때의 스트레인 값을 보여주는 막대 그래프이며, 도 5d는 대상체에 스트레스가 왼쪽보다 오른쪽에 더 세게 인가되었을 때의 스트레인 값을 보여준다. 한편, 도 5e는 각 스캔라인 별로 스트레인을 계산하는 과정에서 오차가 발생하는 경우의 불규칙적인 모양으로 나타나는 스트레인 값을 보여준다. 도 5e에 보이는 바와 같이, 스트레인 그래프를 통하여 탄성 영상을 얻기 위한 영상 처리가 적절히 처리되고 있는 지를 확인 할 수 있다.
본 발명의 일실시예에서는 각 스캔라인의 스트레인을 나타내는 막대 그래프의 아래쪽을 동일하게 고정하고 위쪽으로 스트레인의 크기를 표시하였지만, 본 발명의 다른 실시예에서는 도 6a 내지 6d에 도시된 바와 같이 스트레인 값을 나타내는 막대의 위쪽을 동일하게 고정하고 아래쪽으로 스트레인의 크기를 표시하여 사용 자가 초음파 프로브를 이용하여 인가하는 스트레스 크기를 직관적으로 인식하도록 할 수 있다.
한편, 초음파 프로브를 이용하여 대상체를 압축할 경우 초음파 프로브의 압축 면적이 상대적으로 좁아 대상체를 균일하게 압축하기 어려울 수 있다. 즉, 프로브가 유한한 면적을 가지고 있기 때문에 프로브의 가운데 부분에 비해 바깥쪽에 대응하는 대상체에 대해서 압축 전달이 제대로 이루어지지 않을 수 있다. 이런 경우 대상체에 균일한 압력이 가해지지 않아 탄성 영상에 영향을 줄 수 있으므로 각 스캔라인의 스트레인 값을 이용하여 각 스트레인 값들을 평준화(normalize)해주어 프로브의 유한 플레이트(plate)에 의한 압력의 불균일성을 보상할 수 있다.
또한, 탄성 영상은 실시간으로 보여주는 것이므로 대상체에 스트레스가 균일하게 인가되지 못할 경우 화질이 떨어지는 탄성 영상이 출력될 수 있다. 종래의 탄성 영상에서는 탄성 영상의 화질이 떨어지는 프레임을 삭제하여 보여주지 않거나 이전 프레임을 연속으로 보여주었다. 하지만 이런 경우 화면이 깜박이거나 B-모드 영상과 동기가 맞지 않는 문제가 있다. 따라서, 본 발명의 일실시예에서는 수학식2를 이용한 IIR(Infinite Impulse Response) 필터를 이용하여 탄성영상을 출력한다.
YN = (1-P)YN -1 + PXN
YN은 현재 디스플레이되는 탄성 영상의 프레임을 나타내고, YN -1은 이전에 디스플레이된 탄성 영상의 프레임을 나타내며, XN은 현재 스트레인이 계산된 프레임을 나타내며, P는 퍼시스턴스(persistence) 값이다. P값이 커지면 현재 스트레인이 계산된 프레임에 대한 비중이 커져 보이는 것이고, P값이 작아지면 현재 스트레인 값이 계산된 프레임이 반영되는 부분이 작아져 이전에 디스플레이된 탄성 영상이 거의 그대로 디스플레이 되게 된다. 프레임당 평균 스트레인의 값을 계산하여 이전에서 적정한 스트레인의 범위(0.5~3.0%)에 평균 스트레인 값이 포함되면 P의 값을 크게하여 현재 스트레인이 계산된 프레임의 비중을 크게 하여 출력하고, 평균 스트레인 값이 적정한 스트레인의 범위를 벗어나면 P의 값을 작게 하여 이전에 디스플레이된 프레임의 비중을 크게 하여 출력하여 화질이 개선된 실시간 탄성 영상을 디스플레이할 수 있다.
도 7은 본 발명에 따라 탄성 영상을 디스플레이하는 방법을 보여주는 흐름도 이다. 도 7을 참조하며, 대상체에 스트레스를 인가하기 전에 초음파 신호를 송수신하여 제 1 초음파 데이터를 얻고, 대상체에 스트레스를 인가하여 압축한 후 초음파 신호를 송수신하여 제 2 초음파 데이터를 얻는다(S710). 제 1 초음파 데이터와 제 2 초음파 데이터를 샘플링 간격 단위로 비교하여 제 1 변위(coarse displacement)를 계산하고, 제 1 초음파 데이터와 제 2 초음파 데이터를 보간한 후 샘플링 간격 단위 내에서 데이터를 비교하여 제 2 변위(fine displacement)를 계산한다(S720). 제 1 변위와 제 2 변위를 합하여 최종 변위를 계산하고, 계산된 최종 변위를 이용하여 대상체의 조직들의 스트레인을 계산한다(S730). 이후, 계산된 최종 변위 중 최대 변위를 이용하여 탄성 영상의 프레임을 형성하기 위한 스캔라인 별로 스트레인을 계산한다(S740).
대상체에 인가되는 스트레스의 불균일성으로 인한 탄성 영상의 영향을 줄이기 위해서 스트레인 값을 이용하여 각 스캔라인 별로 스트레인 값을 평준화(normalization)를 한다(S750). 한 프레임을 구성하는 스트레인 값의 평균값을 계산하고, 계산된 스트레인 평균값이 소정의 범위, 예를 들어 0.5% 내지 3% 내에 해당하는 지를 확인한다(S760). 스트레인 평균값이 소정의 범위 내에 해당하지 않으면 대상체에 인가되는 스트레스를 조절하여 스트레인 평균값이 소정의 범위 내에 해당하도록 한다. 스트레인 별로 컬러를 미리 설정하고, 평균 스트레인 값이 소정 범위에 해당하면 각 스트레인 값에 대응하는 칼라를 초음파 데이터 및 스트레인 값을 표시하는 그래프에 매핑(mapping)을 한다(S770).
계속해서, 스트레스가 균일하게 인가되지 못하여 화질이 떨어지는 탄성 영상이 출력되는 것을 막기 위해서 IIR 필터의 P값을 조절하여 탄성 영상을 필터링한다(S780). 이후, 본 발명에 따른 탄성 영상과 스트레인을 표시하는 그래프를 디스플레이한다(S790).
한편, 본 발명의 일실시예에서는 각 스캔라인 별 스트레인을 나타내는 막대 그래프를 2차원의 그래프로 표시하였지만, 본 발명의 다른 실시예에서는 도 8에 도시된 것과 같이 각 스캔라인 별로 스트레인을 시간축 상의 3차원으로 표시할 수 있다. 도 8과 같이 스트레인을 표시하여 시간 변화에 따른 압축의 크기와 그에 따른 스트레인의 크기를 볼 수 있다. 도 9는 초음파 프로브에 2차원 어레이를 사용하는 경우 압축의 크기를 구성되는 모든 스캔라인에 대해서 3차원으로 보여주는 그래프이다. 본 발명의 일실시예에서는 탄성 영상에서 각 스캔라인별 변위의 크기 또는 스트레인 크기를 막대 그래프를 일 예로 이용하여 표시하였지만 선 그래프 등 변위 또는 스트레인 크기를 표시할 수 있는 임의의 그래프를 이용할 수 있다.
상술한 실시예는 본 발명의 원리를 응용한 다양한 실시예의 일부를 나타낸 것에 지나지 않음을 이해해야 한다. 본 기술 분야에서 통상의 지식을 가진 자는 본 발명의 본질로부터 벗어남이 없이 여러 가지 변형이 가능함을 명백히 알 수 있을 것이다.
전술한 바와 같이, 본 발명에 따라 탄성 영상을 구성하기 위한 스트레인을 계산하여 탄성영상과 함께 각 스캔라인의 최대변위에 대한 스트레인을 그래프로 표시함으로써, 각 영상프레임의 평균 스트레인을 계산할 수 있으므로 현재 인가되고 있는 스트레스가 적정한 범위에 들어가는 것인지 판단하여 사용자가 압축의 크기를 조절할 수 있고, 평균 스트레인의 값이 적정 범위에 들어가는지 벗어나는지에 따라 퍼시스턴스(persistent) 값을 변경하여 최적의 탄성 영상 출력할 수 있다.
또한, 스트레인의 그래프 모양을 가지고 대략적인 대상체의 탄성 특성을 알 수 있고, 한쪽에 치우침없는 균일한 압력을 가하도록 유도할 수 있다. 그리고 각 스캔라인의 스트레인 값으로 스캔라인에 따라 다른 평준화를 취함으로써 유한 압축 플레이트에 의한 압축의 불균일을 보상할 수 있다.

Claims (10)

  1. 탄성 영상 디스플레이 방법에 있어서,
    a) 복수의 조직들을 포함하는 대상체에 스트레스(stress) 인가 없이 초음파 신호를 다수의 스캔라인을 따라 송수신하여 제 1 초음파 데이터를 얻는 단계;
    b) 대상체에 소정의 스트레스를 인가하고 초음파 신호를 상기 다수의 스캔라인을 따라 송수신하여 제 2 초음파 데이터를 얻는 단계;
    c) 상기 제 1 초음파 데이터 및 상기 제 2 초음파 데이터를 이용하여 상기 복수의 조직들의 이동 변위(displacement)를 계산하는 단계;
    d) 상기 계산된 변위를 이용하여 탄성 영상을 위한 스트레인을 계산하는 단계;
    e) 상기 스트레인을 이용하여 탄성 영상을 형성하고 상기 스캔라인 별로 상기 조직들의 변화를 표시하는 그래프를 형성하는 단계; 및
    f) 상기 탄성 영상 및 상기 그래프를 디스플레이하는 단계
    를 포함하는 탄성 영상 디스플레이 방법.
  2. 제 1 항에 있어서,
    상기 그래프는 상기 계산된 변위 중 상기 각 스캔라인 별로 최대 변위를 표시하는 변위 그래프인 탄성 영상 디스플레이 방법.
  3. 제 1 항에 있어서,
    상기 그래프는 상기 계산된 변위를 이용하여 계산된 각 스캔라인 별로 스트레인을 표시하는 스트레인 그래프인 탄성 영상 디스플레이 방법.
  4. 제 1 항에 있어서, 상기 단계 c)는,
    c1) 상기 제 1 초음파 데이터와 상기 제 2 초음파 데이터로부터 샘플링 간격 단위로 제 1 변위를 계산하는 단계;
    c2) 상기 제 1 초음파 데이터 및 상기 제 2 초음파 데이터를 보간하는 단계;
    c3) 상기 보간된 제 1 초음파 데이터 및 상기 제 2 초음파 데이터로부터 샘플링 간격 단위 내에서의 제 2 변위를 계산하는 단계; 및
    c4) 상기 제 1 변위와 상기 제 2 변위를 합하여 상기 이동 변위를 계산하는 단계
    를 포함하는 탄성 영상 디스플레이 방법.
  5. 제 4 항에 있어서, d1) 상기 단계 d)에서 계산된 스트레인을 이용하여 상기 각 스캔라인 별로 스트레인을 평준화(normalization) 하는 단계를 더 포함하는 탄성 영상 디스플레이 방법.
  6. 제 5 항에 있어서,
    d2) 상기 각 스캔라인의 스트레인의 평균값을 계산하는 단계;
    d3) 상기 계산된 평균 스트레인이 소정의 범위에 해당하는 지를 판단하는 단계; 및
    d4) 상기 평균 스트레인이 소정의 범위에 해당하지 않으면 상기 대상체에 인가되는 스트레스를 조절하는 단계를
    더 포함하는 탄성 영상 디스플레이 방법.
  7. 제 2 항에 있어서, 상기 단계 e)는 스트레인 별로 사전 설정된 컬러에 따라서 상기 계산된 최대 변위에 해당하는 컬러를 매핑하여 탄성 영상을 형성하는 탄성 영상 디스플레이 방법.
  8. 제 3 항에 있어서, 상기 단계 e)는 스트레인 별로 사전 설정된 컬러에 따라서 상기 계산된 스트레인에 해당하는 컬러를 매핑하여 탄성 영상을 형성하는 탄성 영상 디스플레이 방법.
  9. 제 1 항에 있어서, 상기 그래프는 막대 또는 선 그래프인 탄성 영상 디스플레이 방법.
  10. 제 9 항에 있어서, 상기 탄성 영상을 디스플레이하기 전 스트레인 범위에 따라 다른 값의 P(Persistence)값을 적용한 IIR(Infinite Impulse Response) 필터를 이용하여 디스플레이하는 단계를 더 포함하는 탄성 영상 디스플레이 방법.
KR1020060088547A 2006-09-13 2006-09-13 탄성 영상 디스플레이 방법 Expired - Fee Related KR100908248B1 (ko)

Priority Applications (6)

Application Number Priority Date Filing Date Title
KR1020060088547A KR100908248B1 (ko) 2006-09-13 2006-09-13 탄성 영상 디스플레이 방법
EP07017841A EP1900328B1 (en) 2006-09-13 2007-09-12 Method of displaying an elastic image
DE602007004543T DE602007004543D1 (de) 2006-09-13 2007-09-12 Verfahren zur Anzeige eines elastischen Bildes
US11/854,202 US8133179B2 (en) 2006-09-13 2007-09-12 Method of displaying an elastic image
JP2007238204A JP5596899B2 (ja) 2006-09-13 2007-09-13 弾性映像ディスプレイ方法
JP2014130256A JP5855708B2 (ja) 2006-09-13 2014-06-25 弾性映像ディスプレイ方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020060088547A KR100908248B1 (ko) 2006-09-13 2006-09-13 탄성 영상 디스플레이 방법

Publications (2)

Publication Number Publication Date
KR20080024327A KR20080024327A (ko) 2008-03-18
KR100908248B1 true KR100908248B1 (ko) 2009-07-20

Family

ID=38657119

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020060088547A Expired - Fee Related KR100908248B1 (ko) 2006-09-13 2006-09-13 탄성 영상 디스플레이 방법

Country Status (5)

Country Link
US (1) US8133179B2 (ko)
EP (1) EP1900328B1 (ko)
JP (2) JP5596899B2 (ko)
KR (1) KR100908248B1 (ko)
DE (1) DE602007004543D1 (ko)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101085220B1 (ko) 2009-05-14 2011-11-21 삼성메디슨 주식회사 탄성 영상 구현 장치 및 방법
KR101097651B1 (ko) * 2010-10-01 2011-12-22 삼성메디슨 주식회사 전역 균일 신장에 기초하여 탄성영상을 제공하는 초음파 시스템 및 방법

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100908248B1 (ko) * 2006-09-13 2009-07-20 주식회사 메디슨 탄성 영상 디스플레이 방법
US8100831B2 (en) * 2006-11-22 2012-01-24 General Electric Company Direct strain estimator for measuring elastic properties of tissue
KR100983770B1 (ko) * 2007-03-23 2010-09-28 주식회사 메디슨 초음파 탄성영상을 형성하기 위한 시스템 및 방법
EP2107389B1 (en) * 2008-04-02 2014-05-14 Samsung Medison Co., Ltd. Formation of an elastic image in an ultrasound system
EP2296552A1 (en) * 2008-05-29 2011-03-23 Koninklijke Philips Electronics N.V. Tissue strain analysis
JP5400773B2 (ja) * 2008-06-16 2014-01-29 株式会社日立メディコ 超音波診断装置、超音波画像表示方法及び超音波診断プログラム
KR101014564B1 (ko) * 2008-06-26 2011-02-16 주식회사 메디슨 탄성 영상을 형성하기 위한 초음파 시스템 및 방법
KR101014558B1 (ko) * 2008-07-16 2011-02-16 주식회사 메디슨 스트레인 정규화를 통한 탄성영상 형성 방법 및 그를 위한초음파 시스템
KR101100498B1 (ko) 2008-08-05 2011-12-29 삼성메디슨 주식회사 컬러맵을 형성하는 초음파 시스템 및 방법
KR101134842B1 (ko) * 2008-11-20 2012-04-13 삼성메디슨 주식회사 탄성영상에 적응적으로 퍼시스턴스를 적용하는 초음파 시스템
WO2011052400A1 (ja) * 2009-10-28 2011-05-05 株式会社 日立メディコ 超音波診断装置及び画像構成方法
WO2011102401A1 (ja) * 2010-02-17 2011-08-25 株式会社 日立メディコ 弾性画像の画質評価方法及び超音波診断装置
JP5677757B2 (ja) * 2010-03-16 2015-02-25 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波診断装置
EP2598907B1 (en) * 2010-07-29 2017-03-01 B-K Medical ApS Motion-compensated processing
KR101390185B1 (ko) 2010-08-23 2014-04-29 삼성메디슨 주식회사 탄성영상을 제공하는 초음파 시스템 및 방법
CN104135937B (zh) 2012-02-21 2017-03-29 毛伊图像公司 使用多孔超声确定材料刚度
US10675007B2 (en) * 2016-04-19 2020-06-09 Siemens Medical Solutions Usa, Inc. Frequency compounding in elasticity imaging
CN113397600A (zh) * 2021-06-18 2021-09-17 复旦大学 基于超声射频信号弹性重建的肝纤维化评估方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20010067091A (ko) * 1999-08-23 2001-07-12 추후보정 초음파 이미징에 있어 조직 변형의 실시간 계산 및디스플레이를 제공하는 방법 및 장치
KR20040066795A (ko) * 2001-10-17 2004-07-27 스펙트라소닉스 이미징 인코포레이티드 신체 조직의 기계적 강도 특성을 나타내는 장치 및 방법
KR20070013986A (ko) * 2005-07-27 2007-01-31 주식회사 메디슨 탄성영상 형성 및 디스플레이를 위한 초음파 시스템

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6270459B1 (en) * 1998-05-26 2001-08-07 The Board Of Regents Of The University Of Texas System Method for estimating and imaging of transverse displacements, transverse strains and strain ratios
US6558324B1 (en) 2000-11-22 2003-05-06 Siemens Medical Solutions, Inc., Usa System and method for strain image display
US6508768B1 (en) * 2000-11-22 2003-01-21 University Of Kansas Medical Center Ultrasonic elasticity imaging
US7297116B2 (en) * 2003-04-21 2007-11-20 Wisconsin Alumni Research Foundation Method and apparatus for imaging the cervix and uterine wall
JP2005013283A (ja) * 2003-06-23 2005-01-20 Takeshi Shiina 超音波探触子及び超音波診断装置
JP4455003B2 (ja) * 2003-10-14 2010-04-21 株式会社日立メディコ 超音波診断装置
JP4515799B2 (ja) * 2004-03-24 2010-08-04 株式会社日立メディコ 超音波診断装置
US7118532B2 (en) * 2004-05-06 2006-10-10 General Electric Company Data dependent color wall filters
CN100518662C (zh) 2004-06-09 2009-07-29 株式会社日立医药 超声波诊断装置
US8734351B2 (en) * 2004-08-05 2014-05-27 Hitachi Medical Corporation Method of displaying elastic image and diagnostic ultrasound system
WO2006022238A1 (ja) * 2004-08-25 2006-03-02 Hitachi Medical Corporation 超音波診断装置
WO2006040967A1 (ja) * 2004-10-08 2006-04-20 Hitachi Medical Corporation 超音波診断装置
US9389203B2 (en) * 2005-06-17 2016-07-12 Wisconsin Alumni Research Foundation Automated ultrasonic elasticity image formation with quality measure
KR100898946B1 (ko) * 2006-01-24 2009-05-25 주식회사 메디슨 탄성영상 형성을 위한 초음파 진단 시스템 및 방법
KR100908248B1 (ko) * 2006-09-13 2009-07-20 주식회사 메디슨 탄성 영상 디스플레이 방법

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20010067091A (ko) * 1999-08-23 2001-07-12 추후보정 초음파 이미징에 있어 조직 변형의 실시간 계산 및디스플레이를 제공하는 방법 및 장치
KR20040066795A (ko) * 2001-10-17 2004-07-27 스펙트라소닉스 이미징 인코포레이티드 신체 조직의 기계적 강도 특성을 나타내는 장치 및 방법
KR20070013986A (ko) * 2005-07-27 2007-01-31 주식회사 메디슨 탄성영상 형성 및 디스플레이를 위한 초음파 시스템

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101085220B1 (ko) 2009-05-14 2011-11-21 삼성메디슨 주식회사 탄성 영상 구현 장치 및 방법
KR101097651B1 (ko) * 2010-10-01 2011-12-22 삼성메디슨 주식회사 전역 균일 신장에 기초하여 탄성영상을 제공하는 초음파 시스템 및 방법

Also Published As

Publication number Publication date
DE602007004543D1 (de) 2010-03-18
JP5855708B2 (ja) 2016-02-09
JP2008068094A (ja) 2008-03-27
JP5596899B2 (ja) 2014-09-24
EP1900328B1 (en) 2010-01-27
JP2014168721A (ja) 2014-09-18
EP1900328A1 (en) 2008-03-19
KR20080024327A (ko) 2008-03-18
US8133179B2 (en) 2012-03-13
US20080064956A1 (en) 2008-03-13

Similar Documents

Publication Publication Date Title
KR100908248B1 (ko) 탄성 영상 디스플레이 방법
KR101132531B1 (ko) 서로 마주 보는 트랜스듀서를 구비하는 초음파 진단 장치
JP4455003B2 (ja) 超音波診断装置
US8043216B2 (en) Method of displaying elastic image and diagnostic ultrasound system
US7628754B2 (en) Ultrasonographic device
US6558324B1 (en) System and method for strain image display
EP1747757B1 (en) Ultrasound system for displaying an elastic image
JP2005066041A (ja) 超音波探触子及び超音波診断装置
US20070165925A1 (en) Image processing system and method of enhancing the quality of an ultrasound image
JP2005066041A5 (ko)
EP3513734A1 (en) Ultrasonic imaging apparatus and control method thereof
JP5473527B2 (ja) 超音波診断装置
EP2157442B1 (en) Formation of an elastic image in an ultrasound system
KR101014558B1 (ko) 스트레인 정규화를 통한 탄성영상 형성 방법 및 그를 위한초음파 시스템
KR101286331B1 (ko) 탄성 영상을 제공하는 초음파 시스템 및 방법
JP5623609B2 (ja) 超音波診断装置
KR100847796B1 (ko) 탄성영상 형성 및 디스플레이를 위한 초음파 시스템
KR101107392B1 (ko) 가이드 정보를 제공하는 초음파 시스템 및 방법
KR101562209B1 (ko) 서로 마주 보는 트랜스듀서를 구비하는 초음파 진단 장치

Legal Events

Date Code Title Description
PA0109 Patent application

St.27 status event code: A-0-1-A10-A12-nap-PA0109

A201 Request for examination
P11-X000 Amendment of application requested

St.27 status event code: A-2-2-P10-P11-nap-X000

P13-X000 Application amended

St.27 status event code: A-2-2-P10-P13-nap-X000

PA0201 Request for examination

St.27 status event code: A-1-2-D10-D11-exm-PA0201

D13-X000 Search requested

St.27 status event code: A-1-2-D10-D13-srh-X000

PG1501 Laying open of application

St.27 status event code: A-1-1-Q10-Q12-nap-PG1501

R18-X000 Changes to party contact information recorded

St.27 status event code: A-3-3-R10-R18-oth-X000

D14-X000 Search report completed

St.27 status event code: A-1-2-D10-D14-srh-X000

E701 Decision to grant or registration of patent right
PE0701 Decision of registration

St.27 status event code: A-1-2-D10-D22-exm-PE0701

GRNT Written decision to grant
PR0701 Registration of establishment

St.27 status event code: A-2-4-F10-F11-exm-PR0701

PR1002 Payment of registration fee

St.27 status event code: A-2-2-U10-U11-oth-PR1002

Fee payment year number: 1

PG1601 Publication of registration

St.27 status event code: A-4-4-Q10-Q13-nap-PG1601

PN2301 Change of applicant

St.27 status event code: A-5-5-R10-R13-asn-PN2301

St.27 status event code: A-5-5-R10-R11-asn-PN2301

PN2301 Change of applicant

St.27 status event code: A-5-5-R10-R13-asn-PN2301

St.27 status event code: A-5-5-R10-R11-asn-PN2301

PN2301 Change of applicant

St.27 status event code: A-5-5-R10-R13-asn-PN2301

St.27 status event code: A-5-5-R10-R11-asn-PN2301

PR1001 Payment of annual fee

St.27 status event code: A-4-4-U10-U11-oth-PR1001

Fee payment year number: 4

FPAY Annual fee payment

Payment date: 20130521

Year of fee payment: 5

PR1001 Payment of annual fee

St.27 status event code: A-4-4-U10-U11-oth-PR1001

Fee payment year number: 5

FPAY Annual fee payment

Payment date: 20140703

Year of fee payment: 6

PR1001 Payment of annual fee

St.27 status event code: A-4-4-U10-U11-oth-PR1001

Fee payment year number: 6

FPAY Annual fee payment

Payment date: 20150629

Year of fee payment: 7

PR1001 Payment of annual fee

St.27 status event code: A-4-4-U10-U11-oth-PR1001

Fee payment year number: 7

FPAY Annual fee payment

Payment date: 20160704

Year of fee payment: 8

PR1001 Payment of annual fee

St.27 status event code: A-4-4-U10-U11-oth-PR1001

Fee payment year number: 8

FPAY Annual fee payment

Payment date: 20170704

Year of fee payment: 9

PR1001 Payment of annual fee

St.27 status event code: A-4-4-U10-U11-oth-PR1001

Fee payment year number: 9

FPAY Annual fee payment

Payment date: 20180703

Year of fee payment: 10

PR1001 Payment of annual fee

St.27 status event code: A-4-4-U10-U11-oth-PR1001

Fee payment year number: 10

FPAY Annual fee payment

Payment date: 20190701

Year of fee payment: 11

PR1001 Payment of annual fee

St.27 status event code: A-4-4-U10-U11-oth-PR1001

Fee payment year number: 11

PR1001 Payment of annual fee

St.27 status event code: A-4-4-U10-U11-oth-PR1001

Fee payment year number: 12

PC1903 Unpaid annual fee

St.27 status event code: A-4-4-U10-U13-oth-PC1903

Not in force date: 20210711

Payment event data comment text: Termination Category : DEFAULT_OF_REGISTRATION_FEE

PC1903 Unpaid annual fee

St.27 status event code: N-4-6-H10-H13-oth-PC1903

Ip right cessation event data comment text: Termination Category : DEFAULT_OF_REGISTRATION_FEE

Not in force date: 20210711