WO2014016894A1 - 形状データ生成方法及び装置 - Google Patents

形状データ生成方法及び装置 Download PDF

Info

Publication number
WO2014016894A1
WO2014016894A1 PCT/JP2012/068634 JP2012068634W WO2014016894A1 WO 2014016894 A1 WO2014016894 A1 WO 2014016894A1 JP 2012068634 W JP2012068634 W JP 2012068634W WO 2014016894 A1 WO2014016894 A1 WO 2014016894A1
Authority
WO
WIPO (PCT)
Prior art keywords
shape
input
target
vertex
deformation
Prior art date
Application number
PCT/JP2012/068634
Other languages
English (en)
French (fr)
Inventor
真智子 中川
耕平 畠中
久田　俊明
清了 杉浦
巧 鷲尾
岡田　純一
Original Assignee
富士通株式会社
国立大学法人東京大学
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 富士通株式会社, 国立大学法人東京大学 filed Critical 富士通株式会社
Priority to JP2014526631A priority Critical patent/JP5850380B2/ja
Priority to EP12881557.8A priority patent/EP2876607B1/en
Priority to PCT/JP2012/068634 priority patent/WO2014016894A1/ja
Publication of WO2014016894A1 publication Critical patent/WO2014016894A1/ja
Priority to US14/604,258 priority patent/US9607423B2/en

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T19/00—Manipulating 3D models or images for computer graphics
    • G06T19/20—Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/00—3D [Three Dimensional] image rendering
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00—Geometric image transformations in the plane of the image
    • G06T3/10—Selection of transformation methods according to the characteristics of the input images
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00—Geometric image transformations in the plane of the image
    • G06T3/18—Image warping, e.g. rearranging pixels individually
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • G06T7/12—Edge-based segmentation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • G06T7/149—Segmentation; Edge detection involving deformable models, e.g. active contour models
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30048—Heart; Cardiac
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30101—Blood vessel; Artery; Vein; Vascular
    • G06T2207/30104—Vascular flow; Blood flow; Perfusion
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2219/00—Indexing scheme for manipulating 3D models or images for computer graphics
    • G06T2219/20—Indexing scheme for editing of 3D models
    • G06T2219/2021—Shape modification

Definitions

  • This technology relates to shape data generation technology.
  • the transformation is performed by associating the points in the standard shape with the points in the target shape, but the transformation is not performed successfully unless the points to be associated are appropriately set. There is a problem.
  • an object of the present technology is to provide a novel technology for generating highly accurate three-dimensional shape data as one aspect.
  • an input shape having a simple shape having the same topology as the target shape is obtained with respect to the target shape that is the shape of the deformation target specified by the image data.
  • identifying a first vertex satisfying a predetermined condition including a condition that a normal has an intersection with a target shape from a plurality of vertices of the input shape and
  • the identified first vertex The input shape is deformed so that the first distance shorter than the distance to the target shape is moved in the direction of the normal to the first vertex, and (D) the deformed input shape is input to the processing target. It includes a process of changing to a shape and specifying the first vertex and changing the input shape a predetermined number of times.
  • High-precision 3D shape data can be generated.
  • FIG. 1 is a functional block diagram of the shape data generation apparatus according to the first embodiment.
  • FIG. 2 is a diagram illustrating an example of segment image data.
  • FIG. 3 is a diagram illustrating a sphere that is an example of an input shape.
  • FIG. 4 is a diagram showing a main processing flow in the first embodiment.
  • FIG. 5 is a diagram illustrating an initial arrangement of input shapes.
  • FIG. 6 is a diagram illustrating a processing flow of deformation processing.
  • FIG. 7 is a diagram illustrating a process flow of the landmark setting process.
  • FIG. 8 is a diagram illustrating a processing flow of boundary point search processing.
  • FIG. 9 is a diagram for explaining the relationship between the position of the vertex v and the luminance value.
  • FIG. 9 is a diagram for explaining the relationship between the position of the vertex v and the luminance value.
  • FIG. 10 is a diagram for explaining the relationship between the position of the vertex v and the luminance value.
  • FIG. 11 is a diagram illustrating an example of a case where the search point penetrates the shape before deformation.
  • FIG. 12 is a diagram illustrating an example of a case where the search point penetrates the shape before deformation.
  • FIG. 13 is a diagram for explaining the search for boundary points.
  • FIG. 14 is a diagram showing a conventional problem.
  • FIG. 15 is a diagram for explaining the deformation process in the present embodiment.
  • FIG. 16 is a diagram for explaining Deloni triangulation.
  • FIG. 17 is a diagram for explaining Deloni triangulation.
  • FIG. 18 is a diagram for explaining Deloni triangulation.
  • FIG. 19 is a diagram for explaining Deloni triangulation.
  • FIG. 11 is a diagram illustrating an example of a case where the search point penetrates the shape before deformation.
  • FIG. 13 is a diagram for explaining the search for boundary points.
  • FIG. 20 is a diagram illustrating an example of shape data when mesh subdivision is not performed.
  • FIG. 21 is a diagram illustrating an example of shape data when mesh subdivision is performed.
  • FIG. 22 is a diagram illustrating an example of a target shape having an opening.
  • FIG. 23 is a diagram illustrating an example of a disk-shaped input shape.
  • FIG. 24 is a diagram illustrating a main processing flow according to the second embodiment.
  • FIG. 25 is a diagram illustrating an example in which a source landmark is set as an input shape.
  • FIG. 26 is a diagram illustrating an example in which a target landmark is set in the target shape.
  • FIG. 27 is a diagram for explaining the first deformation process in the second embodiment.
  • FIG. 28 is a diagram illustrating a state after performing the first deformation process in the second embodiment.
  • FIG. 29 is a diagram illustrating an example of a processing result according to the second embodiment.
  • FIG. 30 is a diagram illustrating a target shape in the third embodiment.
  • FIG. 31 is a functional block diagram of the shape data generation apparatus according to the third embodiment.
  • FIG. 32 is a diagram illustrating a main processing flow according to the third embodiment.
  • FIG. 33 is a diagram for explaining a state before the second deformation process in the third embodiment.
  • FIG. 34 is a diagram for explaining a state before the second deformation process in the third embodiment.
  • FIG. 35 is a diagram for explaining a state before the second deformation process according to the third embodiment.
  • FIG. 36 is a diagram illustrating a processing result of the third embodiment.
  • FIG. 37 is a diagram illustrating a processing result of the third embodiment.
  • FIG. 38 is a functional block diagram of a computer.
  • FIG. 1 shows a functional block diagram of the shape data generation apparatus according to the first embodiment of the present technology.
  • the shape data generation apparatus 100 includes an image data storage unit 101, an input shape data storage unit 102, a deformation processing unit 103, a landmark data storage unit 104, a deformed shape data storage unit 105, an output processing unit 106, A display unit 107 and an input unit 108 are provided.
  • the deformation processing unit 103 includes an initial setting unit 1031, a landmark setting unit 1032, a deformation unit 1033, and a mesh processing unit 1034.
  • segment image data is stored.
  • the segment image data is obtained by performing a process of filling the inside of the boundary with different luminance values for each region on a CT (Computed Tomography) image of the heart of a specific patient or the like.
  • CT Computerputed Tomography
  • FIG. 2 three-dimensional data of a target shape, which is a target shape of deformation, is obtained.
  • the layers are stacked in the order of (a), (b), (c), and (d) from the bottom.
  • the input shape data storage unit 102 stores three-dimensional data of an input shape that is a shape to be deformed.
  • a standard shape such as a heart is not used as an input shape, but a simple shape having the same topology as the shape to be expressed and having the least features is used as an input shape.
  • Such an input shape is a spherical input shape if it is a three-dimensional closed surface, and a disk-like input shape if it is a three-dimensional surface having an opening.
  • the target shape is a three-dimensional closed surface
  • the input shape is a sphere.
  • a sphere whose surface is meshed as shown in FIG. 3 is used.
  • the deformation processing unit 103 performs a process of deforming the input shape into the target shape. Specifically, the initial setting unit 1031 of the deformation processing unit 103 executes processing for arranging the input shape with respect to the target shape.
  • the landmark setting unit 1032 executes processing for setting a landmark for the input shape or the input shape after deformation.
  • the deformation unit 1033 executes a deformation process by TPS Warp. For TSPSarp processing, see Fred L. Bookstein "Principal Warps: Thin-Plate Splines and the Decomposition of Deformation", IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE. VOL. 11, NO. 585, J. Since details are described, detailed description is omitted here.
  • the mesh processing unit 1034 executes mesh subdivision processing according to, for example, Delaunay triangulation technology.
  • the Delaunay triangulation is well known and will not be described in detail.
  • the landmark data storage unit 104 stores landmark data set by the landmark setting unit 1032 and used in the changing process in the deformation unit 1033.
  • the deformed shape data storage unit 105 stores shape data during the deformation process and shape data after the deformation process is completed.
  • the output processing unit 106 generates data for displaying the deformed shape data stored in the deformed shape data storage unit 105 on the display unit 107 and outputs the generated data to the display unit 107.
  • the input unit 108 receives an instruction from the user to the deformation processing unit 103 and outputs the instruction to the deformation processing unit 103.
  • the initial setting unit 1031 of the deformation processing unit 103 stores the image stored in the image data storage unit 101.
  • the target shape is identified from the data
  • the spherical input shape data is read from the input shape data storage unit 102
  • the input shape is enlarged or reduced in accordance with the target shape, so that the center position of the target shape matches the center position of the input shape. (FIG. 4: Step S1).
  • the smallest rectangular solid (bounding box) surrounding the target shape is specified, and the average value of the lengths of the x-axis, y-axis, and z-axis sides of the bounding box is calculated.
  • the sphere which is an input shape is expanded or reduced so that it may enter the cube which has the length of the said average value.
  • the scale conversion may be performed on the sphere at different enlargement or reduction ratios in three directions so as to enter the bounding box.
  • FIG. 5 shows a state in which the sphere 11 is arranged so that the centroids thereof coincide with the target shape 10 which is the right atrium, for example.
  • the deformation processing unit 103 executes the deformation process, and stores the shape data after the completion of the deformation process in the deformed shape data storage unit 105 (step S3).
  • the deformation process will be described in detail later.
  • the output processing unit 106 performs processing for displaying the three-dimensional shape data after the completion of the deformation processing on the display unit 107 (step S5).
  • m is a variable for counting the number of processed vertices.
  • the landmark setting process will be described with reference to FIGS.
  • the landmark setting unit 1032 specifies one vertex v at random from the data of the input shape or the input shape after the previous deformation (FIG. 7: step S41). Since there is no point in selecting the same vertex, for example, an unselected vertex is selected while the value of t does not change.
  • Step S43 is a process performed to arrange the vertices v as evenly as possible in the input shape or the input shape after the previous deformation. In step S43, it is determined whether or not the following expression is satisfied.
  • d (v, v i ) represents the Euclidean distance between the point v and the point v i .
  • v i is a source landmark. However, since the source landmark is not set for the first time, it is determined that the Euclidean distance exceeds the threshold value D.
  • step S43: Yes route When it is determined that the minimum of the Euclidean distances between the vertex v and each source landmark is equal to or less than the threshold D (step S43: Yes route), the process returns to the caller process. On the other hand, when it is determined that the smallest Euclidean distance between the vertex v and each source landmark is larger than the threshold D (step S43: No route), the landmark setting unit 1032 executes boundary point search processing. (Step S45). The boundary point search process will be described with reference to FIGS.
  • the landmark setting unit 1032 calculates a unit normal vector n (v) of the vertex v (FIG. 8: Step S61).
  • n (v) is a unit normal vector with respect to the surface H at the vertex v ( ⁇ H).
  • the unit normal vector is a normal vector having a length of 1.
  • H ( ⁇ V) represents the input shape or the shape surface of the input shape after the previous deformation
  • V ( ⁇ R 3 ) is a voxel space specified by the segment image data.
  • R 3 represents a real space.
  • the segment image data takes a binary value of 0 or 1, but may take a value other than 0 or 1, or take a multi-value of 2 or more. It may be.
  • a voxel indicates a unit of the smallest cube (regular lattice) in the three-dimensional representation of digital data.
  • the landmark setting unit 1032 determines whether the vertex v exists inside the target shape (step S63). In step S63, it is determined whether or not the following expression is satisfied.
  • mapping f V ⁇ R from the voxel space V to the real number space R is defined as follows.
  • the elements of the segment image data included in the voxel space V are associated with the real number space R.
  • I is the luminance value of the voxel including the point p ( ⁇ V).
  • step S63 will be described with reference to FIGS.
  • the luminance value f (v) of the voxel space corresponding to the vertex v is larger than 0, the vertex v exists inside the target shape. Therefore, by setting the coefficient k to be incremented by 1 in the process of step S75 described later, the boundary point is searched in the direction from the inside to the outside of the target shape.
  • the luminance value f (v) of the voxel space corresponding to the vertex v becomes 0, the vertex v exists outside the target shape. Therefore, the boundary point is searched in the direction from the outside to the inside of the target shape by setting the coefficient k to be decremented by 1 in the process of step S89 described later.
  • the landmark setting unit 1032 determines whether a search point exists in the voxel space specified by the tomographic image data (step S69). In step S69, it is determined whether the following expression is satisfied.
  • step S69 When it is determined that there is no search point in the voxel space specified by the tomographic image data (step S69: No route), the process returns to the caller process. This is because the search point has reached the outside of the voxel space, so that it can be determined that the normal with respect to the vertex v and the target shape do not have an intersection.
  • step S69 when it is determined that the search point exists in the voxel space specified by the tomographic image data (step S69: Yes route), the landmark setting unit 1032 determines that the search point is a shape before deformation (input shape or after previous deformation). It is determined whether the input shape has been penetrated (step S71). In step S71, it is determined whether the following expression is satisfied.
  • mapping g V ⁇ R 3 is defined as follows. This mapping g, the elements of the segment image data included in the voxel space V is associated with the real space R 3.
  • step S71 will be described with reference to FIGS. If the search point penetrates the pre-deformation shape (the input shape or the input shape after the previous deformation) before reaching the boundary point, the search for the boundary point may not be performed properly.
  • the case where the search point penetrates the pre-deformation shape before reaching the boundary point may be, for example, the case shown in FIG. 11 or the case shown in FIG. That is, there may be a case where no boundary point exists in the search direction depending on the degree of deformation of the shape before deformation with respect to the target shape. In any case, there is a possibility that the boundary point cannot be detected or the boundary point is detected at an inappropriate position.
  • step S71 the inner product of the normal vector for the vertex v and the normal vector for the search point is calculated, and the inner product is smaller than 0 (that is, the angle formed by the normal vector is larger than 90 degrees). In addition, it is determined that the search point has penetrated the shape before deformation.
  • step S ⁇ b> 71: Yes route when it is determined that the search point has penetrated the pre-deformation shape (step S ⁇ b> 71: Yes route), an appropriate boundary point cannot be detected, and the process returns to the caller process.
  • step S71: No route when it is determined that the search point does not penetrate the pre-deformation shape (step S71: No route), the landmark setting unit 1032 displays the luminance value of the voxel space corresponding to the search point and the voxel corresponding to the vertex v. The brightness value of the space is compared, and it is determined whether the brightness value has changed significantly, that is, more than the allowable value (step S73). In step S73, it is determined whether the following expression is satisfied.
  • step S73 when it is determined that the luminance value has changed significantly (step S73: Yes route), the landmark setting unit 1032 sets the search point as a boundary point (step S77).
  • search point data for example, the value of k
  • a storage device such as a main memory.
  • step S63 No route
  • the processing in this case is only different in the search direction from the processing described above, the contents of the basic processing are as described above. That is, the process of step S79 is the same as the process of step S65, the process of step S81 is the same as the process of step S67, the process of step S83 is the same as the process of step S69, and the process of step S85 is the step. It is the same as the process of S71, and the process of step S87 is the same as the process of step S73. Therefore, a detailed description of the processing in steps S79 to S87 is omitted.
  • the search point is moved by one voxel in the normal direction from the outside to the inside of the target shape.
  • step S91 the process of step S91 is the same as the process of step S77.
  • the landmark setting unit 1032 determines whether a boundary point is detected in the boundary point search process (step S47). If it is determined that no boundary point has been detected (step S47: No route), the process returns to the caller process to process the next vertex.
  • the landmark setting unit 1032 uses the internal dividing point of the line connecting the vertex v and the boundary point v + kn (v) as the target landmark.
  • the target landmark is set at the following points.
  • the landmark setting unit 1032 sets the vertex v as a source landmark (step S51).
  • vertex identifier data is stored in the landmark data storage unit 104
  • coordinate data is stored in the landmark data storage unit 104.
  • the source landmark coordinate data is read from the input shape data storage unit 102 when the input shape is processed using the source landmark identifier, and the deformed shape data is used when the input shape after the previous deformation is processed. It is read from the storage unit 105 and used.
  • the source landmark is used in step S43, but the source landmark set in the past is the processing target.
  • the source landmark and target landmark pairs set in steps S49 and S51 are used only in the next step S31.
  • the internal dividing point of the line segment connecting the vertex in the shape before deformation (input shape or input shape after previous deformation) and the boundary point in the target shape is set as the target landmark. Will be able to.
  • the deformation processing unit 103 determines whether m ⁇ N for the variable m (step S29).
  • N is an integer set in advance. If it is determined that m ⁇ N (step S29: Yes route), the process returns to step S25 to process the next vertex.
  • step 29: No route when it is determined that m ⁇ N is not satisfied for the variable m (step 29: No route), the deformation unit 1033 of the deformation processing unit 103 is set for the same t stored in the landmark data storage unit 104. Further, the deformation processing by TPS Warp is executed according to the data of the pair of the source landmark and the target landmark, and the deformed shape data is stored in the deformed shape data storage unit 105 (step S31).
  • the source landmark is arranged in the shape before deformation
  • the target landmark is arranged at the intersection of the normal line of the source landmark and the target shape, and the deformation processing by TPS Warp is performed.
  • TPS Warp the deformation processing by TPS Warp
  • the source landmark (circle with dots by dots) arranged in the shape before the deformation, the normal of the source landmark, and the target shape Target landmarks (circles with cross-hatching) are arranged at points that internally divide the line connecting the intersections (points on the alternate long and short dash line), and deformation processing by TPS Warp is performed.
  • the source landmark (circle with dots) is set again, and the line connecting the intersection of the source landmark normal and the target shape is internally divided A target landmark is arranged, and deformation processing by TPS Warp is performed. By repeatedly performing such deformation processing, the target shape is gradually approached.
  • the mesh processing unit 1034 of the deformation processing unit 103 performs mesh subdivision processing on the shape data after the current deformation stored in the deformation shape data storage unit 105, and the processing result is converted into the deformation shape data storage unit 105. (Step S33).
  • the shape of the mesh element is biased when deforming. Therefore, in this embodiment, mesh elements whose mesh element area exceeds a predetermined threshold and the surrounding mesh elements are repeatedly meshed while being deformed. A mesh re-division process is performed as a process of subdividing. By subdividing the mesh elements, it is possible to prevent such a phenomenon that the aspect ratio is deteriorated due to the uneven shape of the mesh and the smooth shape cannot be maintained. In addition, when drawing a curved surface partitioned by mesh elements, the accuracy depends on the fineness of mesh elements (also called mesh resolution), so arbitrary shapes can be expressed with high accuracy by performing mesh subdivision processing. I can do it.
  • Delaunay triangulation is well known, but here is a brief description as an example.
  • a new point a is added inside a mesh element whose area exceeds a predetermined threshold.
  • the point a is set regardless of the area, but a new point a is set inside the mesh element whose area exceeds a predetermined threshold.
  • the circumscribed circle of each mesh element is drawn, and the mesh element that includes the point a is specified.
  • three mesh elements with hatching are specified.
  • Polygons are generated from these mesh elements (FIG. 18), and the polygons are divided into triangles at the vertices and points a of the polygons (FIG. 19). If such a process is performed, a triangular mesh including triangles of approximately uniform size is generated.
  • the transformation processing unit 103 determines whether t ⁇ T for the variable t (step S35). If it is determined that t ⁇ T (step S33: Yes route), the process returns to step S23 in order to perform further deformation processing.
  • step S35 No route
  • step S33 For example, if the mesh subdivision process in step S33 is not performed, three-dimensional shape data as shown in FIG. 20 is obtained. As can be seen from FIG. 20, three-dimensional shape data with nonuniform mesh element sizes and low accuracy is obtained. On the other hand, if the mesh subdivision process in step S33 is performed, three-dimensional shape data as shown in FIG. 21 is obtained. It can be seen that highly accurate three-dimensional shape data is obtained by the fine triangular elements.
  • a three-dimensional curved surface having an opening is assumed as the target shape. Although there are some parts that are difficult to understand, a target shape as shown in FIG. 22 is assumed. A portion surrounded by a dotted line in this target shape is an opening.
  • FIG. 23 a disk-shaped input shape as shown in FIG. 23 is used.
  • FIG. 23 it seems that the meshing is not performed, but it is assumed that the meshing is performed.
  • the shape data generation apparatus 100 shown in FIG. 1 performs processing as shown in FIGS.
  • the initial setting unit 1031 of the deformation processing unit 103 stores a disk stored in the input shape data storage unit 102.
  • Input shape data is read out, and source landmarks are set at equal intervals on the outer peripheral edge of the disk (FIG. 24: step S101).
  • a spherical source landmark is set on the outer periphery of the disc.
  • the vertex identifier of the source landmark is stored in the landmark data storage unit 104.
  • the source landmark set in this step is handled as a fixed point in the following process.
  • the initial setting unit 1031 identifies the target shape from the image data stored in the image data storage unit 101, and the number of source landmarks set as the input shape at equal intervals on the outer peripheral edge of the opening of the target shape.
  • the same number of target landmarks are set (step S103).
  • the vertex coordinate data is stored in the landmark data storage unit 104.
  • spherical target landmarks are set at equal intervals and the same number on the outer peripheral edge of the opening.
  • FIG. 27 schematically shows this deformation process.
  • the left side of FIG. 27 shows a state before the deformation, and the source landmarks arranged at equal intervals on the outer peripheral edge 21 of the input shape and the targets arranged at equal intervals on the outer peripheral edge 22 of the opening of the target shape. Deforms once toward the landmark. Then, as shown on the right side of FIG.
  • the source landmark moves to the position of the target landmark, and the outer peripheral edge 21 of the other input shape is deformed so as to substantially coincide with the outer peripheral edge 22 of the opening of the target shape. Is done.
  • the state schematically shown in FIG. 28 is obtained.
  • the filled portion on the upper surface corresponds to the input shape.
  • step S107 The basic processing flow is the processing shown in FIGS.
  • the difference from the first embodiment is that the transformation process has already been performed once, the source landmark has already been set, and the source landmark is a fixed point.
  • the difference is that data is used.
  • the point that the source landmark has already been set is different in that the Euclidean distance is calculated even in the first time in step S43 as a processing flow.
  • the point that the source landmark is a fixed point is different in that the coordinate data of the fixed point is not changed in the deformation process in step S31.
  • the output processing unit 106 performs a process for displaying the deformed shape data stored in the deformed shape data storage unit 105 on the display unit 107 (step S109).
  • FIG. 29 a processing result as shown in FIG. 29 is obtained.
  • shape data formed of approximately uniform mesh elements is obtained.
  • atria are connected to blood vessels such as pulmonary veins and vena cava and have a complicated shape as compared with the ventricles. Therefore, as shown in FIG. 30, an example will be described in which processing is performed using a shape in which the vena cava 31 extends from the right atrium as an example of a target shape.
  • FIG. 31 shows a functional block diagram of the shape data generation apparatus 200 according to the present embodiment.
  • the shape data generation device 200 includes an image data storage unit 101b, an input shape data storage unit 102, a deformation processing unit 103b, a landmark data storage unit 104, a fixed point data storage unit 201, and a deformed shape data storage unit 105.
  • the image data storage unit 101b stores both first segment image data for only the right atrium and second segment image data for the right atrium and blood vessels.
  • the fixed point data storage unit 201 stores fixed point data set at the outer peripheral edge at the boundary between the post-deformation input shape and the blood vessel obtained by performing the first deformation process as described below.
  • the output processing unit 106b may output the post-deformation input shape obtained by performing the first deformation process to the display unit 107 in response to an instruction from the deformation processing unit 103b.
  • the deformation processing unit 103b includes an initial setting unit 1031, a landmark setting unit 1032, a deformation unit 1033, a mesh processing unit 1034, and a second initial setting unit 1035.
  • the second initial setting unit 1035 executes a setting process that is performed before the second deformation process.
  • the initial setting unit 1031 of the deformation processing unit 103b starts from the first segment image data stored in the image data storage unit 101b.
  • a target shape such as the right atrium not including a blood vessel or the like is specified, spherical input shape data is read from the input shape data storage unit 102, and the input shape is enlarged or reduced according to the target shape, and the center of gravity position of the target shape is obtained. It arrange
  • the deformation processing unit 103b executes the first deformation process, and stores the shape data after the completion of the deformation process in the deformed shape data storage unit 105 (step S203).
  • the first deformation process is the same as step S3 in FIG. 4, and the process described with reference to FIGS. 6 to 8 is executed.
  • the second initial setting unit 1035 stores the post-deformation input shape, which is the processing result of the first deformation process stored in the deformed shape data storage unit 105, in the second segment stored in the image data storage unit 101b. It arrange
  • the second initial setting unit 1035 sets a fixed point at the outer peripheral edge of the boundary surface between the post-deformation input shape and the second target shape, and uses the fixed point data as the fixed point data.
  • Store in the storage unit 201 step S207).
  • fixed point vertex identifier data is stored in the fixed point data storage unit 201.
  • FIG. 34 shows an enlarged view of a portion surrounded by a circle in FIG.
  • the boundary surface 40 between the post-deformation input shape and the second target shape is specified.
  • a fixed point of a circle is set on the outer peripheral edge of the boundary surface 40, for example, at an equal distance.
  • the processing result of the first deformation process and the second target shape stored in the deformed shape data storage unit 105 are displayed on the display unit 107, and the user is allowed to set a fixed point via the input unit 108. You may do it.
  • An example of the post-deformation input shape in which a fixed point is set is shown in FIG.
  • the fixed points are represented by spheres, but are highlighted in FIG.
  • a triangular mesh element is also formed on the inner surface surrounded by the sphere.
  • transformation process part 103b implements a 2nd deformation
  • step S31 since the fixed point stored in the fixed point data storage unit 201 is not moved, the position of the vertex whose identifier is registered in the fixed point data storage unit 201 is not moved in step S31.
  • the output processing unit 106b executes processing for displaying the deformed shape data stored in the deformed shape data storage unit 105 on the display unit 107 (step S211).
  • FIG. 36 a processing result as shown in FIG. 36 is obtained.
  • FIG. 36 still shows the highlighted fixed points, but the fixed points are not actually visible. That is, three-dimensional shape data as shown in FIG. 37 is obtained as post-deformation input shape data in the second deformation process.
  • FIGS. 1 and 31 The embodiment of the present technology has been described above, but the present technology is not limited to such an embodiment.
  • the functional block diagrams shown in FIGS. 1 and 31 described above are merely examples, and may not necessarily match the actual program module configuration.
  • the order of the steps may be changed, or a plurality of steps may be executed in parallel.
  • the heart particularly the right atrium, has been described as an example.
  • the shape data generation devices 100 and 200 described above are computer devices, and display control connected to a memory 2501, a CPU 2503, a hard disk drive (HDD) 2505, and a display device 2509 as shown in FIG.
  • a unit 2507, a drive device 2513 for a removable disk 2511, an input device 2515, and a communication control unit 2517 for connecting to a network are connected by a bus 2519.
  • An operating system (OS: Operating System) and an application program for performing the processing in this embodiment are stored in the HDD 2505, and are read from the HDD 2505 to the memory 2501 when executed by the CPU 2503.
  • OS Operating System
  • the CPU 2503 controls the display control unit 2507, the communication control unit 2517, and the drive device 2513 according to the processing content of the application program, and performs a predetermined operation. Further, data in the middle of processing is mainly stored in the memory 2501, but may be stored in the HDD 2505. In an embodiment of the present technology, an application program for performing the above-described processing is stored in a computer-readable removable disk 2511 and distributed, and installed from the drive device 2513 to the HDD 2505. In some cases, the HDD 2505 may be installed via a network such as the Internet and the communication control unit 2517. Such a computer apparatus realizes various functions as described above by organically cooperating hardware such as the CPU 2503 and the memory 2501 described above and programs such as the OS and application programs. .
  • an input shape having a simple shape having the same topology as the target shape is set for the target shape that is the shape of the deformation target specified by the image data.
  • the input shape is deformed so that the first distance shorter than the distance to the target shape is moved in the direction of the normal line with respect to the first vertex, and (D) the deformed input shape is the input shape to be processed. And a process of executing the specification of the first vertex and the change of the input shape a predetermined number of times.
  • shape data having the same shape as the target shape can be generated with high accuracy.
  • the mesh subdivision process may be executed for mesh elements whose area exceeds the threshold and surrounding mesh elements among a plurality of mesh elements in the deformed shape.
  • the spherical input shape is reduced or enlarged according to the target shape and arranged at the center of gravity of the target shape. Also good. In this way, the deformation process can be appropriately executed even for a three-dimensional closed surface.
  • the input shape of the flat plate may be deformed according to the opening.
  • the edge portion of the input shape deformed in accordance with the opening may be fixed. In this way, the deformation process can be appropriately executed even for a three-dimensional curved surface having an opening.
  • the present shape data generation method includes (E) specifying the first vertex and the input shape with respect to the second target shape specified by the second image data and having an additional portion added to the target shape.
  • Second input shape that is a second input shape that is a shape after a predetermined number of changes have been executed and is not moved to the outer peripheral edge of the boundary surface with the additional portion of the second target shape (F) from the plurality of vertices of the second input shape, specify a second vertex that satisfies a predetermined condition including a condition that the normal has an intersection with the second target shape, and (G) The second input shape excluding the fixed point so that the identified second vertex is moved in the direction of the normal to the second vertex by a second distance shorter than the distance to the second target shape. And (H) the second input shape after the deformation is the second to be processed. Change of the specific and the second input shape of the second vertex to change the input shape may be a second through a certain number of executions.
  • the luminance value at the destination point is compared with the luminance value at the target vertex to determine whether the luminance value has changed, and (X4) the luminance value changes If it is determined that (X5) When it is determined that the destination point is not included in the voxel space, or (X5) When it is determined that the point has penetrated the input shape or the second input shape, it is determined that the normal of the target vertex has an intersection with the input shape or the second input shape, and (X6 ) When it is determined that the luminance value has not changed, the processing after the determination of whether the destination point is included in the voxel space specified by the image data or the second image data is executed again for the destination point. You may make it do.
  • a program for causing a computer to perform the processing according to the above method can be created.
  • the program can be a computer-readable storage medium such as a flexible disk, a CD-ROM, a magneto-optical disk, a semiconductor memory, or a hard disk. It is stored in a storage device.
  • the intermediate processing result is temporarily stored in a storage device such as a main memory.

Landscapes

  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Software Systems (AREA)
  • Computer Graphics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Architecture (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Quality & Reliability (AREA)
  • Processing Or Creating Images (AREA)
  • Image Generation (AREA)

Abstract

　本形状データ生成方法は、画像データにより特定される変形目標の形状である目標形状に対して、当該目標形状と同じトポロジを有する単純な形状を有する入力形状を設定し、入力形状の複数の頂点から、法線が目標形状と交点を有するという条件を含む所定の条件を満たす第１の頂点を特定し、特定された第１の頂点を、当該第１の頂点についての法線の方向に、目標形状までの距離より短い第１の距離移動させるように、入力形状を変形し、変形後の入力形状を処理対象の入力形状に変更して第１の頂点の特定及び入力形状の変更を所定回数実行する処理を含む。

Description

形状データ生成方法及び装置
　本技術は、形状データ生成技術に関する。
　近年、ＨＰＣ（High Performance Computing）技術による計算機の演算能力の向上に伴い、人体の心臓等の臓器の動きを再現する生体シミュレーションが注目されている。この生体シミュレーションを実行するために、複雑な内部情報を有する３次元の臓器等の形状データを用いる場合がある。
　臓器等の形状データを生成するために、臓器を表す標準的な形状を、画像データが表す臓器の形状に変形する方法が知られている。
　この方法では、例えば、標準的な形状における点と目標となる形状における点との対応付けを行うことにより変形を行うが、対応付けを行う点を適切に設定しなければ変形がうまく行われないという問題がある。
　また、心臓などの臓器の場合、標準形状と個々の患者についての目標形状とに大きな差があることもあり、血管位置等の差などから、標準形状から精度良く変形できない場合があることも分かってきた。
特開２００２－３２９２１６号公報 特開２０１０－６１４３１号公報 特開２００１－３４７７４号公報 特開２００７－９８０２８号公報
Fred L. Bookstein "Principal Warps: Thin-Plate Splines and the Decomposition of Deformation", IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE. VOL. 11, NO. 6, PP. 567-585, JUNE 1989 "Laplacian surface editing", SGP '04 Proceedings of the 2004 Eurographics/ACM SIGGRAPH symposium on Geometry processing, O. Sorkine Tel Aviv University D. Cohen-Or Tel Aviv University Y. Lipman Tel Aviv University M. Alexa Darmstadt University of Technology C. Roessl Max-Planck Institut fuer Informatik, Saarbruecken H.-P. Seidel Max-Planck Institut fuer Informatik, Saarbruecken
　従って、本技術の目的は、一側面として、高精度の三次元形状データを生成するための新規な技術を提供することである。
　本技術の一態様に係る形状データ生成方法は、（Ａ）画像データにより特定される変形目標の形状である目標形状に対して、当該目標形状と同じトポロジを有する単純な形状を有する入力形状を設定し、（Ｂ）入力形状の複数の頂点から、法線が目標形状と交点を有するという条件を含む所定の条件を満たす第１の頂点を特定し、（Ｃ）特定された第１の頂点を、当該第１の頂点についての法線の方向に、目標形状までの距離より短い第１の距離移動させるように、入力形状を変形し、（Ｄ）変形後の入力形状を処理対象の入力形状に変更して第１の頂点の特定及び入力形状の変更を所定回数実行する処理を含む。
　高精度の３次元形状データを生成することができるようになる。
図１は、第１の実施の形態に係る形状データ生成装置の機能ブロック図である。 図２は、セグメント画像データの一例を示す図である。 図３は、入力形状の一例である球を示す図である。 図４は、第１の実施の形態におけるメインの処理フローを示す図である。 図５は、入力形状の初期配置を示す図である。 図６は、変形処理の処理フローを示す図である。 図７は、ランドマーク設定処理の処理フローを示す図である。 図８は、境界点探索処理の処理フローを示す図である。 図９は、頂点ｖの位置と輝度値との関係について説明するための図である。 図１０は、頂点ｖの位置と輝度値との関係について説明するための図である。 図１１は、探索点が変形前形状を貫通するケースの一例を示す図である。 図１２は、探索点が変形前形状を貫通するケースの一例を示す図である。 図１３は、境界点の探索を説明するための図である。 図１４は、従来の問題点を示す図である。 図１５は、本実施の形態における変形処理を説明するための図である。 図１６は、デローニ三角形分割を説明するための図である。 図１７は、デローニ三角形分割を説明するための図である。 図１８は、デローニ三角形分割を説明するための図である。 図１９は、デローニ三角形分割を説明するための図である。 図２０は、メッシュ再分割を行わない場合の形状データの一例を示す図である。 図２１は、メッシュ再分割を行った場合の形状データの一例を示す図である。 図２２は、開口部を有する目標形状の一例を示す図である。 図２３は、円板状の入力形状の一例を示す図である。 図２４は、第２の実施の形態に係るメインの処理フローを示す図である。 図２５は、入力形状にソースランドマークを設定した例を示す図である。 図２６は、目標形状にターゲットランドマークを設定した例を示す図である。 図２７は、第２の実施の形態における最初の変形処理を説明するための図である。 図２８は、第２の実施の形態における最初の変形処理を実施した後の状態を示す図である。 図２９は、第２の実施の形態の処理結果の一例を示す図である。 図３０は、第３の実施の形態における目標形状を示す図である。 図３１は、第３の実施の形態に係る形状データ生成装置の機能ブロック図である。 図３２は、第３の実施の形態に係るメインの処理フローを示す図である。 図３３は、第３の実施の形態における第２の変形処理の前の状態を説明するための図である。 図３４は、第３の実施の形態における第２の変形処理の前の状態を説明するための図である。 図３５は、第３の実施の形態における第２の変形処理の前の状態を説明するための図である。 図３６は、第３の実施の形態の処理結果を示す図である。 図３７は、第３の実施の形態の処理結果を示す図である。 図３８は、コンピュータの機能ブロック図である。
［実施の形態１］
　図１に、本技術の第１の実施の形態に係る形状データ生成装置の機能ブロック図を示す。形状データ生成装置１００は、画像データ格納部１０１と、入力形状データ格納部１０２と、変形処理部１０３と、ランドマークデータ格納部１０４と、変形形状データ格納部１０５と、出力処理部１０６と、表示部１０７と、入力部１０８とを有する。
　変形処理部１０３は、初期設定部１０３１と、ランドマーク設定部１０３２と、変形部１０３３と、メッシュ処理部１０３４とを有する。
　画像データ格納部１０１には、セグメント画像データが格納されている。セグメント画像データは、特定の患者の心臓等のＣＴ（Computed Tomography）画像等に対して、部位毎に境界の内部を異なる輝度値で塗りつぶす処理を施すことにより得られる。例えば図２に模式的に示すようなセグメント画像データを積層することにより、変形の目標となる形状である目標形状の三次元のデータが得られる。図２の例では、下から（ａ）（ｂ）（ｃ）（ｄ）の順番で積層される。
　入力形状データ格納部１０２には、変形される形状である入力形状の三次元データが格納される。なお、本実施の形態では、心臓などの標準形状を入力形状には使わずに、表現したい形状と同様なトポロジを有する最も特徴の無い単純な形状を入力形状として用いる。このような入力形状は、三次元閉曲面であれば球状の入力形状であり、開口部を有する三次元曲面であれば、円板状の入力形状である。本実施の形態では、目標形状が三次元閉曲面である例を説明するので、入力形状は球となる。例えば、図３に示すような表面がメッシングされた球が用いられる。
　変形処理部１０３は、入力形状を目標形状に変形する処理を行う。具体的には、変形処理部１０３の初期設定部１０３１は、目標形状に対して入力形状を配置する処理を実行する。ランドマーク設定部１０３２は、入力形状又は変形後の入力形状に対してランドマークを設定する処理を実行する。変形部１０３３は、ＴＰＳＷａｒｐによる変形処理を実行する。ＴＰＳＷａｒｐ処理については、Fred L. Bookstein "Principal Warps: Thin-Plate Splines and the Decomposition of Deformation", IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE. VOL. 11, NO. 6, PP. 567-585, JUNE 1989に詳細が説明されているので、ここでは詳細な説明を省略する。
　メッシュ処理部１０３４は、例えばDelaunay（デローニ）三角形分割技術に従って、メッシュ再分割処理を実行する。Delaunay三角形分割については、よく知られているので詳細な説明は省略する。
　ランドマークデータ格納部１０４は、ランドマーク設定部１０３２により設定され且つ変形部１０３３での変更処理において用いられるランドマークのデータを格納する。
　変形形状データ格納部１０５は、変形処理途中の形状データ及び変形処理完了後の形状データを格納する。出力処理部１０６は、変形形状データ格納部１０５に格納されている変形処理完了後の形状データを、表示部１０７に表示するためのデータを生成し、表示部１０７に出力する。
　入力部１０８は、ユーザから変形処理部１０３に対する指示などを受け付け、変形処理部１０３に出力する。
　次に、図４乃至図２１を用いて、形状データ生成装置１００の処理内容を説明する。
　例えば、ユーザから目標形状が右心房などの三次元閉曲面である旨の指示を入力部１０８から受け付けると、変形処理部１０３の初期設定部１０３１は、画像データ格納部１０１に格納されている画像データから目標形状を特定し、入力形状データ格納部１０２から球状の入力形状のデータを読み出して目標形状に合わせて入力形状を拡大又は縮小し、目標形状の重心位置と入力形状の重心位置が一致するように配置する（図４：ステップＳ１）。
　より詳しくは、目標形状を囲う最小の直方体（バウンディングボックス）を特定し、当該バウンディングボックスのｘ軸、ｙ軸及びｚ軸の辺の長さの平均値を算出する。そして、当該平均値の長さを有する立方体に入るように、入力形状である球を拡大又は縮小する。但し、バウンディングボックスに入るように、３方向で異なる拡大又は縮小率で球に対してスケール変換を実行するようにしても良い。
　なお、等スケールで拡大又は縮小する場合には、例えば図５に示すように初期配置がなされる。図５は、例えば右心房である目標形状１０に対して球１１が、それらの重心が一致するように配置されている状態を示している。
　次に、変形処理部１０３は、変形処理を実行し、変形処理完了後の形状データを変形形状データ格納部１０５に格納する（ステップＳ３）。変形処理については、後に詳しく述べる。
　その後、出力処理部１０６は、変形処理完了後の３次元形状データを表示部１０７に表示するための処理を行う（ステップＳ５）。
　次に、図６乃至図２１を用いて変形処理について説明する。
　まず、変形処理部１０３は、変形回数をカウントするための変数ｔの初期値をｔ＝０と設定する（図６：ステップＳ２１）。次に、変形処理部１０３は、変数ｔをｔ＝ｔ＋１と増分することにより変形回数をカウントするとともに、また変数ｍの初期値をｍ＝０と設定する（ステップＳ２３）。ｍは、処理した頂点数をカウントするための変数である。
　そして、変形処理部１０３は、変数ｍをｍ＝ｍ＋１と増分し（ステップＳ２５）、変形処理部１０３のランドマーク設定部１０３２は、ランドマーク設定処理を実行する（ステップＳ２７）。ランドマーク設定処理については、図７乃至図１３を用いて説明する。
　まず、ランドマーク設定部１０３２は、入力形状又は前回変形後の入力形状のデータから、頂点ｖをランダムに１つ特定する（図７：ステップＳ４１）。同じ頂点を選択する意味はないので、例えばｔの値が変化しない間で未選択の頂点を選択する。
　そして、ランドマーク設定部１０３２は、ランドマークデータ格納部１０４に格納されている各ソースランドマークと頂点ｖとの間のユークリッド距離をそれぞれ算出する。そして、ランドマーク設定部１０３２は、頂点ｖと各ソースランドマークとのユークリッド距離のうち最小のものが、閾値Ｄ以下であるか判断する（ステップＳ４３）。ステップＳ４３は、頂点ｖを入力形状又は前回変形後の入力形状において可能な限り均等に配置するために行われる処理である。ステップＳ４３においては、以下の式を満たすか否かを判断する。
Figure JPOXMLDOC01-appb-M000001
　ここで、ｄ（ｖ、ｖi）は、点ｖと点ｖiとの間のユークリッド距離を表す。ｖiは、ソースランドマークである。但し、初回は、ソースランドマークが設定されていないので、ユークリッド距離は閾値Ｄを超えると判断するものとする。
　頂点ｖと各ソースランドマークとのユークリッド距離のうち最小のものが閾値Ｄ以下であると判断された場合（ステップＳ４３：Ｙｅｓルート）、呼出元の処理に戻る。一方、頂点ｖと各ソースランドマークとのユークリッド距離のうち最小のものが閾値Ｄより大きいと判断された場合（ステップＳ４３：Ｎｏルート）、ランドマーク設定部１０３２は、境界点探索処理を実行する（ステップＳ４５）。境界点探索処理については、図８乃至図１３を用いて説明する。
　まず、ランドマーク設定部１０３２は、頂点ｖの単位法線ベクトルｎ（ｖ）を算出する（図８：ステップＳ６１）。ここで、ｎ（ｖ）は、頂点ｖ（∈Ｈ）での面Ｈに対する単位法線ベクトルである。単位法線ベクトルとは、長さが１の法線ベクトルをいう。なお、Ｈ（⊂Ｖ）は、入力形状又は前回変形後の入力形状の形状面を表し、Ｖ（⊂Ｒ3）はセグメント画像データにより特定されるボクセル空間である。また、Ｒ3は、実数空間を表す。ここでは簡単化のため、セグメント画像データは、０又は１の２値を取るものとするが、０又は１以外の値を取るものであっても良いし、２値以上の多値を取るものであっても良い。また、ボクセルは、デジタルデータの立体表現において、その最小の立方体（正規格子）の単位を示す。
　また、ランドマーク設定部１０３２は、頂点ｖが目標形状の内側に存在するか判断する（ステップＳ６３）。ステップＳ６３においては、以下の式を満たすか否かを判断する。
Figure JPOXMLDOC01-appb-M000002
　ここで、ボクセル空間Ｖから実数空間Ｒへの写像ｆ：Ｖ→Ｒは以下のように定義される。この写像ｆにより、ボクセル空間Ｖに含まれるセグメント画像データの要素が実数空間Ｒに対応付けられる。
Figure JPOXMLDOC01-appb-M000003
　ここで、Ｉは点ｐ（∈Ｖ）を含むボクセルの輝度値である。
　図９及び図１０を用いて、ステップＳ６３の処理について説明する。図９に示すように、頂点ｖに対応するボクセル空間の輝度値ｆ（ｖ）が０より大きい場合には、頂点ｖが目標形状の内側に存在する。そこで、後に説明するステップＳ７５の処理において係数ｋを１ずつインクリメントするように設定することで、目標形状の内から外へ向かう方向に境界点を探索する。一方、図１０に示すように、頂点ｖに対応するボクセル空間の輝度値ｆ（ｖ）が０になる場合には、頂点ｖが目標形状の外側に存在する。そこで、後に述べるステップＳ８９の処理において係数ｋを１ずつデクリメントするように設定することで、目標形状の外から内へ向かう方向に境界点を探索する。
　そして、頂点ｖが目標形状の内側に存在すると判断された場合（ステップＳ６３：Ｙｅｓルート）、ランドマーク設定部１０３２は、係数ｋをｋ＝０と設定する（ステップＳ６５）。また、ランドマーク設定部１０３２は、境界点であるかどうかを判定する点（以下、探索点と呼ぶ）を以下のように設定する（ステップＳ６７）。
Figure JPOXMLDOC01-appb-M000004
　そして、ランドマーク設定部１０３２は、断層画像データにより特定されるボクセル空間内に探索点が存在するか判断する（ステップＳ６９）。ステップＳ６９においては、以下の式を満たすか判断する。
Figure JPOXMLDOC01-appb-M000005
　断層画像データにより特定されるボクセル空間内に探索点が存在しないと判断された場合（ステップＳ６９：Ｎｏルート）、呼出元の処理に戻る。探索点がボクセル空間外に達したため、頂点ｖについての法線と目標形状とが交点を有さないと判断できるためである。
　一方、断層画像データにより特定されるボクセル空間内に探索点が存在すると判断された場合（ステップＳ６９：Ｙｅｓルート）、ランドマーク設定部１０３２は、探索点が変形前形状（入力形状又は前回変形後の入力形状）を貫通したか判断する（ステップＳ７１）。ステップＳ７１においては、以下の式を満たすか判断する。
Figure JPOXMLDOC01-appb-M000006
　ここで、写像ｇ：Ｖ→Ｒ3は以下のように定義される。この写像ｇにより、ボクセル空間Ｖに含まれるセグメント画像データの要素が実数空間Ｒ3に対応付けられる。
Figure JPOXMLDOC01-appb-M000007
　ここで、写像ｇの制限ｇ|Ｈはｎ(ｖ)になることに注意する。
　図１１及び図１２を用いて、ステップＳ７１の処理について説明する。探索点が境界点に達する前に変形前形状（入力形状又は前回変形後の入力形状）を貫通するような場合には、境界点の探索を適切に行えていない可能性がある。探索点が境界点に達する前に変形前形状を貫通するような場合とは、例えば、図１１に示したような場合と、図１２に示したような場合とが考えられる。すなわち、目標形状に対する変形前形状の変形度合いにより探索方向に境界点が存在しない場合が考えられる。いずれのケースにおいても、境界点を検出できないか、又は適切でない位置に境界点を検出してしまう可能性がある。そこで、ステップＳ７１においては、頂点ｖについての法線ベクトルと探索点についての法線ベクトルとの内積を算出し、内積が０より小さい（すなわち、法線ベクトルのなす角が９０度より大きい）場合に、探索点が変形前形状を貫通したと判断している。
　図８の説明に戻り、探索点が変形前形状を貫通したと判断された場合（ステップＳ７１：Ｙｅｓルート）、適切な境界点を検出できないので、呼出元の処理に戻る。一方、探索点が変形前形状を貫通していないと判断された場合（ステップＳ７１：Ｎｏルート）、ランドマーク設定部１０３２は、探索点に対応するボクセル空間の輝度値と頂点ｖに対応するボクセル空間の輝度値とを比較し、輝度値が有意に、すなわち許容値以上に変化したかを判断する（ステップＳ７３）。ステップＳ７３においては、以下の式を満たすか判断する。
Figure JPOXMLDOC01-appb-M000008
　そして、輝度値が有意には変化していないと判断された場合（ステップＳ７３：Ｎｏルート）、ランドマーク設定部１０３２は、係数ｋをｋ＝ｋ＋１とインクリメントし（ステップＳ７５）、ステップＳ６７の処理に戻る。
　このようにすると、図１３に示すように、頂点ｖから法線方向へ１ボクセル分ずつ探索点を移動させつつ、境界点が存在するかを判定することができる。
　一方、輝度値が有意に変化したと判断された場合（ステップＳ７３：Ｙｅｓルート）、ランドマーク設定部１０３２は、探索点を境界点に設定する（ステップＳ７７）。ステップＳ７７においては、探索点のデータ（例えばｋの値）をメインメモリ等の記憶装置に格納しておく。そして呼出元の処理に戻る。
　これに対して、ステップＳ６３において、頂点ｖが目標形状の外側に存在すると判断された場合（ステップＳ６３：Ｎｏルート）に行われる処理について説明する。この場合の処理は、上で述べた処理とは探索方向が異なるだけであるから、基本的な処理の内容は上で述べたとおりである。すなわち、ステップＳ７９の処理はステップＳ６５の処理と同様であり、ステップＳ８１の処理はステップＳ６７の処理と同様であり、ステップＳ８３の処理はステップＳ６９の処理と同様であり、ステップＳ８５の処理はステップＳ７１の処理と同様であり、ステップＳ８７の処理はステップＳ７３の処理と同様である。従って、ステップＳ７９乃至Ｓ８７の処理の詳細な説明は省略する。
　そして、ステップＳ８９において、ランドマーク設定部１０３２は、係数ｋをｋ＝ｋ－１とデクリメントし（ステップＳ８９）、ステップＳ８１の処理に戻る。これにより、探索点が、目標形状の外側から内側へ向かう法線方向に１ボクセル分移動される。また、ステップＳ９１の処理は、ステップＳ７７の処理と同様である。
　以上のような処理を実施することにより、頂点ｖについての法線と目標形状との交点（すなわち境界点）を検出することができるようになる。
　図７の説明に戻り、ランドマーク設定部１０３２は、境界点探索処理において境界点が検出されたか判断する（ステップＳ４７）。境界点が検出されていないと判断された場合（ステップＳ４７：Ｎｏルート）、次の頂点について処理するため、呼出元の処理に戻る。
　一方、境界点が検出されたと判断された場合（ステップＳ４７：Ｙｅｓルート）、ランドマーク設定部１０３２は、頂点ｖと境界点ｖ＋ｋｎ（ｖ）とを結ぶ線分の内分点をターゲットランドマークに設定する（ステップＳ４９）。具体的には、ターゲットランドマークは以下のような点に設定される。
Figure JPOXMLDOC01-appb-M000009
　そして、ランドマーク設定部１０３２は、頂点ｖをソースランドマークに設定する（ステップＳ５１）。ソースランドマークについては頂点識別子のデータをランドマークデータ格納部１０４に格納し、対応するターゲットランドマークについては座標データをランドマークデータ格納部１０４に格納する。ソースランドマークの座標データについては、ソースランドマークの識別子を用いて、入力形状を処理する場合には入力形状データ格納部１０２から読み出し、前回変形後の入力形状を処理する場合には変形形状データ格納部１０５から読み出して用いる。ソースランドマークについては、ステップＳ４３で用いられるが、過去に設定されたソースランドマークが処理対象となる。ステップＳ４９及びＳ５１で設定されたソースランドマーク及びターゲットランドマークの対については、次のステップＳ３１でのみ用いられる。
　以上のような処理を実行することにより、変形前の形状（入力形状又は前回変形後の入力形状）における頂点と目標形状における境界点とを結ぶ線分の内分点を、ターゲットランドマークに設定することができるようになる。
　図６の説明に戻り、変形処理部１０３は、変数ｍについてｍ＜Ｎであるか判断する（ステップＳ２９）。ここで、Ｎは、予め設定された整数である。ｍ＜Ｎであると判断された場合（ステップＳ２９：Ｙｅｓルート）、次の頂点について処理するため、ステップＳ２５の処理に戻る。
　一方、変数ｍについてｍ＜Ｎではないと判断された場合（ステップ２９：Ｎｏルート）、変形処理部１０３の変形部１０３３は、ランドマークデータ格納部１０４に格納されている、同じｔについて設定された、ソースランドマーク及びターゲットランドマークの対のデータに従い、ＴＰＳＷａｒｐによる変形処理を実行し、変形後の形状データを変形形状データ格納部１０５に格納する（ステップＳ３１）。
　図１４に示すように、ソースランドマークを変形前の形状に配置し、ソースランドマークの法線と目標形状との交点にターゲットランドマークを配置し、ＴＰＳＷａｒｐによる変形処理を行うことが考えられる。しかしながら、図１４に示すように、法線が交差するような状況が発生すると、変形後の形状には、目標形状とは異なる不自然な形状が生じてしまうことがある。
　そこで、本実施の形態における変形処理においては、図１５に示すように、変形前の形状に配置されたソースランドマーク（ドットによるハッチング付き丸）と、ソースランドマークの法線と目標形状との交点とを結ぶ線分を内分する点（一点鎖線上の点）に、ターゲットランドマーク（クロスハッチング付き丸）を配置し、ＴＰＳＷａｒｐによる変形処理を行う。この変形処理で得られた変形後形状においてソースランドマーク（ドットによるハッチング付き丸）を再度設定して、ソースランドマークの法線と目標形状との交点とを結ぶ線分を内分する点にターゲットランドマークを配置して、ＴＰＳＷａｒｐによる変形処理を行う。このような変形処理を繰り返し行うことにより、徐々に目標形状に近付けていく。これにより、変形後の形状に不自然な形状が生じにくくなり、また法線の方向が本来目標とすべき部分に向きやすくなる。なお、内分比は徐々に大きくなるが、ソースランドマークとターゲットランドマークとの距離も短くなるので、適切な変形が行われる。
　そして、変形処理部１０３のメッシュ処理部１０３４は、変形形状データ格納部１０５に格納された、今回の変形後の形状データについて、メッシュ再分割処理を実行し、処理結果を変形形状データ格納部１０５に格納する（ステップＳ３３）。
　目標形状が複雑であるほど変形の際メッシュ要素の形状に偏りが生じるため、本実施の形態では、変形を繰り返しながらメッシュ要素の面積が所定の閾値を超えるメッシュ要素及び周辺のメッシュ要素についてはメッシュの細分化を行う処理としてメッシュ再分割処理を行う。メッシュ要素の細分化を行うことで、メッシュの形状の偏りのためアスペクト比が悪化して滑らかな形状が保てなくなる等の現象を防ぐことが出来る。また、メッシュ要素で区画される曲面を描画する場合には、その精度はメッシュ要素の細かさ（メッシュ解像度とも呼ぶ）に依存するため、メッシュ再分割処理を行うことで、任意形状を精度良く表現することが出来る。
　特に、特徴の無い単純な形状を用いているため変形に偏りが生じ、結果としてメッシュの不均等が起こりやすい。よって、例えば変形前のメッシュ要素のサイズよりも大きくなったメッシュ要素については、当該メッシュ要素等を再分割することによって不都合を解消する。
　メッシュ再分割処理におけるメッシュ分割手法には、八分木メッシュ分割法やDelaunay三角形分割等様々な手法がある。Delaunay三角分割はよく知られているが、一例として、ここで簡単に説明しておく。まず、面積が所定の閾値を超えるメッシュ要素の内部に新たな点ａを追加する。図１６の例では、面積とは関係なく点ａが設定されているが、面積が所定の閾値を超えるメッシュ要素内部に新たな点ａを設定する。次に、各メッシュ要素の外接円を描き、点ａを内包するメッシュ要素を特定する。図１７の例では、ハッチングが付された３つのメッシュ要素が特定される。これらのメッシュ要素で多角形を生成し（図１８）、この多角形の頂点と点ａとで、多角形を三角形に分割する（図１９）。このような処理を行えば、おおよそ均一サイズの三角形が含まれる三角形メッシュが生成されるようになる。
　そして、変形処理部１０３は、変数ｔについてｔ＜Ｔであるか判断する（ステップＳ３５）。ｔ＜Ｔであると判断された場合（ステップＳ３３：Ｙｅｓルート）、さらに変形処理を行うため、ステップＳ２３の処理に戻る。なお、Ｔは総変形回数であり、予め管理者等により設定される（例えばＴ＝５００）。
　一方、変数ｔについてｔ＜Ｔでないと判断された場合（ステップＳ３５：Ｎｏルート）、Ｔ回の変形が終了したため、呼出元の処理に戻る。
　以上のような処理を実行することにより、精度の高い３次元形状データを得ることができるようになる。
　例えば、ステップＳ３３のメッシュ再分割処理を実施しないと、図２０に示すような三次元形状データが得られる。図２０から分かるように、メッシュ要素のサイズは不均一で精度の低い三次元形状データが得られている。一方、ステップＳ３３のメッシュ再分割処理を実施すれば、図２１に示すような三次元形状データが得られる。細かい三角形要素により精度の高い三次元形状データが得られていることが分かる。
［実施の形態２］
　本実施の形態では、目標形状として開口部を有する三次元曲面を想定する。わかりにくい部分もあるが、図２２に示すような目標形状を想定する。この目標形状において点線で囲った部分が開口部である。
　この場合、本実施の形態では、例えば図２３に示すような円板状の入力形状を用いる。図２３ではメッシングされていないように見えるが、メッシングされているものとする。
　このような場合に、図１に示した形状データ生成装置１００は、図２４乃至図２９で示すような処理を実施する。
　まず、ユーザから目標形状が開口を有する三次元曲面である旨の指示を入力部１０８から受け付けると、変形処理部１０３の初期設定部１０３１は、入力形状データ格納部１０２に格納されている円板状の入力形状のデータを読み出し、その円板の外周エッジに等間隔でソースランドマークを設定する（図２４：ステップＳ１０１）。例えば、図２５に模式的に示すように、円板の外周上に球形のソースランドマークが設定される。なお、ソースランドマークについては、その頂点識別子が、ランドマークデータ格納部１０４に格納されている。なお、本ステップで設定されたソースランドマークは、固定点として以下の処理で扱われる。
　さらに、初期設定部１０３１は、画像データ格納部１０１に格納されている画像データから目標形状を特定し、当該目標形状の開口部の外周エッジに等間隔に、入力形状に設定したソースランドマーク数と同数のターゲットランドマークを設定する（ステップＳ１０３）。ターゲットランドマークについては、その頂点座標のデータが、ランドマークデータ格納部１０４に格納される。図２６に模式的に示すように、開口部の外周エッジに、球形のターゲットランドマークが等間隔及び同数設定される。
　そうすると、初期設定部１０３１は、変形部１０３３に対して、ソースランドマークをターゲットランドマークの位置に移動するように、ＴＰＳＷａｒｐによる変形処理を実行させ、変形後の形状データを変形形状データ格納部１０５に格納する（ステップＳ１０５）。図２７に、本変形処理を模式的に示す。図２７の左側は、変形前の状態を示しており、入力形状の外周エッジ２１に等間隔に配置されたソースランドマークを、目標形状の開口部の外周エッジ２２に等間隔に配置されたターゲットランドマークに向けて１回で変形する。そうすると、図２７の右側に示すように、ソースランドマークがターゲットランドマークの位置に移動して、その他の入力形状の外周エッジ２１はおおよそ目標形状の開口部の外周エッジ２２に一致するように変形される。図２５及び図２６の例では、図２８に模式的に示すような状態が得られる。図２８の例では、上面において塗りつぶされた部分が入力形状に相当する。
　そして、変形処理部１０３は、変形処理を実行する（ステップＳ１０７）。基本的な処理の流れは、図６乃至図８の処理を行うことになる。但し、第１の実施の形態と異なる点は、既に１回変形処理が行われている点、既にソースランドマークが設定されている点、当該ソースランドマークは固定点である点が異なる。既に１回変形処理が行われている点については、処理の流れとしては、ｔ＝１の場合においてもランドマーク設定処理において変形形状データ格納部１０５に格納されている前回変形後の入力形状のデータを用いる点が異なる。
　既にソースランドマークが設定されている点は、処理の流れとしてはステップＳ４３において初回でもユークリッド距離が算出される点で異なる。ソースランドマークが固定点である点については、ステップＳ３１における変形処理において固定点の座標データは変更されないという点で異なる。
　それ以外の処理については第１の実施の形態と同様である。
　そして、出力処理部１０６は、変形形状データ格納部１０５に格納されている変形後の形状データを表示部１０７に表示するための処理を実行する（ステップＳ１０９）。
　例えば上で述べた例では、図２９に示すような処理結果が得られる。図２９の例では、おおよそ均等なメッシュ要素で形成された形状データが得られていることが分かる。
［実施の形態３］
　心房は、心室に比して、肺静脈や大静脈などの血管に繋がっており複雑な形状を有している。そこで、図３０に示すように、右心房から大静脈３１が伸びている形状を目標形状の一例として処理を行う例を説明する。
　本実施の形態に係る形状データ生成装置２００の機能ブロック図を図３１に示す。なお、図１に示した形状データ生成装置１００とほぼ同様の機能を有する場合には同じ参照符号を付している。形状データ生成装置２００は、画像データ格納部１０１ｂと、入力形状データ格納部１０２と、変形処理部１０３ｂと、ランドマークデータ格納部１０４と、固定点データ格納部２０１と、変形形状データ格納部１０５と、出力処理部１０６ｂと、表示部１０７と、入力部１０８とを有する。
　画像データ格納部１０１ｂには、右心房のみについての第１のセグメント画像データと、右心房及び血管についての第２のセグメント画像データとが両方格納されている。
　固定点データ格納部２０１は、以下で述べるように１回目の変形処理を実施することで得られる変形後入力形状と血管との境界部分における外周エッジに設定される固定点のデータを格納する。
　出力処理部１０６ｂは、変形処理部１０３ｂからの指示に応じて、１回目の変形処理を実施することで得られる変形後入力形状を表示部１０７に出力する場合もある。
　さらに、変形処理部１０３ｂは、初期設定部１０３１と、ランドマーク設定部１０３２と、変形部１０３３と、メッシュ処理部１０３４と、第２初期設定部１０３５とを有する。第２初期設定部１０３５は、２回目の変形処理の前に行われる設定処理を実行する。
　次に、図３２乃至図３７を用いて、本実施の形態に係る処理を説明する。
　例えば、ユーザから本実施の形態に係る処理の開始指示を入力部１０８から受け付けると、変形処理部１０３ｂの初期設定部１０３１は、画像データ格納部１０１ｂに格納されている第１のセグメント画像データから血管等を含まない右心房などの目標形状を特定し、入力形状データ格納部１０２から球状の入力形状のデータを読み出して目標形状に合わせて入力形状を拡大又は縮小し、目標形状の重心位置と入力形状の重心位置が一致するように配置する（図３２：ステップＳ２０１）。本ステップは、図４のステップＳ１と同様の処理である。
　次に、変形処理部１０３ｂは、第１変形処理を実行し、変形処理完了後の形状データを変形形状データ格納部１０５に格納する（ステップＳ２０３）。第１変形処理については、図４のステップＳ３と同様であり、図６乃至図８を用いて説明された処理を実行する。
　次に、第２初期設定部１０３５は、変形形状データ格納部１０５に格納されている第１変形処理の処理結果である変形後入力形状を、画像データ格納部１０１ｂに格納されている第２セグメント画像データから特定される右心房及び血管を含む第２の目標形状に重畳させるように配置する（ステップＳ２０５）。重畳した状態はわかりにくいが、図３３に示すような状態となる。図３３の例では丸印で囲われた部分が、血管の部分であり、この部分は第１の変形処理による変形後入力形状との差分に相当する。
　そして、第２初期設定部１０３５は、変形後入力形状において、当該変形後入力形状と第２の目標形状との境界面の外周エッジに固定点を設定し、当該固定点のデータを固定点データ格納部２０１に格納する（ステップＳ２０７）。例えば、固定点の頂点識別子のデータを、固定点データ格納部２０１に格納する。
　図３３の丸印で囲まれた部分を拡大すると図３４に示すようになる。図３４の例では、変形後入力形状と第２の目標形状との境界面４０が特定されるようになる。そうすると、この境界面４０の外周エッジに、例えば等距離で丸印の固定点を設定する。なお、変形形状データ格納部１０５に格納されている第１の変形処理の処理結果及び第２の目標形状とを表示部１０７に表示して、ユーザに入力部１０８を介して固定点を設定させるようにしても良い。固定点が設定された変形後入力形状の一例を図３５に示す。固定点は球で表されているが、図３５では強調表示されている。なお、球で囲まれた内部の面にも三角形メッシュ要素が形成されている。
　そして、変形処理部１０３ｂは、第２変形処理を実施し、処理結果を変形形状データ格納部１０５に格納する（ステップＳ２０９）。基本的には、第１の実施の形態に係る変形処理と同様である。
　但し、固定点データ格納部２０１に格納されている固定点については移動させないので、ステップＳ３１において、固定点データ格納部２０１に識別子が登録されている頂点についてはその位置を移動させないようにする。
　そして、出力処理部１０６ｂは、変形形状データ格納部１０５に格納されている変形後の形状データを表示部１０７に表示するための処理を実行する（ステップＳ２１１）。
　例えば上で述べた例では、図３６に示すような処理結果が得られる。図３６は強調表示した固定点を示したままであるが、固定点は実際には見えない。すなわち、図３７に示すような三次元形状データが第２変形処理の変形後入力形状のデータとして得られる。
　このように複雑な形状であっても精度良く形状データを生成することができるようになる。
　以上本技術の実施の形態を説明したが、本技術はこのような実施の形態に限定されるものではない。例えば、上で説明した図１及び図３１に示された機能ブロック図は一例であって、必ずしも実際のプログラムモジュール構成とは一致しない場合もある。
　また、処理フローについても、処理結果が変わらない限り、ステップの順番を入れ替えたり、複数のステップを並列実行するようにしても良い。
　さらに、上では心臓特に右心房を一例に説明したが、他の心房や心室、その他の臓器でも同様に処理できる。
　なお、上で述べた形状データ生成装置１００及び２００は、コンピュータ装置であって、図３８に示すように、メモリ２５０１とＣＰＵ２５０３とハードディスク・ドライブ（ＨＤＤ）２５０５と表示装置２５０９に接続される表示制御部２５０７とリムーバブル・ディスク２５１１用のドライブ装置２５１３と入力装置２５１５とネットワークに接続するための通信制御部２５１７とがバス２５１９で接続されている。オペレーティング・システム（ＯＳ：Operating System）及び本実施例における処理を実施するためのアプリケーション・プログラムは、ＨＤＤ２５０５に格納されており、ＣＰＵ２５０３により実行される際にはＨＤＤ２５０５からメモリ２５０１に読み出される。ＣＰＵ２５０３は、アプリケーション・プログラムの処理内容に応じて表示制御部２５０７、通信制御部２５１７、ドライブ装置２５１３を制御して、所定の動作を行わせる。また、処理途中のデータについては、主としてメモリ２５０１に格納されるが、ＨＤＤ２５０５に格納されるようにしてもよい。本技術の実施例では、上で述べた処理を実施するためのアプリケーション・プログラムはコンピュータ読み取り可能なリムーバブル・ディスク２５１１に格納されて頒布され、ドライブ装置２５１３からＨＤＤ２５０５にインストールされる。インターネットなどのネットワーク及び通信制御部２５１７を経由して、ＨＤＤ２５０５にインストールされる場合もある。このようなコンピュータ装置は、上で述べたＣＰＵ２５０３、メモリ２５０１などのハードウエアとＯＳ及びアプリケーション・プログラムなどのプログラムとが有機的に協働することにより、上で述べたような各種機能を実現する。
　以上述べた本実施の形態をまとめると以下のようになる。
　本実施の形態に係る形状データ生成方法は、（Ａ）画像データにより特定される変形目標の形状である目標形状に対して、当該目標形状と同じトポロジを有する単純な形状を有する入力形状を設定し、（Ｂ）入力形状の複数の頂点から、法線が目標形状と交点を有するという条件を含む所定の条件を満たす第１の頂点を特定し、（Ｃ）特定された第１の頂点を、当該第１の頂点についての法線の方向に、目標形状までの距離より短い第１の距離移動させるように、入力形状を変形し、（Ｄ）変形後の入力形状を処理対象の入力形状に変更して第１の頂点の特定及び入力形状の変更を所定回数実行する処理を含む。
　このような入力形状を用いることで、精度良く目標形状と同様の形状を有する形状データを生成できるようになる。
　なお、上で述べた入力形状の変形において、変形後の形状における複数のメッシュ要素のうち、面積が閾値を超えるメッシュ要素及び周辺のメッシュ要素についてメッシュ再分割処理を実行するようにしても良い。このようにメッシュ再分割処理を実施することで、おおよそ均一なサイズを有するメッシュ要素が生成されることになるので、滑らかな形状を生成できるようになる。特に、上で述べたような入力形状の場合には、不均一なメッシュ要素が生成されるケースが多くなるので、このような処理は精度の良い形状データを生成するために有効である。
　さらに、上で述べた入力形状の設定において、目標形状が三次元閉曲面であれば、球状の入力形状を、目標形状に合わせて縮小又は拡大して目標形状の重心位置に配置するようにしても良い。このようにすれば、三次元閉曲面であっても適切に変形処理を実行できるようになる。
　また、上で述べた入力形状の設定において、目標形状が開口部を有する三次元曲面であれば、平板の入力形状を、開口部に合わせて変形するようにしても良い。この場合、上で述べた入力形状の変形において、開口部に合わせて変形された入力形状のエッジ部分を固定させるようにしても良い。このようにすれば、開口部を有するような三次元曲面であっても適切に変形処理を実行できるようになる。
　さらに、本形状データ生成方法は、（Ｅ）第２の画像データにより特定され且つ目標形状に対して追加部分が付加された第２の目標形状に対して、第１の頂点の特定及び入力形状の変更を所定回数実行させた後の形状である第２の入力形状であって第２の目標形状の追加部分との境界面の外周エッジに移動させない固定点が設定された第２の入力形状を配置させ、（Ｆ）第２の入力形状の複数の頂点から、法線が第２の目標形状と交点を有するという条件を含む所定の条件を満たす第２の頂点を特定させ、（Ｇ）特定された第２の頂点を、当該第２の頂点についての法線の方向に、第２の目標形状までの距離より短い第２の距離移動させるように、固定点を除き第２の入力形状を変形させ、（Ｈ）変形後の第２の入力形状を処理対象の第２の入力形状に変更して第２の頂点の特定及び第２の入力形状の変更を第２の所定回数実行させるようにしても良い。
　複雑な形状を有する形状データを生成する場合には、２段階で変形処理を実行すれば精度の良い形状データを生成できるようになる。
　さらに、上で述べた第１の頂点の特定又は第２の頂点の特定において、（Ｘ１）着目した頂点を当該着目した頂点についての法線の方向へ第３の距離移動させ、移動先の点が前記画像データ又は前記第２の画像データにより特定されるボクセル空間に含まれるかを判断し、（Ｘ２）移動先の点がボクセル空間に含まれると判断された場合、着目した頂点についての法線ベクトルと移動先の点についての法線ベクトルとの内積に基づき、移動先の点が入力形状又は第２の入力形状を貫通したか判断し、（Ｘ３）移動先の点が入力形状又は第２の入力形状を貫通していないと判断された場合、移動先の点における輝度値と着目した頂点における輝度値とを比較し、輝度値が変化したか判断し、（Ｘ４）輝度値が変化したと判断された場合、着目した頂点についての法線が目標形状又は第２の目標形状と交点を有するという条件を満たすと判定し、（Ｘ５）移動先の点がボクセル空間に含まれないと判断された場合又は移動先の点が入力形状又は第２の入力形状を貫通したと判断された場合、着目した頂点についての法線が入力形状又は第２の入力形状と交点を有するという条件を満たさないと判定し、（Ｘ６）輝度値が変化していないと判断された場合、移動先の点が画像データ又は第２の画像データにより特定されるボクセル空間に含まれるかの判断以降の処理を移動先の点について再度実行するようにしても良い。
　このようにすれば、第１の頂点についての法線が第２の形状と交点を有するか否かを適切に判定することができるようになる。
　なお、上記方法による処理をコンピュータに行わせるためのプログラムを作成することができ、当該プログラムは、例えばフレキシブルディスク、ＣＤ－ＲＯＭ、光磁気ディスク、半導体メモリ、ハードディスク等のコンピュータ読み取り可能な記憶媒体又は記憶装置に格納される。尚、中間的な処理結果はメインメモリ等の記憶装置に一時保管される。

Claims (8)

  1. 　コンピュータに、
    　画像データにより特定される変形目標の形状である目標形状に対して、当該目標形状と同じトポロジを有する単純な形状を有する入力形状を設定させ、
    　前記入力形状の複数の頂点から、法線が前記目標形状と交点を有するという条件を含む所定の条件を満たす第１の頂点を特定させ、
    　特定された前記第１の頂点を、当該第１の頂点についての法線の方向に、前記目標形状までの距離より短い第１の距離移動させるように、前記入力形状を変形させ、
    　変形後の入力形状を処理対象の入力形状に変更して前記第１の頂点の特定及び前記入力形状の変更を所定回数実行させる
    　形状データ生成プログラム。
  2. 　前記入力形状の変形において、
    　変形後の形状における複数のメッシュ要素のうち、面積が閾値を超えるメッシュ要素及び周辺のメッシュ要素についてメッシュ再分割処理を実行する
    　請求項１記載の形状データ生成プログラム。
  3. 　前記入力形状の設定において、
    　前記目標形状が三次元閉曲面であれば、球状の入力形状を、前記目標形状に合わせて縮小又は拡大して前記目標形状の重心位置に配置する
    　請求項１又は２記載の形状データ生成プログラム。
  4. 　前記入力形状の設定において、
    　前記目標形状が開口部を有する三次元曲面であれば、平板の入力形状を、前記開口部に合わせて変形し、
    　前記入力形状の変形において、
    　前記開口部に合わせて変形された前記入力形状のエッジ部分を固定させる
    　請求項１又は２記載の形状データ生成プログラム。
  5. 　前記コンピュータに、さらに、
    　第２の画像データにより特定され且つ前記目標形状に対して追加部分が付加された第２の目標形状に対して、前記第１の頂点の特定及び前記入力形状の変更を所定回数実行させた後の形状である第２の入力形状であって前記第２の目標形状の前記追加部分との境界面の外周エッジに移動させない固定点が設定された前記第２の入力形状を配置させ、
    　前記第２の入力形状の複数の頂点から、法線が前記第２の目標形状と交点を有するという条件を含む前記所定の条件を満たす第２の頂点を特定させ、
    　特定された前記第２の頂点を、当該第２の頂点についての法線の方向に、前記第２の目標形状までの距離より短い第２の距離移動させるように、前記固定点を除き前記第２の入力形状を変形させ、
    　変形後の第２の入力形状を処理対象の第２の入力形状に変更して前記第２の頂点の特定及び前記第２の入力形状の変更を第２の所定回数実行させる
    　請求項１乃至４のいずれか１つ記載の形状データ生成プログラム。
  6. 　前記第１の頂点の特定又は前記第２の頂点の特定において、
    　着目した頂点を当該着目した頂点についての法線の方向へ第３の距離移動させ、移動先の点が前記画像データ又は前記第２の画像データにより特定されるボクセル空間に含まれるかを判断し、
    　前記移動先の点が前記ボクセル空間に含まれると判断された場合、前記着目した頂点についての法線ベクトルと前記移動先の点についての法線ベクトルとの内積に基づき、前記移動先の点が前記入力形状又は前記第２の入力形状を貫通したか判断し、
    　前記移動先の点が前記入力形状又は前記第２の入力形状を貫通していないと判断された場合、前記移動先の点における輝度値と前記着目した頂点における輝度値とを比較し、輝度値が変化したか判断し、
    　前記輝度値が変化したと判断された場合、前記着目した頂点についての法線が前記目標形状又は前記第２の目標形状と交点を有するという条件を満たすと判定し、
    　前記移動先の点が前記ボクセル空間に含まれないと判断された場合又は前記移動先の点が前記入力形状又は前記第２の入力形状を貫通したと判断された場合、前記着目した頂点についての法線が前記入力形状又は前記第２の入力形状と交点を有するという条件を満たさないと判定し、
    　前記輝度値が変化していないと判断された場合、前記移動先の点が前記画像データ又は前記第２の画像データにより特定されるボクセル空間に含まれるかの判断以降の処理を前記移動先の点について再度実行する
    　請求項１乃至５のいずれか１つ記載の形状データ生成方法。
  7. 　コンピュータが、
    　画像データにより特定される変形目標の形状である目標形状に対して、当該目標形状と同じトポロジを有する単純な形状を有する入力形状を設定し、
    　前記入力形状の複数の頂点から、法線が前記目標形状と交点を有するという条件を含む所定の条件を満たす第１の頂点を特定し、
    　特定された前記第１の頂点を、当該第１の頂点についての法線の方向に、前記目標形状までの距離より短い第１の距離移動させるように、前記入力形状を変形し、
    　変形後の入力形状を処理対象の入力形状に変更して前記第１の頂点の特定及び前記入力形状の変更を所定回数実行する
    　処理を実行する形状データ生成方法。
  8. 　画像データにより特定される変形目標の形状である目標形状に対して、当該目標形状と同じトポロジを有する単純な形状を有する入力形状を設定する設定部と、
    　前記入力形状の複数の頂点から、法線が前記目標形状と交点を有するという条件を含む所定の条件を満たす第１の頂点を特定する特定部と、
    　特定された前記第１の頂点を、当該第１の頂点についての法線の方向に、前記目標形状までの距離より短い第１の距離移動させるように、前記入力形状を変形する変形部と、
    　変形後の入力形状を処理対象の入力形状に変更して前記特定部及び前記変形部を所定回数動作させる処理部と、
    　を有する形状データ生成装置。
PCT/JP2012/068634 2012-07-23 2012-07-23 形状データ生成方法及び装置 WO2014016894A1 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2014526631A JP5850380B2 (ja) 2012-07-23 2012-07-23 形状データ生成方法及び装置
EP12881557.8A EP2876607B1 (en) 2012-07-23 2012-07-23 Shape data generation method and device
PCT/JP2012/068634 WO2014016894A1 (ja) 2012-07-23 2012-07-23 形状データ生成方法及び装置
US14/604,258 US9607423B2 (en) 2012-07-23 2015-01-23 Shape data generation method and apparatus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2012/068634 WO2014016894A1 (ja) 2012-07-23 2012-07-23 形状データ生成方法及び装置

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US14/604,258 Continuation US9607423B2 (en) 2012-07-23 2015-01-23 Shape data generation method and apparatus

Publications (1)

Publication Number Publication Date
WO2014016894A1 true WO2014016894A1 (ja) 2014-01-30

Family

ID=49996731

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2012/068634 WO2014016894A1 (ja) 2012-07-23 2012-07-23 形状データ生成方法及び装置

Country Status (4)

Country Link
US (1) US9607423B2 (ja)
EP (1) EP2876607B1 (ja)
JP (1) JP5850380B2 (ja)
WO (1) WO2014016894A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2018501549A (ja) * 2014-11-18 2018-01-18 セント・ジュード・メディカル，カーディオロジー・ディヴィジョン，インコーポレイテッド 幾何学的構造のパッチ表面モデルを生成するための方法およびシステム

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015092842A1 (ja) * 2013-12-16 2015-06-25 株式会社日立製作所 解析用メッシュデータ生成方法および解析用メッシュデータ生成装置
JP6762305B2 (ja) * 2015-01-28 2020-09-30 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． 解剖学的構造の有限要素モデリング
JP6536169B2 (ja) * 2015-05-21 2019-07-03 カシオ計算機株式会社 画像処理装置、画像処理方法及びプログラム

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001034774A (ja) 1999-07-21 2001-02-09 Ricoh Co Ltd 自由曲線の近似変換方法、装置及び記録媒体
JP2002329216A (ja) 2001-03-09 2002-11-15 Koninkl Philips Electronics Nv 対象物に含まれる三次元画像をセグメント化する方法
JP2007098028A (ja) 2005-10-07 2007-04-19 Konica Minolta Medical & Graphic Inc モデリング装置、モデリング方法、領域抽出装置、およびプログラム
JP2007159927A (ja) * 2005-12-15 2007-06-28 Konica Minolta Medical & Graphic Inc モデリング装置、モデリング方法、領域抽出装置、領域抽出方法及びプログラム
JP2010061431A (ja) 2008-09-04 2010-03-18 National Institute Of Advanced Industrial Science & Technology 衣類状態推定方法及びプログラム
JP2011224143A (ja) * 2010-04-20 2011-11-10 Hitachi Aloka Medical Ltd 超音波データ処理装置

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004536374A (ja) * 2001-03-09 2004-12-02 コーニンクレッカ　フィリップス　エレクトロニクス　エヌ　ヴィ 画像セグメンテーション
US7221787B2 (en) * 2002-12-10 2007-05-22 Eastman Kodak Company Method for automated analysis of digital chest radiographs
WO2004110309A2 (en) * 2003-06-11 2004-12-23 Case Western Reserve University Computer-aided-design of skeletal implants
US7602965B2 (en) * 2004-10-28 2009-10-13 Siemens Medical Solutions Usa, Inc. Object detection using cross-section analysis
US7620226B2 (en) * 2004-10-28 2009-11-17 Siemens Medical Solutions Usa, Inc. Semi-local active contour for a faint region detection in knee CT images
US7557804B1 (en) * 2006-03-06 2009-07-07 Adobe Systems Inc. Methods and apparatus for three-dimensional isographic navigation
JP4994737B2 (ja) * 2006-07-27 2012-08-08 オリンパスメディカルシステムズ株式会社 医療用画像処理装置及び医療用画像処理方法
WO2008063494A2 (en) * 2006-11-16 2008-05-29 Vanderbilt University Apparatus and methods of compensating for organ deformation, registration of internal structures to images, and applications of same
US8160345B2 (en) * 2008-04-30 2012-04-17 Otismed Corporation System and method for image segmentation in generating computer models of a joint to undergo arthroplasty
US8108189B2 (en) * 2008-03-25 2012-01-31 Align Technologies, Inc. Reconstruction of non-visible part of tooth
JP4806706B2 (ja) * 2008-12-18 2011-11-02 株式会社日立製作所 解析モデルの変形方法および計算機
JP5221436B2 (ja) * 2009-04-02 2013-06-26 トヨタ自動車株式会社 顔特徴点検出装置及びプログラム
JP5546230B2 (ja) * 2009-12-10 2014-07-09 キヤノン株式会社 情報処理装置、情報処理方法、及びプログラム
US20120253170A1 (en) * 2011-03-29 2012-10-04 Samsung Electronics Co., Ltd. Method and apparatus for generating medical image of body organ by using 3-d model

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001034774A (ja) 1999-07-21 2001-02-09 Ricoh Co Ltd 自由曲線の近似変換方法、装置及び記録媒体
JP2002329216A (ja) 2001-03-09 2002-11-15 Koninkl Philips Electronics Nv 対象物に含まれる三次元画像をセグメント化する方法
JP2007098028A (ja) 2005-10-07 2007-04-19 Konica Minolta Medical & Graphic Inc モデリング装置、モデリング方法、領域抽出装置、およびプログラム
JP2007159927A (ja) * 2005-12-15 2007-06-28 Konica Minolta Medical & Graphic Inc モデリング装置、モデリング方法、領域抽出装置、領域抽出方法及びプログラム
JP2010061431A (ja) 2008-09-04 2010-03-18 National Institute Of Advanced Industrial Science & Technology 衣類状態推定方法及びプログラム
JP2011224143A (ja) * 2010-04-20 2011-11-10 Hitachi Aloka Medical Ltd 超音波データ処理装置

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
0. SORKINE; D. COHEN-OR; Y. LIPMAN; M. ALEXA; C. ROESSL; H.-P. SEIDEL: "Laplacian surface editing", SGP '04 PROCEEDINGS OF THE 2004 EUROGRAPHICS/ACM SIGGRAPH SYMPOSIUM ON GEOMETRY PROCESSING
FRED L. BOOKSTEIN: "Principal Warps: Thin-Plate Splines and the Decomposition of Deformations", IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE, vol. 11, no. 6, June 1989 (1989-06-01)
FRED L. BOOKSTEIN: "Principal Warps: Thin-Plate Splines and the Decomposition of Deformations", IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE, vol. 11, no. 6, June 1989 (1989-06-01), pages 567 - 585
See also references of EP2876607A4

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2018501549A (ja) * 2014-11-18 2018-01-18 セント・ジュード・メディカル，カーディオロジー・ディヴィジョン，インコーポレイテッド 幾何学的構造のパッチ表面モデルを生成するための方法およびシステム
US10861246B2 (en) 2014-11-18 2020-12-08 St. Jude Medical, Cardiology Division, Inc. Methods and systems for generating a patch surface model of a geometric structure

Also Published As

Publication number Publication date
US9607423B2 (en) 2017-03-28
JPWO2014016894A1 (ja) 2016-07-07
EP2876607A1 (en) 2015-05-27
EP2876607B1 (en) 2019-07-31
US20150131914A1 (en) 2015-05-14
EP2876607A4 (en) 2016-01-06
JP5850380B2 (ja) 2016-02-03

Similar Documents

Publication Publication Date Title
CN104794758B (zh) 一种三维图像的裁剪方法
JP5120926B2 (ja) 画像処理装置、画像処理方法およびプログラム
JP5866177B2 (ja) 画像処理装置および画像処理方法
US8493389B2 (en) 3D connected shadow mouse pointer
EP2080170B1 (en) Combined intensity projection
CN110168619B (zh) 在虚拟现实或增强现实环境中渲染复杂数据的系统和方法
EP2402910B1 (en) Seamless fracture generation in a graphic pipeline
JP7038377B2 (ja) 生体モデル生成装置、生体モデル生成方法、および生体モデル生成プログラム
US20040036674A1 (en) Apparatus and method for associating voxel information with display positions
JP5850380B2 (ja) 形状データ生成方法及び装置
US20150130804A1 (en) Shape data generation method and apparatus
RU2604674C2 (ru) Системы и способы для создания трехмерного текстурного атласа
JP5721225B2 (ja) 形状データ生成方法、プログラム及び装置
CN109509259B (zh) 一种医学图像三维重建等值面网格提取方法
CN108090953B (zh) 感兴趣区域重建方法、系统以及计算机可读存储介质
EP3726477A1 (en) Chamber reconstruction from a partial volume
Glueck et al. Multiscale 3D reference visualization
US11776207B2 (en) Three-dimensional shape data processing apparatus and non-transitory computer readable medium
JP6106293B2 (ja) 画像処理装置、画像処理システム、および画像処理方法
JP2022171068A (ja) データ処理装置、データ処理方法、及びプログラム
CN117152378A (zh) 血管的建模方法、装置、计算机设备和存储介质
JP2019101966A (ja) 画像生成装置および画像生成方法
Yang et al. A Feature Preserved Mesh Subdivision Framework for Biomedical Mesh
Kedia et al. Layer based 3D clipping
Chen Volume Visualization

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 12881557

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2014526631

Country of ref document: JP

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 2012881557

Country of ref document: EP

NENP Non-entry into the national phase

Ref country code: DE