KR20010101040A - 여과된 히스토그램을 이용한 디지털 엑스레이 검출기를위한 인터레스트 영역의 자동검출 방법 - Google Patents

여과된 히스토그램을 이용한 디지털 엑스레이 검출기를위한 인터레스트 영역의 자동검출 방법 Download PDF

Info

Publication number
KR20010101040A
KR20010101040A KR1020017006365A KR20017006365A KR20010101040A KR 20010101040 A KR20010101040 A KR 20010101040A KR 1020017006365 A KR1020017006365 A KR 1020017006365A KR 20017006365 A KR20017006365 A KR 20017006365A KR 20010101040 A KR20010101040 A KR 20010101040A
Authority
KR
South Korea
Prior art keywords
histogram
digital values
filtered
digital
integer
Prior art date
Application number
KR1020017006365A
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
Application filed by 추후제출, 디렉트 라디오그래피 코포레이션 filed Critical 추후제출
Publication of KR20010101040A publication Critical patent/KR20010101040A/ko

Links

Classifications

    • H—ELECTRICITY
    • H05—ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
    • H05G—X-RAY TECHNIQUE
    • H05G1/00—X-ray apparatus involving X-ray tubes; Circuits therefor
    • H05G1/02—Constructional details
    • 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/10—Segmentation; Edge detection
    • G06T7/11—Region-based segmentation
    • 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/20—Special algorithmic details
    • G06T2207/20021—Dividing image into blocks, subimages or windows
    • 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/20—Special algorithmic details
    • G06T2207/20112—Image segmentation details
    • G06T2207/20132—Image cropping
    • 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/30061—Lung

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Quality & Reliability (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)
  • Solid State Image Pick-Up Elements (AREA)
  • Facsimile Image Signal Circuits (AREA)

Abstract

본 발명은 대상물의 라디오그래프를 디스플레이하는데 사용하기 위해 대상물의 라디오그래픽 노출을 표시하기 위한 다수의 화소에 대응하는 제1세트의 디지털값으로부터 여과된 디지털값의 서브세트를 자동으로 형성하는 방법에 있어서, 제1세트의 디지털값내의 제1히스토그램을 발생시키는 단계와, 상기 제1히스토그램의 정수를 표시하는 정수곡선을 발생시키는 단계와, 제1히스토그램의 정수상에서 변곡점을 확인하고 그 대응의 디지털값(Vpeak1)을 확인하는 단계와, 상기 정수곡선은 가장 큰 디지털값으로부터 감소되는 방향으로 경사도를 먼저 변화시키며, Vtest = Vpeak1 - AR을 연산하는 단계와, 제1세트의 디지털값을 다수의 블럭 서브세트의 디지털값으로 공간적으로 분할하는 단계와, 각각의 블럭 서브세트내에서 디지털값의 발생주파수를 표시하는 블럭 서브세트 히스토그램을 구성하는 단계와, 각각의 블럭 서브세트상에서 블럭 서브세트 히스토그램의 설정의 낮은 포인트에 대응하는 낮은 포인트(BVlow)를 결정하는 단계와, 각각의 BVlow과 Vtest를 비교하는 단계와, Vtest 보다 큰 BVlow 를 갖는 모든 블럭 서브세트의 디지털값으로부터 여과된 디지털값의 서브세트를 형성하는 단계를 포함하며, 상기 정수곡선은 가장 높은 디지털값의 포인트로부터 시작하여 디지털값이 감소되는 방향으로 횡단할 때 먼저 경사도를 떨어뜨리며, 상기 AR은 대상물의 해부학적 범위를 나타내며, 상기 각각의 블럭 서브세트는 다수의 설정된 화소로 구성된 인접 화소로부터 디지털값을 포함하는 방법을 제공한다.

Description

여과된 히스토그램을 이용한 디지털 엑스레이 검출기를 위한 인터레스트 영역의 자동검출 방법{METHOD FOR AUTOMATIC DETECTION OF REGION OF INTEREST FOR DIGITAL X-RAY DETECTORS USING A FILTERED HISTOGRAM}
디지털 의료영상 분야에 있어서, 특히 디렉트 라디오그래피에 있어서, 입사방사선의 일부를 포획하여 이를 전자전하로 변환하는 방사선 검출기는 종종 영상획득을 위해 사용된다. 이러한 장치의 작동범위는 X레이 필름을 위한 40:1의 범위와 비교할 때 10,000:1 이상으로 선형으로 이루어질 수 있다. 검출가능한 노출의 넓은 범위때문에, 노출소자의 부적절한 선택으로 인한 환자를 재영상할 필요성이 없게 된다. 그러나, 디스플레이 장치로부터 획득 영상의 디커플링(decoupling)으로 인해, 때로 룩업테이블(look-up table: LUT)의 형상을 취하는 맵핑 함수(mapping function)를 전개시킬 필요가 있으며, 이것은 필요로 하는 출력 매체[CRT디스플레이, 필름, 반사프린트(reflection print)]상에서 영상을 최적으로 되게 한다.
최적의 그레이 스케일(gray scale)을 전개함에 있어서, 최종 영상의 허용가능성에 일부 주요한 요소가 관련되어 있음을 알 수 있다. 먼저, 검은 배경(black ground)으로서 약화되지 않은 X레이 노출영역을 제공하는 것이 바람직하다. 둘째로, 특정한 광학밀도로 조직의 다른 부분을 디스플레이할 필요가 있는데, 예를 들면 여러 실험자의 피부선(skin-line)을 인지할 수 있는 것이 바람직하다. 이러한 요구사항에 성공적으로 부응하기 위해, 배경 화소의 세기 범위를 인식하는 것처럼, 디스플레이할 필요한 영상특성의 전부 또는 일부를 표시하는 디지털값의 범위를 인식할 필요가 있다. 이러한 작업이 실행되었다면, 그레이 스케일의 특정한 실험자는 두 기준에 부응하는 필요한 조정을 자동으로 수행할 수 있다.
배경 검출의 문제는 의료용 영상분할(image segmentation)의 영역에 속한다. 이에 대해서는 2개의 일반적인 방법 즉, 1)히스토그램 분석과, 2)공간 조직분석법(히스토그램에 의해 제공된 그레이 레벨 세기정보와 조합될 수 있는)이 있다. 제1등급의 방법으로서, 아제웰 등에 허여된 미국특허 5.046.118호에는 히스토그램을 배경영역과 비배경영역으로 분할하기 위한 부분적 엔트로피라는 개념을 사용하는 방법이 서술되어 있다. 카포지 등에 허여된 미국특허 5.164.993호에는 후자의 방법이 개시되어 있으며; 이에 따르면, 배경포인트가 정점의 상부에서 발견되었을 때의 상태를 위한 일부의 부가물을 제공하므로써 이를 선형 및 로가리즘 히스토그램상에서 이용하고 있다. 다나카 등에게 1997년 1월 21일자로 허여된 미국특허 5.596.654호에는 자동임계장치에 의해 다수의 섹션으로 분할된 히스토그램을 이용하여 배경과 전경(콜리메이터 블레이드와 같은 방사선 제한장치의 사용에 따라 매우 소량의 방사선을 수용하는 영상영역)을 찾아내는 방법이 개시되어 있다. 미국특허 5.596.654호에서는, 실험자 타입과 노출기법과 디스플레이될 필요로 하는 몸체 등에 대한 정보와 결합된 판별분석(discrimanant analysis)이 사용되어; 전경과 물체와 배경영역에 대한 필요로 하는 범위가 발견될 때까지 상기 섹션들 사이의 분리포인트를 조정하는 방법을 사용하고 있다.
두번째 방법으로서, 장(Jang) 등에게 허여된 미국특허 5.268.967호에는 형태학적 엣지 검출단계와, 블럭 분류단계와, 블럭 정제단계 등을 포함하는 4단계방법이 개시되어 있다.
상기 미국특허 5.046.118호와 5.164.993호 및 5.596.654호에 개시된 방법은 배경이 불균일한 방법으로 변화되거나 또는 복합적인 배경 정점이 영상의 히스토그램으로 존재할 때는 유용하지 않다. 미국특허 5.268.967호에 개시된 방법은 상술한 바와 같은 용도에서는 유용하지만, 시간이 소요되는 복합적 처리스테이지가 포함된다. 조직분석방법은 소프트웨어로 실행시와 양호한 신뢰성에서는 매우 느리다. 특히 시준된 노출이 불필요한 영상을 생성할 때 배경이 분산되는 경우, 방사선 노출에 응답하여 전기 신호를 생성하는 방사선 검출기에 의해 포획된 라디오그래피의 정확한 그레이 스케일을 자동으로 지속적으로 재생하기 위한 실질적인 방법이 아직도 요망되고 있다.
본 발명은 X레이 세기 분포 라디오그래프를 표시하는 디지털값의 범위를 자동으로 확인하기 위한 방법에 관한 것으로서, 특히 유용한 진단 영상정보를 표시하는 디지털값의 범위를 확인하기 위한 방법에 관한 것이다.
도1은 본 발명에 따라 여과를 실시하기 전에 시준된 X레이에 노출된 손에 대한 전형적인 히스토그램과 그 정수를 도시한 도면.
도2는 필요한 실험 특정 알고리즘(exam specific algorithm)을 인가한 후 도1의 히스토그램에 의해 표시된 데이터를 디스플레이하므로써 생성된 영상을 도시한 도면.
도3은 해부학적 범위 상수(anatomical range constant)의 이탈을 도시하는 예시적인 히스토그램.
도4는 도1의 히스토그램을 구성하는데 사용된 동일 데이터로부터의 전형적인 블럭 히스토그램을 도시한 도면.
도5는 본 발명에 따라 데이터 처리를 한 후 도1의 히스토그램에 대응하는 최종적인 여과된 히스토그램을 도시한 도면.
도6은 도2의 경우에서처럼 필요로 하는 상기 실험 특정 알고리즘을 적용시킨 후 도5의 여과된 히스토그램에 의해 표시된 데이터를 디스플레이하여 최종적으로 생성된 영상을 도시한 도면.
도7은 본 발명의 예시적인 방법의 흐름도를 도시한 도면.
따라서, 본 발명의 목적은 실험-특정 처리(exem-specific processing)를 하기 전에 진단 디스플레이에 사용된 인터레스트의 해부학적 영역을 표시하는 유용한 디지털값의 범위를 자동으로 확인하는 방법을 제공하는 것이다.
본 발명은 대상물의 라디오그래프를 디스플레이하는데 사용하기 위해 대상물의 라디오그래픽 노출을 표시하기 위한 다수의 화소에 대응하는 제1세트의 디지털값으로부터 여과된 디지털값의 서브세트를 자동으로 형성하는 방법에 있어서, 제1세트의 디지털값내의 제1히스토그램을 발생시키는 단계와, 상기 제1히스토그램의 정수를 표시하는 정수곡선(integral curve)을 발생시키는 단계와, 제1히스토그램의 정수상에서 변곡점을 확인하고, 상기 변곡점에 대한 디지털값(Vpeak1)을 얻는 단계와, Vtest = Vpeak1 - AR을 연산하는 단계와, 제1세트의 디지털값을 다수의 블럭 서브세트의 디지털값으로 공간적으로 분할하는 단계와, 각각의 블럭 서브세트내에서 디지털값의 발생주파수를 표시하는 블럭 서브세트 히스토그램을 구성하는 단계와, 각각의 블럭 서브세트상에서 블럭 서브세트 히스토그램의 설정의 낮은 포인트에 대응하는 낮은 포인트(BVlow)를 결정하는 단계와, 각각의 BVlow과 Vtest를 비교하는 단계와, Vtest 보다 큰 BVlow 를 갖는 모든 블럭 서브세트의 디지털값으로부터 여과된 디지털값의 서브세트를 형성하는 단계를 포함하며, 상기 정수곡선은 가장 높은 디지털값의 포인트로부터 시작하여 디지털값이 감소되는 방향으로 횡단할 때 먼저 경사도를 떨어뜨리며, 상기 AR은 대상물의 해부학적 범위를 나타내며, 상기 각각의 블럭 서브세트는 다수의 설정된 화소로 구성된 인접 화소로부터 디지털값을 포함하는 방법을 제공한다.
상기 해부학적 범위(AR)는 설정된 상수이거나 변수 일수도 있으며, 상기 AR은 AR = [100+{(Vpeak-Vedge)-ARCref)}] 이며; ARCref는 실험형태의 함수로서 Vpeak 와 Vedge 사이의 표준편차를 나타내는 실험적으로 유도된 상수이며; Vedge는노출이 증가되는 방향으로 Vlow로부터 시작되는 히스토그램을 따라 횡단할 때, 히스토그램 발생주파수가 이전의 정점 발생주파수의 설정된 퍼센트로 먼저 낙하되는 포인트를 의미한다.
디지털값의 여과된 서브세트는 디스플레이 매체상에 다음과 같은 동작에 의해 라디오그래피를 나타내는 영상을 디스플레이하기 위해 사용된다.
Ⅰ) 디스플레이 매체를 위한 밀도 이송 함수를 표시하는 제1표준화 룩업테이블(look-up tabel: LUT)을 메모리에 저장하는 단계.
Ⅱ) 디스플레이될 형태의 실험적으로 결정된 다수의 상수 특성을 함유하는 제2LUT를 메모리에 저장하는 단계.
Ⅲ) 여과된 히스토그램을 형성하는 단계와, 그 정수를 결정하는 단계와, 상기 정수의 곡선을 생성하는 단계와, 상기 여과된 디지털값의 서브세트로부터 디지털값을 이용하여 여과된 히스토그램 정수곡선상에서 여과된 히스토그램 정수의 설정의 낮은 퍼센트에 대응하는 낮은 포인트(Vlow)를 결정하는 단계와; 여과될 히스토그램상에서 노출이 증가되는 방향으로 Vlow로부터 시작되는 히스토그램을 따라 횡단할 때, 히스토그램 발생주파수가 이전의 정점 발생주파수의 설정된 퍼센트로 먼저 낙하되는 포인트(Vedge)를 확인하는 단계와; 여과된 히스토그램 정수곡선상에서 노출이 감소되는 방향으로 가장 높은 노출 포인트로부터 시작되는 여과된 히스토그램 정수곡선을 따라 횡단할 때, 여과된 히스토그램 정수곡선이 경사도를 먼저 변화시키는 포인트(Vpeak2)를 확인하는 단계와; 상기 경험적으로 결정된 상수(constant)와 Vlow, Vpeak2 및 Vedge를 하나이상 포함하는 특정실험식(exam-specific equation)을 사용하여 시작값(Vmin) 및 종료값(Vmax)을 연산하는 단계.
Ⅳ) Vmin 보다 낮은 여과된 서브시트에서의 모든 디지털값을 Vmin으로 교체하고, Vmax 보다 큰 모든 디지털값을 Vmax로 교체하는 단계.
Ⅴ) 표준화된 디지털값과 영상을 디스플레이할 제1LUT를 이용하여 여과된 서브시트에서의 디지털값을 Vmax - Vmin의 편차를 표시하는 새로운 세트값으로 표준화하는 단계.
상기 방법의 단계는 프로그램 저장장치에서 디지털 라디오그래피 시스템으로 구성되고 상기 방법의 단계를 실행하는 장치에 의해 판독가능하고 실행가능한 명령프로그램으로 사용될 수 있다.
전형적인 직접 영상 포획 셋업에서는 X레이 방사선과, 환자, 및 방사선 검출기일 수도 있는 영상 방사선소스가 있다. 환자는 방사선의 통로에 위치되며, 검출기는 환자를 횡단한 후 방사선을 차단하는 위치에서 방사선 통로에 위치된다.
직접적인 방사선 영상 포획장치는 전형적으로 2차원 어레이로 배치된 다수의 센서를 포함한다. 상기 센서는 전기신호를 발생시키는데 이러한 신호의 크기는 센서에 입사된 방사선의 지속기간과 그 세기에 직접적으로 비례한다.
각각의 화소(PIXEL)로부터 신호전체는 서술된 시컨스로 판독되며; 증폭과 노이즈 필터링에 이어, 이러한 포인트에서 전형적으로 아날로그 신호로 형성된 신호는 각각의 센서의 상대적 노출을 표시하는 디지털값으로 항상 전환된다. 디지털값은 데이터뱅크(data bank)로 언급되는 데이터 저장매체에 저장된다. 저장은 어레이에 포획된 신호분배의 정확한 2차원 재구성을 허용하는 방식으로 이루어진다.
데이터뱅크내의 정보는 본래의 신호분배의 시각적 영상을 2차원 분배의 그레이 스케일값의 형태로 디스플레이하는데 사용되며, 이들 각각의 값은 어레이의 각각의 센서에 의해 생성된 신호에 대응한다. 이러한 각각의 디지털값은 디스플레이된 영상을 구성하는 화소를 형성한다. 이러한 전형적인 방사선 검출기는 본원에 참조인용되고, 1997년 7월 5일자로 리(Lee) 등에 허여된 미국특허 5.648.660호에 개시되어 있다.
영상을 생성하기 위해서는 특정의 디스플레이 매체가 선택되어야만 한다. 대부분의 통상적인 디스플레이 매체는 감광성 필름과 음극선튜브이다. 감광성 필름은 영상방향으로 세기변조된 스캐닝 광비임에 노출될 수도 있으며, 포커스된 전자의 형광체 자극 스캐닝비임의 세기를 영상방향 형태로 변조시키므로써 CRT디스플레이가 실행될 수도 있다. 이러한 디스플레이 형태를 위한 기법은 본 기술분야에 널리 공지되어 있다.
영상을 디스플레이할 때, 디지털값을 표시하는 신호를 수신하고 디스플레이 구동신호를 생성할 수 있는 장치에서 전자비임이나 광비임을 세기를 조절하기 위해 디지털값이 사용된다. 현존의 기법은 디지털 입력값을 발광막(luminous film)의 비임 세기(또는 노출시간)을 변조시키는 변조전압이나, 가시 영상을 생성하기 위해 CRT형광체상에 충돌하는 전자비임 또는 노출비임으로 전환하는데 널리 사용된다. 중요한 것은 최종적인 영상이 필요로 하는 광학 밀도 외관을 가질 수 있도록 각각의 신호강도로 디지털값의 전환이다. 이러한 처리는 사용된 디스플레이 매체의 그레이 스케일 이송(GST) 함수상에 포획된 영상을 제공하는 디지털값의 맵핑으로 언급된다. 상기 GST함수는 디스플레이 매체의 최종 그레이 스케일을 표준화된 그리고 전형적으로 LUT로 저장된 입력 자극(excitation)의 함수로 제공한다.
의료용 진단영상 디스플레이 용도에서는 전통적인 라디오그래피의 외관을 가능한한 정확하게 재생하는 것이 바람직하다. 전형적으로, 전통적인 HeNe 레이저 필름 라디오그래피는 약 3.0의 최대 광학밀도와, 0.2의 포그 레벨(fog level)을 갖고 있다. 따라서, 유용한 밀도범위는 약 2.8이다. 필름의 광학밀도는 전형적으로 로가리즘으로 표시되며, 약 2.8의 역전 로그는 약 600이므로, 필요한 범위의 사용한 그레이 레벨은 600 내지 1 또는 최대 1000 내지 1 이다. 이것은 디스플레이 그레이 스케일 범위가 600/1 또는 1000/1 의 광학밀도 범위를 재생해야 한다는 것을 의미한다. 또는 다른 의미로서, 디스플레이에서 그레이의 제1검출가능 세이드(shade)를 제공하는데 사용될 디지털값(Vmin)과, 디스플레이에서 최대 블랙(maximum black) 또는 포화포인트를 디스플레이하는데 사용될 디지털값(Vmax)이 결정되어야 한다는 것을 의미한다. 또한, Vmin과 Vmax 사이의 나머지 디지털값은 수행될 실험의 형태에 기초하여 그레이 범위를 압축하거나 팽창하면서, 필요로 하는 그레이 범위에서 가장 유용한 시각적 정보를 제공하는 방식으로 펼쳐져야 한다. 이것은 제3디지털값(Vedge)을 이용하여 얻어지는데, 이러한 디지털값은 GST상의 설정포인트에서 특정한 밀도의 값을 대체하는데 사용된다.
본 발명의 양수인에게 양도되고 1998년 2월 20일자로 출원된 계류중인 PCT/US98/03249호에는 방법과, 관련의 실험결정된 실험-특정 알고리즘과, 이러한 알고리즘에 사용하기 위한 관련의 정수가 개시되어 있다. 이러한 방법은 Vmin, Vmax, Vedge를 결정하고 라디오그래피를 디스플레이하기 위해 디지털 데이터에 적용될 수 있다.
상술한 계류중인 특허출원서에 개시된 실험-특정 알고리즘은 임계값과 그 관관계를 확인하기 위해 해부학적 랜드마크로 히스토그램 데이터의 과도한 보정을 이용하고 있으며; 또한, 라디오그래프의 경우 필요로 하는 실험의 형태와 포획된 영상의 특성에 의해 결정되는 방식으로 디지털값을 맵핑하기 위해 이러한 값을 사용한다. 이러한 해부학적 랜드마크는 모든 형태의 영상과 실험에 공통적인 3개의 기준포인트를 사용하여 실험-특정 연산에 의해 결정된다. 도1에는 전형적인 히스토그램(10)과, 여과하기 전에 상기 3개의 기준포인트를 도시하는데 사용되는 손의 시준된 X레이 노출을 위한 정수(12)가 도시되어 있다.
제1기준포인트는 영상 디스플레이 맵핑에 사용될 최소값(Vmin)을 연산하기 위해 사용되는 시작되는 기준포인트(Vlow)이다. 시작 포인트의 확인은 만일 히스토그램 발생주파수값만 사용된다면 불명확해질 수 있다. 히스토그램 값의 정수를 사용하면 보다 명확한 값을 생성할 수 있다. 1% 내지 10%, 양호하기로는 전체 정수의 5%에 대응하는 디지털값(Vlow)이, 디스플레이를 위한 초기 디지털값(Vmin)을 연산하는데 신뢰성있는 기준치로 선택된다.
제2기존포인트는 가장 낮은 노출로부터 가장 높은 노출까지(도1에서 좌측으로부터 우측까지) Vlow로부터 횡단하면서, 히스토그램 발생주파수가 히스토그램상에서 이전 정점값의 20% 내지 30%, 양호하기로는 25% 이하로 떨어지는 Vedge 이다. 예를 들어, PA 가슴영상에서, Vedge는 폐에서 X레이의 투과성이 가장 높은 부분인 것을 확인하며, 하드 카피영상에서 특정화된 광학 밀도를 생성하기 위해 사용된다. 도시된 바와 같이, 폐에 대한 설정의 광학 밀도와 Vmin의 조합은 Vmax값의 연산을 허용한다. 가슴을 측방향으로부터 투시하는 경우, Vedge는 심장영역을 나타내며,측방향 가슴투시에서 Vmax의 연산에 사용된다.
제3기준포인트는 노출값에 대응하는 Vpeak로서, 가장 높은 노출로부터 가장 낮은 노출까지(도1에서 우측으로부터 좌측까지) 횡단하면서, 이러한 포인트에서 히스토그램 정수곡선은 경사도를 변화시킨다. Vpeak는 히스토그램의 가장 높은 불포화된 정점값을 나타낸다. 이러한 값은 피부선을 보고자 하는 영상에서 유용하게 사용된다.
히스토그램 데이터로부터 유도된 기준값(Vlow)을 사용하였을 때, Vmin은 실험-특정 정수(A)를 사용하여 연산된다. 상기 정수(A)는 영상 디스플레이에 포함된 가장 낮은 디지털값에 대해 Vlow 하방으로의 거리를 결정하는데 사용된다. 상기 정수(A)와 그 관련의 쌍으로 이루어진 정수(B)는 실험적으로 결정된 정수로서, 이들은 다수의 검사자로부터 그 어떠한 제한없이 가슴, 손발, 두개골, 복부 등을 포함하는 다른 형태의 라디오그래픽 영상으로부터 유도된다.
Vmin = Vlow - SF ×Log(A)
스케일 팩터(SF)는 히스토그램 구성과 관련이 있다. 히스토그램을 구성하는 양호한 모드에서는 각각의 디지털 노출값을 사용하는 것이 아니라, 다수의 연속값을 포함하는 비인(bin)이 사용된다. 예를 들어, 각각의 비인이 20개의 연속한 디지털값을 포함하고 있다면, 20이 비인의 크기가 된다. 따라서, 디지털 시스템이 12비트 시스템인 경우에는 4096개의 분리된 값이 있으며; 이들은 범위를 20개의 연속한 값의 세그먼트(비인)로 분할하므로써, 또한 비인내에서의 그 어떤 값의 발생 주파수를 단일값으로 하므로써, 205개로 감축될 수 있다.
상기 스케일 팩터는 디스플레이 매체 작동범위로 맞추어진 센서 작동범위로부터 생데이터(raw data)의 로가리즘 전환과도 관련이 있다. 예를 들어, HeNe 레이저 필름 시스템에서, 필름 노출의 유용한 광학 밀보 작동범위는 베이스 플러스 포그를 고려했을 때 600 내지 1 이다. 스케일 팩터(SF)는 다음과 같이 한정된다.
SF = [(유용한 값의 갯수) / (작동범위의 로그)]/ 비인크기
이것은 예를 들어 다음과 같이 주어진다.
SF = [4096 / log(600)]/20 또는 SF = 73.7
일단 Vmin값이 결정되었다면, Vmax는 선택된 노출의 형태에 따라 하기의 관련된 것들중 하나를 사용하여 연산된다.
경우Ⅰ: 일반적인 라디오그래픽 영상(디폴트값)
디지털값의 설정된 실험 특정 작동범위는 디스플레이 데이터에 포함되는 것으로 서술된다. 이러한 모드는 그 감광계 특성의 필요한 작동범위를 갖는 필름을 사용하는 것과 동일한 것이다. Vmin값의 히스토그램 결정은 노출변화를 보상하며, 그레이 스케일 이송특성은 Vmin에 대해 고정될 것이다. Vmax는 다음과 같이 연산된다.
(1) Vmax = Vmin + SF ×Log(B)
B는 A와 쌍으로 이루어진 상수이다.
경우Ⅱ: (대상물 콘트라스트를 위한 보상)
경우Ⅱ에 있어서, Vedge와 Vlow 사이의 편차는 디스플레이에 사용될 디지털 데이터의 작동범위를 결정하는데 사용된다. 경우Ⅰ에서 실험-특정 작동범위변수는 이러한 값에서 실험-특정 표준 편차에 대한 관찰된 편차(Vedge - Vlow)의 비율로 수정된다. 따라서, Vmax는 다음과 같이 연산된다.
(2) Vmax = Vmin + SF ×Log(B) - △범위
△범위 = Vedge - Vlow
여기서, Vedge는 Vmin에서 시작해서 좌측에서 우측으로 스캔되었을 때 이전 가장 높은 정점의 20-30%, 양호하기로는 25%로 떨어지는 포인트이다.
경우Ⅲ: (해부학적 구조와 연관된 디스플레이값)
경우Ⅲ에 대해, Vedge값은 해부학적 구조와 연관되어 있으며, 이러한 구조는 진단 디스플레이에서 상술의 광학밀도(또는 소프트 디스플레이 세기)로 디스플레이될 것이다. 이를 달성하기 위해, 선택된 그레이 스케일 이송곡선이 참고되어야 한다.
예를 들어, 만일 Vedge가 폐에서 최대 투과부에 대응한다면, 상기 Vedge값은 1.7의 광학 밀도(optical density: OD)로 할당되어야 하는 것으로 서술되어 있다. 그레이 스케일 이송곡선으로부터, 주어진 OD = 1.7을 얻기 위해서는 55%로 결정되어야 한다. GST함수는 Vedge 와 Vmin 사이의 거리가 Vmax와 Vmin 사이의 디지털값 범위의 55%에 대응하도록 표준화된다. 편리함을 위해, 이러한 실험-특정 노출 비율변수는 "L"로 표시되었다.
하기의 식은 Vedge가 하드 카피 영상에서 필요로 하는 광학 밀도를 생성할 수 있도록 Vmax를 생성한다.
(3) Vmax = Vmin +(1/L) ×(Vedge - Vmin)
경우Ⅳ: (인터레스트의 영역이 손발인 경우)
예를 들어 손발의 라디오그래프인 경우, 피부선의 가시성이 필요할 때, Vmax 는 다음과 같이 연산된다.
(4) Vmax = Vmin + (1/L) ×(Vpeak - Vmin)
이 경우 L은 그레이 스케일 이송곡선으로 언급되는 실험-특정 변수이며, 피부선에서 요구되는 광학 밀도를 생성하는데 필요한 퍼센트 노출을 나타낸다. 이런 상태에서, 해부학적 영상을 둘러싸고 있는 배경은 피부선을 위해 선택된 동일한 광학밀도에 있게 된다.
상술한 바는 Vmin을 결정하는 방법과, Vmax를 연산하기 위한 다양한 방법에 관한 것이다. 이러한 값들은 진단학적 디스플레이 영상에서 표시될 디지털값의 범위를 나타낸다. 이러한 범위내의 디지털값은 프로그램된 룩업테이블(LUT)을 통해 전송되어 필요로 하는 그레이 스케일 이송특성을 제공하므로써, 진단 데이터의 디스플레이를 최적화시킨다.
이를 달성하기 위해, 본 발명은 변수 테이블에 저장되고 영상과 실험 형태작업자 선택에 의해 인덱스된 다수의 실험 특정 변수를 사용하였다.
이러한 변수는 다음과 같다.
A: Vmin을 결정하기 위한 것.
B: Vmax를 결정하기 위한 것.
L: 필요로 하는 광학 밀도를 생성하기 위한 것.
또한 실험 및 영상 형태에 의존하여, Vmax를 연산하기 위한 4개의 알고리즘이 작업자에 의해 저장되어 선택적으로 사용된다. 마지막으로, 특정의 디스플레이 매체 그레이 이송함수는 출력밀도를 노출의 퍼센트로 도시하는 표준화된 포맷으로도 저장된다. 통상적인 경우에서처럼 하나의 디스플레이가 사용되는 것보다, 하나의 GST 함수 곡선이 저장된다.
상술한 특허출원서에 서술된 바와 같이, 전형적인 시스템은 상이한 실험자에 대한 영상을 디스플레이하기 위해 작업자에 의해 선택되는 다음과 같은 변수를 포함한다.
실험/영상 형태 입력 1 가슴-PA/AP
2 가슴-측방향
3 손발
하기의 표1은 상술한 경우에 있어서 선형의 14비트 디코더 출력으로부터 로가리즘 12비트 도메인까지 전환 맵핑을 따르는 로그 디지털값의 발생주파수의 히스토그램을 이용하는 예시적인 정수값을 제공한다.
표1
실험형태 A B 등식 번호 L GST
1 5 100 3 0.57 F(ilm)
2 5 100 3 0.47 F
3 2 70 2 0.55 F
4 3 100 4 0.55 F
표1및 표2에 사용된 "등식 번호"는 본 발명에서 번호가 매겨진 등식을 의미한다. 표2는 실험-특징 등식 계수 및 변수를 위한 예시적인 범위 및 디폴트값을 제공한다.
표2
변수 MIN MIX DEFAULT
A 0.1 10.0 2.0
B 0/0 100 100
L 0.0 1.55 0.5
콘트라스트 0.5 1.50 1.0
밀도 -0.25 1.55 0.5
범위 0 200 53
해부학적 범위 20 200 100
등식 번호 1 4 1
도2에는 도1의 히스토그램에의해 제공된 데이터를 디스플레이하여 형성된 영상(20)이 도시되어 있다. 히스토그램(10)은 바람직하지 않은 데이터를 포함하는 라디오그래픽 영상(20)을 나타내는 본래의 디지털값(로그 변형에 이어 비인 도메인당-12비트 20-값으로 맵핑되는)의 발생주파수를 도시하고 있다. 이러한 바람직하지 않은 데이터는, 14×17인치 방사선 검출기의 일부에만 노출하기 위해 시준된 방사선 비임이 사용될 때, 차폐는 검출기의 완전노출을 방지하지만, 분산된 방사선이 도2에 도시된 바와 같이 바람직하지 않은 배경 노출(22)을 생성하는 경우에서, 예를 들면 실제 노출영역의 외부로부터의 데이터를 나타낸다.
이러한 바람직하지 않은 노출을 나타내는 데이터가 히스토그램의 구성과 Vmin, Vmax, Vedge의 연산에 포함되었을 때, 그 결과는 바람직하지 않게 된다. 도2에 도시된 바와 같은 영상(20)은 품질이 조악하다. 분산된 방사선 영상에 의해, 유용한 많은 정보가 불명확하게 되었으며, 그 부분은 필요로 하는 영상영역(24)을 둘러싸고 있는 그레이 영역(22)으로 도시되어 있다.
본 발명에 따르면, 이러한 본래의 데이터는 상술한 3개의 값(Vlow, Vpeak, Vedge)을 추출하기 위해 사용되었다. 또한, 본 발명은 제4의 기준값인 Vtest를 사용하는데, 이러한 값은 배경분산이 아닌 예견된 최소한의 데이터값이다. 해부학적 데이터에 대응하는 디지털값은 한정된 해부학적 범위(AR)을 갖는데, 이러한 범위는 일반적으로 약 100비인(각각의 비인이 20개의 디지털값을 포함하는 상술한 히스토그램 구성에 연속하는)이다. 따라서, 배경분산이 아닌 최소한의 데이터값은 다음과 같이 연산될 수 있다.
Vtest = Vpeak - AR (여기서, AR은 100이다)
그러나, 일부 실험자는 100:1 보다 크거나 작은 비정상치(AR)를 가질 수 있다. AR에서의 이러한 편차는 각각의 실험자 형태에 있어서 정상 대 비정상 실험자를 위한 범위(Vpeak - Vedge)의 값에서의 편차로 표시된다. 이러한 범위는 도3에 도시된 바와 같이 ARC(anatomic range compensation: 해부학적 범위 보상)로 언급될 수 있다. 각각의 실험자는 ARCref값(기준 ARC)을 가지며, 이러한 값은 수많은 임상학적 데이터베이스로부터의 히스토그램의 넓은 분포로부터 추정될 수 있으며, ARCact(실제 ARC)는 처리될 특정 히스토그램을 위한 실제값(Vpeak - Vedge)에 대응한다.
따라서, AR은 다음과 같이 표현된다.
AR = 고정범위 + (ARCact - ARCref)
ARCref = Vpeak - Vedge
따라서, 해부학적 범위(AR) 보상을 고려하면, 다음과 같다.
Vtest = Vpeak - [고정범위 + {(Vpeak - Vedge) - ARCref}]
이러한 처리를 도시하기 위해 사용된 양호한 실시예에서, 고정범위값은 실험적으로 결정되어 100이 된다. 히스토그램 스케일에 따라서, 상기 고정범위는 특정의 실험 형태에서 예견할 수 있는 전형적인 유용한 값의 범위(그리고 이에 따른 그레이 스케일 범위)인 실험적으로 결정된 값의 범위를 나타내기 위해 조정될 수도있다.
100의 고정범위에 대해, 이전의 식은 다음과 같이 변화된다.
Vtest = Vpeak - [100 + {(Vpeak - Vedge) - ARCref}]
본 발명에 따르면, 여과된 히스토그램 즉, 라디오그래픽을 디스플레이하는데 사용된 데이터값을 확인하기 위해 사용될 값의 히스토그램은 예비적인 완전영상 히스토그램과 그 대응의 정수곡선을 먼저 편집하므로써 구성된다. 본 발명은 정수곡선상에서 포인트(Vpeak)를 결정하고 Vtest를 얻었다. 이어서, 상기 완전영상은 2차원 정방향 또는 사각형 블럭의 픽셀의 형태를 취하는 다수의 서브영상으로 분할된다. 그후, 각각의 블럭에 대한 히스토그램과 블럭 히스토그램의 정수곡선이 구성된다. 블럭 정수의 1% 내지 10%, 양호하기로는 5%에 대응하는 기준포인트(VBlow)가 각각의 블럭을 위해 생성되어, Vtest와 비교된다. 그후, 여과된 완전영상 히스토그램은 BVlow이 Vtest 보다 큰 블럭만으로 모든 데이터의 총합을 구성한다. 이와 같이 생성된 여과된 히스토그램은 상술한 바와 같이 Vlow, Vpeak2, Vedge를 결정하기 위해 사용되며, 이러한 값들로부터 Vmin과 Vmax가 영상을 디스플레이하기 위해 결정된다. Vpeak1, 및 Vpeak2의 표현은 여과된 완전영상 히스토그램(Vpeak2)으로부터 예비 완전영상 히스토그램(Vpeak1)상에 결정된 Vpeak를 식별하기 위해 사용된다.
일단 여과된 히스토그램의 Vmin, Vedge, Vmax가 결정되었다면, 작업자의 입력이 디스플레이될 영상과 디스플레이될 유용한 값의 범위를 선택한다. 각각의 값이 디스플레이를 위해 이동됨에 따라, 이것은 노출값의 퍼센트로 전환된다. 이러한 퍼센트값은 필요로 하는 광학밀도를 생성하는데 필요한 최종 노출을 위한 값을 발생시키기 위해, 표준화된 노출의 함수로서 디스플레이 매체 밀도에 대응하는 저장된 룩업테이블(LUT)을 이용하여 전환된다. 그리고, 이러한 값은 만일 매체가 감광막일 경우 매체를 디스플레이하기 위해 사용되거나, 또는 소프트 디스플레이가 사용된 경우 CRT 디스플레이 장치를 구동시키기 위해 사용된다.
도4는 본래의 영상이 서브분할되는 다수의 블럭을 위한 디지털값의 발생주파수의 전형적인 히스토그램을 도시하고 있다. 이러한 각각의 서브 히스토그램은 본래의 히스토그램으로부터 Vtest에 대해 테스트되어 얻어진 대응의 Vlow를 갖는다.
도5는 BVlow가 Vtest 보다 큰 블럭 서브세트 히스토그램으로부터의 디지털값으로부터 구성하므로써 얻어진 여과된 히스토그램을 도시한다. 도5에 도시된 바와 같이, x축상에서 비인(85) 이하의 비인에 대응하는 값은 더 이상 포함되지 않는다. 도6은 도2에서 영상을 표시하는데 사요된 동일한 실험-특정 알고리즘(exam-specific algorithm: ESA)의 적용에 이어서, 도5의 여과된 히스토그램에 의해 제공된 데이터의 사용으로부터 파생된 손의 영상을 도시하고 있다.
하드웨어 또는 다양한 작동 스테이지에서 작업자의 입력에 의해 전환과 연산이 실행되지만, 이러한 처리의 연산부는 적절히 프로그램된 컴퓨터를 사용하면 최적으로 실행될 수 있다. 따라서, 상술한 방법의 단계를 실행하기 위해, 특히 이러한 장치가 디지털 라디오그래피 시스템의 일부분일 때, 본 발명은 장치에 의해 실행가능한 명령 프로그램을 실행하며 장치에 의해 판독가능한 프로그램 저장장치를포함한다.
도7에는 이러한 프로그램의 컴퓨터 로직에 사용될 수 있는 예시적인 흐름도가 도시되어 있다. 도7에 도시된 바와 같이, 처리는 본래의 디지털 라디오그래피 데이터를 수집한 단계(100)에서부터 시작된다. 이러한 데이터의 히스토그램 정수에서의 히스토그램은 단계(110)에서 생성되며, Vlow, Vpeak/, Vedge 기준값이 결정된다. 이어서, 단계(130)에서는 실험-특정 알고리즘(ESA)(125)(Vtest)이 연산된다.
단계(140)에서, 본래의 디지털 데이터는 블럭으로 서브분할된다. 단계(150)에서, 특정의 블럭을 위한 데이터가 회수되며; 단계(160)에서는 히스토그램과 이러한 블럭을 위한 히스토그램 정수가 생성된다. 단계(170)에서는 블럭(Vlow)(BVlow)이 생성되며, 단계(180)에서는 Vtest와 비교된다. 만일 Vlow가 Vtest 보다 작거나 이와 같다면, 단계(190)는 데이터가 무시된다. 만일 Vlow가 Vtest 보다 크다면, 데이터는 단계(200)에서 데이터의 여과된 서브세트에 세이브/추가된다. 단계(180)의 결과에 따라, 단계(190) 또는 단계(200)를 통하는 단계(150)는 데이터의 블럭이 회수되어 평가된다.
일단 데이터의 완전 여과 서브세트가 생성된 후에는, 처리는 단계(210)로 진행되며; 이러한 단계에서, 여과된 히스토그램과 상기 데이터의 여과된 서브세트의 여과된 히스토그램 정수가 생성된다. 단계(220)에서는 Vlow, Vpeak2, Vedge가 생성되며, ESA(225)는 Vmin 및 Vmax를 연산하기 위해 사용된다. 디스플레이 데이터는 단계(240)에서 Vmin 이하로 여과된 서브세트에서의 값을 Vmin 으로 대체하고, Vmax 이상의 값을 Vmax로 대체하므로써 생성된다. 표준화된 디스플레이 데이터는 단계(270)에서 투시가능한 디스플레이를 생성하기 위해 단계(260)에서 그레이 스케일 이송(GST) 룩업테이블(LUT)을 통해 처리된다.
본 기술분야의 숙련자라면 특히 상기 주어진 실시예에서 사용된 실제 수치값에 대해 상술한 바를 참조하여 다양한 수정이 가능함을 인식해야 한다, 이러한 수정은 작업자 입력단계를 포함할 수 있으며; 이러한 단계에서 작업자는 필요한 ESA 를 선택하거나 또는 Vedge, Vlow, Vpeak 또는 환자의 특성이나 시행될 실험 형태와 관련있는 기타 다른 입력값을 한정한다. 그 어떠한 수정이라도 하기에 청구될 본 발명의 범주내에 포함되는 것을 인식해야 한다.

Claims (12)

  1. 대상물의 라디오그래프를 디스플레이하는데 사용하기 위해 대상물의 라디오그래픽 노출을 표시하기 위한 다수의 화소에 대응하는 제1세트의 디지털값으로부터 여과된 디지털값의 서브세트를 자동으로 형성하는 방법에 있어서,
    제1세트의 디지털값내에서 디지털값의 발생주파수를 표시하는 제1히스토그램을 발생시키는 단계와,
    상기 제1히스토그램의 정수를 표시하는 정수곡선을 발생시키는 단계와,
    제1히스토그램 정수곡선상에서 변곡점에 대응하는 디지털값(Vpeak1)을 구하는 단계와,
    Vtest = Vpeak1 - AR을 연산하는 단계와,
    제1세트의 디지털값을 다수의 블럭 서브세트의 디지털값으로 공간적으로 분할하는 단계와,
    각각의 블럭 서브세트내에서 디지털값의 발생주파수를 표시하는 블럭 서브세트 히스토그램을 구성하는 단계와,
    각각의 블럭 서브세트상에서 블럭 서브세트 히스토그램의 설정의 낮은 포인트에 대응하는 디지털값(BVlow)을 결정하는 단계와,
    각각의 BVlow과 Vtest를 비교하는 단계와,
    Vtest 보다 큰 BVlow 를 갖는 모든 블럭 서브세트의 디지털값으로부터 여과된 디지털값의 서브세트를 형성하는 단계를 포함하며,
    상기 정수곡선은 먼저 가장 큰 디지털값의 포인트로부터 시작되어 디지털값이 감소되는 방향으로 정수곡선을 따라 횡단하는 경사도를 변화시키며; 상기 AR은 대상물의 해부학적 범위를 나타내며; 상기 각각의 블럭 서브세트는 다수의 설정된 화소로 구성된 인접 화소로부터 디지털값을 포함하는 것을 특징으로 하는, 여과된 디지털값의 서브세트를 자동으로 형성하는 방법.
  2. 제1항에 있어서, 디지털값의 각각의 블럭 서브세트상에서 낮은 포인트를 결정하는 단계는 상기 각각의 블럭 서브세트 히스토그램에 대한 블럭 서브세트 히스토그램 정수곡선을 얻는 단계와, BVlow를 상기 정수곡선의 1% 내지 10%의 범위 사이로 설정하는 단계를 포함하는 것을 특징으로 하는, 여과된 디지털값의 서브세트를 자동으로 형성하는 방법.
  3. 제1항에 있어서, Vtest 연산단계에서 실험-특정 알고리즘이 디지털값의 세트에 적용되는 것을 특징으로 하는, 여과된 디지털값의 서브세트를 자동으로 형성하는 방법.
  4. 제1항에 있어서, Vtest 연산단계에서 다수의 실험-특정 알고리즘중에서 선택된 하나가 디지털값의 세트에 적용되는 것을 특징으로 하는, 여과된 디지털값의 서브세트를 자동으로 형성하는 방법.
  5. 제1항에 있어서,BVlow 연산단계에서 실험 특정 알고리즘이 디지털값의 서브세트에 적용되는 것을 특징으로 하는, 여과된 디지털값의 서브세트를 자동으로 형성하는 방법.
  6. 제1항에 있어서, Vlow로부터 시작하는 히스토그램을 따라 노출증가방향으로 횡단할 때 발생주파수가 발생주파수의 이전 정점주파수의 설정비율로 떨어지는 제1히스토그램상의 엣지포인트에 대응하는 디지털값(Vedge)을 얻는 단계를 부가로 포함하며, AR = [100 + {(Vmax - Vedge) - ARCref)}] 이고, 상기 ARCref는 실험적으로 유도된 정수인 것을 특징으로 하는, 여과된 디지털값의 서브세트를 자동으로 형성하는 방법.
  7. 제6항에 있어서, 상기 ARCref는 설정된 실험 형태에 대해 Vmax - Vedge의 표준값을 표시하며, 상기 표준값은 대응의 임상학적 데이터베이스로부터 대응하는 히스토그램으로부터 유도된 값인 것을 특징으로 하는, 여과된 디지털값의 서브세트를 자동으로 형성하는 방법.
  8. 제1항에 있어서, 다음과 같은 단계에 의해 라디오그래피를 표시하는 영상을 디스플레이 매체상에 디스플레이하는 단계를 포함하는 것을 특징으로 하는 여과된 디지털값의 서브세트를 자동으로 형성하는 방법.
    Ⅰ) 디스플레이 매체를 위한 밀도 이송 함수를 표시하는 제1표준화 룩업테이블(LUT)을 메모리에 저장하는 단계.
    Ⅱ) 디스플레이될 형태의 실험적으로 결정된 다수의 상수 특성을 함유하는 제2LUT를 메모리에 저장하는 단계.
    Ⅲ) 여과된 히스토그램을 형성하는 단계와, 그 정수를 결정하는 단계와, 상기 정수의 곡선을 생성하는 단계와, 상기 여과된 디지털값의 서브세트로부터 디지털값을 이용하여 여과된 히스토그램 정수곡선상에서 여과된 히스토그램 정수의 설정의 낮은 퍼센트에 대응하는 낮은 포인트(Vlow)를 결정하는 단계와; 여과될 히스토그램상에서 노출이 증가되는 방향으로 Vlow로부터 시작되는 히스토그램을 따라 횡단할 때, 히스토그램 발생주파수가 이전의 정점 발생주파수의 설정된 퍼센트로 먼저 낙하되는 포인트(Vedge)를 확인하는 단계와; 여과된 히스토그램 정수곡선상에서 노출이 감소되는 방향으로 가장 높은 노출 포인트로부터 시작되는 여과된 히스토그램 정수곡선을 따라 횡단할 때, 여과된 히스토그램 정수곡선이 경사도를 먼저 변화시키는 포인트(Vpeak2)를 확인하는 단계와; 상기 경험적으로 결정된 상수와 Vlow, Vpeak2 및 Vedge를 하나이상 포함하는 특정실험식(exam-specific equation)을 사용하여 시작값(Vmin) 및 종료값(Vmax)을 연산하는 단계.
    Ⅳ) Vmin 보다 낮은 여과된 서브시트에서의 모든 디지털값을 Vmin으로 교체하고, Vmax 보다 큰 모든 디지털값을 Vmax로 교체하는 단계.
    Ⅴ) 표준화된 디지털값과 영상을 디스플레이할 제1LUT을 이용하여 여과된 서브시트에서의 디지털값을 Vmax - Vmin의 편차를 표시하는 새로운 세트값으로 표준화하는 단계.
  9. 대상물의 라디오그래프를 디스플레이하는데 사용하기 위해 대상물의 라디오그래픽 노출을 표시하기 위한 다수의 화소에 대응하는 디지털값의 세트로부터 여과된 디지털값의 서브세트를 자동으로 형성하는 방법을 실행하는 장치에 의해 판독가능한 명령 프로그램을 사용하는 장치에 의해 판독될 수 있는 프로그램 저장장치에 있어서,
    제1세트의 디지털값내의 제1히스토그램을 발생시키는 단계와,
    상기 제1히스토그램의 정수를 표시하는 정수곡선을 발생시키는 단계와,
    제1히스토그램의 정수상에서 변곡점을 확인하고, 상기 변곡점에 대한 디지털값(Vpeak1)을 얻는 단계와,
    Vtest = Vpeak1 - AR을 연산하는 단계와,
    제1세트의 디지털값을 다수의 블럭 서브세트의 디지털값으로 공간적으로 분할하는 단계와,
    각각의 블럭 서브세트내에서 디지털값의 발생주파수를 표시하는 블럭 서브세트 히스토그램을 구성하는 단계와,
    각각의 블럭 서브세트상에서 블럭 서브세트 히스토그램의 설정의 낮은 포인0트에 대응하는 낮은 포인트(BVlow)를 결정하는 단계와,
    각각의 BVlow과 Vtest를 비교하는 단계와,
    Vtest 보다 큰 BVlow 를 갖는 모든 블럭 서브세트의 디지털값으로부터 여과된 디지털값의 서브세트를 형성하는 단계를 포함하며,
    상기 정수곡선은 가장 높은 디지털값의 포인트로부터 시작하여 디지털값이 감소되는 방향으로 횡단할 때 먼저 경사도를 떨어뜨리며, 상기 AR은 대상물의 해부학적 범위를 나타내며, 상기 각각의 블럭 서브세트는 다수의 설정된 화소로 구성된 인접 화소로부터 디지털값을 포함하는 것을 특징으로 하는 프로그램 저장장치.
  10. 제9항에 있어서, 다음과 같은 단계에 의해 라디오그래피를 나타내는 영상을 디스플레이 매체상에 디스플레이하는 단계를 부가로 포함하는 것을 특징으로 하는 프로그램 저장장치.
    Ⅰ) 디스플레이 매체를 위한 밀도 이송 함수를 표시하는 제1표준화 룩업테이블(LUT)을 메모리에 저장하는 단계.
    Ⅱ) 디스플레이될 형태의 실험적으로 결정된 다수의 상수 특성을 함유하는 제2LUT를 메모리에 저장하는 단계.
    Ⅲ) 여과된 히스토그램을 형성하는 단계와, 그 정수를 결정하는 단계와, 상기 정수의 곡선을 생성하는 단계와, 상기 여과된 디지털값의 서브세트로부터 디지털값을 이용하여 여과된 히스토그램 정수곡선상에서 여과된 히스토그램 정수의 설정의 낮은 퍼센트에 대응하는 낮은 포인트(Vlow)를 결정하는 단계와; 여과될 히스토그램상에서 노출이 증가되는 방향으로 Vlow로부터 시작되는 히스토그램을 따라 횡단할 때, 히스토그램 발생주파수가 이전의 정점 발생주파수의 설정된 퍼센트로 먼저 낙하되는 포인트(Vedge)를 확인하는 단계와; 여과된 히스토그램 정수곡선상에서 노출이 감소되는 방향으로 가장 높은 노출 포인트로부터 시작되는 여과된 히스토그램 정수곡선을 따라 횡단할 때, 여과된 히스토그램 정수곡선이 경사도를 먼저 변화시키는 포인트(Vpeak2)를 확인하는 단계와; 상기 경험적으로 결정된 상수와 Vlow, Vpeak2 및 Vedge를 하나이상 포함하는 특정실험식(exam-specific equation)을 사용하여 시작값(Vmin) 및 종료값(Vmax)을 연산하는 단계.
    Ⅳ) Vmin 보다 낮은 여과된 서브시트에서의 모든 디지털값을 Vmin으로 교체하고, Vmax 보다 큰 모든 디지털값을 Vmax로 교체하는 단계.
    Ⅴ) 표준화된 디지털값과 영상을 디스플레이할 제1LUT을 이용하여 여과된 서브시트에서의 디지털값을 Vmax - Vmin의 편차를 표시하는 새로운 세트값으로 표준화하는 단계.
  11. 제9항에 있어서, 상기 장치는 디지털 라디오그래피 시스템의 한 부품인 것을 특징으로 하는 프로그램 저장장치.
  12. 제10항에 있어서, 상기 장치는 디지털 라디오그래피 시스템의 한 부품인 것을 특징으로 하는 프로그램 저장장치.
KR1020017006365A 1998-11-19 1999-11-12 여과된 히스토그램을 이용한 디지털 엑스레이 검출기를위한 인터레스트 영역의 자동검출 방법 KR20010101040A (ko)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US09/196,391 1998-11-19
US09/196,391 US6118892A (en) 1998-11-19 1998-11-19 Method for automatic detection of region of interest for digital x-ray detectors using a filtered histogram
PCT/US1999/026703 WO2000032016A2 (en) 1998-11-19 1999-11-12 Method for automatic detection of region of interest for digital x-ray detectors using a filtered histogram

Publications (1)

Publication Number Publication Date
KR20010101040A true KR20010101040A (ko) 2001-11-14

Family

ID=22725218

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020017006365A KR20010101040A (ko) 1998-11-19 1999-11-12 여과된 히스토그램을 이용한 디지털 엑스레이 검출기를위한 인터레스트 영역의 자동검출 방법

Country Status (6)

Country Link
US (1) US6118892A (ko)
EP (1) EP1131776A4 (ko)
JP (1) JP2002530180A (ko)
KR (1) KR20010101040A (ko)
CA (1) CA2349224A1 (ko)
WO (1) WO2000032016A2 (ko)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101685821B1 (ko) * 2015-06-18 2016-12-13 인하대학교 산학협력단 흉부 Ｘ－ｒａｙ 영상에서 몸체 및 관심영역 분할 방법 및 시스템

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3237582B2 (ja) * 1997-08-28 2001-12-10 富士ゼロックス株式会社 画像処理装置及び画像処理方法
US6430637B1 (en) * 1999-02-04 2002-08-06 Micron Technology, Inc. Method for multiplexing bus interfaces on a computer expansion bus
US6381351B1 (en) * 1999-11-24 2002-04-30 Direct Radiography Corp. Weighted inverse topography method for digital x-ray image data processing
US6771802B1 (en) * 2000-04-13 2004-08-03 Photon Imaging, Inc. Method and apparatus for imaging and localizing radiation
US20020196965A1 (en) * 2001-06-22 2002-12-26 Wallace Edward S. Image transformation and analysis system and method
US6795526B2 (en) 2002-03-04 2004-09-21 Ge Medical Systems Global Technology Co., Llc Automatic exposure control for a digital image acquisition system
US20040179651A1 (en) * 2003-03-12 2004-09-16 Canon Kabushiki Kaisha Automated quality control for digital radiography
WO2005057493A1 (en) * 2003-12-10 2005-06-23 Agency For Science, Technology And Research Methods and apparatus for binarising images
US7860308B2 (en) * 2005-12-01 2010-12-28 Yahoo! Inc. Approach for near duplicate image detection
US8873818B1 (en) 2013-01-11 2014-10-28 E. Theodore Ostermann System and method for image analysis with characteristic curves

Family Cites Families (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS59151940A (ja) * 1983-02-18 1984-08-30 株式会社東芝 Ｘ線診断装置
US4914295A (en) * 1985-06-25 1990-04-03 Fuji Photo Film Co., Ltd. Radiation image read-out and image signal storing apparatus
US4833625A (en) * 1986-07-09 1989-05-23 University Of Arizona Image viewing station for picture archiving and communications systems (PACS)
JPH0642882B2 (ja) * 1987-04-20 1994-06-08 富士写真フイルム株式会社 所望画像信号範囲決定方法
US4868651A (en) * 1988-05-17 1989-09-19 S&S Inficon, Inc. Digital radiography with image brightness and contrast normalization
US5029226A (en) * 1989-10-10 1991-07-02 Unisys Corporation Method and apparatus for effecting spot/void filtering of image data
US5046118A (en) * 1990-02-06 1991-09-03 Eastman Kodak Company Tone-scale generation method and apparatus for digital x-ray images
US5283736A (en) * 1990-11-29 1994-02-01 Konica Corporation Radiographic image processing apparatus
DE69227008T2 (de) * 1991-11-14 1999-05-20 Agfa-Gevaert N.V., Mortsel Verfahren und Vorrichtung zum Herstellen eines Histogramms in der digitalen Bildverarbeitung mittels statistischer Pixelabtastung
US5164993A (en) * 1991-11-25 1992-11-17 Eastman Kodak Company Method and apparatus for automatic tonescale generation in digital radiographic images
US5268967A (en) * 1992-06-29 1993-12-07 Eastman Kodak Company Method for automatic foreground and background detection in digital radiographic images
JPH06154198A (ja) * 1992-11-27 1994-06-03 Toshiba Corp Ｘ線診断装置
US5588071A (en) * 1994-10-26 1996-12-24 Minnesota Mining And Manufacturing Company Identifying an area of interest using histogram data arranged in predetermined sequence
JPH08235355A (ja) * 1995-02-28 1996-09-13 Minolta Co Ltd 画像処理装置
US5641267A (en) * 1995-06-06 1997-06-24 General Electric Company Controlled leakage shroud panel
US5574764A (en) * 1995-06-06 1996-11-12 General Electric Company Digital brightness detector
US5757022A (en) * 1995-10-19 1998-05-26 Fuji Photo Film Co., Ltd. Image processing apparatus
US5633511A (en) * 1995-12-22 1997-05-27 Eastman Kodak Company Automatic tone scale adjustment using image activity measures
US5606587A (en) * 1996-03-21 1997-02-25 Eastman Kodak Company Determination of direct x-ray exposure regions in digital medical imaging
US5982916A (en) * 1996-09-30 1999-11-09 Siemens Corporate Research, Inc. Method and apparatus for automatically locating a region of interest in a radiograph
EP1008279A4 (en) * 1997-02-21 2002-04-17 Direct Radiography Corp IMAGE DATA PROCESSING FOR X-RAY RADIATION DETECTORS
US5903660A (en) * 1997-07-16 1999-05-11 The Regents Of The University Of California Automatic background recognition and removal (ABRR) in projection digital radiographic images (PDRI)
US6018590A (en) * 1997-10-07 2000-01-25 Eastman Kodak Company Technique for finding the histogram region of interest based on landmark detection for improved tonescale reproduction of digital radiographic images

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101685821B1 (ko) * 2015-06-18 2016-12-13 인하대학교 산학협력단 흉부 Ｘ－ｒａｙ 영상에서 몸체 및 관심영역 분할 방법 및 시스템

Also Published As

Publication number Publication date
US6118892A (en) 2000-09-12
EP1131776A2 (en) 2001-09-12
CA2349224A1 (en) 2000-06-02
JP2002530180A (ja) 2002-09-17
WO2000032016A3 (en) 2000-09-08
WO2000032016A2 (en) 2000-06-02
EP1131776A4 (en) 2005-03-16

Similar Documents

Publication Publication Date Title
US5471987A (en) Method of compressing a dynamic range for a radiation image
US5046118A (en) Tone-scale generation method and apparatus for digital x-ray images
JP3467285B2 (ja) 放射線画像処理方法
JPH0614193A (ja) ディジタル放射線画像における自動階調スケール生成方法
JP3538244B2 (ja) ラジオグラフィー像の表示において飽和画素を捜し出す方法
US6744849B2 (en) Image processing apparatus, image processing method, program, and storage medium
KR20010101040A (ko) 여과된 히스토그램을 이용한 디지털 엑스레이 검출기를위한 인터레스트 영역의 자동검출 방법
US6408046B2 (en) Method for dynamic range management of digital radiographic images
JP2001512607A (ja) ディジタル式ｘ線検出器に使用する画像データ処理法
US5513273A (en) Method for obtaining information about interstitial patterns of the lungs
JP3814864B2 (ja) 放射線画像の画像処理条件決定装置及び画像処理装置並びに放射線画像の画像処理条件決定方法及び画像処理方法
JPH0792588A (ja) 光刺激性蛍光スクリーンに記録された放射線画像を読み取るための装置の感度をモニターする方法
JPH06292008A (ja) 放射線画像のダイナミックレンジ圧縮処理装置
JPH09106448A (ja) 異常陰影候補の検出方法
JP4129598B2 (ja) 画像処理装置及び画像処理方法
JP3455566B2 (ja) 放射線画像のダイナミックレンジ圧縮方法
JP2864162B2 (ja) 放射線画像処理装置
JP2000023950A (ja) 放射線画像の画像処理装置
JP3765320B2 (ja) 放射線画像の画像処理条件決定装置及び画像処理装置
JPH06339025A (ja) 放射線画像のダイナミックレンジ圧縮方法
JP2000157518A (ja) 放射線画像処理装置
JP3316630B2 (ja) 腹部放射線画像の画像処理装置
JP3185105B2 (ja) 放射線画像処理装置
JP3637049B2 (ja) 放射線画像のダイナミックレンジ圧縮方法
EP0527524B1 (en) Image quality enhancing method

Legal Events

Date Code Title Description
PA0105 International application

Patent event date: 20010519

Patent event code: PA01051R01D

Comment text: International Patent Application

PG1501 Laying open of application
A201 Request for examination
PA0201 Request for examination

Patent event code: PA02012R01D

Patent event date: 20041112

Comment text: Request for Examination of Application

E902 Notification of reason for refusal
PE0902 Notice of grounds for rejection

Comment text: Notification of reason for refusal

Patent event date: 20060726

Patent event code: PE09021S01D

E601 Decision to refuse application
PE0601 Decision on rejection of patent

Patent event date: 20061122

Comment text: Decision to Refuse Application

Patent event code: PE06012S01D

Patent event date: 20060726

Comment text: Notification of reason for refusal

Patent event code: PE06011S01I