JP2005261440A - Method for extracting imaged area of image for medical purpose - Google Patents

Method for extracting imaged area of image for medical purpose Download PDF

Info

Publication number
JP2005261440A
JP2005261440A JP2004052302A JP2004052302A JP2005261440A JP 2005261440 A JP2005261440 A JP 2005261440A JP 2004052302 A JP2004052302 A JP 2004052302A JP 2004052302 A JP2004052302 A JP 2004052302A JP 2005261440 A JP2005261440 A JP 2005261440A
Authority
JP
Japan
Prior art keywords
image
region
contrast
low
blood vessel
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.)
Withdrawn
Application number
JP2004052302A
Other languages
Japanese (ja)
Inventor
Jun Masumoto
潤 桝本
Shiyougo Azemoto
将吾 畦元
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Canon Medical Systems Corp
Original Assignee
Aze Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Aze Ltd filed Critical Aze Ltd
Priority to JP2004052302A priority Critical patent/JP2005261440A/en
Publication of JP2005261440A publication Critical patent/JP2005261440A/en
Withdrawn legal-status Critical Current

Links

Images

Landscapes

  • Measurement Of Radiation (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Nuclear Medicine (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)
  • Conversion Of X-Rays Into Visible Images (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To easily and excellently extract only an imaged area on the basis of a low imaging image picked up in a time zone in which an imaging effect is relatively low and a high imaging image picked up in a time zone in which the imaging effect is relatively high. <P>SOLUTION: An dilation processing for multilevel images is executed to the low imaging image 11, and an image 13 after the dilation processing in which a bony part area 15 is dilated is obtained. By obtaining a difference image 14 by executing the subtraction processing of the obtained image 13 after the dilation processing and the high imaging image 12, a cerebrovascular area 17 is extracted. <P>COPYRIGHT: (C)2005,JPO&NCIPI

Description

本発明は主に、CT(computed tomography),MRI(magnetic resonance imaging),核医学,CR(computed radiography),DR(real-time radiography)等の放射線診断システムを用いて撮像された医療用画像の画像処理技術に関し、詳しくは、例えば頭部や頸部等における造影剤注入[投与]前後の医療用画像に基づいて、造影された血管等の領域(被造影領域)を抽出して画像化する医療用画像の被造影領域抽出方法に関する。   The present invention mainly relates to medical images captured using radiation diagnostic systems such as CT (computed tomography), MRI (magnetic resonance imaging), nuclear medicine, CR (computed radiography), and DR (real-time radiography). More specifically, for example, based on medical images before and after contrast medium injection [administration] in the head, neck, etc., the image processing technique is extracted and imaged by contrasting regions such as blood vessels (contrast region). The present invention relates to a method for extracting a contrast-enhanced region of a medical image.

従来、X線血管造影検査において、CTやMRI等の放射線診断システムを用いて撮影された医療用画像に基づき、造影された静脈や動脈血管のみを抽出して表示するDSA(digital subtraction angiography)と称される手法が知られている。   Conventionally, in X-ray angiography, DSA (digital subtraction angiography) which extracts and displays only contrasted veins and arterial blood vessels based on medical images taken using a radiation diagnostic system such as CT or MRI The technique called is known.

一般にこのDSAにおいては、造影効果が相対的に低い時間帯に撮像された低造影画像(通常は、造影前すなわち造影剤注入前や造影剤注入直後に撮影された画像、マスク画像とも称される)と、造影効果が相対的に高い時間帯に撮像された高造影画像(通常は、造影後すなわち造影剤が血管内に満たされる時間帯に撮影された画像、コントラスト画像とも称される)との減算処理(サブトラクション処理)によって差分画像を求めることにより、造影された血管領域のみを抽出するようになっている(下記特許文献1参照)。   In general, in this DSA, a low-contrast image captured in a time zone with a relatively low contrast effect (usually also referred to as an image captured before contrasting, that is, before contrast agent injection or immediately after contrast agent injection, or a mask image) ) And a high-contrast image captured in a time zone in which the contrast effect is relatively high (usually referred to as a contrast image or an image taken in a time zone after contrast enhancement, that is, in which a contrast medium is filled in a blood vessel) Thus, only the contrasted blood vessel region is extracted by obtaining a difference image by subtraction processing (subtraction processing) (see Patent Document 1 below).

一方、コンピュータ・グラフィックス等の画像処理技術は、近年著しい進展を見せており、外科的処置を伴わずに組織の構造変化や機能障害、病変部の有無等を検知することを可能にする技術として、医療分野への適用が経年的に拡大する傾向にある。
例えば、最新の画像処理技術をDSA等により得られた原画像に適用することによって、抽出された血管領域の画像をカラー化したり、抽出された血管領域の全体を様々な方向から観察できるように画像を3次元化したりすることが可能となっている。
On the other hand, image processing technologies such as computer graphics have made remarkable progress in recent years, and it is possible to detect structural changes, functional disorders, presence or absence of lesions, etc. without surgical treatment. As a result, application to the medical field tends to expand over time.
For example, by applying the latest image processing technology to the original image obtained by DSA or the like, the image of the extracted blood vessel region can be colored, or the entire extracted blood vessel region can be observed from various directions. It is possible to make an image three-dimensional.

このような画像処理技術の一手法として、グレイスケールダイレーションと称される手法が知られている。このグレイスケールダイレーション手法は、従来原画像を2値化処理してから行なわれていた拡張処理を、多値画像に対して直接的に適用できるようにしたものであり、2値化処理する場合には必要となる閾値の設定が不要となるなどの利点がある。この手法は、画像処理の分野においては以前より知られているものであるが(下記非特許文献1参照)、これを医療分野に適用した例は知られていない。   As one method of such image processing technology, a method called gray scale dilation is known. In this gray scale dilation method, an extension process that has been performed after binarizing the original image can be directly applied to the multi-valued image. In this case, there is an advantage that it is not necessary to set a necessary threshold value. This technique has been known for a long time in the field of image processing (see Non-Patent Document 1 below), but no example of applying this technique to the medical field is known.

特開平10−322597号公報Japanese Patent Laid-Open No. 10-322597 S. Sternberg, "Grayscale Morphology", Computer Vision, Graphics, and Image Processing, vol. 35, pp. 333-355, 1986.S. Sternberg, "Grayscale Morphology", Computer Vision, Graphics, and Image Processing, vol. 35, pp. 333-355, 1986.

これまで、頭部を撮影した医療用画像から脳血管の領域のみを抽出することは困難であると考えられていた。しかし近年、計算機の処理能力や得られる画像の解像度が大幅に向上していることを受け、DSA等の手法を頭部に適用することにより、脳血管の領域を抽出して画像化する技術の提案が試験的になされている。   Until now, it has been considered difficult to extract only the cerebral blood vessel region from the medical image obtained by photographing the head. However, in recent years, the processing power of computers and the resolution of images obtained have been greatly improved, and a technique for extracting and imaging a region of the cerebral blood vessels by applying a technique such as DSA to the head. Proposals have been made on a trial basis.

これまで提案されている手法は、脳血管の領域を抽出する際、造影前の画像と造影後の画像とを単純に減算処理して差分画像を求めるものである。しかし、このような手法では、造影前後の画像の位置合せを高精度に行なったとしても、パーシャルボリューム効果等の影響によって、差分画像上に不要な領域がゴミのように残ってしまう。   In the method proposed so far, when a cerebral blood vessel region is extracted, a difference image is obtained by simply subtracting the pre-contrast image and the post-contrast image. However, with such a technique, even if the images before and after contrast are aligned with high accuracy, an unnecessary area remains on the difference image as dust due to the influence of the partial volume effect or the like.

このような不要領域は、頭蓋骨と他の組織や液体部分との境界部分に連続的に残存するため、そのままでは脳血管の状態を観察することが困難となり診断に支障が生じてしまう。したがって、これまでの提案手法では、差分画像に残存するこれらの不要領域を除去する画像処理を行なっている。しかし、脳全体を覆う頭蓋骨の領域に沿って連続的に残存する不要領域を脳血管の領域と峻別することは容易ではなく、また残存する範囲も広範囲に及ぶために、不要領域除去の画像処理を行なった画像に対する信頼性が十分に得られず、またそのような画像処理に要する計算量や時間も甚大になるという問題がある。   Such an unnecessary region remains continuously at the boundary portion between the skull and other tissues or liquid portions, so that it is difficult to observe the state of the cerebral blood vessel as it is, and this causes a trouble in diagnosis. Therefore, in the proposed methods so far, image processing for removing these unnecessary areas remaining in the difference image is performed. However, it is not easy to distinguish the unnecessary area that remains continuously along the skull area that covers the entire brain from the cerebral blood vessel area, and the remaining area covers a wide range. However, there is a problem that sufficient reliability is not obtained with respect to the image subjected to the image processing, and the calculation amount and time required for such image processing are enormous.

このような問題は、頭蓋骨が脳全体を覆う頭部において顕著に生じるものであるが、頸部や脊椎部、あるいは腰部や大腿部、胸部など、骨部の領域と血管等の被造影領域とが近接する部位においても程度の差はあるものの、同様に生じるものである。   Such a problem occurs remarkably in the head where the skull covers the entire brain, but the region of the bone and the region to be imaged such as the blood vessel such as the neck, spine, lumbar, thigh, and chest. Even in a region where the and are close to each other, there is a difference in degree, but the same occurs.

本発明は、このような事情に鑑みなされたものであり、低造影画像と高造影画像とに基づいて、被造影領域のみを容易かつ良好に抽出することが可能な医療用画像の被造影領域抽出方法を提供することを目的とする。   The present invention has been made in view of such circumstances, and a contrast-enhanced region of a medical image capable of easily and satisfactorily extracting only the contrast-enhanced region based on a low-contrast image and a high-contrast image. An object is to provide an extraction method.

上記課題解決のため本発明方法においては、低造影画像に多値画像用の拡張処理を施して骨部領域を拡張し、これをマスク画像として高造影画像との減算処理を行なって差分画像を求め、この差分画像に基づいて被造影領域を抽出する。   In order to solve the above problem, in the method of the present invention, a low contrast image is subjected to multivalue image expansion processing to expand the bone region, and this is used as a mask image to perform subtraction processing from the high contrast image to obtain a difference image. The contrast-enhanced region is extracted based on the difference image.

すなわち、本発明に係る医療用画像の被造影領域抽出方法は、造影効果が相対的に低い時間帯に撮像された低造影画像と、造影効果が相対的に高い時間帯に撮像された高造影画像とに基づいて、被造影領域を抽出する医療用画像の被造影領域抽出方法において、多値画像用の拡張処理を前記低造影画像に施して、該低造影画像における骨部領域を拡張し、前記高造影画像と前記拡張処理後の前記低造影画像との減算処理によって差分画像を求めることにより、前記被造影領域を抽出することを特徴とする。   That is, the method for extracting a contrast-enhanced region of a medical image according to the present invention includes a low-contrast image captured in a time zone with a relatively low contrast effect and a high-contrast image captured in a time zone with a relatively high contrast effect. In a method for extracting a contrast-enhanced region of a medical image based on the image, the low-contrast image is subjected to a multi-value image expansion process to expand a bone region in the low-contrast image. The contrast region is extracted by obtaining a difference image by subtracting the high contrast image and the low contrast image after the expansion process.

本発明において、前記差分画像において抽出された前記被造影領域のうち縮減された必要部分に対しては、領域の復元処理を行なうことが好ましい。また、前記拡張処理としては、上記非特許文献1に記載されたグレイスケールダイレーション手法を用いることができる。   In the present invention, it is preferable to perform a region restoration process on a reduced necessary portion of the contrasted region extracted in the difference image. As the expansion process, the gray scale dilation method described in Non-Patent Document 1 can be used.

本発明は、頭部,頸部,脊椎部,腰部,大腿部,胸部など、骨部の領域と血管等の被造影領域とが近接する部位を撮影した医療用画像に対して適用することが有効であるが、特に頭蓋骨が脳全体を覆う頭部、すなわち被造影領域が脳血管の領域(腫瘍を含む)である場合の有効性が顕著である。   The present invention is applied to a medical image obtained by imaging a region where a bone region and a contrasted region such as a blood vessel are close to each other, such as the head, neck, spine, waist, thigh, and chest. However, the effectiveness is particularly remarkable when the skull covers the entire brain, that is, the region to be imaged is a cerebral blood vessel region (including a tumor).

なお、上記「造影効果が相対的に低い時間帯」は、造影効果が全く得られない時間帯、例えば造影剤注入または投与前の時間帯も含む。   The “time zone in which the contrast effect is relatively low” includes a time zone in which no contrast effect is obtained, for example, a time zone before contrast agent injection or administration.

本発明に係る医療用画像の被造影領域抽出方法によれば、多値画像用の拡張処理を施して骨部領域を拡張した低造影画像と、高造影画像との減算処理を行なって差分画像を求めることにより、差分画像に残存する不要領域を低減させることができるので、被造影領域のみを容易かつ良好に抽出することが可能となる。   According to the method for extracting a contrast-enhanced region of a medical image according to the present invention, a difference image is obtained by performing a subtraction process between a high-contrast image and a low-contrast image obtained by performing a multi-value image expansion process to expand a bone region. Since it is possible to reduce unnecessary areas remaining in the difference image, it is possible to easily and satisfactorily extract only the contrast target area.

以下、本発明の一実施形態に係る医療用画像の被造影領域抽出方法(以下「本実施形態方法」と称する)について、図面を参照しながら説明する。なお、以下の説明は、造影剤を用いてCT撮影された頭部における2時相の画像系列(低造影画像および高造影画像)のデータに基づき、脳血管領域を抽出して画像化する場合を例にとって行なう。   Hereinafter, a method for extracting a contrast region of a medical image according to an embodiment of the present invention (hereinafter referred to as “method of the present embodiment”) will be described with reference to the drawings. In the following explanation, a cerebral blood vessel region is extracted and imaged based on data of a two-phase image sequence (low contrast image and high contrast image) in the head taken by CT using a contrast agent. For example.

本実施形態方法の実施手順は、大きく次の3つの手順に分けることができる。すなわち、(1)画像間の位置合せ手順、(2)多値画像用の拡張処理を用いた減算処理手順、(3)縮減された必要部分の領域復元処理手順の3つの手順である。   The implementation procedure of the method of this embodiment can be roughly divided into the following three procedures. That is, there are three procedures: (1) a registration procedure between images, (2) a subtraction processing procedure using an expansion process for a multi-valued image, and (3) an area restoration processing procedure for a reduced necessary part.

(1)画像間の位置合せ手順は、造影効果が相対的に低い時間帯(例えば、造影剤注入前の時間帯)に撮像された低造影画像と、造影効果が相対的に高い時間帯(例えば、造影剤が脳血管全体に行き渡る時間帯)に撮像された高造影画像との間の座標上の位置合せを行なう手順であり、本実施形態方法では剛体レジストレーション手法を用いて行なう。   (1) The registration procedure between images includes a low-contrast image captured in a time zone in which the contrast effect is relatively low (for example, a time zone before contrast agent injection) and a time zone in which the contrast effect is relatively high ( For example, this is a procedure for performing coordinate alignment with a high-contrast image captured in a time zone in which the contrast agent spreads over the entire cerebral blood vessel. In the method of this embodiment, a rigid registration method is used.

剛体レジストレーション手法とは、位置合せの対象となる物体(この例では頭部)が剛体、すなわち物体の向きや位置は変化するが物体の形状や大きさは変化しないと仮定して、平行移動と回転運動の6つの運動パラメータを推定して位置合せを行なう手法である。本実施形態方法で用いる剛体レジストレーション手法においては、低造影画像と高造影画像との2つの画像系列間で、理論的にはCT値(濃度値)が変化しない領域(脳組織,脂肪,骨)の相関係数が最大となるように上記運動パラメータを求めるようになっている。   Rigid registration method means that the object to be aligned (in this example, the head) is a rigid body, that is, the direction and position of the object will change, but the shape and size of the object will not change. This is a method for performing alignment by estimating six motion parameters of rotational motion. In the rigid registration method used in the method of the present embodiment, a region in which the CT value (concentration value) does not change theoretically between two image series of a low contrast image and a high contrast image (brain tissue, fat, bone) The motion parameters are obtained so that the correlation coefficient of

(2)多値画像用の拡張処理を用いた減算処理手順は、位置合せが行なわれた低造影画像と高造影画像との減算処理により差分画像を求めるものであり、減算処理の際に低造影画像に多値画像用の拡張処理を行ない、骨部領域を拡張させる点に特徴がある。また、本実施形態方法では、多値画像用の拡張処理として前掲の非特許文献1に記載されたグレイスケールダイレーション手法を用いる。   (2) The subtraction processing procedure using the extended processing for multi-valued images is to obtain a difference image by subtracting the low contrast image and the high contrast image that have been aligned, and is low during the subtraction processing. The contrast image is characterized in that a multi-value image expansion process is performed to expand the bone region. In the method of the present embodiment, the gray scale dilation method described in Non-Patent Document 1 described above is used as the expansion processing for multi-valued images.

この手順の特徴について、図1および図2を用いてより詳細に説明する。図1は多値画像用の拡張処理を用いた減算処理手順を模式的に示す図、図2はグレイスケールダイレーション手法を模式的に示す図である。   The feature of this procedure will be described in more detail with reference to FIGS. FIG. 1 is a diagram schematically showing a subtraction processing procedure using an expansion process for multi-valued images, and FIG. 2 is a diagram schematically showing a gray scale dilation method.

図1に示すように低造影画像11においては、主に頭蓋骨の領域に相当する骨部領域15の濃度値が非常に高く、脳組織や脂肪、脳血管といった非骨領域16の濃度値は低くなっている。一方、高造影画像12においては、上記骨部領域15に加えて上記非骨領域16のうち造影された脳血管領域17の濃度値も高くなっている。   As shown in FIG. 1, in the low-contrast image 11, the density value of the bone region 15 mainly corresponding to the skull region is very high, and the density value of the non-bone region 16 such as brain tissue, fat, and cerebral blood vessel is low. It has become. On the other hand, in the high contrast image 12, in addition to the bone part region 15, the concentration value of the contrasted cerebral blood vessel region 17 in the non-bone region 16 is also high.

ここで、低造影画像11に上記グレイスケールダイレーション手法を用いた拡張処理を行ない、拡張処理後画像13を得る。このグレイスケールダイレーション手法は、後述するように、画像上の各局所において高い濃度値を有する画素の当該濃度値を、その周辺の画素に拡散させる作用を奏する。これにより、低造影画像11において高い濃度値を有する骨部領域15が、拡張処理後画像13においてはその全体に占める割合が増大した状態、すなわち拡張された状態で描出される。   Here, the low contrast image 11 is subjected to the expansion process using the gray scale dilation method, and the expanded image 13 is obtained. As will be described later, this gray scale dilation method has an effect of diffusing the density value of a pixel having a high density value in each local area on the image to surrounding pixels. As a result, the bone region 15 having a high density value in the low contrast image 11 is depicted in a state in which the ratio of the entire bone region 15 in the expanded image 13 is increased, that is, in an expanded state.

次に、高造影画像12と拡張処理後画像13との減算処理を行ない、差分画像14を求める。拡張処理後画像13においては骨部領域15が拡張されているため、単純に高造影画像12と低造影画像11との減算処理を行なう従来手法では頭蓋骨に沿って多く残存することになる不要領域が、本手法により得られた差分画像14にほとんど存在せず、脳血管領域17が良好に抽出される。   Next, the high contrast image 12 and the post-expansion image 13 are subtracted to obtain a difference image 14. Since the bone region 15 is expanded in the post-expansion image 13, an unnecessary region that remains much along the skull in the conventional method in which the high contrast image 12 and the low contrast image 11 are simply subtracted. However, it hardly exists in the difference image 14 obtained by this method, and the cerebral blood vessel region 17 is well extracted.

上記グレイスケールダイレーション手法は、具体的には以下のように行なわれる。すなわち図2に示すように、拡張処理を行なう第1画像41上を、所定の大きさ(図2では3×3正方)の構造化要素(カーネル,オペレータ)42を1通り順次移動させる。その際、構造化要素42の中心に位置する画素と、それと上下左右で隣接する4つの画素(図中の構造化要素42において斜線を付した4つの領域に位置する画素)とがそれぞれに有する濃度値を比較し、そのうちで最も高い濃度値を、中心の画素が有する濃度値と置き換える。このような拡張処理を画像全体に対して行なうことにより、拡張処理後の第2画像43上においては高い濃度値を有する範囲が周囲に拡散された状態となる。   Specifically, the gray scale dilation method is performed as follows. That is, as shown in FIG. 2, a structured element (kernel, operator) 42 having a predetermined size (3 × 3 square in FIG. 2) is sequentially moved one by one on the first image 41 to be expanded. At that time, each pixel has a pixel located in the center of the structuring element 42 and four pixels adjacent to it in the vertical and horizontal directions (pixels located in the four regions shaded in the structuring element 42 in the figure). The density values are compared, and the highest density value among them is replaced with the density value of the central pixel. By performing such an expansion process on the entire image, a range having a high density value is diffused around the second image 43 after the expansion process.

(3)縮減された必要部分の領域復元処理手順は、上述したグレイスケールダイレーション手法を用いた減算処理によって求められた差分画像において、脳血管領域のうち実際の大きさよりも縮減されて描出された部分(必要部分)の領域を復元する手順であり、このようなことが生じる場合に行なわれる。   (3) The region restoration processing procedure of the reduced necessary portion is depicted by reducing the actual size of the cerebral blood vessel region in the difference image obtained by the subtraction processing using the gray scale dilation method described above. This is a procedure for restoring the area of the necessary part (necessary part), and is performed when such a situation occurs.

脳血管の領域が実際の大きさよりも縮減されて描出される虞があるのは、主に脳血管と頭蓋骨とが近接している部分においてである。以下では、領域復元処理を、(3a)脳血管が頭蓋骨を貫通している部分に対する処理と、(3b)脳血管が頭蓋骨に沿うように接している部分に対する処理とに分けて説明する。なお、図3は脳血管と頭蓋骨とが近接している部分を例示する図で、同図(a)は脳血管が頭蓋骨を貫通している部分(図中円内の部分)を示し、同図(b)は脳血管が頭蓋骨に沿うように接している部分(図中楕円内の部分)を示している。   There is a possibility that the region of the cerebral blood vessel is drawn with a size smaller than the actual size mainly in a portion where the cerebral blood vessel and the skull are close to each other. In the following, the region restoration process will be described by dividing it into (3a) a process for a part where the cerebral blood vessel penetrates the skull and (3b) a process for a part where the cerebral blood vessel is in contact with the skull. FIG. 3 is a diagram illustrating a portion where the cerebral blood vessel and the skull are close to each other. FIG. 3A shows a portion where the cerebral blood vessel penetrates the skull (the portion in the circle in the figure). FIG. (B) shows a portion where the cerebral blood vessel is in contact with the skull (the portion within the ellipse in the figure).

(3a)脳血管が頭蓋骨を貫通している部分に対する処理は、図4に示すような手順で行なわれる。図4は脳血管が頭蓋骨を貫通している部分に対する処理手順を模式的に示す図である。   (3a) The process for the portion where the cerebral blood vessel penetrates the skull is performed according to the procedure shown in FIG. FIG. 4 is a diagram schematically showing a processing procedure for a portion where the cerebral blood vessel penetrates the skull.

図4に示す低造影画像21においては、濃度値が非常に高い骨部領域25の中に濃度値が低い非骨領域26の存在が確認される。この非骨領域26が骨部領域25を貫通する脳血管領域27である場合、高造影画像22においては、この脳血管領域27が造影されるために高い濃度値を有することになる。   In the low contrast image 21 shown in FIG. 4, the existence of a non-bone region 26 having a low density value is confirmed in the bone region 25 having a very high density value. When the non-bone region 26 is a cerebral blood vessel region 27 penetrating the bone region 25, the high contrast image 22 has a high density value because the cerebral blood vessel region 27 is imaged.

ここで、低造影画像21に上記拡張処理を行なうと、拡張処理後画像23においては脳血管領域27が縮減されて描出される。この状態で、高造影画像22と拡張処理後画像23との減算処理を行なうと、差分画像24においては、脳血管領域27が実際の大きさよりも縮減されて小さく描出されることになる。そこで、脳血管領域27の周辺拡大画像28に示すように、描出された脳血管領域27aの周囲に存在する縮減された脳血管領域27bの範囲を特定し、その領域を復元する処理を行なう。   Here, when the expansion process is performed on the low contrast image 21, the cerebral blood vessel region 27 is reduced in the post-expansion image 23 and is rendered. In this state, when the high contrast image 22 and the post-expansion image 23 are subtracted, the cerebral blood vessel region 27 is drawn smaller than the actual size in the difference image 24. Therefore, as shown in the enlarged peripheral image 28 of the cerebral vascular region 27, the range of the reduced cerebral vascular region 27b existing around the rendered cerebral vascular region 27a is specified, and the region is restored.

縮減された脳血管領域27bの特定および復元は、例えば差分画像24上において所定の濃度値αよりも大きい濃度値Vを有する領域(V>α)を縮減された疑いのある領域として特定し、その特定された領域については、低造影画像21と高造影画像22との減算処理を行なうことにより、その領域を復元する方法を用いることができる。 For example, the reduced cerebrovascular region 27b is identified and restored by specifying a region (V 1 > α) having a density value V 1 larger than a predetermined density value α on the difference image 24 as a region suspected of being reduced. And about the specified area | region, the method of decompress | restoring the area | region by performing the subtraction process of the low contrast image 21 and the high contrast image 22 can be used.

(3b)脳血管が頭蓋骨に沿うように接している部分に対する処理は、図5に示すような手順で行なわれる。図5は脳血管が頭蓋骨に沿うように接している部分に対する処理手順を模式的に示す図である。   (3b) The processing for the portion where the cerebral blood vessel is in contact with the skull is performed according to the procedure shown in FIG. FIG. 5 is a diagram schematically showing a processing procedure for a portion where the cerebral blood vessel is in contact with the skull.

図5に示す低造影画像31においては、濃度値が非常に高い骨部領域35に沿って濃度値が低い非骨領域36の存在が確認される。この非骨領域36が骨部領域35に沿うように接している脳血管領域37を含む場合、高造影画像32においては、この脳血管領域37が造影されて高い濃度値を有することになる。   In the low contrast image 31 shown in FIG. 5, the existence of the non-bone region 36 having a low density value is confirmed along the bone region 35 having a very high density value. When the non-bone region 36 includes a cerebral blood vessel region 37 that is in contact with the bone region 35, the high blood contrast image 32 has a high density value in the high contrast image 32.

ここで、低造影画像31に上記拡張処理を行なうと、拡張処理後画像33においては脳血管領域37が縮減されて描出される。この状態で、高造影画像32と拡張処理後画像33との減算処理を行なうと、差分画像34においては、脳血管領域37が実際の大きさよりも縮減されて細く描出されることになる。そこで、描出された脳血管領域37aの傍らに存在する縮減された脳血管領域37bの範囲を特定し、その領域を復元する処理を行なう。   Here, when the expansion process is performed on the low contrast image 31, the cerebral blood vessel region 37 is reduced in the post-expansion image 33 and is rendered. When the high contrast image 32 and the post-expansion image 33 are subtracted in this state, the cerebral blood vessel region 37 is reduced in size and drawn finely in the difference image 34. Therefore, the range of the reduced cerebral vascular region 37b existing beside the rendered cerebral vascular region 37a is specified, and processing for restoring the region is performed.

縮減された脳血管領域37bの特定は、例えば拡張処理後画像33上において所定の濃度値βよりも大きい濃度値Vを有し、かつ低造影画像31上においても所定の濃度値γより大きい濃度値Vを有する領域(V>βかつV>γ)を縮減された疑いのある領域として特定する方法を用いることができる。また、その特定された領域の復元については、例えば低造影画像21と高造影画像22との減算処理を行なった結果と、差分画像34との結果を所定の割合で加算処理する方法を用いることができる。 Specific reduction brain blood vessel region 37b, for example greater than a predetermined density value γ even extended on processed image 33 has a larger density value V 2 than the predetermined density value beta, and low contrast images on 31 A method of specifying a region having the density value V 3 (V 2 > β and V 3 > γ) as a region suspected of being reduced can be used. In addition, for the restoration of the specified region, for example, a method of adding the result of the subtraction process between the low contrast image 21 and the high contrast image 22 and the result of the difference image 34 at a predetermined ratio is used. Can do.

これらの処理を行なった後、画像上に残る細かい不要領域を除去するとともに、所定体積以上の連結領域を最終的に脳血管の領域として確定する。なお、図6は脳血管領域の確定処理前(同図(a))と処理後(同図(b))の画像を示す図である。   After performing these processes, a fine unnecessary region remaining on the image is removed, and a connection region having a predetermined volume or more is finally determined as a cerebral blood vessel region. FIG. 6 is a diagram showing images before and after the cerebrovascular region determination process (FIG. 6A) and after the process (FIG. 6B).

また、図7は従来手法により得られる差分画像(同図(a))と本実施形態方法により得られる差分画像(同図(b))とを示す図である。この図7に示すように、造影前の画像と造影後の画像とを単純に減算処理して差分画像を求める従来手法では、造影前後の画像の位置合せを高精度に行なったとしても、パーシャルボリューム効果等の影響によって、差分画像上に不要な領域がゴミのように残ってしまう。一方、造影前の画像に対してグレイスケールダイレーション手法を用いた拡張処理を行なって骨部領域を拡張させ、この拡張処理後の画像と造影後の画像との減算処理により差分画像を求める本実施形態方法では、差分画像上に不要な領域がほとんど残らず、脳血管の領域を良好に抽出して画像化することが可能となる。   FIG. 7 is a diagram showing a difference image (FIG. (A)) obtained by a conventional method and a difference image (FIG. (B)) obtained by the method of the present embodiment. As shown in FIG. 7, in the conventional method in which the difference image is obtained by simply subtracting the pre-contrast image and the post-contrast image, even if the images before and after the contrast are aligned with high accuracy, Due to the effect of the volume effect or the like, an unnecessary area remains on the difference image like dust. On the other hand, the bone region is expanded by performing an expansion process using the gray scale dilation method on the pre-contrast image, and a difference image is obtained by subtracting the image after the expansion process and the post-contrast image. In the embodiment method, there is almost no unnecessary area remaining on the difference image, and the cerebral blood vessel area can be satisfactorily extracted and imaged.

以上、本発明の実施形態について説明したが、本発明は上記実施形態に限定されるものではなく、種々に態様を変更することが可能である。   As mentioned above, although embodiment of this invention was described, this invention is not limited to the said embodiment, A various aspect can be changed.

例えば上記実施形態では、造影剤を用いたCT撮影により得られた頭部の画像データに基づき、脳血管の領域を抽出しているが、撮影手法はCTに限られずまた本発明の適用部位は頭部に限られない。撮影手法としてはCT以外に、MRI,核医学,DR等の種々の放射線診断システムを用いることが可能である。他の撮影手法で得られた画像データに本発明を適用するのに際して特段の困難はない。同一領域に関する造影前後の画像が得られさえすれば、上述した本発明の手順を撮影手法に関係なく略同様に適用することができる。   For example, in the above embodiment, the cerebral blood vessel region is extracted based on the image data of the head obtained by CT imaging using a contrast agent. However, the imaging method is not limited to CT, and the application site of the present invention is as follows. It is not limited to the head. In addition to CT, various radiological diagnosis systems such as MRI, nuclear medicine, and DR can be used as imaging methods. There is no particular difficulty in applying the present invention to image data obtained by other imaging methods. As long as images before and after contrast enhancement relating to the same region are obtained, the above-described procedure of the present invention can be applied in substantially the same manner regardless of the imaging technique.

本発明の適用部位に関しては、頭部以外に頸部,脊椎部,腰部,大腿部,胸部など、骨部の領域と血管等の被造影領域とが近接する種々の部位を、その対象とすることが可能である。   Regarding the application site of the present invention, in addition to the head, various sites where the bone region and the contrasted region such as a blood vessel are close to each other, such as the neck, spine, waist, thigh, and chest, are considered as targets. Is possible.

また、上記実施形態では多値画像用の拡張処理手法としてグレイスケールダイレーション手法を用いており、その際の構造化要素として3×3正方のものを例示しているが、構造化要素の大きさや形、あるいは構造化要素内での濃度値の置換方法を変えることが可能であり、また、上記非特許文献1に記載されたグレイスケールダイレーション手法に替え、多値画像において濃度値の高い範囲を周囲に拡張し得る他の多値画像用の拡張処理手法を用いることも可能である。   In the above embodiment, the gray scale dilation method is used as an expansion processing method for multi-valued images, and a 3 × 3 square is exemplified as a structuring element at that time. It is possible to change the method of replacing the density value in the sheath or the structuring element, and in addition to the gray scale dilation method described in Non-Patent Document 1, the density value is high in the multi-value image. It is also possible to use other expansion processing techniques for multi-valued images that can expand the range to the periphery.

多値画像用の拡張処理を用いた減算処理手順を模式的に示す図The figure which shows typically the subtraction process procedure using the extended process for multi-valued images グレイスケールダイレーション手法を模式的に示す図Diagram showing gray scale dilation technique 脳血管と頭蓋骨とが近接している部分を例示する図The figure which illustrates the part where the cerebral blood vessel and the skull are close 脳血管が頭蓋骨を貫通する部分の処理手順を模式的に示す図The figure which shows the processing procedure of the part where the cerebral blood vessel penetrates the skull typically 脳血管が頭蓋骨に沿う部分の処理手順を模式的に示す図The figure which shows typically the processing procedure of the part where a cerebral blood vessel follows a skull. 脳血管領域の最終的な確定処理前後の画像を示す図The figure which shows the image before and after the final decision processing of the cerebral blood vessel region 従来手法と本発明方法とにより得られる差分画像の比較図Comparison diagram of difference images obtained by the conventional method and the method of the present invention

符号の説明Explanation of symbols

11,21,31 低造影画像
12,22,32 高造影画像
13,23,33 拡張処理後画像
14,24,34 差分画像
15,25,35 骨部領域
16,26,36 非骨領域
17,27,37 脳血管領域
27a,37a 抽出された脳血管領域
27b,37b 縮減された脳血管領域
28 周辺拡大画像
41 第1画像
42 構造化要素(カーネル,オペレータ)
43 第2画像
11, 21, 31 Low contrast image 12, 22, 32 High contrast image 13, 23, 33 Extended image 14, 24, 34 Difference image 15, 25, 35 Bone region 16, 26, 36 Non-bone region 17, 27, 37 Cerebrovascular region 27a, 37a Extracted cerebrovascular region 27b, 37b Reduced cerebrovascular region 28 Peripheral enlarged image 41 First image 42 Structuring element (kernel, operator)
43 Second image

Claims (4)

造影効果が相対的に低い時間帯に撮像された低造影画像と、造影効果が相対的に高い時間帯に撮像された高造影画像とに基づいて、被造影領域を抽出する医療用画像の被造影領域抽出方法において、
多値画像用の拡張処理を前記低造影画像に施して、該低造影画像における骨部領域を拡張し、
前記高造影画像と前記拡張処理後の前記低造影画像との減算処理によって差分画像を求めることにより、前記被造影領域を抽出することを特徴とする医療用画像の被造影領域抽出方法。
Based on a low-contrast image captured in a time zone with a relatively low contrast effect and a high-contrast image captured in a time zone with a relatively high contrast effect, In the contrast area extraction method,
The low contrast image is subjected to expansion processing for a multivalued image, and a bone region in the low contrast image is expanded,
A method for extracting a contrast region of a medical image, wherein the contrast region is extracted by obtaining a difference image by a subtraction process between the high contrast image and the low contrast image after the expansion process.
前記差分画像において抽出された前記被造影領域のうち縮減された必要部分に対して、領域の復元処理を行なうことを特徴とする請求項1記載の医療用画像の被造影領域抽出方法。   The method for extracting a contrast region of a medical image according to claim 1, wherein a region restoration process is performed on a reduced necessary portion of the contrast region extracted in the difference image. 前記被造影領域が、脳血管の領域であることを特徴とする請求項1または2記載の医療用画像の被造影領域抽出方法。   3. The method for extracting a contrast-enhanced region of a medical image according to claim 1, wherein the region to be contrasted is a cerebral blood vessel region. 前記拡張処理が、グレイスケールダイレーション手法を用いるものであることを特徴とする請求項1〜3までのいずれか1項記載の医療用画像の被造影領域抽出方法。   The method for extracting a contrast-enhanced region of a medical image according to any one of claims 1 to 3, wherein the expansion processing uses a gray scale dilation method.
JP2004052302A 2004-02-19 2004-02-26 Method for extracting imaged area of image for medical purpose Withdrawn JP2005261440A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2004052302A JP2005261440A (en) 2004-02-19 2004-02-26 Method for extracting imaged area of image for medical purpose

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2004043601 2004-02-19
JP2004052302A JP2005261440A (en) 2004-02-19 2004-02-26 Method for extracting imaged area of image for medical purpose

Publications (1)

Publication Number Publication Date
JP2005261440A true JP2005261440A (en) 2005-09-29

Family

ID=35086591

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2004052302A Withdrawn JP2005261440A (en) 2004-02-19 2004-02-26 Method for extracting imaged area of image for medical purpose

Country Status (1)

Country Link
JP (1) JP2005261440A (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008043564A (en) * 2006-08-17 2008-02-28 Fujifilm Corp Domain definer
JP2008073301A (en) * 2006-09-22 2008-04-03 Toshiba Corp Medical imaging diagnostic apparatus and medical image processor
JP2010011980A (en) * 2008-07-02 2010-01-21 Toshiba Corp Image processor and image diagnostic apparatus
JP2012020174A (en) * 2011-10-17 2012-02-02 Toshiba Corp Medical diagnostic imaging apparatus and medical image processor
JP2014008414A (en) * 2012-07-02 2014-01-20 Toshiba Corp Motion correction apparatus and method
JP2014237005A (en) * 2007-04-20 2014-12-18 メディシム・ナムローゼ・フエンノートシャップ Method for deriving shape information
WO2015015777A1 (en) * 2013-07-31 2015-02-05 富士フイルム株式会社 Medical image display control device, method of operation for same, and medical image display control program
JP2018166657A (en) * 2017-03-29 2018-11-01 ザイオソフト株式会社 Medical image processing device, medical image processing method, and medical image processing program
CN109350089A (en) * 2018-10-25 2019-02-19 杭州医学院 A kind of thyroid gland areas automatic division method based on ct images
JP2020039507A (en) * 2018-09-07 2020-03-19 株式会社日立製作所 Magnetic resonance imaging device, image processing device, and image processing method

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008043564A (en) * 2006-08-17 2008-02-28 Fujifilm Corp Domain definer
JP2008073301A (en) * 2006-09-22 2008-04-03 Toshiba Corp Medical imaging diagnostic apparatus and medical image processor
US9439608B2 (en) 2007-04-20 2016-09-13 Medicim Nv Method for deriving shape information
JP2014237005A (en) * 2007-04-20 2014-12-18 メディシム・ナムローゼ・フエンノートシャップ Method for deriving shape information
JP2010011980A (en) * 2008-07-02 2010-01-21 Toshiba Corp Image processor and image diagnostic apparatus
JP2012020174A (en) * 2011-10-17 2012-02-02 Toshiba Corp Medical diagnostic imaging apparatus and medical image processor
JP2014008414A (en) * 2012-07-02 2014-01-20 Toshiba Corp Motion correction apparatus and method
JP2015029518A (en) * 2013-07-31 2015-02-16 富士フイルム株式会社 Medical image display control apparatus, method of operating the same, and medical image display control program
WO2015015777A1 (en) * 2013-07-31 2015-02-05 富士フイルム株式会社 Medical image display control device, method of operation for same, and medical image display control program
US10062184B2 (en) 2013-07-31 2018-08-28 Fujifilm Corporation Medical image display control device, method of operation for same, and medical image display control program
JP2018166657A (en) * 2017-03-29 2018-11-01 ザイオソフト株式会社 Medical image processing device, medical image processing method, and medical image processing program
JP2020039507A (en) * 2018-09-07 2020-03-19 株式会社日立製作所 Magnetic resonance imaging device, image processing device, and image processing method
JP7125312B2 (en) 2018-09-07 2022-08-24 富士フイルムヘルスケア株式会社 MAGNETIC RESONANCE IMAGING APPARATUS, IMAGE PROCESSING APPARATUS, AND IMAGE PROCESSING METHOD
US11701024B2 (en) 2018-09-07 2023-07-18 Fujifilm Healthcare Corporation Magnetic resonance imaging apparatus, image processor, and image processing method
CN109350089A (en) * 2018-10-25 2019-02-19 杭州医学院 A kind of thyroid gland areas automatic division method based on ct images
CN109350089B (en) * 2018-10-25 2022-03-08 杭州医学院 Automatic thyroid region segmentation method based on CT (computed tomography) image

Similar Documents

Publication Publication Date Title
US20070016108A1 (en) Method for 3D visualization of vascular inserts in the human body using the C-arm
US5832134A (en) Data visualization enhancement through removal of dominating structures
JP3486615B2 (en) Region extraction method for medical images
JP6448972B2 (en) Medical image processing apparatus and medical image processing method
Patil et al. Preprocessing to be considered for MR and CT images containing tumors
JP2013141603A (en) Medical image processing apparatus, medical image processing method and medical image processing program
US7330573B2 (en) Visualization of medical images
JP2011139891A (en) Medical image processing apparatus and medical image processing method
CN109350089B (en) Automatic thyroid region segmentation method based on CT (computed tomography) image
JP2003153082A (en) Image aligner and image processor
US20130057569A1 (en) 2d-2d fusion for interventional guidance in trans-catheter aortic valve implantation
CN105761217A (en) Image reconstruction method and device
JP2005261440A (en) Method for extracting imaged area of image for medical purpose
CN110876627B (en) X-ray imaging apparatus and X-ray image processing method
US20220138936A1 (en) Systems and methods for calcium-free computed tomography angiography
JP2010213760A (en) Image processing device, and image processing program
Kroon et al. Coherence filtering to enhance the mandibular canal in cone-beam CT data
US12039718B2 (en) Inference apparatus, medical apparatus, and program
Sepehrian et al. Individual teeth segmentation in CBCT and MSCT dental images using watershed
JPH0759763A (en) Three-dimensional object measuring system
JP2015177972A (en) Image processor and x-ray diagnostic apparatus
van Andel et al. Interobserver variability in the detection of cerebral venous thrombosis using CT venography with matched mask bone elimination
Hiroshima et al. Digital subtraction angiogram registration method with local distortion vectors to decrease motion artifact
CN111862312A (en) Brain blood vessel display device and method
Fukuda et al. Three-dimensional visualization of aortic arch morphology via rearrangement of non-enhanced thoracic CT images

Legal Events

Date Code Title Description
A300 Application deemed to be withdrawn because no request for examination was validly filed

Free format text: JAPANESE INTERMEDIATE CODE: A300

Effective date: 20070501