WO2006135024A1

WO2006135024A1 - Image processing method, image processing device, and image processing program

Info

Publication number: WO2006135024A1
Application number: PCT/JP2006/312066
Authority: WO
Inventors: Yoshinori Arai
Original assignee: Matsumoto Dental University
Priority date: 2005-06-16
Filing date: 2006-06-15
Publication date: 2006-12-21
Also published as: JP2006350724A; JP4710030B2

Abstract

Provided are an image processing method, an image processing device and an image processing program for displaying a two-dimensional image cut out from three-dimensional image data. The image processing method, the image processing device and the image processing program can intuitively recognize and precisely set the positional relation between the first three-dimensional image data and the second three-dimensional image data. There are displayed the two-dimensional image cut out from the first three-dimensional image data, the two-dimensional image cut out from the second three-dimensional data, and a comparison image, in which the two-dimensional image cut out from the first three-dimensional image data and the two-dimensional image cut out from the second three-dimensional image data are overlapped. The relational position between the first three-dimensional image data and the second three-dimensional image data is set such that the first three-dimensional image data and the second three-dimensional image data are substantially identical on the basis of the two-dimensional image cut out from the first three-dimensional image data displayed, the two-dimensional image cut out from the second three-dimensional image data, and the comparison image.

Description

Image processing method, image processing apparatus, and image processing program

The present invention relates to an image processing method, an image processing apparatus, and an image processing program, and more particularly to an image processing method, an image processing apparatus, and an image processing program for displaying a 2D image cut out from 3D image data. .

Background art

[0002] In the medical and medical fields, it is desired that one individual can be observed over a long period of time in order to eliminate individual differences among experimental animals and perform experiments with as few individuals as possible. In particular, when observing the internal state of an individual, it is currently the case that a large number of equivalent individuals are bred under the same conditions and observed using different solids for each observation. It was impossible to observe the individual for a long time.

[0003] Micro X-ray CT has been developed to observe the internal structure of an individual without damaging the individual. Micro X-ray CT captures a number of two-dimensional images while an X-ray generator and an X-ray imaging unit arranged opposite to the periphery of the imaging target rotate around the imaging target. By constructing a three-dimensional image from the above, a three-dimensional image of the internal structure of the imaging target can be acquired.

[0004] A three-dimensional image acquired by micro X-ray CT can be cut at a desired plane to obtain a necessary cross section. Therefore, it is possible to freely observe the internal structure of the imaging target. For this reason, it becomes possible to observe one individual over a long period of time.

Disclosure of the invention

Problems to be solved by the invention

[0005] However, three-dimensional images captured at different timings by micro X-ray CT naturally have different coordinate positions. When observing a single individual, it is necessary to observe images such as cut surfaces of the same part side by side. However, since the coordinate positions are different, it was difficult to display the cut surface of the same part. In addition, since the extraction of the cut surface relies on the observer's feeling, it is applied with accuracy. [0006] In addition, the processing for positioning the positions of the two three-dimensional images requires a large-sized, dedicated computer and the like, and it is difficult to easily introduce the processing.

[0007] Therefore, it has been desired to be able to set the positions of three-dimensional images of two identical individuals with different imaging timings with simple processing and accurately.

[0008] The present invention has been made in view of the above points, and image processing capable of intuitively recognizing and accurately setting the positional relationship between the first three-dimensional image data and the second three-dimensional image data. It is an object to provide a method, an image processing apparatus, and an image processing program.

Means for solving the problem

[0009] The present invention provides a two-dimensional image cut out from the first three-dimensional image data, a two-dimensional image cut out from the second three-dimensional image data, and the first three-dimensional image data. A two-dimensional image cut out from the displayed first three-dimensional image data is displayed by displaying the cut-out two-dimensional image and a comparison image obtained by superimposing the two-dimensional image cut out from the second three-dimensional image data. The first 3D image data and the second 3D image data are substantially matched based on the image, the 2D image cut out from the second 3D image data, and the comparison image. It is characterized in that the relative position between the 3D image data and the second 3D image data is set.

[0010] Further, in the comparison image, a two-dimensional image cut out from the first three-dimensional image and a two-dimensional image cut out from the second three-dimensional image are distinguished from each other in a portion that is the same and a portion that is different from each other. It is displayed so that it can be performed.

[0011] Further, the comparison image is different between the two-dimensional image cut out from the first three-dimensional image and the two-dimensional image cut out from the second three-dimensional image. It is characterized in that the colors are displayed differently for each part.

[0012] In the present invention, the two-dimensional image cut out from the second three-dimensional image and the portion where the values of the two-dimensional image cut out from the two-dimensional image cut out from the first three-dimensional image coincide with each other. It is characterized by converging the positions of the first 3D image and the second 3D original image so as to increase.

The invention's effect

[0013] According to the present invention, a two-dimensional image cut out from the first three-dimensional image data, 2D image cut out from 2D 3D image data, 2D image cut out from 1st 3D image data, and 2D image cut out from 2D 3D image data 2D image cut out from the displayed first 3D image data, 2D image cut out from the second 3D image data, and comparison image Based on the above, the relationship position between the first 3D image data and the second 3D image data is set so that the first 3D image data and the second 3D image data substantially coincide with each other. Thus, the positional relationship between the first three-dimensional image data and the second three-dimensional image data can be intuitively recognized by the comparison image, and the processing can be easily performed by comparing the two-dimensional images. Exactly, first 3D image data and second It is possible to set the relationship between the position of the three-dimensional image data.

Brief Description of Drawings

FIG. 1 is a system configuration diagram of an embodiment of the present invention.

2 is a block configuration diagram of the image processing apparatus 112. FIG.

3 is an image processing flowchart of the image processing apparatus 112. FIG.

FIG. 4 is a diagram showing an example of a comparative image screen.

FIG. 5 is a diagram showing a screen displayed on display device 133.

FIG. 6 is a diagram showing a screen displayed on display device 133.

FIG. 7 is a diagram showing a screen displayed on display device 133.

FIG. 8 is a diagram showing a screen displayed on display device 133.

FIG. 9 is an explanatory diagram of repositioning operation.

FIG. 10 shows a flowchart of the repositioning process.

FIG. 11 is an axis repositioning process flowchart in the repositioning process. Explanation of symbols

[0015] 100 X-ray CT, 112 image processing device

121 stages, 122 imaging holes, 123 X-ray generator

124 X-ray imaging unit, 125 rotation drive mechanism, 126 control unit, 127 communication unit

131 Computer, 132 Input device, 133 Display device

141 interface, 142 CPU, 143 hard disk drive 144 replaceable disk drive, 145 memory

146 graphics controller, 147 communication unit

150 Display screen

151, 154 Z— X plane image, 152, 155 Z— Y plane image

153, 156 X—Y plane image, 157 Z—X plane comparison image

158 Z—Y side comparison image

159 X-Y comparison image, 160, 161, 162 grayscale image

BEST MODE FOR CARRYING OUT THE INVENTION

[0016] [Configuration]

FIG. 1 shows a system configuration diagram of an embodiment of the present invention.

The system 100 according to the present embodiment includes an X-ray CT 111 and an image processing device 112.

The X-ray CT 111 includes a stage 121, an imaging hole 122, an X-ray generation unit 123, an X-ray imaging unit 124, a rotation drive mechanism 125, a control unit 126, and a communication unit 127.

[0019] On the stage 121, an imaging target 128 is mounted. The imaging target 128 is inserted into the imaging hole 122 by the stage 121. The X-ray generator 123 generates X-ray Rx and radiates it in the direction of the imaging hole 122. The X-ray imaging unit 124 is provided opposite to the X-ray generation unit 123 with the imaging hole 122 interposed therebetween, and is a two-dimensional image of the X-ray Rx emitted from the X generation unit 123 and transmitted through the imaging target 128. Image.

[0020] The X-ray generation unit 123 and the X-ray imaging unit 124 are disposed so as to face each other with the imaging hole 122 interposed therebetween, and sequentially rotate while rotating around the imaging target 128 by the rotation drive mechanism 125. The captured two-dimensional image is supplied to the image processing apparatus 112 via the control unit 126 and the communication unit 127. The image processing apparatus 112 constructs a three-dimensional image from a plurality of two-dimensional images captured by the X-ray CT 111 and supplied.

[0021] The control unit 126 controls the X-ray generation unit 123, the X-ray imaging unit 124, and the rotation drive mechanism 125 to image the imaging target 128 from a predetermined angle. The communication unit 127 communicates with the image processing device 112 and transmits the captured two-dimensional image to the image processing device 112.

[Image processing apparatus 112] FIG. 2 is a block diagram of the image processing apparatus 112.

The image processing apparatus 112 is configured with a computer system such as a personal computer, and mainly includes a computer main body 131, an input device 132, and a display device 133.

The computer main body 131 includes an interface 141, a CPU 142, a hard disk drive device 143, a replaceable disk drive device 144, a memory 145, a graphics controller 146, and a communication unit 147.

The interface 141 is equipped with a USB interface or the like, and interfaces with the input device 132 or the like. The CPU 142 performs image analysis processing based on the image analysis program installed in the hard disk drive device 143.

The hard disk drive device 143 stores an image analysis program and image data including first 3D image data and second 3D image data to be analyzed.

[0027] In the present embodiment, the first 3D image data and the second 3D image data are 3D image data of the same imaging target with different imaging dates and times.

[0028] The removable disk drive device 144 is a device that performs recording and Z or reproduction of a removable disk such as a CD-ROM, DVD-ROM, etc., and an image analysis program installed in the hard disk drive device 143. Is provided by a replaceable disc that can be played back by the replaceable disc drive device 144, for example.

The memory 145 is composed of a semiconductor storage device and is used as a working storage area for the CPU 142. The graphics controller 146 generates an image signal for displaying an image on the display device 133 from the image data supplied from the CPU 142, and supplies the image signal to the display device 133.

[0030] The communication unit 147 communicates with the X-ray CT111. The display device 113 displays an image corresponding to the image signal supplied from the graphics controller 146.

The image processing device 112 constructs a three-dimensional image from a plurality of two-dimensional images supplied from the X-ray CT 111 to the communication unit 147 and stores the three-dimensional image in the hard disk drive device 147. Note that a 3D image is constructed from a plurality of 2D images on the X-ray CT 111 side and sent to the image processing device 112. Make sure to believe.

[Image processing]

FIG. 3 shows an image processing flowchart of the image processing apparatus 112.

[0033] Here, an image obtained by imaging an animal bone by X-ray CT111 will be described.

[0034] The CPU 142 reads out the first three-dimensional image and the second three-dimensional image stored in the hard disk drive 147 in step S1-1, and in step S1-2, the first three-dimensional image and the second three-dimensional image. A two-dimensional image is cut out from the original image. The two-dimensional image to be cut out is the two-dimensional image force of the three surfaces of Z-X plane, Z-Y plane, and XY plane orthogonal to each other of the first three-dimensional image and the second three-dimensional image. Constructed!

[0035] In step S1-3, the CPU 142 extracts a Z-X plane image, a Z-Y plane image, a three-dimensional two-dimensional image of the X-Y plane, and a second plane clipped from the first three-dimensional image. The Z-X plane image, Z-Y plane image, and X-Y plane image cut out from the three-dimensional image of the above are binary-valued. In the binarized image, a white display indicates a site where a bone exists, and a black display indicates a site where no bone exists.

[0036] In step S1-4, the CPU 142 determines that the binary image of the Z-X plane image cut out from the first three-dimensional image and the Z-X plane cut out from the second three-dimensional image camera. The binary image of the image is compared, and a comparison image corresponding to the comparison result and its gray scale image are generated and displayed.

FIG. 4 is a diagram showing an example of a comparative image.

[0038] The comparison image includes a binary image of the Z-X plane image cut out from the first three-dimensional image and a binary image of the Z-X plane image cut out from the second three-dimensional image force. The image of the matching part is displayed in green. In addition, the comparison image is a Z-X plane image cut out from the first three-dimensional image, that is, a portion that was white, that is, a portion where bones were present, was cut out from the second three-dimensional image. If the surface image is black, that is, if there are no bones, that portion is displayed in red. In addition, the comparison image is the black part of the Z-X plane image cut out from the first three-dimensional image, that is, the part where the bone was not present and the second three-dimensional image force was also cut out. If the Z-X plane image is white, that is, if a bone is present, that portion is displayed in blue. 5 to 8 are diagrams showing screens displayed on the display device 133. FIG. Figure 4 shows Z-X plane image 151, Z-Y plane image 152, X-Y plane image 153, and second Z-dimensional image Z extracted from the first three-dimensional image. — X plane image 154, Z— Y plane image 155, X—Y plane image 156 in grayscale display, Figure 5 shows the Z—X plane image 151, Z—Y plane cut out from the first 3D image Image 152, X—Y plane image 153, and second 3D image power Z-X plane image 154, Z—Y plane image 155, and X—Y plane image 156 binarized images Fig. 6 shows the screen before X-Y plane image 156 is displaced in the Z direction, and Fig. 7 shows the screen in Fig. 6 after X-Y plane image 156 is displaced in the Z direction. It shows.

[0040] The display screen 150 displayed on the screen of the display device 133 includes a Z-X plane image 151, a Z-Y plane image 152, an XY plane image 153, and the like extracted from the first three-dimensional image. Z-X image 154, Z-Y image 155, X-Y image 156, Z-X image 157, Z-Y image 158, X-Y comparison image 159, ZX comparison image 157 grayscale image 160, Z—Y comparison image 158 grayscale image 161, X—Y comparison image 159 grayscale image 162 Yes.

[0041] At this time, the Z-X plane image 151, the Z-Y plane image 152, the XY plane image 153, and the second three-dimensional image force that are also cut out from the first three-dimensional image force. The display position of each of the —X plane image 154, the Z —Y plane image 155, and the X—Y plane image 156 can be manipulated by operating the input device 132. At this time, if the position of one of the first three-dimensional image forces extracted Z-X plane image 151, Z-Y plane image 152, and X-Y plane image 153 is adjusted, the other two images The position of is also adjusted according to the position of the adjusted single plane image. Similarly, adjusting the position of one of the Z-X plane image 154, Z-Y plane image 155, and X-Y plane image 156 from which the second three-dimensional image force is also extracted, the other two The position of the image is also adjusted according to the position of the adjusted single plane image.

[0042] Further, according to the displacement of the position of the adjusted image, Z-X plane comparison image 157, Z-Y plane comparison image 158, X-Y plane comparison image 159, and Z-X plane comparison image 157 grayscale image 160, Z—Y comparison image 158 grayscale image 161, X—Y comparison image 1 59 grayscale images 162 also change.

[0043] The CPU 142 operates the input device 132 in step S1-5 to operate the Z-X plane image 151, the Z-Y plane image 152, the X-Y plane image 153, which are cut out from the first three-dimensional image. When the display position of each of the Z-X plane image 154, Z-Y plane image 155, and X-Y plane image 156, from which the second three-dimensional image force is cut out, is manipulated in step S1-6. As described above, the position of the adjusted plane image of the three-dimensional image is adjusted, and the display of the image changes with the displacement of the position.

[0044] In addition, when repositioning is instructed by operating the input device 132 in step S1-7, the CPU 142 selects a Z-X plane image 151 cut out from the first three-dimensional image in step S1-8. Z—X plane image 154 cut out from the second 3D image, Z—Y plane image 152 cut out from the first 3D image, and Z—cut out from the second 3D image Y-plane image 155, X-Y plane image 153 cut out from the first three-dimensional image, and second cubic so that the X-Y plane image 156 cut out from the second three-dimensional image matches. The position of the original image is adjusted. At this time, for example, the CPU 142 adjusts the positions of the Z-X plane image 154, the Z-Y plane image 155, and the XY plane image 156 of the second three-dimensional image, respectively, and compares the Z-X plane comparison image. 157, Z—Y surface comparison image 158, X—Y surface comparison image 159 converge to a position where the green part is as much as possible.

[0045] The position of the second three-dimensional image can be made to substantially coincide with the first three-dimensional image by repositioning.

Here, the repositioning process will be described.

FIG. 9 is a diagram for explaining the operation of repositioning. FIG. 9A shows a first three-dimensional image, and FIG. 9B shows a second three-dimensional image.

In this embodiment, the first and second 3D images are as shown in FIG. 9A. The first 3D image is an area of XXYXZ = 512 dots × 512 dots × 384 dots Al 1 force. Among them, the area A12 of XXYX Z = 200 dots X 200 dots X 200 dots is used for repositioning. In addition, as shown in Fig. 9 (B), the second 3D image is composed of A 51 = 51 2 dots x 512 dots x 384 dots area A21 force, of which ΧΧΥΧ Ζ = 200 dots x 200 dots X 200 dot area Α22 is used for repositioning. Lipo In shifting, find the Y, X, and Y surfaces of area A12 and the YX surface of area Α22, and maximize the difference between the corresponding surfaces. In FIG. 9, 01 and 02 indicate imaging targets.

FIG. 10 shows a flowchart of the surface repositioning process in the repositioning process. Here, repositioning of the heel and Z 'planes is explained. Here, the saddle surface indicates the X-plane of the region A12 set in the first three-dimensional image, and the Z 'plane indicates the X-plane of the region Α22 set in the second three-dimensional image. . S (Z,) is the coincidence area between the Z plane and the plane, S (Z'-l) is the coincidence area between the Z plane and the plane when shifted by 1 dot in the Z direction, and S (Z '+ l ) Shows the coincidence area between the Z plane and the plane when shifted by 1 dot in the + Z direction.

[0050] CPU 142 determines in step S2-1 whether S (Ζ'-1) is maximum. If S (Ζ'-1) is the maximum in step S2-1, CPU 142 sets Ζ '= (Ζ'-1) in step S2-2 and returns to the processing in step S2-1. On the other hand, when S (Z, -1) is not the maximum in step S2-1, the CPU 142 determines whether S (Z '+ l) is the maximum in step S2-3.

If S (Z ′ + l) is the maximum in step S2—3, CPU 142 sets Z ′ = (Z ′ + 1) in step S2—4 and returns to the process of step S2-1. In addition, CPU 142 can determine that S (Z ') is the maximum when S (Z' + l) is not the maximum in step S2-3, so the position in the Z direction at that time is determined on the Z plane. Match surface.

[0052] Repositioning is performed by performing the above process on the X plane corresponding to the Y-Z plane and the Y plane corresponding to the X-Z plane.

In addition, the origin is moved and the axis is converted by repositioning.

The axis conversion will be described.

[0055] FIG. 11 shows a flowchart of the axis repositioning process in the repositioning process.

[0056] The total area of the X—Y, X—Z, and Y—Z planes of region A12 and region A22 when the angle of the Z-plane axis is p and degrees is S (p,), and the axis angle Is the total coincidence area of the X — Y plane, X — Z plane, and Y — Z plane of region A12 and region A22 when is (p, 1 1) degrees, and the axis angle is (ρ Let S (p '+ l) be the total coincidence area of Χ-Χ surface, Χ-Ζ surface, and Υ- Ζ surface of region A12 and region Α22 at' + 1) degrees. The CPU 142 determines whether or not S (p′-l) is maximum in step S3-1. If S (ρ′−1) is the maximum in step S3-1, CPU 142 sets ρ ′ = (ρ′−1) in step S3-2 and returns to the processing in step S3-1. On the other hand, if S (p, -1) is not the maximum in step S3-1, the CPU 142 determines whether S (p '+ l) is the maximum force in step S3-3.

If S (p ′ + l) is the maximum in step S3—3, CPU 142 sets p ′ = (p ′ + 1) in step S3—4 and returns to the process of step S3-1. Further, the CPU 142 can determine that S (p ′) is the maximum when S (p ′ + l) is not the maximum in step S3-3, and therefore the angle of the axis at that time is determined as the axis coincidence angle. To do. Repositioning is performed by performing the above processing on the axis q orthogonal to the X plane corresponding to the Y-Z plane and the axis r orthogonal to the Y plane corresponding to the X-Z plane.

[0059] By the above processing, S (Z '), S (X'), S (Y '), S (p'), S (q '), and S (r') are the maximum positions. The position of area A22 converges with respect to area A12. In addition, if any one area is the maximum when shifted by 1 dot or 1 degree, the calculation is performed again.

[0060] At this time, in this embodiment, the first three-dimensional image and the second three-dimensional image are obtained by the Z-X plane comparison image 157, the Z-Y plane comparison image 158, and the XY plane comparison image 159. The positional relationship between the two is intuitively easy to recognize. Therefore, while referring to the Z-X plane comparison image 157, the Z-Y plane comparison image 158, and the X-Y plane comparison image 159, the Z-X plane image 151, in which the first three-dimensional image force was manually cut out, Adjust the position of the Z—Y plane image 152, the X—Y plane image 153, or the Z—X plane image 154, the Z—Y plane image 155, and the X—Y plane image 156 of the second three-dimensional image. Alignment can be performed. This makes the repositioning process easy and accurate.

[0061] Furthermore, the position at which it is desired to manually observe the positions of the X-Y plane image 153 from which the first three-dimensional image force has been cut out and the X-Y plane image 156 from which the second three-dimensional image force has also been cut out is manually centered It can also be adjusted so that As a result, it is possible to obtain an image in which the image to be observed is recognized.

[0062] CPU 142 ends the process when the end operation is performed in step S1-9.

[0063] According to the present embodiment, the Z-X plane image 151, Z-Y plane cut out from the first three-dimensional image The two-dimensional image power of image 152, X—Y plane image 153, or Z—X plane image 154 of the second three-dimensional image, Z—Y plane image 155, and X—Y plane image 156 is also the Z—X plane comparison image 157. Z-Y comparison image 158, X-Y comparison image 159, Z-X comparison image 157 grayscale image 160, Z-Y comparison image 158 grayscale image 161, X-Y Since the grayscale image 162 of the surface comparison image 159 is generated and displayed, the processing in the image processing device 112 can be simplified.

[0064] Further, according to the present embodiment, since the same imaging target can be observed without damaging, the internal structure of the same individual can be observed over a long period of time. This allows observations without individual differences. In addition, the number of individuals required for observation can be reduced. In addition, this can reduce the speed and cost of the experiment.

[0065] Needless to say, the present invention is not limited to the above-described embodiments, and various modifications can be made without departing from the gist of the present invention!

[0066] In addition, this international application claims priority based on Japanese Patent Application 2005-176600 filed on June 16, 2005. The entire contents of 2005-176600 are claimed in this international application. Incorporated into.

Claims

The scope of the claims

[1] A two-dimensional image cut out from the first three-dimensional image data, a two-dimensional image cut out from the second three-dimensional image data, and a cut-out from the first three-dimensional image data And a comparison image obtained by superimposing the two-dimensional image cut out from the second three-dimensional image data,

The two-dimensional image cut out from the displayed first three-dimensional image data, the two-dimensional image cut out from the second three-dimensional image data, and the first image based on the comparison image A relational position between the first 3D image data and the second 3D image data is set so that the 3D image data and the second 3D image data substantially match. Image processing method.

[2] The comparison image includes a portion in which the two-dimensional image cut out from the first three-dimensional image and the two-dimensional image cut out from the second three-dimensional image match each other and each other! 2. The image processing method according to claim 1, wherein different parts are displayed so as to be distinguishable.

[3] The comparison image is the same as the two-dimensional image extracted from the first three-dimensional image and the two-dimensional image extracted from the second three-dimensional image. 3. The image processing method according to claim 2, wherein the image is displayed so that the color is different from that of the portion different from the eyelid.

[4] The two-dimensional image cut out from the second three-dimensional image and the two-dimensional image cut out from the first three-dimensional image so that there are many portions where the values of the two values coincide with each other. The image processing method according to claim 1, wherein the positions of the first three-dimensional image and the second three-dimensional original image are converged.

[5] An image processing apparatus for displaying a two-dimensional image cut out from the first three-dimensional image data and a two-dimensional image cut out from the second three-dimensional image data,

An input unit for instructing repositioning of the two-dimensional image;

2D image cut out from the first 3D image data, 2D image cut out from the second 3D image data, and 2D cut out from the first 3D image data Secondary extracted from the image and the second 3D image data A comparison image obtained by superimposing the original image is displayed, and the first three-dimensional image data and the second three-dimensional image data substantially match based on the comparison image in accordance with an instruction from the input unit. An image processing apparatus comprising a processing unit that sets a relative position between the first three-dimensional image data and the second three-dimensional image data.

[6] The processing unit may be configured to match the comparison image between a two-dimensional image cut out from the first three-dimensional image and a two-dimensional image cut out from the second three-dimensional image. 6. The image processing apparatus according to claim 5, wherein different parts are displayed so as to be distinguished from each other.

[7] The processing unit may be configured such that the comparison image includes a portion where a two-dimensional image cut out from the first three-dimensional image and a two-dimensional image cut out from the second three-dimensional image coincide with each other. 7. The image processing apparatus according to claim 6, wherein the display is performed so that the colors are different from each other in different parts.

[8] The processing unit binarizes a two-dimensional image cut out from the second three-dimensional image and a two-dimensional image cut out from the first three-dimensional image according to a repositioning instruction from the input unit. 6. The image according to claim 5, wherein the positions of the first three-dimensional image and the second three-dimensional original image are converged so that there are more portions where the values of the matched images match. Processing equipment.

[9] On the computer,

A two-dimensional image cut out from the first three-dimensional image data, a two-dimensional image cut out from the second three-dimensional image data, and two pieces cut out from the first three-dimensional image data Display a comparison image obtained by superimposing a two-dimensional image and a two-dimensional image cut out from the second three-dimensional image data,

Based on the two-dimensional image cut out from the first three-dimensional image data displayed on the display device, the two-dimensional image cut out from the second three-dimensional image data, and the comparison image A relational position between the first 3D image data and the second 3D image data is set so that the first 3D image data and the second 3D image data substantially coincide with each other. An image processing program readable by a computer.

[10] The comparison image is a portion in which the two-dimensional image cut out from the first three-dimensional image and the two-dimensional image cut out from the second three-dimensional image match each other! 10. The computer-readable image processing program according to claim 9, wherein different parts are displayed so as to be distinguishable.

[11] The comparison image is the same as a portion where the two-dimensional image cut out from the first three-dimensional image and the two-dimensional image cut out from the second three-dimensional image match each other! 11. The computer readable image processing program according to claim 10, wherein the image processing program is displayed so that a color is different from that of a portion different from the eyelid.

[12] The two-dimensional image cut out from the second three-dimensional image and the two-dimensional image cut out from the first three-dimensional image have more portions where the values of the two values coincide with each other. 10. The computer-readable image processing program according to claim 9, wherein the positions of the first three-dimensional image and the second three-dimensional original image are converged.