JPH0760088B2 - Target area volume measuring device - Google Patents
Target area volume measuring deviceInfo
- Publication number
- JPH0760088B2 JPH0760088B2 JP63185441A JP18544188A JPH0760088B2 JP H0760088 B2 JPH0760088 B2 JP H0760088B2 JP 63185441 A JP63185441 A JP 63185441A JP 18544188 A JP18544188 A JP 18544188A JP H0760088 B2 JPH0760088 B2 JP H0760088B2
- Authority
- JP
- Japan
- Prior art keywords
- image information
- contour
- distance value
- slice
- volume
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 claims description 14
- 238000012545 processing Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 206010028980 Neoplasm Diseases 0.000 description 2
- 210000000988 bone and bone Anatomy 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000002512 chemotherapy Methods 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000001959 radiotherapy Methods 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 230000001225 therapeutic effect Effects 0.000 description 1
Landscapes
- Measurement Of Length, Angles, Or The Like Using Electric Or Magnetic Means (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Length-Measuring Devices Using Wave Or Particle Radiation (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Description
【発明の詳細な説明】 ［発明の目的］ （産業上の利用分野） 本発明は被検体内骨、腫瘍等の対象物領域の体積を正確
に計測する被検体対象領域体積計測装置に関する。DETAILED DESCRIPTION OF THE INVENTION (Industrial field of application) The present invention relates to an object target area volume measuring device for accurately measuring the volume of an object area such as bone or tumor in an object.
（従来の技術） 従来、骨、腫瘍等の対象物に対する手術、放射線療法、
化学療法等の治療の効果を評価するため、Ｘ線CT（Comp
uted Tomography）あるいは磁気共鳴イメージング装置
を用い治療前後における対象物の画像を撮影し、対象物
の体積変化を調べることが行われている。(Conventional technology) Conventionally, surgery, radiation therapy for an object such as bone or tumor,
In order to evaluate the effect of treatment such as chemotherapy, X-ray CT (Comp
uted tomography) or a magnetic resonance imaging apparatus, images of an object before and after treatment are taken to examine the volume change of the object.
（発明が解決しようとする課題） しかしながらこのような従来技術には次に述べるような
問題がある。(Problems to be Solved by the Invention) However, such a conventional technique has the following problems.
すなわち、CTあるいはMRIにより得られるものは対象物
の断面に関する情報であるにすぎず、対象物の体積に関
する情報を直接与えるものではない。That is, what is obtained by CT or MRI is only information about the cross section of the object, and does not directly give information about the volume of the object.
本発明は治療効果を正確に評価することができるよう
に、対象物の体積を直接正確に計測する被検体対象領域
体積計測装置を提供することを目的とするものである。It is an object of the present invention to provide an object target area volume measuring device that directly and accurately measures the volume of an object so that the therapeutic effect can be accurately evaluated.
［発明の構成］ （課題を解決するための手段） 上記目的を達成するため、本発明の被検体対象領域体積
計測装置においては、 被検体に関するマルチスライス画像情報を入力し、隣
接するスライス画像情報間に挿入すべきスライス画像情
報を補間演算により作成する手段と、各スライス画像情
報および補間画像情報における被検体の対象領域の輪郭
を認識する輪郭認識手段と、輪郭認識手段の認識結果に
基づき被検体の対象領域の体積の演算を行う体積演算手
段とを有する。[Structure of the Invention] (Means for Solving the Problems) In order to achieve the above-mentioned object, in a subject object region volume measuring apparatus of the present invention, multi-slice image information relating to a subject is input, and adjacent slice image information is input. A means for creating slice image information to be inserted between them by interpolation calculation, a contour recognition means for recognizing the contour of the target area of the subject in each slice image information and the interpolated image information, and a contour recognition means based on the recognition result of the contour recognition means. And a volume calculation means for calculating the volume of the target region of the sample.
被検体に関するマルチスライス画像情報を入力し、各
スライス画像情報における被検体の対象領域の輪郭を認
識する輪郭認識手段と、前記各スライス画像の各画素
と、前記輪郭認識手段により認識された対象領域の輪郭
部位との距離を求め、その距離に応じた画素値を有する
距離値画像情報を作成する距離値画像情報作成手段と、
距離値画像情報を補間して補間距離値画像情報を作成す
る補間距離値画像情報作成手段と、輪郭認識手段におけ
る認識結果と補間距離値画像情報とに基づき被検体の対
象領域の体積の演算を行う体積演算手段とを有する。Inputting multi-slice image information about a subject, a contour recognizing unit that recognizes a contour of a target region of the subject in each slice image information, each pixel of each slice image, and a target region recognized by the contour recognizing unit A distance value image information creating means for creating a distance value image information having a pixel value according to the distance,
Interpolation distance value image information creating means for creating the interpolation distance value image information by interpolating the distance value image information, and calculation of the volume of the target region of the subject based on the recognition result in the contour recognition means and the interpolation distance value image information. And a volume calculating means for performing the calculation.
（作用） 上記のように構成された被検体対象領域体積計測装置
において、補間画像情報作成手段は最終的なスライス間
隔が各スライス画像情報の空間分解能に等しくなるよう
に補間画像情報の作成を行い、輪郭認識手段は各スライ
ス画像情報および補間画像情報において被検体の対象領
域の輪郭を認識し、体積演算手段は輪郭認識手段の認識
結果にもとづき被検体の対象領域の体積を演算する。(Operation) In the subject object area volume measuring device configured as described above, the interpolation image information creating unit creates the interpolation image information so that the final slice interval becomes equal to the spatial resolution of each slice image information. The contour recognition means recognizes the contour of the target area of the subject in each slice image information and the interpolated image information, and the volume calculation means calculates the volume of the target area of the subject based on the recognition result of the contour recognition means.
あるいは、上記のように構成された被検体対象領域体
積計測装置において、輪郭認識手段は各スライス画像情
報において被検体の対象領域の輪郭を認識し、距離値画
像情報作成手段は輪郭認識手段の認識結果にもとづき距
離値画像情報を作成し、補間距離値画像情報作成手段は
距離値画像情報を補間する補間距離値画像情報を作成
し、体積演算手段は輪郭認識手段における認識結果と補
間距離値画像情報とにもとずき被検体の対象領域の体積
を演算する。Alternatively, in the object target area volume measuring device configured as described above, the contour recognition means recognizes the contour of the target area of the object in each slice image information, and the distance value image information creation means recognizes the contour recognition means. Distance value image information is created based on the result, interpolation distance value image information creating means creates interpolated distance value image information that interpolates distance value image information, and volume calculation means is the recognition result in the contour recognition means and the interpolated distance value image. The volume of the target region of the subject is calculated based on the information.
（実施例） 以下図面を参照しながら本発明の実施例について説明を
行う。Embodiments Embodiments of the present invention will be described below with reference to the drawings.
第１図に本発明の第１の実施例の対象領域体積計測装置
のブロック構成を示す。第２図は同装置の動作を示すフ
ローチャートである。FIG. 1 shows a block configuration of a target area volume measuring apparatus according to a first embodiment of the present invention. FIG. 2 is a flowchart showing the operation of the device.
本実施例にあってはスライス画像情報から後述の距離値
画像情報を作成するものであり、第３図は同装置におけ
る距離値画像情報および補間距離値画像情報を示す説明
図である。In the present embodiment, distance value image information described later is created from slice image information, and FIG. 3 is an explanatory diagram showing distance value image information and interpolated distance value image information in the same apparatus.
第４図は第３図において太線で示される各画像上の対象
領域輪郭情報を３次元的に表現したものである。FIG. 4 is a three-dimensional representation of the target area contour information on each image indicated by the bold line in FIG.
第１図に示す被検体対象領域体積計測装置１は、Ｘ線CT
装置や磁気共鳴診断装置（MRI）等から出力される被検
体のマルチスライス画像（連続断層像）情報を取込ん
で、従来公知の方法により各スライス画像情報における
被検体の対象領域の輪郭の抽出、輪郭の境界認識を行う
輪郭認識部２と、この輪郭認識部２により認識された被
検体の対象領域の輪郭情報をもとに、前記各スライス画
像情報を距離値画像情報に変換する距離値画像情報作成
部３と、この距離値画像情報作成部３から出力される各
距離値画像情報を取り込んで後述する線型補間処置によ
り各距離値画像の間に挿入される補間距離値画像情報を
求める線型補間部４と、距離値画像情報と補間距離値画
像情報とを用い被検体の対象領域体積を演算する体積演
算部５と、この装置の動作プログラムを格納したプログ
ラムメモリ６と、この装置で処理される前記各情報を記
憶する記憶部７と、前記動作プログラムをもとにこの装
置の各要素の動作を制御するCPU8とを有している。The object target area volume measuring device 1 shown in FIG. 1 is an X-ray CT.
The multi-slice image (continuous tomographic image) information of the subject output from the device or magnetic resonance diagnostic apparatus (MRI) is captured, and the contour of the target region of the subject in each slice image information is extracted by a conventionally known method. , A contour recognition unit 2 for recognizing a contour boundary, and a distance value for converting each slice image information into distance value image information based on the contour information of the target region of the subject recognized by the contour recognition unit 2. The image information creating unit 3 and the distance value image information output from the distance value image information creating unit 3 are fetched and interpolated distance value image information inserted between the distance value images is obtained by a linear interpolation process described later. A linear interpolation unit 4, a volume calculation unit 5 for calculating a target region volume of the subject using the distance value image information and the interpolation distance value image information, a program memory 6 storing an operation program of this apparatus, A storage unit 7 for storing each of the information processed by the location, and a CPU8 for controlling the operation of each element of the device on the basis of the operation program.
前記輪郭認識部２は、例えば第６図に示すように対象領
域を含むスライス領域のCT値情報を包含するマルチスラ
イス画像情報のうち、スライス番号i,i＋１の各スライ
ス画像情報における対象領域の輪郭情報（第６図におい
て太線で示す）を抽出する輪郭抽出部11と、この輪郭抽
出部11の抽出した輪郭情報を基にCT値０とCT値100とに
２値化する２値化処理部12とを具備している。すなわ
ち、輪郭認識部２において、被検体対象領域内の画素の
値は一定値100とされ、領域外の画素の値は一定値０と
される。For example, as shown in FIG. 6, the contour recognition unit 2 defines the contour of the target area in each slice image information of slice number i, i + 1 in the multi-slice image information including the CT value information of the slice area including the target area. A contour extraction unit 11 for extracting information (thick line in FIG. 6), and a binarization processing unit for binarizing CT value 0 and CT value 100 based on the contour information extracted by the contour extraction unit 11 It has 12 and. That is, in the contour recognition unit 2, the value of the pixel inside the subject area is set to the constant value 100, and the value of the pixel outside the area is set to the constant value 0.
前記距離値画像情報作成部３は、CPU8の制御のもとに前
記輪郭認識部２で認識した輪郭に沿った各画素の値（CT
値で100）を基準として対象領域に属する側の各画素に
対しては輪郭境界から上下、左右方向に一つ離れる毎に
CT値を＋１し、対象領域に属さない側の各画素に対して
は輪郭境界から上下、左右方向に一つ離れる毎にCT値を
−１していくように作用する。The distance value image information creation unit 3 controls the value (CT) of each pixel along the contour recognized by the contour recognition unit 2 under the control of the CPU 8.
For each pixel on the side that belongs to the target area with a value of 100) as a reference,
The CT value is incremented by 1, and the CT value is decremented by -1 for each pixel on the side that does not belong to the target area every time the pixel is separated from the contour boundary in the vertical and horizontal directions.
前記線型補間部４は、前記距離値画像情報作成部３で作
成した各距離値画像情報を取り込み、これらに対し次に
述べる線型補間処理を実行する。The linear interpolation unit 4 takes in each distance value image information created by the distance value image information creation unit 3 and executes a linear interpolation process described below on them.
スライス番号i,i＋１の２枚のスライスCT画像から得ら
れたスライス番号i,i＋１の距離値画像の間に線形補間
により補間距離値画像を作成する場合、補間距離値画像
の挿入位置は、スライス画像の空間分解能とスライス厚
によって決定される。When an interpolation distance value image is created by linear interpolation between the distance value images of slice numbers i, i + 1 obtained from two slice CT images of slice number i, i + 1, the insertion position of the interpolation distance value image is Determined by image spatial resolution and slice thickness.
例えば、空間分解能0.82mm、スライス厚1.64mmである場
合、補間距離値画像は両距離値画像の中間に１枚挿入す
ればよい。すなわち一般に補間距離値画像も含めた最終
的なスライス厚が空間分解能に等しくなるように補間が
行なわれる。For example, when the spatial resolution is 0.82 mm and the slice thickness is 1.64 mm, one interpolation distance value image may be inserted in the middle of both distance value images. That is, in general, interpolation is performed so that the final slice thickness including the interpolated distance value image becomes equal to the spatial resolution.
補間距離値画像の画素値Ｔ（x,y）すなわち補間距離値
画像情報は次の式（１）によって得られる： ここに、 Si（x,y）,Si＋１（x,y）：スライス番号ｉとｉ＋１の
距離値画像上位置（x,y）の画素の有する値、 l,m:スライス番号ｉの距離値画像および補間距離値画像
間の距離と、スライス番号ｉ＋１の距離値画像および補
間距離値画像間の距離との比を示す数、 である。The pixel value T (x, y) of the interpolated distance value image, that is, the interpolated distance value image information is obtained by the following equation (1): Where Si (x, y), Si + 1 (x, y): the value of the pixel at the position (x, y) on the distance value image of slice numbers i and i + 1, l, m: the distance value image of slice number i And a number indicating the ratio between the distance between the interpolated distance value images and the distance between the distance value image of slice number i + 1 and the interpolated distance value image.
次に上記構成の装置の作用を、第２図に示すこの装置の
動作手順を示すフローチャートと第３図および第４図に
示す距離値画像および補間距離値画像の説明図をも参照
して説明する。Next, the operation of the apparatus having the above-described configuration will be described with reference to the flowchart showing the operation procedure of this apparatus shown in FIG. 2 and the explanatory views of the distance value image and the interpolated distance value image shown in FIGS. 3 and 4. To do.
この装置にマルチスライスCT画像情報を入力すると（ST
1）、輪郭抽出部11は例えば第６図に示されるスライス
番号i,i＋１の各スライス画像情報における対象領域の
輪郭情報（第６図において太線を付して示す）を抽出す
る（ST2）。When multi-slice CT image information is input to this device (ST
1) The contour extracting unit 11 extracts, for example, the contour information (shown by thick lines in FIG. 6) of the target area in each slice image information of slice numbers i and i + 1 shown in FIG. 6 (ST2).
次に２値化処理部12は、輪郭抽出部11が抽出した輪郭情
報を含む前記各スライス画像情報を取込み、各スライス
画像情報の各画素の値を輪郭の外部において０とし、輪
郭の内部において一定値100とするように２値化し（ST
3）、これらをCPU8を介して距離値画像情報作成部３へ
送出する。Next, the binarization processing unit 12 takes in each slice image information including the contour information extracted by the contour extraction unit 11, sets the value of each pixel of each slice image information to 0 outside the contour, and inside the contour. Binarize so that the constant value is 100 (ST
3), these are sent to the distance value image information creation unit 3 via the CPU 8.
距離値画像情報作成部３は、CPU8の制御のもとに前記輪
郭情報における輪郭に沿った画素の値を基準として対象
領域に属する側の各画素に対して輪郭から上下、左右方
向に一つ離れる毎に＋１の値を加え、また、対象領域に
属しない側の各画素の値に対して輪郭から上下、左右方
向に一つ離れる毎に−１の値を加えていく（ST4,ST
5）。このような動作が各２値化スライス画像情報の全
ての画素について行われ、これにより、第３図の最下部
及び最上部に示すスライス番号i,i＋１等の各距離値画
像情報が得られる（ST6）。画素の値は、その画素と対
象領域輪郭との距離を示す。Under the control of the CPU 8, the distance value image information creation unit 3 uses the values of the pixels along the contour in the contour information as a reference for each pixel on the side belonging to the target area in the vertical and horizontal directions from the contour. A value of +1 is added each time the distance from the contour is increased, and a value of -1 is added to the value of each pixel on the side that does not belong to the target area in the vertical and horizontal directions from the contour (ST4, ST
Five). Such an operation is performed for all the pixels of each binarized slice image information, whereby each distance value image information such as slice numbers i, i + 1 shown at the bottom and top of FIG. 3 is obtained ( ST6). The pixel value indicates the distance between the pixel and the target area contour.
各距離値画像情報は、CPU8の制御のもとに線型補間部４
に送られ、ここで、既述した（１）式を用い、線型補間
処理が実行される（ST7）。これにより、第３図の中央
部分に示す補間距離値画像情報を求めることができる。
すなわち、今の場合（１）式においてｌ＝ｍ＝１である
から補間距離値画像上位置（x,y）の画素の値はスライ
ス番号ｉおよびｉ＋１の距離値画像上同じ位置（x,y）
の２つの画素の値の平均として求められる。Each distance value image information is stored in the linear interpolation unit 4 under the control of the CPU 8.
Then, the linear interpolation processing is executed using the above-mentioned equation (1) (ST7). As a result, the interpolated distance value image information shown in the central portion of FIG. 3 can be obtained.
That is, in this case, since l = m = 1 in the equation (1), the pixel value at the position (x, y) on the interpolated distance value image is the same position (x, y) on the distance value image at slice numbers i and i + 1. )
Is calculated as the average of the values of the two pixels.
この補間距離値画像情報における各画素の値（距離値）
は、スライス番号i,i＋１の各距離値画像情報の各画素
の値を平均化したものであることから、画素値が100で
ある画素の集合として定まるその輪郭は第３図に太線で
示し、また第４図に３次元的に視覚化して示すように各
距離値画像情報の輪郭の平均位置を占めることになる。
このようにして得られる距離値画像情報および補間距離
値画像情報における輪郭情報の全体は、実際の被検体対
象領域の形状に関し極めて正確な情報を与える。Value of each pixel in this interpolated distance value image information (distance value)
Is a value obtained by averaging the values of the respective pixels of the distance value image information of slice numbers i, i + 1, the contour determined as a set of pixels having a pixel value of 100 is shown by a bold line in FIG. Further, as shown three-dimensionally in FIG. 4, the average position of the contour of each distance value image information is occupied.
The entire contour information in the distance value image information and the interpolated distance value image information obtained in this way provides extremely accurate information regarding the actual shape of the subject region of interest.
このようにして求めた各補間距離値画像情報及び前記各
距離値画像情報は体積演算部５に送られ、ここで各画像
情報においてCT値100以上の画素の数が計数され、それ
らの総和が求められる。The interpolation distance value image information and the distance value image information thus obtained are sent to the volume calculation unit 5, where the number of pixels having a CT value of 100 or more in each image information is counted, and the sum of them is calculated. Desired.
その総和値に１つの画素の表示する実際の面積すなわち
本実施例にあっては0.82mm×0.82mmに、実効的なスライ
ス厚すなわち距離値画像と補間距離値画像との間の距離
0.82mmを乗じた値すなわちボクセルの体積を乗じること
により、被検体対象領域の体積を演算する（ST8）。す
なわち、被検体対象領域を多数のボクセル＝小立方体に
分割し、対象領域に含まれるボクセルの数を計数するこ
とによって体積を計測するボクセル法を用いる。The actual area displayed by one pixel in the total value, that is, 0.82 mm × 0.82 mm in this embodiment, is the effective slice thickness, that is, the distance between the distance value image and the interpolated distance value image.
The volume of the subject region is calculated by multiplying the value multiplied by 0.82 mm, that is, the voxel volume (ST8). That is, the voxel method is used in which the subject target area is divided into a large number of voxels = small cubes, and the volume is measured by counting the number of voxels included in the target area.
このように本実施例にあってはCTから入力されたスライ
ス画像情報から距離値画像情報を作成し、それを用いて
補間距離値画像情報を作成し、補間を行うので、対象領
域を実際にスライスして得たスライス画像の輪郭の段差
を被検体対象領域の実際の形状に近似し滑らかに補間す
ることができ、対象領域体積の正確な計測が可能となる
ものである。この方法は特に微小な対象領域の体積を計
測するときに有効である。In this way, in this embodiment, distance value image information is created from slice image information input from CT, interpolation distance value image information is created using it, and interpolation is performed, so that the target area is actually The contour step of the slice image obtained by slicing can be approximated to the actual shape of the subject region of interest and smoothly interpolated, and the volume of the subject region can be accurately measured. This method is particularly effective when measuring the volume of a minute target area.
もちろん、対象領域の体積の計測を行うのに、必らずし
も距離値画像を作成し、この距離値画像にもとずいて補
間を行う必要はなく、CT値を用いて補間を行い、スライ
ス画像情報と補間画像情報から輪郭認識手段を用い被写
体の領域の輪郭を抽出し、体積演算手段において各スラ
イス画像情報および補間画像情報の輪郭内に含まれる画
素の数を計数し、それらの総和をとることによって前記
実施例の場合と同様に対象領域の体積を計測することが
できる。Of course, in order to measure the volume of the target area, it is not always necessary to create a distance value image, and it is not necessary to perform interpolation based on this distance value image, but interpolation is performed using CT values, The contour of the subject area is extracted from the slice image information and the interpolation image information by using the contour recognition means, and the volume calculation means counts the number of pixels included in the contour of each slice image information and the interpolation image information, and sums them. By taking the above, the volume of the target region can be measured as in the case of the above embodiment.
次にこの方法をとる本発明の第２の実施例について述べ
る。本実施例にあっては、前記第１の実施例の場合と異
なり距離値画像に変換することなく、CTから入力される
スライスCT画像情報を用い線形補間により補間画像情報
を求めるものである。Next, a second embodiment of the present invention which adopts this method will be described. In the present embodiment, unlike the case of the first embodiment, interpolation image information is obtained by linear interpolation using slice CT image information input from CT without conversion to distance value images.
第５図に本実施例における処理手順を、第６図に前記第
１の実施例と同様のスライス画像情報および本実施例の
方法で得られる補間画像情報を、第７図に第６図に示さ
れる各画像情報における対象領域各輪郭の位置関係を示
す。FIG. 5 shows the processing procedure in this embodiment, FIG. 6 shows the same slice image information as in the first embodiment and interpolation image information obtained by the method of this embodiment, and FIG. The positional relationship of each outline of the target area in each image information shown is shown.
本実施例の被検体対象領域体積計測装置は補間画像情報
作成部と、輪郭認識部と、体積演算部等を有し、特許請
求の範囲（１）で述べた構成を有する。次に、第５図を
参照しながら本実施例における処理手順を説明する。The subject target area volume measuring apparatus according to the present embodiment has an interpolation image information creation unit, a contour recognition unit, a volume calculation unit, and the like, and has the configuration described in claim (1). Next, a processing procedure in this embodiment will be described with reference to FIG.
補間画像情報作成部ではCTから入力されたスライスCT画
像情報にもとずき、補間画像情報が作成される（ST1お
よびST2）。このとき、補間としては（１）式と同様な
線形補間が用いられ、ただ第１の実施例の場合と異な
り、画素のCT値について平均がとられ、その値が補間画
像上画素の値とされる。次に各スライス画像情報および
補間画像情報にもとずき輪郭認識部において被検体対象
領域の輪郭に対応する画素の抽出が行われる（ST3）。
この結果にもとずき、体積演算部において、前記第１の
実施例の場合と同様に、各画像上の輪郭内にある画素数
が各画像について計数され、さらにそれらの総和が求め
られる。こうして得られた総和値、すなわち被検体対象
領域内ボクセル数に、各ボクセルの体積、すなわち0.82
mm×0.82mm×0.82mmを乗じることによって被検体対象領
域の体積が得られる（ST4）。The interpolated image information creation unit creates interpolated image information based on the slice CT image information input from the CT (ST1 and ST2). At this time, the linear interpolation similar to the equation (1) is used as the interpolation, and unlike the case of the first embodiment, the CT values of the pixels are averaged, and the average is taken as the value of the pixel on the interpolated image. To be done. Next, based on the slice image information and the interpolated image information, the contour recognition unit extracts pixels corresponding to the contour of the subject region (ST3).
Based on this result, in the volume calculation unit, as in the case of the first embodiment, the number of pixels within the contour on each image is counted for each image, and the sum of them is obtained. The total value thus obtained, that is, the number of voxels in the object region of interest, the volume of each voxel, namely 0.82
The volume of the subject area is obtained by multiplying by mm × 0.82 mm × 0.82 mm (ST4).
ただこの方法では第７図に示されるようにスライス画像
情報の内容によっては補間がやや不自然な結果を与える
ことがある。すなわち、実際の対象領域の輪郭はなめら
かにつながっているのにもかかわらず、補間画像の輪郭
はスライス画像の輪郭間の段差を小さくするものとはな
っていない。したがってこの方法では前記第１の実施例
の場合に比べると被検体対象領域体積の計測に関し、や
や不正確な結果を与えるが、しかしいずれにしてもその
ような誤差は特に小さな領域を正確に計測する場合のほ
かは大きな問題となることはなく、本実施例の方法で多
くの場合十分な効果が得られる。However, in this method, the interpolation may give a somewhat unnatural result depending on the contents of the slice image information, as shown in FIG. That is, although the contours of the actual target area are smoothly connected, the contour of the interpolation image does not reduce the step between the contours of the slice image. Therefore, this method gives a slightly inaccurate result with respect to the measurement of the subject object region volume as compared with the case of the first embodiment, but in any case, such an error is accurately measured particularly in a small region. There is no major problem other than the case, and the method of this embodiment can obtain a sufficient effect in many cases.
また、本発明の実施態様としては上述した実施例に限定
されるものではなくその要旨の範囲内で種々の変形が可
能である。Further, the embodiment of the present invention is not limited to the above-described embodiments, and various modifications can be made within the scope of the gist thereof.
例えば、上述した実施例ではＸ線CT装置によるマルチス
ライス画像情報を用いる場合について説明したが、これ
に限らず、磁気共鳴診断装置によるマルチスライス画像
情報を用いても同様に実施できる。For example, in the above-described embodiment, the case of using the multi-slice image information by the X-ray CT apparatus has been described, but the present invention is not limited to this, and the same can be performed by using the multi-slice image information by the magnetic resonance diagnostic apparatus.
［発明の効果］ 以上詳述した本発明によれば、直接に被検体対象領域の
体積を計測可能な被検体対象領域体積計測装置を提供す
ることができる。[Effects of the Invention] According to the present invention described in detail above, it is possible to provide a subject target region volume measuring apparatus capable of directly measuring the volume of the subject target region.
また、スライス画像を距離値画像に変換した後画像補間
を行うようにすれば、より正確な被検体対象領域の体積
計測が可能となる。Further, if the image interpolation is performed after converting the slice image into the distance value image, more accurate volume measurement of the subject region of interest can be performed.
第１図は本発明の第１の実施例の被検体対象領域体積計
測装置の構成を示すブロック図、第２図は同じく第１の
実施例装置の動作を示すフローチャート、第３図は同じ
く第１の実施例装置における距離値画像情報及び補間距
離値画像情報を示す説明図、第４図は第３図の画像情報
に含まれる被検体対象領域の輪郭情報を３次元的に表現
した説明図、第５図は第２の実施例の被検体対象領域体
積計測装置の動作を示すフローチャート、第６図は同じ
く第２の実施例装置におけるスライスCT画像情報および
補間画像情報を示す説明図、第７図は第６図の画像情報
に含まれる被検体対象領域の輪郭情報を３次元的に表現
した説明図である。 １……対象領域体積計測装置 ２……輪郭認識部 ３……距離値画像情報作成部 ４……線型補間部FIG. 1 is a block diagram showing the configuration of a subject object region volume measuring apparatus according to the first embodiment of the present invention, FIG. 2 is a flow chart showing the operation of the first embodiment apparatus, and FIG. FIG. 4 is an explanatory view showing distance value image information and interpolated distance value image information in the apparatus of the first embodiment, and FIG. 4 is an explanatory view three-dimensionally expressing the contour information of the subject target area included in the image information of FIG. FIG. 5 is a flow chart showing the operation of the subject object region volume measuring apparatus of the second embodiment, and FIG. 6 is an explanatory diagram showing slice CT image information and interpolation image information in the same apparatus of the second embodiment. FIG. 7 is an explanatory diagram in which the contour information of the subject target area included in the image information of FIG. 6 is three-dimensionally expressed. 1 ... Target area volume measuring device 2 ... Contour recognition unit 3 ... Distance value image information creation unit 4 ... Linear interpolation unit
───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.6 識別記号 庁内整理番号 ＦＩ 技術表示箇所 Ｇ０１Ｂ 7/00 Ｚ 15/00 Ａ ─────────────────────────────────────────────────── ─── Continuation of front page (51) Int.Cl. 6 Identification code Office reference number FI technical display location G01B 7/00 Z 15/00 A
Claims (3)
入力し、隣接するスライス画像情報間に挿入すべきスラ
イス画像情報を補間演算により作成する補間画像情報作
成手段と、 各スライス画像情報および補間画像情報における被検体
の対象領域の輪郭を認識する輪郭認識手段と、 輪郭認識手段の認識結果に基づき被検体の対象領域の体
積の演算を行う体積演算手段を有することを特徴とする
被検体対象領域体積計測装置。1. Interpolation image information creating means for inputting multi-slice image information about an object and creating slice image information to be inserted between adjacent slice image information by interpolation calculation, and each slice image information and interpolation image information. The target object area volume characterized by having contour recognizing means for recognizing the contour of the target area of the subject and volume calculating means for calculating the volume of the target area of the subject based on the recognition result of the contour recognizing means. Measuring device.
入力し、各スライス画像情報における被検体の対象領域
の輪郭を認識する輪郭認識手段と、前記各スライス画像
の各画素と、前記輪郭認識手段により認識された対象領
域の輪郭部位との距離を求め、その距離に応じた画素値
を有する距離値画像情報を作成する距離値画像情報作成
手段と、距離値画像情報を補間して補間距離値画像情報
を作成する補間距離値画像情報作成手段と、輪郭認識手
段における認識結果と補間距離値画像情報とに基づき被
検体の対象領域の体積の演算を行う体積演算手段を有す
ることを特徴とする被検体対象領域体積計測装置。2. A contour recognizing unit for inputting multi-slice image information about a subject and recognizing a contour of a target region of the subject in each slice image information, each pixel of each slice image, and the contour recognizing unit. Distance value image information creating means for creating a distance value image information having a pixel value according to the distance, and an interpolated distance value image by interpolating the distance value image information. An interpolation distance value image information generating means for generating information, and a volume calculation means for calculating the volume of the target region of the subject based on the recognition result in the contour recognition means and the interpolation distance value image information. Specimen target area volume measuring device.
像の空間分解能Ｂに基づき各スライス画像間の補間画像
枚数をA/Bにて求めることを特徴とする特許請求の範囲
第１項または第２項記載の被検体対象領域体積計測装
置。3. The method according to claim 1, wherein the interval A between the slice images is recognized, and the number of interpolated images between the slice images is obtained by A / B based on the spatial resolution B of the images. Alternatively, the subject object region volume measuring apparatus according to the second item.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP63185441A JPH0760088B2 (en) | 1988-07-27 | 1988-07-27 | Target area volume measuring device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP63185441A JPH0760088B2 (en) | 1988-07-27 | 1988-07-27 | Target area volume measuring device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH0236305A JPH0236305A (en) | 1990-02-06 |
| JPH0760088B2 true JPH0760088B2 (en) | 1995-06-28 |
Family
ID=16170846
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP63185441A Expired - Fee Related JPH0760088B2 (en) | 1988-07-27 | 1988-07-27 | Target area volume measuring device |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH0760088B2 (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH03273118A (en) * | 1990-03-23 | 1991-12-04 | Daitoo:Kk | Method and system for measuring sediment loadage on vehicle |
| JPH06154226A (en) * | 1992-11-25 | 1994-06-03 | Matsushita Electric Ind Co Ltd | Ultrasonic diagnosing apparatus |
| CN110309561A (en) * | 2019-06-14 | 2019-10-08 | 吉旗物联科技（上海）有限公司 | Goods space volume measuring method and device |
-
1988
- 1988-07-27 JP JP63185441A patent/JPH0760088B2/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH0236305A (en) | 1990-02-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7822246B2 (en) | Method, a system and a computer program for integration of medical diagnostic information and a geometric model of a movable body | |
| US8423124B2 (en) | Method and system for spine visualization in 3D medical images | |
| EP1744671B1 (en) | Medical imaging system for accurate measurement evaluation of changes in a target lesion | |
| US6373998B2 (en) | Image processing electronic device for detecting dimensional variations | |
| Baert et al. | Three-dimensional guide-wire reconstruction from biplane image sequences for integrated display in 3-D vasculature | |
| CN107563998B (en) | Method for processing heart image in medical image | |
| US20090074276A1 (en) | Voxel Matching Technique for Removal of Artifacts in Medical Subtraction Images | |
| US20110075896A1 (en) | Computer readable medium, systems and methods for medical image analysis using motion information | |
| US8055044B2 (en) | Flexible 3D rotational angiography and computed tomography fusion | |
| WO2007079099A2 (en) | Cross-time and cross-modality medical diagnosis | |
| US20080205724A1 (en) | Method, an Apparatus and a Computer Program For Segmenting an Anatomic Structure in a Multi-Dimensional Dataset | |
| US20200051257A1 (en) | Scan alignment based on patient-based surface in medical diagnostic ultrasound imaging | |
| AU2018301580B2 (en) | Three-dimensional ultrasound image display method | |
| JP5296981B2 (en) | Automatic registration of medical volume images in modalities using affine transformation | |
| CN116703882A (en) | Non-contact type heart three-dimensional mapping method and system | |
| US8009910B2 (en) | Method, a system, a computer program product and a user interface for segmenting image sets | |
| US20200286240A1 (en) | Method of segmenting a 3d object in a medical radiation image | |
| CN107292351B (en) | A kind of matching process and device of tubercle | |
| JPH0760088B2 (en) | Target area volume measuring device | |
| US8165375B2 (en) | Method and system for registering CT data sets | |
| WO2008002325A2 (en) | Cross-time inspection method for medical diagnosis | |
| Veisterä et al. | Reconstructing 3D boundary element heart models from 2D biplane fluoroscopy | |
| Tsai et al. | Computer-based measurement of right and left ventricular volumes in magnetic resonance imaging | |
| Veistera et al. | Creating 3D boundary element models of the heart from 2D projections | |
| Patra et al. | Medical Image Processing in Nuclear Medicine and Bone Arthroplasty |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| LAPS | Cancellation because of no payment of annual fees |