KR101946577B1 - Method and apparatus for generating 3d volume panorama - Google Patents
Method and apparatus for generating 3d volume panorama Download PDFInfo
- Publication number
- KR101946577B1 KR101946577B1 KR1020180078251A KR20180078251A KR101946577B1 KR 101946577 B1 KR101946577 B1 KR 101946577B1 KR 1020180078251 A KR1020180078251 A KR 1020180078251A KR 20180078251 A KR20180078251 A KR 20180078251A KR 101946577 B1 KR101946577 B1 KR 101946577B1
- Authority
- KR
- South Korea
- Prior art keywords
- volume
- image
- volume image
- images
- function
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 230000009466 transformation Effects 0.000 claims abstract description 143
- 230000006870 function Effects 0.000 claims description 453
- 238000006243 chemical reaction Methods 0.000 claims description 333
- 238000005457 optimization Methods 0.000 claims description 127
- 230000000877 morphologic effect Effects 0.000 claims description 43
- 239000000523 sample Substances 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 8
- 230000002194 synthesizing effect Effects 0.000 claims description 7
- 239000002131 composite material Substances 0.000 description 28
- 239000011159 matrix material Substances 0.000 description 20
- 239000013598 vector Substances 0.000 description 16
- 238000010586 diagram Methods 0.000 description 13
- 230000008569 process Effects 0.000 description 12
- 238000002604 ultrasonography Methods 0.000 description 8
- 238000002591 computed tomography Methods 0.000 description 4
- 210000003754 fetus Anatomy 0.000 description 4
- 238000002595 magnetic resonance imaging Methods 0.000 description 4
- 238000013519 translation Methods 0.000 description 4
- 238000003745 diagnosis Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 238000012285 ultrasound imaging Methods 0.000 description 3
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- FKOQWAUFKGFWLH-UHFFFAOYSA-M 3,6-bis[2-(1-methylpyridin-1-ium-4-yl)ethenyl]-9h-carbazole;diiodide Chemical compound [I-].[I-].C1=C[N+](C)=CC=C1C=CC1=CC=C(NC=2C3=CC(C=CC=4C=C[N+](C)=CC=4)=CC=2)C3=C1 FKOQWAUFKGFWLH-UHFFFAOYSA-M 0.000 description 1
- 230000004931 aggregating effect Effects 0.000 description 1
- 210000000601 blood cell Anatomy 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 210000002381 plasma Anatomy 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/08—Volume rendering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/41—Medical
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computer Graphics (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Image Processing (AREA)
Abstract
볼륨 파노라마 영상 생성 방법 및 장치에 따르면, 대상체 내부의 관찰 영역을 3차원적으로 나타내는 복수의 볼륨 영상들의 영상 데이터들을 입력받고, 입력된 볼륨 영상들 중 제1 볼륨 영상에 제1 볼륨 영상과 공통된 영역을 갖는 제2 볼륨 영상이 정합되도록 제2 볼륨 영상을 이동시키는 변환 함수를 결정하고, 적어도 결정된 변환 함수를 포함하는 복수의 변환 함수들에 기초하여 복수의 볼륨 영상들의 통합 변환 특성을 결정하고, 결정된 통합 변환 특성에 기초하여 복수의 볼륨 영상들의 영상 데이터로부터 볼륨 파노라마 영상의 영상 데이터를 생성한다. According to the method and apparatus for generating a volume panorama image, a plurality of volume images of a volume image that three-dimensionally represents an observation region inside a target object are input, and a first volume image of the input volume images Determining a transform function to move the second volume image so that the second volume image having the predetermined transformation function is matched with the second volume image having the second volume image having the determined transformation function, And generates image data of the volume panoramic image from the image data of the plurality of volume images based on the integrated transformation property.
Description
의료 영상을 생성하는 방법 및 장치에 관한 것으로 볼륨 영상들간의 변환 함수에 기초하여 볼륨 파노라마 영상을 생성하는 방법 및 장치에 관한 것이다. And more particularly, to a method and apparatus for generating a volume panorama image based on a transformation function between volume images.
환자를 진단하기 위한 다양한 의료 장비들이 사용 또는 개발 중에 있다. 환자 진단 과정에서의 환자의 편의, 환자 진단 결과의 신속성 등으로 인하여 초음파 영상 장치, CT(Computed Tomography), MRI(Magnetic Resonance Imaging) 등과 같이 인체 내부 단면의 모습을 영상으로 보여주는 의료 장비들의 중요성이 부각되고 있다. 초음파 영상 장치는 대상체의 내부의 소정 지점을 향하여 초음파 신호를 전달하고, 대상체 내부에서 반사된 초음파 신호의 정보를 이용하여 대상체 내부에 관한 이미지를 얻는 장치이다. 이와 같은 초음파 영상 장치는 소형이고, 저렴하며, 실시간으로 표시 가능하고, X선 등의 피폭이 없어 안전성이 높은 장점을 가지고 있다. Various medical devices for diagnosing patients are in use or under development. Because of the convenience of the patient in the diagnosis process and the promptness of the patient's diagnosis results, the importance of medical devices such as ultrasound imaging, CT (Computed Tomography) and MRI (Magnetic Resonance Imaging) . An ultrasound imaging apparatus is an apparatus that acquires an image related to an inside of a target object by transmitting an ultrasound signal toward a predetermined point inside the target object and using information of the ultrasound signal reflected inside the target object. Such an ultrasound imaging apparatus is small, inexpensive, can be displayed in real time, and has high safety because there is no exposure such as X-rays.
한편, 의료 장비들의 급속한 발전에 따라 인체 내부 단면을 나타내는 2차원 영상을 넘어서 3차원 영상을 출력하는 의료 장비들이 등장하고 있다. 나아가 보다 넓은 관찰 영역을 확보하기 위해 복수의 3차원 영상들을 합성하여 3차원 파노라마 영상을 생성하는 의료 장비들이 소개되고 있다. On the other hand, as the medical equipment rapidly develops, there are emerging medical devices that output three-dimensional images beyond the two-dimensional image showing the internal cross-section of the human body. Furthermore, medical devices that synthesize a plurality of three-dimensional images to generate a three-dimensional panoramic image are introduced to secure a wider viewing area.
복수의 3차원 볼륨 영상들간의 복수의 변환 함수들을 통합적으로 고려하여 보다 정확하게 볼륨 영상들을 정합하는 볼륨 파노라마 영상을 생성하는 방법 및 장치를 제공하는데 있다. 복수의 변환 함수들로부터 결정된 통합 변환 정보를 갱신하여 많은 개수의 볼륨 영상들을 효율적으로 처리하는 볼륨 파노라마 영상을 생성하는 방법 및 장치를 제공하는데 있다. 또한, 상기 방법들 각각을 컴퓨터에서 실행시키기 위한 프로그램을 기록한 컴퓨터로 읽을 수 있는 기록 매체를 제공하는데 있다. 본 실시예가 이루고자 하는 기술적 과제는 상기된 바와 같은 기술적 과제들로 한정되지 않으며, 또 다른 기술적 과제들이 존재할 수 있다. And a method and apparatus for generating a volume panorama image that more accurately matches volume images by integrally considering a plurality of conversion functions between a plurality of three-dimensional volume images. And a method and apparatus for generating a volume panorama image that efficiently processes a large number of volume images by updating integrated transformation information determined from a plurality of transformation functions. The present invention also provides a computer-readable recording medium having recorded thereon a program for causing a computer to execute each of the above methods. The technical problem to be solved by this embodiment is not limited to the above-described technical problems, and other technical problems may exist.
본 발명의 일측에 따르면, 볼륨 파노라마 영상 생성 방법은 상기 볼륨 영상들 중 제1 볼륨 영상과 상기 제1 볼륨 영상과 공통된 영역을 갖는 제2 볼륨 영상간의 변환 관계를 나타내는 변환 함수를 입력받는 단계, 적어도 상기 변환 함수를 포함하는 복수의 변환 함수들을 통합적으로 고려하여 상기 복수의 변환 함수들 각각으로부터 최적화 변환 함수를 생성하는 단계 및 상기 최적화 변환 함수에 기초하여 상기 볼륨 파노라마 영상을 생성하는 단계를 포함한다. According to an aspect of the present invention, there is provided a method of generating a volume panorama image, the method comprising: receiving a transformation function indicating a transformation relationship between a first volume image of the volume images and a second volume image having an area common to the first volume image; Generating an optimization conversion function from each of the plurality of conversion functions by integrally considering a plurality of conversion functions including the conversion function, and generating the volume panorama image based on the optimization conversion function.
본 발명의 다른 측면에 따라 상기 볼륨 파노라마 영상 생성 방법을 컴퓨터에서 실행시키기 위한 프로그램을 기록한 컴퓨터로 읽을 수 있는 기록매체가 제공된다. According to another aspect of the present invention, there is provided a computer-readable recording medium having recorded thereon a program for causing a computer to execute the volume panorama image generating method.
본 발명의 또 다른 측면에 따르면, 볼륨 파노라마 영상 생성 장치는 상기 복수의 볼륨 영상들의 영상 데이터들을 입력받는 입력부, 복수의 볼륨 영상들간의 변환 관계를 수식적으로 나타내는 복수의 변환 함수들을 통합적으로 고려하여 상기 복수의 변환 함수들 중 어느 하나로부터 최적화 변환 함수를 생성하고, 상기 최적화 변환 함수에 기초하여 상기 복수의 볼륨 영상들로부터 볼륨 파노라마 영상을 생성하는 영상 프로세서 및 상기 생성된 볼륨 파노라마 영상을 출력하는 출력부를 포함한다. According to another aspect of the present invention, an apparatus for generating a volume panorama image includes an input unit for receiving image data of the plurality of volume images, a plurality of conversion functions expressively representing a conversion relation between a plurality of volume images, An image processor for generating an optimization conversion function from any one of the plurality of conversion functions and generating a volume panorama image from the plurality of volume images based on the optimization conversion function and an output for outputting the generated volume panorama image .
본 발명의 또 다른 측면에 따르면, 입력부는 복수의 변환 함수들을 입력받고, 영상 프로세서는 상기 볼륨 영상들 중 제1 볼륨 영상과 상기 제1 볼륨 영상과 공통된 영역을 갖는 제2 볼륨 영상간의 변환 관계를 나타내는 변환 함수를 포함하는 복수의 변환 함수들을 통합적으로 고려하여 상기 복수의 변환 함수들 각각으로부터 최적화 변환 함수를 생성하는 최적화 변환 함수 생성부 및 상기 최적화 변환 함수에 기초하여 복수의 볼륨 영상들로부터 볼륨 파노라마 영상을 생성하는 볼륨 파노라마 영상 생성부를 포함한다. According to another aspect of the present invention, an input unit receives a plurality of conversion functions, and an image processor converts a first volume image of the volume images and a second volume image having an area common to the first volume image An optimization conversion function generation unit for generating an optimization conversion function from each of the plurality of conversion functions by integrally considering a plurality of conversion functions including a conversion function indicating a conversion function representing a volume panorama from the plurality of volume images based on the optimization conversion function, And a volume panorama image generation unit for generating an image.
본 발명의 또 다른 측면에 따르면, 영상 프로세서는 상기 볼륨 영상들 중 제1 볼륨 영상과 상기 제1 볼륨 영상과 공통된 영역을 갖는 제2 볼륨 영상간의 변환 관계를 나타내는 변환 함수를 결정하는 변환 함수 결정부, 상기 변환 함수를 포함하는 복수의 변환 함수들을 통합적으로 고려하여 상기 복수의 변환 함수들 각각으로부터 최적화 변환 함수를 생성하는 최적화 변환 함수 생성부 및 상기 최적화 변환 함수에 기초하여 복수의 볼륨 영상들로부터 볼륨 파노라마 영상을 생성하는 볼륨 파노라마 영상 생성부를 포함한다. According to another aspect of the present invention, an image processor includes a transform function determining unit determining a transform function indicating a transform relationship between a first volume image of the volume images and a second volume image having an area common to the first volume image, An optimization conversion function generation unit for collectively considering a plurality of conversion functions including the conversion function to generate an optimization conversion function from each of the plurality of conversion functions, And a volume panorama image generation unit for generating a panorama image.
복수의 3차원 볼륨 영상들간의 복수의 변환 함수들을 통합적으로 고려하여 통합 변환 특성을 생성함으로써, 보다 정확하게 볼륨 영상들을 정합할 수 있는 볼륨 파노라마 영상을 생성하는 방법 및 장치를 제공할 수 있다. 또한, 복수의 변환 함수들 각각을 통합 변환 정보의 파라미터로 결정하고, 결정된 통합 변환 특성을 최적화함으로써, 볼륨 파노라마 영상의 품질을 향상시키는 볼륨 파노라마 영상을 생성하는 방법 및 장치를 제공할 수 있다. 또한, 복수의 변환 함수들로부터 결정된 통합 변환 특성을 갱신함으로써, 많은 개수의 볼륨 영상들간의 변환 함수들을 효율적으로 관리하는 볼륨 파노라마 영상을 생성하는 방법 및 장치를 제공할 수 있다. It is possible to provide a method and apparatus for generating a volume panorama image that can more accurately match volume images by integrally considering a plurality of conversion functions between a plurality of three-dimensional volume images. It is also possible to provide a method and apparatus for generating a volume panorama image that improves the quality of a volume panorama image by determining each of the plurality of conversion functions as parameters of the integrated transformation information and optimizing the determined integrated transformation characteristics. Also, it is possible to provide a method and apparatus for generating a volume panorama image that efficiently manages conversion functions between a large number of volume images by updating the integrated transformation property determined from the plurality of conversion functions.
도 1은 본 발명의 일 실시예에 따른 의료 영상 시스템의 구성도이다.
도 2는 도 1에 도시된 볼륨 파노라마 영상 생성 장치(20)의 일 예에 따른 구성도이다.
도 3은 도 2에 도시된 영상 프로세서(22)에서 복수의 볼륨 영상들로부터 볼륨 파노라마 영상을 생성하는 모습을 도시한 도면이다.
도 4는 도 2에 도시된 최적화 변환 함수 생성부(221)에서 입력된 통합 변환 정보와 볼륨 영상들의 영상 데이터에 기초하여 최적화 변환 함수를 생성하는 과정을 설명하기 위한 도면이다.
도 5는 도 2에 도시된 최적화 변환 함수 생성부(221)의 일 예에 따른 구성도이다.
도 6은 부분 영역의 형태적인 특성을 결정하는 과정을 설명하기 위한 도면이다.
도 7은 도 6의 평활화된 구 형상의 영역(64)의 일 영역에 포함된 복셀들의 강도들의 변화하는 정도를 나타내는 지표를 나타낸 도면이다.
도 8은 도 2의 합성 영상 데이터 생성부(222)에서 합성되는 볼륨 영상의 영상 데이터를 생성하는 과정을 나타낸 동작 흐름도이다.
도 9는 합성되는 볼륨 영상으로부터 적어도 하나의 지역 볼륨 영상을 분할하는 일 실시예를 나타낸 도면이다.
도 10은 발명의 다른 실시예에 따른 볼륨 파노라마 영상 생성 장치(100)의 구성도이다.
도 11은 본 발명의 일 실시예에 따른 볼륨 파노라마 영상 생성 방법의 동작 흐름도이다.
도 12는 본 발명의 다른 실시예에 따른 볼륨 파노라마 영상 생성 방법의 동작 흐름도이다. 1 is a block diagram of a medical image system according to an embodiment of the present invention.
FIG. 2 is a block diagram illustrating an example of the volume panorama
FIG. 3 is a view showing a volume panorama image generated from a plurality of volume images in the
4 is a diagram for explaining a process of generating an optimization conversion function based on the integrated conversion information and the image data of the volume images input by the optimization conversion
5 is a block diagram illustrating an example of the optimization
6 is a diagram for explaining a process of determining the morphological characteristics of a partial area.
FIG. 7 is a graph showing the degree of variation of the intensities of voxels contained in one area of the smoothed
FIG. 8 is a flowchart illustrating a process of generating image data of a volume image synthesized by the composite image
9 is a view showing an embodiment of dividing at least one local volume image from a synthesized volume image.
FIG. 10 is a block diagram of a volume panoramic
11 is a flowchart illustrating an operation of a volume panorama image generating method according to an embodiment of the present invention.
12 is a flowchart illustrating an operation of a volume panorama image generating method according to another embodiment of the present invention.
이하에서는 도면을 참조하여 본 발명의 실시예들을 상세히 설명한다. Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings.
도 1은 본 발명의 일 실시예에 따른 의료 영상 시스템의 구성도이다. 도 1을 참조하면, 도 1에 도시된 실시예에 따른 의료 진단 시스템은 3차원 볼륨 영상 생성 장치(10), 볼륨 파노라마 영상 생성 장치(20), 영상 표시 장치(30)로 구성된다. 3차원 볼륨 영상 생성 장치(10)는 대상체(40) 내부의 관찰 영역을 3차원적으로 나타내는 볼륨 영상들의 영상 데이터들을 생성한다. 이 때, 이와 같은 3차원 볼륨 영상 생성 장치(10)의 일 예에는 초음파 진단 장치, CT(Computed Tomography), MRI(Magnetic Resonance Imaging) 등과 같이 대상체 내부 모습을 영상으로 보여주는 의료 장비들이 포함된다. 특히, 3차원 볼륨 영상 생성 장치(10)가 초음파 진단 장치인 경우, 3차원 볼륨 영상 생성 장치(10)는 이것에 장착된 프로브(11)로부터 발생된 소스 신호(source signal)가 의사 등과 같은 의료 전문가가 진단하고자 하는 대상체(40) 내부의 관찰 영역에 전달됨으로써 발생되는 반응 신호를 이용하여 이러한 관찰 영역을 3차원적으로 나타내는 볼륨 영상들의 영상 데이터들을 생성한다. 이 때, 소스 신호는 초음파, X선 등 여러 종류의 신호가 될 수 있다. 이하에서는 설명의 편의를 위하여 3차원 볼륨 영상 생성 장치(10)가 초음파(ultrasound)를 이용하여 대상체(40, 예를 들어 환자의 신체)로부터 3차원 볼륨 영상을 검출하는 초음파 진단 장치인 경우를 예로 들어 설명하도록 한다. 그러나, 본 발명의 복수의 볼륨 영상들이 초음파 진단 장치에 의해 생성되는 것으로만 한정 해석되지는 않는다. 1 is a block diagram of a medical image system according to an embodiment of the present invention. Referring to FIG. 1, the medical diagnosis system according to the embodiment shown in FIG. 1 includes a 3D volume
초음파 진단 장치에서의 프로브(11)는 일반적으로 적어도 하나의 변환기의 배열로 구성된다. 3차원 볼륨 영상 생성 장치(10)의 프로브(11)로부터 초음파 신호가 대상체(40)의 내부의 관찰 영역으로 전달되면, 이러한 초음파 신호는 여러 다른 조직들(tissue) 사이의 계층들로부터 부분적으로 반사된다. 특히, 초음파 신호는 대상체(40) 내부에서의 밀도 변화가 있는 곳, 예를 들어, 혈장(blood plasma) 내의 혈구들(blood cells), 장기들(organs) 내의 작은 조직들(structures) 등에서 반사된다. 이와 같이 반사된 초음파 신호들은 프로브(11)의 변환기들을 진동시키고, 변환기들은 이 진동들에 따른 전기적 펄스들(electrical pulses)을 출력한다. 이와 같은 전기적 펄스들이 3차원 볼륨 영상으로 변환될 수 있다. The
3차원 볼륨 영상 생성 장치(10)는 대상체(40) 위에서 프로브(11)의 위치(location)와 방향(orientation)을 변화시키면서, 대상체(40)에 대한 3차원 볼륨 영상을 검출한다. 예를 들어, 3차원 볼륨 영상 생성 장치(10)는 대상체(40)의 관찰 영역으로 복수 회의 초음파 신호들을 송신함에 따라 대상체(40)의 특정 부위에 대한 다수의 단면 영상들을 검출하고, 이와 같은 단면 영상들을 축척하여 대상체(40) 내부의 관찰 영역을 3차원적으로 나타내는 3차원 볼륨(volume) 영상의 영상 데이터를 생성한다. 이와 같이 단면 영상들을 축척하여 3차원 볼륨 영상의 영상 데이터를 생성하는 방식을 MPR(Multiplanar Recunstruction) 방식이라고 한다. 그러나, 이하에서 설명될 실시예들의 특징은 3차원 볼륨 영상을 생성에 있는 것이 아니라, 3차원 볼륨 영상들로부터 대상체(40) 내부의 보다 넓은 관찰 영역에 대한 볼륨 파노라마 영상을 얻는 것에 있다. 따라서, 앞서 설명된 3차원 볼륨 영상의 생성 과정은 하나의 예일 뿐, 다른 방식들로 생성된 3차원 볼륨 영상에 이하에서 설명될 실시예들이 적용될 수 있다. 예를 들어, x 축, y 축, z 축의 위치 데이터를 포함하는 3차원 수신 신호를 프로브(11)의 변환기들로 입력받아, 이러한 3차원 수신 신호로부터 3차원 볼륨 영상들의 영상 데이터를 생성하는 방식으로 생성된 3차원 볼륨 영상에 이하에서 설명될 실시예들이 적용될 수도 있는 것이다. The 3D volume
그런데, 이와 같이 생성된 3차원 볼륨 영상에는 한 번에 관찰 가능한 관찰 영역 (Field of View)의 크기에 제한이 수반될 수 있다. 특히, 초음파 신호로부터 생성된 초음파 3차원 볼륨 영상은 프로브의 종류, 변환기의 배열 형태, 변환기의 개수 등에 의해서 한 번에 관찰 가능한 관찰 영역의 크기가 제한될 수 있다. 이 때, 한 번에 관찰 가능한 관찰 영역은 프로브(11)를 대상체(40)의 소정 위치로 위치시키고, 위치의 이동 없이 상기 위치에서 획득하는 초음파 영상을 의미한다. 예를 들어, 본 발명의 실시예에 따른 3차원 볼륨 영상 생성 장치(10)가 대상체(40)의 피부로부터 깊이 15cm, 시야각 60 내지 100도 정도의 관찰 영역을 한 번에 관찰 가능한 경우, 3차원 볼륨 영상 생성 장치(10)로부터 출력되는 3차원 볼륨 영상은 대상체(40) 내부의 장기 기관이나 또는 태아 전체를 한 번에 관찰하기 위한 용도로 활용함에 있어 제한을 가질 수 있다. 이에 따라, 보다 넓은 시역의 확보를 위해서는 순차적으로 획득한 복수의 3차원 볼륨 영상들을 합성하여, 볼륨 파노라마 영상을 생성하는 것이 요구된다. However, the 3D volume image thus generated may be accompanied by a limitation on the size of a field of view that can be observed at one time. Particularly, the size of the observation region that can be observed at a time can be limited by the type of the probe, the arrangement type of the transducer, the number of the transducers, and the like in the ultrasonic three-dimensional volume image generated from the ultrasonic signal. At this time, the observable observation region at one time means an ultrasound image obtained by positioning the
한편, 복수의 3차원 볼륨 영상들을 합성하여 볼륨 파노라마 영상을 생성함에 있어서, 복수의 3차원 볼륨 영상들간의 정합이 요구된다. 일반적으로, 이러한 정합은 볼륨 영상들간의 변환 함수에 의하여 수행된다. 예를 들어, 복수의 볼륨 영상들 중 제1 볼륨 영상과 제2 볼륨 영상간의 정합은 제1 볼륨 영상과 제2 볼륨 영상간의 변환 관계를 나타내는 변환 함수에 의하여 수행될 수 있다. 이러한 변환 함수는 제1 볼륨 영상에 제2 볼륨 영상을 정합하기 위하여, 제2 볼륨 영상에 포함된 복셀들의 위치와 방향을 이동시켜 제1 볼륨 영상에 매칭시키는 것을 의미할 수 있다. 그런데 복수의 볼륨 영상들 중 적어도 둘간의 변환 함수는 다른 볼륨 영상들간의 변환 함수에 오차를 야기할 수 있다. 예를 들어, 순차적으로 획득되는 제1 볼륨 영상, 제2 볼륨 영상 및 제3 볼륨 영상을 합성하여 볼륨 파노라마 영상을 생성함에 있어서, 제1 볼륨 영상과 제2 볼륨 영상간의 변환 함수에 발생하는 오차는 제2 볼륨 영상과 제3 볼륨 영상간의 변환 함수에 오차를 야기할 수 있다. 이에 따라, 적어도 셋 이상의 볼륨 영상들을 합성하여 볼륨 파노라마 영상을 생성함에 있어서, 복수의 볼륨 영상들간의 변환 함수들을 통합적으로 고려하여 변환 함수들 각각으로부터 최적화 변환 함수들을 생성하는 것이 요구된다. 이러한 최적화 변환 함수들은 복수의 볼륨 영상들을 합성하여 볼륨 파노라마 영상을 생성함에 있어 오차를 최소화할 수 있다. 이하에서 설명될 실시예들은 3차원 볼륨 영상들간의 복수의 변환 함수들을 통합적으로 고려하여 복수의 변환 함수들 각각으로부터 최적화 변환 함수를 생성하고, 생성된 최적화 변환 함수에 기초하여 복수의 볼륨 영상들로 이루어진 볼륨 파노라마 영상을 생성하는 방식을 제시한다. On the other hand, in generating a volume panorama image by combining a plurality of three-dimensional volume images, matching between a plurality of three-dimensional volume images is required. In general, such matching is performed by a transform function between volume images. For example, the matching between the first volume image and the second volume image among the plurality of volume images may be performed by a conversion function indicating a conversion relationship between the first volume image and the second volume image. This conversion function may mean that the position and direction of the voxels included in the second volume image are moved to match the first volume image to match the second volume image to the first volume image. However, the conversion function between at least two of the plurality of volume images may cause an error in the conversion function between different volume images. For example, in generating a volume panorama image by sequentially synthesizing the first volume image, the second volume image, and the third volume image, the error generated in the conversion function between the first volume image and the second volume image is It is possible to cause an error in the conversion function between the second volume image and the third volume image. Accordingly, in generating a volume panorama image by synthesizing at least three volume images, it is required to integrally consider the conversion functions between a plurality of volume images to generate optimization conversion functions from each of the conversion functions. Such optimization conversion functions can minimize errors in generating a volume panorama image by synthesizing a plurality of volume images. The embodiments described below collectively consider a plurality of conversion functions between three-dimensional volume images to generate an optimization conversion function from each of the plurality of conversion functions, and generate a plurality of volume images based on the generated optimization conversion function A volume panorama image is generated.
도 2는 도 1에 도시된 볼륨 파노라마 영상 생성 장치(20)의 일 예에 따른 구성도이다. 도 2를 참조하면, 도 1에 도시된 볼륨 파노라마 영상 생성 장치(20)는 입력부(21), 영상 프로세서(22), 저장부(23) 및 출력부(24)로 구성된다. 다만, 도 2에 도시된 볼륨 파노라마 영상 생성 장치(20)는 본 발명의 하나의 구현 예에 불과하며, 도 2에 도시된 구성 요소들을 기초로 하여 여러 가지 변형이 가능함을 본 발명의 일 실시예가 속하는 기술 분야에서 통상의 지식을 가진 자라면 이해할 수 있다. 예를 들어, 볼륨 파노라마 영상 생성 장치(20)는 사용자 인터페이스를 더 포함할 수 있다. 이와 같은 사용자 인터페이스는 의료 전문가 등과 같은 사용자로부터 어떤 명령 내지 정보를 입력받기 위한 인터페이스이다. 사용자 인터페이스는 일반적으로 키보드, 마우스 등과 같은 입력 장치가 될 수도 있으나, 영상 표시 장치(30)에 표현되는 그래픽 유저 인터페이스(GUI, Graphical User interface)가 될 수도 있다. FIG. 2 is a block diagram illustrating an example of the volume panorama
입력부(21)는 3차원 볼륨 영상 생성 장치(10)로부터 볼륨 영상의 영상 데이터를 입력받는다. 이 때, 볼륨 영상은 대상체(40) 내부의 관찰 영역을 3차원적으로 나타낸다. 일반적으로, 입력부(21)는 3차원 볼륨 영상 생성 장치(10)로부터 복수의 3차원 볼륨 영상의 영상 데이터들을 입력받는다. 이 때, 복수의 볼륨 영상들 각각은 서로 다른 관찰 영역을 갖는다. 예를 들어, 복수의 볼륨 영상들 중 어느 하나는 대상체(40) 내부의 태아의 머리 부분의 관찰 영역을 갖고, 복수의 볼륨 영상들 중 다른 하나는 대상체(40) 내부의 태아의 몸통 부분의 관찰 영역을 가질 수 있다. 이와 같은 입력부(21)는 3차원 볼륨 영상 생성 장치(10)로부터 입력받은 복수의 볼륨 영상들의 영상 데이터들을 영상 프로세서(22)에 전달한다. The
본 발명의 일 실시예에 따르면, 입력부(21)는 3차원 볼륨 영상 생성 장치(10)로부터 볼륨 영상들간의 변환 함수를 입력받는다. 이 때, 볼륨 영상들간의 변환 함수는 볼륨 영상들 중 적어도 둘을 정합하기 위한 함수를 의미한다. 일반적으로, 볼륨 영상들을 정합한다는 것은 복수의 볼륨 영상들 중 어느 하나에 포함된 복셀들의 위치와 방향을 이동시켜 볼륨 영상들 중 다른 하나에 매칭시키는 것을 의미한다. 이와 같이, 볼륨 영상들간의 변환 함수는 볼륨 영상들 각각에 대응하는 복셀들간의 변환 함수를 의미할 수 있다. 예를 들어, 제1 볼륨 영상과 제2 볼륨 영상의 변환 함수는 제1 볼륨 영상에 대응하는 복셀들과 제2 볼륨 영상에 대응하는 복셀들간의 변환 함수를 의미할 수 있다. 이 때, 제1 볼륨 영상에 대응하는 복셀들과 제2 볼륨 영상에 대응하는 복셀들간의 변환 함수는 제1 볼륨 영상에 대응하는 복셀들을 기준으로 제2 볼륨 영상에 대응하는 복셀들을 제1 볼륨 영상에 대응하는 복셀들로 정합시키기 위한 제2 볼륨 영상에 대응하는 복셀들의 변환 함수를 의미한다. 제1 볼륨 영상에 대응하는 복셀들은 제1 볼륨 영상에 포함된 복셀들을 의미할 수 있다. 이와 같은 맥락으로, 제2 볼륨 영상에 대응하는 복셀들은 제2 볼륨 영상에 포함된 복셀들을 의미할 수 있다. 다만, 이에 한정되지 않는다. 예를 들어, 제1 볼륨 영상에 대응하는 복셀들은 제1 볼륨 영상에 포함된 복셀들 중 임계값 이상의 강도를 갖는 복셀들만을 의미할 수도 있다. According to an embodiment of the present invention, the
본 발명의 일 실시예에 따르면, 입력부(21)는 3차원 볼륨 영상 생성 장치(10)로부터 볼륨 영상들간의 복수의 변환 함수들을 입력받는다. 예를 들어, 입력부(21)는 3차원 볼륨 영상 생성 장치(10)로부터 볼륨 영상들 중 제1 볼륨 영상과 제2 볼륨 영상간의 변환 함수를 입력받고, 이와 비교하여 볼륨 영상들 중 제2 볼륨 영상과 제3 볼륨 영상간의 변환 함수를 입력받는다. 이와 같은, 제1 볼륨 영상, 제2 볼륨 영상 및 제3 볼륨 영상은 순차적으로 획득된 볼륨 영상들일 수도 있고, 시간의 선후에 상관 없이 무작위로 획득된 볼륨 영상들일 수도 있다. 일반적으로, 제1 볼륨 영상과 제2 볼륨 영상간에는 공통된 영역이 존재한다. 이와 마찬가지로, 제2 볼륨 영상과 제3 볼륨 영상간에는 공통된 영역이 존재한다. 다만, 제1 볼륨 영상과 제2 볼륨 영상간의 공통된 영역과 제2 볼륨 영상과 제3 볼륨 영상간의 공통된 영역은 서로 다른 영역일 수 있다. 또한, 이러한 공통된 영역은 관찰 영역을 달리하는 복수의 볼륨 영상들 각각의 관찰 영역간에 공통된 영역을 의미할 수 있다. According to an embodiment of the present invention, the
지금까지 입력부(21)가 3차원 볼륨 영상 생성 장치(10)로부터 복수의 변환 함수들간의 변환 함수를 수신하는 실시예에서 설명하였으나, 본 발명의 다른 실시예에 따르면, 입력부(21)는 3차원 볼륨 영상 생성 장치(10)로부터 복수의 볼륨 영상들의 영상 데이터들만을 입력받고, 영상 프로세서(22)에서 이와 같은 복수의 볼륨 영상들의 영상 데이터들에 기초하여 복수의 볼륨 영상들간의 변환 함수를 결정할 수도 있다. 이와 같은 다른 실시예에 대한 구체적인 설명 및 복수의 볼륨 영상들간의 변환 함수들에 관한 구체적인 설명은 이하에서 이루어진다. Although the
출력부(24)는 영상 프로세서(22)에 의해 복수의 볼륨 영상들의 영상 데이터들로부터 합성된 볼륨 파노라마 영상의 영상 데이터를 영상 표시 장치(30)로 출력한다. 이와 같은, 출력부(24)와 앞서 설명된 입력부(21) 각각은 영상 프로세서(22)와 영상 표시 장치(30)를 연결하고, 3차원 볼륨 영상 생성 장치(10)와 영상 프로세서(22)를 연결하기 위한 일종의 인터페이스(interface)이다. 영상 표시 장치(30)는 출력부(24)로부터 수신한 볼륨 파노라마 영상의 영상 데이터를 이용하여 볼륨 파노라마 영상을 표시한다. 이와 같은 영상 표시 장치(30)의 일 예에는 볼륨 파노라마 영상을 스크린 또는 종이 위에 디스플레이하는 장치가 포함된다. 다만, 이에 한정되지 아니한다. The
저장부(23)에는 영상 프로세서(22)에서 수행되는 영상 프로세싱 과정에서 발생되는 다양한 데이터들이 저장된다. 예를 들어, 저장부(23)에는 입력부(21)로부터 입력된 복수의 볼륨 영상들의 영상 데이터들이 저장되고, 복수의 볼륨 영상들간의 변환 함수들이 저장되고, 출력부(24)로 전달되는 볼륨 파노라마 영상의 영상 데이터가 저장될 수 있다. 또한, 본 발명의 다양한 실시예들에 따르면, 저장부(23)에는 이하에서 설명되는 파라미터, 통합 변환 정보 등의 데이터들이 저장될 수도 있다. 이러한 저장부(23)의 일 예에는 하드디스크드라이브, ROM(Read Only Memory), RAM(Random Access Memory), 플래쉬메모리 및 메모리카드 등이 포함된다. The
영상 프로세서(22)는 입력부(21)에 입력된 복수의 볼륨 영상들의 영상 데이터들을 이용하여 복수의 볼륨 영상들로 이루어진 볼륨 파노라마 영상의 영상 데이터를 생성한다. 이 때, 볼륨 파노라마 영상은 복수의 볼륨 영상들 각각의 관찰 영역보다 넓은 관찰 영역을 갖는 볼륨 영상을 의미한다. 도 3은 도 2에 도시된 영상 프로세서(22)에서 복수의 볼륨 영상들로부터 볼륨 파노라마 영상을 생성하는 모습을 도시한 도면이다. 도 3을 통하여 예시하면, 영상 프로세서(22)는 입력부(21)에 입력된 제1 볼륨 영상(31)의 영상 데이터와 제2 볼륨 영상(32)의 영상 데이터를 이용하여 제1 볼륨 영상(31)의 관찰 영역 또는 제2 볼륨 영상(32)의 관찰 영역보다 넓은 관찰 영역을 갖는 볼륨 파노라마 영상(34)의 영상 데이터를 생성한다. The
영상 프로세서(22)는 볼륨 영상들간의 변환 관계를 나타내는 변환 함수에 기초하여 볼륨 영상들로부터 볼륨 파노라마 영상을 생성한다. 이 때, 변환 함수는 볼륨 파노라마 영상을 생성하기 위하여 합성되는 복수의 볼륨 영상들 중 어느 하나를 기준으로 복수의 볼륨 영상들 중 다른 하나를 이동시키는 것을 의미한다. 이와 같은 변환 함수는 볼륨 영상들 상호간을 정합하기 위한 함수를 의미할 수도 있다. 일반적으로, 영상 프로세서(22)는 볼륨 영상들간의 변환 관계를 나타내는 변환 함수를 결정하고, 결정된 변환 함수에 기초하여 볼륨 영상들을 정합하고, 정합된 볼륨 영상들을 합성하여 볼륨 파노라마 영상을 생성한다. 도 3을 통해 예시하면, 영상 프로세서(22)는 복수의 볼륨 영상들 중 제1 볼륨 영상(31)과 제2 볼륨 영상(32)간의 변환 함수를 결정하고, 결정된 변환 함수에 기초하여 제2 볼륨 영상(32)으로부터 제1 볼륨 영상과 합성되는 볼륨 영상(33)을 생성하고, 생성된 볼륨 영상(33)과 제1 볼륨 영상(31)을 합성하여 볼륨 파노라마 영상(34)을 생성할 수 있다. The
영상 프로세서(22)는 최적화 변환 함수에 기초하여 복수의 볼륨 영상들로부터 볼륨 파노라마 영상을 생성한다. 이 때, 최적화 변환 함수는 볼륨 영상들간의 형태적인 특성의 유사도를 최대화시키는 함수를 의미한다. 도 3을 통해 예시하면, 최적화 변환 함수는 제1 볼륨 영상(31)과 제2 볼륨 영상(32)간의 형태적인 특성의 유사도를 최대화하는 제1 볼륨 영상(31)과 제2 볼륨 영상(32)간의 변환 관계를 나타내는 변환 함수를 의미한다. 다만, 본 발명의 다양한 실시예들에 따르면, 이와 같은 최적화 변환 함수는 다양하게 정의될 수 있다. 도 3을 통해 예시하면, 최적화 변환 함수는 제1 볼륨 영상(31)과 합성되는 볼륨 영상(33)간의 형태적인 특성의 유사도를 최대화하는 제1 볼륨 영상(31)과 제2 볼륨 영상(32)간의 변환 관계를 나타내는 변환 함수를 의미할 수도 있다. 또한, 본 발명의 일 실시예에 따른 영상 프로세서(22)는 입력부(21)에 입력된 복수의 볼륨 영상들간의 복수의 변환 함수들을 통합적으로 고려하여, 복수의 변환 함수들 각각으로부터 최적화 변환 함수를 생성하고, 생성된 최적화 변환 함수에 기초하여 볼륨 파노라마 영상의 영상 데이터를 생성한다. 이와 같은 영상 프로세서(22)의 동작에 대해서는 이하 보다 구체적으로 설명된다. The
도 2를 참조하면, 영상 프로세서(22)는 최적화 변환 함수 생성부(221), 합성 영상 데이터 생성부(222), 볼륨 파노라마 영상 생성부(223)로 구성된다. 이러한 영상 프로세서(22)는 상기된 바와 같은 구성 요소들의 기능을 수행하는 전용 칩(chip)들로 제작될 수도 있고, 범용 CPU와 저장부(23)에 저장된 전용 프로그램으로 구현될 수도 있다. 2, the
최적화 변환 함수 생성부(221)는 볼륨 영상들간의 변환 관계를 나타내는 변환 함수를 입력부(21)로부터 입력받는다. 이 때, 변환 함수는 볼륨 영상들을 합성하는 볼륨 파노라마 영상을 생성함에 있어서, 볼륨 영상들간의 정합을 위해 결정되는 함수를 의미한다. 도 3을 통해 예시하면, 볼륨 영상들 중 제1 볼륨 영상(31)과 제2 볼륨 영상(32)간의 변환 함수는 제1 볼륨 영상(31)을 기준으로 제2 볼륨 영상(32)을 제1 볼륨 영상(31)에 정합시키기 위해서 제2 볼륨 영상(32)에 적용되는 함수를 의미할 수 있다. 또한, 볼륨 영상들간의 변환 함수는 볼륨 영상들 각각에 대응하는 복셀들간의 변환 함수를 의미할 수 있다. 도 3을 통하여 예시하면, 복수의 볼륨 영상들 중 제1 볼륨 영상(31)과 제2 볼륨 영상(32)간의 변환 함수는 제1 볼륨 영상(31)에 대응하는 복셀들과 제2 볼륨 영상(32)에 대응하는 복셀들간의 변환 함수를 의미할 수 있다. 일반적으로, 볼륨 영상들 각각에 대응하는 복셀들은 볼륨 영상들 각각에 포함된 복셀들을 의미한다. 다만, 본 발명의 다양한 실시예에 따르면, 이와 같은 볼륨 영상들 각각에 대응하는 복셀들은 다양한 형태로 결정될 수도 있다. 도 3을 통해 예시하면, 제1 볼륨 영상(31)에 대응하는 복셀들은 제1 볼륨 영상(31) 내 소정 영역에 포함된 복셀들 그리고 그 주변부의 복셀들을 의미할 수도 있고, 제1 볼륨 영상(31)에 포함된 복셀들 중 관찰 대상(예를 들어, 태아)에 포함되는 복셀들을 의미할 수도 있는 것이다. The optimization conversion
최적화 변환 함수 생성부(221)는 복수의 볼륨 영상들 중 어느 하나의 볼륨 영상과 이러한 어느 하나의 볼륨 영상과 공통된 영역을 갖는 복수의 볼륨 영상들 중 다른 하나의 볼륨 영상간의 변환 함수를 입력받는다. 또한, 최적화 변환 함수 생성부(221)는 이러한 변환 함수와 다른 제2 변환 함수를 더 입력받는다. 도 4는 도 2에 도시된 최적화 변환 함수 생성부(221)에서 입력된 통합 변환 정보와 볼륨 영상들의 영상 데이터에 기초하여 최적화 변환 함수를 생성하는 과정을 설명하기 위한 도면이다. 도 4를 통하여 예시하면, 최적화 변환 함수 생성부(221)는 복수의 볼륨 영상들 중 제1 볼륨 영상(41)과 이러한 제1 볼륨 영상(41)과 공통된 영역을 갖는 제2 볼륨 영상(42)간의 제1 변환 함수를 입력받고, 복수의 볼륨 영상들 중 제2 볼륨 영상(42)과 이러한 제2 볼륨 영상(42)과 앞서 언급된 공통된 영역과 다른 공통된 영역을 갖는 제3 볼륨 영상(43)간의 제2 변환 함수를 입력받을 수 있다. 이 때, 제1 볼륨 영상(41), 제2 볼륨 영상(42) 및 제3 볼륨 영상(43) 각각은 입력부(21)로부터 순차적으로 입력되고, 제1 변환 함수와 제2 변환 함수 역시 순차적으로 입력될 수 있다. 다시 말하면, 최적화 변환 함수 생성부(221)는 순차적으로 획득된 볼륨 영상들 각각의 순서를 기준으로, 제1 볼륨 영상(41)과 제2 볼륨 영상(42)간의 제1 변환 함수를 입력받고, 제2 볼륨 영상(42)과 제3 볼륨 영상(43)간의 제2 변환 함수를 입력받을 수 있다. 다만, 본 발명의 다른 실시예에 따르면, 최적화 변환 함수 생성부(221)는 복수의 볼륨 영상들 중 어느 하나를 제1 볼륨 영상으로 결정하고, 결정된 제1 볼륨 영상과 다른 볼륨 영상들 각각과의 변환 함수들을 볼륨 영상들 각각의 획득 또는 입력 순서와 상관없이 입력받을 수도 있다. The optimization conversion
최적화 변환 함수 생성부(221)는 복수의 변환 함수들을 통합적으로 고려하여 복수의 변환 함수들 각각으로부터 최적화 변환 함수를 생성한다. 일반적으로, 복수의 변환 함수들을 통합적으로 고려한다는 것은 복수의 변환 함수들 중 어느 하나의 변환 함수로부터 최적화 변환 함수를 생성함에 있어서, 이러한 어느 하나의 변환 함수와 복수의 변환 함수들 중 다른 변환 함수들을 함께 고려하여 어느 하나의 변환 함수로부터 최적화 변환 함수를 생성하는 것을 의미한다. 이 때, 최적화 변환 함수는 변환 함수들 각각으로부터 변경된 것을 의미한다. 일반적으로, 최적화 변환 함수는 복수의 변환 함수들간의 형태적인 특성의 유사도를 최대화하는 변환 함수를 의미할 수 있으나, 이와 같은 정의에 한정 해석되지는 않는다. The optimization conversion
도 4를 참조하면, 최적화 변환 함수 생성부(221)는 복수의 변환 함수들을 통합적으로 고려하여 제1 변환 함수로부터 제1 최적화 변환 함수를 생성하고, 제2 변환 함수로부터 제2 최적화 변환 함수를 생성한다. 이 때, 제1 변환 함수는 앞서 설명된 바와 같이 제1 볼륨 영상(41)과 제2 볼륨 영상(42)간의 변환 관계를 나타내는 변환 함수이고, 제2 변환 함수는 제2 볼륨 영상(42)과 제3 볼륨 영상(43)간의 변환 관계를 나타내는 변환 함수로 가정된다. 일반적으로, 최적화 변환 함수 생성부(221)는 복수의 볼륨 영상들간의 유사도를 결정하고, 복수의 유사도들을 통합적으로 고려하여 최적화 변환 함수를 생성한다. 도 4를 통해 예시하면, 최적화 변환 함수 생성부(221)는 복수의 볼륨 영상들 중 제1 볼륨 영상(41)과 제2 볼륨 영상(42)간의 제1 변환 함수에 기초하여 제1 볼륨 영상(41)과 제2 볼륨 영상(42)간의 유사도를 결정하고, 복수의 볼륨 영상들 중 제2 볼륨 영상(42)과 제3 볼륨 영상(43)간의 제2 변환 함수에 기초하여 제2 볼륨 영상(42)과 제3 볼륨 영상(43)간의 유사도를 결정하고, 복수의 유사도들을 통합적으로 고려하여 제1 변환 함수로부터 제1 최적화 변환 함수를 생성하고, 제2 변환 함수로부터 제2 최적화 변환 함수를 생성한다. 이 때, 앞서 설명된 바와 같이, 제1 볼륨 영상(41)과, 제2 볼륨 영상(42)간의 유사도는 제1 볼륨 영상(41)과 제1 변환 함수에 기초하여 제2 볼륨 영상(42)으로부터 생성된 볼륨 영상간의 유사도를 의미할 수도 있다. Referring to FIG. 4, the optimization conversion
일반적으로, 볼륨 영상들간의 유사도의 합은 볼륨 영상들간의 형태적인 특성의 유사도를 의미한다. 이와 같은 볼륨 영상들 각각의 형태적인 특성은 볼륨 영상들 각각에 대응하는 복셀들에 의해 결정된다. 예를 들어, 볼륨 영상들 각각의 형태적인 특성은 볼륨 영상들 각각에 대응하는 복셀들의 위치 정보 및 정보량에 따라 정의될 수 있다. 이 때, 정보량의 일 예는 복셀들 각각의 강도(intensity)이다. 또한, 볼륨 영상들간의 형태적인 특성의 유사도는 볼륨 영상들간의 상호 정보관계 (Mutual Information)로 정의될 수 있다. 이 때, 상호 정보관계는 평활화된 상호 정보관계 (Normalized Mutual Information)를 의미할 수도 있다. 다만, 이와 같은 상호 정보관계는 볼륨 영상들간의 형태적인 특성의 유사도를 결정하는 일 예에 불과하며, 본 발명의 다양한 실시예들에 따르면 이러한 유사도는 다양하게 결정될 수 있다. 예를 들어, 볼륨 영상들 각각에 대응하는 복셀들간의 강도 분포의 유사도를 의미할 수도 있고, 볼륨 영상들 각각에 포함된 부분 영역들간의 유사도를 의미할 수도 있고, 볼륨 영상들 각각에 포함된 에지를 구성하는 복셀들간의 유사도를 의미할 수도 있다. 이에 대해서는 이하 자세하게 설명된다. Generally, the sum of similarities between volume images means the similarity of morphological characteristics between volume images. The morphological characteristics of each of these volume images are determined by the voxels corresponding to the volume images. For example, the morphological characteristic of each volume image may be defined according to the position information and the information amount of the voxels corresponding to each volume image. An example of the amount of information is the intensity of each of the voxels. In addition, the similarity of morphological characteristics between volume images can be defined as mutual information between volume images. In this case, the mutual information relationship may mean a normalized mutual information relationship. However, such a mutual information relationship is merely an example of determining the similarity of morphological characteristics between volume images, and according to various embodiments of the present invention, such similarity can be determined in various ways. For example, it may mean the similarity of the intensity distribution between voxels corresponding to each of the volume images, the similarity between the partial regions included in each volume image, May be used to refer to the similarity between the voxels. This will be described in detail below.
일반적으로, 최적화 변환 함수 생성부(221)는 복수의 유사도들의 합을 최대화하도록 복수의 변환 함수들 각각을 변경하고, 변경된 변환 함수로부터 최적화 변환 함수를 생성한다. 도 4를 통해 예시하면, 최적화 변환 함수 생성부(221)는 제1 볼륨 영상(41)과 제2 볼륨 영상(42)간의 형태적인 특성의 유사도와 제2 볼륨 영상(42)과 제3 볼륨 영상(43)간의 형태적인 특성의 유사도의 합이 최대화되도록, 제1 변환 함수와 제2 변환 함수를 변경하고, 변경된 제1 변환 함수 및 제2 변환 함수 각각을 제1 최적화 변환 함수 및 제2 최적화 변환 함수로 결정할 수 있다. 다만, 본 발명의 다양한 실시예들에 따르면, 최적화 변환 함수 생성부(221)는 복수의 유사도들의 합이 기 결정된 임계값에 근사하도록 복수의 변환 함수들 각각을 변경할 수도 있고, 유사도들을 대신하여 다른 파라미터들의 합이 최대화 또는 최소화되도록 복수의 변환 함수들 각각을 변경할 수도 있다. Generally, the optimization conversion
최적화 변환 함수 생성부(221)로부터 출력되는 최적환 변환 함수는 복수의 볼륨 영상들로부터 볼륨 파노라마 영상을 생성하기 위해 이용된다. 도 4를 통해 예시하면, 최적화 변환 함수 생성부(221)로부터 출력되는 최적환 변환 함수는 제1 볼륨 영상(41), 제2 볼륨 영상(42) 및 제3 볼륨 영상(43)으로부터 볼륨 파노라마 영상(44)을 생성하기 위하여 이용될 수 있다. 이에 대해서는 이하에서 보다 상세하게 설명된다. The optimal conversion function output from the optimization
도 5는 도 2에 도시된 최적화 변환 함수 생성부(221)의 일 예에 따른 구성도이다. 도 5를 참조하면, 최적화 변환 함수 생성부(221)는 통합 변환 정보 생성부(2211) 및 통합 변환 정보 최적화부(2212)로 구성된다. 최적화 변환 함수 생성부(221)는 복수의 변환 함수들을 통합적으로 나타내는 통합 변환 정보를 생성한다. 이 때, 통합 변환 정보는 복수의 변환 함수들 각각으로부터 추출된 적어도 하나의 파라미터로 이루어진 벡터를 포함한다. 또한, 이러한 파라미터의 일 예에는 볼륨 영상들간의 방향의 변환 관계를 나타내는 파라미터 및 볼륨 영상들간의 위치의 변환 관계를 나타내는 파라미터 중 적어도 하나가 포함된다. 일반적으로, 볼륨 영상들간의 방향의 변환 관계를 나타내는 파라미터 및 볼륨 영상들간의 위치의 변환 관계를 나타내는 파라미터 각각은 볼륨 영상들 각각에 대응하는 복셀들간의 방향의 변환 관계를 나타내는 파라미터 및 볼륨 영상들 각각에 대응하는 복셀들간의 위치의 변환 관계를 나타내는 파라미터를 의미할 수 있다. 5 is a block diagram illustrating an example of the optimization
통합 변환 정보 생성부(2211)는 복수의 변환 함수들을 이용하여 통합 변환 정보를 생성한다. 앞서 설명된 바와 같이, 복수의 변환 함수들 각각은 볼륨 영상들간의 변환 관계를 나타내며, 본 발명의 일 실시예에 따르면, 볼륨 영상들간의 변환 관계는 볼륨 영상들 각각에 대응하는 복셀들간의 변환 관계를 나타낸다. 예를 들어, 제1 볼륨 영상과 제2 볼륨 영상간의 제1 변환 함수는 제1 볼륨 영상에 포함된 복셀들 중 어느 하나가 제2 볼륨 영상에 포함된 복셀들 중 어느 하나로 변환하는 관계를 나타낼 수 있다. 이를 수식적으로 표현하면, 수학식 1과 같다. 구체적으로, 수학식 1은 N개의 볼륨 영상들과 N-1개의 볼륨 영상들간의 변환 함수들이 통합 변환 정보 생성부(2211)로 입력되는 경우, N-1번째 볼륨 영상과 N번째 볼륨 영상간의 변환 함수를 정의한다. 다만, N은 2 이상의 정수이다. 이 때, 앞서 설명된 바와 같이, 볼륨 영상들간의 변환 함수는 볼륨 영상들 각각에 대응하는 복셀들간의 변환 함수를 의미할 수 있다. 따라서, 수학식 1을 다시 해석하면, 수학식 1은 N개의 볼륨 영상들 중 N번째 볼륨 영상에 대응하는 임의의 복셀 xn으로부터 볼륨 영상들 중 N-1번째 볼륨 영상에 대응하는 임의의 복셀 xn-1으로의 변환 관계를 나타낸 것이다. 이 때, An,n-1은 제2 볼륨 영상에 대응하는 복셀들로부터 제1 볼륨 영상에 대응하는 복셀들로의 방향의 변환 관계를 나타내는 파라미터이고, Tn,n-1은 제1 볼륨 영상에 대응하는 복셀들로의 위치의 변환 관계를 나타내는 파라미터를 의미할 수 있다. The integrated conversion
[수학식 1][Equation 1]
수학식 1에서 N-1번째 볼륨 영상과 N번째 볼륨 영상간의 변환 함수를 정의하는 (An,n-1, Tn,n-1)을 입력으로 가정하는 경우, (An,m, Tn,m)의 An,m과 Tn,m 각각은 수학식 2와 같이 정의된다. If you assume the N-1 first volume image and the N second volume image converter type the (A n, n-1, T n, n-1) that defines the function between the Equation 1, (A n, m, T n, m ), A n, m and T n, m, respectively, are defined as in Equation (2).
[수학식 2]&Quot; (2) "
각각의 변환 함수는 일반적으로 복수 개의 파라미터로 정의된다. 이 때, 파라미터의 일 예에는 볼륨 영상들간의 방향의 변환 관계를 나타내는 파라미터 및 볼륨 영상들간의 위치의 변환 관계를 나타내는 파라미터 중 적어도 하나가 포함된다. 예를 들어, 제1 볼륨 영상과 제2 볼륨 영상간의 변환 함수는 제1 볼륨 영상과 제2 볼륨 영상간의 방향의 변환 관계를 나타내는 파라미터 및 제1 볼륨 영상과 제2 볼륨 영상간의 위치의 변환 관계를 나타내는 파라미터로 정의될 수 있다. 또한, 볼륨 영상들간의 위치 또는 방향의 변환 관계를 나타내는 파라미터 각각은 볼륨 영상들 각각에 대응하는 복셀들간의 위치 또는 방향의 변환 관계를 나타내는 파라미터 각각일 수 있음은 앞서 설명된 바와 같다. Each transform function is generally defined by a plurality of parameters. At this time, one example of the parameters includes at least one of a parameter indicating a conversion relationship between directions of volume images and a parameter indicating a conversion relationship of positions between volume images. For example, the conversion function between the first volume image and the second volume image may include a parameter indicating a conversion relationship between directions of the first volume image and the second volume image, and a conversion relationship between positions of the first volume image and the second volume image And the like. Each of the parameters indicating the positional relationship or the directional relationship between the volume images may be a parameter indicating the positional relationship or the directional relationship between voxels corresponding to each volume image, as described above.
일반적으로, 변환 함수를 정의하는 (An,n-1, Tn,n-1)은 복수 개의 파라미터로 정의된다. 이 때, 변환 함수가 리지드 변환 함수(Rigid Transformation)인 경우, (An,n-1, Tn,n-1)은 6 개 또는 7 개의 파라미터로 표현될 수 있다. 예를 들어, 6 개의 파라미터들로 표현되는 경우, 6 개의 파라미터들 중 3 개의 파라미터들은 방향의 변환을 정의하는 파라미터들이고, 나머지 3 개의 파라미터들은 위치 이동을 정의하는 파라미터들일 수 있다. 이 때, 방향의 변환을 정의하는 파라미터들의 일 예는 3 개의 오일러 각(Euler Angle)들이고, 위치 이동을 정의하는 파라미터들의 일 예는 3 개의 전환 벡터(Translation Vector)들이다. 다른 예를 들어, 7 개의 파라미터들로 표현되는 경우, 7 개의 파라미터들 중 4 개의 파라미터들은 방향의 변환을 정의하는 파라미터들이고, 나머지 3 개의 파라미터들은 위치 이동을 정의하는 파라미터들일 수 있다. 이 때, 방향의 변환을 정의하는 파라미터들의 일 예는 4개의 쿼터니언(quaternion) 구성(element)들이고, 위치 이동을 정의하는 파라미터들의 일 예는 3 개의 전환 벡터(Translation Vector)들이다. 또한, 본 발명의 다른 실시예에 따라, 변환 함수가 아핀 변환 함수(Affine Transformation)인 경우, (An,n-1, Tn,n-1)은 6 개 또는 7 개의 파라미터로 표현될 수도 있다. 일반적으로, 리지드 변환 함수는 이동변환과 회전변환을 나타내며 물체(예를 들어, 볼륨 영상) 자체의 모습은 변화가 없는 변환을 의미한다. 다시 말하면, 리지드 변환 함수는 유클리드 공간(Euclidean Space) 상에서 모든 점들 사이의 거리를 보전하는 변환을 의미할 수도 있다. 또한, 아핀 변환 함수는 n 차원의 공간 내 점들로부터 변환된 점들로의 변환 관계를 1차식으로 나타낸 변환 함수를 의미할 수 있다. 다만, 이러한 정의에 한정 해석되는 것은 아니다. Generally, (A n, n-1 , T n, n-1 ) defining the transform function is defined by a plurality of parameters. In this case, (A n, n-1 , T n, n-1 ) can be expressed by six or seven parameters when the transform function is a rigid transformation function. For example, when represented by six parameters, three of the six parameters are parameters defining the directional transformation, and the remaining three parameters may be parameters defining the positional movement. One example of parameters defining the directional transformation is three Euler angles, and an example of the parameters defining the translation is three translation vectors. In another example, when represented by seven parameters, four of the seven parameters are parameters defining the directional transformation, and the remaining three parameters may be parameters defining the position movement. One example of the parameters defining the directional transformation is four quaternion elements, and one example of the parameters defining the translation is three translation vectors. According to another embodiment of the present invention, when the transform function is an affine transform function, (A n, n-1 , T n, n-1 ) may be expressed by six or seven parameters have. In general, the rigid transformation function refers to a transformation that represents a motion transformation and a rotation transformation, and the appearance of an object (e.g., a volume image) itself is unchanged. In other words, the rigid transformation function may mean a transformation that conserves the distance between all points on the Euclidean space. In addition, the affine transformation function may mean a transformation function that linearly expresses a transformation relation from n-dimensional in-space points to transformed points. However, the present invention is not limited to these definitions.
통합 변환 정보 생성부(2211)는 복수의 변환 함수들을 이용하여 통합 변환 정보를 생성한다. 이 때, 통합 변환 정보 생성부(2211)는 복수의 변환 함수들 각각을 정의하는 파라미터들을 이용하여 통합 변환 정보를 생성한다. 예를 들어, 통합 변환 정보 생성부(2211)는 제1 변환 함수를 나타내는 6 개의 파라미터들과 제2 변환 함수를 나타내는 6 개의 파라미터들을 이용하여 통합 변환 정보를 생성할 수 있다. 이 때, 제1 변환 함수는 제1 볼륨 영상과 제2 볼륨 영상간의 변환 함수를 나타내고, 제2 변환 함수는 제2 볼륨 영상과 제3 볼륨 영상간의 변환 함수를 나타냄은 앞서 설명된 예시와 같다. 일반적으로, 통합 변환 정보 생성부(2211)는 복수의 변환 함수들 각각의 파라미터들을 벡터로 정의하고, 정의된 벡터 형태의 파라미터들을 이용하여 통합 변환 정보를 생성한다. 따라서, 이와 같은 통합 변환 정보는 복수개의 벡터들로부터 정의된 벡터가 될 수 있다. 이러한 통합 변환 정보의 일 예를 수식적으로 표현하면, 수학식 3과 같다. 이 때, 통합 변환 정보 v를 구성하는 v2, 1 내지 vN,N-1각각은 복수의 변환 함수들 각각으로부터 추출된 파라미터들을 나타내는 벡터를 의미할 수 있다. 예를 들어, n=2 내지 N인 경우, vn,n-1은 앞서 설명된 (An,n-1, Tn,n-1)을 나타내는 복수 개의 파라미터들을 벡터로 표현한 것을 의미할 수 있다. 지금까지의 설명은 복수의 변환 함수들이 기 정의되어 통합 변환 정보 생성부(2211)로 입력되는 것을 가정하여 설명되었으나, 앞서 언급된 바와 같이 본 발명의 다른 실시예에 따르면, 통합 변환 정보 생성부(2211)는 복수의 볼륨 영상들의 영상 데이터들만을 입력받고 스스로 복수의 볼륨 영상들간의 변환 함수를 정의할 수도 있다. The integrated conversion
[수학식 3]&Quot; (3) "
통합 변환 정보 최적화부(2212)는 생성된 통합 변환 정보에 기초하여 최적화 변환 함수를 생성한다. 이 때, 통합 변환 정보는 앞서 살펴본 바와 같이 복수의 변환 함수들 각각을 나타내는 파라미터들을 모두 포함한다. 따라서, 통합 변환 정보 최적화부(2212)에 의해 통합 변환 정보가 고려된다는 것은 통합 변환 정보 최적화부(2212)에 의해 복수의 변환 함수들이 통합적으로 고려된다는 것을 의미할 수 있다. 이와 같이, 통합 변환 정보 최적화부(2212)는 복수의 변환 함수들을 통합적으로 고려하여 최적화 변환 함수를 생성한다. The integrated transformation
통합 변환 정보 최적화부(2212)는 생성된 통합 변환 정보에 기초하여 복수의 변환 함수들 각각으로부터 최적화 변환 함수를 생성한다. 앞서 설명된 바와 같이, 통합 변환 정보는 복수의 변환 함수들의 정보들을 모두 포함한다. 따라서, 통합 변환 정보 최적화부(2212)에 의하여 통합 변환 정보에 기초하여 복수의 변환 함수들 각각으로부터 최적화 변환 함수를 생성한다는 것은 통합 변환 정보에 포함된 복수의 변환 함수들 각각의 정보를 변경하여, 복수의 변환 함수들 각각으로부터 변경된 변환 함수들 각각을 생성하고, 생성된 변환 함수들 각각을 최적화 변환 함수로서 결정한다는 것을 의미할 수 있다. 이 때, 복수의 변환 함수들 각각의 정보를 변경하는 기준은 앞서 설명된 바와 같이 볼륨 영상들간의 유사도들의 합을 최대화시키 위한 복수의 변환 함수들 각각의 정보를 결정하는 것을 의미할 수 있다. 결론적으로, 통합 변환 정보 최적화부(2212)는 생성된 통합 변환 정보에 포함된 복수의 변환 함수들 각각의 정보를 볼륨 영상들간의 유사도들의 합이 최대화되도록 변경하고, 변경된 정보에 기초하여 각각의 변환 함수로부터 각각의 최적화 변환 함수를 생성하는 것이다. The integrated conversion
볼륨 영상들간의 유사도를 수식적으로 나타내면 수학식 4와 같다. 수학식 4에서 sn은 복수의 볼륨 영상들 중 n 볼륨 영상을 기준으로, n 볼륨 영상과 복수의 볼륨 영상들 중 다른 영상들인 m 볼륨 영상들 각각간의 유사도들의 합을 나타낸다. 예를 들어, 제1 볼륨 영상, 제2 볼륨 영상, 제3 볼륨 영상이 입력되는 경우, sn은 제1 볼륨 영상과 제2 볼륨 영상간의 유사도와 제1 볼륨 영상과 제3 볼륨 영상간의 유사도의 합을 의미할 수 있다. 다른 예를 들어, 제1 볼륨 영상, 제2 볼륨 영상, 제3 볼륨 영상이 순차적으로 입력되는 경우, sn은 중간에 입력된 제2 볼륨 영상과 최초 입력된 제1 볼륨 영상간의 유사도와 중간에 입력된 제2 볼륨 영상과 최후 입력된 제3 볼륨 영상간의 유사도의 합을 의미할 수도 있다. The similarity degree between the volume images can be expressed by Equation (4). In Equation (4), s n represents the sum of similarities between the n-volume image and the m-volume images, which are different images among the plurality of volume images, with reference to the n-volume image among the plurality of volume images. For example, when a first volume image, a second volume image, and a third volume image are input, s n is a degree of similarity between the first volume image and the second volume image, and a degree of similarity between the first volume image and the third volume image It can mean sum. In another example, when the first volume image, the second volume image, and the third volume image are sequentially input, s n is a degree of similarity between the second volume image input in the middle and the first volume image input first, May mean the sum of similarities between the input second volume image and the last input third volume image.
[수학식 4]&Quot; (4) "
또한, 수학식 4로부터 N 개의 볼륨 영상들에 대해서 볼륨 영상들간의 유사도들의 합을 표현하는 경우, 수학식 5와 같이 표현될 수 있다. 이 때, 볼륨 영상들간의 유사도들의 합인 S는 복수의 볼륨 영상들에 의해 가능한 볼륨 영상 쌍들의 유사도들을 모두 합한 값을 의미할 수 있다. Also, from Equation (4), when the sum of similarities between volume images is represented with respect to N volume images, it can be expressed as Equation (5). In this case, S, which is the sum of the similarities between the volume images, may be a sum of the similarities of the volume image pairs that are possible by the plurality of volume images.
[수학식 5]&Quot; (5) "
통합 변환 정보 최적화부(2212)는 복수의 볼륨 영상들간의 유사도들을 결정하고, 복수의 유사도들을 통합적으로 고려하여 최적화 변환 함수를 생성한다. 일반적으로, 통합 변환 정보 최적화부(2212)는 복수의 유사도들의 합이 최대화되도록 통합 변환 정보를 갱신하고, 갱신된 통합 변환 정보에 기초하여 각각의 변환 함수들로부터 최적화 변환 함수를 생성한다. 수학식 5를 통해 예시하면, 통합 변환 정보 최적화부(2212)는 볼륨 영상들간의 유사도들의 합인 S가 최대화되도록 통합 변환 정보를 갱신하고, 갱신된 통합 변환 정보에 기초하여 각각의 변환 함수들로부터 최적화 변환 함수를 생성한다. 이 때, 통합 변환 정보를 갱신한다는 것은 통합 변환 정보에 포함된 복수의 변환 함수들 각각의 파라미터들을 갱신하는 것을 의미하고, 갱신된 통합 변환 정보에 기초하여 각각의 변환 함수들로부터 최적화 변환 함수를 생성하는 것은 갱신된 파라미터들에 기초하여 복수의 변환 함수들 각각으로부터 최적화 변환 함수를 생성하는 것을 의미할 수 있다. The integrated transformation
앞서 설명된 바와 같이, 볼륨 영상들간의 유사도는 볼륨 영상들간의 형태적인 특성의 유사도를 의미한다. 이와 같은 볼륨 영상들 각각의 형태적인 특성은 볼륨 영상들 각각에 대응하는 복셀들에 의해 결정된다. 예를 들어, 볼륨 영상들 각각의 형태적인 특성은 볼륨 영상들 각각에 대응하는 복셀들의 위치 정보 및 정보량에 따라 정의될 수 있다. 이 때, 정보량의 일 예는 복셀들 각각의 강도(intensity)이다. 또한, 볼륨 영상들간의 형태적인 특성의 유사도는 볼륨 영상들간의 상호 정보관계 (Mutual Information)로 정의될 수 있다. 이 때, 상호 정보관계는 평활화된 상호 정보관계 (Normalized Mutual Information)를 의미할 수도 있다. 다만, 이와 같은 상호 정보관계는 볼륨 영상들간의 형태적인 특성의 유사도를 결정하는 일 예에 불과하며, 본 발명의 다양한 실시예들에 따르면 이러한 유사도는 다양하게 결정될 수 있다. As described above, the similarity between volume images means the similarity of morphological characteristics between volume images. The morphological characteristics of each of these volume images are determined by the voxels corresponding to the volume images. For example, the morphological characteristic of each volume image may be defined according to the position information and the information amount of the voxels corresponding to each volume image. An example of the amount of information is the intensity of each of the voxels. In addition, the similarity of morphological characteristics between volume images can be defined as mutual information between volume images. In this case, the mutual information relationship may mean a normalized mutual information relationship. However, such a mutual information relationship is merely an example of determining the similarity of morphological characteristics between volume images, and according to various embodiments of the present invention, such similarity can be determined in various ways.
통합 변환 정보 최적화부(2212)는 볼륨 영상들간의 복수의 변환 함수들에 최적화 알고리즘을 적용하여 복수의 유사도들의 합이 최대화시키는 최적화 변환 함수를 결정한다. 다시 말하면, 통합 변환 정보 최적화부(2212)는 최적화 알고리즘에 기초하여 복수의 유사도들의 합을 최대화시키도록 볼륨 영상들간의 복수의 변환 함수들 각각을 대표하는 통합 변환 정보의 파라미터들을 갱신하고, 갱신된 파라미터들을 이용하여 복수의 변환 함수들 각각에 대응하는 최적화 변환 함수를 결정할 수 있다. 예를 들어, 통합 변환 정보가 제1 볼륨 영상과 제2 볼륨 영상간의 제1 변환 함수로부터 추출된 제1 파라미터와 제2 볼륨 영상과 제3 볼륨 영상간의 제2 변환 함수로부터 추출된 제2 파라미터로 구성되는 경우, 통합 변환 정보 최적화부(2212)는 최적화 알고리즘을 적용하여 제1 볼륨 영상과 제2 볼륨 영상간의 유사도와 제2 볼륨 영상과 제3 볼륨 영상간의 유사도의 합을 최대화하는 제1 파라미터와 제2 파라미터를 계산하고, 계산된 제1 파라미터에 기초하여 제1 변환 함수에 대응하는 제1 최적화 변환 함수를 생성하고, 계산된 제2 파라미터에 기초하여 제2 변환 함수에 대응하는 제2 최적화 변환 함수를 생성한다. 이 때, 최적화 알고리즘의 일 예는 Downhill Simplex 알고리즘이다. 다만, 이와 같은 최적화 알고리즘은 본 발명의 다양한 실시예들에 따라 다양하게 선택될 수 있다. 예를 들어, 최적화 알고리즘은 Downhill simplex 알고리즘뿐만 아니라, Conjugate Gradient 알고리즘, Powell 알고리즘 등이 선택될 수도 있고, 복수의 최적화 알고리즘들을 함께 선택될 수도 있다. The integrated transformation
통합 변환 정보 최적화부(2212)는 볼륨 영상들간의 유사도들의 합이 최대화되도록 통합 변환 정보를 갱신하고, 갱신된 통합 변환 정보로부터 최적화 변환 함수를 생성한다. 예를 들어, 통합 변환 정보 최적화부(2212)는 복수의 변환 함수들로부터 추출된 파라미터들을 구성 요소로 하는 통합 변환 정보인 벡터 v를 볼륨 영상들간의 유사도들의 합이 최대화되도록 갱신하여 갱신된 통합 변환 정보인 v*를 생성하고, 갱신된 통합 변환 정보인 v*를 구성하는 파라미터들을 이용하여 최적화 변환 함수를 생성한다. 수학식 3 및 5를 통해 구체화하면, 통합 변환 정보 최적화부(2212)는 수학식 5의 S를 최대화시키도록 수학식 3의 v에 포함된 v2, 1 내지 vN,N-1각각을 갱신하여 v*를 생성하고, 생성된 v*를 구성하는 파라미터들을 이용하여 최적화 변환 함수를 생성하는 것이다. The integrated transformation
본 발명의 다른 실시예에 따르면, 볼륨 영상들간의 유사도는 볼륨 영상들 각각에 포함된 부분 영역들간의 유사도에 기초하여 결정될 수 있다. 이 경우, 부분 영역들 각각의 유사도는 부분 영역들 각각의 형태적인 특성의 유사도를 의미한다. 예를 들어, 제1 볼륨 영상과 제2 볼륨 영상간의 유사도는 제1 볼륨 영상에 포함된 제1 부분 영역과 제2 볼륨 영상에 포함된 제2 부분 영역간의 유사도로부터 결정되고, 이 때, 제1 부분 영역과 제2 부분 영역간의 유사도는 제1 부분 영역의 형태적인 특성과 제2 부분 영역의 형태적인 특성간의 유사도를 의미할 수 있다. 본 발명의 일 실시예에 따르면, 부분 영역들 각각의 형태적인 특성은 부분 영역들 각각을 구(sphere) 형상의 영역으로 평활화(normalize)함에 기초하여 결정된다. 이하에서는 본 발명의 일 실시예에 따라 부분 영역의 형태적인 특성을 결정하는 과정을 설명한다. According to another embodiment of the present invention, the similarity between the volume images can be determined based on the similarity between the partial regions included in each volume image. In this case, the degree of similarity of each of the partial regions means the similarity of the morphological characteristics of each of the partial regions. For example, the similarity between the first volume image and the second volume image is determined from the similarity between the first partial region included in the first volume image and the second partial region included in the second volume image, The similarity degree between the partial region and the second partial region may mean the similarity between the morphological characteristic of the first partial region and the morphological characteristic of the second partial region. According to one embodiment of the present invention, the morphological characteristics of each of the subregions are determined based on normalizing each of the subregions into a sphere-shaped region. Hereinafter, a process of determining a morphological characteristic of a partial region according to an embodiment of the present invention will be described.
도 6은 부분 영역의 형태적인 특성을 결정하는 과정을 설명하기 위한 도면이다. 이와 같은 과정은 통합 변환 정보 최적화부(2212)에 의해 수행된다. 도 6을 참조하면, 통합 변환 정보 최적화부(2212)는 복수의 볼륨 영상들 중 제1 볼륨 영상에 포함된 부분 영역(61)을 타원체(ellipsoid) 형상의 영역(62)으로 변환하고, 변환된 타원체 형상의 영역(62)을 구(sphere) 형상의 영역(63)으로 변환하고, 구(sphere) 형상의 영역(63)을 평활화(normalize)하여 평활화된(normalized) 영역(64)을 생성하고, 평활화된(normalized) 영역(64)에 기초하여 부분 영역(61)의 형태적인 특성을 결정할 수 있다. 이하에서 보다 구체적으로 설명한다. 6 is a diagram for explaining a process of determining the morphological characteristics of a partial area. This process is performed by the integrated conversion
도 6을 참조하면, 통합 변환 정보 최적화부(2212)는 부분 영역(61)을 타원체 형상의 영역(62)으로 변환한다. 일반적으로, 통합 변환 정보 최적화부(2212)는 부분 영역(61)에 포함된 복수의 복셀들 중 어느 하나를 타원체 형상의 영역(62)의 중심 복셀로 결정하고, 결정된 중심 복셀을 이용하여 타원체 형상의 영역(62)을 정의할 수 있다. 예를 들어, 통합 변환 정보 최적화부(2212)는 타원체 형상의 영역(62)을 부분 영역(61)에 대응하는 복셀들 각각인 x에 대하여 수학식 6과 같이 정의할 수 있다. 이 때, 부분 영역(61)에 대응하는 복셀들은 부분 영역(61)에 포함된 복셀들 또는 부분 영역(61)에 포함된 복셀과 부분 영역(61)의 주변의 복셀들을 모두를 의미할 수 있다. 또한, c는 타원체 형상의 영역(62)에 포함된 복셀들 중 중심 복셀을 의미하고, 
[수학식 6]&Quot; (6) "
한편, 부분 영역(61)은 일반적으로 볼륨 영상들 각각에 포함된 적어도 하나의 복셀로 구성된 소정 영역을 의미한다. 일반적으로, 이러한 부분 영역은 3차원적으로 표현되나 이에 한정되는 것은 아니다. 즉, 부분 영역은 2차원으로도 표현될 수 있다. 또한, 이러한 부분 영역에는 둘 이상의 복수의 복셀들이 포함된다. 예를 들어, 부분 영역은 볼륨 영상들 각각에 포함된 복셀들 중 20 개의 복셀들로 구성된 소정의 3차원 영역을 의미할 수 있다. 또한, 부분 영역은 둘 이상의 복수의 복셀들로 구성된 3차원 볼륨 세그먼트라 표현될 수도 있다. 또한, 부분 영역(61)은 일반적으로 볼륨 영상들 각각의 복셀들의 강도(intensity)에 기초하여 복수의 볼륨 영상들 각각으로부터 추출한다. 예를 들어, 부분 영역(61)은 볼륨 영상에 포함된 복수의 강도들을 비교함에 기초하여, 복수의 복셀들 중 상호간에 유사한 강도를 갖는 복셀들의 집합으로서 결정될 수 있다. 이와 같이, 상호간에 유사한 강도를 갖는 복셀들을 추출하는 일 실시예는 Maximally stable extremal regions 방식(J. Matas et al., "Robust wide baseline stereo from maximally stable extremal regions," BMVC 2002)을 3차원적으로 이용한 방식이 포함된다. 다만, 이러한 실시예는 본 발명에 따른 하나의 실시예에 불과하므로, 본원 발명의 권리 범위가 이러한 실시예에 한정 해석되지는 않는다. 예를 들어, 부분 영역(61)은 볼륨 영상에 포함된 복셀들 중 어느 하나의 복셀을 임의로 결정하고, 결정된 복셀의 주변 복셀들의 강도를 비교한 후, 유사한 강도를 갖는 복셀들의 집합으로부터 추출될 수도 있고, 볼륨 영상에 포함된 모든 복셀들 각각의 위치와 강도를 기준으로, 상호간에 유사한 강도를 갖는 복셀들의 집합으로부터 추출될 수도 있다. Meanwhile, the
도 6을 참조하면, 통합 변환 정보 최적화부(2212)는 타원체 형상의 영역(62)을 구 형상의 영역(63)으로 변환한다. 이 때, 구 형상 영역(63)은 앞서 설명된 타원체 형상의 영역(62)과 마찬가지로 부분 영역(61)에 대응하는 복셀들 각각인 x에 대한 수학식 7과 같이 정의한다. 구체적으로, 통합 변환 정보 최적화부(2212)는 구 형상의 영역(63)을 수학식 6에 포함된 공분산(covariance) 행렬의 역행렬인 
[수학식 7]&Quot; (7) "
도 6을 참조하면, 통합 변환 정보 최적화부(2212)는 구 형상의 영역(63)을 평활화하여 평활화된 구 형상의 영역(64)으로 변환한다. 이 때, 평활화된 구 형상의 영역(64)은 구 형상의 영역(63)으로부터 구의 중심 복셀을 기준으로 상호 직교하는 3 개의 벡터 성분들에 의해서 정의된다. 예를 들어, 통합 변환 정보 최적화부(2212)는 평활화된 구 형상의 영역(64)을 수학식 7에 로테이션 행렬 R을 적용하여 부분 영역(61)에 대응하는 복셀들 각각인 x에 대하여 수학식 8과 같이 정의한다. Referring to FIG. 6, the integrated transformation
[수학식 8]&Quot; (8) "
도 6을 참조하면, 통합 변환 정보 최적화부(2212)는 수학식 9를 이용하여 부분 영역(61)에 대응하는 복셀들을 변환한 후, 수학식 9에 의해 변환된 복셀들의 강도를 이용하여 로테이션 행렬 R을 결정할 수 있다. 이 때, 로테이션 행렬 R은 일반적으로 3차원을 구성하는 3 개의 벡터 성분을 구성요소로서 포함한다. 따라서, 통합 변환 정보 최적화부(2212)는 수학식 9를 이용하여 부분 영역(61)에 대응하는 복셀들을 변환한 후, 수학식 9에 의해 변환된 복셀들의 강도를 이용하여 강도의 변화도(gradient)의 크기가 큰 방향을 순차적으로 검출함으로써, 3 개의 벡터 성분을 결정할 수 있다. 예를 들어, 통합 변환 정보 최적화부(2212)는 수학식 9를 이용하여 부분 영역(61)에 대응하는 복셀들을 이동(wraping)하고, 강도의 변화도(gradient)의 크기에 비례하는 가중치로 변화도 방향의 히스토그램을 만든 후, 가장 빈도수가 높은 방향을 찾고 그 방향을 벡터 v1으로 결정하고, v1에 직교한 두 개의 방향 중 보다 빈도수가 높은 방향을 벡터 v2로 결정하고, v1과 v2에 모두 직교하는 방향을 벡터 v3로 결정함으로써, 로테이션 행렬 R을 결정할 수 있다. 수학식 10은 수학식 8의 로테이션 행렬 R을 나타낸다. Referring to FIG. 6, the integrated transformation
[수학식 9]&Quot; (9) "
[수학식 10]&Quot; (10) "
도 6을 참조하면, 통합 변환 정보 최적화부(2212)는 평활화된 구 형상의 영역(64)에 기초하여 부분 영역(61)의 형태적인 특성을 결정할 수 있다. 구체적으로, 통합 변환 정보 최적화부(2212)는 수학식 8을 이용하여 부분 영역(61)에 대응하는 복셀들을 변환한 후, 변환된 복셀들의 강도를 이용하여 평활화된 구 형상의 영역(64)의 일 영역별로 강도의 변화하는 정도를 나타내는 지표를 생성하여 이러한 지표들을 하나의 벡터로 집적(aggregation)시킨 형태적인 특성을 결정할 수 있다. 이와 같은 형태적인 특성은 불변하는 특징 식별자(Invariant Feature Descriptor)로 표현될 수 있다. 또한, 변화하는 정도를 나타내는 지표의 일 예에는 Intensity Gradient Orientation Histogram이 포함된다. Referring to FIG. 6, the integrated transformation
도 7은 평활화된 구 형상의 영역(64)의 일 영역에 포함된 복셀들의 강도들의 변화하는 정도를 나타내는 지표를 나타낸 도면이다. 도 6 및 7을 참조하면, 통합 변환 정보 최적화부(2212)는 부분 영역(61)에 대응하는 복셀들로부터 평활화된 구 형상의 영역(64)의 일 영역(631)에 포함된 복셀들의 강도의 변화하는 정도를 나타내는 지표(632)를 생성한다. 또한, 통합 변환 정보 최적화부(2212)는 평활화된 구 형상의 영역(64)의 다른 영역들 각각에 포함된 복셀들의 강도의 변화하는 정도를 나타내는 지표들을 생성하고, 생성된 복수의 지표들을 이용하여 부분 영역(61)에 대한 형태적인 특성을 결정할 수 있다. FIG. 7 is a diagram showing an index indicating the degree of variation of intensities of voxels included in one area of the smoothed
이와 같이, 통합 변환 정보 최적화부(2212)는 볼륨 영상들 각각의 부분 영역들로부터 볼륨 영상들 각각의 형태적인 특성을 결정할 수 있다. 또한, 통합 변환 정보 최적화부(2212)는 볼륨 영상들 각각의 부분 영역들로부터 결정된 볼륨 영상들 각각의 형태적인 특성을 비교하여, 볼륨 영상들간의 유사도를 결정할 수 있다. 볼륨 영상에 포함된 부분 영역들로부터 이러한 볼륨 영상의 형태적인 특성을 결정한다는 것은 볼륨 영상에 포함된 부분 영역들의 형태적인 특성의 평균을 볼륨 영상의 형태적인 특성으로 결정한다는 것을 의미할 수도 있고, 볼륨 영상에 포함된 적어도 하나의 부분 영역 각각의 형태적인 특성의 셋을 볼륨 영상의 형태적인 특성으로 결정한다는 것을 의미할 수도 있다. 이와 같이, 볼륨 영상의 형태적인 특성을 결정함에 있어서, 이러한 볼륨 영상에 포함된 부분 영역들의 형태적인 특성을 이용하는 실시예들은 다양하게 결정될 수 있다. 또한, 통합 변환 정보 최적화부(2212)는 부분 영역들간의 유사도에 기초하여 볼륨 영상들간의 유사도를 결정할 수도 있다. 예를 들어, 통합 변환 정보 최적화부(2212)는 볼륨 영상들 중 제1 볼륨 영상에 포함된 적어도 하나의 제1 부분 영역 각각과 볼륨 영상들 중 제2 볼륨 영상에 포함된 적어도 하나의 제2 부분 영역 각각간의 유사도에 기초하여 제1 볼륨 영상과 제2 볼륨 영상간의 유사도를 결정할 수 있다. In this manner, the integrated transformation
또한, 본 발명의 다른 실시예에 따르면, 통합 변환 정보 최적화부(2212)는 볼륨 영상들간의 목적 함수를 극소화시키도록 통합 변환 정보를 갱신하고, 갱신된 통합 변환 정보에 기초하여 최적화 변환 함수를 생성할 수도 있다. 예를 들어, 통합 변환 정보 최적화부(2212)는 볼륨 영상들 각각에 포함된 복셀들을 이용하여 볼륨 영상들 각각의 에지 응답(Edge Response)을 생성하고, 볼륨 영상들 각각의 에지 응답의 크기 및 방향의 유사도를 에지 응답의 크기 및 방향 차이로 대체하여 목적 함수를 만들고, 이와 같은 목적 함수에 기초하여 최적화 변환 함수를 생성할 수도 있다. According to another embodiment of the present invention, the integrated transformation
또한, 본 발명의 다른 실시예에 따르면, 통합 변환 정보 최적화부(2212)는 복수의 볼륨 영상들 중 일부의 볼륨 영상들간의 유사도를 최대화하도록 최적화 변환 함수를 결정할 수도 있다. 예를 들어, 통합 변환 정보 최적화부(2212)는 제1 볼륨 영상과 제2 볼륨 영상간의 제1 변환 함수, 제2 볼륨 영상과 제3 볼륨 영상간의 제2 변환 함수 및 제3 볼륨 영상과 제4 볼륨 영상간의 제3 변환 함수 각각으로부터 추출된 파라미터들에 의해 통합 변환 정보가 결정된 경우, 제1 변환 함수, 제2 변환 함수 및 제3 변환 함수로부터 추출된 파라미터를 모두 이용하여, 제1 볼륨 영상과 제2 볼륨 영상간의 유사도를 제외한 제2 볼륨 영상과 제3 볼륨 영상간의 유사도와 제3 볼륨 영상과 제4 볼륨 영상간의 유사도의 합이 최대화되도록 최적화 변환 함수를 결정할 수도 있다. 또한, 다른 예를 들어, 통합 변환 정보 최적화부(2212)는 제1 볼륨 영상과 제2 볼륨 영상간의 제1 변환 함수, 제2 볼륨 영상과 제3 볼륨 영상간의 제2 변환 함수 및 제3 볼륨 영상과 제4 볼륨 영상간의 제3 변환 함수 각각으로부터 추출된 파라미터들에 의해 통합 변환 정보가 결정된 경우, 제1 변환 함수로부터 추출된 파라미터를 제외한 제2 변환 함수와 제3 변환 함수로부터 추출된 파라미터들만을 이용하여, 제1 볼륨 영상과 제2 볼륨 영상간의 유사도, 제2 볼륨 영상과 제3 볼륨 영상간의 유사도 및 제3 볼륨 영상과 제4 볼륨 영상간의 유사도의 합이 최대화되도록 최적화 변환 함수를 결정할 수도 있다. Also, according to another embodiment of the present invention, the integrated transformation
앞서 설명된 바와 같이, 영상 프로세서(22)는 최적화 변환 함수에 기초하여 볼륨 파노라마 영상을 생성한다. 구체적으로, 영상 프로세서(22)는 최적화 변환 함수에 기초하여 복수의 볼륨 영상들의 영상 데이터로부터 볼륨 파노라마 영상을 나타내는 영상 데이터를 생성한다. 이하에서는 영상 프로세서(22)에 포함된 합성 영상 데이터 생성부(222) 및 볼륨 파노라마 영상 생성부(223)에 의한 동작 설명을 통해 볼륨 파노라마 영상을 생성하는 과정을 더욱 상세하게 설명하도록 한다. As described above, the
합성 영상 데이터 생성부(222)는 복수의 변환 함수들 각각으로부터 생성된 최적화 변환 함수에 기초하여 복수의 볼륨 영상들의 영상 데이터들로부터 합성되는 볼륨 영상들의 영상 데이터들을 생성한다. 도 3을 통해 예시하면, 합성 영상 데이터 생성부(222)는 제1 볼륨 영상(31)과 제2 볼륨 영상(32)간의 제1 변환 함수로부터 생성된 제1 최적화 변환 함수에 기초하여 제2 볼륨 영상(32)의 영상 데이터로부터 제1 볼륨 영상(31)과 합성되는 볼륨 영상(33)의 영상 데이터를 생성한다. 이 때, 볼륨 영상(33)은 제2 볼륨 영상(32)에 제1 최적화 변환 함수가 반영된 것을 의미한다. 따라서, 볼륨 영상(33)은 제1 볼륨 영상(31)을 기준으로 제2 볼륨 영상(32)을 제1 볼륨 영상(31)에 정합시킨 영상을 의미할 수 있다. 또한, 합성 영상 데이터 생성부(222)는 제2 볼륨 영상과 제3 볼륨 영상간의 제2 변환 함수로부터 생성된 제2 최적화 변환 함수에 기초하여 제3 볼륨 영상의 영상 데이터로부터 제2 볼륨 영상과 합성되는 다른 볼륨 영상의 영상 데이터를 생성할 수 있다. 또한, 합성 영상 데이터 생성부(222)는 일반적으로 최적화 변환 함수에 기초하여 제1 볼륨 영상에 포함된 복셀들을 기준으로 제2 볼륨 영상에 포함된 복셀들을 이동하여 합성되는 볼륨 영상의 복셀들을 생성한다. 다만, 이에 한정되지는 않는다. The composite image
도 8은 도 2의 합성 영상 데이터 생성부(222)에서 합성되는 볼륨 영상의 영상 데이터를 생성하는 과정을 나타낸 동작 흐름도이다. 도 8의 단계 81 내지 단계 83은 합성 영상 데이터 생성부(222)에 의해 수행된다. 단계 81에서 합성 영상 데이터 생성부(222)는 전역 변환 함수에 기초하여 제1 볼륨 영상과 합성되는 볼륨 영상의 영상 데이터를 생성한다. 단계 82에서 합성 영상 데이터 생성부(222)는 합성되는 볼륨 영상으로부터 분할된 적어도 하나의 지역 볼륨 영상에 기초하여 지역 변환 함수를 결정한다. 구체적으로, 합성 영상 데이터 생성부(222)는 합성되는 볼륨 영상을 복수개의 지역 볼륨 영상들로 분할하고, 분할된 지역 볼륨 영상들 각각에 기초하여 지역 볼륨 영상들 각각에 대한 지역 변환 함수를 결정한다. FIG. 8 is a flowchart illustrating a process of generating image data of a volume image synthesized by the composite image
도 9는 합성되는 볼륨 영상으로부터 적어도 하나의 지역 볼륨 영상을 분할하는 일 실시예를 나타낸 도면이다. 도 9를 통해 예시하면, 합성 영상 데이터 생성부(222)는 합성되는 볼륨 영상(91)을 복수개의 지연 볼륨 영상들(92)로 분할하고, 분할된 지연 볼륨 영상들(92) 각각에 대한 지역 변환 함수를 결정한다. 본 발명의 일 실시예에 따르면, 합성 영상 데이터 생성부(222)는 지역 볼륨 영상들(92) 각각에 대응하는 복셀들과 제2 볼륨 영상들에 대응하는 복셀들간의 변환 함수에 기초하여 지역 볼륨 영상들(92) 각각에 대한 지역 변환 함수를 결정한다. 이 때, 합성 영상 데이터 생성부(222)는 지역 볼륨 영상들(92) 각각에 대응하는 복셀들과 제2 볼륨 영상들에 대응하는 복셀들간의 변환 함수를 초기값으로 하여 최적화 알고리즘을 적용함으로써, 지역 볼륨 영상들(92) 각각에 대한 지역 변환 함수를 결정할 수 있다. 예를 들어, 합성 영상 데이터 생성부(222)는 지역 볼륨 영상들(92) 각각의 초기 지역 변환 함수 (I, O)와 초기 지역 변환 함수의 주변에서 샘플링 한 지역 변혼 특성에 기초하여 최적화 알고리즘을 적용할 수 있다. 이 때, I는 3행 3열의 단위 행렬을 의미하고, O는 3 차원 0 벡터를 의미한다. 9 is a view showing an embodiment of dividing at least one local volume image from a synthesized volume image. 9, the composite image
도 9를 참조하면, 본 발명의 일 실시예에 따른 합성 영상 데이터 생성부(222)는 지역 볼륨 영상들(92) 중 어느 하나의 지역 볼륨 영상의 지역 변환 함수를 결정하기 위해, 지역 볼륨 영상들(92) 중 어느 하나의 지역 볼륨 영상의 주변 지역 볼륨 영상의 지역 변환 함수를 이용할 수도 있다. 예를 들어, 합성 영상 데이터 생성부(222)는 지역 볼륨 영상들(92) 중 어느 하나의 지역 볼륨 영상의 주변 지역 볼륨 영상의 지역 변환 함수를 이용하여 지역 볼륨 영상들(92) 중 어느 하나의 지역 볼륨 영상의 지역 변환 함수에 대한 보간(Interpolation)을 수행할 수도 있다. 9, the composite image
도 9를 참조하면, 본 발명의 일 실시예에 따른 합성 영상 데이터 생성부(222)는 합성되는 볼륨 영상들(91)을 계층적으로 분할한다. 예를 들어, 합성 영상 데이터 생성부(222)는 합성되는 볼륨 영상들(91)을 4 개의 영역으로 분할하여 지역 볼륨 영상들(92)을 생성하고, 생성된 지역 볼륨 영상들(92) 각각을 분할하여 지역 볼륨 영상들(93)을 생성할 수 있다. 일반적으로, 지역 볼륨 영상들 중 텍스쳐가 많은 지역 볼륨 영상은 보다 작은 영역으로 분할할수록 보다 정확한 지연 변환 함수가 획득된다. 따라서, 본 발명의 일 실시예에 따른 합성 영상 데이터 생성부(222)는 지역 볼륨 영상들 각각에 포함된 텍스쳐의 양을 고려하여 지역 볼륨 영상들 각각마다 얼마나 작은 영역으로 분할할지를 적응적으로 결정할 수도 있다. Referring to FIG. 9, the composite image
단계 83에서 합성 영상 데이터 생성부(222)는 결정된 지역 변환 함수에 기초하여 제1 볼륨 영상과 합성되는 볼륨 영상의 영상 데이터를 갱신한다. 구체적으로, 합성 영상 데이터 생성부(222)는 합성되는 볼륨 영상으로부터 분할된 지역 볼륨 영상들 각각에 지역 볼륨 영상들 각각의 지역 변환 함수를 적용하여, 합성되는 볼륨 영상의 영상 데이터를 갱신할 수 있다. In
볼륨 파노라마 영상 생성부(223)는 복수의 볼륨 영상들의 영상 데이터와 합성되는 볼륨 영상들의 영상 데이터에 기초하여 볼륨 파노라마 영상을 나타내는 영상 데이터를 생성한다. 도 4를 통해 예시하면, 볼륨 파노라마 영상 생성부(223)는 제1 최적화 변환 함수에 기초하여 제2 볼륨 영상(42)의 영상 데이터로부터 생성된 제1 볼륨 영상(41)과 합성되는 볼륨 영상의 영상 데이터, 제2 최적화 변환 함수에 기초하여 제3 볼륨 영상(43)의 영상 데이터로부터 생성된 제2 볼륨 영상(42)과 합성되는 다른 볼륨 영상의 영상 데이터 및 제1 볼륨 영상(41)의 영상 데이터에 기초하여 볼륨 파노라마 영상(44)을 나타내는 영상 데이터를 생성한다. The volume panorama
볼륨 파노라마 영상 생성부(223)는 제1 볼륨 영상에 포함된 복셀들, 합성되는 볼륨 영상에 포함된 복셀들 및 합성되는 다른 볼륨 영상에 포함된 복셀들을 합성하여 볼륨 파노라마 영상을 생성한다. 일반적으로, 최적화 변환 함수에 기초하여 제1 볼륨 영상으로부터 생성된 합성되는 볼륨 영상에 포함된 복셀들 각각은 제1 볼륨 영상에 포함된 복셀들 각각에 대응한다. 그럼에도 불구하고, 합성되는 볼륨 영상에 포함된 복셀들 각각의 강도와 합성되는 볼륨 영상에 포함된 복셀들 각각에 대응하는 제1 볼륨 영상에 포함된 복셀들 각각의 강도는 다를 수 있다. 이러한 강도의 차이는 일반적으로 초음파 신호의 쉐도우(shadow) 효과에 의해서 나타날 수 있다. 이 경우, 볼륨 파노라마 영상 생성부(223)는 제1 볼륨 영상의 복셀들 중 어느 하나의 강도, 합성되는 볼륨 영상의 복셀들 중 어느 하나의 강도 및 합성되는 다른 볼륨 영상의 복셀들 중 어느 하나의 강도에 기초하여 볼륨 파노라마 영상의 복셀들 중 어느 하나의 강도를 결정할 수 있다. 예를 들어, 볼륨 파노라마 영상 생성부(223)는 제1 볼륨 영상의 복셀들 중 어느 하나의 강도, 합성되는 볼륨 영상의 복셀들 중 어느 하나의 강도 및 합성되는 다른 볼륨 영상의 복셀들 중 어느 하나의 강도 중 가장 낮은 강도 또는 가장 큰 강도를 볼륨 파노라마 영상의 복셀들 중 어느 하나의 강도로 결정할 수 있다. 또한, 본 발명의 일 실시예에 따르면, 볼륨 파노라마 영상 생성부(223)는 제1 볼륨 영상의 복셀들 중 어느 하나의 강도, 합성되는 볼륨 영상의 복셀들 중 어느 하나의 강도 및 합성되는 다른 볼륨 영상의 복셀들 중 어느 하나의 강도의 평균을 볼륨 파노라마 영상의 복셀들 중 어느 하나의 강도로 결정할 수도 있다. The volume panorama
도 10은 발명의 다른 실시예에 따른 볼륨 파노라마 영상 생성 장치(100)의 구성도이다. 도 10의 볼륨 파노라마 영상 생성 장치(100)의 영상 프로세서(22)는 도 2의 영상 프로세서(22)와 비교하여 변환 함수 결정부(224)를 더 포함한다. 이와 같은 변환 함수 결정부(224)는 입력부(21)로부터 입력된 볼륨 영상들의 영상 데이터에 기초하여 볼륨 영상들간의 변환 함수를 생성하고, 생성된 변환 함수를 최적화 변환 함수 생성부(221)로 전달한다. 변환 함수 결정부(224)에 의해 생성되는 변환 함수들은 앞서 설명된 복수의 볼륨 영상들간의 복수의 변환 함수들에 대해 설명된 것에 의한다. FIG. 10 is a block diagram of a volume panoramic
본 발명의 일 실시예에 따르면, 변환 함수 결정부(224)는 볼륨 영상들 각각에 포함된 부분 영역들간의 부분 변환 함수에 기초하여 볼륨 영상들간의 변환 함수를 결정할 수 있다. 예를 들어, 변환 함수 결정부(224)는 복수의 볼륨 영상들 중 제1 볼륨 영상의 제1 부분 영역과 복수의 볼륨 영상들 중 제2 볼륨 영상의 제2 부분 영역간의 부분 변환 함수를 결정하고, 결정된 제1 부분 영역과 제2 부분 영역간의 부분 변환 함수에 기초하여 제1 볼륨 영상과 제2 볼륨 영상간의 변환 함수를 결정할 수 있다. 이하에서 구체적으로 설명한다. According to an embodiment of the present invention, the transform
변환 함수 결정부(224)는 제1 볼륨 영상의 제1 부분 영역을 구 형상의 영역으로 평활화하여 평활화된 구 형상의 영역을 생성한다. 이 때, 앞서 설명된 바와 같이, 변환 함수 결정부(224)는 수학식 6 내지 10을 이용하여, 제1 볼륨 영상의 제1 부분 영역을 타원체 형상의 영역으로 변환하고, 변환된 타원체 형상의 영역을 구 형상의 영역으로 변환하고, 변환된 구 형상의 영역을 평활화하여 평활화된 구 형상의 영역을 생성할 수 있다. 이와 같은 맥락으로, 변환 함수 결정부(224)는 제2 볼륨 영상의 제2 부분 영역을 구 형상의 영역으로 평활화하여 평활화된 구 형상의 영역을 생성한다. The transform
변환 함수 결정부(224)는 제1 부분 영역과 제2 부분 영역 각각을 구 형상의 영역으로 변환하는 적어도 하나의 파라미터를 결정하고, 결정된 파라미터에 기초하여 제1 부분 영역과 제2 부분 영역간의 부분 변환 함수를 결정한다. 예를 들어, 변환 함수 결정부(224)는 제1 부분 영역에 대응하는 복셀들 각각인 x1에 대하여 수학식 8로부터 변형된 수학식 11과 같이 정의할 수 있다. 이 때, c1는 타원체 형상의 영역에 포함된 복셀들 중 중심 복셀을 의미하고, 
[수학식 11]&Quot; (11) "
[수학식 12]&Quot; (12) "
변환 함수 결정부(224)는 결정된 파라미터에 기초하여 제1 부분 영역과 제2 부분 영역간의 부분 변환 함수를 결정한다. 이 때, 결정된 파라미터는 제1 부분 영역에 대한 제1 파라미터와 제2 부분 영역에 대한 제2 파라미터를 포함한다. 또한, 앞서 설명된 바와 같이, 제1 파라미터에는 제1 부분 영역에 대응하는 복셀들의 위치 이동을 나타내는 제1 파라미터 및 제1 부분 영역에 대응하는 복셀들의 방향 변환을 나타내는 제1 파라미터 중 적어도 하나가 포함되고, 제2 파라미터에는 제2 부분 영역에 대응하는 복셀들의 위치 이동을 나타내는 제2 파라미터 및 제2 부분 영역에 대응하는 복셀들의 방향 변환을 나타내는 제2 파라미터 중 적어도 하나가 포함된다. 또한, 위치 이동을 나타내는 제1 파라미터 및 제2 파라미터 각각은 공분산 행렬을 의미하고, 방향 변환을 나타내는 제1 파라미터 및 제2 파라미터 각각은 로테이션 행렬을 의미할 수 있다. 또한, 본 발명의 일 실시예에 따르면, 제1 파라미터와 제2 파라미터는 앞서 설명된 통합 변환 정보를 생성하기 위해 변환 함수들로부터 추출된 파라미터에 대응한다. The conversion
변환 함수 결정부(224)는 제1 파라미터와 제2 파라미터에 기초하여 제1 부분 영역과 제2 부분 영역간의 부분 변환 함수를 결정한다. 예를 들어, 변환 함수 결정부(224)는 수학식 11 및 수학식 12로부터 변형된 수학식 13과 같이 정의할 수 있다. 수학식 13을 참조하면, 제1 부분 영역과 제2 부분 영역간의 변환 함수는 제2 부분 영역에 대응하는 복셀들 각각을 나타내는 x2가 제1 부분 영역에 대응하는 복셀들 각각을 나타내는 x1으로 변환되는 관계로 정의될 수 있다. The conversion
[수학식 13]&Quot; (13) "
변환 함수 결정부(224)는 제1 부분 영역과 제2 부분 영역간의 부분 변환 함수에 기초하여 제1 볼륨 영상과 제2 볼륨 영상간의 변환 함수를 결정한다. 이 때, 일반적으로 제1 볼륨 영상과 제2 볼륨 영상간의 변환 함수는 제1 볼륨 영상에 대응하는 복셀들과 제2 볼륨 영상에 대응하는 복셀들간의 변환 함수를 의미할 수 있다. 이 때, 제1 볼륨 영상에 대응하는 복셀들과 제2 볼륨 영상에 대응하는 복셀들간의 변환 함수는 제1 볼륨 영상에 대응하는 복셀들을 기준으로 제2 볼륨 영상에 대응하는 복셀들을 제1 볼륨 영상에 대응하는 복셀들로 정합시키기 위한 제2 볼륨 영상에 대응하는 복셀들의 변환 함수를 의미한다. 제1 볼륨 영상에 대응하는 복셀들은 제1 볼륨 영상에 포함된 복셀들을 의미할 수 있다. 이와 같은 맥락으로, 제2 볼륨 영상에 대응하는 복셀들은 제2 볼륨 영상에 포함된 복셀들을 의미할 수 있다. 다만, 이에 한정되지 않는다. 예를 들어, 제1 볼륨 영상에 대응하는 복셀들은 제1 볼륨 영상에 포함된 복셀들 중 임계값 이상의 강도를 갖는 복셀들만을 의미할 수도 있다. 따라서, 제1 부분 영역과 제2 부분 영역간의 변환 함수에 기초하여 제1 볼륨 영상과 제2 볼륨 영상간의 변환 함수를 결정하는 것은 제1 부분 영역과 제2 부분 영역간의 변환 함수를 제1 볼륨 영상에 포함된 복셀들로부터 제2 볼륨 영상에 포함된 복셀들로의 변환을 의미할 수 있다. 예를 들어, 변환 함수 결정부(224)는 제1 부분 영역과 제2 부분 영역간의 변환 함수를 나타내는 수학식 13을 이용하여 제1 볼륨 영상에 포함된 복셀들로부터 제2 볼륨 영상에 포함된 복셀들로의 변환을 수행할 수 있다. The conversion
본 발명의 일 실시예에 따르면, 변환 함수 결정부(224)는 복수의 대응쌍들을 이루는 복수의 제1 부분 영역들 각각과 복수의 제2 부분 영역들 각각에 대한 부분 변환 함수를 결정한다. 예를 들어, 변환 함수 결정부(224)는 복수의 제1 부분 영역들 중 어느 하나의 제1 부분 영역과 이러한 어느 하나의 제1 부분 영역에 대응하는 제2 부분 영역들 중 다른 하나의 제2 부분 영역간의 제1 부분 변환 함수를 결정하고, 복수의 제1 부분 영역들 중 다른 하나의 제1 부분 영역과 이러한 다른 하나의 제1 부분 영역에 대응하는 제2 부분 영역들 중 다른 하나의 제2 부분 영역간의 제2 부분 변환 함수를 결정할 수 있다. 또한, 변환 함수 결정부(224)는 복수의 부분 변환 함수들에 기초하여 제1 볼륨 영상과 제2 볼륨 영상간의 변환 함수를 결정할 수도 있다. 또한, 변환 함수 결정부(224)는 복수의 부분 변환 함수들 중 적어도 하나의 부분 변환 함수를 선택하고, 선택된 부분 변환 함수에 기초하여 제1 볼륨 영상과 제2 볼륨 영상간의 변환 함수를 결정할 수도 있다. According to one embodiment of the present invention, the transform
변환 함수 결정부(224)는 복수의 부분 변환 함수들 중 적어도 하나를 선택함에 있어서, 부분 변환 함수들 각각에 기초하여 제2 볼륨 영상을 제1 볼륨 영상을 기준으로 이동(wraping)하고, 이동한 결과들을 비교하여 적어도 하나의 부분 변환 함수를 선택할 수 있다. 예를 들어, 변환 함수 결정부(224)는 제1 부분 변환 함수에 기초하여 제2 볼륨 영상을 이동시킨 결과와 제2 부분 변환 함수에 기초하여 제2 볼륨 영상을 이동시킨 결과를 비교한 후, 비교 결과에 따라 제1 부분 변환 함수 및 제2 부분 변환 함수 중 어느 하나를 선택할 수 있다. 일반적으로, 변환 함수 결정부(224)는 제1 부분 변환 함수에 기초하여 제2 볼륨 영상을 이동시킨 결과와 제2 부분 변환 함수에 기초하여 제2 볼륨 영상을 이동시킨 결과를 비교함에 있어서, 볼륨 간 유사도를 이용한다. 이 때, 볼륨 간 유사도는 제1 부분 변환 함수에 기초하여 제2 볼륨 영상을 이동시킨 결과와 제1 볼륨 영상간의 유사도를 의미한다. 이와 마찬가지로, 볼륨 간 유사도는 제2 부분 변환 함수에 기초하여 제2 볼륨 영상을 이동시킨 결과와 제1 볼륨 영상간의 유사도를 의미할 수도 있다. 따라서, 변환 함수 결정부(224)는 제1 부분 변환 함수에 기초하여 제2 볼륨 영상을 이동시킨 결과와 제1 볼륨간의 제1 유사도를 계산하고, 제2 부분 변환 함수에 기초하여 제2 볼륨 영상을 이동시킨 결과와 제1 볼륨간의 제2 유사도를 계산한 후, 제1 유사도와 제2 유사도를 비교한 결과 보다 높은 유사도를 갖는 제1 유사도를 선택하고, 제1 유사도에 대응하는 제1 부분 변환 함수를 선택할 수 있다. 이 때, 볼륨 간 유사도는 제1 볼륨 영상과 이동한 제2 볼륨 영상간의 복셀들의 강도의 분포 유사도, 동일 위치에 대응하는 복셀들간의 강도의 변화(gradient)의 크기 및 방향의 유사도가 이용될 수 있다. 이 때, 복셀들의 강도의 분포 유사도의 일 예는 Normalized Mutual Information이다. In selecting at least one of the plurality of partial conversion functions, the conversion
변환 함수 결정부(224)는 제1 부분 영역과 제2 부분 영역간의 부분 변환 함수에 기초하여 제1 볼륨 영상과 제2 볼륨 영상간의 변환 함수를 결정한다. 이 때, 제1 부분 영역과 제2 부분 영역간의 부분 변환 함수는 앞서 설명된 바와 같이 복수의 유사도들에 기초하여 복수의 변환 함수들 중 선택된 적어도 하나의 부분 변환 함수를 의미할 수도 있다. 예를 들어, 변환 함수 결정부(224)는 복수의 부분 변환 함수들 중 M 개의 부분 변환 함수를 선택하고, 이와 같은 M 개의 부분 변환 함수들에 최적화 알고리즘을 적용하여, 제1 볼륨 영상과 제2 볼륨 영상간의 유사도를 최대화하는 부분 변환 함수를 결정할 수 있다. 이 때, 최적화 알고리즘의 일 예는 Downhill Simplex이다. 다만, 이와 같은 최적화 알고리즘은 본 발명의 다양한 실시예들에 따라 다양하게 선택될 수 있다. 예를 들어, 최적화 알고리즘은 Downhill simplex 알고리즘뿐만 아니라, Conjugate Gradient 알고리즘, Powell 알고리즘 등이 선택될 수도 있고, 복수의 최적화 알고리즘들을 함께 선택될 수도 있다. 또한, 본 발명의 일 실시예에 따르면, 변환 함수 결정부(224)는 복수의 부분 변환 함수들로부터 선택된 부분 변환 함수들이 N (M>N)개 존재하는 경우, 부족한 L 개(L=M-N)의 부분 변환 함수들은 N 개의 부분 변환 함수들 각각의 주변에서 샘플링하여 생성할 수 있다. The conversion
본 발명의 일 실시예에 따르면, 변환 함수 결정부(224)는 최적화 알고리즘의 적용 없이, 복수의 부분 변환 함수들 중 적어도 하나의 부분 변환 함수를 그대로 이용하여 제1 볼륨 영상과 제2 볼륨 영상간의 변환 함수를 결정할 수 있다. 예를 들어, 변환 함수 결정부(224)는 제1 부분 영역과 제2 부분 영역간의 부분 변환 함수를 나타내는 수학식 13을 변환 함수로 결정할 수도 있다. According to an embodiment of the present invention, the transform
본 발명의 일 실시예에 따르면, 변환 함수 결정부(224)는 결정된 변환 함수에 대한 미세 보정(Refinement)을 수행한다. 이를 위해, 변환 함수 결정부(224)는 결정된 변환 함수를 제2 볼륨 영상에 적용하고, 적용된 제2 볼륨 영상과 제2 볼륨 영상간의 변환 함수를 샘플링한 후, 샘플링 한 변환 함수에 최적화 알고리즘을 다시 적용함으로써, 변환 함수에 대한 미세 보정을 수행할 수 있다. 이와 같이 미세 보정은 변환 함수를 갱신하는 것을 의미한다. According to an embodiment of the present invention, the transform
이와 같은 도 10의 볼륨 파노라마 영상 생성 장치(100)에 대하여 설명되지 아니한 사항은 앞서 도 2의 볼륨 파노라마 영상 생성 장치(20)에 대해 설명된 내용과 동일 및 설명된 내용으로부터 당업자에 의해 용이하게 유추 가능한 것으로 이하 설명을 생략하도록 한다. 10 are the same as those described above for the volume panorama
도 11은 본 발명의 일 실시예에 따른 볼륨 파노라마 영상 생성 방법의 동작 흐름도이다. 도 11에 도시된 실시예에 따른 볼륨 파노라마 영상 생성 방법은 도 2에 도시된 볼륨 파노라마 영상 생성 장치(20)에서 시계열적으로 처리되는 단계들로 구성된다. 따라서, 이하 생략된 내용이라고 하더라도 도 2에 도시된 볼륨 파노라마 영상 생성 장치(20)에 관하여 이상에서 기술된 내용은 도 11에 도시된 실시예에 따른 볼륨 파노라마 영상 생성 방법에도 적용된다. 11 is a flowchart illustrating an operation of a volume panorama image generating method according to an embodiment of the present invention. The volume panorama image generating method according to the embodiment shown in FIG. 11 is composed of the steps of time-series processing in the volume panorama
단계 111에서 입력부(21)는 볼륨 영상들 중 제1 볼륨 영상과 제1 볼륨 영상과 공통된 영역을 갖는 제2 볼륨 영상간의 변환 관계를 나타내는 변환 함수를 입력받는다. 단계 112에서 영상 프로세서(22)는 복수의 변환 함수들을 통합적으로 고려하여 복수의 변환 함수들 각각으로부터 최적화 변환 함수를 생성한다. 단계 113에서 영상 프로세서(22)는 최적화 변환 함수에 기초하여 볼륨 파노라마 영상을 생성한다. In
도 12는 본 발명의 다른 실시예에 따른 볼륨 파노라마 영상 생성 방법의 동작 흐름도이다. 도 12에 도시된 실시예에 따른 볼륨 파노라마 영상 생성 방법은 도 10에 도시된 볼륨 파노라마 영상 생성 장치(100)에서 시계열적으로 처리되는 단계들로 구성된다. 따라서, 이하 생략된 내용이라고 하더라도 도 10에 도시된 볼륨 파노라마 영상 생성 장치(100)에 관하여 이상에서 기술된 내용은 도 12에 도시된 실시예에 따른 볼륨 파노라마 영상 생성 방법에도 적용된다. 12 is a flowchart illustrating an operation of a volume panorama image generating method according to another embodiment of the present invention. The volume panorama image generation method according to the embodiment shown in FIG. 12 is composed of the steps of time-series processing in the volume panorama
단계 111에서 입력부(21)는 볼륨 영상들의 영상 데이터들을 입력받는다. 단계 112에서 영상 프로세서(22)는 볼륨 영상들의 영상 데이터들에 기초하여 볼륨 영상들 중 제1 볼륨 영상과 상기 제1 볼륨 영상과 공통된 영역을 갖는 제2 볼륨 영상간의 변환 관계를 나타내는 변환 함수를 결정한다. 단계 113에서 영상 프로세서(22)는 복수의 변환 함수들을 통합적으로 고려하여 복수의 변환 함수들 각각으로부터 최적화 변환 함수를 생성한다. 단계 114에서 영상 프로세서(22)는 최적화 변환 함수에 기초하여 볼륨 파노라마 영상을 생성한다. In
도 11 및 12를 각각을 통해 설명된 실시예에 따른 볼륨 파노라마 영상 생성 방법은 컴퓨터에서 실행될 수 있는 프로그램으로 작성 가능하고, 컴퓨터로 읽을 수 있는 기록매체를 이용하여 상기 프로그램을 동작시키는 범용 디지털 컴퓨터에서 구현될 수 있다. 상기 컴퓨터로 읽을 수 있는 기록매체는 마그네틱 저장매체(예를 들면, 롬, 플로피 디스크, 하드 디스크 등), 광학적 판독 매체(예를 들면, 시디롬, 디브이디 등)와 같은 저장매체를 포함한다. 11 and 12, a volume panorama image generating method according to the embodiment described above can be realized by a general-purpose digital computer that can be created as a program that can be executed by a computer and operates the program using a computer- Can be implemented. The computer-readable recording medium includes a storage medium such as a magnetic storage medium (e.g., ROM, floppy disk, hard disk, etc.), optical reading medium (e.g., CD ROM,
이제까지 본 발명에 대하여 그 바람직한 실시예들을 중심으로 살펴보았다. 본 발명이 속하는 기술 분야에서 통상의 지식을 가진 자는 본 발명이 본 발명의 본질적인 특성에서 벗어나지 않는 범위에서 변형된 형태로 구현될 수 있음을 이해할 수 있을 것이다. 그러므로 개시된 실시예들은 한정적인 관점이 아니라 설명적인 관점에서 고려되어야 한다. 본 발명의 범위는 전술한 설명이 아니라 특허청구범위에 나타나 있으며, 그와 동등한 범위 내에 있는 모든 차이점은 본 발명에 포함된 것으로 해석되어야 할 것이다. The present invention has been described with reference to the preferred embodiments. It will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention as defined by the appended claims. Therefore, the disclosed embodiments should be considered in an illustrative rather than a restrictive sense. The scope of the present invention is defined by the appended claims rather than by the foregoing description, and all differences within the scope of equivalents thereof should be construed as being included in the present invention.
20 ... 볼륨 파노라마 영상 생성 장치
22 ... 영상 프로세서
221 ... 최적화 변환 함수 생성부
222 ... 합성 영상 데이터 생성부
223 ... 볼륨 파노라마 영상 생성부20 ... Volume panoramic image generation device
22 ... image processor
221 ... optimization conversion function generation unit
222 ... composite image data generation unit
223 ... volume panorama image generating section
Claims (19)
상기 복수의 변환 함수들로부터 상기 복수의 3차원 볼륨 영상들 사이의 유사도에 기초하여 최적화 변환 함수를 생성하고, 상기 최적화 변환 함수를 이용하여 변환된 복수의 3차원 볼륨 영상들을 이용하여 볼륨 파노라마 영상을 생성하는 하나 이상의 프로세서; 및
상기 볼륨 파노라마 영상을 출력하는 출력부를 포함하고,
상기 유사도는, 상기 복수의 3차원 볼륨 영상들 중 제1 볼륨 영상의 형태적인 특성과, 상기 획득된 변환 함수들 중 하나에 기초하여 이동된 제2 볼륨 영상의 형태적인 특성 간의 유사도를 포함하고,
상기 하나 이상의 프로세서는, 상기 순차적인 볼륨 영상들의 대응하는 복셀들 각각의 위치 정보 및 강도(intensity)에 기초하여 상기 유사도를 결정하는, 볼륨 파노라마 생성 장치.Dimensional volume image of a plurality of three-dimensional volume images, which is a function for matching at least two of the plurality of three-dimensional volume images and a plurality of three-dimensional volume images of a sequential order, from an external apparatus that generates three- An input unit for receiving a plurality of conversion functions between the plurality of conversion functions;
Dimensional volume images from the plurality of transformation functions, and generates a volume panorama image using a plurality of transformed three-dimensional volume images using the optimization transformation function One or more processors; And
And an output unit for outputting the volume panorama image,
Wherein the degree of similarity includes a degree of similarity between morphological characteristics of a first volume image of the plurality of three-dimensional volume images and morphological characteristics of a second volume image moved based on one of the obtained conversion functions,
Wherein the one or more processors determine the similarity based on location information and intensity of each of the corresponding voxels of the sequential volume images.
상기 입력부는, 상기 외부 장치로부터 상기 복수의 3차원 볼륨 영상들과 함께 프로브의 위치 데이터를 수신하는, 볼륨 파노라마 생성 장치.The method according to claim 1,
Wherein the input unit receives the position data of the probe together with the plurality of three-dimensional volume images from the external apparatus.
상기 하나 이상의 프로세서는, 복수의 변환 함수들에 기초하여, 상기 순차적인 볼륨 영상들의 형태적인 특성의 유사도를 결정하고, 상기 결정된 유사도에 기초하여 상기 최적화 변환 함수를 생성하는, 볼륨 파노라마 생성 장치.The method according to claim 1,
Wherein the one or more processors determine a similarity of the morphological characteristics of the sequential volume images based on the plurality of conversion functions and generate the optimization conversion function based on the determined similarity.
상기 유사도들의 합을 최대화하도록 상기 복수의 변환 함수들을 변경하고, 상기 변경된 변환 함수들에 기초하여 상기 최적화 변환 함수를 생성하는, 볼륨 파노라마 생성 장치.2. The apparatus of claim 1,
Change the plurality of transformation functions to maximize a sum of the similarities, and generate the optimization transform function based on the modified transformation functions.
상기 생성된 최적화 변환 함수에 기초하여 제 2 볼륨 영상의 영상 데이터로부터 제 1 볼륨 영상과 합성되는 볼륨 영상의 영상 데이터를 생성하고, 상기 생성된 상기 제 1 볼륨 영상과 합성되는 볼륨 영상의 영상 데이터와 상기 제 1 볼륨 영상의 영상 데이터를 합성하는, 볼륨 파노라마 생성 장치.2. The apparatus of claim 1,
Generating a volume image of the volume image to be synthesized with the first volume image from the image data of the second volume image based on the generated optimization conversion function and generating image data of the volume image synthesized with the generated first volume image And synthesizes the image data of the first volume image.
상기 제 1 볼륨 영상과 합성된 볼륨 영상으로부터 분할된 지역 볼륨 영상들에 기초하여 지역 변환 함수를 결정하고, 상기 결정된 지역 변환 함수에 기초하여 상기 생성된 상기 제 1 볼륨 영상과 합성되는 볼륨 영상의 영상 데이터를 갱신하는, 볼륨 파노라마 생성 장치.7. The apparatus of claim 6,
Determining a local transformation function based on the local volume images segmented from the volume images synthesized with the first volume image, and generating a volume image of the volume image synthesized with the generated first volume image based on the determined local transformation function A volume panorama generator for updating data.
상기 복수의 변환 함수들 중 적어도 하나의 변환 함수는 제1 볼륨 영상의 부분 영역과 제2 볼륨 영상의 부분 영역간의 변환 관계를 나타내는, 볼륨 파노라마 생성 장치.The method according to claim 1,
Wherein at least one transformation function of the plurality of transformation functions represents a transformation relationship between a partial region of the first volume image and a partial region of the second volume image.
상기 획득된 변환 함수들 중 상기 적어도 하나의 변환 함수는, 상기 제1 볼륨 영상의 부분 영역 및 상기 제2 볼륨 영상의 부분 영역을 구 형상의 영역들로 평활화하는 적어도 하나의 파라미터에 기초하여 결정되는, 볼륨 파노라마 생성 장치.9. The method of claim 8,
Wherein the at least one transform function of the obtained transform functions is determined based on at least one parameter that smoothes the partial region of the first volume image and the partial region of the second volume image into spherical regions , Volume panorama generation device.
상기 복수의 변환 함수들로부터 상기 복수의 3차원 볼륨 영상들 사이의 유사도에 기초하여 최적화 변환 함수를 생성하는 단계;
상기 최적화 변환 함수를 이용하여 변환된 복수의 3차원 볼륨 영상들을 이용하여 볼륨 파노라마 영상을 생성하는 단계; 및
상기 볼륨 파노라마 영상을 출력하는 단계를 포함하고,
상기 유사도는, 상기 복수의 3차원 볼륨 영상들 중 제1 볼륨 영상의 형태적인 특성과, 상기 획득된 변환 함수들 중 하나에 기초하여 이동된 제2 볼륨 영상의 형태적인 특성 간의 유사도를 포함하고,
상기 볼륨 파노라마 생성 방법은, 상기 순차적인 볼륨 영상들의 대응하는 복셀들 각각의 위치 정보 및 강도(intensity)에 기초하여 상기 유사도를 결정하는 단계를 더 포함하는, 볼륨 파노라마 생성 방법.Dimensional volume image of a plurality of three-dimensional volume images, which is a function for matching at least two of the plurality of three-dimensional volume images and a plurality of three-dimensional volume images of a sequential order, from an external apparatus that generates three- Receiving a plurality of transform functions between the plurality of transform functions;
Generating an optimization transform function based on the similarity between the plurality of three-dimensional volume images from the plurality of transform functions;
Generating a volume panorama image using a plurality of transformed three-dimensional volume images using the optimization conversion function; And
And outputting the volume panorama image,
Wherein the degree of similarity includes a degree of similarity between morphological characteristics of a first volume image of the plurality of three-dimensional volume images and morphological characteristics of a second volume image moved based on one of the obtained conversion functions,
Wherein the volume panorama generation method further comprises determining the similarity based on location information and intensity of each of the corresponding voxels of the sequential volume images.
상기 외부 장치로부터 상기 복수의 3차원 볼륨 영상들과 함께 프로브의 위치 데이터를 수신하는 단계를 더 포함하는, 볼륨 파노라마 생성 방법.11. The method of claim 10,
And receiving position data of the probe along with the plurality of three-dimensional volume images from the external device.
상기 복수의 변환 함수들에 기초하여, 상기 순차적인 볼륨 영상들의 형태적인 특성의 유사도를 결정하는 단계; 및
상기 결정된 유사도에 기초하여 상기 최적화 변환 함수를 생성하는 단계를 더 포함하는 볼륨 파노라마 생성 방법.11. The method of claim 10,
Determining a similarity of the morphological characteristics of the sequential volume images based on the plurality of conversion functions; And
And generating the optimization transform function based on the determined similarity.
상기 유사도들의 합을 최대화하도록 상기 복수의 변환 함수들을 변경하고, 상기 변경된 변환 함수들에 기초하여 상기 최적화 변환 함수를 생성하는 단계를 더 포함하는, 볼륨 파노라마 생성 방법.11. The method of claim 10,
Modifying the plurality of transformation functions to maximize the sum of the similarities, and generating the optimization transform function based on the modified transformation functions.
상기 생성된 최적화 변환 함수에 기초하여 제 2 볼륨 영상의 영상 데이터로부터 제 1 볼륨 영상과 합성되는 볼륨 영상의 영상 데이터를 생성하는 단계; 및
상기 생성된 상기 제 1 볼륨 영상과 합성되는 볼륨 영상의 영상 데이터와 상기 제 1 볼륨 영상의 영상 데이터를 합성하는 단계를 더 포함하는, 볼륨 파노라마 생성 방법.11. The method of claim 10,
Generating image data of a volume image to be synthesized with a first volume image from image data of a second volume image based on the generated optimization conversion function; And
And synthesizing the image data of the volume image and the image data of the first volume image, which are combined with the generated first volume image.
상기 제 1 볼륨 영상과 합성된 볼륨 영상으로부터 분할된 지역 볼륨 영상들에 기초하여 지역 변환 함수를 결정하는 단계; 및
상기 결정된 지역 변환 함수에 기초하여 상기 생성된 상기 제 1 볼륨 영상과 합성되는 볼륨 영상의 영상 데이터를 갱신하는 단계를 더 포함하는, 볼륨 파노라마 생성 방법.16. The method of claim 15,
Determining a region transform function based on the local volume images segmented from the volume image combined with the first volume image; And
And updating the image data of the volume image to be synthesized with the generated first volume image based on the determined local conversion function.
상기 복수의 변환 함수들 중 적어도 하나의 변환 함수는 제1 볼륨 영상의 부분 영역과 제2 볼륨 영상의 부분 영역간의 변환 관계를 나타내는, 볼륨 파노라마 생성 방법.11. The method of claim 10,
Wherein at least one transformation function of the plurality of transformation functions represents a transformation relationship between a partial region of the first volume image and a partial region of the second volume image.
상기 획득된 변환 함수들 중 상기 적어도 하나의 변환 함수는, 상기 제1 볼륨 영상의 부분 영역 및 상기 제2 볼륨 영상의 부분 영역을 구 형상의 영역들로 평활화하는 적어도 하나의 파라미터에 기초하여 결정되는, 볼륨 파노라마 생성 방법.18. The method of claim 17,
Wherein the at least one transform function of the obtained transform functions is determined based on at least one parameter that smoothes the partial region of the first volume image and the partial region of the second volume image into spherical regions , Volume panorama generation method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| KR1020180078251A KR101946577B1 (en) | 2018-07-05 | 2018-07-05 | Method and apparatus for generating 3d volume panorama |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| KR1020180078251A KR101946577B1 (en) | 2018-07-05 | 2018-07-05 | Method and apparatus for generating 3d volume panorama |
Related Parent Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| KR1020110079153A Division KR101880634B1 (en) | 2011-08-09 | 2011-08-09 | Method and apparatus for generating 3d volume panorama |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| KR20180083282A KR20180083282A (en) | 2018-07-20 |
| KR101946577B1 true KR101946577B1 (en) | 2019-02-11 |
Family
ID=63103596
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| KR1020180078251A KR101946577B1 (en) | 2018-07-05 | 2018-07-05 | Method and apparatus for generating 3d volume panorama |
Country Status (1)
| Country | Link |
|---|---|
| KR (1) | KR101946577B1 (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20200109883A (en) * | 2019-03-15 | 2020-09-23 | 오스템임플란트 주식회사 | Method and Apparatus for generating x-ray image, computer-readable recording medium |
| CN116778127B (en) * | 2023-07-05 | 2024-01-05 | 广州视景医疗软件有限公司 | Panoramic view-based three-dimensional digital scene construction method and system |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060239571A1 (en) * | 2005-03-29 | 2006-10-26 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Method of volume-panorama imaging processing |
-
2018
- 2018-07-05 KR KR1020180078251A patent/KR101946577B1/en active IP Right Grant
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060239571A1 (en) * | 2005-03-29 | 2006-10-26 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Method of volume-panorama imaging processing |
Non-Patent Citations (1)
| Title |
|---|
| J Matas, et al., ‘Robust wide-baseline stereo form maximally stable extremal regions’, Image and Vision Computing, Vol. 22, Issue 10, 2004. 9.1. |
Also Published As
| Publication number | Publication date |
|---|---|
| KR20180083282A (en) | 2018-07-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR101805619B1 (en) | Apparatus and method for creating optimal 2-dimensional medical image automatically from 3-dimensional medical image | |
| RU2663649C2 (en) | Segmentation of large objects from multiple three-dimensional views | |
| EP3338625B1 (en) | Medical image display device and medical image processing method | |
| US9179890B2 (en) | Model-based positioning for intracardiac echocardiography volume stitching | |
| US8116548B2 (en) | Method and system for detecting 3D anatomical structures using constrained marginal space learning | |
| US6500123B1 (en) | Methods and systems for aligning views of image data | |
| KR101932721B1 (en) | Method and Appartus of maching medical images | |
| KR101805624B1 (en) | Method and apparatus for generating organ medel image | |
| US10026016B2 (en) | Tracking and representation of multi-dimensional organs | |
| KR101880634B1 (en) | Method and apparatus for generating 3d volume panorama | |
| US11013495B2 (en) | Method and apparatus for registering medical images | |
| US10945708B2 (en) | Method and apparatus for registration of medical images | |
| US10417777B2 (en) | Image processing apparatus, image processing method, and non-transitory computer-readable storage medium | |
| US20100123715A1 (en) | Method and system for navigating volumetric images | |
| KR101783000B1 (en) | Method and apparatus for generating 3d volume panorama based on a plurality of 3d volume images | |
| JP7267928B2 (en) | Volume rendered ultrasound image | |
| JP2000185036A (en) | Medical image display device | |
| EP3108456B1 (en) | Motion adaptive visualization in medical 4d imaging | |
| EP2601637B1 (en) | System and method for multi-modality segmentation of internal tissue with live feedback | |
| US20100274132A1 (en) | Arranging A Three-Dimensional Ultrasound Image In An Ultrasound System | |
| KR101946577B1 (en) | Method and apparatus for generating 3d volume panorama | |
| US9449425B2 (en) | Apparatus and method for generating medical image | |
| JP6501796B2 (en) | Acquisition Orientation Dependent Features for Model-Based Segmentation of Ultrasound Images | |
| KR100466409B1 (en) | System and method for displaying a virtual endoscopy and computer-readable recording medium having virtual endoscopy displaying program recorded thereon | |
| JP7313392B2 (en) | Optimal ultrasound organ segmentation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A107 | Divisional application of patent | ||
| A201 | Request for examination | ||
| E902 | Notification of reason for refusal | ||
| E701 | Decision to grant or registration of patent right | ||
| GRNT | Written decision to grant |