JPH0716240A - Within-image position transformation device - Google Patents
Within-image position transformation deviceInfo
- Publication number
- JPH0716240A JPH0716240A JP5190725A JP19072593A JPH0716240A JP H0716240 A JPH0716240 A JP H0716240A JP 5190725 A JP5190725 A JP 5190725A JP 19072593 A JP19072593 A JP 19072593A JP H0716240 A JPH0716240 A JP H0716240A
- Authority
- JP
- Japan
- Prior art keywords
- image
- coordinate system
- coordinates
- coordinate
- virtual
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000009466 transformation Effects 0.000 title abstract description 14
- 238000006243 chemical reaction Methods 0.000 claims description 57
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 abstract description 3
- 239000013598 vector Substances 0.000 description 25
- 239000003550 marker Substances 0.000 description 19
- 239000000523 sample Substances 0.000 description 14
- 238000001356 surgical procedure Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 8
- 238000000034 method Methods 0.000 description 8
- 210000004556 brain Anatomy 0.000 description 6
- 230000003902 lesion Effects 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 2
- 239000004698 Polyethylene Substances 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 238000007428 craniotomy Methods 0.000 description 1
- 238000012631 diagnostic technique Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- -1 polyethylene Polymers 0.000 description 1
- 229920000573 polyethylene Polymers 0.000 description 1
- 210000003625 skull Anatomy 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
Description
【0001】[0001]
【産業上の利用分野】この発明は、画像撮影装置で撮影
されたある画像が属する画像座標系の座標を、別の画像
座標系の座標に変換する画像内位置変換装置に係り、例
えば、X線CT装置や核磁気断層撮影装置(MRI装
置)等で得られた画像上に写し出された病変部位と実際
に手術を施す位置との対応関係をモニタ上で確認できる
ようにした定位的脳手術支援装置等に用いる画像内位置
変換装置に関する。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an intra-image position conversion device for converting coordinates of an image coordinate system to which an image photographed by an image photographing device belongs to coordinates of another image coordinate system. Stereotactic brain surgery that allows the operator to check the correspondence between the lesion site projected on the image obtained with a CT-ray machine or a nuclear magnetic tomography apparatus (MRI apparatus) and the actual surgery position on the monitor The present invention relates to an in-image position conversion device used for a support device or the like.
【0002】[0002]
【従来の技術】近年、X線CT装置やMRI装置等(以
下、これらを総称する場合には、単にCT装置ともい
う)の画像診断技術の進展に伴い、例えば、頭蓋内の病
変部位が高精度で決定できるようになった。しかしなが
ら、決定された病変部位に手術操作を加える段階では、
解剖学的知識に基づいてCT装置の画像と実際の手術位
置とのおおまかな対応関係を考えつつ、開頭位置の決
定、および病変部位の抽出等が行なわれる。このような
実際の手術操作位置とCT装置の画像上の病変部位との
正確な対応関係を付けるために、例えば、特開昭62-327
号や特開平3-267054号公報等に開示されている装置があ
る。2. Description of the Related Art In recent years, with the progress of image diagnostic techniques such as an X-ray CT apparatus and an MRI apparatus (hereinafter collectively referred to as a CT apparatus), for example, a lesion site in the skull has become high. You can now decide with precision. However, at the stage of applying a surgical operation to the determined lesion site,
The craniotomy position is determined and the lesion site is extracted while considering the rough correspondence between the CT image and the actual surgical position based on anatomical knowledge. In order to establish such an accurate correspondence between the actual surgical operation position and the lesion site on the CT image, for example, Japanese Patent Laid-Open No. 62-327
And Japanese Patent Application Laid-Open No. 3-267054.
【0003】この内、特開平3-267054号公報等に開示さ
れている定位的脳手術支援装置では、複数個の指標点
(マーカー)を頭部に貼付してCT装置で前記マーカー
を写し込んだ画像を撮影し、画像内に写し込まれたマー
カーのCT座標系(画像座標系)での座標と、磁場ソー
スと手術用プローブに取り付けられた磁場センサーと3
次元デジタイザー等からなる3次元座標測定装置で測定
した前記マーカーの実空間上(MS座標系)の位置(座
標)とに基づいて、MS座標系から画像座標系への変換
式を求め、手術用プローブの先端のMS座標上での位置
を3次元座標測定装置で測定し、その位置を前記変換式
を用いて画像座標系に変換し、プローブの先端位置をモ
ニタに表示されている画像に重畳表示して、手術操作位
置の特定等に役立てる。Among these, in the stereotactic brain surgery support device disclosed in Japanese Patent Application Laid-Open No. 3-267054, a plurality of index points (markers) are attached to the head and the marker is imaged by a CT device. Image of the marker, the coordinates of the marker embedded in the image in the CT coordinate system (image coordinate system), the magnetic field source, and the magnetic field sensor attached to the surgical probe.
Based on the position (coordinates) of the marker in the real space (MS coordinate system) measured by a three-dimensional coordinate measuring device such as a three-dimensional digitizer, a conversion formula from the MS coordinate system to the image coordinate system is obtained, and the conversion formula is used for surgery. The position of the tip of the probe on the MS coordinates is measured by a three-dimensional coordinate measuring device, the position is converted into an image coordinate system using the conversion formula, and the tip position of the probe is superimposed on the image displayed on the monitor. It is displayed to help identify the operation position of surgery.
【0004】[0004]
【発明が解決しようとする課題】しかしながら、このよ
うな構成を有する従来例の場合には、次のような問題が
ある。従来装置では、MS座標系から単一のCT装置で
単一の方向から得られた画像が属する画像座標系への変
換のみが可能であり、例えば、あるCT装置で撮影され
た画像上に重畳表示したプローブの先端の位置を、画像
座標系が異なる別のCT装置で撮影した画像上に重畳表
示することができないという問題がある。However, the conventional example having such a structure has the following problems. In the conventional apparatus, only the conversion from the MS coordinate system to the image coordinate system to which an image obtained from a single direction with a single CT apparatus belongs can be performed. For example, the image is superimposed on an image captured by a certain CT apparatus. There is a problem that the displayed tip position of the probe cannot be superimposed and displayed on an image taken by another CT apparatus having a different image coordinate system.
【0005】しかし、単一のCT装置により単一の方向
から得られた画像上に対するプローブの先端の表示のみ
では、手術中のオリエンテーションの認識が難しい場合
がある。However, it may be difficult to recognize orientation during surgery only by displaying the tip of the probe on an image obtained from a single direction by a single CT apparatus.
【0006】この発明は、このような事情に鑑みてなさ
れたものであって、ある画像が属する画像座標系の座標
を、別の画像座標系の座標に変換することができる画像
内位置変換装置を提供することを目的とする。The present invention has been made in view of the above circumstances, and is an image position conversion device capable of converting the coordinates of an image coordinate system to which an image belongs to the coordinates of another image coordinate system. The purpose is to provide.
【0007】[0007]
【課題を解決するための手段】この発明は、このような
目的を達成するために、次のような構成をとる。すなわ
ち、この発明は、画像撮影装置で撮影されたある画像が
属する画像座標系の座標を、別の画像座標系の座標に変
換する画像内位置変換装置であって、前記画像撮影装置
で、空間内に配置された複数個の指標点を写し込んで撮
影された異なる画像座標系に属する複数枚の画像を記憶
する画像記憶手段と、前記画像記憶手段に記憶された各
画像内に写し込まれた各指標点の、前記各画像が属する
画像座標系における座標を特定する画像座標特定手段
と、仮想座標系における前記各指標点の座標を特定する
仮想座標特定手段と、前記画像座標特定手段で特定され
た前記各画像座標系における前記各指標点の座標と、前
記仮想座標特定手段で特定された前記各指標点の座標と
に基づいて、前記各画像座標系と前記仮想座標系との各
変換式を特定する画像座標系/仮想座標系座標変換式特
定手段と、前記各画像座標系/仮想座標系座標変換式特
定手段で特定された各変換式に基づいて、各画像座標系
間の変換式を特定する画像座標系間座標変換式特定手段
とを備えたものである。The present invention has the following constitution in order to achieve such an object. That is, the present invention is an intra-image position conversion device that converts the coordinates of an image coordinate system to which an image captured by an image capturing device belongs to the coordinates of another image coordinate system, wherein Image storage means for storing a plurality of images belonging to different image coordinate systems captured by imprinting a plurality of index points arranged therein, and imprinted in each image stored in the image storage means The image coordinate specifying means for specifying the coordinates of each index point in the image coordinate system to which each image belongs, the virtual coordinate specifying means for specifying the coordinates of each index point in the virtual coordinate system, and the image coordinate specifying means. Based on the coordinates of each of the index points in each of the identified image coordinate systems and the coordinates of each of the index points identified by the virtual coordinate identifying means, each of the image coordinate system and the virtual coordinate system Specify conversion formula Based on each image coordinate system / virtual coordinate system coordinate conversion formula specifying unit and each conversion formula specified by each image coordinate system / virtual coordinate system coordinate conversion formula specifying unit. And an image coordinate system coordinate conversion formula specifying means.
【0008】[0008]
【作用】この発明の作用は次のとおりである。画像記憶
手段には、画像撮影装置で、空間内に配置された複数個
の指標点を写し込んで撮影された異なる画像座標系に属
する複数枚の画像が記憶されている。画像座標特定手段
は、画像記憶手段に記憶された各画像内に写し込まれた
各指標点の、前記各画像が属する画像座標系における座
標を特定する。仮想座標特定手段は、仮想座標系におけ
る前記各指標点の座標を特定する。画像座標系/仮想座
標系座標変換式特定手段は、画像座標特定手段で特定さ
れた各画像座標系における各指標点の座標と、仮想座標
特定手段で特定された各指標点の座標とに基づいて、各
画像座標系と仮想座標系との(各画像座標系ごとの)各
変換式を特定する。画像座標系間座標変換式特定手段
は、各画像座標系/仮想座標系座標変換式特定手段で特
定されたある画像座標系に対する変換式と別の画像座標
系に対する変換式とに基づいて、それら画像座標系間の
変換式を特定する。他の画像座標系間の変換式も同様に
特定する。すなわち、各画像座標系間の変換式を任意に
設定される仮想座標系を仲介として特定するものであ
る。The operation of the present invention is as follows. The image storage means stores a plurality of images belonging to different image coordinate systems, which are captured by the image capturing apparatus by capturing a plurality of index points arranged in the space. The image coordinate specifying means specifies the coordinates of each index point projected in each image stored in the image storage means in the image coordinate system to which each image belongs. The virtual coordinate specifying means specifies the coordinates of each index point in the virtual coordinate system. The image coordinate system / virtual coordinate system coordinate conversion formula specifying means is based on the coordinates of each index point in each image coordinate system specified by the image coordinate specifying means and the coordinates of each index point specified by the virtual coordinate specifying means. Then, each conversion formula (for each image coordinate system) between each image coordinate system and the virtual coordinate system is specified. The coordinate conversion formula specifying means between the image coordinate systems is based on the conversion formula for one image coordinate system specified by each image coordinate system / virtual coordinate system coordinate conversion specifying device and the conversion formula for another image coordinate system. Specify the conversion formula between image coordinate systems. Conversion formulas between other image coordinate systems are similarly specified. That is, the conversion formula between the image coordinate systems is specified by using a virtual coordinate system that is arbitrarily set as an intermediary.
【0009】[0009]
【実施例】以下、図面を参照してこの発明の一実施例を
説明する。図1は、この発明の一実施例に係る画像内位
置変換装置を組み込んだ定位的脳手術支援装置の構成を
示すブロック図である。DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS An embodiment of the present invention will be described below with reference to the drawings. FIG. 1 is a block diagram showing the configuration of a stereotactic brain surgery support device incorporating an intra-image position conversion device according to an embodiment of the present invention.
【0010】図中、符号1は各種の画像撮影装置であ
り、この画像撮影装置1は、例えば、X線CT装置やM
RI装置、ポジトロンCT装置等で構成されている。各
画像撮影装置1では、図2に示すように、被検者Mの任
意の位置に貼付された指標点(マーカー)M1 〜M4 を
写し込むように単数枚または複数枚の画像2が撮影され
る。例えば、画像撮影装置11 では、図3(a)に示す
ような方向で断層画像が撮影され、画像撮影装置12 で
は、図3(b)に示すような方向で断層画像が撮影さ
れ、画像撮影装置13 では、図3(c)に示すような方
向で断層画像が撮影される。撮影された画像を図4に示
す。図4(a)〜(c)に示すように、各画像21 、2
2 、23 にそれぞれマーカーM1 〜M4 が写し込まれる
ように画像が撮影される。撮影された各画像21 〜2n
は、例えば、固体撮影素子(CCD)カメラやイメージ
スキャナ等の画像読み取り装置(図示せず)で読み取ら
れて画像記憶手段としてのメモリ3に記憶される。ま
た、画像撮影装置1から直接画像データを転送してメモ
リ3に記憶するようにしてもよい。そして、メモリ3に
記憶された各画像2は、表示制御部4によりモニタ5に
表示される。なお、マーカーM1 〜M4 は、例えば、ポ
リエチレンチューブに造影剤を注入して構成されてい
る。In the figure, reference numeral 1 denotes various image capturing apparatuses, and the image capturing apparatus 1 is, for example, an X-ray CT apparatus or an M
It is composed of an RI apparatus, a positron CT apparatus, and the like. In each image capturing apparatus 1, as shown in FIG. 2, a single image or a plurality of images 2 are formed so that the index points (markers) M 1 to M 4 attached to arbitrary positions of the subject M are imaged. To be photographed. For example, the image capturing apparatus 1 1, the tomographic image is captured in a direction as shown in FIG. 3 (a), the image capturing apparatus 1 2 and the tomographic image in the direction as shown in FIG. 3 (b) is captured, in the image capturing apparatus 1 3, tomographic images are taken in the direction as shown in Figure 3 (c). The captured image is shown in FIG. As shown in FIGS. 4A to 4C, each of the images 2 1 , 2
Images are taken so that the markers M 1 to M 4 are imprinted on 2 and 2 3 , respectively. Each image taken 2 1 to 2 n
Is read by an image reading device (not shown) such as a solid-state imaging device (CCD) camera or an image scanner, and stored in the memory 3 as an image storage means. Alternatively, the image data may be directly transferred from the image capturing device 1 and stored in the memory 3. Then, each image 2 stored in the memory 3 is displayed on the monitor 5 by the display control unit 4. The markers M 1 to M 4 are formed by injecting a contrast agent into a polyethylene tube, for example.
【0011】図1の一点鎖線で示す部分は、この発明に
係る画像内位置変換装置であり、この画像内位置変換装
置10の各部の構成を、画像内位置変換手順に従って以
下に詳述する。The portion indicated by the alternate long and short dash line in FIG. 1 is the intra-image position conversion device according to the present invention, and the configuration of each part of this intra-image position conversion device 10 will be described in detail below in accordance with the intra-image position conversion procedure.
【0012】画像座標特定部11は、各画像21 〜2n
が属する画像座標系をそれぞれ特定し、それらの画像座
標系におけるマーカーM1 〜M4 の座標を特定する。ま
ず、各画像21 〜2n が属する画像座標系をそれぞれ特
定する。ここで、画像座標系とは、画像撮影装置11 〜
1n が持つ座標系を基準として特定してもよいし、モニ
タ5の座標系を基準にして特定してもよい。例えば、画
像21 について、断層面をX−Y平面とし、断層方向を
Z軸とし、特定の断層画像(例えば、図3(a)、図4
(a)に示す断層画像211)の中心点を原点とすること
により、画像21 が属する画像座標系を特定する。な
お、原点を含む断層画像の設定は入力装置6から行なわ
れる。また、画像22 および23 についても、断層面を
X−Y平面とし、断層方向をZ軸とし、図3(b)、図
4(b)および図3(c)図4(c)に示す断層画像2
21、231の中心点を原点とすることにより、画像22 、
23 が属する画像座標系を特定する。このように、各画
像21 〜2n が属する画像座標系をそれぞれ特定する。
なお、上記では各画像座標系を直交座標系として特定し
たが、斜交座標系で特定してもよい。The image coordinate specifying unit 11 is provided for each of the images 2 1 to 2 n.
Specify the image coordinate systems to which each belongs and specify the coordinates of the markers M 1 to M 4 in those image coordinate systems. First, the image coordinate system to which each of the images 2 1 to 2 n belongs is specified. Here, the image coordinate system refers to the image capturing devices 11 to.
The coordinate system of 1 n may be specified as a reference, or the coordinate system of the monitor 5 may be specified as a reference. For example, the image 2 1, the fault plane and the X-Y plane, the fault direction and the Z-axis, the specific tomographic image (e.g., FIG. 3 (a), FIG. 4
By using the center point of the tomographic image 2 11 ) shown in (a) as the origin, the image coordinate system to which the image 2 1 belongs is specified. The input device 6 sets the tomographic image including the origin. Further, the image 2 2 and 2 3 also, the fault plane and the X-Y plane, the fault direction and the Z-axis, in FIG. 3 (b), FIG. 4 (b) and FIG. 3 (c) FIG. 4 (c) Tomographic image 2
By setting the center points of 21 and 2 31 as the origins, the image 2 2 ,
Identify the image coordinate system to which 2 3 belongs. In this way, the image coordinate system to which each of the images 2 1 to 2 n belongs is specified.
Although each image coordinate system is specified as an orthogonal coordinate system in the above, it may be specified as an oblique coordinate system.
【0013】次に、メモリ3に記憶されている各画像2
1 〜2n を読み出し、各座標系におけるマーカーM1 〜
M4 の座標を特定する。例えば、画像21 に写されたマ
ーカーM3 については、図5に示すように、上述の画像
21 が属する画像座標系において、断層画像21 内で
X、Y座標を特定し、原点Oを含む断層画像から何枚目
の断層画像にマーカーM3 が写っているかと、断層厚さ
等の情報によりZ座標を特定する。このとき、原点Oを
含む断層画像から何枚目の断層画像にマーカーM3 が写
っているかの情報は、原点Oを含む断層画像(画像座標
系)が特定されているので、容易に検出することができ
る。また、断層厚さは、入力装置6から設定される。Next, each image 2 stored in the memory 3
1 to 2 n are read and the marker M 1 to in each coordinate system is
Specify the coordinates of M 4 . For example, the image for the 2 1 marker M 3 which photographed, as shown in FIG. 5, in the image 2 1 is an image coordinate system belonging above to identify X, Y coordinates in the tomographic image 2 within 1, the origin O The Z-coordinate is specified from the tomographic image including the marker M 3 in the tomographic image including the information and the information such as the tomographic thickness. At this time, the information about how many tomographic images from the tomographic image including the origin O the marker M 3 appears in is easy to detect because the tomographic image (image coordinate system) including the origin O is specified. be able to. Further, the slice thickness is set from the input device 6.
【0014】このようにして、各画像21 〜2n に属す
る画像座標系において、各画像21〜2n にそれぞれ写
っているマーカーM1 〜M4 の座標(xij,yij,
xij)(但し、iは各画像に対応し、i=1〜n、jは
各マーカーに対応し、j=1〜4)を特定する。このと
き、例えば、同じマーカーM1 であっても、画像21 の
属する画像座標系における座標(x11,y11,z11)
と、画像22 の属する画像座標系における座標(x21,
y21,z21)とは、それぞれの座標系が異なるので、一
致しない。特定された座標(xij,yij,xij)は、画
像座標系/仮想座標系座標変換式特定部12に与えられ
る。[0014] Thus, in the image coordinate system belonging to each image 2 1 to 2 n, of the marker M 1 ~M 4 that is reflected to each image 2 1 to 2 n coordinates (x ij, y ij,
x ij ) (where i corresponds to each image, i = 1 to n, j corresponds to each marker, and j = 1 to 4) are specified. At this time, for example, even with the same marker M 1 , the coordinates (x 11 , y 11 , z 11 ) in the image coordinate system to which the image 2 1 belongs
And the coordinates in the image coordinate system to which the image 2 2 belongs (x 21 ,
Since y 21 and z 21 ) have different coordinate systems, they do not match. The specified coordinates (x ij , y ij , x ij ) are given to the image coordinate system / virtual coordinate system coordinate conversion formula specifying unit 12.
【0015】仮想座標特定部13は、仮想的な座標系を
決め、その仮想座標系におけるマーカーM1 〜M4 の座
標を特定する。まず、仮想的な座標系を決める。ここ
で、仮想座標系は、実空間内において何を基準にして決
めてもよい。例えば、従来例(特開平3-267054号公報)
に開示されている、磁場ソース、手術用プローブに取り
付けられた磁場センサー、3次元デジタイザ等からなる
3次元座標測定装置7における磁場ソースの中心点を原
点としたMS座標系を仮想的な座標系として決めたとす
る。The virtual coordinate specifying unit 13 determines a virtual coordinate system and specifies the coordinates of the markers M 1 to M 4 in the virtual coordinate system. First, a virtual coordinate system is determined. Here, the virtual coordinate system may be determined based on anything in the real space. For example, a conventional example (Japanese Patent Laid-Open No. 3-267054)
Of the magnetic field source, the magnetic field sensor attached to the surgical probe, the three-dimensional coordinate measuring device 7 including the three-dimensional digitizer, etc. Suppose you decide.
【0016】このとき、その仮想座標系におけるマーカ
ーM1 〜M4 の座標は、固定された被検者に貼付された
マーカーM1 〜M4 をプローブの先端で指し、仮想的座
標系(MS座標系)における各マーカーM1 〜M4 の座
標(Xj ,Yj ,Zj )(但し、jは各マーカーに対応
し、j=1〜4)を、磁場センサーが受信した磁場ソー
スからの磁場の受信強度等に基づいて演算部8が特定す
る。特定された座標(Xj ,Yj ,Zj )は、画像座標
系/仮想座標系座標変換式特定部12に与えられる。At this time, the coordinates of the markers M 1 to M 4 in the virtual coordinate system refer to the markers M 1 to M 4 attached to the fixed subject by the tip of the probe, and the virtual coordinate system (MS Coordinates (X j , Y j , Z j ) (where j corresponds to each marker, j = 1 to 4) of each marker M 1 to M 4 in the coordinate system) from the magnetic field source received by the magnetic field sensor. The calculation unit 8 specifies the received intensity of the magnetic field. The specified coordinates (X j , Y j , Z j ) are given to the image coordinate system / virtual coordinate system coordinate conversion formula specifying unit 12.
【0017】画像座標系/仮想座標系座標変換式特定部
12は、画像座標特定部11で特定された各画像座標系
における各マーカーM1 〜M4 の座標(xij,yij,x
ij)と、仮想座標特定部13で特定された各マーカーM
1 〜M4 の座標(Xj ,Yj,Zj )とに基づいて、各
画像座標系と仮想座標系との各変換式を特定する。例え
ば、仮想座標系から画像21 が属する画像座標系への変
換を、各座標系で特定したマーカーM1 〜M4 の座標
(x1j,y1j,x1j)と(Xj ,Yj ,Zj)を用いて
4次元Affine(アフィン)変換で表すと、次式の
ようになる。The image coordinate system / virtual coordinate system coordinate conversion formula specifying unit 12 determines the coordinates (x ij , y ij , x) of the markers M 1 to M 4 in each image coordinate system specified by the image coordinate specifying unit 11.
ij ) and each marker M identified by the virtual coordinate identifying unit 13
Based on the coordinates (X j , Y j , Z j ) of 1 to M 4 , each conversion formula between each image coordinate system and the virtual coordinate system is specified. For example, the coordinates (x 1j , y 1j , x 1j ) and (X j , Y j of the markers M 1 to M 4 specified in each coordinate system for conversion from the virtual coordinate system to the image coordinate system to which the image 2 1 belongs are converted. , Z j ) is represented by a four-dimensional Affine transformation, the following equation is obtained.
【0018】[0018]
【数1】 [Equation 1]
【0019】ここで、T1 は、変換行列を表し、この変
換行列T1 は、上式(1)から次式のようになる。Here, T 1 represents a conversion matrix, and this conversion matrix T 1 is expressed by the following expression from the above expression (1).
【0020】[0020]
【数2】 [Equation 2]
【0021】各座標(x1j,y1j,x1j)、(Xj ,Y
j ,Zj )は、画像座標特定部11、仮想座標特定部1
3により特定されているので、変換行列T1 を求めるこ
とができる。Each coordinate (x 1j , y 1j , x 1j ), (X j , Y
j , Z j ) is the image coordinate specifying unit 11 and the virtual coordinate specifying unit 1
Since it is specified by 3, the conversion matrix T 1 can be obtained.
【0022】同様にして、仮想座標系から画像22 〜2
n が属する画像座標系への変換式を4次元Affine
変換で表し、各変換行列T2 〜Tn を求める。求めた各
変換行列T1 〜Tn は、画像座標系間座標変換式特定部
14に与えられる。Similarly, from the virtual coordinate system, the images 2 2 to 2 2
A four-dimensional Affine conversion formula to the image coordinate system to which n belongs
The transformation matrix is used to obtain each transformation matrix T 2 to T n . The obtained transformation matrices T 1 to T n are provided to the inter-image coordinate system coordinate transformation formula specifying unit 14.
【0023】画像座標系間座標変換式特定部14は、各
画像座標系/仮想座標系座標変換式特定部12で特定さ
れた各変換式に基づいて、各画像座標系間の変換式を特
定する。例えば、仮想座標系での任意の点の座標(X,
Y,Z)を、画像2k (k=1〜n)が属する画像座標
系の座標(xk ,yk ,zk )に変換する場合には、上
式(1)から次式のように表すことができる。 (xk ,yk ,zk ,1)=(X,Y,Z,1)*Tk …… (3)The inter-image coordinate system coordinate conversion formula specifying unit 14 specifies a conversion formula between the image coordinate systems based on each conversion formula specified by each image coordinate system / virtual coordinate system coordinate conversion formula specifying unit 12. To do. For example, the coordinates (X,
When converting Y, Z) to the coordinates (x k , y k , z k ) of the image coordinate system to which the image 2 k (k = 1 to n) belongs, the following equation is obtained from the above equation (1). Can be expressed as (X k , y k , z k , 1) = (X, Y, Z, 1) * T k (3)
【0024】同様に、仮想座標系での任意の点の座標
(X,Y,Z)を、画像2l (l=1〜n)が属する画
像座標系の座標(xl ,yl ,zl )に変換する場合に
は次式で表すことができる。 (xl ,yl ,zl ,1)=(X,Y,Z,1)*Tl …… (4)Similarly, the coordinates (X, Y, Z) of an arbitrary point in the virtual coordinate system are converted into the coordinates (x l , y l , z) of the image coordinate system to which the image 2 l (l = 1 to n) belongs. When converted into l ), it can be expressed by the following equation. (X l , y l , z l , 1) = (X, Y, Z, 1) * T l (4)
【0025】ここで、(X,Y,Z,1)は、(4)式
から次式で表すことができる。 (X,Y,Z,1)=(xl ,yl ,zl ,1)*Tl -1 ……(5) (3)式の(X,Y,Z,1)に(5)式を代入してや
ると、次式となる。 (xk ,yk ,zk ,1)=((xl ,yl ,zl ,1)*Tl -1)*Tk =(xl ,yl ,zl ,1)*(Tl -1*Tk ) ……(6) ここで、Tk 、Tl は、既に求められており、T
l -1は、Tl から求めることができるので、画像2l が
属する画像座標系での任意の点の座標(xl ,yl ,z
l )を、画像2k が属する画像座標系の座標(xk ,y
k ,zk )に変換することができる。Here, (X, Y, Z, 1) can be expressed by the following equation from the equation (4). (X, Y, Z, 1) = (x l , y l , z l , 1) * T l -1 (5) In (X, Y, Z, 1) of the equation (3), (5) Substituting an expression gives the following expression. (X k , y k , z k , 1) = ((x l , y l , z l , 1) * T l −1 ) * T k = (x l , y l , z l , 1) * ( T l -1 * T k ) (6) Here, T k and T l have already been calculated, and T
Since l −1 can be obtained from T l , the coordinates (x l , y l , z) of an arbitrary point in the image coordinate system to which the image 2 l belongs
l ) is the coordinates (x k , y) of the image coordinate system to which the image 2 k belongs.
can be converted to k , z k ).
【0026】例えば、3次元座標測定装置7によりMS
座標系におけるプローブの先端位置の座標が特定される
と、演算部8により、その座標が、画像21 が属する画
像座標系の座標(x1j,y1j,z1j)に変換され、表示
制御部4により、モニタ5に表示されている画像21 に
プローブの先端位置が重畳表示される。このとき、従来
例では、他の画像22 〜2n には、プローブの先端位置
は表示されない。しかし、この実施例では、画像21 が
属する画像座標系におけるプローブの座標(x1j,
y1j,z1j)に基づいて、他の画像22 〜2n が属する
画像座標系におけるプローブ先端位置の座標(xij,y
ij,zij)を求めることができので、図6に示すよう
に、各画像21 〜2n (図では、21 〜23 を示してい
る)上にそれぞれプローブの先端位置S(「×」で示
す)を重畳表示することができる。For example, by using the three-dimensional coordinate measuring device 7, the MS
When the coordinates of the tip position of the probe in the coordinate system are specified, the calculation unit 8 converts the coordinates into the coordinates (x 1j , y 1j , z 1j ) of the image coordinate system to which the image 2 1 belongs, and the display control is performed. the part 4, the tip position of the probe in image 2 1 displayed on the monitor 5 is superimposed and displayed. In this case, in the conventional example, the other image 2 2 to 2 n, the tip position of the probe is not displayed. However, in this embodiment, the probe in the image coordinate system of the image 2 1 belongs coordinates (x 1j,
y 1j , z 1j ) based on the coordinates (x ij , y) of the probe tip position in the image coordinate system to which the other images 2 2 to 2 n belong.
ij, since it is possible to determine the z ij), as shown in FIG. 6, in each image 2 1 to 2 n (figure 2 1 to 2 probes each 3 on the show) tip position S ( " (Indicated by “X”) can be displayed in a superimposed manner.
【0027】ところで、上述では、仮想的な座標系とし
てMS座標系を用いたが、他の座標系を仮想的な座標系
としてもよい。以下に、仮想的な座標系を特定する方法
をいくつか説明する。まず、各画像21 〜2n に写され
た4個のマーカーM1 〜M4 を用いて、仮想的な座標系
を作成する場合を説明する。同一平面上にないマーカー
M1 〜M4 の画像座標系での座標(xij,yij,zij)
を用い、マーカーM1 を原点とし、M1 からM2 、M1
からM3 、M1 からM4 への方向性を有する線形独立な
ベクトルi、j、kを作成し、各ベクトルで張られる座
標空間を仮想的な座標系とする。By the way, although the MS coordinate system is used as the virtual coordinate system in the above description, another coordinate system may be used as the virtual coordinate system. Hereinafter, some methods for specifying the virtual coordinate system will be described. First, a case will be described in which a virtual coordinate system is created using the four markers M 1 to M 4 imaged on each of the images 2 1 to 2 n . Coordinates (x ij , y ij , z ij ) in the image coordinate system of the markers M 1 to M 4 that are not on the same plane
Using the marker M 1 as the origin, M 1 to M 2 , M 1
To M 3 and M 1 to M 4 , linearly independent vectors i, j, and k having directions are created, and the coordinate space spanned by each vector is used as a virtual coordinate system.
【0028】具体的には、例えば、画像21 に写された
マーカーM1 〜M4 の画像21 が属する画像座標系での
座標(x1j,y1j,z1j)を用いて各ベクトルi、j、
kを表すと、次式のようになる。[0028] Specifically, for example, coordinates in the image coordinate system which the image 2 1 markers M 1 ~M 4 was photographed in the image 2 1 belongs (x 1j, y 1j, z 1j) using the vectors i, j,
When k is expressed, the following equation is obtained.
【0029】[0029]
【数3】 [Equation 3]
【0030】これらのベクトルi、j、kを仮想的な座
標系における単位ベクトルとすると、 原点(x11,y11,z11)=(0,0,0) i=(aX ,aY ,aZ )=(1,0,0) j=(bX ,bY ,bZ )=(0,1,0) k=(cX ,cY ,cZ )=(0,0,1) となり、画像座標系での各マーカーM1 〜M4 の位置
は、 M1 =(x11,y11,z11)=(0,0,0) M2 =(x12,y12,z12)=(1,0,0) M3 =(x13,y13,z13)=(0,1,0) M4 =(x14,y14,z14)=(0,0,1) となる。従って、仮想的な座標系から画像21 が属する
画像座標系への変換を次式で表すことができる。Assuming that these vectors i, j, k are unit vectors in a virtual coordinate system, the origin (x 11 , y 11 , z 11 ) = (0,0,0) i = (a X , a Y , A Z ) = (1,0,0) j = (b X , b Y , b Z ) = (0,1,0) k = (c X , c Y , c Z ) = (0,0, 1), and the positions of the markers M 1 to M 4 in the image coordinate system are M 1 = (x 11 , y 11 , z 11 ) = (0,0,0) M 2 = (x 12 , y 12 , Z 12 ) = (1,0,0) M 3 = (x 13 , y 13 , z 13 ) = (0,1,0) M 4 = (x 14 , y 14 , z 14 ) = (0, 0, 1). Therefore, the conversion from the virtual coordinate system to the image coordinate system to which the image 2 1 belongs can be expressed by the following equation.
【0031】[0031]
【数4】 [Equation 4]
【0032】また、画像22 に写されたマーカーM1 〜
M4 の画像22 が属する画像座標系での座標(x2j,y
2j,z2j)を用いて作成した各ベクトルi、j、kで張
られる座用空間を仮想座標系として、仮想的な座標系か
ら画像21 が属する画像座標系への変換式を表すと、次
式のようになる。[0032] In addition, image 2 2 photographed the markers M 1 ~
The coordinates (x 2j , y in the image coordinate system to which the image 2 2 of M 4 belongs
2j , z 2j ), a transformation space from the virtual coordinate system to the image coordinate system to which the image 2 1 belongs is represented by using the sitting space stretched by each vector i, j, k created as a virtual coordinate system. , Is as follows.
【0033】[0033]
【数5】 [Equation 5]
【0034】画像21 が属する画像座標系と、画像22
が属する画像座標系とは一致しないが、各画像座標系で
のマーカーM1 〜M4 の座標を用いて作成した各ベクト
ルi、j、kは同じ方向性を示しているので、それらの
ベクトルi、j、kで張られる座標空間(仮想座標系)
は同じ座標系を示すことになる。The image coordinate system to which the image 2 1 belongs and the image 2 2
The vector i, j, k created by using the coordinates of the markers M 1 to M 4 in each image coordinate system does not match the image coordinate system to which the vector belongs. A coordinate space defined by i, j, and k (virtual coordinate system)
Indicate the same coordinate system.
【0035】次に、画像21 〜2n に写された3点のマ
ーカーM1 〜M3 を用いて、仮想的な座標系を作成する
場合を説明する。同一直線上にないマーカーM1 〜M3
の画像座標系での座標(xik,yik,zik)(但し、k
はマーカーに対応し、k=1〜3)を用い、マーカーM
1 を原点とし、以下のように線形独立なベクトルi、
j、kを作成し、各ベクトルで張られる座標空間を仮想
的な座標系とする。Next, a case will be described in which a virtual coordinate system is created using the three markers M 1 to M 3 imaged on the images 2 1 to 2 n . Markers are not on the same straight line M 1 ~M 3
In the image coordinate system of (x ik , y ik , z ik ) (where k is
Corresponds to a marker, and k = 1 to 3) is used, and the marker M
With 1 as the origin, linearly independent vector i,
j and k are created, and the coordinate space defined by each vector is used as a virtual coordinate system.
【0036】まず、M1 からM2 、M1 からM3 への方
向性を有するベクトルl、sと、ベクトルlとベクトル
sの外積で求められるベクトルn、ベクトルnとベクト
ルlの外積で求められるベクトルmを作成する。これを
次式に示す。First, the vectors l and s having directivity from M 1 to M 2 and M 1 to M 3 and the vector n obtained by the outer product of the vector l and the vector s, and the outer product of the vector n and the vector 1 are obtained. Create a vector m that is This is shown in the following equation.
【0037】[0037]
【数6】 [Equation 6]
【0038】次に、ベクトルl、m、nをそれぞれの長
さで割ったベクトルをi、j、kとする。これを以下に
示す。 i=l/|l| j=m/|m| k=n/|n| 原点をM1 とし、線形独立なベクトルi、j、kで張ら
れる直交座標空間を仮想的な座標系とする。なお、上述
のi、j、kを作成する手順を図7に示す。図7(a)
は、ベクトルl、sを示し、図7(b)は、ベクトル
l、sから作成されたベクトルnを示し、図7(c)
は、ベクトルn、lから作成されたベクトルmを示し、
図7(d)は、ベクトルi、j、kを示す。Next, the vectors l, m, and n are divided by their respective lengths to be i, j, and k. This is shown below. i = l / | l | j = m / | m | k = n / | n | The origin is M 1 and the orthogonal coordinate space spanned by linearly independent vectors i, j, and k is a virtual coordinate system. . The procedure for creating the above i, j, and k is shown in FIG. Figure 7 (a)
Shows the vectors l and s, FIG. 7B shows the vector n created from the vectors l and s, and FIG.
Indicates a vector m created from the vectors n and l,
FIG. 7D shows the vectors i, j, and k.
【0039】なお、仮想的な座標系は、各画像座標系間
の変換式を特定するための仲介としての機能、すなわ
ち、各画像座標系からみて同じ座標系であれば、上述以
外の方法で特定してもよい。The virtual coordinate system has a function as an intermediary for specifying a conversion formula between the image coordinate systems, that is, if the coordinate system is the same as seen from the image coordinate systems, a method other than the above is used. May be specified.
【0040】また、仮想的な座標系から画像座標系への
変換は、上述の方法以外であってもよい。Further, the conversion from the virtual coordinate system to the image coordinate system may be performed by a method other than the above method.
【0041】さらに、同じ断層撮影装置で撮影された画
像であっても、撮影された方向が違うこと等により、そ
れぞれの画像が属する画像座標系が違う場合の各画像座
標系間の変換も行なうことができる。Further, even in the case of images taken by the same tomographic imaging apparatus, conversion between the image coordinate systems is performed when the image coordinate systems to which the images belong are different due to the different shooting directions. be able to.
【0042】また、上述の実施例では、定位的脳手術支
援装置に組み込んだ場合について説明したが、その他の
画像座標系間の座標の変換を行なう場合に適用すること
もできる。Further, in the above-mentioned embodiment, the case where the device is incorporated in the stereotactic brain surgery support device has been described, but the present invention can be applied to the case of converting the coordinates between other image coordinate systems.
【0043】[0043]
【発明の効果】以上の説明から明らかなように、この発
明によれば、ある画像の属する画像座標系と別の画像が
属する画像座標系とを仮想的な座標系を介して変換可能
に構成したことにより、この画像内位置変換装置を、例
えば、定位的脳手術支援装置に組み込めば、ある画像座
標系の画像上に表示された任意点に対応する対応点を他
の画像座標系の画像上に重畳表示することできるので、
手術中のオリエンテーションの認識を容易にすることが
できる。As is apparent from the above description, according to the present invention, an image coordinate system to which an image belongs and an image coordinate system to which another image belongs can be converted via a virtual coordinate system. Therefore, if this intra-image position conversion device is incorporated into, for example, a stereotactic brain surgery support device, the corresponding point corresponding to an arbitrary point displayed on the image of a certain image coordinate system becomes an image of another image coordinate system. Since it can be superimposed on top,
It can facilitate the recognition of orientation during surgery.
【図1】この発明の一実施例に係る画像内位置変換装置
を組み込んだ定位的脳手術支援装置の構成を示すブロッ
ク図である。FIG. 1 is a block diagram showing the configuration of a stereotactic brain surgery support device incorporating an intra-image position conversion device according to an embodiment of the present invention.
【図2】被検者に貼付されたマーカーの位置を示す図で
ある。FIG. 2 is a diagram showing positions of markers attached to a subject.
【図3】各画像撮影装置で撮影される画像の断層方向を
示す図である。FIG. 3 is a diagram showing a tomographic direction of an image captured by each image capturing device.
【図4】各画像撮影装置で撮影された画像を示す図であ
る。FIG. 4 is a diagram showing an image captured by each image capturing device.
【図5】画像座標系におけるマーカーの座標を特定する
手順を説明する図である。FIG. 5 is a diagram illustrating a procedure of identifying coordinates of a marker in an image coordinate system.
【図6】各画像にプローブの先端位置を重畳表示した状
態を示す図である。FIG. 6 is a diagram showing a state in which the tip position of the probe is displayed in an overlapping manner on each image.
【図7】3個のマーカーから仮想座標系を特定する手順
を示す図である。FIG. 7 is a diagram showing a procedure for identifying a virtual coordinate system from three markers.
11 〜1n … 画像撮影装置 21 〜2n … 画像 3 … メモリ(画像記憶手段) 11 … 画像座標特定部(画像座標特定手段) 12 … 仮想座標特定部(仮想座標特定手段) 13 … 画像座標系/仮想座標系座標変換式特定部 (画像座標系/仮想座標系座標変換式特定手段) 14 … 画像座標系間座標変換式特定部 (画像座標系間座標変換式特定手段) M1 〜M4 … マーカー(指標点)1 1 to 1 n ... Image photographing device 2 1 to 2 n ... Image 3 ... Memory (image storage means) 11 ... Image coordinate identification unit (image coordinate identification means) 12 ... Virtual coordinate identification unit (virtual coordinate identification means) 13 ... Image coordinate system / virtual coordinate system coordinate conversion formula specifying unit (image coordinate system / virtual coordinate system coordinate conversion formula specifying unit) 14 ... Image coordinate system coordinate conversion formula specifying unit (image coordinate system coordinate conversion formula specifying unit) M 1 ~ M 4 ... Marker (index point)
───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.6 識別記号 庁内整理番号 FI 技術表示箇所 G06T 1/00 3/00 8420−5L G06F 15/66 360 ─────────────────────────────────────────────────── ─── Continuation of the front page (51) Int.Cl. 6 Identification code Internal reference number FI Technical display location G06T 1/00 3/00 8420-5L G06F 15/66 360
Claims (1)
する画像座標系の座標を、別の画像座標系の座標に変換
する画像内位置変換装置であって、前記画像撮影装置
で、空間内に配置された複数個の指標点を写し込んで撮
影された異なる画像座標系に属する複数枚の画像を記憶
する画像記憶手段と、前記画像記憶手段に記憶された各
画像内に写し込まれた各指標点の、前記各画像が属する
画像座標系における座標を特定する画像座標特定手段
と、仮想座標系における前記各指標点の座標を特定する
仮想座標特定手段と、前記画像座標特定手段で特定され
た前記各画像座標系における前記各指標点の座標と、前
記仮想座標特定手段で特定された前記各指標点の座標と
に基づいて、前記各画像座標系と前記仮想座標系との各
変換式を特定する画像座標系/仮想座標系座標変換式特
定手段と、前記各画像座標系/仮想座標系座標変換式特
定手段で特定された各変換式に基づいて、各画像座標系
間の変換式を特定する画像座標系間座標変換式特定手段
とを備えたことを特徴とする画像内位置変換装置。1. An in-image position conversion device for converting the coordinates of an image coordinate system to which an image captured by an image capturing device belongs to the coordinates of another image coordinate system, wherein Image storage means for storing a plurality of images belonging to different image coordinate systems, which are photographed by imprinting a plurality of index points arranged in the image storage means, and imprinted in each image stored in the image storage means. Image coordinate specifying means for specifying coordinates of each index point in the image coordinate system to which each image belongs, virtual coordinate specifying means for specifying coordinates of each index point in the virtual coordinate system, and specification by the image coordinate specifying means Based on the coordinates of each of the index points in each of the image coordinate systems and the coordinates of each of the index points identified by the virtual coordinate identifying means, each conversion between each of the image coordinate systems and the virtual coordinate system. Image Zodiac that identifies an expression An image that specifies the conversion formula between the image coordinate systems based on the standard system / virtual coordinate system coordinate conversion formula specifying unit and each conversion formula specified by the image coordinate system / virtual coordinate system coordinate conversion formula specifying unit An intra-image position conversion device, comprising: an inter-coordinate system coordinate conversion expression specifying means.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP5190725A JPH0716240A (en) | 1993-07-01 | 1993-07-01 | Within-image position transformation device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP5190725A JPH0716240A (en) | 1993-07-01 | 1993-07-01 | Within-image position transformation device |
Publications (1)
Publication Number | Publication Date |
---|---|
JPH0716240A true JPH0716240A (en) | 1995-01-20 |
Family
ID=16262775
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP5190725A Pending JPH0716240A (en) | 1993-07-01 | 1993-07-01 | Within-image position transformation device |
Country Status (1)
Country | Link |
---|---|
JP (1) | JPH0716240A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001224586A (en) * | 2000-02-11 | 2001-08-21 | Toshiba Corp | Image mapping method and system |
WO2007102509A1 (en) * | 2006-03-09 | 2007-09-13 | Imagnosis Inc. | Medical image processing method |
JP2008080131A (en) * | 2006-09-28 | 2008-04-10 | Accuray Inc | Radiation treatment plan using four-dimensional imaging data |
-
1993
- 1993-07-01 JP JP5190725A patent/JPH0716240A/en active Pending
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001224586A (en) * | 2000-02-11 | 2001-08-21 | Toshiba Corp | Image mapping method and system |
JP2011078822A (en) * | 2000-02-11 | 2011-04-21 | Toshiba Corp | Image mapping method and system |
WO2007102509A1 (en) * | 2006-03-09 | 2007-09-13 | Imagnosis Inc. | Medical image processing method |
JPWO2007102509A1 (en) * | 2006-03-09 | 2009-07-23 | イマグノーシス株式会社 | Medical image processing method |
JP2008080131A (en) * | 2006-09-28 | 2008-04-10 | Accuray Inc | Radiation treatment plan using four-dimensional imaging data |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4822634B2 (en) | A method for obtaining coordinate transformation for guidance of an object | |
US7634122B2 (en) | Registering intraoperative scans | |
JP2950340B2 (en) | Registration system and registration method for three-dimensional data set | |
US6351661B1 (en) | Optically coupled frameless stereotactic space probe | |
US6259943B1 (en) | Frameless to frame-based registration system | |
US6049622A (en) | Graphic navigational guides for accurate image orientation and navigation | |
US7760909B2 (en) | Video tracking and registering | |
JP5027922B2 (en) | Ultrasonic diagnostic equipment | |
EP1467317B1 (en) | Combining first and second image data of an object | |
KR19990023201A (en) | Image display method and image display device | |
JP2002330958A (en) | Method and device for selecting and displaying medical image data | |
US5509084A (en) | Method and apparatus for postprocessing medical images | |
US20090080742A1 (en) | Image display device and image display program storage medium | |
DE102011114333B4 (en) | A method of registering an X-ray volume with a position sensing system using a registration phantom | |
US10078906B2 (en) | Device and method for image registration, and non-transitory recording medium | |
JP2013220326A (en) | Transcranial brain function analyzing method | |
US7965882B2 (en) | Image display apparatus and computer-readable image display program storage medium | |
US20100063420A1 (en) | Method for verifying the relative position of bone structures | |
WO2001057805A2 (en) | Image data processing method and apparatus | |
JP2001017422A (en) | Image processing device and marker member for the same | |
JPH0716240A (en) | Within-image position transformation device | |
JP3634416B2 (en) | Surgical instrument position display device | |
JPH03231643A (en) | Image display device | |
JPH087080A (en) | Image diagnostic device | |
JP2018061844A (en) | Information processing apparatus, information processing method, and program |