KR20080025055A - 고속의 변형 가능한 포인트에 기초한 영상화에 관한 포인트하위선택 - Google Patents

고속의 변형 가능한 포인트에 기초한 영상화에 관한 포인트하위선택 Download PDF

Info

Publication number
KR20080025055A
KR20080025055A KR1020077028448A KR20077028448A KR20080025055A KR 20080025055 A KR20080025055 A KR 20080025055A KR 1020077028448 A KR1020077028448 A KR 1020077028448A KR 20077028448 A KR20077028448 A KR 20077028448A KR 20080025055 A KR20080025055 A KR 20080025055A
Authority
KR
South Korea
Prior art keywords
vertices
vertex
imaged object
pair
distance
Prior art date
Application number
KR1020077028448A
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 KR20080025055A publication Critical patent/KR20080025055A/ko

Links

Images

Classifications

    • 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
    • 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16Measuring radiation intensity
    • G01T1/24Measuring radiation intensity with semiconductor detectors

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • Radiology & Medical Imaging (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Radiation-Therapy Devices (AREA)

Abstract

영상화된 물체에 대한 변형 가능한 등록을 수행하기 위해 정점을 선택하는 방법이 제공된다. 선택된 정점은 대응하는 쌍을 형성하고 각 쌍은 제 1 영상화된 물체로부터의 정점과 제 2 영상화된 물체로부터의 정점을 포함한다. 대응하는 정점 쌍은 대응하는 정점 쌍을 구성하는 정점들 사이의 거리의 순서대로 분류된다. 가장 큰 거리를 지닌 대응하는 정점 쌍이 맨 위의 우선 순위가 주어진다. 선택된 대응하는 정점 쌍으로부터의 선택된 거리 내에 있는 대응하는 정점 쌍이 버려진다. 이러한 식으로, 영상화된 물체의 변형 가능한 등록을 위해 사용된 정점 쌍의 개수가 감소하고 따라서 임상적으로 받아들일 수 있는 처리 시간을 허용한다.

Description

고속의 변형 가능한 포인트에 기초한 영상화에 관한 포인트 하위선택{POINT SUBSELECTION FOR FAST DEFORMABLE POINT-BASED IMAGING}
방사선 치료는 X-선 방사선과 같은 방사선으로 암 종양과 같은 질병을 치료하는 방법이다. 병든 세포 조직을 방사선으로 치료하는 과정에서, 일부 건강한 세포 조직 또한 방사선에 노출된다. 건강한 세포 조직이 방사선에 노출되는 것은 치료 관련 합병증을 야기할 수 있다. 이와 같이, 방사선이 주로 병든 세포 조직에 조사되고 그 주변의 건강한 세포 조직에는 최소로 조사되도록, 정확하고 정밀하게 병든 구역의 윤곽을 그리는 것이 바람직하다.
치료된 구역의 정확하고 정밀한 윤곽선(계획 목표 분량, 즉 PTV: planning target volume)은 세분화된 치료 동안 목표물의 움직임을 통합시킨다. 움직임은 환자의 물리적인 움직임(셋업 에러)이나 심장, 호흡기 및 소화기 계통 또는 치료 반응과 같은 생리적 기능에 의해 야기된, 병든 세포 조직을 포함하는 내부 세포 조직의 움직임과 변형일 수 있다. 종래의 치료 계획에서는, PTV가 환자 모집단의 통계치에 기초하여, 너무 크거나 부정확한 목표 영역을 초래한다. 환자 특정 운동을 평가하기 위해, 병든 조직 세포와 그 둘레 기관의 기하학적인 변경의 3차원 설명을 얻기 위해 일정 시간에 걸쳐 일련의 영상이 취해진다. 시간적 샘플링(temporal sampling)은, 예컨대 4차원의 게이팅된(gated) 영상화(imaging)를 사용하는 호흡 움직임을 감시하기 위한 초단위 또는 4차원 게이팅된 영상화 기술을 사용하는 주단위 영상화와 같은 날(day)과 주(week) 단위 또는 이들의 조합일 수 있다. 방사선 치료에서 3차원 영상의 시간 샘플을 통합하는 것을, 일반적으로 영상-안내된 방사선 치료(IGRT: image-guided radiotherapy) 또는 적응성 방사선치료라고 부른다. 그러한 조정은 더 정밀한 목표 구역에 방사선을 적용하는 것을 허용한다.
치료하는 동안 발생하는 변형 과정의 양적인 측정이, IGRT에서의 많은 응용에 있어 필수적이다. 투약 분량 히스토그램(DVH: dose volume histogram)이 치료 품질을 평가하는 중요한 수단으로, 이는 목표 영역과 위험한 장기에서의 투약 분포를 통계학적으로 설명한다. DVH는 종래에는 단일 3차원 영상에 기초하고 있었다. 4차원의 투약 히스토그램을 모으기 위해서는, 목표 구역의 분량 요소와 그 둘레 장기 사이의 공간 대응이 모든 영상에서 계산될 필요가 있다. 이는 변형 가능한 등록에 의해 행해질 수 있다. 창자 가스(gas)와 대변과 같은 운동 및 대량의 운동에 의해 야기된 변화로 인해, 복셀(voxel)에 기초한 등록 방법이 통상 비효율적인데, 이는 그것이 영상에서의 대응하는 포인트들을 매칭하는 것이 너무 어렵기 때문이다. 대안적으로, 기하학적 대응(예컨대, 포인트 또는 표면 사이의)에 기초한 등록 방법이 사용될 수 있다. 모델에 기초한 분할 방법은 해당 물체의 표면으로 3각형 표면 메시를 자동으로 변형시킴으로써, 관심 구역에 있는 구조체를 분할하기 위해 성공적으로 사용되었다. 만약 이러한 기술이 모든 영상에서의 모든 구조체에 적용된다면, 적응된 메시의 정점들은 기준 물체 표면과 영상화된 물체 표면 사이의 대응하 는 경계표(landmark)를 한정한다. 이러한 경계표는 정렬되고, 그 영상에서의 모든 복셀에 관한 변형 필드가 변형 가능한 포인트에 기초한 등록 방법으로 추정될 수 있어, 그 영상에서의 각 복셀의 운동을 설명하는 분량 측정 변형 필드를 초래한다.
변형 가능한 모델은 생물학적 구조체의 상당한 변이성을 수용할 수 있지만, 그러한 모델은 복잡하다. 통상, 해부학적 물체에 관해 한정된 정점의 양은 수천 개에 이른다. 수천 개의 정점 또는 경계표에 기초한 포인트-기초의 등록 기술을 사용하여 분량 측정 변형 필드를 계산하는 것은, 상당한 처리 시간을 필요로 하고, 임상 환경에서는 실행 가능하지 않다. 이와 같이, 임상적으로 받아들일 수 있는 시간 기간 내에서 해부학적 물체에 관한 추정된 분량 측정 변형을 제공할 수 있는 변형 절차에 관한 필요성이 존재한다.
본 발명은 임상적으로 받아들일 수 있는 시간 기간 내에서 분량 측정 물체를 생성하기 위한 변형 모델을 지향한다. 이러한 변형 모델은 변형 가능한 등록에서 사용될 감소된 개수의 정점을 체계적으로 선택한다. 사용된 감소된 개수의 정점은 받아들일 수 있는 처리 시간을 허용한다. 사용될 정점의 체계적인 선택은 영상화된 물체의 전체 표면에 걸쳐 선택된 정점의 분포를 제공하고, 가장 가까운 정점에 의해 야기된 지형적(topographical) 오차를 감소시킨다.
일 실시예에서, 선택 거리에 기초하여 대응하는 정점 쌍의 개수를 감소시키는 방법이 제공된다. 이 방법은 정점 쌍의 정점들 사이의 거리에 기초하여 내림 차순으로 대응하는 정점 쌍을 분류한다. 가장 큰 거리를 지닌 대응하는 정점 쌍에 가장 높은 우선 순위가 주어진다. 선택된 거리 내에 있는 대응하는 정점 쌍이 버려지고, 이를 통해 변형 가능한 등록에서 사용된 대응하는 정점 쌍의 개수를 감소시킨다. 일부 실시예에서, 선택된 거리는 선택된 대응하는 정점 쌍의 정점들 사이의 거리이다. 다른 실시예에서는, 선택된 거리가 미리 결정된다. 일부 실시예에서는, 그 거리가 변화될 수 있다.
본 명세서에 통합되고 본 명세서의 일부를 구성하는 첨부한 도면에서, 위에서 주어진 본 발명의 일반적인 설명과 아래에 주어진 상세한 설명과 함께 본 발명의 원리를 예시하는 역할을 하는 본 발명의 실시예가 예시되어 있다. 당업자라면 이들 예시적인 실시예가 본 발명을 제한하려는 것이 아니라 본 발명의 원리를 통합하는 예들을 제공한다는 점을 깨닫게 된다.
도 1은 정점 서브세트를 지닌 3각형 표면 모델로서의 제 1 영상화된 방광(bladder)을 도시하는 도면.
도 2는 정점 서브세트를 지닌 3각형 표면 모델로서의 제 2 영상화된 방광을 도시하는 도면.
도 3은 정정 서브세트에 기초한 제 1 영상화된 방광과 제 2 영상화된 방광 사이의 대응을 도시하는 도면.
도 4는 물체의 3각형 메시 정점들의 선택된 서브세트에 기초한 영상화된 물체의 분량 측정 변형을 계산하기 위한 예시적인 방법 흐름도.
도 5는 대응하는 정점 쌍의 선택된 서브세트를 결정하는 예시적인 방법을 도 시하는 도면.
본 명세서에서 개시된 변형 모델은 임상적으로 받아들일 수 있는 시간 기간 내에서 분량 측정 영상 몰체의 생성을 제공하는 것이다. 이러한 변형 모델은 운동을 설명하는 영상화된 물체의 재구성을 허용한다. 변형 모델은 관심 목표의 더 정확한 윤곽 형성을 허용하여, 더 정밀한 방사선 치료법을 허용한다. 비록 본 개시물이 영상 안내된 방사선 치료에 본 발명을 적용하는 가능성을 논의하지만, 당업자라면 개시된 변형 모델이 다른 영상화 절차나 방식, 특히 시간 간격에 걸쳐 관심 구역의 운동을 설명하는 것에 동등하게 적용 가능하다는 점을 알아야 한다.
개시된 변형 모델은 변형에서 사용된 정점의 개수를 감소시키는 방법을 제공하여, 임상 사용을 위한 적당한 처리 시간을 허용한다. 정점의 감소는 영상화 정확도의 감소를 일으키지 않고 처리 시간의 감소를 허용하도록 되어야 한다. 이러한 점에서, 본 명세서에서 개시된 변형 모델은, 선택된 정점이 영상화된 물체의 표면에 걸쳐 분포하도록 정점의 감소를 제공한다. 선택된 정점은 표면의 모든 부분을 나타내고, 가깝게 위치한 경계표의 탄성(elastic) 포인트에 기초한 등록에 의해 야기된 위상적인 위반을 실질적으로 감소시키도록 충분히 공간적으로 떨어져 배치된다.
본 발명의 일 실시예에서, 변형 모델은 도 4의 10에서 시작한다. 도 1과 도 2에 도시된 방광과 같은 대응하는 물체의 한 쌍의 3각형 표면 메시가 20에서 입력된다. 3각형 표면 메시의 각 정점의 위치는 도 4의 30에서 도시된 정점 인덱스에 기록된다. 비록 많은 경우에 물체들이 동일한 개수의 정점을 가지더라도, 2개의 물체는 이러한 방법이 효과적이 되도록 동일한 개수의 정점을 가질 필요는 없다. 35에서, 2개의 3각형 표면 메시의 각 정점은 매칭되어 대응하는 정점 쌍의 세트를 만들어낸다. 만약 2개의 메시가 동일한 토폴로지를 가진다면, 대응하는 정점 쌍은 정점의 인덱스에 의해 결정된다. 그렇지 않으면 대응하는 정점이, 예컨대 유클리드 거리 개념에서 가장 가까운 것으로 또는 정점 위치에서 수직인 표면을 따라 가장 가까운 것으로 또는 일부 다른 방법으로 결정될 수 있다.
정점의 각각의 대응하는 쌍 사이의 유클리드 거리는 38에서 계산되고, 이후 40에서 그 거리가 내림차순으로 분류된다. 이와 같이, 가장 큰 거리를 가지는 대응하는 정점 쌍에 가장 높은 우선 순위가 주어진다. 가장 큰 거리를 지닌 대응하는 정점 쌍은 45에서 선택된다. 선택된 대응하는 정점 쌍의 반경 내에 있는 임의의 대응하는 정점 쌍은 이후 48에서 버려진다. 이는 도 5를 고려함으로써 가장 잘 이해될 수 있고, 도 5는 소수의 대응하는 정점 쌍의 오버레이(overlay)를 나타내며, 이 경우 제 1 물체로부터의 정점은 실선으로 나타나고 제 2 물체로부터의 정점은 점선으로 나타나 있다. 포인트 A와 B가 가장 큰 거리를 지닌 대응하는 정점 쌍으로 결정된다면, 선택된 대응하는 정점 쌍 사이의 거리를 나타내는 반경(R)이 결정된다. 포인트(A) 둘레에서 반경(R)만큼 스쳐 지나가면, 포인트 C와 포인트 D가 반경 R 내에 있고, 포인트 E는 반경 R의 밖에 있다는 것을 알 수 있다. 단계(48)는 제 2 영상화된 물체와 대응하는 정점 쌍 목록으로부터 포인트 C, 포인트 D 및 그것들의 대응하는 정점을 버린다. 버려지는 포인트는 변형 가능한 등록에서 사용되지 않는다.
일단 선택된 정점 쌍의 반경 내에 있는 대응하는 정점 쌍이 버려지면, 본 발명의 방법은 임의의 추가적인 대응하는 정점 쌍이 정점 쌍 목록에 남아있는지 여부가 결정되는 단계(50)로 진행한다. 만약 남아있지 않다면, 즉 모든 대응하는 정점 쌍이 선택되거나 버려졌다면, 본 발명의 방법은 단계(60)로 진행하고, 변형 가능한 등록이 수행된다. 대응하는 정점 쌍 목록에 대응하는 정점 쌍이 남아있다면, 본 발명의 방법은 단계(45)로 되돌아가 가장 큰 거리를 지닌 다음 대응하는 정점 쌍을 결정하는데, 이러한 정점 쌍은 이전에 선택되거나 버려지지 않은 것이다. 이후 본 발명의 방법은 단계(48)로 진행하여, 새롭게 선택된 대응하는 정점 쌍의 반경 내의 대응하는 포인트가 버려진다. 이러한 과정은 모든 대응하는 정점 쌍이 선택되거나 버려질 때가지 반복된다.
이러한 과정의 결과, 통상 수천 개인 대응하는 정점 쌍의 개수는 대략 10부터 대략 100의 범위 또는 더 많은 범위를 가지는 인자에 의해 체계적으로 감소된다. 예컨대, 수천 개의 대응하는 정점 쌍을 가지는 해부학적 물체는 50 또는 100개의 대응하는 정점 쌍으로 감소될 수 있다. 선택된 정점 쌍들의 개수는, 물체의 크기와 영상 사이의 상대적인 기하학적인 변화량을 포함하는 다수의 인자에 따라 변한다. 감소된 개수의 정점은 영상화된 물체의 변형 가능한 등록에 관한 임상적으로 받아들일 수 있는 처리 시간을 허용한다. 본 발명의 방법은 또한 물체의 전체 표면에 걸쳐 퍼져나가는 정점을 선택하기 위해 제공한다. 이는 물체의 모든 부분이 변형 가능한 등록에서 나타나는 것을 보장하고, 위상 오차를 초래할 수 있는 가장 가까운 정점 쌍의 사용을 방지한다.
변형 가능한 등록 방법의 또 다른 실시예에서, 버려지는 대응하는 정점 쌍은 선택된 대응하는 정점 쌍의 반경에 의해 결정되지 않고, 대신 미리 결정된 거리에 기초한다. 일부 실시예에서, 미리 결정된 거리는 임상의에 의해 바뀔 수 있다. 일부 실시예에서, 임상의에게 거리를 입력할 것을 요구하는 단계가 추가될 수 있다. 이후 본 발명의 방법이 실행되고, 선택된 대응하는 정점 쌍의 개수가 결정된다. 이후 본 발명의 방법은 임상의에게 변형 가능한 등록을 수행하기 전에 선택된 대응하는 쌍의 개수가 바라는 범위 내에 있는지를 물을 수 있다. 선택된 대응하는 정점 쌍의 개수가 바라는 범위 내에 있지 않으면, 임상의는 미리 결정된 거리를 바꾸고, 본 발명의 방법이 새로운 개수의 선택된 대응하는 정점 쌍을 결정하는 것을 허용할 수 있다. 대안적으로, 선택된 정점 쌍의 개수는 임계 범위를 가질 수 있다. 선택된 대응하는 정점 쌍의 개수가 임계 범위 내에 있다면, 변형 가능한 등록이 수행된다. 선택된 대응하는 정점 쌍의 개수가 임계 범위 밖에 있다면, 선택된 대응하는 정점 쌍의 개수가 임계 범위 내에 있을 때까지 미리 결정된 거리가 바뀐다.
당업자라면, 위에서 설명한 방법이 어떤 영상에서의 단일 물체나 어떤 영상에서의 다수의 물체의 변형 가능한 등록을 수행하기 위해 구현될 수 있다는 사실을 알아야 한다. 다수의 물체가 그 영상 내에 있다면, 그러한 모든 물체의 정점은 각 물체마다 하나씩 대응하는 정점 구름(cloud)으로 결합된다. 본 발명의 방법은 이후, 각각의 대응하는 정점 구름에 관해 수행된다. 영상화된 물체가 동일한 개수의 정점을 가질 필요는 없다는 사실이 또한 주목되어야 한다. 정점의 양이 대응하는 물체에서 동일하지 않다면, 더 작은 개수의 정점을 지닌 물체에 기초한 대응이 한 정될 수 있다. 더 많은 개수의 정점을 지닌 물체는 여분의 정점을 가지게 되는데, 즉 이들은 버려진 대응하는 정점 쌍을 초래하지 않는다.
본 발명이 하나 이상의 바람직한 실시예를 참조하여 설명되었다. 분명히, 본 명세서를 읽고 이해하는 다른 당업자에게, 수정과 대안이 떠오르게 될 것이다. 모든 그러한 수정과 대안은 그것들이 첨부된 청구항 또는 그 등가물의 범주 내에 있는 한 포함되는 것으로 의도된다.
전술한 바와 같이, 본 발명은 방사선 치료에서 임상적으로 받아들일 수 있는 시간 기간 내에서 해부학적 물체에 관한 판단된 분량 측정 변형을 제공할 수 있는 변형 절차에 이용 가능하다.

Claims (19)

  1. 영상화된 물체에 대한 변형 가능한 등록을 수행하기 위해 정점(vertices)을 선택하는 방법으로서,
    제 1 영상화된 물체와 제 2 영상화된 물체 위의 복수의 정점의 위치를 결정하는 단계(20, 30),
    각각의 대응하는 정점 쌍이 제 1 영상화된 물체로부터의 정점과 제 2 영상화된 물체로부터의 정점을 포함하는, 대응하는 정점 쌍(A-B,C,D)의 세트를 결정하는 단계(35),
    각각의 대응하는 정점 쌍(A-B,C,D) 사이의 거리(R)를 결정하는 단계(38),
    각각의 대응하는 정점 쌍(A-B,C,D) 사이의 거리(R)에 기초한 대응하는 정점 쌍(A-B,C,D)의 세트를 분류하는 단계(40),
    가장 큰 거리를 지닌 대응하는 정점 쌍(A-B)을 선택하는 단계(45),
    선택된 정점 쌍(A-B)의 선택된 거리 내에 있는 대응하는 정점 쌍(C, D)을 버리는 단계(48), 및
    모든 대응하는 정점 쌍(A-B,C,D)이 선택되거나 버려질 때까지 선택하는 단계(45)와 버리는 단계(48)를 반복하는 단계를
    포함하는, 영상화된 물체에 대한 변형 가능한 등록을 수행하기 위해 정점을 선택하는 방법.
  2. 제 1항에 있어서, 상기 선택된 거리는 선택된 정점 쌍(A-B)의 정점 사이의 거리(R)인, 영상화된 물체에 대한 변형 가능한 등록을 수행하기 위해 정점을 선택하는 방법.
  3. 제 1항에 있어서, 상기 선택된 거리는 미리 결정된 값인, 영상화된 물체에 대한 변형 가능한 등록을 수행하기 위해 정점을 선택하는 방법.
  4. 제 1항에 있어서, 상기 선택된 거리는 변화될 수 있는, 영상화된 물체에 대한 변형 가능한 등록을 수행하기 위해 정점을 선택하는 방법.
  5. 제 1항에 있어서, 선택된 대응하는 정점 쌍(A-B,C,D)의 개수를 결정하는 단계를 더 포함하는, 영상화된 물체에 대한 변형 가능한 등록을 수행하기 위해 정점을 선택하는 방법.
  6. 제 5항에 있어서, 선택된 정점 쌍의 개수가 미리 결정된 임계 범위 내에 있지 않다면, 상기 선택된 거리를 수정하는 단계를 더 포함하는, 영상화된 물체에 대한 변형 가능한 등록을 수행하기 위해 정점을 선택하는 방법.
  7. 제 1항에 있어서, 상기 제 1 영상화된 물체에서의 정점의 개수는, 상기 제 2 영상화된 물체에서의 정점의 개수와 같지 않은, 영상화된 물체에 대한 변형 가능한 등록을 수행하기 위해 정점을 선택하는 방법.
  8. 제 1항에 있어서, 상기 제 1 영상화된 물체와 제 2 영상화된 물체 각각은 다수의 해부학적(anatomical) 물체로 이루어지는, 영상화된 물체에 대한 변형 가능한 등록을 수행하기 위해 정점을 선택하는 방법.
  9. 제 1항에 있어서, 상기 제 1 영상화된 물체와 제 2 영상화된 물체 위의 복수의 정점은 3각형 표면 메시(mesh)로부터 유도되는, 영상화된 물체에 대한 변형 가능한 등록을 수행하기 위해 정점을 선택하는 방법.
  10. 영상화된 물체에 대한 변형 가능한 등록을 수행하기 위해 정점을 선택하는 장치로서,
    제 1 영상화된 물체와 제 2 영상화된 물체 위의 복수의 정점의 위치를 결정하는 수단(20, 30),
    각각의 대응하는 정점 쌍이 제 1 영상화된 물체로부터의 정점과 제 2 영상화된 물체의 정점을 포함하는, 대응하는 정점 쌍(A-B,C,D)의 세트를 결정하는 수단(35),
    각각의 대응하는 정점 쌍(A-B,C,D) 사이의 거리(R)에 기초한 대응하는 정점 쌍(A-B,C-D)의 세트를 분류하는 수단(40),
    가장 큰 거리를 지닌 대응하는 정점 쌍(A-B)을 선택하는 수단(45),
    선택된 정점 쌍(A-B)의 선택된 거리 내에 있는 대응하는 정점 쌍(C, D)을 버리는 수단(48), 및
    다음으로 가장 큰 거리를 지닌 대응하는 정점 쌍을 선택하고, 새롭게 선택된 정점 쌍의 상기 선택된 거리 내에 있는 대응하는 정점 쌍을 버리며, 모든 대응하는 정점 쌍이 선택되거나 버려질 때까지, 대응하는 정점 쌍을 선택하는 것(45)과 버리는 것(48)의 반복(50)을 제공하는 수단을
    포함하는, 영상화된 물체에 대한 변형 가능한 등록을 수행하기 위해 정점을 선택하는 장치.
  11. 제 10항에 있어서, 상기 선택된 거리는 선택된 정점 쌍(A-B)의 정점 사이의 거리(R)인, 영상화된 물체에 대한 변형 가능한 등록을 수행하기 위해 정점을 선택하는 장치.
  12. 제 10항에 있어서, 상기 선택된 거리는 미리 결정된 값인, 영상화된 물체에 대한 변형 가능한 등록을 수행하기 위해 정점을 선택하는 장치.
  13. 제 10항에 있어서, 상기 선택된 거리는 변화될 수 있는, 영상화된 물체에 대한 변형 가능한 등록을 수행하기 위해 정점을 선택하는 장치.
  14. 제 10항에 있어서, 선택된 대응하는 정점 쌍(A-B,C,D)의 개수를 결정하는 수 단을 더 포함하는, 영상화된 물체에 대한 변형 가능한 등록을 수행하기 위해 정점을 선택하는 장치.
  15. 제 14항에 있어서, 선택된 정점 쌍의 개수가 미리 결정된 임계 범위 내에 있지 않다면, 상기 선택된 거리를 수정하는 수단을 더 포함하는, 영상화된 물체에 대한 변형 가능한 등록을 수행하기 위해 정점을 선택하는 장치.
  16. 제 10항에 있어서, 상기 제 1 영상화된 물체에서의 정점의 개수는, 상기 제 2 영상화된 물체에서의 정점의 개수와 같지 않은, 영상화된 물체에 대한 변형 가능한 등록을 수행하기 위해 정점을 선택하는 장치.
  17. 제 10항에 있어서, 상기 제 1 영상화된 물체와 제 2 영상화된 물체 각각은 다수의 해부학적 물체로 이루어지는, 영상화된 물체에 대한 변형 가능한 등록을 수행하기 위해 정점을 선택하는 장치.
  18. 제 1항에 있어서, 상기 제 1 영상화된 물체와 제 2 영상화된 물체 위의 복수의 정점은 3각형 표면 메시로부터 유도되는, 영상화된 물체에 대한 변형 가능한 등록을 수행하기 위해 정점을 선택하는 장치.
  19. 영상화된 물체에 대한 변형 가능한 등록을 수행하는 방법으로서,
    제 1 영상화된 물체와 제 2 영상화된 물체 위의 복수의 정점의 위치를 결정하는 단계(20, 30),
    각각의 대응하는 정점 쌍이 제 1 영상화된 물체로부터의 정점과 제 2 영상화된 물체로부터의 정점을 포함하는, 대응하는 정점 쌍(A-B,C,D)의 세트를 결정하는 단계(35),
    각각의 대응하는 정점 쌍(A-B,C,D) 사이의 거리를 결정하는 단계(38),
    각각의 대응하는 정점 쌍(A-B,C,D) 사이의 거리에 기초한 대응하는 정점 쌍(A-B,C,D)의 세트를 분류하는 단계(40),
    가장 큰 거리를 지닌 대응하는 정점 쌍(A-B)을 선택하는 단계(45),
    선택된 정점 쌍의 선택된 거리 내에 있는 대응하는 정점 쌍(C, D)을 버리는 단계(48),
    모든 대응하는 정점 쌍이 선택되거나 버려질 때까지 선택하는 단계(45)와 버리는 단계(48)를 반복하는 단계,
    영상화된 물체에 대한 변형 가능한 등록을 수행하기 위해, 선택된 대응하는 정점 쌍을 사용하는 단계(60)를
    포함하는, 영상화된 물체에 대한 변형 가능한 등록을 수행하는 방법.
KR1020077028448A 2005-06-08 2006-05-31 고속의 변형 가능한 포인트에 기초한 영상화에 관한 포인트하위선택 KR20080025055A (ko)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US59512205P 2005-06-08 2005-06-08
US60/595,122 2005-06-08

Publications (1)

Publication Number Publication Date
KR20080025055A true KR20080025055A (ko) 2008-03-19

Family

ID=37027053

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020077028448A KR20080025055A (ko) 2005-06-08 2006-05-31 고속의 변형 가능한 포인트에 기초한 영상화에 관한 포인트하위선택

Country Status (6)

Country Link
US (1) US8094895B2 (ko)
EP (1) EP1894163B1 (ko)
JP (1) JP5058985B2 (ko)
KR (1) KR20080025055A (ko)
CN (1) CN101194288B (ko)
WO (1) WO2006131848A2 (ko)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007133932A2 (en) 2006-05-11 2007-11-22 Koninklijke Philips Electronics, N.V. Deformable registration of images for image guided radiation therapy
WO2011073820A1 (en) * 2009-12-16 2011-06-23 Koninklijke Philips Electronics N.V. Use of collection of plans to develop new optimization objectives
EP2678826B1 (en) 2011-02-23 2018-09-19 Koninklijke Philips N.V. Automatic projection of landmarks to generate additional correspondences in image registration
EP2686831A1 (en) 2011-03-15 2014-01-22 Koninklijke Philips N.V. Correlated image mapping pointer
EP2721580B1 (en) 2011-06-16 2019-04-10 Koninklijke Philips N.V. Hybrid point-based registration
US9654173B2 (en) 2014-11-20 2017-05-16 At&T Intellectual Property I, L.P. Apparatus for powering a communication device and methods thereof
US10650940B2 (en) 2015-05-15 2020-05-12 At&T Intellectual Property I, L.P. Transmission medium having a conductive material and methods for use therewith
CN110547766B (zh) * 2019-08-22 2023-04-28 苏州佳世达光电有限公司 口扫机的操作方法

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5689577A (en) * 1994-10-14 1997-11-18 Picker International, Inc. Procedure for the simplification of triangular surface meshes for more efficient processing
US5886702A (en) 1996-10-16 1999-03-23 Real-Time Geometry Corporation System and method for computer modeling of 3D objects or surfaces by mesh constructions having optimal quality characteristics and dynamic resolution capabilities
US6106466A (en) 1997-04-24 2000-08-22 University Of Washington Automated delineation of heart contours from images using reconstruction-based modeling
US6301496B1 (en) * 1998-07-24 2001-10-09 Biosense, Inc. Vector mapping of three-dimensionally reconstructed intrabody organs and method of display
US6633686B1 (en) * 1998-11-05 2003-10-14 Washington University Method and apparatus for image registration using large deformation diffeomorphisms on a sphere
DE19916978C1 (de) * 1999-04-15 2001-04-26 Bock Orthopaed Ind Verfahren zum Vermessen eines Körperbereichs
JP4159226B2 (ja) * 2000-03-21 2008-10-01 住友重機械工業株式会社 患部トラッキング方法、装置、及び、これを用いた放射線照射方法、装置、並びに放射線治療装置
EP1433141A1 (en) * 2001-01-30 2004-06-30 Koninklijke Philips Electronics N.V. Image processing method for displaying an image sequence of a deformable 3-d object with indications of the object wall motion
JP2002298141A (ja) * 2001-03-29 2002-10-11 Nec Corp パターン照合装置とそのパターン照合方法、及びパターン照合プログラム
JP3903783B2 (ja) * 2001-12-14 2007-04-11 日本電気株式会社 顔メタデータ生成方法および装置、並びに顔類似度算出方法および装置
US6937751B2 (en) * 2003-07-30 2005-08-30 Radiological Imaging Technology, Inc. System and method for aligning images
JP4356831B2 (ja) * 2003-08-28 2009-11-04 株式会社Aze 複数画像間の非剛体レジストレーション方法
US7379598B2 (en) * 2003-08-29 2008-05-27 The Johns Hopkins University Distance sorting algorithm for matching patterns
FR2875146B1 (fr) * 2004-09-15 2007-09-07 Degremont Sa Installation de traitement d'effluents, et procede de clarification et de filtration utilisant cette installation

Also Published As

Publication number Publication date
US20090279739A1 (en) 2009-11-12
WO2006131848A3 (en) 2007-03-15
CN101194288A (zh) 2008-06-04
EP1894163B1 (en) 2015-09-30
CN101194288B (zh) 2010-06-09
WO2006131848A2 (en) 2006-12-14
JP5058985B2 (ja) 2012-10-24
JP2008545498A (ja) 2008-12-18
US8094895B2 (en) 2012-01-10
EP1894163A2 (en) 2008-03-05

Similar Documents

Publication Publication Date Title
US10300303B2 (en) Therapy control using motion prediction based on cyclic motion model
CN107072624B (zh) 用于自动治疗计划的***和方法
CN101443816B (zh) 用于图像引导下辐射治疗的图像可变形配准
RU2629235C2 (ru) Вычисление дозы на уровне сегментов пучка и отслеживание движения во времени для адаптивного планирования лечения
JP6732340B2 (ja) 4次元画像情報を構成するための磁気共鳴投影
RU2603606C2 (ru) Изучение дозиметрического воздействия движения на формирование адаптивных границ для конкретного пациента при планировании наружной дистанционной лучевой терапии
US10376714B2 (en) Method for improved surface tracking-based motion management and dynamic planning in adaptive external beam radiation therapy
US7596207B2 (en) Method of accounting for tumor motion in radiotherapy treatment
KR20080025055A (ko) 고속의 변형 가능한 포인트에 기초한 영상화에 관한 포인트하위선택
US9498645B2 (en) Dose deformation error calculation method and system
US9919163B2 (en) Methods, systems and computer readable storage media for determining optimal respiratory phase for treatment
McClelland Estimating internal respiratory motion from respiratory surrogate signals using correspondence models
Ertan et al. Assessment of the Reproducibility of Deep Inspiration Breath Hold Technique During Left-Sided Breast Cancer Radiotherapy

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