JP6937187B2 - Breast shape change prediction method, breast shape change prediction system, and breast shape change prediction program - Google Patents

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.

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.

Japanese Unexamined Patent Publication No. 2004-209019

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.

It is a block diagram which shows one Embodiment of the breast shape change prediction system according to this invention. It is a flow chart which shows one Embodiment of the breast shape change prediction method according to this invention. It is a figure which shows an example of a breast model and an affected part model. It is an enlarged view of FIG. 3C, and is the figure which shows the operating point and the fixed point. It is a figure which shows the center of gravity of the affected part model. It is a figure which shows the fixed vector and the moving direction vector. It is a figure which shows the movement direction line and the point of action. (A) is a diagram showing a tension line and an intersection of the tension line and the affected area model, and (b) is an enlarged view thereof. It is a figure explaining the tension line total length and tension line part length. It is a figure which shows the operating point after movement and the outer surface of a breast model. It is a figure which shows the breast model three-dimensionally. It is a figure which shows the affected part model three-dimensionally. It is a figure which shows the breast model and the affected part model shown in FIG. 11 and FIG. It is a figure which shows the fixed vector which connects the center of gravity and a fixed point. It is a figure which shows the operating point in a breast part. It is a figure which shows the movement direction line and the point of action. It is a figure which shows the outer surface of the breast model after the affected part was removed. It is a figure explaining the embodiment using one breast model obtained by measuring the breast part and the chest part of a patient. It is a figure explaining the embodiment using another breast model obtained by measuring the breast part and the chest part of a patient. It is a figure which shows the operating point and the fixed point in the state which set a plurality of affected part models. It is a figure which shows the center of gravity of the affected part model in the state which set a plurality of affected part models. It is a figure which shows the fixed vector and the moving direction vector in the state which set a plurality of affected part models. It is a figure which shows the movement direction line and the point of action in the state which set a plurality of affected part models. It is a figure which shows the tension line in the state which set a plurality of affected part models, and the intersection of a tension line and an affected part model. In the enlarged view of FIG. 24, FIG. 24A is a diagram showing the periphery of the first affected area model, and FIG. 24B is a diagram showing the periphery of the second affected area model. It is a figure which showed the tension line vector and the composite tension line vector. It is a figure which shows the operating point after movement and the outer surface of a breast model in the state which set a plurality of affected part models. It is a figure which shows a more preferable state in setting a plurality of affected part models.

Hereinafter, an embodiment of the present invention will be described with reference to the drawings.

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 change prediction system 1 can be configured by using, for example, a computer. The breast shape change prediction system 1 of the present embodiment includes a control unit 2, a storage unit 3, an input unit 4, a display unit 5, and a bus 6 connecting them.

The control unit 2 can be realized by, for example, a central processing unit (CPU) inside a computer.

The storage unit 3 can be realized by using a hard disk, ROM, or RAM inside (or outside) the computer. The storage unit 3 stores various programs including a breast shape change prediction program. It also stores data showing the breast model before removal of the affected area and data showing the affected area model. As the data of the breast model or the affected part model, data showing a simple shape of the breast or the affected part, data showing the average shape of a plurality of people, data showing a shape peculiar to a patient, or the like can be used.

The input unit 4 can be realized by using a mouse, a keyboard, or the like. Further, it can be realized by using a communication interface for inputting information by communicating with a database (not shown). Further, in order to input the data of the breast model or the affected part model and store it in the storage unit 3, it can be realized by using an interface for connecting to a device such as CT, MRI, or a three-dimensional measuring instrument.

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.

The breast shape change prediction system 1 stores the data of the breast model and the affected area model acquired from the input unit 4 in the storage unit 3 together with the breast shape change prediction program, and executes this breast shape change prediction program by the control unit 2. Then, the breast shape change prediction method is realized. Hereinafter, a breast shape change prediction method using the breast shape change prediction system 1 will be described with reference to FIGS. 1 to 9.

First, the control unit 2 shown in FIG. 1 is a breast model including a breast part and a breast part as shown in FIG. 2 based on the breast shape change prediction program stored in the storage unit 3 and the data of the breast model and the affected part model. And the affected part model included in the breast model (step S1, breast model / affected part model identification step).

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).

These steps S1 and S2 will be described in more detail with reference to FIG. In FIG. 3, reference numeral 11 indicates an example of the breast model, and reference numeral 12 indicates an example of the affected part model included in the breast model 11. 3 (a) and 3 (b) are plan views (views taken from the front when a person takes a standing posture). 3 (c) is a cross-sectional view taken along the line AA shown in FIG. 3 (b), and FIG. 3 (d) is a cross-sectional view taken along the line BB shown in FIG. 3 (b). The x-axis indicates the left-right direction when the person in the standing posture is used as a reference, the y-axis indicates the up-down direction, and the z-axis indicates the front-back direction. Here, the breast model 11 is a combination of a hemispherical breast portion 11a and a flat breast portion 11b, and the affected portion model 12 has a spherical shape. In specifying these models, for example, data showing the average shape of a plurality of people may be used, or data showing a patient-specific shape may be used. Further, the actual breast is mainly composed of a large number of branched mammary glands and fat surrounding the mammary glands, but both the breast model 11 and the affected area model 12 are assumed to have an even solid state.

Further, FIG. 3A shows a state in which the outer surface of the breast model 11 is divided into a plurality of regions. In the present embodiment, the sections are rectangularly divided so as to be parallel to the x-axis and the y-axis. In dividing into these areas, it is preferable to set the center of gravity of the affected area model 12, which will be described later, at the center. Then, in the breast portion 11a, an operating point is set at the boundary of the partitioned region. In the present embodiment, the operating point is set at the corner of the rectangular region. Further, a fixed point is set at the boundary between the breast portion 11a and the chest portion 11b.

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 breast portion 11a and the chest portion 11b. be. The following description of the breast model 11 and the affected area model 12 will be performed using a two-dimensional cross-sectional view as shown in FIG.

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 portion model 12 forming a circle is the center of gravity 12a.

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 gravity 12a and the fixed point f1 and the fixed vector V2 connecting the center of gravity 12a and the fixed point f2, these composite vectors are used. A certain movement direction vector V3 can be derived.

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 gravity 12a and is parallel to the moving direction vector V3 described with reference to FIG. Further, the symbol e indicates an action point which is an intersection of the movement direction line T and the boundary of the affected part model 12.

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 affected part model 12 are the affected parts. It is considered that it hardly moves even when it is removed. Therefore, as shown in FIG. 8B, the tension lines L5 to L11 intersecting the affected area model 12 are adopted, and the intersection points p5 to p11 between the tension lines L5 to L11 and the boundary of the affected area model 12 are obtained. Then, as shown in FIG. 9, for each tension line L5 to L11, the total length of the tension line from each operating point m5 to m11 to the action point e is calculated, and the tension from each intersection p5 to p11 to the action point e is calculated. Calculate the line part length respectively.

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.

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.

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.

The steps described so far relate to a two-dimensional cross-sectional view of the breast model 11 in which the outer surface is divided into a plurality of regions in the region along BB shown in FIG. 3 (b). By expanding this to all areas of the breast model 11, the shape of the breast after the affected part is removed can be displayed three-dimensionally.

11 to 17 are views showing steps S1 to S8 described above using a three-dimensional breast model. As shown in FIG. 11, the breast portion 11a of the present embodiment is hemispherical, and the chest portion 11b is planar. The affected area model 12 is columnar. The X-axis, Y-axis, and Z-axis, like the x-axis, y-axis, and z-axis described above, indicate the left-right direction, the up-down direction, and the front-back direction, respectively, when a person in a standing posture is used as a reference. ing.

Also in this embodiment, the breast shape change prediction system 1 identifies the breast model including the breast portion 11a and the breast portion 11b and the affected part model 12 included in the breast model as shown in FIGS. 11 to 13 (step). S1, breast model / affected area model identification step). Then, the outer surface of the breast model is divided into a plurality of regions (triangular in the present embodiment), and two or more operating points m are set at the boundary between these regions in the breast portion 11a (see FIG. 15). Two or more fixed points f are set at the boundary between 11a and the chest portion 11b (see FIGS. 13 and 14. Step S2, operating point / fixed point setting step). Further, as shown in FIG. 12, the center of gravity 12a of the affected area model 12 is calculated (step S3, the center of gravity calculation step), and two or more fixed vectors V connecting each of the center of gravity 12a and the fixed point f are calculated (see FIG. 14). ), Further, a moving direction vector which is a composite vector of these fixed vectors V is calculated (not shown. Step S4, fixed vector / moving direction vector calculation step), and as shown in FIG. 16, the moving direction passes through the center of gravity 12a. The intersection of the moving direction line T, which is a straight line parallel to the vector, and the boundary of the affected area model 12 is set as the action point e (step S5, action point setting step).

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 area model 12 to the action point e. The tension line partial length is calculated (not shown. Step S6, tension line total length / tension line partial length calculation step), and the operating point m is set along the tension line according to the ratio of the tension line partial length to the total tension line. (Step S7, movement point movement step), and the outer surface formed by the regions connecting the movement points after movement is specified as the outer surface of the breast model after the affected part is removed (FIG. 17). See step S8, outer surface identification step). In this way, the change in the shape of the breast accompanying the removal of the affected part can be shown three-dimensionally.

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.

Here, when the shape of the breast obtained by the breast shape change prediction system 1 was compared with the shape of the breast from which the affected part was actually removed, a high correlation was observed. That is, according to the present invention, it is possible to predict the shape of the breast after removing the affected portion with high accuracy.

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.

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).

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.

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.

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.

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.

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.

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 change prediction system 1 having such high performance. I found. In addition, when connecting three affected part models in a string, the case is not limited to one in which the affected part models A, B, and C having the same shape are connected in series as shown in FIG. 28 (a), but as shown in FIG. 28 (b). , The affected part models A, B, and C may have different radii, or may be connected so as to bend in the middle.

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: Breast shape change prediction system 2: Control unit 3: Storage unit 4: Input unit 5: Display unit 6: Bus 11: Breast model 11a: Breast part 11b: Chest part 12: Affected part model 12a: Center of gravity A: First affected part Model Aa: Center of gravity of the first affected area model B: Second affected area model Ba: Center of gravity of the second affected area model L1 to L11: Tension line T: Movement direction lines V1, V2, V: Fixed vector V3: Movement direction vector e: Action Points f1, f2, f: Fixed points m1 to m11, m: Operating points

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. 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.

Images