KR20000052875A - 기준 이미지와 x-레이 이미지와의 정합용 장치 - Google Patents

기준 이미지와 x-레이 이미지와의 정합용 장치 Download PDF

Info

Publication number
KR20000052875A
KR20000052875A KR1019990703726A KR19997003726A KR20000052875A KR 20000052875 A KR20000052875 A KR 20000052875A KR 1019990703726 A KR1019990703726 A KR 1019990703726A KR 19997003726 A KR19997003726 A KR 19997003726A KR 20000052875 A KR20000052875 A KR 20000052875A
Authority
KR
South Korea
Prior art keywords
dpis
dsis
image
generating
images
Prior art date
Application number
KR1019990703726A
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
Family has litigation
First worldwide family litigation filed litigation Critical https://patents.darts-ip.com/?family=24973119&utm_source=google_patent&utm_medium=platform_link&utm_campaign=public_patent_search&patent=KR20000052875(A) "Global patent litigation dataset” by Darts-ip is licensed under a Creative Commons Attribution 4.0 International License.
Application filed by 유니버시티 오브 피츠버그 오브 더 코몬웰츠 시스템 오브 하이어 에듀케이션 filed Critical 유니버시티 오브 피츠버그 오브 더 코몬웰츠 시스템 오브 하이어 에듀케이션
Publication of KR20000052875A publication Critical patent/KR20000052875A/ko

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/33Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/74Image or video pattern matching; Proximity measures in feature spaces
    • G06V10/75Organisation of the matching processes, e.g. simultaneous or sequential comparisons of image or video features; Coarse-fine approaches, e.g. multi-scale approaches; using context analysis; Selection of dictionaries
    • G06V10/757Matching configurations of points or features
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048Monitoring, verifying, controlling systems and methods
    • A61N5/1064Monitoring, verifying, controlling systems and methods for adjusting radiation treatment in response to monitoring
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10081Computed x-ray tomography [CT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10STECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S378/00X-ray or gamma ray systems or devices
    • Y10S378/901Computer tomography program or processor

Abstract

방사선요법 포털 이미지와 시뮬레이션 이미지와 같은 X-레이 이미지(53)는, 비교용으로 디스플레이될 수 있는 디지털화된 정합 신호를 생성하기위해 이미지를 디지털화하고 디지털화된 신호를 자동적으로 처리하는 장치에 의해 정합된다. 디지털화된 이미지는, 환자위에 위치한 X-레이 불투명 기점의 검출 그리고 식별(120)을 통해, 또는 2개의 이미지로부터 상호작용하여 선택된 시드점(110)으로부터 생성된 변환을 사용하여 먼저 대략 정렬이 이루어진다. 이후 2개 이미지(150)의 영역을 우선 교차시키고 상기 영역(154)을 향상시킴으로써 미세 정렬(35)은 수행된다. 두번째로, 연속적인 향상 레벨의 해상도에서 상기 영역내의 로버스트 이동흐름을 사용하여 갱신된 변환이 발생한다(160). 이후 갱신된 변환은 비교를 위해 디스플레이되는 이미지(167)를 정렬하는데 사용된다. 방사선요법 장비(140)를 제어하기위해 갱신된 변환이 또한 사용될 수 있다.

Description

기준 이미지와 X-레이 이미지와의 정합용 장치{APPARATUS FOR MATCHING X-RAY IMAGES WITH REFERENCE IMAGES}
응용 의학에는 X-레이 이미지의 정합이 필요하다. 예를 들어, 컴퓨터 제어 방사선요법에서, 종양을 둘러싸는 건강한 조직이 최소한 노출되면서 종양에 대한 방사선치료를 최대화하기위해 고 에너지 방사의 치료광선이 다수의 방향으로부터 종양으로 향해진다. 이러한 방사선요법 치료는 특히 시뮬레이션 위상, 및 실제 치료 위상인 2개의 특별한 위상을 갖는다. 시뮬레이션 위상에서, 고 에너지 방사선이 발생하지 않는 장비를 제외하고 환자는 치료 장비와 유사한 장비위에 위치하게 된다. 종양 치료 전문가에 의해 상기 설명되는 치료광선의 시퀀스 전달을 시뮬레이션하기위해 시뮬레이션 장비는 연속적으로 위치된다. 따라서 요구되는 치료광선을 전달하도록 그리고 장비와 환자사이에서 또는 장비의 이동성 구성요소사이에서 충돌없이 한개의 치료광선으로부터 그다음 치료광선으로 이동하도록 장비는 위치될 수 있다. 이러한 프로시저 동안, 시뮬레이션 이미지라 불리는 저 조사량 X-레이 이미지가 이용된다. (킬로전자 볼트범위에서) 사용되는 엑스레인선의 저 에너지 때문에 일반적으로 좋은 대조 및 세부를 갖는 상기 시뮬레이션 이미지로 인해 종양치료 전문가는 종양의 위치를 지정할 수 있고 이에따라 연속적인 치료광선을 전달하기위해 장비 구성요소의 위치를 확립한다.
실제 치료위상동안, 환자를 치료하기위해 전형적으로 메가전자볼트 범위의 정규 조사량 X-레이 방사선이 사용되기전에 시뮬레이션에서와 같이 환자는 장비위에서 정확히 동일한 위치에 배치된다. 상기 위상동안, 또다른 X-레이 이미지가 취해지고, 포털 이미지라 불린다.
치료의 완료후, 시뮬레이션 이미지와 포털 이미지는 시뮬레이션 이미지에서 식별된 바와같은 종양이 포털 이미지에서 방사선으로 적절히 치료되었는지를 결정하기위해 전문가에 의해 비교된다. 치료가 완전하지 않다면, 환자는 정정 치료를 받게된다.
현재 허용되는 치료프로시저에는 포털 이미지와 시뮬레이션 이미지와의 수동 비교가 포함된다. 상이한 장비에서 그리고 상이한 방사선 노출 정도로 2개의 X-레이가 흔히 이용된다는 사실에서 볼 때 정밀한 수동 비교가 매우 요구된다. 후자의 경우는 종양 영역이 흔히 포털 X-레이에서 보이지 않는다는 것을 의미하며, 따라서 포털 이미지를 시뮬레이션 이미지와 정합하는 것은 선명하게 보이지 않을 해부용 표시 차원인 수동 평가에 의존해야 한다.
종래에, 포털 이미지는 개발될 필요가 있는 X-레이 필름을 사용하여 생성되었다. 오직 한개 또는 몇개의 치료광선만이 이용된다는 것은 심각한 단점이 아니다. 그러나, 상기 X-레이 필름은 다수의 치료광선이 사용되는 컴퓨터 제어 방사선요법에서 상당한 제한요소이다. 전자 디스플레이 디바이스에 디스플레이될 수 있는 디지털화된 이미지를 생성하는 전자 포털 이미저가 개발되었다. 불행히도, 전자적으로 생성된 포털 이미지의 대조와 선명도에서도 같은 문제가 존재한다.
포털 이미지를 시뮬레이션 이미지와 정합시키는 문제에는 스큐, 스케일링 차이, 회전, 번역 및 차이로 인한 정위의 차이를 상기 이미지가 갖는다는 점이 있다.
정위 방사선학에서, 상기 2개의 프레임에 나타나는 공지된 고정 차원의 스틸 프레임으로부터 얻어진 스케일링을 가함으로써 디지털화된 컴퓨터 토모그래피 X-레이 이미지 및 자기 공명 이미지(MRI)가 자동적으로 정합되었다. 이러한 공지된 차원의 고정 표시는 종래의 방사선요법 이미지에서 이용할 수 없다.
따라서, 방사선요법에서 시뮬레이션 이미지를 X-레이 이미지와 자동적으로 정합시키고 특히 포털 이미지를 시뮬레이션 이미지와 정합시키는 장치가 필요하다.
다중 치료광선을 위한 시뮬레이션 이미지와 포털 이미지를 정합시킬 수 있는 이러한 장치가 필요하다.
스큐, 회전, 스케일링, 만곡, 또는 원근으로 인한 매우 상이한 대조, 선명도, 및 차이를 갖는 시뮬레이션 이미지와 포털 이미지를 정합시킬 수 있는 장치가 또한 필요하다.
컴퓨터 제어 방사선요법동안 환자의 정렬을 얻고 유지하며 또는 정렬이 허용가능하지 않다면 방사를 종료시키는 장치가 또한 필요하다.
발명의 개요
이러한 요구는, X-레이 이미지와 기준 이미지를 자동적으로 정합시키며, 방사선요법 치료가 적절한 것이었는지를 결정하기위해 또는 방사선요법 치료의 동작을 제어하도록 연속 포털 이미지를 정합시키기 위해, 특히 포털 이미지와 시뮬레이션 이미지를 정합시키는 장치에 대한 본 발명에 의해 해결된다. 이미지를 정합하는데 있어서, 디지털화 수단은 제 1 세트의 디지털 이미지 신호 즉 포털 이미지인 경우 디지털 포털 이미지 신호(DPIS)를 생성하기위해 포털 이미지와 같은 X-레이 이미지를 디지털화한다. 디지털화 수단은 또한 제 2 디지털 이미지 신호 즉 디지털 시뮬레이션 신호(DSIS)를 생성하기위해 시뮬레이션 이미지와 같은 기준 이미지를 디지털화한다. 처리 수단은 이러한 디지털 이미지를 처리하여 정합된 디지털 이미지 신호를 생성한다. 상기 처리는 이미지에서 어떠한 특징의 물리적 차원의 종래 지식없이 수행된다. 출력 수단은 예를 들어 정합된 디지털 이미지 신호로부터 디스플레이를 생성하며 그리고/또는 치료/진단 장비를 제어한다.
처리 수단은 우선 디지털 포털 이미지와 디지털 시뮬레이션 이미지 신호사이의 대략 정렬에 영향을 주는 대략 정렬 수단을 포함한다. DPIS로 표시되는 포털 이미지 그리고 DSIS로 표시되는 시뮬레이션 이미지내의 시드점을 선택함으로써 대략 정렬은 시작된다. 2개 이미지의 디스플레이위에 상응하는 점이 무엇인가를 선택하기위해 마우스와 같은 포인팅 디바이스를 사용하여, 또는 자동적으로 환자위에 배치된 X-레이 불투명 기점을 사용함으로써 시드점 선택은 상호작용식으로 행해질 수 있다. 어느 경우에서든지, 2개의 이미지사이의 변환을 계산하기위해 시드점이 사용된다. 이후 한 이미지에 있는 점을 다른 쪽 이미지의 좌표로 변환하기위해 디지털 이미지 신호 세트중 한개에 상기 변환을 적용하며, 이에따라 대략 정렬된 DPIS와 DSIS가 발생하는 수단이 사용된다.
대략 정렬 이후에, 미세 정렬이 수행된다. 상기 미세 정렬을 행하는데 있어서, 우선 교차하거나 중첩하는 이미지의 영역을 위해, 그리고 이미지의 교차영역내에서 직각과 같이 일반적인 형태의 영역을 위해 선택된 DPIS와 DSIS를 선택함으로써 대략 정렬된 DPIS와 DSIS가 준비된다. 이후 이러한 영역을 위한 디지털 이미지 신호는 향상되어 유사한 동적 범위와 픽셀 명암도를 갖는 준비된 이미지가 발생한다. 미세 정렬 수단은, 정합된 DPIS와 DSIS를 생성하기위해, 준비된 DPIS와 DSIS로부터 갱신된 변환을 생성하는 수단, 및 대략 DPIS와 DSIS 또는 준비된 DPIS와 DSIS에 갱신된 변환을 적용하는 수단을 포함한다.
갱신된 변환을 생성하는 수단은 준비된 DPIS와 DSIS로부터 이동흐름 성분을 생성하는 수단, 및 이동흐름 성분으로부터 갱신된 변환을 계산하는 계산 수단을 포함한다. 바람직하게, 이동흐름 성분을 생성하는 수단은 이동흐름 변화도 구성요소를 생성하고 계산 수단은 갱신된 변환을 계산하기위해 로버스트 최적화를 적용하는 수단을 포함한다. 갱신된 변환을 생성하는 수단은 갱신된 변환을 생성하기위해 준비된 DPIS와 DSIS의 연속적인 향상레벨의 해상도를 사용한다.
추적 모드에서, 연속 세트의 디지털 포털 이미지 신호사이의 움직임을 추적하기위해 갱신된 변환이 사용된다. 환자 움직임이 특정 한계를 초과한다면 방사를 종료시키기위해 추적이 사용될 수 있고, 또는 치료받을 영역이 적절한 정렬상태로 방사선을 유지하도록 환자 위치지정 어셈블리를 동작시키기위해 사용될 수 있다.
본 발명은, 또한 또다른 X-레이 이미지 또는 또다른 형태의 이미지일 수 있는 기준 이미지와 X-레이 이미지를 정합시키는데 사용될 수 있다.
본 발명은, 유사한 X-레이 이미지의 정합에 관한 것이며, 판단을 위해, 그리고/또는 치료 장비를 제어하기위해 원하는 타겟이 실질적으로 방사선치료를 받는 것을 측정하도록 시뮬레이션 기계에서 치료전에 생성된 시뮬레이션 이미지와 치료기계의 방사선요법 치료동안 생성된 포털 이미지를 온라인상에서 자동적으로 정합시키는 컴퓨터로 제어되는 방사선요법 장치에 특별히 응용된다.
본 발명을 전체 이해하는 것은 다음에 따르는 바람직한 실시예의 상세한 설명과 첨부된 도면으로 가능하다.
도 1은 본 발명을 실현하는 장치의 개략적인 도이다.
도 2a는 본 발명이 적용될 수 있는 시뮬레이션 이미지의 간략한 도이다.
도 2b는 본 발명이 적용될 수 있는 포털 이미지의 간략한 도이다.
도 2c는 본 발명을 이용하는 도 2a와 도 2b의 시뮬레이션 이미지와 포털 이미지를 중첩하는 디스플레이의 간략한 도이다.
도 3 내지 도 11은 도 1 장치에서 본 발명을 실현하는데 이용되는 소프트웨어의 흐름도이다.
본 발명은 X-레이 이미지를 기준 이미지와 정합시키는 것이며 컴퓨터 제어 방사선요법에서 생성된 포털 이미지를 시뮬레이션 이미지와 정합하는데 적용될 설명될 것이다. 그러나, 본 발명은 예를 들어 진단에서 사용된 것과 같이 다른 X-레이 이미지를 정합하는 데 있어서 응용 폭이 넓다는 것을 이해할 것이다. 명백하듯이, 본 발명은 또한 방사선의 게이팅 또는 환자의 위지지정을 제어하기위한 것과 같은 연속 포털 이미지로 움직임을 추적하는 응용을 갖는다.
도 1에서, 치료받을 환자(3)내부의 종양과 같은 영역의 위치를 결정하기위해 그리고 치료광선의 시퀀스를 설정하기위해 시뮬레이션 설정(1)이 사용된다. 설정 장비는 기계 베이스(9)에 의해 지지되는 수평 피봇(7)에 대해 회전용으로 장착된 갠트리(5)를 포함한다. 킬로전자볼트 범위의 저 에너지인, X-레이 광선(1)은 갠트리(5)위에 장착된 콜리메이터(13)에 의해 피봇(7)의 연장을 통해 대각선으로 연장되는 경로를 향한다.
환자(3)는 지지대에 대하여 3차원 번역을 위해 지지대(19)위의 소파(17)를 포함하는 환자 위치지정 어셈블리(15)위에서 지지된다. 지지대(19)는 회전테이블(21)위에 장착된다. 소파(17)의 번역을 통해, 피봇(7)에 대한 갠트리(5)의 회전과 회전테이블(21)의 회전, 복수의 치료광선이 시뮬레이션될 수 있다. 시뮬레이션 장비(1)를 연속 광선을 생성하기위해 필요한 위치에 배열함으로써, 필요한 모든 광선이 얻어질 수 있는지 여부 그리고 장비와 환자사이 또는 장비의 구성요소들 사이의 충돌을 피하기위해 광선 사이의 장비의 움직임이 조절될 수 있는지 여부가 결정될 수 있다.
콜리메이터(13)로부터 환자(3)의 다른 쪽에 X-레이 필름(23)을 X-레이 선(11)과 정렬하여 배치함으로써 시뮬레이션 이미지를 생성하기위해 저 에너지 X-레이 선(11)이 사용된다. 피봇 축(7)의 투영과 광선(11)의 교차인 설정의 등선량중심점에서 종양과 같은 치료받을 환자의 영역을 위치지정하는데 상기 시뮬레이션 이미지가 사용된다.
시뮬레이션의 완료후에, 환자(3)는 치료 설정(1')으로 이송된다. 도시된 바와같이, X-레이 광선(11')이 메가전자볼트 범위내라는 것을 제외하고 치료 설정(1')은 시뮬레이션 설정(1)과 유사하다. 포털 이미지는 치료 설정(1')에 의해 생성된다. 상기 포털 이미지는 시뮬레이션 설정에서와 같이 X-레이 필름에 의해 캡쳐될 수 있다. 그러나, 전자 포털 이미저(25)가 사용되는 것이 바람직하다. 가능하다면, 시뮬레이션 설정(1)에서 X-레이 필름(23)대신에 전자 이미저가 또한 사용될 수 있다.
상기에서와 같이, 시뮬레이션 이미지 및 포털 이미지는 상당히 상이할 수 있다. 이렇게 상이한 주요 이유중 한가지는 광선(11 및 11')의 에너지 차이에 있다. 치료 조사량이 적절한 치료영역에 전달되었는지를 결정하기위해 시뮬레이션과 포털 이미지를 정합시키는 본 발명이 사용될 수 있다. 움직임이 적정 한계를 초과한다면 X-레이 광선(11')의 생성을 종료시키기위해, 또는 장비를 조정하여 적절한 정렬을 유지하기위해 치료받는동안 환자 움직임을 검출하는데 본 발명이 또한 사용될 수 있다.
이미지 정합 시스템(27)은, X-레이 필름(23)에서 발생된 것과 같은 시뮬레이션 이미지, 및 전자 포털 이미저(25)에 의해 생성된 것과 같은 포털 이미지를 디지털화하는 디지타이저(29)를 포함한다. 보다 일반적인 의미에서, 정합 시스템(27)은 포털 이미지와 같은 X-레이 이미지를 시뮬레이션 이미지와 같은 기준 이미지와 정합시킨다.
이미지 정합 시스템(27)은 미세 정렬(35)용 모듈이 뒤따르는 대략 정렬(33)용 모듈을 포함하는 프로세서(31)를 더 포함한다. 프로세서의 출력은 디스플레이 디바이스(37)에 나타나는 정합된 포털(X-레이) 및 시뮬레이션(기준) 이미지일 수 있다. 디스플레이 디바이스(37)와 관련된 것으로는 마우스 또는 트랙볼과 같은 포인팅 디바이스(43)와 키보드(41)를 포함할 수 있는 인터페이스 디바이스(39)가 있다.
본 발명이 포털 X-레이 이미지를 시뮬레이션 기준 이미지와 정합시키는 것을 도 2a 내지 도 2c는 도시한다. 도 2a는 시뮬레이션 설정(1)을 사용하여 생성된 시뮬레이션 이미지(45)를 나타낸다. 상기 이미지를 위해 사용되는 저 에너지 X-레이는 대조와 세부를 갖는 이미지를 발생시키며, 따라서 종양(51)뿐만 아니라 환자와 뼈 구조(49)의 윤곽이 도시된다. 도 2b에서는, 보다 높은 에너지 치료 광선을 갖는 포털 이미지가 균일하게 어두운 부분과 같은 치료 영역(55)을 도시한다. 광선(11')을 종양의 형태에 대체적으로 맞추기위해 치료 영역(55)의 불규칙한 모서리는 콜리메이터(13)에서 사용되는 리브에 의해 발생한다. 포털 이미지(55)의 나머지는 매우 적은 세부를 도시하며 뼈의 위치를 나타내지 않는다. 명백하듯이, 2개의 이미지(45및 53)는 서로에 대하여 번역될 수 있고, 상이하게 스케일링될 수 있으며, 스큐되고 (예를 들어 90도로) 회전할 수 있다. 또한 2개의 이미지는 투시적이며 만곡되어 상이할 수 있다.
대략 정렬 모듈(33)은 2개 이미지의 일반적인 정렬을 발생시키며, 이후 미세 정렬 모듈(35)은 로버스트 이동흐름을 사용하여 완벽한 이미지 정합을 빠르고도 정밀하게 행한다. 디스플레이 디바이스(37)는 상이한 방식으로 정합된 이미지를 표시한다. 일실시예에서, 디스플레이(37)는 약 6 내지 20Hz로, 보통 약 12Hz로, 2개의 이미지사이에서 교체되고, 도 2c에 도시된 바와같이, 관찰자는 합성 이미지(59)처럼 중첩된 이미지를 보게된다. 예에서 도시된 바와같이, 정합된 포털 이미지의 치료 영역(55')은 정합된 시뮬레이션 이미지내의 종양(51')과 중첩된다. (도시되지 않은) 디스플레이의 또다른 형태에서, 포털 이미지로부터의 치료 영역의 윤곽은 처리되는 시뮬레이션 이미지위로 투영되고, 따라서 타겟인 종양이 실제로 치료받는다면 보일 것이다.
대략 정렬을 수행하는데 있어서, 대략 변환이 디지털화된 X-레이 또는 포털 이미지 신호(DPIS)에 적용되어 디지털 기준 또는 시뮬레이션 이미지 신호(DSIS)의 좌표 시스템으로 변환된다. 명백하듯이, 상기 변환을 생성하는데 필요한 정보는, 상호작용하여 2개의 이미지에서 상응하는 점으로 보이는 것을 포인터 디바이스(43)를 사용하는, 또는 시뮬레이션 설정 및 치료 설정(도 1 참조)에서 환자위에 배치된 X-레이 불투명 기점(61)을 자동적으로 사용하는 조작자에 의해 선택함으로써 생성될 수 있다. 이렇게 생성된 점은 시드점이라 불린다. 포털 이미지 좌표로부터 시뮬레이션 좌표로의 대략 변환(H)은 수학식 1과 같다.
(x y)벡터는 시뮬레이션 이미지와 상응하는 포털 이미지에 있는 각 시드점의 중심의 행과 열 좌표를 나타낸다. 행렬에서 4개의 로트스큐스케일 성분은 이미지를 대략적으로 정렬하는데 필요한 전체 아핀 변환을 설명한다. 이 단계에서, 그 다음 단계가 상당히 작은 편향을 조절할 수 있기 때문에 시드점의 상호작용적 선택 또는 기점의 배치는 정밀할 필요가 없다.
대략 정렬의 결과를 사용하여, 포털 이미지는 시뮬레이션 이미지를 향해 휘게 된다. 이후, 2개 이미지의 중첩 영역은 컴퓨터로 향상되어 상응하는 명암도 레벨이 유사하게 된다. 마지막으로, 이동흐름, 또는 미세 스케일 변환이 계산되며 포털 이미지는 시뮬레이션 이미지를 향해 차이 변화도로 미끄러진다. 이 단계에서, 입력 위치 벡터가 수학식 2로 표시되는 보다 포괄적인 변환 모델이 사용된다.
그리고 변환 행렬이 수학식 3으로 표시된다.
따라서 결과는 수학식 4이다.
여기서, Δ포털이고,이다. 파라미터(α0내지 α5)는 대략 정렬에서처럼 아핀 변환을 포함하며, 파라미터(P0, P1)는 투시적 변환을 포함하고, c는 예를 들어 호흡에 의해 발생될 수 있는 변형을 다룬다.
벡터()의 파라미터를 복구하기위해, 이동흐름의 결과로서 이미지 차이를, 또는 2개 이미지사이의 거리를 공식화한다.
여기서 ∀xf이고, f는 변형이 계산되는 이미지의 영역이다. 수학식 5에서, I(x)는 점(x)에서의 명암도 함수이고, t+1에서의 이미지는 포털 이미지이며, t에서는 시뮬레이션 이미지이다. 다양한 미분 기술을 사용하여, 수학식 6과 같은 변화도(또는 차이 변화도)를 사용하여 이동흐름이 공식화된다.
이 단계에서, (수학식 3 참조) Q 원소를 계산하기위해 비구속 최적화를 사용하여 로버스트 복귀 방법이 사용된다. 이러한 기술로서 2개의 이미지사이의 잉여 오차 차이뿐만 아니라 대략 정렬 단계로부터의 상당히 작은 편차에 대처할 수 있다. 로버스트 기술을 사용함으로써 노이즈와 잉여 오차에 의한 영향을 받을 수 있는 위험없이 주된 변환만이 복구될 것이다.
도 3 내지 도 11은 본 발명을 실현하는 소프트웨어의 흐름도이다. 도 3은, 블록(120)에서 자동적으로 또는 블록(110)에서 상호작용적으로 대략 정렬을 수행하는 것을 포함하는 주요 루틴(100)을 도시한다. 어느 경우에서든지, 수학식 1을 사용하여 포털 이미지와 시뮬레이션 이미지사이에서 변환의 대략적인 근사화가 계산된다. 이후 사용자는 대략적인 근사로 인해 130에서 이미지의 만족할만한 정렬이 제공되었는지 여부를 결정할 수 있다. 제공되었다면, 상기 프로시저는 완료된다. 그렇지 않다면, 미세 정렬이 수행된다. 설명된 바와같이, 환자 움직임을 추적하는데 본 발명이 또한 사용될 수 있으며, 이 경우 기점의 갱신된 위치를 대략적으로 결정하기위해 140에서 2개의 이미지사이의 변환이 이용된다. 추적동안 그리고 이미지 정합에서 사용자에 의한 요구가 있다면, 150에서 미세 정렬을 위한 이미지가 준비된다. 이후 160에서 정밀한 이미지 변환이 계산되며 이미지 정합 모드가 170에서 결정된 바와같이 선택된다면, 변환이 달성되며 180에서 상기 설명된 방식으로 이미지가 디스플레이된다. 190에서 추적 모드가 선택되었다면, 그다음 위치를 생성하기위해 200에서 140으로 루틴이 복귀된다. 다시 사용자는 이미지 정합이 만족할만지 여부를 결정하기위해 200에서 최종 결정을 한다. 만족스럽지 않다면, 루틴은 110으로 복귀되고 대략적인 계산이 다시 시작된다.
도 3의 블록(110)에서 상호작용적으로 요구되는 변환의 대략적인 근사를 계산하는 프로시저가 도 4에서 상세히 도시된다. 사용자는, 예를 들어 111에서 표시된 바와같이 마우스(43)를 사용하여 시뮬레이션 이미지와 포털 이미지의 영역 또는 상응하는 시드점을 선택한다. 선택된 영역 또는 점은 이후 112에 표시된 바와같이 프로시저(A)를 호출함으로써 포털 이미지와 시뮬레이션 이미지사이의 대략적인 변환을 계산하는데 사용된다. 이러한 대략적인 변환은 113에서 표시된 바와같이 프로시저(B)를 호출함으로써 포털 이미지를 시뮬레이션 이미지 좌표로 변환하는데 사용된다. 이후 이미지는 114에 표시된 바와같이 모니터(37)에 디스플레이된다.
대략적인 변환을 계산하는데 사용되는 프로시저(A)의 상세한 사항은 도 5에 도시된다. 사용자가 A1에서 결정된 것과 같은 영역을 표시했다면, 시스템은 A2에 표시된 것과 상응하는 영역 내부로부터 자동적으로 임의 점을 선택한다. 이후, 사용자가 영역이 아닌 점을 선택했다면, A3에 표시된 바와같은 최소 스퀘어(LSQ)방법을 사용하여 변환 파라미터를 계산하기위해 상응하는 영역 쌍이 사용된다.
포털을 시뮬레이션 좌표로 변환하는 프로시저(B)의 세부사항은 도 5에 도시된다. 먼저, 변환된 포털 이미지의 행과 열 극한은 수학식 1의 역인 변환 행렬(H)을 사용하여 B1에서 결정된다. 이후 결과로 나타나는 포털 이미지는 B2에서 래스터 주사되고, 각 픽셀을 위해 변환을 사용하여 위치가 결정된다. 초기 포털 이미지에서 주위의 픽셀 위치사이의 선형 보간법을 사용하여 각 픽셀을 위한 명암도값이 다음에 계산된다.
환자위의 기점을 사용하여 대략 정렬을 자동적으로 수행하는 루틴(120)은 도 7에 도시된다. 121에서 포털과 시뮬레이션 이미지에서 X-레이 불투명 기점(61)이 검출되고 상응하는 표시가 122에서 식별된다. 이후 도 5의 프로시저(A) 및 시드점인 각 표시의 도심을 사용하여 이미지 변환이 123에서 계산된다. 이후 도 6의 프로시저(B)와 계산된 변환을 사용하여 124에서 포털 이미지가 시뮬레이션 좌표로 변환된다. 125에서 결정된 바와같이 정합 모드일때, 이미지는 도 2a 내지 2c와 함께 상기 설명된 방식으로 126에서 디스플레이된다.
미세 정렬을 위해 대략 정렬된 디지털 이미지 신호를 준비하는 루틴(150)이 도 8에 도시된다. 우선, 포털과 시뮬레이션 이미지사이의 중첩에 대한 교차 영역은 수학식 1의 변환을 사용하여 151에서 계산된다. 이후, 교차 영역에 꼭 맞는 가장 큰 정사각형 영역이 152에서 계산된다. 사각형 등과 같이 다른 일반적인 기하학적 형태가 직사각형대신에 사용될 수 있다. 포털과 시뮬레이션 이미지의 직사각형 교차 영역을 나타내는 새로운 이미지는 153에서 형성된다. 이렇게 나타나는 이미지는 이후 준비된 디지털 이미지 신호를 생성하기위해 154에서 향상된다. 유사한 동적 범위과 픽셀 명암도를 갖는 준비된 이미지를 발생시키기위해 막대그래프 균일화, 가우스의 라플라시안, 고주파 필터링 및 다른 기술과 같이 향상된 다양한 형태가 사용된다.
도 9는 미세 정렬을 위해 갱신된 변환을 계산하는 루틴(160)을 도시한다. 이 프로세스는, 본 예에서 1/8 해상도인 최저 해상도로 시작되는 디지털 이미지 신호의 여러 레벨의 해상도에서 수행된다. 따라서, 161에서 준비된 포털과 시뮬레이션 이미지를 위한 최저 해상도에서의 이미지가 형성된다. 이러한 이미지는 가장 최근의 갱신된 변환 파라미터, 즉, 162에서 이전 레벨의 해상도에서 계산된 변환 파라미터를 사용하여 갱신된다. 본 발명의 중요한 부분은 미세 정렬을 수행하기위해 로버스트 이동흐름이 사용된다는 것이다. 특히, 이동흐름 변화도 구성요소는 163에서 생성된다. 변화도 구성요소를 사용하는 이동흐름의 응용은, 1993년 5월 독일 베를린 컴퓨터 비전(ICCV 93)의 제 4 차 국제회의에서 발표된 "A Framework for the robust estimation of optical flow"라는 M.J.Black과 P.Anandan의 논문에 설명되어 있다. 이동흐름은 한개 이미지의 픽셀이 다른 이미지의 상응하는 픽셀과 정렬되기위해 필요한 움직임에 가해진다. 로버스트 움직임은 이동한 대부분의 픽셀이 유사하게 이동한 만큼 움직임에 가해지는 한편, 상이한 움직임을 나타내는 다른 것이 있을 수 있다. 이후 갱신된 이미지 변환 파라미터가 로버스트 최적화를 사용하여 164에서 계산된다. 해상도의 상한이 165에서 결정된 바와같이 도달되지 않았다면, 해상도는 166에서 증가하며 갱신된 변환 파라미터는 새로운 레벨의 해상도에서 다시 계산된다.
최고 레벨의 해상도가 165에서 도달되었을 때, 최종 변환행렬(Q)이 167에서 생성된다. 도 9의 블록(164)의 로버스트 최적화를 사용하여 갱신된 이미지 변환 파라미터를 계산하는 루틴의 세부사항은 도 10에 도시된다. 상기 설명된 Black과 Anandan에 의한 논문에서 설명된 바와같이, 로버스트 움직임은 인라이어라 불리는 데이터 포인트에 의해 표시된다. 이렇게 존재하는 다른 움직임은 아웃라이너로서 식별된다. 본 발명에서, 데이터 포인트는 픽셀값이다. 픽셀은 일정한 이동흐름 벡터에 대한 정도에 근거하여 인라이어와 아웃라이어로 연속적으로 분리된다. 주요 이동흐름을 계산하기위해 인라이어 세트의 픽셀이 사용되며, 상기 벡터에 대한 정도는 로버스트 최적화동안 계산되는 가중 인자에 의존한다.
도 10에서, 개별적인 가중 인자를 사용하여 각 인라이어 포인트가 표시되는 164.1에서 루프가 시작된다. 처음에, 픽셀의 가중 인자는 모두 1로 설정되어 모든 픽셀이 인라이어로 된다. 164.2에서, 아웃라이어의 프로시저의 감도를 결정하는 최적화 파라미터(σ)가 설정된다. 가중 인자는 상기 파라미터(σ)에 의존한다. σ값이 작아질수록, 보다 많은 포인트가 인라이어로 제거되며 인라이어는 이동흐름 벡터의 현재 예측으로 보다 가깝게 된다. 따라서, 모든 포인트가 포함되도록 큰 σ이 처음에 사용된다. 연속적인 루프에서, σ는 σ스케쥴링이라 불린다. σ스케쥴링은 조심스럽게 행해져야 한다. σ이 너무 빨리 작아진다면, 해결책이 없을 수 있고, 반면에, σ을 너무 천천히 작게한다면 처리 시간이 길어지게 된다. 본 발명에 따라서, σ는 이동흐름 파라미터의 최대 에러에 의존하여 작아진다. 이후에, 164.4에서 Q행렬의 변환 파라미터를 위한 갱신된 값의 계산에서 각 인라이어 데이터 포인트가 사용되는 164.3에서 또다른 루프가 시작된다. 하강 변화도가 또한 사용될 수 있지만, 164.4에서 사용되는 수학식은 바람직하게 켤레 변화도를 사용하여 유도된다. 또한, 변환 파라미터를 결정하기위한 수학식을 유도하는데 이동흐름과 로버스트 통계가 사용된다. σ뿐만 아니라 최종 계산으로부터의 변화인 변환 파라미터의 에러는 픽셀을 위한 가중 인자를 조절하기위해 164.5에서 사용된다. 모든 인라이어 데이터 포인트/픽셀이 164.3에서 결정된 바와같이 사용되었을때, 해법이 원하는 정도까지 수렴되었는지 여부를 결정하기위해 164.6에서 검사가 행해진다. 그렇지 않다면, 루틴은 164.1로 복귀하며 갱신된 가중 인자를 사용하여 인라이어 데이터 포인트는 다시 표시된다.
도 11은 140에서의 추적 루틴을 도시한다. 141에서 나타난 바와같이, 증가 갱신 및 변환(H 그리고/또는 Q)은 결합되어 변환이 항상 초기 시뮬레이션 또는 기준 이미지와 다시 연관된다. 추적 루틴을 통한 초기 통과시에, 필요하다면 현재의 포털 이미지는 시뮬레이션 이미지를 대신하고, 이후 새로운 포털 이미지가 143에서 얻어진다. 추적이 계속됨에 따라, 연속적인 포털 이미지는 그다음 선행하는 포털 이미지와 정합되어 갱신된 변환이 생성된다. 144에서 나타난 바와같이, 광선을 온오프하고 그리고/또는 환자 위치지정 어셈블리를 구동하는 것과 같이 방사선요법 장비를 제어하도록 추적 신호를 생성하기위해 연속 포털 이미지로부터의 기점 패턴에 있어서의 변화 또는 기점의 연속 위치가 사용된다.
본 발명의 특정 실시예가 상세하 설명되었지만, 당 기술에 숙련된 당업자는 본 명세서의 전체 교시라는 관점에서 다양한 수정 및 변경이 발생할 수 있다는 것을 이해할 것이다. 따라서, 설명된 특정 장비는 단순히 도시된 것이며 본 발명의 사상을 제한하는 것이 아니고 본 발명에는 청구범위 및 청구범위와 등가인 것이 제공된다.

Claims (28)

  1. 포털 이미지(53)를 시뮬레이션 이미지(45)와 자동적으로 정합시키는 장치(27)에 있어서,
    디지털 포털 이미지 신호(DPIS)와 디지털 시뮬레이션 이미지 신호(DSIS)를 각각 생성하기위해 시뮬레이션 이미지(45)와 상기 포털 이미지(53)를 디지털화하는 수단(29);
    정합된 DPIS와 DSIS를 생성하기위해 상기 DPIS와 DSIS를 처리하는 처리 수단(31); 및
    상기 정합된 DPIS와 DSIS로부터 출력을 생성하는 출력 수단(37, 140)을 포함하는 것을 특징으로 하는 장치.
  2. 제 1 항에 있어서, 상기 처리 수단(31)은, 상기 DPIS와 DSIS로부터 대략 정렬된 DPIS와 DSIS를 생성하는 대략 정렬 수단(33), 및 상기 시뮬레이션 이미지와 포털 이미지의 영역을 중첩하기위해 상기 대략 정렬된 DPIS와 DSIS로부터 상기 정합된 DPIS와 DSIS를 생성하는 미세 정렬 수단(35)을 포함하는 것을 특징으로 하는 장치(27).
  3. 제 2 항에 있어서, 상기 대략 정렬 수단(33)은, 상기 DSIS에 의해 표시되는 상기 시뮬레이션 이미지(45)와 상기 DPIS로 표시되는 상기 포털 이미지(53)에서 상응하는 시드점을 선택하는 수단(111, 121, 122), 상기 상응하는 시드점으로부터 상기 시뮬레이션 이미지(45)와 상기 포털 이미지(53)사이의 변환을 계산하는 수단(112, 123), 및 상기 대략 정렬된 DPIS와 DSIS를 상기 DPIS와 상기 DSIS중 한 개와 함께 생성하기위해, 상기 변환을 상기 DPIS와 상기 DSIS중 한개에 적용하는 수단(113, 124)을 포함하는 것을 특징으로 하는 장치(27).
  4. 제 3 항에 있어서, 상응하는 시드점을 선택하는 상기 수단(111)은 상기 DPIS와 DSIS로부터 생성된 디스플레이에서 상응하는 점을 선택하는 상호작용 수단(43)을 포함하는 것을 특징으로 하는 장치(27).
  5. 제 3 항에 있어서, 상응하는 시드점을 선택하는 상기 수단(121, 122)은, 상기 DPIS와 DSIS에서 X-레이 불투명 기점을 검출하는 수단(121), 및 상기 DPIS와 DSIS의 상응하는 기점을 상기 상응하는 시드점으로 식별하는 수단(122)을 포함하는 것을 특징으로 하는 장치(27).
  6. 제 3 항에 있어서, 상기 미세 정렬 수단(33)은, 상기 대략 정렬된 DPIS와 DSIS로부터 준비된 DPIS와 DSIS를 생성하는 수단(150), 상기 준비된 DPIS와 DSIS로부터 갱신된 변환을 생성하는 수단(161 내지 165), 및 상기 정합된 DPIS와 DSIS를 생성하기위해 준비된 대략 DPIS와 DSIS중 한 개에 상기 갱신된 변환을 적용하는 수단(167)을 포함하는 것을 특징으로 하는 장치(27).
  7. 제 2 항에 있어서, 상기 미세 정렬 수단(33)은, 상기 대략 정렬된 DPIS와 DSIS로부터 준비된 DPIS와 DSIS를 생성하는 수단(150), 상기 준비된 DPIS와 DSIS로부터 갱신된 변환을 생성하는 수단(161 내지 165), 및 상기 정합된 DPIS와 DSIS를 생성하기위해 준비된 대략 DPIS와 DSIS중 한 개에 상기 갱신된 변환을 적용하는 수단(167)을 포함하는 것을 특징으로 하는 장치(27).
  8. 제 7 항에 있어서, 상기 준비된 DPIS와 DSIS를 생성하는 상기 수단(150)은, 교차하는 상기 DPIS와 DSIS에 의해 표시되는 이미지의 영역을 위해 선택된 DPIS와 DSIS를 선택하는 수단(151, 152)을 포함하는 것을 특징으로 하는 장치(27).
  9. 제 8 항에 있어서, 상기 준비된 DPIS와 DSIS를 생성하는 상기 수단(150)은 상기 선택된 DPIS와 DSIS를 향상시키는 수단을 더 포함하는 것을 특징으로 하는 장치(27).
  10. 제 9 항에 있어서, 상기 선택된 DPIS와 상기 선택된 DSIS를 선택하는 상기 수단(151, 152)은, 소정의 규칙적인 형태를 가지며 상기 DPIS와 DSIS에 의해 표시되는 이미지의 일부 영역내의 DPIS와 DSIS를 선택하는 수단(152)을 더 포함하는 것을 특징으로 하는 장치(27).
  11. 제 7 항에 있어서, 상기 갱신된 변형을 생성하는 상기 수단(161 내지 165)은, 상기 준비된 DPIS와 DSIS로부터 이동흐름 성분(163)를 생성하는 수단, 상기 이동흐름 성분로부터 상기 갱신된 변환을 계산하는 계산 수단(164)을 포함하는 것을 특징으로 하는 장치(27).
  12. 제 11 항에 있어서, 이동흐름 성분을 생성하는 상기 수단(163)은, 이동흐름 변화도 구성요소를 생성하고, 상기 계산 수단(164)은 상기 갱신된 변환을 계산하기위해 로버스트 최적화를 적용하는 수단을 포함하는 것을 특징으로 하는 장치(27).
  13. 제 12 항에 있어서, 상기 갱신된 변환을 생성하는 상기 수단(161 내지 165)은, 상기 준비된 DPIS와 DSIS의 연속적인 향상 레벨의 해상도를 반복적으로 사용하는 상기 계산 수단과 이동흐름 변화도 구성요소를 생성하는 상기 수단(163)을 이용하는 수단(161, 162, 166)을 포함하는 것을 특징으로 하는 장치(27).
  14. 제 7 항에 있어서, 상기 갱신된 변환을 생성하는 상기 수단(161 내지 165)은 상기 갱신된 변환을 생성하기위해 상기 준비된 DPIS와 DSIS의 연속적인 향상 레벨의 해상도를 사용하는 수단(161, 162, 166)을 포함하는 것을 특징으로 하는 장치(27).
  15. 제 7 항에 있어서, 상기 갱신된 변환을 생성하는 상기 수단(161 내지 166)은 로버스트 이동흐름을 상기 준비된 DPIS와 DSIS에 적용하는 수단(164)을 포함하는 것을 특징으로 하는 장치(27).
  16. 제 15 항에 있어서, 로버스트 이동흐름을 상기 준비된 DPIS와 DSIS에 적용하는 상기 수단(164)은, 상기 DPIS와 DSIS의 연속적인 향상 레벨의 해상도에 로버스트 이동흐름을 적용하는 것을 특징으로 하는 장치(27).
  17. 제 1 항에 있어서, 상기 출력 수단(37, 140)은 상기 정합된 DPIS와 DSIS로부터 디스플레이(59)를 생성하는 디스플레이 수단을 포함하는 것을 특징으로 하는 장치(27).
  18. 제 1 항에 있어서, 상기 출력 수단(37, 140)은 상기 DPIS에 의해 표시되는 상기 이미지에서의 움직임을 추적하는 추적 수단(140)을 포함하는 것을 특징으로 하는 장치(27).
  19. 제 18 항에 있어서, 상기 출력 수단(37, 140)은, 상기 포털 이미지(53)를 생성하는 방사선에 대한 환자위치를 정하는 위치지정 수단(15'), 및 상기 추적 수단(140)에 의해 추적되는 움직임에 반응하여 상기 위치지정 수단(15')을 제어하는 수단을 더 포함하는 것을 특징으로 하는 장치(27).
  20. 제 18 항에 있어서, 상기 출력 수단(37, 140)은, 상기 추적 수단(140)에 의해 추적되는 움직임에 반응하여 상기 포털 이미지(53)를 발생시키는 방사선(11)의 생성을 제어하는 수단(9')을 제어하는 수단을 포함하는 것을 특징으로 하는 장치(27).
  21. 방사선요법 진단장비를 제어하기위해 포털 이미지를 정합시키는 장치에 있어서,
    연속적인 디지털 포털 이미지 신호(DPIS)의 세트를 생성하기위해 연속적인 포털 이미지를 디지털화하는 수단(29); 및
    연속적인 DPIS 세트사이의 움직임을 추적하는 추적 수단(140, 150, 160)을 포함하는 것을 특징으로 하는 장치(27).
  22. 제 21 항에 있어서, 상기 추적 수단(140, 150, 160)은, 상기 연속적인 DPIS 세트에 로버스트 이동흐름을 적용함으로써 연속적인 포털 이미지(53)사이의 갱신된 변환을 생성하는 수단(160), 및 상기 연속적인 DPIS 세트사이의 상기 움직임을 추적하기위해 상기 갱신된 변환을 사용하는 수단(140)을 포함하는 것을 특징으로 하는 장치(27).
  23. 제 22 항에 있어서, 상기 갱신된 변환을 생성하는 상기 수단(160)은, 상기 연속적인 DPIS 세트로부터 이동흐름 성분을 생성하는 수단, 및 상기 이동흐름 성분을 사용하여 연속적인 포털 이미지사이의 상기 갱신된 변환을 계산하는 수단(164)을 포함하는 것을 특징으로 하는 장치(27).
  24. 제 23 항에 있어서, 이동흐름 성분을 생성하는 상기 수단(163)은 이동흐름 변화도 구성요소를 생성하고, 상기 계산 수단(164)은 상기 갱시된 변환을 계산하기위해 로버스트 최적화를 적용하는 수단을 포함하는 것을 특징으로 하는 장치(27).
  25. 제 24 항에 있어서, 상기 갱신된 변환을 생성하는 상기 수단(160)은, 상기 연속적인 DPIS 세트의 연속적인 향상 레벨의 해상도를 반복적으로 사용하는 상기 계산 수단(164)과 이동흐름 변화도 구성요소를 생성하는 상기 수단(163)을 이용하는 수단(166)을 포함하는 것을 특징으로 하는 장치(27).
  26. X-레이 이미지(53)를 기준 이미지(45)와 자동적으로 정합시키는 장치(27)에 있어서,
    제 1 디지털 이미지 신호와 제 2 디지털 이미지 신호를 생성하기위해 각각 상기 X-레이 이미지(53)와 기준 이미지(45)를 디지털화하는 수단(29);
    정합된 디지털 이미지 신호를 생성하기위해 상기 이미지내에서 어떠한 물리적 차원의 특성 입력없이 상기 제 1 그리고 제 2 디지털 신호를 처리하는 처리 수단(31); 및
    상기 정합된 디지털 이미지 신호로부터 디스플레이(59)를 생성하는 디스플레이 수단(37)을 포함하는 것을 특징으로 하는 장치(27).
  27. 제 26 항에 있어서, 상기 처리 수단(31)은, 상기 제 1 그리고 제 2 디지털 이미지 신호로부터 대략 정렬된 디지털 이미지 신호를 생성하는 대략 정렬 수단(33), 로버스트 이동흐름을 사용하여 상기 X-레이와 기준 이미지의 영역을 중첩시키기위해 상기 대략 정렬된 디지털 이미지 신호사이의 변환을 생성하는 미세 정렬 수단(35, 160), 및 상기 정합된 디지털 이미지 신호를 생성하기위해 상기 대략 정렬된 디지털 이미지 신호중 한 개에 상기 변환을 적용하는 수단(167)을 포함하는 것을 특징으로 하는 장치(27).
  28. 제 27 항에 있어서, 상기 미세 정렬 수단(35)은, 유사한 동적 범위와 조명도를 갖는 준비된 대략 정렬 이미지 신호를 생성하기위해 상기 대략 정렬된 디지털 이미지 신호를 향상시키는 수단(154), 및 로버스트 이동흐름을 이용하여 상기 준비된 대략 정렬 디지털 이미지 신호사이의 상기 변환을 생성하는 수단(160)을 포함하는 것을 특징으로 하는 장치(27).
KR1019990703726A 1996-10-29 1997-10-28 기준 이미지와 x-레이 이미지와의 정합용 장치 KR20000052875A (ko)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US08/739,622 US5784431A (en) 1996-10-29 1996-10-29 Apparatus for matching X-ray images with reference images
US8/739,622 1996-10-29

Publications (1)

Publication Number Publication Date
KR20000052875A true KR20000052875A (ko) 2000-08-25

Family

ID=24973119

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1019990703726A KR20000052875A (ko) 1996-10-29 1997-10-28 기준 이미지와 x-레이 이미지와의 정합용 장치

Country Status (8)

Country Link
US (1) US5784431A (ko)
EP (1) EP0934573A1 (ko)
JP (1) JP2001503176A (ko)
KR (1) KR20000052875A (ko)
CN (1) CN1235684A (ko)
AU (1) AU5002997A (ko)
CA (1) CA2270357A1 (ko)
WO (1) WO1998019272A1 (ko)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101440826B1 (ko) * 2013-02-05 2014-09-23 광주과학기술원 이미지 정합 방법 및 이미지 정합 장치
KR102110136B1 (ko) * 2018-11-07 2020-05-13 한국항공대학교산학협력단 특징점 누적 기반 호모그래피를 이용한 고정형 비디오 스티칭 장치 및 방법

Families Citing this family (116)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1998002091A1 (en) * 1996-07-11 1998-01-22 The Board Of Trustees Of The Leland Stanford Junior University High-speed inter-modality image registration via iterative feature matching
US6597818B2 (en) 1997-05-09 2003-07-22 Sarnoff Corporation Method and apparatus for performing geo-spatial registration of imagery
US6512857B1 (en) * 1997-05-09 2003-01-28 Sarnoff Corporation Method and apparatus for performing geo-spatial registration
GB2360684B (en) * 1997-11-15 2001-12-19 Elekta Ab Analysis of radiographic images
GB9724110D0 (en) * 1997-11-15 1998-01-14 Elekta Ab Analysis of radiographic images
US8788020B2 (en) * 1998-10-23 2014-07-22 Varian Medical Systems, Inc. Method and system for radiation application
US7158610B2 (en) * 2003-09-05 2007-01-02 Varian Medical Systems Technologies, Inc. Systems and methods for processing x-ray images
US6621889B1 (en) 1998-10-23 2003-09-16 Varian Medical Systems, Inc. Method and system for predictive physiological gating of radiation therapy
US6980679B2 (en) * 1998-10-23 2005-12-27 Varian Medical System Technologies, Inc. Method and system for monitoring breathing activity of a subject
US6937696B1 (en) 1998-10-23 2005-08-30 Varian Medical Systems Technologies, Inc. Method and system for predictive physiological gating
US6973202B2 (en) * 1998-10-23 2005-12-06 Varian Medical Systems Technologies, Inc. Single-camera tracking of an object
US6279579B1 (en) * 1998-10-23 2001-08-28 Varian Medical Systems, Inc. Method and system for positioning patients for medical treatment procedures
WO2000024467A1 (en) 1998-10-23 2000-05-04 Varian Medical Systems, Inc. Method and system for physiological gating of radiation therapy
DE69943026D1 (de) * 1998-12-16 2011-01-20 Michael I Miller Verfahren und gerät zur verarbeitung von bildern mit zielobjekte repräsentierenden regionen
US6754374B1 (en) 1998-12-16 2004-06-22 Surgical Navigation Technologies, Inc. Method and apparatus for processing images with regions representing target objects
US6694057B1 (en) 1999-01-27 2004-02-17 Washington University Method and apparatus for processing images with curves
US6778850B1 (en) * 1999-03-16 2004-08-17 Accuray, Inc. Frameless radiosurgery treatment system and method
US6501981B1 (en) * 1999-03-16 2002-12-31 Accuray, Inc. Apparatus and method for compensating for respiratory and patient motions during treatment
US6587601B1 (en) 1999-06-29 2003-07-01 Sarnoff Corporation Method and apparatus for performing geo-spatial registration using a Euclidean representation
WO2001095109A2 (en) * 2000-06-02 2001-12-13 Koninklijke Philips Electronics N.V. Method and apparatus for merging images
US6909792B1 (en) * 2000-06-23 2005-06-21 Litton Systems, Inc. Historical comparison of breast tissue by image processing
KR100453509B1 (ko) * 2000-08-25 2004-10-20 (주)바이오메디시스 치과용 디지털 엑스레이 이미지의 메모 작성 및 그 전자우편 전송 방법과 그 프로그램이 저장된 기록매체
US6422750B1 (en) * 2000-12-22 2002-07-23 Ge Medical Systems Global Technology Company, Llc Digital x-ray imager alignment method
US7769430B2 (en) * 2001-06-26 2010-08-03 Varian Medical Systems, Inc. Patient visual instruction techniques for synchronizing breathing with a medical procedure
US6810107B2 (en) * 2001-11-02 2004-10-26 Siemens Medical Solutions Usa, Inc. System and method for measuring beam quality and dosimetry using electronic portal imaging
US6810108B2 (en) 2001-11-02 2004-10-26 Siemens Medical Solutions Usa, Inc. System and method for positioning an electronic portal imaging device
US6839404B2 (en) * 2001-11-02 2005-01-04 Siemens Medical Solutions Usa, Inc. System and method for positioning an electric portal imaging device
US20030086596A1 (en) * 2001-11-07 2003-05-08 Medical Metrics, Inc. Method, computer software, and system for tracking, stabilizing, and reporting motion between vertebrae
JP4430937B2 (ja) * 2001-12-07 2010-03-10 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ ノイズの多い画像中で構造を空間的に強調する医用ビューイングシステム及び方法
DE10160611A1 (de) * 2001-12-11 2003-06-26 Siemens Ag Bildgebende medizinische Untersuchungseinrichtung
ES2217210T3 (es) * 2002-02-22 2004-11-01 Brainlab Ag Instrumento de calibrado de reducida altura.
ATE357190T1 (de) 2002-03-27 2007-04-15 Brainlab Ag Medizinische navigation bzw. prä-operative behandlungsplanung mit unterstützung durch generische patientendaten
EP1348394B1 (de) 2002-03-27 2006-02-22 BrainLAB AG Planungs- bzw. Navigationsunterstützung durch generische und erfasste Patientendaten mit zweidimensionaler Anpassung
DK1519681T3 (da) * 2002-07-09 2007-04-10 Aecc Entpr Ltd Fremgangsmåde til billeddannelse af den relative bevægelse af skeletstykker
US20040068182A1 (en) * 2002-09-18 2004-04-08 Misra Satrajit Chandra Digitally reconstruced portal image and radiation therapy workflow incorporating the same
EP1571999B1 (en) * 2002-12-11 2007-11-28 Philips Intellectual Property & Standards GmbH Method of tomographic imaging
US7412029B2 (en) * 2003-06-25 2008-08-12 Varian Medical Systems Technologies, Inc. Treatment planning, simulation, and verification system
CA2473963A1 (en) * 2003-07-14 2005-01-14 Sunnybrook And Women's College Health Sciences Centre Optical image-based position tracking for magnetic resonance imaging
US7873403B2 (en) * 2003-07-15 2011-01-18 Brainlab Ag Method and device for determining a three-dimensional form of a body from two-dimensional projection images
US20050053267A1 (en) * 2003-09-05 2005-03-10 Varian Medical Systems Technologies, Inc. Systems and methods for tracking moving targets and monitoring object positions
US8571639B2 (en) * 2003-09-05 2013-10-29 Varian Medical Systems, Inc. Systems and methods for gating medical procedures
FI20031835A (fi) * 2003-12-15 2005-06-16 Instrumentarium Corp Menetelmä ja järjestelmä referenssimerkin paikantamiseksi digitaalisista projektiokuvista
WO2005081178A1 (en) * 2004-02-17 2005-09-01 Yeda Research & Development Co., Ltd. Method and apparatus for matching portions of input images
JP2006014931A (ja) * 2004-07-01 2006-01-19 Fuji Photo Film Co Ltd 画像読影支援システム、および、画像位置合わせ装置、画像出力装置、ならびにプログラム
EP1790203B1 (en) 2004-07-21 2015-12-30 Mevion Medical Systems, Inc. A programmable radio frequency waveform generator for a synchrocyclotron
US20060074305A1 (en) * 2004-09-30 2006-04-06 Varian Medical Systems Technologies, Inc. Patient multimedia display
DE102005016256B3 (de) * 2005-04-08 2006-06-08 Siemens Ag Verfahren zum Darstellen präoperativ aufgenommener dreidimensionaler Bilddaten bei Aufnahme zweidimensionaler Röntgenbilder
EP1874388A4 (en) 2005-04-13 2011-09-07 Univ Maryland TECHNIQUES FOR COMPENSATING THE MOVEMENT OF A TREATMENT TARGET IN A PATIENT
US9119541B2 (en) * 2005-08-30 2015-09-01 Varian Medical Systems, Inc. Eyewear for patient prompting
EP2389978B1 (en) 2005-11-18 2019-03-13 Mevion Medical Systems, Inc. Charged particle radiation therapy
US7564950B2 (en) * 2006-02-17 2009-07-21 Siemens Medical Solutions Usa, Inc. Multi-leaf collimator based field size clipping for automatic adaptation to allowed image area
US8676293B2 (en) * 2006-04-13 2014-03-18 Aecc Enterprises Ltd. Devices, systems and methods for measuring and evaluating the motion and function of joint structures and associated muscles, determining suitability for orthopedic intervention, and evaluating efficacy of orthopedic intervention
WO2007130313A2 (en) * 2006-05-02 2007-11-15 Telesis Technologies, Inc. Laser safety system
CN101453950A (zh) * 2006-06-01 2009-06-10 皇家飞利浦电子股份有限公司 分层运动估计
EP1868157A1 (en) 2006-06-14 2007-12-19 BrainLAB AG Shape reconstruction using X-ray images
US7848592B2 (en) * 2006-07-31 2010-12-07 Carestream Health, Inc. Image fusion for radiation therapy
US8160364B2 (en) * 2007-02-16 2012-04-17 Raytheon Company System and method for image registration based on variable region of interest
WO2008139374A1 (en) * 2007-05-11 2008-11-20 Philips Intellectual Property & Standards Gmbh Method for planning 2d x-ray examinations
US7953247B2 (en) * 2007-05-21 2011-05-31 Snap-On Incorporated Method and apparatus for wheel alignment
WO2009012577A1 (en) * 2007-07-20 2009-01-29 Resonant Medical Inc. Methods and systems for compensating for changes in anatomy of radiotherapy patients
EP3281671A1 (en) 2007-09-13 2018-02-14 ProCure Treatment Centers, Inc. Patient positioner system
EP2197547B1 (en) * 2007-09-13 2014-06-11 Toby D. Henderson Imaging positioning system having robotically positioned d-arm
US7933380B2 (en) * 2007-09-28 2011-04-26 Varian Medical Systems International Ag Radiation systems and methods using deformable image registration
US20090099481A1 (en) 2007-10-10 2009-04-16 Adam Deitz Devices, Systems and Methods for Measuring and Evaluating the Motion and Function of Joints and Associated Muscles
US8581523B2 (en) 2007-11-30 2013-11-12 Mevion Medical Systems, Inc. Interrupted particle source
US8933650B2 (en) 2007-11-30 2015-01-13 Mevion Medical Systems, Inc. Matching a resonant frequency of a resonant cavity to a frequency of an input voltage
JP2011506032A (ja) * 2007-12-18 2011-03-03 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 一貫性指標に基づく画像レジストレーション
US8825136B2 (en) * 2008-03-14 2014-09-02 Baylor Research Institute System and method for pre-planning a radiation treatment
US8170319B2 (en) * 2008-09-05 2012-05-01 Siemens Medical Solutions Usa, Inc. Motion detection by direct imaging during radiotherapy
US10667727B2 (en) * 2008-09-05 2020-06-02 Varian Medical Systems, Inc. Systems and methods for determining a state of a patient
US20100061596A1 (en) * 2008-09-05 2010-03-11 Varian Medical Systems Technologies, Inc. Video-Based Breathing Monitoring Without Fiducial Tracking
US8394007B2 (en) 2008-10-31 2013-03-12 Toby D Henderson Inclined beamline motion mechanism
JP2010183976A (ja) * 2009-02-10 2010-08-26 Mitsubishi Heavy Ind Ltd 放射線治療装置制御装置および放射線照射方法
US8437577B2 (en) * 2009-03-05 2013-05-07 Tektronix, Inc. Methods and systems for image registration
US9138163B2 (en) 2009-09-25 2015-09-22 Ortho Kinematics, Inc. Systems and devices for an integrated imaging system with real-time feedback loop and methods therefor
US8758263B1 (en) 2009-10-31 2014-06-24 Voxel Rad, Ltd. Systems and methods for frameless image-guided biopsy and therapeutic intervention
US10657613B2 (en) * 2010-06-17 2020-05-19 Koninklijke Philips N.V. Identity matching of patient records
WO2012012664A2 (en) * 2010-07-21 2012-01-26 Moehrle Armin E Image reporting method
EP3649937A1 (en) 2010-12-13 2020-05-13 Statera Spine, Inc. Methods, systems and devices for clinical data reporting and surgical navigation
US8837791B2 (en) 2010-12-22 2014-09-16 Kabushiki Kaisha Toshiba Feature location method and system
US9235892B2 (en) 2011-03-31 2016-01-12 Denise De Andrade Castro Method and device for comparing radiographic images
JP5693388B2 (ja) * 2011-06-10 2015-04-01 三菱電機株式会社 画像照合装置、患者位置決め装置及び画像照合方法
US9533173B2 (en) * 2012-02-21 2017-01-03 Koninklijke Philips N.V. Adaptive radiotherapy with spectral tissue of interest imaging and tracking
CN102697561A (zh) * 2012-05-17 2012-10-03 深圳市一体医疗科技股份有限公司 一种无创体外固定标记点肿瘤定位系统及定位方法
CN102670234B (zh) * 2012-05-17 2013-11-20 西安一体医疗科技有限公司 一种伽玛辐射摆位验证装置及方法
EP2901820B1 (en) 2012-09-28 2021-02-17 Mevion Medical Systems, Inc. Focusing a particle beam using magnetic field flutter
US9545528B2 (en) 2012-09-28 2017-01-17 Mevion Medical Systems, Inc. Controlling particle therapy
TW201424466A (zh) 2012-09-28 2014-06-16 Mevion Medical Systems Inc 磁場再生器
WO2014052718A2 (en) 2012-09-28 2014-04-03 Mevion Medical Systems, Inc. Focusing a particle beam
US10254739B2 (en) 2012-09-28 2019-04-09 Mevion Medical Systems, Inc. Coil positioning system
WO2014052721A1 (en) 2012-09-28 2014-04-03 Mevion Medical Systems, Inc. Control system for a particle accelerator
JP6523957B2 (ja) 2012-09-28 2019-06-05 メビオン・メディカル・システムズ・インコーポレーテッド 磁場を変更するための磁性シム
TW201424467A (zh) 2012-09-28 2014-06-16 Mevion Medical Systems Inc 一粒子束之強度控制
EP3581242B1 (en) 2012-09-28 2022-04-06 Mevion Medical Systems, Inc. Adjusting energy of a particle beam
US8791656B1 (en) 2013-05-31 2014-07-29 Mevion Medical Systems, Inc. Active return system
US9730308B2 (en) 2013-06-12 2017-08-08 Mevion Medical Systems, Inc. Particle accelerator that produces charged particles having variable energies
WO2015048468A1 (en) 2013-09-27 2015-04-02 Mevion Medical Systems, Inc. Particle beam scanning
WO2015085252A1 (en) * 2013-12-06 2015-06-11 Sonitrack Systems, Inc. Radiotherapy dose assessment and adaptation using online imaging
US9962560B2 (en) 2013-12-20 2018-05-08 Mevion Medical Systems, Inc. Collimator and energy degrader
US10675487B2 (en) 2013-12-20 2020-06-09 Mevion Medical Systems, Inc. Energy degrader enabling high-speed energy switching
US9661736B2 (en) 2014-02-20 2017-05-23 Mevion Medical Systems, Inc. Scanning system for a particle therapy system
US9950194B2 (en) 2014-09-09 2018-04-24 Mevion Medical Systems, Inc. Patient positioning system
EP3034000A1 (en) * 2014-12-16 2016-06-22 Agfa Healthcare Motion correction method in dual energy radiography
US20160354161A1 (en) 2015-06-05 2016-12-08 Ortho Kinematics, Inc. Methods for data processing for intra-operative navigation systems
CA3021538C (en) * 2015-06-11 2023-09-26 University Of Pittsburgh-Of The Commonwealth System Of Higher Education Systems and methods for finding regions of interest in hematoxylin and eosin (h&e) stained tissue images and quantifying intratumor cellular spatial heterogeneity in multiplexed/hyperplexed fluorescence tissue images
US10786689B2 (en) 2015-11-10 2020-09-29 Mevion Medical Systems, Inc. Adaptive aperture
US10925147B2 (en) 2016-07-08 2021-02-16 Mevion Medical Systems, Inc. Treatment planning
US10881466B2 (en) 2016-08-29 2021-01-05 Covidien Lp Systems, methods, and computer-readable media of providing distance, orientation feedback and motion compensation while navigating in 3D
US11103730B2 (en) 2017-02-23 2021-08-31 Mevion Medical Systems, Inc. Automated treatment in particle therapy
US11478662B2 (en) 2017-04-05 2022-10-25 Accuray Incorporated Sequential monoscopic tracking
US10653892B2 (en) 2017-06-30 2020-05-19 Mevion Medical Systems, Inc. Configurable collimator controlled using linear motors
US11471217B2 (en) 2017-12-11 2022-10-18 Covidien Lp Systems, methods, and computer-readable media for improved predictive modeling and navigation
CN109656502A (zh) * 2018-11-30 2019-04-19 贵州电网有限责任公司 Gis设备数字射线图像双屏灰阶度自动对比诊断方法及装置
US10806339B2 (en) 2018-12-12 2020-10-20 Voxel Rad, Ltd. Systems and methods for treating cancer using brachytherapy
WO2020185543A1 (en) 2019-03-08 2020-09-17 Mevion Medical Systems, Inc. Collimator and energy degrader for a particle therapy system
CN111161297B (zh) * 2019-12-31 2023-06-16 上海联影医疗科技股份有限公司 限束器边缘的确定方法、装置及x射线系统

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5099846A (en) * 1988-12-23 1992-03-31 Hardy Tyrone L Method and apparatus for video presentation from a variety of scanner imaging sources
US4995068A (en) * 1989-10-02 1991-02-19 S&S Inficon, Inc. Radiation therapy imaging apparatus
DE4207632C2 (de) * 1992-03-11 1995-07-20 Bodenseewerk Geraetetech Vorrichtung und Verfahren zur Positionierung eines Körperteils für Behandlungszwecke

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101440826B1 (ko) * 2013-02-05 2014-09-23 광주과학기술원 이미지 정합 방법 및 이미지 정합 장치
KR102110136B1 (ko) * 2018-11-07 2020-05-13 한국항공대학교산학협력단 특징점 누적 기반 호모그래피를 이용한 고정형 비디오 스티칭 장치 및 방법

Also Published As

Publication number Publication date
AU5002997A (en) 1998-05-22
CA2270357A1 (en) 1998-05-07
JP2001503176A (ja) 2001-03-06
CN1235684A (zh) 1999-11-17
WO1998019272A1 (en) 1998-05-07
EP0934573A1 (en) 1999-08-11
US5784431A (en) 1998-07-21

Similar Documents

Publication Publication Date Title
KR20000052875A (ko) 기준 이미지와 x-레이 이미지와의 정합용 장치
Gilhuijs et al. Automatic on‐line inspection of patient setup in radiation therapy using digital portal images
Penney et al. A comparison of similarity measures for use in 2-D-3-D medical image registration
US6415171B1 (en) System and method for fusing three-dimensional shape data on distorted images without correcting for distortion
US7480399B2 (en) Apparatus and method for determining measure of similarity between images
CN104025119B (zh) 用于手术和介入性医疗过程中的成像系统和方法
JP4271941B2 (ja) 患者の断層撮影投影画像を増強するための方法
EP0840253B1 (en) Methods and apparatus for digital subtraction angiography
US6301495B1 (en) System and method for intra-operative, image-based, interactive verification of a pre-operative surgical plan
EP2992819B1 (en) Weighted surface-to-surface mapping
EP2823463B1 (en) Method and system to assist 2d-3d image registration
US20100160836A1 (en) Determination of indicator body parts and pre-indicator trajectories
Hristov et al. A grey‐level image alignment algorithm for registration of portal images and digitally reconstructed radiographs
GB2418495A (en) Image processing system for use with a patient postioning device
JP2006247268A (ja) 患者位置決めシステム及び患者位置決め方法
JP2023040242A (ja) 医療処置中のエフェクタの画像定位のためのシステム及び方法
Malsch et al. An enhanced block matching algorithm for fast elastic registration in adaptive radiotherapy
Penney Registration of tomographic images to X-ray projections for use in image guided interventions
CN103976751A (zh) 自动视野尺寸计算约束
Aouadi et al. Accurate and precise 2D–3D registration based on X-ray intensity
Bögel et al. Respiratory motion compensation using diaphragm tracking for cone-beam C-arm CT: A simulation and a phantom study
Cai et al. A simple algorithm for planar image registration in radiation therapy
Lei et al. Deep learning‐based fast volumetric imaging using kV and MV projection images for lung cancer radiotherapy: A feasibility study
WO2019198394A1 (ja) 医用画像処理装置、医用画像処理方法、およびプログラム
Weese et al. 2D/3D registration and motion tracking for surgical interventions

Legal Events

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