JP6937187B2 - Breast shape change prediction method, breast shape change prediction system, and breast shape change prediction program - Google Patents
Breast shape change prediction method, breast shape change prediction system, and breast shape change prediction program Download PDFInfo
- Publication number
- JP6937187B2 JP6937187B2 JP2017157284A JP2017157284A JP6937187B2 JP 6937187 B2 JP6937187 B2 JP 6937187B2 JP 2017157284 A JP2017157284 A JP 2017157284A JP 2017157284 A JP2017157284 A JP 2017157284A JP 6937187 B2 JP6937187 B2 JP 6937187B2
- Authority
- JP
- Japan
- Prior art keywords
- breast
- model
- tension line
- vector
- length
- 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
Links
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Description
本発明は、乳がんなどの患部を摘出するにあたって、***の形状がどのように変化するかを予測するための技術に関するものである。 The present invention relates to a technique for predicting how the shape of a breast changes when an affected part such as breast cancer is removed.
現在、乳がんなどの***に係わる病気を患った場合には、例えば患部の大きさなどを考慮して、***のほぼ全てを摘出する***切除術や、***の一部を切除する***温存術(例えば特許文献1参照)によって治療が行われている。***温存術で切除される部分の大きさは医療技術の進歩に伴って以前よりも小さくなっていて、手術前後での***の形状変化も目立ちにくくなっている。 Currently, when suffering from a breast-related disease such as breast cancer, mastectomy to remove almost all of the breast or breast-conserving surgery to remove a part of the breast (for example, considering the size of the affected area) For example, the treatment is performed according to Patent Document 1). The size of the part to be excised by breast-conserving surgery has become smaller than before due to advances in medical technology, and changes in the shape of the breast before and after surgery have become less noticeable.
このように患部摘出に伴う***の形状変化は小さくなってきたものの、患者にとっては***の形状がどの程度変化するのか分かり難いため、手術に不安を感じることがある。また医師にとっても患者への具体的な説明が難しく、負担を抱えている。 Although the change in the shape of the breast due to the removal of the affected area has become smaller in this way, it is difficult for the patient to understand how much the shape of the breast changes, and the surgery may cause anxiety. In addition, it is difficult for doctors to give specific explanations to patients, which is a burden.
本発明は、このような問題点に鑑みてなされたものであり、患部摘出に伴う***の形状変化を予測するための***形状変化予測方法、***形状変化予測システム、及び***形状変化予測プログラムを提供することを目的とする。 The present invention has been made in view of such problems, and provides a breast shape change prediction method, a breast shape change prediction system, and a breast shape change prediction program for predicting a breast shape change due to removal of an affected area. The purpose is to provide.
本発明は、患部摘出に伴う***の形状変化を予測するための***形状変化予測方法であって、
***部分と胸部分とを含む***モデルと、該***モデルに含まれる患部モデルとを特定する***モデル・患部モデル特定ステップと、
前記***モデルの外表面を複数の領域に区画して、前記***部分における該領域の境界に動作点を設定する一方、前記***部分と前記胸部分との境界に二以上の固定点を設定する動作点・固定点設定ステップと、
前記患部モデルの重心を算出する重心算出ステップと、
前記重心と前記固定点のそれぞれとを結ぶ二以上の固定ベクトルを算出するとともにこれらの固定ベクトルの合成ベクトルである移動方向ベクトルを算出する固定ベクトル・移動方向ベクトル算出ステップと、
前記重心を通って前記移動方向ベクトルに平行な直線である移動方向線と前記患部モデルの境界との交点を作用点と設定する作用点設定ステップと、
前記動作点と前記作用点とを直線で結んだ張力線の長さである張力線全長を算出するとともに、該張力線と前記患部モデルの境界との交点から該作用点に至る長さである張力線部分長を算出する張力線全長・張力線部分長算出ステップと、
前記張力線全長に対する前記張力線部分長の割合に応じて前記動作点を前記張力線に沿って移動させる動作点移動ステップとを備える***形状変化予測方法である。
The present invention is a breast shape change prediction method for predicting a breast shape change associated with removal of an affected area.
A breast model / affected area model identification step for identifying a breast model including a breast part and a breast part and an affected part model included in the breast model, and
The outer surface of the breast model is divided into a plurality of regions, and an operating point is set at the boundary of the region in the breast portion, while two or more fixed points are set at the boundary between the breast portion and the chest portion. Operating point / fixed point setting step and
The center of gravity calculation step for calculating the center of gravity of the affected area model and
A fixed vector / moving direction vector calculation step for calculating two or more fixed vectors connecting the center of gravity and each of the fixed points and calculating a moving direction vector which is a composite vector of these fixed vectors.
An action point setting step of setting the intersection of the movement direction line, which is a straight line passing through the center of gravity and parallel to the movement direction vector, and the boundary of the affected area model as the action point,
The total length of the tension line, which is the length of the tension line connecting the operation point and the action point with a straight line, is calculated, and the length from the intersection of the tension line and the boundary of the affected area model to the action point. Tension line total length / tension line part length calculation step to calculate tension line part length,
This is a breast shape change prediction method including an operating point moving step of moving the operating point along the tension line according to the ratio of the tension line portion length to the total length of the tension line.
上記発明にあっては、移動後の前記動作点を相互につないだ領域によって形成される外表面を、患部が摘出された後の***モデルの外表面として特定する外表面特定ステップを備えることが好ましい。 In the above invention, the outer surface specifying step for identifying the outer surface formed by the regions connecting the operating points after movement as the outer surface of the breast model after the affected part is removed may be provided. preferable.
また上記発明において、前記***モデル・患部モデル特定ステップでは、前記患部モデルを複数特定し、
それぞれの患部モデル毎に前記重心算出ステップ、前記固定ベクトル・移動方向ベクトル算出ステップ、前記作用点設定ステップ、前記張力線全長・張力線部分長算出ステップを実行し、
前記動作点移動ステップでは、それぞれの患部モデル毎に、前記張力線全長に対する前記張力線部分長の割合に応じた大きさであって前記張力線に沿う向きの張力線ベクトルを算出し、更にこれら張力線ベクトルの合成ベクトルである合成張力線ベクトルを算出して、該合成張力線ベクトルに基づいて前記動作点を移動させることが好ましい。
Further, in the above invention, in the breast model / affected area model identification step, a plurality of the affected area models are specified.
For each affected area model, the center of gravity calculation step, the fixed vector / moving direction vector calculation step, the action point setting step, and the tension line total length / tension line partial length calculation step are executed.
In the operation point movement step, for each affected portion model, a tension line vector having a size corresponding to the ratio of the tension line portion length to the total length of the tension line and in the direction along the tension line is calculated, and further, these It is preferable to calculate the composite tension line vector, which is a composite vector of the tension line vector, and move the operation point based on the composite tension line vector.
また本発明は、患部摘出に伴う***の形状変化を予測するための***形状変化予測システムであって、
***部分と胸部分とを含む***モデルと、該***モデルに含まれる患部モデルとを特定する***モデル・患部モデル特定手段と、
前記***モデルの外表面を複数の領域に区画して、前記***部分における該領域の境界に動作点を設定する一方、前記***部分と前記胸部分との境界に二以上の固定点を設定する動作点・固定点設定手段と、
前記患部モデルの重心を算出する重心算出手段と、
前記重心と前記固定点のそれぞれとを結ぶ二以上の固定ベクトルを算出するとともにこれらの固定ベクトルの合成ベクトルである移動方向ベクトルを算出する固定ベクトル・移動方向ベクトル算出手段と、
前記重心を通って前記移動方向ベクトルに平行な直線である移動方向線と前記患部モデルの境界との交点を作用点と設定する作用点設定手段と、
前記動作点と前記作用点とを直線で結んだ張力線の長さである張力線全長を算出するとともに、該張力線と前記患部モデルの境界との交点から該作用点に至る長さである張力線部分長を算出する張力線全長・張力線部分長算出手段と、
前記張力線全長に対する前記張力線部分長の割合に応じて前記動作点を前記張力線に沿って移動させる動作点移動手段とを備える***形状変化予測システムである。
Further, the present invention is a breast shape change prediction system for predicting a breast shape change accompanying removal of an affected area.
A breast model including a breast part and a chest part, a breast model for identifying the affected part model included in the breast model, a means for identifying the affected part model, and a means for identifying the affected part model.
The outer surface of the breast model is divided into a plurality of regions, and an operating point is set at the boundary of the region in the breast portion, while two or more fixed points are set at the boundary between the breast portion and the chest portion. Operating point / fixed point setting means and
A means for calculating the center of gravity of the affected area model and a means for calculating the center of gravity.
A fixed vector / moving direction vector calculating means for calculating two or more fixed vectors connecting the center of gravity and each of the fixed points and calculating a moving direction vector which is a composite vector of these fixed vectors.
An action point setting means for setting an intersection of a movement direction line, which is a straight line passing through the center of gravity and parallel to the movement direction vector, and a boundary of the affected area model as an action point.
The total length of the tension line, which is the length of the tension line connecting the operation point and the action point with a straight line, is calculated, and the length from the intersection of the tension line and the boundary of the affected area model to the action point. Tension line total length / tension line part length calculation means for calculating tension line part length,
It is a breast shape change prediction system including an operating point moving means for moving the operating point along the tension line according to the ratio of the tension line portion length to the total length of the tension line.
上記発明にあっては、移動後の前記動作点を相互につないだ領域によって形成される外表面を、患部が摘出された後の***モデルの外表面として特定する外表面特定手段を備えることが好ましい。 In the above invention, it is possible to provide an outer surface specifying means for identifying the outer surface formed by the regions connecting the operating points after movement as the outer surface of the breast model after the affected part is removed. preferable.
また本発明は、患部摘出に伴う***の形状変化を予測するための***形状変化予測プログラムであって、
コンピュータに、
***部分と胸部分とを含む***モデルと、該***モデルに含まれる患部モデルとを特定する***モデル・患部モデル特定ステップと、
前記***モデルの外表面を複数の領域に区画して、前記***部分における該領域の境界に動作点を設定する一方、前記***部分と前記胸部分との境界に二以上の固定点を設定する動作点・固定点設定ステップと、
前記患部モデルの重心を算出する重心算出ステップと、
前記重心と前記固定点のそれぞれとを結ぶ二以上の固定ベクトルを算出するとともにこれらの固定ベクトルの合成ベクトルである移動方向ベクトルを算出する固定ベクトル・移動方向ベクトル算出ステップと、
前記重心を通って前記移動方向ベクトルに平行な直線である移動方向線と前記患部モデルの境界との交点を作用点と設定する作用点設定ステップと、
前記動作点と前記作用点とを直線で結んだ張力線の長さである張力線全長を算出するとともに、該張力線と前記患部モデルの境界との交点から該作用点に至る長さである張力線部分長を算出する張力線全長・張力線部分長算出ステップと、
前記張力線全長に対する前記張力線部分長の割合に応じて前記動作点を前記張力線に沿って移動させる動作点移動ステップとを実行させる***形状変化予測プログラムである。
Further, the present invention is a breast shape change prediction program for predicting a breast shape change associated with removal of an affected area.
On the computer
A breast model / affected area model identification step for identifying a breast model including a breast part and a breast part and an affected part model included in the breast model, and
The outer surface of the breast model is divided into a plurality of regions, and an operating point is set at the boundary of the region in the breast portion, while two or more fixed points are set at the boundary between the breast portion and the chest portion. Operating point / fixed point setting step and
The center of gravity calculation step for calculating the center of gravity of the affected area model and
A fixed vector / moving direction vector calculation step for calculating two or more fixed vectors connecting the center of gravity and each of the fixed points and calculating a moving direction vector which is a composite vector of these fixed vectors.
An action point setting step of setting the intersection of the movement direction line, which is a straight line passing through the center of gravity and parallel to the movement direction vector, and the boundary of the affected area model as the action point,
The total length of the tension line, which is the length of the tension line connecting the operation point and the action point with a straight line, is calculated, and the length from the intersection of the tension line and the boundary of the affected area model to the action point. Tension line total length / tension line part length calculation step to calculate tension line part length,
This is a breast shape change prediction program that executes an operating point moving step of moving the operating point along the tension line according to the ratio of the tension line portion length to the total length of the tension line.
上記発明にあっては、移動後の前記動作点を相互につないだ領域によって形成される外表面を、患部が摘出された後の***モデルの外表面として特定する外表面特定ステップを実行させることが好ましい。 In the above invention, the outer surface identification step of identifying the outer surface formed by the regions connecting the operating points after movement as the outer surface of the breast model after the affected part is removed is executed. Is preferable.
本発明によれば、患部摘出後の***の形状を予測することができるため、手術を控えた患者の不安を取り除くことができ、また医師にとっても患者へ分かりやすく説明することができる。 According to the present invention, since the shape of the breast after removal of the affected part can be predicted, the anxiety of the patient who is about to undergo surgery can be removed, and the doctor can explain to the patient in an easy-to-understand manner.
以下、図面を参照して、本発明の一実施形態について説明する。 Hereinafter, an embodiment of the present invention will be described with reference to the drawings.
図1は、本発明に従う***形状変化予測システムの一実施形態につき、その概略構成を示すブロック図である。***形状変化予測システム1は、例えばコンピュータを用いて構成することが可能である。本実施形態の***形状変化予測システム1は、制御部2、記憶部3、入力部4、表示部5、及びこれらを接続するバス6を備えている。
FIG. 1 is a block diagram showing a schematic configuration of an embodiment of a breast shape change prediction system according to the present invention. The breast shape
制御部2は、例えばコンピュータ内部の中央演算処理装置(CPU)によって実現することができる。
The
記憶部3は、コンピュータ内部(又は外部)のハードディスクやROM又はRAMを用いて実現することができる。記憶部3では、***形状変化予測プログラムを含む各種のプログラムを記憶している。また、患部摘出前の***モデルを示すデータや患部モデルを示すデータも記憶している。***モデルや患部モデルのデータとしては、***や患部の簡易的な形状を示すデータや、複数の人の平均的な形状を示すデータ、或いは患者固有の形状を示すデータなどを用いることができる。
The
入力部4は、マウスやキーボードなどを用いて実現することができる。また、図示しないデータベースとの通信によって情報を入力する通信インターフェースを用いて実現することができる。更に、***モデルや患部モデルのデータを入力して記憶部3に記憶させるべく、CTやMRI、三次元計測器などの装置に接続するためのインターフェースを用いて実現することもできる。
The
表示部5は、モニタなどの表示装置を用いて実現することができる。また、図示しないデータベースに情報を出力するべく、上述した通信インターフェースを用いて実現することができる。更に、レーザープリンタなどの印刷装置や3Dプリンタなどの造形装置に接続するためのインターフェースを用いて実現することもできる。 The display unit 5 can be realized by using a display device such as a monitor. Further, in order to output information to a database (not shown), it can be realized by using the above-mentioned communication interface. Further, it can be realized by using an interface for connecting to a printing device such as a laser printer or a modeling device such as a 3D printer.
***形状変化予測システム1は、入力部4から取得した***モデルや患部モデルのデータを、***形状変化予測プログラムとともに記憶部3に記憶させておき、この***形状変化予測プログラムを制御部2によって実行して、***形状変化予測方法を実現するものである。以下、図1〜図9を参照しながら、***形状変化予測システム1を用いた***形状変化予測方法について説明する。
The breast shape
まず図1に示す制御部2は、記憶部3に記憶した***形状変化予測プログラム、並びに***モデル及び患部モデルのデータに基づいて、図2に示すように***部分と胸部分とを含む***モデルと、***モデルに含まれる患部モデルとを特定する(ステップS1、***モデル・患部モデル特定ステップ)。
First, the
続いて、***モデルの外表面を複数の領域に区画する。更に***部分におけるこれらの領域の境界に二以上の動作点を設定し、また***部分と胸部分との境界に二以上の固定点を設定する(ステップS2、動作点・固定点設定ステップ)。 Subsequently, the outer surface of the breast model is partitioned into multiple regions. Further, two or more operating points are set at the boundary of these regions in the breast portion, and two or more fixed points are set at the boundary between the breast portion and the chest portion (step S2, operating point / fixed point setting step).
これらステップS1、S2について、図3を参照しながらより詳細に説明する。図3において、符号11は***モデルの一例を示していて、符号12は***モデル11に含まれる患部モデルの一例を示している。なお、図3(a)、図3(b)は平面図(人が立位姿勢をとる時の正面視での図)である。また図3(c)は、図3(b)に示すA−Aに沿う断面図であり、図3(d)は、図3(b)に示すB−Bに沿う断面図である。またx軸は、立位姿勢の人を基準とした場合での左右方向、y軸は上下方向、z軸は前後方向を示す。ここで***モデル11は、半球状の***部分11aと平面状の胸部分11bとを組み合わせたものであり、患部モデル12は球状をなすものである。なお、これらのモデルを特定するにあたっては、例えば複数の人の平均的な形状を示すデータを用いてもよく、また患者固有の形状を示すデータを用いてもよい。また実際の***は、主に多数に枝分かれした乳腺と乳腺の周りを囲む脂肪で構成されているが、***モデル11及び患部モデル12は、ともに均等な中実状であるとする。
These steps S1 and S2 will be described in more detail with reference to FIG. In FIG. 3,
また図3(a)は、***モデル11の外表面を複数の領域に区画した状態を示している。本実施形態では、x軸、y軸に平行になるようにして矩形状に区画している。なおこれらの領域に区画するに当たっては、後述する患部モデル12の重心を中心に据えて行うことが好ましい。そして***部分11aにおいては、区画した領域の境界において動作点を設定する。本実施形態では、矩形状をなす領域の角部に動作点を設定している。また、***部分11aと胸部分11bとの境界には固定点を設定する。
Further, FIG. 3A shows a state in which the outer surface of the
図4は、図3(d)の拡大図を示している。ここで符号m1〜m11は、上述したように矩形状をなす領域の角部に設けられる動作点であり、符号f1〜f2は、***部分11aと胸部分11bとの境界に設けられる固定点である。なお以降の***モデル11と患部モデル12の説明は、図4に示す如き2次元断面図を用いて行う。
FIG. 4 shows an enlarged view of FIG. 3 (d). Here, the reference numerals m1 to m11 are operating points provided at the corners of the rectangular region as described above, and the reference numerals f1 to f2 are fixed points provided at the boundary between the
ステップS2を行った後は、図2に示すように、患部モデルの重心を算出する(ステップS3、重心算出ステップ)。図5に示す2次元断面図においては、円形をなす患部モデル12の中心が重心12aになる。
After performing step S2, the center of gravity of the affected area model is calculated as shown in FIG. 2 (step S3, center of gravity calculation step). In the two-dimensional cross-sectional view shown in FIG. 5, the center of the affected
その後、図2に示すように、重心と固定点のそれぞれとを結ぶ二以上の固定ベクトルを算出するとともにこれらの固定ベクトルの合成ベクトルである移動方向ベクトルを算出する(ステップS4、固定ベクトル・移動方向ベクトル算出ステップ)。本実施形態においては、図6に示すように、重心12aと固定点f1とを結ぶ固定ベクトルV1と重心12aと固定点f2とを結ぶ固定ベクトルV2とを算出することにより、これらの合成ベクトルである移動方向ベクトルV3を導くことができる。
After that, as shown in FIG. 2, two or more fixed vectors connecting the center of gravity and each of the fixed points are calculated, and a moving direction vector which is a composite vector of these fixed vectors is calculated (step S4, fixed vector / moving). Direction vector calculation step). In the present embodiment, as shown in FIG. 6, by calculating the fixed vector V1 connecting the center of
次に、図2に示すように、重心を通って移動方向ベクトルに平行な直線である移動方向線と患部モデルの境界との交点を、作用点として設定する(ステップS5、作用点設定ステップ)。図7において符号Tは、重心12aを通り、図6で説明した移動方向ベクトルV3に平行となる移動方向線を示す。また符号eは、移動方向線Tと患部モデル12の境界との交点である作用点を示す。
Next, as shown in FIG. 2, the intersection of the movement direction line, which is a straight line parallel to the movement direction vector through the center of gravity, and the boundary of the affected area model is set as the action point (step S5, action point setting step). .. In FIG. 7, reference numeral T is a moving direction line that passes through the center of
その後、図2に示すように、動作点と作用点とを直線で結んだ張力線の長さである張力線全長を算出するとともに、張力線と患部モデルの境界との交点から作用点に至る長さである張力線部分長を算出する(ステップS6、張力線全長・張力線部分長算出ステップ)。図8においては、各動作点m1〜m11と作用点eとを結ぶ張力線L1〜L11を求める。なお、各動作点m1〜m11は、患部を摘出することによってそれぞれ張力線L1〜L11上を移動すると想定しているが、患部モデル12と交差しない張力線L1〜L11上の動作点は、患部を摘出した際にもほとんど移動しないものと考えられる。このため、図8(b)に示すように患部モデル12と交差する張力線L5〜L11を採用し、張力線L5〜L11と患部モデル12の境界との交点p5〜p11を求めることとする。そして図9に示すように、それぞれの張力線L5〜L11について、各動作点m5〜m11から作用点eに至る張力線全長をそれぞれ算出するとともに、各交点p5〜p11から作用点eに至る張力線部分長をそれぞれ算出する。
After that, as shown in FIG. 2, the total length of the tension line, which is the length of the tension line connecting the operation point and the action point with a straight line, is calculated, and the point of action is reached from the intersection of the tension line and the boundary of the affected area model. The tension line partial length, which is the length, is calculated (step S6, tension line total length / tension line partial length calculation step). In FIG. 8, tension lines L1 to L11 connecting the operating points m1 to m11 and the operating point e are obtained. It is assumed that each operating point m1 to m11 moves on the tension lines L1 to L11 by removing the affected part, but the operating points on the tension lines L1 to L11 that do not intersect with the
そして、図2に示すように、張力線全長に対する張力線部分長の割合に応じて、動作点を張力線に沿って移動させる(ステップS7、動作点移動ステップ)。例えば、張力線L5における張力線全長に対する張力線部分長の割合が10%であり、張力線L9における割合が50%である場合は、動作点m9が動作点m9’へ移動する距離は、動作点m5が動作点m5’へ移動する距離の5倍になる。なお、動作点を移動させる際には、算出した各割合に対して、例えば***の大きさや年齢などを考慮した係数を掛けて全体的に補正を行ってもよい。 Then, as shown in FIG. 2, the operating point is moved along the tension line according to the ratio of the tension line portion length to the total length of the tension line (step S7, operating point moving step). For example, when the ratio of the tension line portion length to the total length of the tension line in the tension line L5 is 10% and the ratio in the tension line L9 is 50%, the distance that the operating point m9 moves to the operating point m9'is the operation. The point m5 is five times the distance moved to the operating point m5'. When moving the operating point, each calculated ratio may be corrected as a whole by multiplying each calculated ratio by a coefficient that considers, for example, the size and age of the breast.
本実施形態では更に、図2に示すように、移動後の動作点を相互につないだ領域によって形成される外表面を、患部が摘出された後の***モデルの外表面として特定する(ステップS8、外表面特定ステップ)。図10では、移動後の動作点m5’〜 m11’を相互につないだ線を描くことによって、患部が摘出された後の***モデルの外表面として特定している。 In the present embodiment, as shown in FIG. 2, the outer surface formed by the regions connecting the operating points after movement is further specified as the outer surface of the breast model after the affected part is removed (step S8). , Outer surface identification step). In FIG. 10, by drawing a line connecting the operating points m5'to m11' after movement to each other, the affected part is identified as the outer surface of the breast model after being removed.
その後は、図10のようにして描いた図を、図1に示した表示部5に表示させる。このように本実施形態では、移動後の動作点、及び移動後の動作点を相互につないだ線を表示させることができるため、患部を摘出した後の***の形状を患者に対して視覚で認識させることができる。 After that, the figure drawn as shown in FIG. 10 is displayed on the display unit 5 shown in FIG. As described above, in the present embodiment, since the operating point after the movement and the line connecting the operating points after the movement can be displayed, the shape of the breast after the affected part is removed can be visually displayed to the patient. Can be recognized.
なお、これまでに説明したステップは、外表面を複数の領域に区画した***モデル11のうち、図3(b)に示すB−Bに沿う領域での2次元断面図についてのものであるが、これを***モデル11の全ての領域に展開することによって、患部を摘出した後の***の形状を3次元的に表示することができる。
The steps described so far relate to a two-dimensional cross-sectional view of the
図11〜図17は、3次元の***モデルを用いて上述したステップS1〜S8を示した図である。なお、図11に示すように本実施形態の***部分11aは半球状であって、胸部分11bは平面状である。また患部モデル12は円柱状である。そしてX軸、Y軸、Z軸は、先に説明したx軸、y軸、z軸と同様に、立位姿勢の人を基準とした場合において、それぞれ左右方向、上下方向、前後方向を示している。
11 to 17 are views showing steps S1 to S8 described above using a three-dimensional breast model. As shown in FIG. 11, the
本実施形態においても***形状変化予測システム1は、図11〜図13に示すように***部分11aと胸部分11bとを含む***モデルと、***モデルに含まれる患部モデル12とを特定する(ステップS1、***モデル・患部モデル特定ステップ)。そして***モデルの外表面を複数の領域に区画し(本実施形態では三角形状)、更に***部分11aにおけるこれらの領域の境界に二以上の動作点mを設定し(図15参照)、***部分11aと胸部分11bとの境界に二以上の固定点fを設定する(図13、図14参照。ステップS2、動作点・固定点設定ステップ)。また、図12に示すように患部モデル12の重心12aを算出し(ステップS3、重心算出ステップ)、重心12aと固定点fのそれぞれとを結ぶ二以上の固定ベクトルVを算出し(図14参照)、更にこれらの固定ベクトルVの合成ベクトルである移動方向ベクトルを算出し(不図示。ステップS4、固定ベクトル・移動方向ベクトル算出ステップ)、図16に示すように、重心12aを通って移動方向ベクトルに平行な直線である移動方向線Tと患部モデル12の境界との交点を、作用点eとして設定する(ステップS5、作用点設定ステップ)。
Also in this embodiment, the breast shape
その後は、動作点mと作用点eとを直線で結んだ張力線の長さである張力線全長を算出するとともに、張力線と患部モデル12の境界との交点から作用点eに至る長さである張力線部分長を算出し(不図示。ステップS6、張力線全長・張力線部分長算出ステップ)、張力線全長に対する張力線部分長の割合に応じて、動作点mを張力線に沿って移動させ(ステップS7、動作点移動ステップ)、移動後の動作点を相互につないだ領域によって形成される外表面を、患部が摘出された後の***モデルの外表面として特定する(図17参照。ステップS8、外表面特定ステップ)。このようにして、患部摘出に伴う***の形状変化を3次元的にも示すことができる。
After that, the total length of the tension line, which is the length of the tension line connecting the operation point m and the action point e with a straight line, is calculated, and the length from the intersection of the tension line and the boundary of the affected
そして、例えば三次元計測器によって患者の***部分及び胸部分を計測し、これを***モデルに使用することも可能である。図18(a)及び図19(a)は、三次元計測器を用いて得られた***モデルに関し、患部を摘出する前の状態を示す図であって、図18(b)及び図19(b)は、患部を摘出した後の状態を示す図である。 Then, for example, it is possible to measure the breast portion and the breast portion of the patient with a three-dimensional measuring instrument and use this in the breast model. 18 (a) and 19 (a) are views showing the state of the breast model obtained by using the three-dimensional measuring instrument before the affected part is removed, and are shown in FIGS. 18 (b) and 19 (a). b) is a figure which shows the state after removing the affected part.
ここで、***形状変化予測システム1で得られる***の形状と実際に患部を摘出した***の形状との比較を行ったところ、高い相関関係が認められた。すなわち本発明によれば、患部を摘出した後の***の形状を高い精度で予測することが可能である。
Here, when the shape of the breast obtained by the breast shape
上述した実施形態では、患部モデルを1つ特定して***の形状変化を予測したが、患部モデルを複数特定して実際に取り出す患部の形状に近づけることによって、更に高い精度で***の形状変化を予測することができる。以下、前述の図4〜図8、図10に対応する図を用いて、患部モデルを複数特定した場合について説明する。 In the above-described embodiment, one affected part model is specified to predict the shape change of the breast, but by specifying a plurality of affected part models and bringing them closer to the shape of the affected part to be actually taken out, the shape change of the breast can be predicted with higher accuracy. Can be predicted. Hereinafter, a case where a plurality of affected area models are specified will be described with reference to the drawings corresponding to FIGS. 4 to 8 and 10 described above.
図20に示すように本実施形態では、図4で示した実施形態と同様にして***モデルを特定するとともに、第1患部モデルAと第2患部モデルBの2つの患部モデルを特定するものとする(ステップS1、***モデル・患部モデル特定ステップ)。また、図4で示した実施形態と同様にして動作点m1〜m11と固定点f1〜f2を設定するものとする(ステップS2、動作点・固定点設定ステップ)。 As shown in FIG. 20, in the present embodiment, the breast model is specified in the same manner as in the embodiment shown in FIG. 4, and two affected part models, the first affected part model A and the second affected part model B, are specified. (Step S1, breast model / affected area model identification step). Further, it is assumed that the operating points m1 to m11 and the fixed points f1 to f2 are set in the same manner as in the embodiment shown in FIG. 4 (step S2, operating point / fixed point setting step).
続いて、図21に示すように、第1患部モデルAの重心Aaと第2患部モデルBの重心Baを算出し(ステップS3、重心算出ステップ)、更に、第1患部モデルA、第2患部モデルB毎に、固定ベクトルと移動方向ベクトルを算出する(ステップS4、固定ベクトル・移動方向ベクトル算出ステップ)。具体的には図22に示すように、第1患部モデルAの重心Aaについては、重心Aaと固定点f1とを結ぶ固定ベクトルVA1と重心Aaと固定点f2とを結ぶ固定ベクトルVA2とを算出することにより、これらの合成ベクトルである移動方向ベクトルVA3を導くこととする。同様にして第2患部モデルBの重心Baについても、固定ベクトルVB1と固定ベクトルVB2から移動方向ベクトルVB3を導くこととする。 Subsequently, as shown in FIG. 21, the center of gravity Aa of the first affected area model A and the center of gravity Ba of the second affected area model B are calculated (step S3, center of gravity calculation step), and further, the first affected area model A and the second affected area are calculated. A fixed vector and a moving direction vector are calculated for each model B (step S4, fixed vector / moving direction vector calculation step). Specifically, as shown in FIG. 22, for the center of gravity Aa of the first affected area model A, a fixed vector VA1 connecting the center of gravity Aa and the fixed point f1 and a fixed vector VA2 connecting the center of gravity Aa and the fixed point f2 are calculated. By doing so, the moving direction vector VA3, which is a composite vector of these, is derived. Similarly, for the center of gravity Ba of the second affected area model B, the moving direction vector VB3 is derived from the fixed vector VB1 and the fixed vector VB2.
次に、第1患部モデルA、第2患部モデルB毎に、それぞれの重心を通って移動方向ベクトルに平行な直線である移動方向線と患部モデルの境界との交点を、作用点として設定する(ステップS5、作用点設定ステップ)。図23に示すように本実施形態においては、第1患部モデルAについては移動方向線TAと作用点eAを設定し、第2患部モデルBについては移動方向線TBと作用点eBを設定するものとする。 Next, for each of the first affected part model A and the second affected part model B, the intersection of the moving direction line, which is a straight line parallel to the moving direction vector through the center of gravity, and the boundary of the affected part model is set as the action point. (Step S5, step of setting the point of action). As shown in FIG. 23, in the present embodiment, the movement direction line TA and the action point eA are set for the first affected part model A, and the movement direction line TB and the action point eB are set for the second affected part model B. And.
その後、第1患部モデルA、第2患部モデルB毎に、動作点と作用点とを直線で結んだ張力線の長さである張力線全長を算出するとともに、張力線と患部モデルの境界との交点から作用点に至る長さである張力線部分長を算出する(ステップS6、張力線全長・張力線部分長算出ステップ)。本実施形態においては、図24に示すように、第1患部モデルAについては、各動作点m1〜m11と作用点eAとを結ぶ張力線LA1〜LA11を求め、第1患部モデルAと交差する張力線上の動作点が患部を摘出した際に実質的に移動するものであるとして、図25(a)に示すように張力線LA5〜LA11と第1患部モデルAの境界との交点pA5〜pA11を求めることとする。そして、図9を参照しながら説明した前述の実施形態のように、それぞれの張力線LA5〜LA11について、各動作点m5〜m11から作用点eAに至る張力線全長をそれぞれ算出するとともに、各交点pA5〜pA11から作用点eAに至る張力線部分長をそれぞれ算出する。同様にして、第2患部モデルBについては、図24に示すように各動作点m1〜m11と作用点eBとを結ぶ張力線LB1〜LB11を求め、更に、図25(b)に示すように張力線LB1〜LB10と第2患部モデルBの境界との交点pB1〜pB10を求めた後、張力線LB1〜LB10について、張力線全長と張力線部分長をそれぞれ算出する。 After that, for each of the first affected part model A and the second affected part model B, the total length of the tension line, which is the length of the tension line connecting the operation point and the action point with a straight line, is calculated, and the boundary between the tension line and the affected part model is calculated. The tension line partial length, which is the length from the intersection of the above to the point of action, is calculated (step S6, the tension line total length / tension line partial length calculation step). In the present embodiment, as shown in FIG. 24, for the first affected part model A, the tension lines LA1 to LA11 connecting the operating points m1 to m11 and the operating point eA are obtained and intersect with the first affected part model A. Assuming that the operating point on the tension line substantially moves when the affected part is removed, as shown in FIG. 25 (a), the intersection point pA5 to pA11 between the tension line LA5 to LA11 and the boundary of the first affected part model A Will be sought. Then, as in the above-described embodiment described with reference to FIG. 9, for each tension line LA5 to LA11, the total length of the tension line from each operating point m5 to m11 to the point of action eA is calculated, and each intersection point. The tension line partial lengths from pA5 to pA11 to the operating point eA are calculated respectively. Similarly, for the second affected area model B, the tension lines LB1 to LB11 connecting the operating points m1 to m11 and the operating point eB are obtained as shown in FIG. 24, and further, as shown in FIG. 25 (b). After obtaining the intersections pB1 to pB10 between the tension lines LB1 to LB10 and the boundary of the second affected part model B, the total length of the tension line and the partial length of the tension line are calculated for the tension lines LB1 to LB10, respectively.
そして、第1患部モデルA、第2患部モデルB毎に、張力線全長に対する張力線部分長の割合に応じた大きさであって張力線に沿う向きの張力線ベクトルを算出し、更にこれら張力線ベクトルの合成ベクトルである合成張力線ベクトルを算出して、合成張力線ベクトルに基づいて動作点を移動させる(ステップS7、動作点移動ステップ)。例えば、動作点m7に対し、第1患部モデルAについては張力線LA7における張力線全長に対する張力線部分長の割合が30%であって、第2患部モデルBについては、張力線LB7における張力線全長に対する張力線部分長の割合が60%であったとする。この場合は図26に示すように、第1患部モデルAについては、張力線全長に対する張力線部分長の割合が30%分の大きさであって張力線LA7に沿う向きの張力線ベクトルTVA7を算出し、第2患部モデルBについては、張力線全長に対する張力線部分長の割合が60%分の大きさであって張力線LB7に沿う向きの張力線ベクトルTVB7を算出する。そして、張力線ベクトルTVA7と張力線ベクトルTVB7の合成ベクトルである合成張力線ベクトルTV7を算出して、合成張力線ベクトルTV7に基づいて動作点m7をm7’へ移動させる。なお、合成張力線ベクトルに基づいて動作点を移動させる際には、例えば***の大きさや年齢などを考慮した係数を掛けて全体的に補正を行ってもよいし、患部モデルの数に従う平均化(本実施形態では患部モデルは2つなので合成張力線ベクトルの長さを1/2にする)を行ってもよい。このようにして、その他の動作点についても合成張力線ベクトルに基づいて移動させるものとする。 Then, for each of the first affected part model A and the second affected part model B, a tension line vector having a size corresponding to the ratio of the tension line portion length to the total length of the tension line and in the direction along the tension line is calculated, and further, these tensions are obtained. A composite tension line vector, which is a composite vector of the line vectors, is calculated, and the operation point is moved based on the composite tension line vector (step S7, operation point movement step). For example, the ratio of the tension line portion length to the total length of the tension line in the tension line LA7 is 30% with respect to the operating point m7 for the first affected part model A, and the tension line in the tension line LB7 for the second affected part model B. It is assumed that the ratio of the tension line portion length to the total length is 60%. In this case, as shown in FIG. 26, for the first affected area model A, the tension line vector TVA7 in which the ratio of the tension line portion length to the total length of the tension line is 30% and is oriented along the tension line LA7 is provided. For the second affected area model B, the tension line vector TVB7 in the direction along the tension line LB7, in which the ratio of the tension line portion length to the total length of the tension line is 60%, is calculated. Then, the combined tension line vector TV7, which is a combined vector of the tension line vector TVA7 and the tension line vector TVB7, is calculated, and the operating point m7 is moved to m7'based on the combined tension line vector TV7. When moving the operating point based on the synthetic tension line vector, for example, it may be corrected by multiplying it by a coefficient considering the size and age of the breast, or it may be averaged according to the number of affected area models. (Since there are two affected area models in this embodiment, the length of the synthetic tension line vector is halved). In this way, the other operating points are also moved based on the synthetic tension line vector.
その後は、移動後の動作点を相互につないだ領域によって形成される外表面を、患部が摘出された後の***モデルの外表面として特定する(ステップS8、外表面特定ステップ)。本実施形態では、図27に示すように、移動後の動作点m1’〜 m11’を相互につないだ線を描くことによって、患部が摘出された後の***モデルの外表面として特定している。そして、図27のようにして描いた図を、図1に示した表示部5に表示させることによって、患部を摘出した後の***の形状を患者に対して視覚で認識させることができる。 After that, the outer surface formed by the regions connecting the operating points after movement is specified as the outer surface of the breast model after the affected part is removed (step S8, outer surface identification step). In the present embodiment, as shown in FIG. 27, by drawing a line connecting the operating points m1'to m11' after movement to each other, the affected part is specified as the outer surface of the breast model after being removed. .. Then, by displaying the figure drawn as shown in FIG. 27 on the display unit 5 shown in FIG. 1, the shape of the breast after the affected part is removed can be visually recognized by the patient.
なお、患部モデルを複数設定するにあたっては、上述した実施形態のように2つに限られるものではなく、更に数を増やしてもよい。また、本実施形態では複数の患部モデルを設定した場合について2次元のモデルで説明したが、3次元のモデルでも適用可能である。ところで、本願発明者が様々な***モデルと患部モデルを用いて検討を重ねたところ、3次元モデルに対して患部モデルを複数設定するにあたっては、それぞれが球状になる患部モデルを数珠つなぎに3つ連ねて実際の患部を模した場合は、***形状変化予測システム1で使用するコンピュータにおいてそれ程高い性能のものを必要とすることなく、患部を摘出した後の***の形状を特に高い精度で予測できることを見出した。なお、患部モデルを数珠つなぎに3つ連ねるにあたっては、図28(a)のように同一形状の患部モデルA、B、Cが直列状に連なるものに限られず、図28(b)のように、患部モデルA、B、Cの半径がそれぞれ異なっていてもよく、また、途中で屈曲するように連なっていてもよい。
When setting a plurality of affected area models, the number is not limited to two as in the above-described embodiment, and the number may be further increased. Further, in the present embodiment, the case where a plurality of affected area models are set has been described with a two-dimensional model, but a three-dimensional model can also be applied. By the way, when the inventor of the present application repeated studies using various breast models and affected area models, when setting a plurality of affected area models for a three-dimensional model, three affected area models, each of which is spherical, are connected in a row. When the actual affected area is imitated in a row, the shape of the breast after the affected area is removed can be predicted with particularly high accuracy without requiring a computer used in the breast shape
以上、具体例を挙げて本発明の実施形態を説明したが、本発明の特許請求の範囲から逸脱しない限りにおいて、あらゆる変形や変更が可能であることは当業者に明らかである。例えば、***モデルの外表面を複数の領域に区画するにあたっては、外表面に対して平行になるように分割して基本的に領域の面積が全て等しくなるように設定してもよい。また領域の形状は矩形状や三角形状に限られず、他の多角形状であってもよい。 Although the embodiments of the present invention have been described above with specific examples, it is clear to those skilled in the art that any modification or modification is possible as long as the invention does not deviate from the claims of the present invention. For example, when the outer surface of the breast model is divided into a plurality of regions, the outer surface may be divided so as to be parallel to the outer surface and set so that the areas of the regions are basically the same. Further, the shape of the region is not limited to a rectangular shape or a triangular shape, and may be another polygonal shape.
1:***形状変化予測システム
2:制御部
3:記憶部
4:入力部
5:表示部
6:バス
11:***モデル
11a:***部分
11b:胸部分
12:患部モデル
12a:重心
A:第1患部モデル
Aa:第1患部モデルの重心
B:第2患部モデル
Ba:第2患部モデルの重心
L1〜L11:張力線
T:移動方向線
V1、V2、V:固定ベクトル
V3:移動方向ベクトル
e:作用点
f1、f2、f:固定点
m1〜m11、m:動作点
1: Breast shape change prediction system 2: Control unit 3: Storage unit 4: Input unit 5: Display unit 6: Bus 11:
Claims (4)
***部分と胸部分とを含む***モデルと、該***モデルに含まれる患部モデルとを特定する***モデル・患部モデル特定ステップと、
前記***モデルの外表面を複数の領域に区画して、前記***部分における該領域の境界に動作点を設定する一方、前記***部分と前記胸部分との境界に二以上の固定点を設定する動作点・固定点設定ステップと、
前記患部モデルの重心を算出する重心算出ステップと、
前記重心と前記固定点のそれぞれとを結ぶ二以上の固定ベクトルを算出するとともにこれらの固定ベクトルの合成ベクトルである移動方向ベクトルを算出する固定ベクトル・移動方向ベクトル算出ステップと、
前記重心を通って前記移動方向ベクトルに平行な直線である移動方向線と前記患部モデルの境界との交点を作用点と設定する作用点設定ステップと、
前記動作点と前記作用点とを直線で結んだ張力線の長さである張力線全長を算出するとともに、該張力線と前記患部モデルの境界との交点から該作用点に至る長さである張力線部分長を算出する張力線全長・張力線部分長算出ステップと、
前記張力線全長に対する前記張力線部分長の割合に応じて前記動作点を前記張力線に沿って移動させる動作点移動ステップと、
移動後の前記動作点を相互につないだ領域によって形成される外表面を、患部が摘出された後の***モデルの外表面として特定する外表面特定ステップと、を備える***形状変化予測方法。 It is a breast shape change prediction method for predicting the shape change of the breast due to the removal of the affected area.
A breast model / affected area model identification step for identifying a breast model including a breast part and a breast part and an affected part model included in the breast model, and
The outer surface of the breast model is divided into a plurality of regions, and an operating point is set at the boundary of the region in the breast portion, while two or more fixed points are set at the boundary between the breast portion and the chest portion. Operating point / fixed point setting step and
The center of gravity calculation step for calculating the center of gravity of the affected area model and
A fixed vector / moving direction vector calculation step for calculating two or more fixed vectors connecting the center of gravity and each of the fixed points and calculating a moving direction vector which is a composite vector of these fixed vectors.
An action point setting step of setting the intersection of the movement direction line, which is a straight line passing through the center of gravity and parallel to the movement direction vector, and the boundary of the affected area model as the action point,
The total length of the tension line, which is the length of the tension line connecting the operation point and the action point with a straight line, is calculated, and the length from the intersection of the tension line and the boundary of the affected area model to the action point. Tension line total length / tension line part length calculation step to calculate tension line part length,
An operating point moving step of moving the operating point along the tension line according to the ratio of the tension line portion length to the total length of the tension line .
A method for predicting a change in breast shape , comprising an outer surface specifying step for identifying an outer surface formed by regions connecting the operating points after movement as an outer surface of a breast model after the affected part is removed.
それぞれの患部モデル毎に前記重心算出ステップ、前記固定ベクトル・移動方向ベクトル算出ステップ、前記作用点設定ステップ、前記張力線全長・張力線部分長算出ステップを実行し、
前記動作点移動ステップでは、それぞれの患部モデル毎に、前記張力線全長に対する前記張力線部分長の割合に応じた大きさであって前記張力線に沿う向きの張力線ベクトルを算出し、更にこれら張力線ベクトルの合成ベクトルである合成張力線ベクトルを算出して、該合成張力線ベクトルに基づいて前記動作点を移動させる請求項1に記載の***形状変化予測方法。 In the breast model / affected area model identification step, a plurality of the affected area models are specified.
For each affected area model, the center of gravity calculation step, the fixed vector / moving direction vector calculation step, the action point setting step, and the tension line total length / tension line partial length calculation step are executed.
In the operation point movement step, for each affected portion model, a tension line vector having a size corresponding to the ratio of the tension line portion length to the total length of the tension line and in the direction along the tension line is calculated, and further, these and calculates the combined tension line vector is a composite vector of the line of tension vector, breast shape change predicting method according to claim 1 to move the operating point based on the composite tension line vector.
***部分と胸部分とを含む***モデルと、該***モデルに含まれる患部モデルとを特定する***モデル・患部モデル特定手段と、
前記***モデルの外表面を複数の領域に区画して、前記***部分における該領域の境界に動作点を設定する一方、前記***部分と前記胸部分との境界に二以上の固定点を設定する動作点・固定点設定手段と、
前記患部モデルの重心を算出する重心算出手段と、
前記重心と前記固定点のそれぞれとを結ぶ二以上の固定ベクトルを算出するとともにこれらの固定ベクトルの合成ベクトルである移動方向ベクトルを算出する固定ベクトル・移動方向ベクトル算出手段と、
前記重心を通って前記移動方向ベクトルに平行な直線である移動方向線と前記患部モデルの境界との交点を作用点と設定する作用点設定手段と、
前記動作点と前記作用点とを直線で結んだ張力線の長さである張力線全長を算出するとともに、該張力線と前記患部モデルの境界との交点から該作用点に至る長さである張力線部分長を算出する張力線全長・張力線部分長算出手段と、
前記張力線全長に対する前記張力線部分長の割合に応じて前記動作点を前記張力線に沿って移動させる動作点移動手段と、
移動後の前記動作点を相互につないだ領域によって形成される外表面を、患部が摘出された後の***モデルの外表面として特定する外表面特定手段と、を備える***形状変化予測システム。 A breast shape change prediction system for predicting breast shape changes associated with removal of the affected area.
A breast model including a breast part and a chest part, a breast model for identifying the affected part model included in the breast model, and a means for identifying the affected part model,
The outer surface of the breast model is divided into a plurality of regions, and an operating point is set at the boundary of the region in the breast portion, while two or more fixed points are set at the boundary between the breast portion and the chest portion. Operating point / fixed point setting means and
A means for calculating the center of gravity of the affected area model and a means for calculating the center of gravity.
A fixed vector / moving direction vector calculating means for calculating two or more fixed vectors connecting the center of gravity and each of the fixed points and calculating a moving direction vector which is a composite vector of these fixed vectors.
An action point setting means for setting an intersection of a movement direction line, which is a straight line passing through the center of gravity and parallel to the movement direction vector, and a boundary of the affected area model as an action point.
The total length of the tension line, which is the length of the tension line connecting the operation point and the action point with a straight line, is calculated, and the length from the intersection of the tension line and the boundary of the affected area model to the action point. Tension line total length / tension line part length calculation means for calculating tension line part length,
An operating point moving means for moving the operating point along the tension line according to the ratio of the tension line portion length to the total length of the tension line .
A breast shape change prediction system including an outer surface identifying means for identifying an outer surface formed by regions connecting the operating points after movement as an outer surface of a breast model after the affected part is removed.
コンピュータに、
***部分と胸部分とを含む***モデルと、該***モデルに含まれる患部モデルとを特定する***モデル・患部モデル特定ステップと、
前記***モデルの外表面を複数の領域に区画して、前記***部分における該領域の境界に動作点を設定する一方、前記***部分と前記胸部分との境界に二以上の固定点を設定する動作点・固定点設定ステップと、
前記患部モデルの重心を算出する重心算出ステップと、
前記重心と前記固定点のそれぞれとを結ぶ二以上の固定ベクトルを算出するとともにこれらの固定ベクトルの合成ベクトルである移動方向ベクトルを算出する固定ベクトル・移動方向ベクトル算出ステップと、
前記重心を通って前記移動方向ベクトルに平行な直線である移動方向線と前記患部モデルの境界との交点を作用点と設定する作用点設定ステップと、
前記動作点と前記作用点とを直線で結んだ張力線の長さである張力線全長を算出するとともに、該張力線と前記患部モデルの境界との交点から該作用点に至る長さである張力線部分長を算出する張力線全長・張力線部分長算出ステップと、
前記張力線全長に対する前記張力線部分長の割合に応じて前記動作点を前記張力線に沿って移動させる動作点移動ステップと、
移動後の前記動作点を相互につないだ領域によって形成される外表面を、患部が摘出された後の***モデルの外表面として特定する外表面特定ステップと、を実行させる***形状変化予測プログラム。 A breast shape change prediction program for predicting breast shape changes associated with removal of the affected area.
On the computer
A breast model / affected area model identification step for identifying a breast model including a breast part and a breast part and an affected part model included in the breast model, and
The outer surface of the breast model is divided into a plurality of regions, and an operating point is set at the boundary of the region in the breast portion, while two or more fixed points are set at the boundary between the breast portion and the chest portion. Operating point / fixed point setting step and
The center of gravity calculation step for calculating the center of gravity of the affected area model and
A fixed vector / moving direction vector calculation step for calculating two or more fixed vectors connecting the center of gravity and each of the fixed points and calculating a moving direction vector which is a composite vector of these fixed vectors.
An action point setting step of setting the intersection of the movement direction line, which is a straight line passing through the center of gravity and parallel to the movement direction vector, and the boundary of the affected area model as the action point,
The total length of the tension line, which is the length of the tension line connecting the operation point and the action point with a straight line, is calculated, and the length from the intersection of the tension line and the boundary of the affected area model to the action point. Tension line total length / tension line part length calculation step to calculate tension line part length,
An operating point moving step of moving the operating point along the tension line according to the ratio of the tension line portion length to the total length of the tension line .
A breast shape change prediction program that executes an outer surface identification step of identifying the outer surface formed by the regions connecting the operating points after movement as the outer surface of the breast model after the affected part is removed.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2016163193 | 2016-08-24 | ||
JP2016163193 | 2016-08-24 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2018033951A JP2018033951A (en) | 2018-03-08 |
JP6937187B2 true JP6937187B2 (en) | 2021-09-22 |
Family
ID=61564604
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2017157284A Active JP6937187B2 (en) | 2016-08-24 | 2017-08-16 | Breast shape change prediction method, breast shape change prediction system, and breast shape change prediction program |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6937187B2 (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR102329034B1 (en) * | 2019-09-10 | 2021-11-19 | 고려대학교 산학협력단 | Method for predicting simulation of breast deformity and apparatus for executing the method |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6119033A (en) * | 1997-03-04 | 2000-09-12 | Biotrack, Inc. | Method of monitoring a location of an area of interest within a patient during a medical procedure |
US11020184B2 (en) * | 2013-09-24 | 2021-06-01 | Koninklijke Philips N.V. | System and method for accurately determining and marking skin incision lines for surgical intervention |
-
2017
- 2017-08-16 JP JP2017157284A patent/JP6937187B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
JP2018033951A (en) | 2018-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
ES2704080T3 (en) | Biometric simulation device, method to control the biometric simulation device, and program to control the biometric simulation device | |
JP4521271B2 (en) | An image processing system for displaying information on the movement of the top of a deformable three-dimensional object | |
Lacko et al. | Evaluation of an anthropometric shape model of the human scalp | |
US9808213B2 (en) | Image processing apparatus, image processing method, medical image diagnostic system, and storage medium | |
KR101693808B1 (en) | Method and apparatus 3d surgery simulation of nasal cavity | |
US20170236279A1 (en) | Image analyzing device, image analyzing method, and computer program product | |
JP2015080720A (en) | Apparatus and method for computer-aided diagnosis | |
US20150073765A1 (en) | System and method for prediction of respiratory motion from 3d thoracic images | |
Wallman et al. | Computational methods to reduce uncertainty in the estimation of cardiac conduction properties from electroanatomical recordings | |
CN106898027A (en) | Mapping of the 3-D view to two dimensional image | |
US20180005378A1 (en) | Atlas-Based Determination of Tumor Growth Direction | |
WO2015010745A1 (en) | Multi-modal segmentation of image data | |
Hachaj et al. | Visualization of perfusion abnormalities with GPU-based volume rendering | |
Wang et al. | Real time 3D simulation for nose surgery and automatic individual prosthesis design | |
US20220108540A1 (en) | Devices, systems and methods for generating and providing image information | |
Wei et al. | A (near) real-time simulation method of aneurysm coil embolization | |
JP6937187B2 (en) | Breast shape change prediction method, breast shape change prediction system, and breast shape change prediction program | |
JP2007312837A (en) | Region extracting apparatus, region extracting method and program | |
JP5721225B2 (en) | Shape data generation method, program, and apparatus | |
JP4910478B2 (en) | Modeling apparatus, modeling method and program | |
Segars et al. | Improved dynamic cardiac phantom based on 4D NURBS and tagged MRI | |
US20180315181A1 (en) | Deformable registration of magnetic resonance and ultrasound images using biomechanical models | |
de Bruin et al. | Interactive 3D segmentation using connected orthogonal contours | |
Quammen et al. | The virtual pediatric airways workbench | |
Wang et al. | Using optimal transport to improve spherical harmonic quantification of complex biological shapes |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AA64 | Notification of invalidation of claim of internal priority (with term) |
Free format text: JAPANESE INTERMEDIATE CODE: A241764 Effective date: 20170828 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20170914 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20200721 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20210528 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20210601 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20210628 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20210824 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20210830 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6937187 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |