JP6877534B2 - Image processing systems and methods for interactive contouring of 3D medical data - Google Patents

Image processing systems and methods for interactive contouring of 3D medical data Download PDF

Info

Publication number
JP6877534B2
JP6877534B2 JP2019512197A JP2019512197A JP6877534B2 JP 6877534 B2 JP6877534 B2 JP 6877534B2 JP 2019512197 A JP2019512197 A JP 2019512197A JP 2019512197 A JP2019512197 A JP 2019512197A JP 6877534 B2 JP6877534 B2 JP 6877534B2
Authority
JP
Japan
Prior art keywords
source data
edge
interest
subset
processing unit
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.)
Active
Application number
JP2019512197A
Other languages
Japanese (ja)
Other versions
JP2019535336A (en
Inventor
ウェイストランド，オラ
Original Assignee
レイサーチ ラボラトリーズ，エービー
レイサーチ ラボラトリーズ，エービー
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by レイサーチ ラボラトリーズ，エービー, レイサーチ ラボラトリーズ，エービー filed Critical レイサーチ ラボラトリーズ，エービー
Publication of JP2019535336A publication Critical patent/JP2019535336A/en
Application granted granted Critical
Publication of JP6877534B2 publication Critical patent/JP6877534B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • G06T7/12—Edge-based segmentation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • G06T7/149—Segmentation; Edge detection involving deformable models, e.g. active contour models
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • G06T7/174—Segmentation; Edge detection involving the use of two or more images
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20092—Interactive image processing based on input by user
    • G06T2207/20096—Interactive definition of curve of interest
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Software Systems (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)
  • Processing Or Creating Images (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Description

本発明は、一般に、医用画像データの情報内容を強化するためのソリューションに関する。より具体的には、本発明は、請求項１のプリアンブルに係る画像処理システム及び対応する方法に関する。本発明はまた、コンピュータプログラム及びプロセッサ可読媒体に関する。 The present invention generally relates to a solution for enhancing the information content of medical image data. More specifically, the present invention relates to the image processing system according to the preamble of claim 1 and the corresponding method. The present invention also relates to computer programs and processor readable media.

どのボクセルが特定の解剖学的構造体を表すかを定義するプロセス、すなわち、いわゆる臓器の描写は、放射線治療計画作成のなかでも最も長々と時間がかかることの１つである。このプロセスは、普通は、簡単な描画ツールを用いて二次元スライスにおいて手作業で輪郭描出することに関係し、計画作成のために用いられる高解像度の三次元データセットにおける関心あるすべての構造体を描写するのに数時間を要する場合がある。 The process of defining which voxels represent a particular anatomical structure, the so-called organ depiction, is one of the longest-running radiotherapy planning processes. This process usually involves manual delineation of 2D slices using simple drawing tools and all structures of interest in high resolution 3D datasets used for planning. It may take several hours to depict.

Ｐｅｋａｒ，Ｖ．，ｅｔ ａｌ．，“Ａｕｔｏｍａｔｅｄ Ｍｏｄｅｌ−Ｂａｓｅｄ Ｏｒｇａｎ Ｄｅｌｉｎｅａｔｉｏｎ ｆｏｒ Ｒａｄｉｏｔｈｅｒａｐｙ Ｐｌａｎｎｉｎｇ ｉｎ Ｐｒｏｓｔａｔｉｃ Ｒｅｇｉｏｎ”， Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｊｏｕｒｎａｌ ｏｆ Ｒａｄｉａｔｉｏｎ Ｏｎｃｏｌｏｇｙ−Ｂｉｏｌｏｇｙ−Ｐｈｙｓｉｃｓ，Ｖｏｌ．６０，Ｎｏ．３，ｐｐ９７３−９８０，２００４は、関心ある解剖学的構造体の境界に３Ｄの変形可能な表面モデルを適応させるための方法を開示する。この適応は、或る画像特徴へ寄せることにより生じるモデルの変形とモデルの形状完全性との兼ね合いに基づいている。自動モデル適応に失敗し得る問題のある領域は、対話型ツールにより補正することができる。 Pekar, V.I. , Et al. , "Automated Model-Based Organization for Radiation therapy Planning in Prostate Region", International Journal of Radiation Oncology-Biology-Physics. 60, No. 3, pp973-980, 2004 discloses a method for adapting a 3D deformable surface model to the boundaries of an anatomical structure of interest. This adaptation is based on the trade-off between model deformation and model shape integrity caused by shifting to certain image features. Problem areas where automatic model adaptation can fail can be corrected with interactive tools.

ＵＳ２０１１／０２６８３３０は、医用画像の組を輪郭描出するためのシステム及び方法を説明している。例示的なシステムは、画像データベース、画像変形エンジン、及び輪郭変換エンジンを含み得る。画像データベースは、医用画像の組を記憶するのに用いられ得る。画像変形エンジンは、画像データベースにおける医用画像の組からソース画像及び対象画像を受信するように構成され、ソース画像及び対象画像に変形アルゴリズムを用いてソース画像から対象画像への１つ以上のオブジェクト間の変化を示す変形フィールドデータを生成するようにさらに構成され得る。輪郭変換エンジンは、ソース画像内の１つ以上のオブジェクトを識別するソース輪郭データを受信するように構成され、変形フィールドデータ及びソース輪郭データを用いて対象画像内の１つ以上のオブジェクトを識別する自動対象輪郭データを生成するようにさらに構成され得る。画像変形エンジン及び輪郭変換エンジンは、１つ以上のメモリデバイスに記憶され、１つ以上のプロセッサにより実行可能な、ソフトウェア命令を備え得る。 US2011 / 0268330 describes a system and method for outlining a set of medical images. An exemplary system may include an image database, an image transformation engine, and a contour transformation engine. The image database can be used to store a set of medical images. The image transformation engine is configured to receive a source image and a target image from a set of medical images in an image database, and between one or more objects from the source image to the target image using a transformation algorithm on the source image and the target image. It may be further configured to generate transformation field data indicating changes in. The contour transformation engine is configured to receive source contour data that identifies one or more objects in the source image and uses deformed field data and source contour data to identify one or more objects in the target image. It may be further configured to generate automatic target contour data. The image transformation engine and contour transformation engine may include software instructions stored in one or more memory devices and executed by one or more processors.

上記のソリューションは、形状の或る範囲内の所定のタイプの臓器、すなわち、関心ある構造体に関して満足な最終結果をもたらすことができる場合がある。しかしながら、それに基づいて描写が行われる記憶された輪郭のライブラリにその形状があてはまらない関心ある構造体全体の描写は、問題がある及び／又は時間がかかる場合がある。 The above solution may be able to give satisfactory end results for certain types of organs within a range of shapes, i.e. structures of interest. However, the depiction of the entire structure of interest whose shape does not apply to the library of stored contours on which the depiction is made can be problematic and / or time consuming.

したがって、本発明の目的は、一方では、比較的軽度の人手の介入を必要とし、他方では、三次元画像データセットにおける改良された幾何学的補間をもたらすソリューションを提供することである。 Therefore, an object of the present invention is to provide a solution that, on the one hand, requires relatively light manual intervention and, on the other hand, results in improved geometric interpolation in a three-dimensional image dataset.

本発明の一態様によれば、目的は、命令の第１のセットを受信した後で処理ユニットが以下のことを行うように構成される、最初に説明した構成により達成される。ソースデータにおける第２の軸方向の位置で関心ある解剖学的構造体を通る第２の画像スライスにおいて構成されるソースデータの第２の二次元グラフィック表現を定義するソースデータの第２のサブセットを選択する。第２の軸方向の位置は、第１の軸方向の位置とは異なる。次いで、処理ユニットは、グラフィカルディスプレイ上に提示するために第２のインターフェースにソースデータの第２のサブセットを出力させるように構成された制御命令を生成するように構成される。第３のインターフェースを介して、処理ユニットは、第２の画像スライスにおける関心ある解剖学的構造体の第２のエッジを識別する命令の第２のセットを受信するように構成される。その後、第１及び第２のエッジとソースデータに基づいて、処理ユニットは、関心ある解剖学的構造体の境界面の近似を表す三次元シェルを計算するように構成される。 According to one aspect of the invention, an object is achieved by the configuration described first, in which the processing unit is configured to do the following after receiving a first set of instructions. A second subset of the source data that defines a second two-dimensional graphic representation of the source data that consists of a second image slice that passes through the anatomical structure of interest at the second axial position in the source data. select. The second axial position is different from the first axial position. The processing unit is then configured to generate control instructions configured to cause the second interface to output a second subset of the source data for presentation on the graphical display. Through the third interface, the processing unit is configured to receive a second set of instructions that identify the second edge of the anatomical structure of interest in the second image slice. Based on the first and second edges and source data, the processing unit is then configured to calculate a 3D shell that represents an approximation of the interface of the anatomical structure of interest.

この構成は、第１及び第２のエッジとソースデータを併せて用いることにより、比較的少ない離散ステップで関心ある解剖学的構造体の高品質の描写に達することが可能となるので有利である。他の利点は、輪郭ライブラリが必要とされないことである。さらに、ソースデータの異なるサブセットは、互いに特に類似している必要はない。言い換えれば、画像スライスは、比較的大きい距離だけ分離することができる。もちろん、これはソリューションをユーザの見地から融通性の高いものにする。 This configuration is advantageous because the combined use of the first and second edges and the source data allows a high quality depiction of the anatomical structure of interest to be achieved with relatively few discrete steps. .. Another advantage is that no contour library is needed. Moreover, the different subsets of source data need not be particularly similar to each other. In other words, image slices can be separated by a relatively large distance. Of course, this makes the solution flexible from the user's point of view.

本発明のこの態様の一実施形態によれば、処理ユニットは、データポイントの第１及び第２のセット間を排他的に補間することによりソースデータにおける第１及び第２の軸方向の位置間の三次元シェルの最初の推定を計算するように構成され、データポイントの第１のセットは、第１のエッジのアウトラインを記述し、データポイントの第２のセットは、第２のエッジのアウトラインを記述する。したがって、以下の近似プロセスのための信頼できる基礎が確立される。 According to one embodiment of this aspect of the invention, the processing unit is between the first and second axial positions in the source data by exclusively interpolating between the first and second sets of data points. The first set of data points describes the outline of the first edge, and the second set of data points outlines the second edge. To describe. Therefore, a reliable basis for the following approximation process is established.

本発明のこの態様の別の実施形態によれば、第１の画像スライスはソースデータにおいて第１の配向を有し、第２の画像スライスは第１の配向と平行に配向される。処理ユニットはさらに、表面再構築アルゴリズム、例えば、ユーザにより描かれた輪郭、例えば凸包又は同様のものを近似する三次元表面を計算するための形状に基づく補間を採用するアルゴリズム又は任意の他のアルゴリズムを用いることにより三次元シェルの最初の推定を計算するように構成される。これにより、近似プロセスは、ほとんどのタイプの解剖学的構造体に関して容認できる結果へ非常に迅速に収束する。 According to another embodiment of this aspect of the invention, the first image slice has a first orientation in the source data and the second image slice is oriented parallel to the first orientation. The processing unit also employs a surface reconstruction algorithm, eg, a shape-based interpolation for computing a three-dimensional surface that approximates a user-drawn contour, eg, a convex hull or the like, or any other algorithm. It is configured to calculate the first estimate of the 3D shell by using an algorithm. This allows the approximation process to converge very quickly to acceptable results for most types of anatomical structures.

したがって、好ましくは、処理ユニットは、三次元シェルの最初の推定とソースデータに基づいて、三次元シェルの第１の後続する推定を計算するように構成される。その後、三次元シェルの先行する推定と第１及び第２のエッジとソースデータに基づいて、任意の後続する推定が計算される。これにより、三次元シェルの推定は、手書きの輪郭と画像データとの両方に幾何学的に緊密に保たれ、これは高品質の最終結果を保証する。 Therefore, preferably, the processing unit is configured to calculate the first subsequent estimate of the 3D shell based on the initial estimate of the 3D shell and the source data. Any subsequent estimates are then calculated based on the preceding estimates of the 3D shell and the first and second edges and source data. This keeps the 3D shell estimation geometrically close to both the handwritten contours and the image data, which guarantees a high quality end result.

本発明のこの態様のさらに別の実施形態によれば、処理ユニットはさらに、第１のエッジの近傍のソースデータの第１のサブセットの第１のデータ部分を解析し、第２のエッジの近傍のソースデータの第２のサブセットの第２のデータ部分を解析するように構成される。これらの解析に基づいて、次に、処理ユニットは、関心ある解剖学的構造体の推定される外周を定義する画像特徴の集合を導出し、例えば調和関数に基づく補間を用いて、画像特徴の集合にさらに基づいて三次元シェルの後続する推定を計算するように構成され、画像特徴の集合は補間プロセスの境界条件を表す。その結果、画像特徴の特性は、関心ある構造体の周りの表面にわたってかなり変化し得るが、その外周は、正確に推定することができる。例えば、この手法は、その骨格が構造体のいくつかの部分で非常に強いエッジを有し、他の部分で非常に弱いエッジを有し得る（そこでは代わりに周囲の構造体が非常に強いエッジを有し得る）大腿骨頭の描写を比較的問題がないものにする。 According to yet another embodiment of this aspect of the invention, the processing unit further analyzes the first data portion of the first subset of source data in the vicinity of the first edge and in the vicinity of the second edge. It is configured to analyze the second data portion of the second subset of the source data of. Based on these analyses, the processing unit then derives a set of image features that define the estimated perimeter of the anatomical structure of interest, using, for example, harmonic function-based interpolation of the image features. It is configured to calculate subsequent estimates of the three-dimensional shell based on the set, and the set of image features represents the boundary conditions of the interpolation process. As a result, the properties of the image features can vary considerably over the surface around the structure of interest, but its perimeter can be estimated accurately. For example, in this technique, the skeleton can have very strong edges in some parts of the structure and very weak edges in other parts (where the surrounding structure is very strong instead). Makes the depiction of the femoral head (which may have edges) relatively benign.

本発明のこの態様のまたさらに別の実施形態によれば、処理ユニットは、三次元シェルを推定される外周に近似する三角形メッシュ構造の第１反復を生成するように構成される。三角形メッシュ構造は、それぞれメッシュ構造のいくつかの三角形が出会う頂点の組を含む。好ましくは、処理ユニットはさらに、第１反復に加えて、三角形メッシュ構造の少なくとも１つのさらなる反復を生成するように構成され、三角形メッシュ構造の少なくとも１つのさらなる反復は、推定される外周と三角形メッシュ構造の頂点との間の全距離測度に関する第１反復に対する改良である、推定される外周の近似を表す。したがって、関心ある解剖学的構造体の外周を、三角形メッシュ構造を通じて正確に表現することができ、これは効率的なグラフィカル処理及び視覚化を可能にする。 According to yet another embodiment of this aspect of the invention, the processing unit is configured to generate a first iteration of a triangular mesh structure that approximates the estimated perimeter of the 3D shell. A triangular mesh structure contains a set of vertices where several triangles of the mesh structure meet. Preferably, the processing unit is further configured to generate at least one additional iteration of the triangular mesh structure in addition to the first iteration, with at least one additional iteration of the triangular mesh structure being the estimated perimeter and triangular mesh. Represents an approximation of the estimated perimeter, which is an improvement over the first iteration for the total distance measurement to the vertices of the structure. Therefore, the perimeter of the anatomical structure of interest can be accurately represented through the triangular mesh structure, which enables efficient graphical processing and visualization.

本発明のこの態様の１つのさらなる実施形態によれば、処理ユニットは、第１及び第２のサブセットに加えて、ソースデータの少なくとも１つのさらなるサブセットを選択するように構成される。少なくとも１つのさらなるサブセットは、ソースデータにおける少なくとも１つのさらなる軸方向の位置で関心ある解剖学的構造体を通る少なくとも１つのさらなる画像スライスにおいて構成されるソースデータの少なくとも１つのさらなる二次元グラフィック表現を定義する。少なくとも１つのさらなる軸方向の位置は、第１及び第２の軸方向の位置とは異なる。少なくとも１つのさらなるサブセットのそれぞれに関して、処理ユニットは、グラフィカルディスプレイ上に提示するために第２のインターフェースにソースデータの少なくとも１つのさらなるサブセットを出力させるように構成された制御命令を生成するように構成される。次いで、第３のインターフェースを介して受信した命令の少なくとも１つのさらなるセットに応答して、処理ユニットは、少なくとも１つのさらなる画像スライスにおける関心ある解剖学的構造体の少なくとも１つのさらなるエッジを識別するように構成される。その後、第１、第２、及び少なくとも１つのさらなるエッジとソースデータに基づいて、処理ユニットは、関心ある解剖学的構造体の表面推定を表す三次元シェルの改良バージョンを計算するように構成される。これは、ユーザが１つ以上の手書きの輪郭を追加することにより自動描写プロセスを都合よく支援し、スピードアップできることを意味する。 According to one further embodiment of this aspect of the invention, the processing unit is configured to select at least one additional subset of source data in addition to the first and second subsets. At least one additional subset provides at least one additional two-dimensional graphic representation of the source data composed of at least one additional image slice passing through the anatomical structure of interest at at least one additional axial position in the source data. Define. The at least one additional axial position is different from the first and second axial positions. For each of the at least one additional subset, the processing unit is configured to generate control instructions configured to cause the second interface to output at least one additional subset of source data for presentation on the graphical display. Will be done. The processing unit then identifies at least one additional edge of the anatomical structure of interest in at least one additional image slice in response to at least one additional set of instructions received via the third interface. It is configured as follows. Based on the first, second, and at least one additional edge and source data, the processing unit is then configured to calculate an improved version of the 3D shell that represents the surface estimation of the anatomical structure of interest. To. This means that the user can conveniently assist and speed up the automatic drawing process by adding one or more handwritten contours.

本発明のこの態様のさらに別の実施形態によれば、第３のインターフェースは、命令の補足的なセットを処理ユニットに送るように構成され、命令の補足的なセットは、前記画像スライスのうちの少なくとも１つにおける関心ある解剖学的構造体の調整されたエッジを識別する。そのうえ、処理ユニットは、調整されたエッジにさらに基づいて三次元シェルを計算するように構成される。したがって、ユーザは、関心ある解剖学的構造体のエッジのアウトラインを直観的な様態で、手動で変更することにより描写プロセスを微調整することもできる。 According to yet another embodiment of this aspect of the invention, the third interface is configured to send a complementary set of instructions to the processing unit, the complementary set of instructions being among the image slices. Identify the adjusted edges of the anatomical structure of interest in at least one of the. Moreover, the processing unit is configured to calculate the 3D shell based further on the adjusted edges. Therefore, the user can also fine-tune the depiction process by manually changing the outline of the edges of the anatomical structure of interest in an intuitive manner.

本発明のこの態様のまたさらに別の実施形態によれば、第３のインターフェースは、例えばコンピュータマウス又は同様のポインティングデバイスを介して生成された、ユーザコマンドの形態の命令のセットを受信することに特化して構成される。 According to yet another embodiment of this aspect of the invention, the third interface is to receive a set of instructions in the form of user commands generated, for example, via a computer mouse or similar pointing device. Specially configured.

本発明の別の態様によれば、この目的は、命令の第１のセットを受信した後で、ソースデータの第２の二次元グラフィック表現を定義するソースデータの第２のサブセットが選択される、最初に説明した方法により達成される。第２のサブセットは、ソースデータにおける第２の軸方向の位置で関心ある解剖学的構造体を通る第２の画像スライスにおいて構成され、第２の軸方向の位置は第１の軸方向の位置とは異なる。第２のインターフェースを介して、ソースデータの第２のサブセットが、グラフィカルディスプレイ上に提示するために出力される。その後、第３のインターフェースを介して受信した命令の第２のセットに応答して、第２の画像スライスにおける関心ある解剖学的構造体の第２のエッジが識別される。最後に、第１及び第２のエッジとソースデータに基づいて、三次元シェルが計算される。三次元シェルは、関心ある解剖学的構造体の境界面の近似を表す。この方法、並びに、その好ましい実施形態の利点は、提案されるシステムに関する上記の説明から明白である。 According to another aspect of the invention, the object is to select a second subset of the source data that defines a second two-dimensional graphic representation of the source data after receiving the first set of instructions. , Achieved by the method described first. The second subset consists of a second image slice through the anatomical structure of interest at the second axial position in the source data, where the second axial position is the first axial position. Is different. Through the second interface, a second subset of the source data is output for presentation on a graphical display. The second edge of the anatomical structure of interest in the second image slice is then identified in response to a second set of instructions received via the third interface. Finally, the 3D shell is calculated based on the first and second edges and the source data. The 3D shell represents an approximation of the interface of the anatomical structure of interest. The advantages of this method, as well as its preferred embodiments, are evident from the above description of the proposed system.

本発明のさらなる態様によれば、この目的は、少なくとも１つのプロセッサのメモリにロード可能であり、且つ、少なくとも１つのプロセッサ上で実行されるときに上記で提案した方法を実施するように適応されたソフトウェアを含む、コンピュータプログラムにより達成される。 According to a further aspect of the invention, this object is loadable to the memory of at least one processor and is adapted to carry out the method proposed above when running on at least one processor. Achieved by computer programs, including software.

本発明の別の態様によれば、この目的は、少なくとも１つのプロセッサにロードされるときに上記で提案した方法を行うように少なくとも１つのプロセッサを制御するためのプログラムが記録されている、プロセッサ可読媒体により達成される。 According to another aspect of the invention, the object is a processor in which a program is recorded to control at least one processor to perform the method proposed above when loaded into at least one processor. Achieved by readable media.

本発明のさらなる利点、有益な特徴、及び用途は、以下の説明及び従属請求項から明白であろう。 Further advantages, beneficial features, and uses of the present invention will be apparent from the following description and dependent claims.

ここで、添付図を参照しながら、例として開示される好ましい実施形態によって本発明をより詳しく説明する。 Here, the present invention will be described in more detail by the preferred embodiments disclosed as examples with reference to the accompanying drawings.

本発明の一実施形態に係るシステムの概要を示す図である。It is a figure which shows the outline of the system which concerns on one Embodiment of this invention. 本発明の実施形態に係る関心ある解剖学的構造体のエッジがソースデータの異なる画像スライスにおいてどのように識別され得るかを例示する図である。FIG. 5 illustrates how the edges of an anatomical structure of interest according to an embodiment of the invention can be identified in different image slices of source data. 本発明の一実施形態に係る生成された三角形メッシュの例を示す図である。It is a figure which shows the example of the generated triangular mesh which concerns on one Embodiment of this invention. 提案される三角形メッシュに関連して距離測度を例示する図である。It is a figure which illustrates the distance measure in relation to the proposed triangular mesh. 本発明の一実施形態に係る描かれたエッジの近傍のソースデータがどのように解析されるかを例示する図である。It is a figure which illustrates how the source data in the vicinity of the drawn edge which concerns on one Embodiment of this invention is analyzed. 本発明に係る一般的方法を流れ図で例示する図である。It is a figure which illustrates the general method which concerns on this invention in the flow chart. 提案される方法の一実施形態を例示する流れ図の一部を示す図である。It is a figure which shows a part of the flow chart which illustrates one Embodiment of the proposed method.

図１は、本発明の一実施形態に係る三次元医用画像データを輪郭描出するための画像処理システム１００の概要を示す。システム１００は、処理ユニット１４０と、第１、第２、及び第３のインターフェース１１０、１２０、及び１３０を含む。好ましくは、処理ユニット１４０上で実行されているときに後述する手順を実行するためのソフトウェアを記憶する、メモリ１５０も含まれる。提示の目的で、図１は、インターフェース１１０、１２０、及び１３０を別個のエンティティとして例示する。しかしながら、実際的な実装では、インターフェースのうちの２つ以上が共通のユニットへ一体化されてよい。 FIG. 1 shows an outline of an image processing system 100 for delineating three-dimensional medical image data according to an embodiment of the present invention. The system 100 includes a processing unit 140 and first, second, and third interfaces 110, 120, and 130. Preferably, a memory 150 is also included, which stores software for performing the procedure described below when running on the processing unit 140. For presentation purposes, FIG. 1 illustrates interfaces 110, 120, and 130 as separate entities. However, in a practical implementation, two or more of the interfaces may be integrated into a common unit.

第１のインターフェース１１０は、関心ある解剖学的構造体を表す三次元医用画像データの形態のソースデータＳＤに処理ユニット１４０がアクセスすることを可能にするように構成される。ソースデータＳＤは、関心ある解剖学的構造体に隣接する組織も含む。提案されるシステム１００の目的は、例えば、関心ある解剖学的構造体の放射線治療を計画するために、この構造体を隣接する組織から区別することができるように、関心ある解剖学的構造体を輪郭描出する、すなわち描写することである。 The first interface 110 is configured to allow the processing unit 140 to access the source data SD in the form of three-dimensional medical image data representing the anatomical structure of interest. Source data SD also includes tissue adjacent to the anatomical structure of interest. The purpose of the proposed system 100 is to allow the anatomical structure of interest to be distinguished from adjacent tissues, for example, in order to plan radiotherapy for the anatomical structure of interest. Is to outline, that is, to depict.

ここで図２も参照すると、処理ユニット１４０は、ソースデータＳＤにおける第１の軸方向の位置ｚ１で関心ある解剖学的構造体を通る第１の画像スライスにおいて構成されるソースデータＳＤの第１の二次元グラフィック表現を定義するソースデータＳＤの第１のサブセットｓｄ１を選択するように構成される。 Also referring to FIG. 2, the processing unit 140 is the first of the source data SD composed of the first image slice passing through the anatomical structure of interest at the first axial position z1 in the source data SD. It is configured to select the first subset sd1 of the source data SD that defines the two-dimensional graphic representation of.

第２のインターフェース１２０は、グラフィカルディスプレイ１７０上に提示するためにソースデータの第１のサブセットｓｄ１を出力するように構成される。ソースデータの第１のサブセットｓｄ１は、処理ユニット１４０からの制御命令に応答して出力される。 The second interface 120 is configured to output a first subset sd1 of source data for presentation on the graphical display 170. The first subset sd1 of the source data is output in response to a control instruction from the processing unit 140.

第３のインターフェース１３０は、命令［Ｉ］ｅｄｇｅの第１のセットを処理ユニット１４０に送るように構成される。命令［Ｉ］ｅｄｇｅの第１のセットは、第１の画像スライスにおける関心ある解剖学的構造体の第１のエッジＥ１を識別し、好ましくは、キーボード及び／又はコンピュータマウス、又は同様のポインティングデバイスを介して入力されるユーザコマンドｕｃｍｄ１及び／又はｕｃｍｄ２に応答して生成される。 The third interface 130 is configured to send a first set of instructions [I] edges to the processing unit 140. The first set of instructions [I] edge identifies the first edge E1 of the anatomical structure of interest in the first image slice, preferably a keyboard and / or computer mouse, or similar pointing device. It is generated in response to the user commands ucmd1 and / or ucmd2 input via.

図２は、第１の画像スライスによって表されるソースデータの第１のサブセットｓｄ１でのソースデータＳＤにおける関心ある解剖学的構造体の形状に第１のエッジＥ１がどのように従うかを例示する。 FIG. 2 illustrates how the first edge E1 follows the shape of the anatomical structure of interest in the source data SD at the first subset sd1 of the source data represented by the first image slice. ..

命令［Ｉ］ｅｄｇｅの第１のセットを受信した後で、処理ユニット１４０は、ソースデータＳＤの第２の二次元グラフィック表現を定義するソースデータＳＤの第２のサブセットｓｄ２を選択するように構成される。第２のサブセットｓｄ２は、ソースデータＳＤにおける第２の軸方向の位置ｚ２で関心ある解剖学的構造体を通る第２の画像スライスにおいて構成される。第２の軸方向の位置ｚ２は、第１の軸方向の位置ｚ１とは異なる。言い換えれば、第２のサブセットｓｄ２は、ソースデータＳＤからのデータの別個の二次元スライスを含み、このデータのスライスは、第１のサブセットｓｄ１によって表されるスライスと平行であってよく、又は平行でなくてもよい。しかしながら、後述するように、第１の画像スライスと第２の画像スライスが実際に互いに平行であれば有利である。 After receiving the first set of instructions [I] edge , the processing unit 140 is configured to select the second subset sd2 of the source data SD that defines the second two-dimensional graphic representation of the source data SD. Will be done. The second subset sd2 is composed of a second image slice through the anatomical structure of interest at the second axial position z2 in the source data SD. The second axial position z2 is different from the first axial position z1. In other words, the second subset sd2 contains a separate two-dimensional slice of data from the source data SD, and the slice of this data may be parallel or parallel to the slice represented by the first subset sd1. It does not have to be. However, as will be described later, it is advantageous if the first image slice and the second image slice are actually parallel to each other.

処理ユニットは、グラフィカルディスプレイ１７０上に提示するために、すなわち、ユーザがこれを視覚的に観察し、第２の画像スライスにおける関心ある解剖学的構造体の第２のエッジＥ２を記述するためのコマンドを生成することができるように、第２のインターフェース１２０にソースデータＳＤの第２のサブセットｓｄ２を出力させるように適応された制御命令を生成するように構成される。 The processing unit is for presentation on the graphical display 170, i.e., for the user to visually observe it and describe the second edge E2 of the anatomical structure of interest in the second image slice. The second interface 120 is configured to generate a control instruction adapted to output a second subset sd2 of the source data SD so that the command can be generated.

したがって、処理ユニットはまた、第３のインターフェース１３０を介して、第２の画像スライスにおける関心ある解剖学的構造体の第２のエッジＥ２を識別する命令［Ｉ］ｅｄｇｅの第２のセットを受信するように構成される。上記と同様に、命令［Ｉ］ｅｄｇｅの第２のセットは、好ましくは、キーボード及び／又はコンピュータマウス、又は同様のポインティングデバイスを介して入力されるユーザコマンドｕｃｍｄ１及び／又はｕｃｍｄ２に応答して生成される。 Therefore, the processing unit also via a third interface 130, receives a second set of instructions [I] edge identifying a second edge E2 of interest anatomical structure in the second image slices It is configured to do. Similar to the above , the second set of instructions [I] edge is preferably generated in response to user commands ucmd1 and / or ucmd2 entered via a keyboard and / or computer mouse, or similar pointing device. Will be done.

最後に、処理ユニット１４０は、第１のエッジＥ１及び第２のエッジＥ２とソースデータＳＤに基づいて、三次元シェル３ＤＳを計算するように構成される。ここで、三次元シェル３ＤＳは、関心ある解剖学的構造体の境界面の近似を表す。三次元シェル３ＤＳは、ソースデータＳＤと第１のエッジＥ１及び第２のエッジＥ２への幾何学的距離との両方を目的関数に組み込む非線形最適化問題を解くことにより反復的に適応される、変形可能表面モデルによって計算されてよい。この手法は、ユーザ入力、すなわち、第１のエッジＥ１及び第２のエッジＥ２を定義するユーザコマンドｕｃｍｄ１及び／又はｕｃｍｄ２に三次元シェル３ＤＳが密接に結び付けられることを保証するので有益である。これはまた、比較的少ない離散ステップで高品質の結果に達する。第１のサブセットｓｄ１及び第２のサブセットｓｄ２が平行な画像スライスをそれぞれ表すならば、三次元シェル３ＤＳの最初の推定は、表面再構築アルゴリズム、例えば、ユーザにより描かれた輪郭、例えば凸包又は同様のものを近似する三次元表面を計算するための形状に基づく補間を採用するアルゴリズム又は任意の他のアルゴリズムを用いることにより計算されてよい。すなわち、これにより、近似プロセスは、ほとんどのタイプの解剖学的構造体に関して容認できる結果へ非常に迅速に収束する。 Finally, the processing unit 140 is configured to calculate the 3D shell 3DS based on the first edge E1 and the second edge E2 and the source data SD. Here, the 3D shell 3DS represents an approximation of the interface of the anatomical structure of interest. The 3D shell 3DS is iteratively adapted by solving a nonlinear optimization problem that incorporates both the source data SD and the geometric distances to the first edge E1 and the second edge E2 into the objective function. It may be calculated by a deformable surface model. This technique is useful because it ensures that the 3D shell 3DS is tightly coupled to the user input, i.e., the user commands ucmd1 and / or ucmd2 that define the first edge E1 and the second edge E2. It also achieves high quality results with relatively few discrete steps. If the first subset sd1 and the second subset sd2 represent parallel image slices, respectively, then the first estimation of the 3D shell 3DS is a surface reconstruction algorithm, eg, a user-drawn contour, eg, a convex hull or It may be calculated by using an algorithm that employs shape-based interpolation to calculate a three-dimensional surface that approximates a similar one, or any other algorithm. That is, this allows the approximation process to converge very quickly to acceptable results for most types of anatomical structures.

遡れば、ユーザは、既に描かれた輪郭、例えば第１のエッジＥ１及び／又は第２のエッジＥ２が、ソースデータＳＤの画像内容とより良好に一致するように多少異なるアウトラインを有するべきであることに気付くかもしれない。したがって、本発明の一実施形態によれば、第３のインターフェース１３０は、命令［Ｉ］ｅｄｇｅの補足的なセットを処理ユニット１４０に送るように構成され、命令［Ｉ］ｅｄｇｅの補足的なセットは、前記画像スライスのうちの少なくとも１つにおける関心ある解剖学的構造体の調整されたエッジを識別する。また、好ましくは、命令［Ｉ］ｅｄｇｅの補足的なセットは、キーボード及び／又はコンピュータマウス、又は同様のポインティングデバイスを介して入力されているユーザコマンドｕｃｍｄ１及び／又はｕｃｍｄ２に応答して生成される。処理ユニット１４０はさらに、調整されたエッジにさらに基づいて三次元シェル３ＤＳを計算するように構成される。 Going back, the user should have a slightly different outline so that the contours already drawn, eg, the first edge E1 and / or the second edge E2, better match the image content of the source data SD. You may notice that. Therefore, according to one embodiment of the invention, the third interface 130 is configured to send a complementary set of instructions [I] edge to the processing unit 140, which is a complementary set of instructions [I] edge. Identify the adjusted edges of the anatomical structure of interest in at least one of the image slices. Also preferably, a complementary set of instructions [I] edge is generated in response to user commands ucmd1 and / or ucmd2 being entered via a keyboard and / or computer mouse, or similar pointing device. .. The processing unit 140 is further configured to calculate the 3D shell 3DS based further on the adjusted edges.

図３は、三次元シェル３ＤＳを関心ある解剖学的構造体の推定される外周Ｐに近似する三角形メッシュＭの例を示す。三角形メッシュ構造Ｍは、それぞれメッシュ構造Ｍのいくつかの三角形が出会う頂点ｖの組を含む。このタイプの三角形メッシュ表現は、データの効率的なグラフィカル処理及び視覚化を可能にするので有益である。 FIG. 3 shows an example of a triangular mesh M that approximates the 3D shell 3DS to the estimated perimeter P of the anatomical structure of interest. The triangular mesh structure M includes a set of vertices v where several triangles of the mesh structure M meet. This type of triangular mesh representation is useful because it allows for efficient graphical processing and visualization of the data.

したがって、処理ユニット１４０は、好ましくは、三角形分割と、形状に基づく補間に基づいて、三角形メッシュ構造Ｍの第１反復を生成するように構成される。次いで、メッシュ構造Ｍは、第１のエッジＥ１及び第２のエッジＥ２のそれぞれ及びソースデータＳＤと離散数のステップで一致するように連続的に適応される。ここで、現在の近似表面が、適応のための拘束表面として用いられる。最初の近似では、この表面は３ＤＳ’で表される。 Therefore, the processing unit 140 is preferably configured to generate a first iteration of the triangular mesh structure M based on triangulation and shape-based interpolation. The mesh structure M is then continuously adapted to match each of the first edge E1 and the second edge E2 and the source data SD in discrete steps. Here, the current approximate surface is used as the constraining surface for adaptation. In the first approximation, this surface is represented by 3DS'.

好ましくは、処理ユニット１４０は、第１反復に加えて、三角形メッシュ構造Ｍの少なくとも１つのさらなる反復を生成するように構成され、三角形メッシュ構造Ｍの少なくとも１つのさらなる反復は、推定される外周Ｐと三角形メッシュ構造Ｍの頂点ｖとの間の全距離測度に関する第１反復に対する改良である、推定される外周Ｐの近似を表す。図４は、三角形メッシュＭ及び推定される外周Ｐに関連してこの距離測度を概略的に示す。 Preferably, the processing unit 140 is configured to generate at least one additional iteration of the triangular mesh structure M in addition to the first iteration, and at least one additional iteration of the triangular mesh structure M is the estimated perimeter P. Represents an approximation of the estimated perimeter P, which is an improvement over the first iteration for the total distance measurement between and the vertex v of the triangular mesh structure M. FIG. 4 schematically shows this distance measure in relation to the triangular mesh M and the estimated perimeter P.

上記で簡潔に述べたように、推定される外周Ｐを定義する画像特徴の特性は、関心ある構造体の周りの表面にわたってかなり変化し得る。難しいシナリオでは、関心ある構造体は、その境界面のいくつかの部分で非常に強いエッジを有し、他の部分で非常に弱いエッジを有する。後者の部分の近くの、周囲の構造体／組織は、代わりに非常に強いエッジを有し得る。この問題に対処するために、図５は、本発明の一実施形態に従ってエッジＥ１及びＥ２の近傍のソースデータＳＤの一部がどのように解析されるかを例示する。 As briefly mentioned above, the properties of image features that define the estimated perimeter P can vary considerably over the surface around the structure of interest. In difficult scenarios, the structure of interest has very strong edges in some parts of its interface and very weak edges in other parts. Peripheral structures / tissues near the latter part may instead have very strong edges. To address this issue, FIG. 5 illustrates how a portion of the source data SD in the vicinity of edges E1 and E2 is analyzed according to one embodiment of the present invention.

本発明のこの実施形態では、処理ユニット１４０は、第１のエッジＥ１の近傍のソースデータＳＤの第１のサブセットｓｄ１の第１のデータ部分Ｂ１を解析し、第２のエッジＥ２の近傍のソースデータＳＤの第２のサブセットｓｄ２の第２のデータ部分Ｂ２を解析するように構成される。これらの解析に基づいて、処理ユニット１４０は、関心ある解剖学的構造体の推定される外周Ｐを定義する画像特徴の集合を導出し、画像特徴の集合にさらに基づいて、三次元シェル３ＤＳの後続する推定を計算するように構成される。 In this embodiment of the invention, the processing unit 140 analyzes the first data portion B1 of the first subset sd1 of the source data SD near the first edge E1 and the source near the second edge E2. It is configured to analyze the second data portion B2 of the second subset sd2 of the data SD. Based on these analyzes, the processing unit 140 derives a set of image features that define the estimated perimeter P of the anatomical structure of interest, and further based on the set of image features, the 3D shell 3DS. It is configured to calculate subsequent estimates.

好ましくは、それぞれ第１のエッジＥ１及び第２のエッジＥ２の周りで局所的に導出されるこれらの画像特徴は、調和関数を用いて三次元シェル３ＤＳの推定上へ補間される。例えば、処理ユニット１４０は、調和関数を用いることにより最初の推定３ＤＳ’にわたって推定される外周Ｐを定義する画像特徴の集合を補間するように構成することができる。ここで、第１のデータ部分Ｂ１及び第２のデータ部分Ｂ２から導出される画像特徴は、補間プロセスの境界条件を表す。 Preferably, these image features locally derived around the first edge E1 and the second edge E2, respectively, are interpolated onto the estimation of the 3D shell 3DS using a harmonic function. For example, the processing unit 140 can be configured to interpolate a set of image features that define the outer circumference P estimated over the first estimated 3DS'by using a harmonic function. Here, the image features derived from the first data portion B1 and the second data portion B2 represent the boundary conditions of the interpolation process.

これにより、大腿骨頭の骨格などの同じく非常に複雑な表面を識別し、正確にモデル化することができる。磁気共鳴画像における構造体も、ここでの強度が構造体の表面にわたってかなり変化する場合であっても、上手く描写することができる。 This makes it possible to identify and accurately model similarly very complex surfaces such as the skeleton of the femoral head. Structures in magnetic resonance imaging can also be well depicted, even if the intensities here vary significantly over the surface of the structure.

ここで図２に戻り、本発明の別の実施形態を説明する。関心ある解剖学的構造体の形状に応じて、構造体の適切な描写を得るために異なる数のエッジを描くことが適切であり得る。一般に、多くの凹面要素を有する構造体は、全体として凸形の形状を有する構造体よりも多くのエッジをユーザが定義することを必要とする。 Here, returning to FIG. 2, another embodiment of the present invention will be described. Depending on the shape of the anatomical structure of interest, it may be appropriate to draw a different number of edges to obtain a proper depiction of the structure. In general, a structure with many concave elements requires the user to define more edges than a structure with an overall convex shape.

したがって、融通性の向上のために、処理ユニット１４０は、好ましくは、第１のサブセットｓｄ１及び第２のサブセットｓｄ２それぞれに加えて、ソースデータＳＤの少なくとも１つのさらなるサブセットを選択するように構成される。図２は、関心ある解剖学的構造体を通る少なくとも１つのさらなる画像スライスにおいて構成されるソースデータＳＤのさらなる二次元グラフィック表現を定義するこのようなサブセットｓｄ３の例を示す。さらなるサブセットｓｄ３は、ソースデータＳＤにおけるさらなる軸方向の位置ｚ３に存在し、さらなる軸方向の位置ｚ３は、第１の軸方向の位置ｚ１及び第２の軸方向の位置ｚ２とは異なる。 Therefore, for increased flexibility, the processing unit 140 is preferably configured to select at least one additional subset of the source data SD in addition to each of the first subset sd1 and the second subset sd2. To. FIG. 2 shows an example of such a subset sd3 that defines a further two-dimensional graphic representation of the source data SD composed of at least one additional image slice passing through the anatomical structure of interest. A further subset sd3 is present at the additional axial position z3 in the source data SD, and the additional axial position z3 is different from the first axial position z1 and the second axial position z2.

図２に示された例では、さらなるサブセットｓｄ３は、第１のサブセットｓｄ１によって表される画像スライスと平行な画像スライスを表す。このような平行な配向は、三次元シェル３ＤＳが表面再構築アルゴリズム、例えば形状に基づく補間を採用するアルゴリズムを用いることにより計算されることを可能にするので有利である。これは、次に、一般に、近似プロセスがほとんどのタイプの解剖学的構造体に関して容認できる結果へ非常に迅速に収束することにつながる。 In the example shown in FIG. 2, the additional subset sd3 represents an image slice parallel to the image slice represented by the first subset sd1. Such parallel orientation is advantageous because it allows the 3D shell 3DS to be calculated by using a surface reconstruction algorithm, eg, an algorithm that employs shape-based interpolation. This in turn leads to, in general, the approximation process converging very quickly to acceptable results for most types of anatomical structures.

ソースデータＳＤのさらなるサブセットｓｄ３に関して（又は１つ以上存在する場合に、各さらなるサブセットに関して）、処理ユニット１４０は、ユーザがこの画像データを視覚的に観察し得るように、グラフィカルディスプレイ１７０上に提示するために第２のインターフェース１２０にさらなるサブセットｓｄ３を出力させるように適応された制御命令を生成するように構成される。 For a further subset of the source data SD (or for each further subset, if one or more), the processing unit 140 is presented on the graphical display 170 so that the user can visually observe this image data. The second interface 120 is configured to generate a control instruction adapted to output a further subset sd3.

処理ユニット１４０は、第３のインターフェース１３０を介して、さらなる画像スライスにおける関心ある解剖学的構造体の少なくとも１つのさらなるエッジＥ３を識別する命令［Ｉ］ｅｄｇｅのさらなるセットを受信するように構成される。上記と同様に、命令［Ｉ］ｅｄｇｅのさらなるセットは、好ましくは、キーボード及び／又はコンピュータマウス、又は同様のポインティングデバイスを介して入力されているユーザコマンドｕｃｍｄ１及び／又はｕｃｍｄ２に応答して生成される。 Processing unit 140, via a third interface 130, configured to receive a further set of instructions [I] edge identifying at least one additional edge E3 of the anatomical structures of interest in the further image slices To. Similar to the above, an additional set of instructions [I] edge is preferably generated in response to user commands ucmd1 and / or ucmd2 being entered via a keyboard and / or computer mouse, or similar pointing device. To.

処理ユニット１４０は、第１、第２、及びさらなるエッジＥ１、Ｅ２、及びＥ３それぞれとソースデータＳＤに基づいて、三次元シェルの改良バージョン３ＤＳｉを計算するように構成される。ここで、三次元シェルの改良バージョン３ＤＳｉは、関心ある解剖学的構造体の更新された表面推定を表す。 The processing unit 140 is configured to calculate an improved version of the 3D shell, 3DSi, based on the first, second, and additional edges E1, E2, and E3, respectively, and the source data SD. Here, the improved version of the 3D shell, 3DSi, represents an updated surface estimation of the anatomical structure of interest.

処理ユニット１４０は、コンピュータプログラムを実行することにより上記の手順を行うように構成されれば一般に有利である。したがって、処理ユニット１４０は、好ましくは、コンピュータプログラム製品を記憶するメモリユニットに通信可能に接続され、コンピュータプログラム製品は、処理ユニット１４０により実行可能な命令を含み、それにより、処理ユニット１４０は、コンピュータプログラム製品が処理ユニット１４０上で実行されるときに前述のアクションを実行するように動作可能である。 It is generally advantageous if the processing unit 140 is configured to perform the above procedure by executing a computer program. Therefore, the processing unit 140 is preferably communicatively connected to a memory unit that stores the computer program product, the computer program product including instructions that can be executed by the processing unit 140, whereby the processing unit 140 is a computer. It is possible to operate to perform the above-mentioned actions when the program product is executed on the processing unit 140.

要約するために、図６での流れ図を参照しながら、ここで、三次元医用画像データを輪郭描出するための本発明に係る一般的方法を説明する。 For the sake of brevity, a general method according to the present invention for delineating three-dimensional medical image data will be described here with reference to the flow chart in FIG.

最初のステップ６１０において、関心ある解剖学的構造体、並びに、関心ある解剖学的構造体に隣接する組織の３Ｄ画像を表す、ソースデータが受信される。 In the first step 610, source data is received that represents the anatomical structure of interest as well as a 3D image of the tissue adjacent to the anatomical structure of interest.

次いで、ステップ６２０において、ソースデータの第１のサブセットが選択され、グラフィカルディスプレイ上に提示される。ソースデータの第１のサブセットは、ソースデータの第１の二次元グラフィック表現を定義し、第１のサブセットは、ソースデータにおける第１の軸方向の位置で関心ある解剖学的構造体を通る第１の画像スライスにおいて構成される。 Then, in step 620, a first subset of the source data is selected and presented on a graphical display. The first subset of the source data defines the first two-dimensional graphic representation of the source data, and the first subset passes through the anatomical structure of interest at the first axial position in the source data. It is composed of one image slice.

その後、ステップ６３０で、命令の第１のセットが受信されているかをチェックする。好ましくは、命令の第１のセットは、ユーザにより例えばコンピュータマウスを介して入力される手動コマンドに応答して生成される。このような命令の第１のセットが受信されている場合、手順はステップ６４０へ続く。そうでなければ、手順はループバックし、ステップ６３０にとどまる。 Then, in step 630, it is checked whether the first set of instructions has been received. Preferably, the first set of instructions is generated in response to a manual command entered by the user, eg, through a computer mouse. If the first set of such instructions has been received, the procedure continues to step 640. Otherwise, the procedure loops back and stays at step 630.

ステップ６４０において、命令の第１のセットに応答して第１の画像スライスにおける関心ある解剖学的構造体の第１のエッジが識別される。好ましくは、第１のエッジを例示するグラフィックスも、ユーザへのフィードバックとしてグラフィカルディスプレイ上に提示される。 In step 640, the first edge of the anatomical structure of interest in the first image slice is identified in response to the first set of instructions. Preferably, the graphics illustrating the first edge are also presented on the graphical display as feedback to the user.

次いで、ステップ６５０において、ソースデータの第２のサブセットが選択され、グラフィカルディスプレイ上に提示される。ソースデータの第２のサブセットは、ソースデータの第２の二次元グラフィック表現を定義し、第２のサブセットは、ソースデータにおける第２の軸方向の位置で関心ある解剖学的構造体を通る第２の画像スライスにおいて構成され、第２の軸方向の位置は第１の軸方向の位置とは異なる。 Then, in step 650, a second subset of the source data is selected and presented on a graphical display. A second subset of the source data defines a second two-dimensional graphic representation of the source data, and the second subset passes through the anatomical structure of interest at the second axial position in the source data. It is composed of two image slices, and the position in the second axial direction is different from the position in the first axial direction.

その後、ステップ６６０で、命令の第２のセットが受信されているかをチェックする。上記と同様に、命令の第２のセットは、好ましくは、ユーザにより例えばコンピュータマウスを介して入力される手動コマンドに応答して生成される。このような命令の第２のセットが受信されている場合、手順はステップ６７０へ続く。そうでなければ、手順はループバックし、ステップ６６０にとどまる。 Then, in step 660, it is checked whether a second set of instructions has been received. Similar to the above, the second set of instructions is preferably generated in response to manual commands entered by the user, eg, through a computer mouse. If a second set of such instructions has been received, the procedure continues to step 670. Otherwise, the procedure loops back and stays at step 660.

ステップ６７０において、命令の第２のセットに応答して第２の画像スライスにおける関心ある解剖学的構造体の第２のエッジが識別される。 In step 670, the second edge of the anatomical structure of interest in the second image slice is identified in response to the second set of instructions.

その後、ステップ６８０において、第１及び第２のエッジとソースデータに基づいて、三次元シェルが計算される。三次元シェルは、関心ある解剖学的構造体の境界面の近似を表し、したがって、三次元シェルは反復的に計算される。図７は、本発明の一実施形態に従ってこれをどのように行うことができるかを例示する流れ図を示す。 Then, in step 680, a 3D shell is calculated based on the first and second edges and the source data. The 3D shell represents an approximation of the interface of the anatomical structure of interest, so the 3D shell is calculated iteratively. FIG. 7 shows a flow chart illustrating how this can be done according to an embodiment of the present invention.

ここで、最初のステップ７１０において、データポイントの第１及び第２のセット間の補間により、ソースデータにおける第１及び第２の軸方向の位置間の三次元シェルの最初の推定が計算される。データポイントの第１のセットは、第１のエッジの第１のアウトラインを記述し、データポイントの第２のセットは、第２のエッジの第２のアウトラインを記述する。 Here, in the first step 710, interpolation between the first and second sets of data points calculates the first estimate of the 3D shell between the first and second axial positions in the source data. .. The first set of data points describes the first outline of the first edge, and the second set of data points describes the second outline of the second edge.

その後、ステップ７２０において、関心ある解剖学的構造体の推定される外周を定義する画像特徴の集合が導出される。好ましくは、画像特徴の集合は、第１及び第２のエッジのそれぞれの近傍のソースデータの解析に基づいて導出される。 Then, in step 720, a set of image features that define the estimated perimeter of the anatomical structure of interest is derived. Preferably, the set of image features is derived based on the analysis of source data in the vicinity of each of the first and second edges.

次いで、ステップ７３０において、三次元シェルの後続する推定が、三次元シェルの先行する推定（すなわち、第１反復での最初の推定）、第１及び第２のエッジ、画像特徴の集合、及びソースデータに基づいて計算される。 Then, in step 730, the subsequent estimation of the 3D shell is the preceding estimation of the 3D shell (ie, the first estimation in the first iteration), the first and second edges, the set of image features, and the source. Calculated based on the data.

その後、ステップ７４０で、品質基準が満たされているかをチェックし、例えば、基準は、画像情報と第１及び第２のエッジへの幾何学的距離との両方を目的関数に組み込む非線形最適化問題の観点で表される。前記基準が満たされている場合、手順は終了する。そうでなければ、手順はステップ７３０へループバックする。 Then, in step 740, it is checked whether the quality criteria are met, for example, the criteria is a nonlinear optimization problem that incorporates both the image information and the geometric distances to the first and second edges into the objective function. It is expressed from the viewpoint of. If the criteria are met, the procedure ends. Otherwise, the procedure loops back to step 730.

上記の図６及び図７を参照して説明したプロセスステップ並びにステップの任意のサブシーケンスのすべては、プログラムされたプロセッサによって制御され得る。そのうえ、図面を参照して前述した本発明の実施形態は、プロセッサと、少なくとも１つのプロセッサにおいて行われるプロセスを含むが、本発明は、コンピュータプログラム、特に、本発明を実施するように適応されるキャリア上又は内のコンピュータプログラムにも拡張される。プログラムは、部分的にコンパイルされた形態又は本発明に係るプロセスの実装に用いるのに適した任意の他の形態などのソースコード、オブジェクトコード、中間ソースコード、及びオブジェクトコードの形態であり得る。プログラムは、オペレーティングシステムの一部、又は別個のアプリケーションであり得る。キャリアは、プログラムを搭載することができる任意のエンティティ又はデバイスであり得る。例えば、キャリアは、フラッシュメモリ、ＲＯＭ（Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）、例えばＤＶＤ（Ｄｉｇｉｔａｌ Ｖｉｄｅｏ／Ｖｅｒｓａｔｉｌｅ Ｄｉｓｋ）、ＣＤ（Ｃｏｍｐａｃｔ Ｄｉｓｃ）、又は半導体ＲＯＭ、ＥＰＲＯＭ（Ｅｒａｓａｂｌｅ Ｐｒｏｇｒａｍｍａｂｌｅ Ｒｅａｄ−Ｏｎｌｙ Ｍｅｍｏｒｙ）、ＥＥＰＲＯＭ（Ｅｌｅｃｔｒｉｃａｌｌｙ Ｅｒａｓａｂｌｅ Ｐｒｏｇｒａｍｍａｂｌｅ Ｒｅａｄ−Ｏｎｌｙ Ｍｅｍｏｒｙ）などの記憶媒体、又は磁気記録媒体、例えばフロッピーディスク又はハードディスクを含み得る。さらに、キャリアは、電気ケーブル又は光ケーブルを介して若しくは無線又は他の手段により伝達され得る電気信号又は光信号などの伝送可能なキャリアであり得る。プログラムが、ケーブル又は他のデバイス又は手段により直接伝達され得る信号で具体化されるとき、キャリアは、このようなケーブル又はデバイス又は手段により構成され得る。代替的に、キャリアは、プログラムが組み込まれる集積回路であってよく、集積回路は、該当するプロセスを行う又は行う際に用いるように適応される。 All of the process steps described with reference to FIGS. 6 and 7 above and any subsequence of steps can be controlled by a programmed processor. Moreover, although the embodiments of the present invention described above with reference to the drawings include a processor and a process performed on at least one processor, the present invention is adapted to implement computer programs, in particular the present invention. It also extends to computer programs on or within the carrier. The program can be in the form of source code, object code, intermediate source code, and object code, such as a partially compiled form or any other form suitable for use in implementing a process according to the invention. The program can be part of an operating system or a separate application. The carrier can be any entity or device that can carry the program. For example, the carrier may be a flash memory, a ROM (Read Only Memory), for example, a DVD (Digital Video / Versaille Disk), a CD (Compact Disk), or a semiconductor ROM, an EPROM (Erasable Programmable Read-Ollemagry Memory). It may include a storage medium such as a Read-Only Memory) or a magnetic recording medium such as a floppy disk or a hard disk. Further, the carrier can be a transmittable carrier such as an electrical or optical signal that can be transmitted via an electrical or optical cable or by radio or other means. When the program is embodied in a signal that can be transmitted directly by a cable or other device or means, the carrier may be composed of such a cable or device or means. Alternatively, the carrier may be an integrated circuit into which the program is incorporated, and the integrated circuit is adapted for use in performing or performing the relevant process.

本明細書で用いられるときの「備える、含む」という用語は、表記された特徴、整数、ステップ、又はコンポーネントの存在を明記するために採用される。しかしながら、該用語は、１つ以上の付加的な特徴、整数、ステップ、又はコンポーネント、又はその群の存在又は追加を除外しない。 The term "preparing, including" as used herein is adopted to specify the presence of a described feature, integer, step, or component. However, the term does not exclude the presence or addition of one or more additional features, integers, steps, or components, or groups thereof.

本発明は、図面で説明した実施形態に限定されず、請求項の範囲内で自由に変更され得る。 The present invention is not limited to the embodiments described in the drawings and may be freely modified within the scope of the claims.

Claims (10)

処理ユニット（１４０）と、第１、第２、及び第３のインターフェース（１１０、１２０、１３０）とを備え、
前記第１のインターフェース（１１０）は、処理ユニット（１４０）が三次元医用画像データのソースデータ（ＳＤ）にアクセスすることを可能にするように構成され、前記ソースデータ（ＳＤ）は輪郭描出されるべき関心ある解剖学的構造体及び前記関心ある解剖学的構造体に隣接する組織を表し、
前記処理ユニット（１４０）が、前記ソースデータ（ＳＤ）における第１の軸方向の位置（ｚ１）で前記関心ある解剖学的構造体を通る第１の画像スライスにおいて構成される前記ソースデータ（ＳＤ）の第１の二次元グラフィック表現を定義する前記ソースデータ（ＳＤ）の第１のサブセット（ｓｄ１）を選択するように構成され、
前記第２のインターフェース（１２０）が、前記処理ユニット（１４０）からの制御命令に応答して、グラフィカルディスプレイ（１７０）上に提示するために前記ソースデータ（ＳＤ）の第１のサブセット（ｓｄ１）を出力するように構成され、
前記第３のインターフェース（１３０）が、命令（［Ｉ］ｅｄｇｅ）の第１のセットを前記処理ユニット（１４０）に送るように構成され、前記命令（［Ｉ］ｅｄｇｅ）の第１のセットは、ユーザコマンド（ｕｃｍｄ１、ｕｃｍｄ２）の形態で前記第３のインターフェース（１３０）に受信され、第１の画像スライスにおける関心ある解剖学的構造体の第１のエッジ（Ｅ１）を識別し、
前記命令（［Ｉ］ｅｄｇｅ）の第１のセットを受信した後で、前記処理ユニット（１４０）が、
前記ソースデータ（ＳＤ）における第２の軸方向の位置（ｚ２）で前記関心ある解剖学的構造体を通る第２の画像スライスにおいて構成される前記ソースデータ（ＳＤ）の第２の二次元グラフィック表現を定義する前記ソースデータ（ＳＤ）の第２のサブセット（ｓｄ２）を選択し、前記第２の軸方向の位置（ｚ２）は前記第１の軸方向の位置（ｚ１）とは異なり、
前記グラフィカルディスプレイ（１７０）上に提示するために前記ソースデータ（ＳＤ）の前記第２のサブセット（ｓｄ２）を前記第２のインターフェース（１２０）に出力させるように構成された制御命令を生成し、
前記第３のインターフェース（１３０）を介して、前記第２の画像スライスにおける前記関心ある解剖学的構造体の第２のエッジ（Ｅ２）を識別する命令（［Ｉ］ｅｄｇｅ）の第２のセットを受信し、
前記第１及び第２のエッジ（Ｅ１、Ｅ２）と前記ソースデータ（ＳＤ）に基づいて三次元シェル（３ＤＳ）を計算する、
ように構成され、前記三次元シェル（３ＤＳ）が前記関心ある解剖学的構造体の境界面の近似を表し、
データポイントの第１及び第２のセット間の補間により前記ソースデータ（ＳＤ）における第１及び第２の軸方向の位置（ｚ１、ｚ２）間の三次元シェル（３ＤＳ）の最初の推定（３ＤＳ’）を計算し、前記データポイントの第１のセットが前記第１のエッジ（Ｅ１）の第１のアウトラインを記述し、前記データポイントの第２のセットが前記第２のエッジ（Ｅ２）の第２のアウトラインを記述し、
前記第１のエッジ（Ｅ１）の近傍の前記ソースデータ（ＳＤ）の前記第１のサブセット（ｓｄ１）の第１のデータ部分（Ｂ１）を解析し、
前記第２のエッジ（Ｅ２）の近傍の前記ソースデータ（ＳＤ）の前記第２のサブセット（ｓｄ２）の第２のデータ部分（Ｂ２）を解析し、これらに基づいて、
前記関心ある解剖学的構造体の推定される外周（Ｐ）を定義する画像特徴の集合を導出し、
前記画像特徴の集合にさらに基づいて前記三次元シェル（３ＤＳ）の後続する推定を計算する、
三次元医用画像データを輪郭描出するための画像処理システム（１００）であって、
前記処理ユニット（１４０）が、調和関数を用いることにより前記最初の推定（３ＤＳ’）にわたって前記推定される外周（Ｐ）を定義する画像特徴の集合を補間するようにさらに構成され、前記第１及び第２のデータ部分（Ｂ１、Ｂ２）から導出される画像特徴が境界条件を表す、
ことを特徴とする、画像処理システム（１００）。 It comprises a processing unit (140) and first, second, and third interfaces (110, 120, 130).
The first interface (110) is configured to allow the processing unit (140) to access the source data (SD) of the three-dimensional medical image data, and the source data (SD) is outlined. Represents the anatomical structure of interest and the tissue adjacent to the anatomical structure of interest.
The source data (SD) in which the processing unit (140) is composed of a first image slice passing through the anatomical structure of interest at a first axial position (z1) in the source data (SD). ) Is configured to select the first subset (sd1) of the source data (SD) that defines the first two-dimensional graphic representation.
A first subset (sd1) of the source data (SD) for the second interface (120) to present on the graphical display (170) in response to a control instruction from the processing unit (140). Is configured to output
The third interface (130) is configured to send a first set of instructions ([I] edge ) to the processing unit (140), and the first set of instructions ([I] edge) is , Received by the third interface (130) in the form of user commands (umdd1, ucmd2) to identify the first edge (E1) of the anatomical structure of interest in the first image slice.
After receiving the first set of instructions ([I] edge ), the processing unit (140)
A second two-dimensional graphic of the source data (SD) composed of a second image slice passing through the anatomical structure of interest at a second axial position (z2) in the source data (SD). A second subset (sd2) of the source data (SD) that defines the representation is selected, and the second axial position (z2) is different from the first axial position (z1).
A control instruction configured to output the second subset (sd2) of the source data (SD) to the second interface (120) for presentation on the graphical display (170) is generated.
A second set of instructions ([I] edge ) that identify the second edge (E2) of the anatomical structure of interest in the second image slice via the third interface (130). Received
A three-dimensional shell (3DS) is calculated based on the first and second edges (E1, E2) and the source data (SD).
The three-dimensional shell (3DS) represents an approximation of the interface of the anatomical structure of interest.
The first estimation (3DS) of a three-dimensional shell (3DS) between the first and second axial positions (z1, z2) in the source data (SD) by interpolation between the first and second sets of data points. ') Is calculated, the first set of the data points describes the first outline of the first edge (E1), and the second set of the data points is of the second edge (E2). Write a second outline,
The first data portion (B1) of the first subset (sd1) of the source data (SD) in the vicinity of the first edge (E1) is analyzed.
The second data portion (B2) of the second subset (sd2) of the source data (SD) in the vicinity of the second edge (E2) is analyzed and based on these, the second data portion (B2) is analyzed.
We derive a set of image features that define the estimated perimeter (P) of the anatomical structure of interest.
Compute subsequent estimates of the 3D shell (3DS) based further on the set of image features.
An image processing system (100) for delineating three-dimensional medical image data.
The processing unit (140) is further configured to interpolate a set of image features defining the estimated outer circumference (P) over the first estimate (3DS') by using a harmonic function. And the image features derived from the second data part (B1, B2) represent the boundary conditions.
An image processing system (100), characterized in that. 前記第１の画像スライスが前記ソースデータ（ＳＤ）において第１の配向を有し、前記第２の画像スライスが前記第１の配向と平行に配向され、
前記処理ユニット（１４０）が、表面再構築アルゴリズムを用いることにより前記三次元シェル（３ＤＳ）の最初の推定（３ＤＳ’）を計算するように構成される、
請求項１に記載の画像処理システム（１００）。 The first image slice has a first orientation in the source data (SD) and the second image slice is oriented parallel to the first orientation.
The processing unit (140) is configured to calculate the first estimate (3DS') of the three-dimensional shell (3DS) by using a surface reconstruction algorithm.
The image processing system (100) according to claim 1. 前記処理ユニット（１４０）が、前記三次元シェル（３ＤＳ）の最初の推定（３ＤＳ’）、前記第１及び第２のエッジ（Ｅ１、Ｅ２）、及び前記ソースデータ（ＳＤ）に基づいて、前記三次元シェル（３ＤＳ）の後続する推定を計算するように構成される、請求項２に記載の画像処理システム（１００）。 The processing unit (140) is based on the first estimation (3DS') of the three-dimensional shell (3DS), the first and second edges (E1, E2), and the source data (SD). The image processing system (100) according to claim 2, configured to calculate subsequent estimates of the three-dimensional shell (3DS). 前記処理ユニット（１４０）がさらに、
前記三次元シェル（３ＤＳ）を推定される前記外周（Ｐ）に近似する三角形メッシュ構造（Ｍ）の第１反復を生成するように構成され、前記三角形メッシュ構造（Ｍ）は、前記三角形メッシュ構造（Ｍ）のいくつかの三角形が出会う頂点（ｖ）の組を含む、請求項１に記載の画像処理システム（１００）。 The processing unit (140) further
The three-dimensional shell (3DS) is configured to generate a first iteration of a triangular mesh structure (M) that approximates the estimated perimeter (P), and the triangular mesh structure (M) is the triangular mesh structure. The image processing system (100) according to claim 1, comprising a set of vertices (v) where several triangles of (M) meet. 前記処理ユニット（１４０）がさらに、前記第１反復に加えて、前記三角形メッシュ構造（Ｍ）の少なくとも１つのさらなる反復を生成するように構成され、前記三角形メッシュ構造（Ｍ）の少なくとも１つのさらなる反復は、前記推定される外周（Ｐ）と前記三角形メッシュ構造（Ｍ）の頂点（ｖ）との間の全距離測度に関する前記第１反復に対する改良である、前記推定される外周（Ｐ）の近似を表す、請求項４に記載の画像処理システム（１００）。 The processing unit (140) is further configured to generate at least one additional iteration of the triangular mesh structure (M) in addition to the first iteration, and at least one additional iteration of the triangular mesh structure (M). The iteration is an improvement over the first iteration with respect to the total distance measurement between the estimated outer circumference (P) and the apex (v) of the triangular mesh structure (M) of the estimated outer circumference (P). The image processing system (100) according to claim 4, which represents an approximation. 前記処理ユニット（１４０）が、
前記第１及び第２のサブセット（ｓｄ１、ｓｄ２）に加えて、前記ソースデータ（ＳＤ）の少なくとも１つのさらなるサブセット（ｓｄ３）を選択し、前記少なくとも１つのさらなるサブセット（ｓｄ３）は、前記ソースデータ（ＳＤ）における少なくとも１つのさらなる軸方向の位置（ｚ３）で前記関心ある解剖学的構造体を通る少なくとも１つのさらなる画像スライスにおいて構成される前記ソースデータ（ＳＤ）の少なくとも１つのさらなる二次元グラフィック表現を定義し、前記少なくとも１つのさらなる軸方向の位置（ｚ３）は、前記第１及び第２の軸方向の位置（ｚ１、ｚ２）とは異なり、
前記少なくとも１つのさらなるサブセット（ｓｄ３）のそれぞれに関して、前記グラフィカルディスプレイ（１７０）上に提示するために前記第２のインターフェース（１２０）に前記ソースデータ（ＳＤ）の少なくとも１つのさらなるサブセット（ｓｄ３）を出力させるように構成された制御命令を生成し、
前記第３のインターフェース（１３０）を介して、前記少なくとも１つのさらなる画像スライスにおける前記関心ある解剖学的構造体の少なくとも１つのさらなるエッジ（Ｅ３）を識別する命令（［Ｉ］ｅｄｇｅ）の少なくとも１つのさらなるセットを受信し、
前記第１、第２、及び少なくとも１つのさらなるエッジ（Ｅ１、Ｅ２、Ｅ３）と前記ソースデータ（ＳＤ）に基づいて、前記関心ある解剖学的構造体の表面推定を表す前記三次元シェルの改良バージョン（３ＤＳｉ）を計算する、
ように構成される、請求項２〜５のいずれか一項に記載の画像処理システム（１００）。 The processing unit (140)
In addition to the first and second subsets (sd1, sd2), at least one additional subset (sd3) of the source data (SD) is selected and the at least one additional subset (sd3) is the source data. At least one additional two-dimensional graphic of the source data (SD) composed of at least one additional image slice passing through the anatomical structure of interest at at least one additional axial position (z3) in (SD). Defining the representation, said at least one additional axial position (z3) is different from said first and second axial positions (z1, z2).
For each of the at least one additional subset (sd3), the second interface (120) is provided with at least one additional subset (sd3) of the source data (SD) for presentation on the graphical display (170). Generates a control command configured to output and
At least one of the instructions ([I] edge ) that identifies at least one additional edge (E3) of the anatomical structure of interest in said at least one additional image slice via the third interface (130). Receive one more set,
Improvements to the three-dimensional shell that represent surface estimates of the anatomical structure of interest based on the first, second, and at least one additional edge (E1, E2, E3) and the source data (SD). Calculate the version (3DSi),
The image processing system (100) according to any one of claims 2 to 5, which is configured as described above. 前記第３のインターフェース（１３０）が、命令（［Ｉ］ｅｄｇｅ）の補足的なセットを前記処理ユニット（１４０）に送るように構成され、前記命令（［Ｉ］ｅｄｇｅ）の補足的なセットは、前記少なくとも１つのさらなる画像スライスにおける前記関心ある解剖学的構造体の調整されたエッジを識別し、
前記処理ユニット（１４０）が、前記調整されたエッジにさらに基づいて三次元シェル（３ＤＳ）を計算するように構成される、
請求項６に記載の画像処理システム（１００）。 The third interface (130) is configured to send a complementary set of instructions ([I] edge ) to the processing unit (140), and the complementary set of instructions ([I] edge) is Identify the adjusted edges of the anatomical structure of interest in the at least one additional image slice.
The processing unit (140) is configured to calculate a three-dimensional shell (3DS) based further on the adjusted edges.
The image processing system (100) according to claim 6. 第１のインターフェース（１１０）を介して、三次元医用画像データのソースデータ（ＳＤ）を受信することであって、前記ソースデータ（ＳＤ）は輪郭描出されるべき関心ある解剖学的構造体及び前記関心ある解剖学的構造体に隣接する組織を表す、受信することと、
前記ソースデータ（ＳＤ）の第１の二次元グラフィック表現を定義する前記ソースデータ（ＳＤ）の第１のサブセット（ｓｄ１）を選択することであって、前記第１のサブセット（ｓｄ１）は、前記ソースデータ（ＳＤ）における第１の軸方向の位置（ｚ１）で前記関心ある解剖学的構造体を通る第１の画像スライスにおいて構成される、第１のサブセット（ｓｄ１）を選択することと、
第２のインターフェース（１２０）を介して、グラフィカルディスプレイ（１７０）上に提示するために前記ソースデータ（ＳＤ）の第１のサブセット（ｓｄ１）を出力することと、
ユーザコマンド（ｕｃｍｄ１、ｕｃｍｄ２）の形態で、第３のインターフェース（１３０）において命令（［Ｉ］ｅｄｇｅ）の第１のセットを受信することと、
第３のインターフェース（１３０）を介して受信した前記命令（［Ｉ］ｅｄｇｅ）の第１のセットに応答して、前記第１の画像スライスにおける前記関心ある解剖学的構造体の第１のエッジ（Ｅ１）を識別することと、
前記命令（［Ｉ］ｅｄｇｅ）の第１のセットを受信した後で、
前記ソースデータ（ＳＤ）の第２の二次元グラフィック表現を定義する前記ソースデータ（ＳＤ）の第２のサブセット（ｓｄ２）を選択することであって、前記第２のサブセット（ｓｄ２）は、前記ソースデータ（ＳＤ）における第２の軸方向の位置（ｚ２）で前記関心ある解剖学的構造体を通る第２の画像スライスにおいて構成され、前記第２の軸方向の位置（ｚ２）は前記第１の軸方向の位置（ｚ１）とは異なる、第２のサブセット（ｓｄ２）を選択することと、
前記第２のインターフェース（１２０）を介して、前記グラフィカルディスプレイ（１７０）上に提示するために前記ソースデータ（ＳＤ）の前記第２のサブセット（ｓｄ２）を出力することと、
前記第３のインターフェース（１３０）を介して受信した前記命令（［Ｉ］ｅｄｇｅ）の第２のセットに応答して、前記第２の画像スライスにおける前記関心ある解剖学的構造体の第２のエッジ（Ｅ２）を識別することと、
前記第１及び第２のエッジ（Ｅ１、Ｅ２）と前記ソースデータ（ＳＤ）に基づいて、前記関心ある解剖学的構造体の境界面の近似を表す三次元シェル（３ＤＳ）を計算することと、
データポイントの第１及び第２のセット間の補間により前記ソースデータ（ＳＤ）における前記第１及び第２の軸方向の位置（ｚ１、ｚ２）間の三次元シェル（３ＤＳ）の最初の推定（３ＤＳ’）を計算することであって、前記データポイントの第１のセットが前記第１のエッジ（Ｅ１）の第１のアウトラインを記述し、前記データポイントの第２のセットが前記第２のエッジ（Ｅ２）の第２のアウトラインを記述する、計算することと、
前記第１のエッジ（Ｅ１）の近傍の前記ソースデータ（ＳＤ）の前記第１のサブセット（ｓｄ１）の第１のデータ部分（Ｂ１）を解析することと、
前記第２のエッジ（Ｅ２）の近傍の前記ソースデータ（ＳＤ）の前記第２のサブセット（ｓｄ２）の第２のデータ部分（Ｂ２）を解析することと、これらに基づいて、
前記関心ある解剖学的構造体の推定される外周（Ｐ）を定義する画像特徴の集合を導出することと、
前記画像特徴の集合にさらに基づいて前記三次元シェル（３ＤＳ）の後続する推定を計算することと、
を含む、三次元医用画像データを輪郭描出するプロセッサにより実施される方法であって、
調和関数を用いることにより前記最初の推定（３ＤＳ’）にわたって前記推定される外周（Ｐ）を定義する画像特徴の集合を補間し、前記第１及び第２のデータ部分（Ｂ１、Ｂ２）から導出される画像特徴が境界条件を表す、
ことを特徴とする方法。 The source data (SD) of the three-dimensional medical image data is received via the first interface (110), the source data (SD) being the anatomical structure of interest to be contoured and Receiving and representing the tissue adjacent to the anatomical structure of interest,
Selecting a first subset (sd1) of the source data (SD) that defines a first two-dimensional graphic representation of the source data (SD), wherein the first subset (sd1) is said. Selecting a first subset (sd1) composed of a first image slice passing through the anatomical structure of interest at a first axial position (z1) in the source data (SD).
To output a first subset (sd1) of said source data (SD) for presentation on a graphical display (170) via a second interface (120).
Receiving a first set of instructions ([I] edge ) on a third interface (130) in the form of user commands (umdd1, ucmd2).
The first edge of the anatomical structure of interest in the first image slice in response to a first set of said instructions ([I] edge) received via a third interface (130). Identifying (E1) and
After receiving the first set of instructions ([I] edge),
Selecting a second subset (sd2) of the source data (SD) that defines a second two-dimensional graphic representation of the source data (SD), wherein the second subset (sd2) is said. The second axial position (z2) in the source data (SD) is composed of a second image slice passing through the anatomical structure of interest, and the second axial position (z2) is said to be the second. Choosing a second subset (sd2) that is different from the axial position of 1 (z1) and
To output the second subset (sd2) of the source data (SD) for presentation on the graphical display (170) via the second interface (120).
A second of the anatomical structures of interest in the second image slice in response to a second set of instructions ([I] edges ) received via the third interface (130). Identifying the edge (E2) and
Based on the first and second edges (E1, E2) and the source data (SD), calculating a three-dimensional shell (3DS) that represents an approximation of the interface of the anatomical structure of interest. ,
The first estimation of the 3D shell (3DS) between the first and second axial positions (z1, z2) in the source data (SD) by interpolation between the first and second sets of data points ( 3DS'), where the first set of data points describes the first outline of the first edge (E1) and the second set of data points is the second. Writing, calculating, and describing the second outline of the edge (E2),
Analyzing the first data portion (B1) of the first subset (sd1) of the source data (SD) in the vicinity of the first edge (E1).
Analyzing the second data portion (B2) of the second subset (sd2) of the source data (SD) in the vicinity of the second edge (E2) and based on these.
To derive a set of image features that define the estimated perimeter (P) of the anatomical structure of interest.
Computing subsequent estimates of the 3D shell (3DS) based further on the set of image features.
A method performed by a processor that outlines three-dimensional medical image data, including
By using the harmonic function, the set of image features that define the estimated outer circumference (P) is interpolated over the first estimation (3DS') and derived from the first and second data portions (B1, B2). The image features that are created represent the boundary conditions,
A method characterized by that. 少なくとも１つの処理ユニット（１４０）のメモリ（１５０）にロード可能であり、且つ、前記少なくとも１つの処理ユニット（１４０）上で実行されるときに請求項８に記載の方法を実行するためのソフトウェアを含む、コンピュータプログラム。 Software for performing the method of claim 8 when loaded into memory (150) of at least one processing unit (140) and executed on said at least one processing unit (140). Including computer programs. 少なくとも１つの処理ユニット（１４０）にロードされるときに前記少なくとも１つの処理ユニット（１４０）に請求項８に記載の方法を実行させるためのプログラムが記録されている、プロセッサ可読媒体（１５０）。 The at least one processing unit (140) to the program for executing the method according to claim 8 is recorded, the processor-readable medium when loaded to one processing unit (140) even without least (150 ).
JP2019512197A 2016-09-05 2017-09-04 Image processing systems and methods for interactive contouring of 3D medical data Active JP6877534B2 (en)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP16187218.9 2016-09-05
EP16187218.9A EP3291175B1 (en) 2016-09-05 2016-09-05 Image processing system and method for interactive contouring of three-dimensional medical data
PCT/EP2017/072135 WO2018042043A1 (en) 2016-09-05 2017-09-04 Image processing system and method for interactive contouring of three-dimensional medical data

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2020123219A Division JP2020195779A (en) 2016-09-05 2020-07-17 Image processing system and method for three-dimensional medical data interactive contour delineation

Publications (2)

Publication Number Publication Date
JP2019535336A JP2019535336A (en) 2019-12-12
JP6877534B2 true JP6877534B2 (en) 2021-05-26

Family

ID=56943318

Family Applications (2)

Application Number Title Priority Date Filing Date
JP2019512197A Active JP6877534B2 (en) 2016-09-05 2017-09-04 Image processing systems and methods for interactive contouring of 3D medical data
JP2020123219A Pending JP2020195779A (en) 2016-09-05 2020-07-17 Image processing system and method for three-dimensional medical data interactive contour delineation

Family Applications After (1)

Application Number Title Priority Date Filing Date
JP2020123219A Pending JP2020195779A (en) 2016-09-05 2020-07-17 Image processing system and method for three-dimensional medical data interactive contour delineation

Country Status (5)

Country Link
US (1) US10453199B2 (en)
EP (1) EP3291175B1 (en)
JP (2) JP6877534B2 (en)
CN (1) CN109716389B (en)
WO (1) WO2018042043A1 (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3291175B1 (en) * 2016-09-05 2018-12-12 RaySearch Laboratories AB Image processing system and method for interactive contouring of three-dimensional medical data
CN115145017B (en) * 2018-01-25 2023-12-22 台湾东电化股份有限公司 Optical system
CN109933862A (en) * 2019-02-26 2019-06-25 中国人民解放军军事科学院国防科技创新研究院 A kind of electromagnetic model construction method and device suitable for magnetic induction spectrum emulation
WO2021124683A1 (en) * 2019-12-16 2021-06-24 富士フイルム株式会社 Contour extracting device, contour extracting method, and contour extracting program
WO2021133811A1 (en) 2019-12-26 2021-07-01 Teva Pharmaceuticals International Gmbh Solid state forms of cenicriviroc and process for preparation thereof
US11455733B2 (en) * 2020-06-16 2022-09-27 Xue Feng System and method to improve model-based organ segmentation with image post-processing
US11974881B2 (en) * 2020-08-26 2024-05-07 GE Precision Healthcare LLC Method and system for providing an anatomic orientation indicator with a patient-specific model of an anatomical structure of interest extracted from a three-dimensional ultrasound volume
KR102667231B1 (en) 2020-12-21 2024-05-20 주식회사 인피니트헬스케어 System and method for contouring a set of medical images based on deep learning algorighm and anatomical characteristics

Family Cites Families (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0670923A (en) * 1992-08-28 1994-03-15 Hitachi Medical Corp Extraction method for specific area of medical image
JPH08336503A (en) * 1995-06-15 1996-12-24 Ge Yokogawa Medical Syst Ltd Medical image diagnostic system
JP3692050B2 (en) * 2001-04-12 2005-09-07 株式会社東芝 Image processing apparatus and method
JP3725442B2 (en) * 2001-04-18 2005-12-14 株式会社東芝 Medical image diagnostic apparatus and method
US8358819B2 (en) * 2005-06-24 2013-01-22 University Of Iowa Research Foundation System and methods for image segmentation in N-dimensional space
US20090003666A1 (en) * 2007-06-27 2009-01-01 Wu Dee H System and methods for image analysis and treatment
US8311306B2 (en) * 2008-04-30 2012-11-13 Otismed Corporation System and method for image segmentation in generating computer models of a joint to undergo arthroplasty
CN101807296A (en) * 2009-02-16 2010-08-18 海信集团有限公司 Segmentation method of medical ultrasonic image three-dimensional target object
US8805035B2 (en) 2010-05-03 2014-08-12 Mim Software, Inc. Systems and methods for contouring a set of medical images
US8983189B2 (en) * 2011-02-08 2015-03-17 Algotec Systems Ltd. Method and systems for error correction for three-dimensional image segmentation
EP2847605B1 (en) * 2012-05-09 2019-01-02 Laboratoires Bodycad Inc. Segmentation of magnetic resonance imaging data
US10827983B2 (en) * 2012-10-30 2020-11-10 The Johns Hopkins University System and method for personalized cardiac arrhythmia risk assessment by simulating arrhythmia inducibility
JP6273291B2 (en) * 2012-12-03 2018-01-31 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． Image processing apparatus and method
CN103903255B (en) * 2012-12-31 2016-08-31 中国科学院深圳先进技术研究院 A kind of ultrasonic image division method and system
JP6396420B2 (en) * 2013-04-03 2018-09-26 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． 3D ultrasound imaging system, corresponding method and computer program for implementing such a method
CN104091365B (en) * 2014-07-12 2017-11-10 大连理工大学 Towards the acetabular bone tissue model reconstruction method of serializing hip joint CT images
US10147185B2 (en) * 2014-09-11 2018-12-04 B.G. Negev Technologies And Applications Ltd., At Ben-Gurion University Interactive segmentation
US10013810B2 (en) * 2014-09-26 2018-07-03 Varian Medical Systems International Ag. Contouring tool having automatic interpolation and extrapolation
EP3234926B1 (en) * 2014-12-18 2020-11-04 Koninklijke Philips N.V. Medical image editing
US9978150B2 (en) * 2015-08-05 2018-05-22 Algotec Systems Ltd. Method and system for spatial segmentation of anatomical structures
CN105787978A (en) * 2016-02-29 2016-07-20 深圳市医诺智能科技发展有限公司 Automatic medical image interlayer sketching method, device and system
EP3291175B1 (en) * 2016-09-05 2018-12-12 RaySearch Laboratories AB Image processing system and method for interactive contouring of three-dimensional medical data

Also Published As

Publication number Publication date
EP3291175A1 (en) 2018-03-07
WO2018042043A1 (en) 2018-03-08
JP2020195779A (en) 2020-12-10
US10453199B2 (en) 2019-10-22
CN109716389A (en) 2019-05-03
CN109716389B (en) 2020-04-14
JP2019535336A (en) 2019-12-12
US20190197692A1 (en) 2019-06-27
EP3291175B1 (en) 2018-12-12

Similar Documents

Publication Publication Date Title
JP6877534B2 (en) Image processing systems and methods for interactive contouring of 3D medical data
US9710880B2 (en) User-guided shape morphing in bone segmentation for medical imaging
EP2810598B1 (en) Surgical support device, surgical support method and surgical support program
US8582854B2 (en) Method and system for automatic coronary artery detection
JP5113077B2 (en) A method to facilitate post-processing of images using deformable meshes
US11132801B2 (en) Segmentation of three-dimensional images containing anatomic structures
CN108038862B (en) Interactive medical image intelligent segmentation modeling method
EP2856428B1 (en) Segmentation highlighter
EP2827301B1 (en) Image generation device, method, and program
JP2009072576A (en) Method, apparatus and computer program for quantitative bifurcation analysis on angiographic image
US9649167B2 (en) Pattern and surgery support set, apparatus, method and program
JP5748636B2 (en) Image processing apparatus and method, and program
JP2022526418A (en) Modeling the area of interest of the anatomy
CN117115150B (en) Method, computing device and medium for determining branch vessels
CN107851337B (en) Interactive grid editing
JP5231210B2 (en) Method, system and computer program for segmenting a surface in a multidimensional data set
JP6678819B2 (en) Image processing system and method for interactive contouring of three-dimensional medical data
JP2008535614A5 (en)
US10467761B2 (en) System and method for handling image data
JP5954846B2 (en) Shape data generation program, shape data generation method, and shape data generation apparatus
Reska et al. Fast 3D segmentation of hepatic images combining region and boundary criteria
JP2007512890A (en) How to determine the structure of a moving object

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190614

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190614

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20190614

A975 Report on accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A971005

Effective date: 20191017

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20191112

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200210

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20200317

C60 Trial request (containing other claim documents, opposition documents)

Free format text: JAPANESE INTERMEDIATE CODE: C60

Effective date: 20200717

C116 Written invitation by the chief administrative judge to file amendments

Free format text: JAPANESE INTERMEDIATE CODE: C116

Effective date: 20200804

C22 Notice of designation (change) of administrative judge

Free format text: JAPANESE INTERMEDIATE CODE: C22

Effective date: 20200804

C22 Notice of designation (change) of administrative judge

Free format text: JAPANESE INTERMEDIATE CODE: C22

Effective date: 20201020

C23 Notice of termination of proceedings

Free format text: JAPANESE INTERMEDIATE CODE: C23

Effective date: 20210222

C03 Trial/appeal decision taken

Free format text: JAPANESE INTERMEDIATE CODE: C03

Effective date: 20210330

C30A Notification sent

Free format text: JAPANESE INTERMEDIATE CODE: C3012

Effective date: 20210330

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20210427

R150 Certificate of patent or registration of utility model

Ref document number: 6877534

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE

Ref document number: 6877534

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250