WO2004003850A1 - 3次元イメージ比較プログラム、3次元イメージ比較方法、および3次元イメージ比較装置 - Google Patents

3次元イメージ比較プログラム、3次元イメージ比較方法、および3次元イメージ比較装置 Download PDF

Info

Publication number
WO2004003850A1
WO2004003850A1 PCT/JP2002/006621 JP0206621W WO2004003850A1 WO 2004003850 A1 WO2004003850 A1 WO 2004003850A1 JP 0206621 W JP0206621 W JP 0206621W WO 2004003850 A1 WO2004003850 A1 WO 2004003850A1
Authority
WO
WIPO (PCT)
Prior art keywords
image
dimensional
difference
poxels
poxel
Prior art date
Application number
PCT/JP2002/006621
Other languages
English (en)
French (fr)
Inventor
Makoto Amakai
Original Assignee
Fujitsu Limited
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 Fujitsu Limited filed Critical Fujitsu Limited
Priority to PCT/JP2002/006621 priority Critical patent/WO2004003850A1/ja
Priority to KR1020047021408A priority patent/KR100639139B1/ko
Priority to EP02736203A priority patent/EP1519318A4/en
Priority to JP2004517224A priority patent/JP4001600B2/ja
Priority to CNB028292359A priority patent/CN100388317C/zh
Publication of WO2004003850A1 publication Critical patent/WO2004003850A1/ja
Priority to US10/989,464 priority patent/US20050068317A1/en

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0004Industrial image inspection
    • G06T7/001Industrial image inspection using an image reference approach
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T19/00Manipulating 3D models or images for computer graphics
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30108Industrial image inspection
    • G06T2207/30164Workpiece; Machine component
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2210/00Indexing scheme for image generation or computer graphics
    • G06T2210/36Level of detail

Definitions

  • 3D image comparison program 3D image comparison method, and 3D image comparison device
  • the present invention relates to a three-dimensional image comparison program, a three-dimensional image comparison method, and a three-dimensional image comparison device for comparing three-dimensional images and displaying differences, and in particular, to compare the three-dimensional image poxel differences with high accuracy.
  • the present invention relates to a three-dimensional image comparison program, a three-dimensional image comparison method, and a three-dimensional image comparison device.
  • the three-dimensional images include, for example, a design image by three-dimensional CAD (Computer Aided Design) in various manufacturing industries, and a three-dimensional image of a diseased part by a three-dimensional ultrasonic diagnostic imaging apparatus in the medical field.
  • CAD Computer Aided Design
  • a three-dimensional image of a diseased part by a three-dimensional ultrasonic diagnostic imaging apparatus in the medical field.
  • the image to be compared is a design model consisting of a three-dimensional free-form surface created by CAD, a product manufactured from the model, or point cloud information obtained by measuring a part or part with a three-dimensional shape measuring instrument, or The cut created from the point cloud information This is information about the screen.
  • the difference was evaluated based on the distance between a point on one surface and a point on the other surface.
  • These points can be measured with a computed tomography (CT) scanner or 3D digitizer, or obtained by dividing a given surface according to certain rules.
  • CT computed tomography
  • the distance between the surfaces is the perpendicular distance from each point to the surface. This distance information is plotted on a point on the surface, and the value is displayed, or the color density or the luminance is proportional to the distance, and the image is rendered visible. Therefore, a specific description will be given below with reference to FIG.
  • FIG. 26 is a diagram showing a three-dimensional image (CAD data: three-dimensional curved surface) and a three-dimensional point cloud (measurement data) of the object.
  • CAD data three-dimensional curved surface
  • measurement data measurement data
  • the distance between a certain point (x., Y., Z.) Is the minimum of ⁇ (xx 0 ) 2 + (yy 0 ) 2 + (z-Z 0 ) 2 ⁇ 1/2 .
  • This calculation, the point (x., Y., Z 0) is performed for all surfaces in the vicinity innumerable, and smallest with by point and Ime distance di things therein.
  • FIG. 27 is a conceptual diagram of a conventional comparison method using a three-dimensional image.
  • the conventional comparison method if the surface and the point cloud are compared with high accuracy, the number of target point clouds is innumerable.
  • the actual CAD model has many small surfaces, and the number of surfaces to be compared is enormous. In other words, in the actual CAD model, there are innumerable combinations of points and surfaces to be calculated. As a result, the following problems appear.
  • the number of measurement points required in the actual manufacturing industry is as shown in Table (1) below. Regardless of the method used, the number of measurement points becomes a bottleneck, making it practically impossible to implement. I have.
  • the number of measurement points is assumed to be five times the surface size (400 x 400 mm) of the above parts. Disclosure of the invention
  • the present invention has been made in view of the above points, and provides a three-dimensional image comparison slogram, a three-dimensional image comparison method, and a three-dimensional image that can perform three-dimensional image comparison with high accuracy and high speed. It is an object to provide a comparison device.
  • a three-dimensional image comparison program as shown in FIG. 1 is provided to solve the above problem.
  • the three-dimensional image comparison program of the present invention is applied when comparing three-dimensional images and displaying differences. Further, the three-dimensional image comparison program of the present invention causes a computer to execute the following processing.
  • the first three-dimensional image previously converted into a poxel is compared with the second three-dimensional image to generate a three-dimensional difference image (step S 1).
  • a fine 3D difference image between the first 3D image and the second 3D image is generated with higher precision than the comparison on the surface of the 3D difference image (step S). 2).
  • comparing the three-dimensional difference minute image of the first three-dimensional image with the three-dimensional difference minute image of the second three-dimensional image, the difference part of the three-dimensional difference image The display format is determined (step S3).
  • the three-dimensional difference image is displayed in the display format as a result of the poxel comparison (step S).
  • the computer first compares the first three-dimensional image and the second three-dimensional image, which have been converted into the poxels, in advance.
  • a three-dimensional difference image is generated.
  • a three-dimensional difference image between the first three-dimensional image and the second three-dimensional image is generated on the surface of the three-dimensional difference image with finer precision than the comparison. Is done.
  • a display format indicating the difference of the poxels is determined for the three-dimensional difference image and the three-dimensional difference minute image. Then, the three-dimensional difference image is displayed in the display format as a result of the poxel comparison.
  • a computer compares a first three-dimensional image previously converted into a poxel with a second three-dimensional image, and generates a three-dimensional difference image. Generate and partially extract the difference part of the three-dimensional difference image, count the number of non-matching poxels in the vertical direction from each face of the three-dimensional shape with respect to the matching surface poxels, and count Determining a display format representing the difference of the pixels on the basis of the number of pixels, and displaying the three-dimensional difference image having the determined display format in the display format.
  • a 3D image comparison program is provided.
  • the computer first compares the first three-dimensional image and the second three-dimensional image, which have been previously converted into poxels, to obtain a three-dimensional difference image. Is generated. Further, the difference part of the three-dimensional difference image is partially extracted, and the number of non-matching poxels is counted in the vertical direction from each face of the three-dimensional shape with respect to the matching surface poxels. Further, a display format indicating the difference of the poxels is determined based on the number of the force points. Then, the three-dimensional difference image for which the display format has been determined is displayed in the display format.
  • a computer compares a first three-dimensional image and a second three-dimensional image that have been previously converted into a poxel to generate a three-dimensional difference image. Generate the ratio of the number of surface poxels in the difference portion of the three-dimensional difference image to the total number of poxels, and display a display format representing the difference of the poxels based on the ratio with respect to the three-dimensional difference image. And determine the display format A three-dimensional image comparison program is provided, wherein the determined three-dimensional difference image is displayed in a predetermined display format corresponding to the ratio to execute the processing.
  • a computer first compares a first three-dimensional image and a second three-dimensional image that have been previously converted into a poxel to generate a three-dimensional difference image. You. Next, the ratio between the number of surface poxels in the different part of the three-dimensional difference image and the total number of poxels is calculated, and a display format representing the difference of the poxels is determined for the three-dimensional difference image based on the ratio. Is done. Then, for the three-dimensional difference image for which the display format has been determined, the difference portion is displayed in a predetermined display format corresponding to the ratio.
  • FIG. 1 is a principle configuration diagram in the first embodiment.
  • FIG. 2 is a conceptual diagram of a method for comparing three-dimensional images according to the first embodiment.
  • FIG. 3 is a diagram illustrating an example of a hardware configuration of a computer capable of executing a three-dimensional image comparison program.
  • FIG. 4 is a functional block diagram showing a functional configuration of the three-dimensional image comparison device.
  • FIG. 5 is a diagram showing an outline of the comparison processing of the three-dimensional image comparison program.
  • FIG. 6 is a diagram showing a procedure for comparing the compared three-dimensional images with finer poxels.
  • FIG. 7 is a diagram showing color determination for a common poxel in a method of setting a color determined by the number of poxel differences.
  • FIG. 8 is a diagram showing color determination for a difference poxel (reference image) in a method of setting a color determined by the number of poxel differences.
  • FIG. 9 is a diagram showing color determination for a difference poxel (comparison target image) in a method of setting a color determined by the number of poxel differences.
  • FIG. 10 is a diagram showing the relationship between element numbers of nodes and nodes.
  • FIG. 11 is a diagram showing a general data structure at a node on a three-dimensional image.
  • FIG. 12 is a diagram showing a data structure of the poxel method at a poxel node.
  • FIG. 13 is a flowchart showing the overall flow of the three-dimensional image comparison processing.
  • FIG. 14 is a flowchart showing comparison information input processing in the three-dimensional image comparison processing of FIG.
  • FIG. 15 is a flowchart showing the poxel processing in the three-dimensional image comparison processing of FIG.
  • FIG. 16 is a flowchart showing the image registration processing in the three-dimensional image comparison processing of FIG.
  • FIG. 17 is a flowchart showing a first image comparison process in the three-dimensional image comparison process of FIG.
  • FIG. 18 is a diagram showing an outline of the comparison processing of the three-dimensional image comparison program in the second embodiment.
  • FIG. 19 is a flowchart showing a second image comparison process in the three-dimensional image comparison process of FIG.
  • FIG. 20 is a diagram showing a difference evaluation (1) of the comparison processing of the three-dimensional image comparison program in the third embodiment.
  • FIG. 21 is a diagram showing colors (light colors) determined by the comparison processing of the three-dimensional image comparison program in the third embodiment.
  • FIG. 22 is a diagram showing a difference evaluation (2) of the comparison processing of the three-dimensional image comparison program in the third embodiment.
  • FIG. 23 is a diagram showing colors (dark colors) determined by the comparison processing of the three-dimensional image comparison program in the third embodiment.
  • FIG. 24 is a diagram showing the difference evaluation (3) of the comparison processing of the three-dimensional image comparison program in the third embodiment.
  • FIG. 25 is a flowchart showing a third image comparison process in the three-dimensional image comparison process of FIG.
  • Figure 26 is a diagram showing a three-dimensional image (CAD data: three-dimensional curved surface) of an object and a three-dimensional point cloud (measurement data).
  • FIG. 27 is a conceptual diagram of a conventional comparison method using a three-dimensional image.
  • the present invention can be applied to, for example, a method for comparing geometric shapes with the following three-dimensional images.
  • two 3D images are both expressed as a uniform cube or a rectangular parallelepiped (that is, a voxel), a Boolean operation is performed for each individual poxel, and the two images are compared. This is done by evaluating the agreement and disagreement.
  • the two images are aligned based on reference points or characteristic shapes.
  • image differences are represented by the number of unmatched poxels.
  • the information of the unmatched pixels is displayed on the surface.
  • the number of difference pixels in the depth direction from each surface (X method, y direction, z direction) of the reference poxel is displayed in a predetermined display format for each surface of the reference poxel.
  • the difference in the three-dimensional image is represented by each surface of the reference poxel (color corresponding to the depth of the number of difference poxels).
  • a three-dimensional image comparison program capable of expressing a difference in a three-dimensional image by a point (color of a poxel), a three-dimensional image Describes the comparison method and 3D image comparison device I will tell.
  • a three-dimensional image comparison program capable of expressing a difference in a three-dimensional image on each surface of a reference poxel (a color corresponding to the depth of the number of difference poxels), The comparison method and the three-dimensional image comparison device are described.
  • a three-dimensional image comparison program capable of expressing a difference of a three-dimensional image by the entire difference part (color corresponding to the ratio) Will be described.
  • three-dimensional images are compared in units of poxels, and differences in three-dimensional images are displayed in a predetermined display format.
  • FIG. 1 is a principle configuration diagram in the first embodiment.
  • the three-dimensional image comparison program according to the first embodiment causes a computer to perform comparison Z difference display of three-dimensional images on a computer.
  • the comparison Z difference display by computer is executed on the following page of Tagawa. This shows a procedure in which a computer compares a 3D image B1 in which a design model is converted into a poxel with a 3D image B2 in which measured data is converted into a poxel, and displays the data in a predetermined display format.
  • the computer generates a three-dimensional difference image VB1 by comparing the three-dimensional image B2 and the three-dimensional image B2 that have been converted into the poxels in advance (step S1).
  • the computer generates, on the surface of the three-dimensional difference image VB1, a three-dimensional difference minute image VB2 of the three-dimensional image B1 and the three-dimensional image B2 with finer accuracy than the comparison in step S1. (Step S2).
  • the computer determines the display format of the difference portion of the three-dimensional difference image VB1 by comparing the three-dimensional difference detail image of the three-dimensional image B1 and the three-dimensional difference detail image of the three-dimensional image B2. Yes (step S3).
  • This display format is an additive color mixture (RGB) based on the three primary colors of red, green, and blue.
  • the display format is subtractive color mixing (C) based on the three primary colors of cyan, magenta, and yellow. MY). If subtractive color mixing is used, black can be added in addition to the three primary colors of cyan, magenta, and yellow.
  • the computer displays the three-dimensional difference image V B1 whose display format has been determined on a display or the like in a predetermined display format as a voxel comparison result (step S4).
  • the poxel comparison result can be printed not only on a display or the like but also on a printer or a plotter.
  • the computer compares the three-dimensional image B1 previously converted into the poxels with the three-dimensional image B2 to generate the three-dimensional difference image VB1. .
  • the computer generates a detailed three-dimensional difference image VB2 of the three-dimensional image B1 and the three-dimensional image B2 on the surface of the three-dimensional difference image VB1 with a finer precision than the comparison of the step S1. Is done.
  • the computer compares the three-dimensional difference minute image of the three-dimensional image B1 with the three-dimensional difference minute image of the three-dimensional image B2 to determine the display format of the difference part of the three-dimensional difference image VB1. Is done.
  • the computer displays the three-dimensional difference image V B 1 whose display format has been determined on a display or the like in a predetermined display format as a result of the poxel comparison. This makes it possible to perform the comparison processing with high accuracy and high speed in comparison of three-dimensional images.
  • FIG. 2 is a conceptual diagram of a method for comparing three-dimensional images according to the first embodiment.
  • the comparison method in the first embodiment as shown in FIG. 2, it is not possible to directly compare the surface data I m1 1 which is a design model with another surface data I m1 2. Therefore, in the comparison method according to the first embodiment, the volumes enclosed by the surface data I ml 1 and I ml 2 are converted into poxels, and the poxels B 1 1 and B 12 are generated, respectively. The comparison between 1 and the poxel B 1 2 is performed. This makes it possible to easily compare two 3D images.
  • FIG. 3 is a diagram illustrating an example of a hardware configuration of a computer that can execute a three-dimensional image comparison program.
  • the entire computer 100 is controlled by a CPU (Central Processing Unit) 101.
  • the CPU 101 is connected to a random access memory (RAM) 102, a hard disk drive (HDD) 103, a graphic processing device 104, an input interface 105, and a communication interface 106 via a bus 107. .
  • RAM random access memory
  • HDD hard disk drive
  • the RAM 102 temporarily stores at least a part of an S (Operating System) program and an application program to be executed by the CPU 101. Further, the RAMI 02 stores various data required for processing by the CPU 101.
  • the HDD 103 stores an OS and application programs such as a three-dimensional image comparison program.
  • a monitor P111 is connected to the graphic processing device 104. Further, the graphic processing device 104 displays an image on the screen of the monitor P111 according to a command from the CPU 101.
  • the input interface 105 is connected with a keypad P112 and a mouse P113. The input interface 105 transmits a signal transmitted from the keyboard P 112 or the mouse P 113 to the CPU 101 via the bus 107.
  • Communication interface 106 is connected to network 110.
  • the communication interface 106 transmits and receives data to and from another computer via the network 110.
  • the computer 100 By causing the computer 100 having the above hardware configuration to execute the three-dimensional image comparison program, the computer 100 functions as a three-dimensional image comparison device.
  • the processing functions of the three-dimensional image comparison device realized by the computer executing the three-dimensional image comparison program will be described below.
  • FIG. 4 is a functional block diagram showing a functional configuration of the three-dimensional image comparison device.
  • the three-dimensional image comparison device includes a comparison information input unit 10 for inputting comparison information, a poxel processing unit 20 for generating a poxel from the input three-dimensional image, and a Image alignment processing unit 30 for aligning dimensional images, Poxel difference evaluation processing for evaluating and displaying poxel differences And a poxel superimposition processing unit 50 for superimposing the displayed poxels.
  • the comparison information input unit 10 apart from its own function, is a measurement accuracy input unit 11, a measurement accuracy storage unit 12, a three-dimensional surface input unit 13, an image storage unit 14, a point cloud It has a processing switching section 15, a point group (surface) input section 16, a point group (surface) storage section 17, a point group (volume) input section 18, and a point group (volume) storage section 19.
  • the poxel processing section 20 has a poxel generating section 21, a poxel generating section 22, and a poxel storage section 23, separately from its own function. The details of these functions are listed below.
  • the comparison information input unit 10 is connected to the poxel processing unit 20, the three-dimensional CAD image DB 11, and the three-dimensional measurement image DB 12. Enter Here, the comparison information input unit 10 accepts a user's instruction to input a three-dimensional image. Further, the comparison information input unit 10 determines the image type of the received three-dimensional image. Then, the comparison information input unit 10 determines whether or not the three-dimensional image is stored in the image storage unit 14, the point cloud (surface) storage unit 17, or the point cloud (volume) storage unit 19. I do.
  • the measurement accuracy input unit 11 is connected to the measurement accuracy storage unit 12 and receives an input of the measurement accuracy by the user.
  • the measurement accuracy can be freely determined by the user, and is the minimum unit of a 3D image to be turned into a poxel.
  • the measurement accuracy it is preferable to use, for example, the required accuracy (0.01 mm or the like) of the product to be manufactured.
  • the measurement accuracy storage unit 12 is connected to the measurement accuracy input unit 11 and stores the measurement accuracy received by the measurement accuracy input unit 11.
  • the three-dimensional surface input unit 13 is connected to the image storage unit 14 and the three-dimensional CAD image DB 11, and if the image type of the three-dimensional image specified by the user is a three-dimensional surface. For example, input of a three-dimensional surface is received from the three-dimensional CAD image DB 11.
  • the image storage unit 14 is connected to the three-dimensional surface input unit 13 and stores the three-dimensional surface received by the three-dimensional surface input unit 13 in an internal storage device.
  • the point cloud processing switching section 15 is connected to the point cloud (surface) input section 16 and the point cloud (volume) input section 18. Connected to switch point cloud processing.
  • the point cloud processing switching unit 15 is a point cloud (surface) input unit 16 and a point cloud (volume). Switches to input section 18.
  • the point cloud (surface) input unit 16 is connected to the point cloud processing switching unit 15, point cloud (surface) storage unit 17, and three-dimensional measurement image DB 12, and accepts input of measurement surface data.
  • the point cloud (surface) input unit 16 receives the input of the measurement surface data from the three-dimensional measurement image DB 12 if the image type of the three-dimensional image specified by the user is the measurement surface data.
  • the point cloud (surface) storage unit 17 is connected to the point cloud (surface) input unit 16 and stores the measured surface data.
  • the point cloud (surface) storage unit 17 stores the measured surface data received by the point cloud (surface) input unit 16 in an internal storage device.
  • the point cloud (volume) input unit 18 is connected to the point cloud processing switching unit 15, the point cloud (volume) storage unit 19, and the three-dimensional measurement image DB 12, and accepts input of measurement cross-section data. Do.
  • the point cloud (volume) input unit 18 accepts input of measured cross-sectional data from the three-dimensional measured image DB 12 if the image type of the three-dimensional image specified by the user is the measured cross-sectional data.
  • the point cloud (volume) storage unit 19 is connected to the point cloud (volume) input unit 18 and stores the measured section data.
  • the point cloud (volume) storage unit 19 stores the measured cross-sectional data received by the point cloud (volume) input unit 18 in an internal storage device.
  • the poxel processing unit 20 is connected to the comparison information input unit 10 and the image registration processing unit 30 and generates an image by a poxel from the three-dimensional image input by the comparison information input unit 10.
  • the poxel processing unit 20 selects the comparison information recorded in the comparison information input unit 10.
  • the poxel processing unit 20 determines the poxel size based on the measurement accuracy recorded in the comparison information input unit 10. Then, the poxel processing unit 20 receives the image type of the three-dimensional image stored in the comparison information input unit 10 and determines the next process.
  • the poxel processing unit 20 selects a three-dimensional surface from the image storage unit 14 of the comparison information input unit 10 if the image type is a three-dimensional surface. Select and read, and pass this 3D surface to the poxel generator 21. If the image type is the measurement surface data, the poxel processing unit 20 selects and reads the measurement surface data from the point cloud (surface) storage unit 17 of the comparison information input unit 10 and reads the measurement surface data. To the poxel generator 22.
  • the poxel processing section 20 selects and reads the measurement section data from the point cloud (volume) storage section 19 of the comparison information input section 10.
  • the measurement section data is passed to the poxel generator 22.
  • the poxel processing section 20 determines whether or not the poxel is generated by the poxel generating section 21 and the poxel generating section 22.
  • the poxel generator 21 generates an image based on the passed three-dimensional surface.
  • the poxel generator 22 generates an image using the poxels based on the passed measurement surface data. Further, the poxel generating unit 22 generates an image using the poxels based on the passed measurement cross-sectional data.
  • the image alignment processing unit 30 is connected to the poxel processing unit 20 and the poxel difference evaluation processing unit 40, and aligns the two images converted into the poxels by the poxel processing unit 20.
  • the image registration processing unit 30 selects a reference point that has been added to the two three-dimensional images in advance.
  • the two three-dimensional images refer to a three-dimensional surface that is CAD data, and measurement surface data or measurement cross-section data that is measurement data. Further, the image registration processing unit 30 calculates an offset from the selected reference point. Further, based on the calculated offset, the image registration processing unit 30 moves the three-dimensional surface converted into the poxels. Then, the image registration processing unit 30 stores the moved three-dimensional surface in an internal storage device.
  • the poxel difference evaluation processing unit 40 is connected to the image alignment processing unit 30 and the poxel overlay processing unit 50, and is aligned by the image alignment processing unit 30.
  • the two images are compared and evaluated, and the two evaluated images are displayed in a predetermined display format.
  • the predetermined display format is, for example, that the depth of each color of RGB is changed according to the difference between images and displayed.
  • the difference evaluation processing unit 40 selects a surface poxel. With this selection, the poxel difference evaluation processing unit 40 determines whether the selected poxel is a surface poxel. Then, if a surface poxel is selected by this determination, the poxel difference evaluation processing unit 40 re-poxels with evaluation accuracy. That is, according to the evaluation accuracy, the surface pixels are made into finer pixels.
  • the evaluation accuracy can be determined in advance by the user. In addition, the evaluation accuracy can be refined stepwise from several times the measurement accuracy.
  • the poxel difference evaluation processing unit 40 compares and counts the difference poxels based on the three-dimensional surface converted into the poxels and the measured surface Z-section data. The details of this comparison processing will be described later.
  • the poxel difference evaluation processing section 40 determines the color depth as a display format based on the counted difference poxels, and colors the poxels.
  • the poxel difference evaluation processing unit 40 determines whether or not the poxels are colored. That is, the poxel difference evaluation processing unit 40 determines whether or not the three-dimensional images converted into the poxels have been compared.
  • the poxel overlay processing section 50 is connected to the poxel difference evaluation processing section 40, and overlays a plurality of images displayed by the poxel difference evaluation processing section 40.
  • the three-dimensional CAD image DB 11 is connected to the comparison information input unit 10 and stores a three-dimensional surface that is CAD data.
  • the three-dimensional measurement image DB12 is connected to the comparison information input unit 10, and stores measurement surface data and measurement cross-section data that are measurement data.
  • the calculation time is reduced by comparing the three-dimensional shapes with the poxels.
  • the image is divided into fine and uniform cubes (poxels), the inside of which is turned on, and the outside is turned off.
  • Poxels fine and uniform cubes
  • the image measured by the X-ray CT scanner pixels with high X-ray density are turned on, and those with low X-ray density are turned off.
  • FIG. 5 is a diagram showing an outline of the comparison processing of the three-dimensional image comparison program.
  • the comparison process converts the two images into a poxel like a three-dimensional image B21 and a three-dimensional image B22, and sets the overlapping part (3D difference image VB21) as a Boolean image.
  • the difference can be expressed visually.
  • the computational complexity of the comparison is several tens of 1Z of the conventional method.
  • a non-matching pixel can be extracted by performing an exclusive OR (XOR) between the pixel of the three-dimensional image B21 and the pixel of the three-dimensional image B22.
  • XOR exclusive OR
  • FIG. 6 is a diagram showing a procedure for comparing the compared three-dimensional images with finer poxels.
  • a coarse poxel (3D difference image VB21) is generated.
  • the poxel including the surface is finely divided (three-dimensional difference fine image VB22), and the difference between the two images is evaluated based on the number of difference poxels.
  • the difference between the two images is "8". This indicates, for example, the difference between the fine poxels in the three-dimensional difference fine image VB22 of the images I21 and I22 in the coarse poxels.
  • the color depth of the poxels is determined based on the number of differences.
  • the overall difference becomes visually understandable.
  • the image I22 can be determined to be a difference.
  • the poxel Bx16 since the color of this portion (for example, the depth of red) is determined, it can be seen that the image I21 is a difference.
  • the color of this part since the color of this part is not determined, it can be determined that there is no difference.
  • the poxels are further subdivided by the poxels having a similarity ratio of 1/4.
  • the difference is up to 64 (4X4X4) poxels.
  • FIG. 7 is a diagram showing color determination for a common poxel in a method of setting a color determined by the number of poxel differences.
  • the brightness is set to the maximum value (for example, 255, 255, 255) for the three RGB colors at the setting position ST1. I do.
  • FIG. 8 is a diagram showing color determination for a difference poxel (reference image) in a method of setting a color determined by the number of poxel differences.
  • FIG. 9 is a diagram showing color determination for a difference poxel (comparison target image) in a method of setting a color determined by the number of poxel differences.
  • the color values are determined in proportion to the number of poxels in each image, as shown in Figs.
  • FIG. 10 is a diagram showing the relationship between the element numbers of nodes and nodes.
  • the poxels converted from the three-dimensional image are, for example, eight poxels BX1 to BX8.
  • the poxel B xl is composed of articulation points (1, 2, 5, 4, 10, 11, 14, 13).
  • the I injection number at the node is counted counterclockwise on each of the upper and lower surfaces.
  • the coordinate information is registered in the poxel data structure for each of such nodes. Have been. Therefore, the details of this coordinate information will be described with reference to FIG. 11 and FIG.
  • FIG. 11 is a diagram showing a general data structure at a node of a three-dimensional image.
  • FIG. 11 shows a data structure in a typical finite element analysis.
  • the data structure D1 of the poxels has a data type area D11 indicating a data type, a data number area D12 indicating a serial number in each data type, and a data type area D12 indicating a serial number in each data type. It comprises a position information area D13 indicating position information.
  • the overnight area D11 is composed of two items, "GRID” and "CHEXA". Where “GR I D” indicates a node and "CHE
  • XA indicates the elements obtained by dividing the 3D image.
  • RID " is assigned a serial number of 1 to 27.
  • the serial numbers 1 to 27 are the node numbers described above with reference to FIG. 10.
  • the data number area D12 for example, As shown in Fig. 11, for the data type "CHEXA",
  • serial numbers 1 to 8 are the element numbers of the poxels described above with reference to FIG.
  • the coordinate information at each node of “0.000000, 0.000000, 50.00000” is set for the data number 1 of “ID.” This is the same as the data number (“GR ID”).
  • the coordinates of x, y, and z are represented by floating-point numbers, which are set as shown in the figure below. For example, as shown in Fig. 11, for the data number 1 of the data type "CHEXA", the node ID of each element of "17, 21, 27, 22, 7, 10, 26, 1" is set. This is the data number ("CHEXA"
  • FIG. 12 is a diagram showing the data structure of the poxel method at the nodes of the poxels.
  • Fig. 12 shows the data structure in the case of expressing three dimensions by the Poxel method. ing.
  • the poxel data structure D2 sets the size of the element and the node of each element for each of the predefined records.
  • the size of the element the length of the side of X, y, z is represented by a floating point. Nodes that define each element are set as nodes of each element. However, the element IDs described above do not need to be defined because they are regularly arranged.
  • the record (RECORD 7) has an element size expressed in exponential notation (floating point notation) “2.500000 e + 001, 2.50000000 e + 001, 2. 50000004 e + 001 "is set.
  • FIG. 13 is a flowchart showing the overall flow of the three-dimensional image comparison processing. This process is a process in which the three-dimensional image comparison device is activated by turning on the power, activating a program, or any predetermined timing, and is executed by the CPU 101. Hereinafter, the processing illustrated in FIG. 13 will be described along the step numbers. The names of the functions in this chart are explained based on Fig. 4.
  • the comparison information input unit 10 performs a comparison information input process of inputting CAD data of a three-dimensional image and measurement data, which are comparison information. The details of the comparison information input processing will be described later with reference to FIG.
  • Step S20 The poxel processing unit 20 performs a poxel process of generating an image using the poxels from the three-dimensional image input in step S10. The details of the poxel processing will be described later with reference to FIG.
  • Step S30 The image alignment processing unit 30 determines in step S20 Performs image alignment processing to align the two images converted to poxels. The details of the image registration process will be described later with reference to FIG.
  • Step S40 The poxel difference evaluation processing unit 40 performs an image comparison process of comparing and evaluating the two images positioned in step S30. The details of the image comparison process will be described later with reference to FIG.
  • the poxel difference evaluation processing unit 40 displays the plurality of images compared and evaluated in step S40 in a predetermined display format.
  • the predetermined display format is, for example, that the depth of each color of RGB is changed according to the difference between images and displayed.
  • the poxel overlay processing unit 50 overlays the plurality of images displayed in step S50.
  • the comparison result is superimposed on the 3D image.
  • a three-dimensional difference image which is a comparison result, is generated by exclusive OR based on the two images that have been subjected to the poxelation, and only the difference part can be visualized.
  • the poxel superimposition processing unit 50 superimposes the two poxelized images, the three-dimensional difference image, or the three-dimensional difference fine image, and displays the original image and the difference part in different colors. Can be done.
  • FIG. 14 is a flowchart showing comparison information input processing in the three-dimensional image comparison processing of FIG.
  • step S10 comparison information input processing
  • processing is performed according to the following flow. This comparison information input processing is performed by the comparison information input unit 10 of the three-dimensional image comparison device.
  • the measurement accuracy input unit 11 of the comparison information input unit 10 accepts the input of the measurement accuracy by the user.
  • Step S102 The measurement accuracy storage unit 12 of the comparison information input unit 10 stores the measurement accuracy received in step S101.
  • the comparison information input unit 10 receives a user's input instruction of a three-dimensional image.
  • Step S104 The comparison information input unit 10 accepts in step S103 Determines the image type of the 3D image.
  • the comparison information input unit 10 proceeds to step S105 if the image type is the three-dimensional surface, and proceeds to step S107 if the image type is the measurement surface data, or If the image type is the measurement section data, the flow advances to step S109. If the image type is measurement surface data or measurement section data, the point group (surface) input section 16 and the point group (surface) input section 16 are switched by the point group processing switching section 15 of the comparison information input section 10. (Volume) Switching to input section 18 is performed.
  • Step S105 The 3D surface input unit 13 of the comparison information input unit 10 uses the 3D CAD image DB 11 since the image type is 3D surface in step S104. Accept input of 3D surface.
  • Step S106 The image storage unit 14 of the comparison information input unit 10 stores the three-dimensional surface received in step S105 in an internal storage device.
  • Step S107 The point group (surface) of the comparison information input unit 10 (surface)
  • the input unit 16 is a 3D measurement image because the image type is the measurement surface data in step S104. Accept input of measurement surface data from DB12.
  • Step S108 The point group (surface) storage unit 17 of the comparison information input unit 10 stores the measurement surface data received in step S107 in an internal storage device.
  • the point group (volume) input part 18 of the comparison information input part 10 is a three-dimensional measurement image because the image type is the measurement section data in step S104. Accepts input of measurement section data from DB12.
  • Step S110 The point cloud (volume) storage unit 19 of the comparison information input unit 10 stores the measured cross-sectional data received in step S109 in an internal storage device.
  • Step S111 The comparison information input unit 10 determines whether a three-dimensional image has been stored in step S106, step S108, or step S110. Here, the comparison information input unit 10 returns to step S104 if a three-dimensional image has not been recorded, and repeats the same processing. If a three-dimensional image has been recorded, this processing ends. Then, the process returns to step S10 in FIG.
  • FIG. 15 is a flowchart showing the poxel processing in the three-dimensional image comparison processing of FIG.
  • step S20 poxel processing
  • Fig. 13 Processing is performed according to the following flow. This poxel processing is performed by the poxel processing unit 20 of the three-dimensional image comparison device.
  • Step S201 The poxel processing unit 20 selects the comparison information stored in step S10.
  • Step S202 The poxel processing unit 20 determines the poxel size based on the measurement accuracy stored in step S10.
  • Step S203 The poxel processing unit 20 receives the type of the three-dimensional image stored in step S10.
  • Step S204 The poxel processing unit 20 determines the image type of the three-dimensional image received in step S203.
  • the poxel processing unit 20 proceeds to step S205 if the image type is a three-dimensional surface, and proceeds to step S207 if the image type is a measurement surface data. If the image type is the measurement section data, the process proceeds to step S209.
  • Step S 205 Since the image type is a three-dimensional surface in step S 204, the poxel processing unit 20 changes the image storage units 14 to 3 of the comparison information input unit 10. Select and read the dimensional surface, and pass this 3D surface to the poxel generator 21.
  • Step S206 The poxel generation unit 21 of the poxel processing unit 20 generates a 3D image by the poxel based on the 3D surface passed in step S205. Generate.
  • the poxel processing unit 20 stores the point group (surface) of the comparison information input unit 10 since the image type is measured surface data in step S204. Select and read measurement surface data from 17 and pass this measurement surface data to the pixel generator 22.
  • Step S208 The poxel generating unit 22 of the poxel processing unit 20 generates a three-dimensional image using the poxels based on the measurement surface data passed in step S207. .
  • the poxel processing unit 20 stores the point group (volume) storage unit of the comparison information input unit 10 because the image type is the measurement section data in step S204. Select and read the measurement cross-section data from 19 and pass this measurement cross-section data to the pixel generator 22.
  • Step S210 The poxel generation unit 22 of the poxel processing unit 20 generates a three-dimensional image using the poxels based on the measurement cross-sectional data passed in step S209. .
  • step S211 the poxel processing section 20 determines whether or not a three-dimensional image is generated by the poxels. Here, if no poxel has been generated, the process returns to step S204 to repeat the same processing. If a poxel has been generated, the poxel processing unit 20 ends this process and returns to the step in FIG. Return to S20.
  • FIG. 16 is a flowchart showing an image registration process in the three-dimensional image comparison process of FIG.
  • step S30 image registration processing
  • processing is performed according to the following flow. Note that this image registration processing is performed by the image registration processing unit 30 of the three-dimensional image comparison device.
  • the image alignment processing unit 30 selects a reference point added to two three-dimensional images in advance.
  • the two 3D images refer to a 3D surface that is CAD data and a surface data or a cross section that is a measurement data.
  • Step S302 The image registration processing unit 30 calculates an offset from the reference point selected in step S301.
  • Step S303 The image alignment processing unit 30 moves the poxel-transformed three-dimensional image based on the offset calculated in step S302.
  • Step S304 The image registration processing unit 30 stores the three-dimensional image moved in step S303 in an internal storage device, ends this processing, and ends the processing in step S3 of FIG. Returns to 0.
  • FIG. 17 is a flowchart showing a first image comparison process in the three-dimensional image comparison process of FIG.
  • step S40 image comparison processing
  • processing is performed according to the following flow. Note that this first image comparison processing is performed by the poxel difference evaluation processing unit 40 of the three-dimensional image comparison device. Also,
  • Poxel difference evaluation processing section 40 selects a surface poxel.
  • the selection of the surface poxels can be easily determined by, for example, that the inside of the surface poxels is on and the outside is off for a three-dimensional image that has been turned into a poxel.
  • Step S 4 1 2 The poxel difference evaluation processing section 40 determines whether the poxel selected in step S 411 is a surface poxel. Here, the poxel difference evaluation processing unit 40 proceeds to step S413 if a surface poxel is selected, and returns to step S411 if a surface poxel is not selected.
  • Step S4 13 The poxel difference evaluation processing section 40 performs re-poxelization with evaluation accuracy because the surface poxel is selected in step S4 12. That is, the surface poxels are further refined according to the evaluation accuracy.
  • the poxel difference evaluation processing unit 40 compares and counts the difference poxels based on the three-dimensional surface converted into the poxels and the measured surface Z section data.
  • Step S4l5 The poxel difference evaluation processing section 40 determines the color depth of the display format based on the difference poxels counted in step S4l4, and colors the poxels. Do.
  • Step S 4 16 The poxel difference evaluation processing section 40 determines whether or not the poxels have been colored in step S 4 15. That is, the poxel difference evaluation processing unit 40 determines whether or not the three-dimensional images converted into the poxels have been compared. Here, if the poxel is not colored, the poxel difference evaluation processing unit 40 returns to step S411 and repeats the same processing. If the poxel is colored, This process ends, and the process returns to step S40 in FIG.
  • the computer compares the three-dimensional image B 21 previously converted into the poxels with the three-dimensional image B 22 to obtain the three-dimensional difference image VB 21. Generated.
  • the computer converts the three-dimensional difference image VB 22 between the three-dimensional image B 21 and the three-dimensional image B 22 to the surface of the three-dimensional difference image VB 21 with higher precision than the previous comparison. Generated.
  • the computer compares the three-dimensional difference detailed image of the three-dimensional image B21 with the three-dimensional difference fine image of the three-dimensional image B22, and displays the difference part of the three-dimensional difference image VB21.
  • the format is determined. In this determination, the color is determined according to the number of differences of the poxels.
  • the computer displays the three-dimensional difference image VB21 whose display format has been determined on a display or the like in a predetermined display format as a result of the poxel comparison.
  • This makes it possible to perform the comparison processing with high accuracy and high speed in comparison of three-dimensional images.
  • the image comparison target a three-dimensional poxel, it is not necessary to compare innumerable point clouds, and it is possible to increase the efficiency of memory resources and reduce the number of operations.
  • the number of operations can be significantly reduced.
  • the 3D image surface coarse pixels are compared to generate a 3D difference fine image, and the comparison is performed. Only the coarse poxels with differences can be compared by generating a 3D difference fine image.
  • the second embodiment differs from the first embodiment in the form of comparison processing. That is, in the second embodiment, three-dimensional images are compared on a per-poxel basis, and the difference between the three-dimensional images is further determined by a predetermined display format corresponding to the depth of the number of different poxels from each surface of the reference poxel. And each surface of the reference poxel is displayed.
  • FIG. 18 is a diagram illustrating an outline of a comparison process of the three-dimensional image comparison program in the second embodiment.
  • the comparison process in the second embodiment uses the poxels created in the first embodiment, and displays the difference in color for each face of the poxels in one of the three-dimensional images.
  • the color is determined by the number of poxels in the vertical direction (depth) from each surface as follows.
  • a three-dimensional image is compared using the poxels created in the first embodiment (step Sib).
  • the mismatched part of the compared three-dimensional difference image VB21 is partially extracted (step S2b).
  • the number of unmatched poxels (vertical direction) from each face (zy plane, zx plane, yx plane) of the three-dimensional shape is determined (step S3b).
  • the computer compares the three-dimensional images using the poxels created in the first embodiment. Also, the mismatched part of the compared three-dimensional difference image VB21 is partially extracted.
  • the computer calculates the number of non-matching poxels (vertical poxels BB 1) in the vertical direction from each of the three-dimensional shapes (zy, zx, and yx) with respect to the matching poxels (reference poxel BB1). Is counted, and the display format is determined.
  • the computer displays each of the matching poxels (reference poxels) in a predetermined display format corresponding to the count number.
  • the functions in the above outline differ only in the function of the poxel difference evaluation processing unit 40 in the function block diagram of FIG. Therefore, in the second embodiment, only new functions (including duplication) of the pixel difference evaluation processing unit 40 will be described.
  • the poxel on the surface adjacent to and coincident with the mismatched portion is described as the reference poxel, but in the following description, the surface adjacent to and coincident with the difference portion is described.
  • a poxel is selected and a range is determined, and the poxel when this range is re-poxelled with evaluation accuracy is set as a reference poxel.
  • a new function of the poxel difference evaluation processing unit 40 is that the image alignment processing unit 30 performs an image comparison process of comparing the aligned two images, and compares the compared two images in a predetermined display format. To display.
  • the predetermined display format is, for example, that the depth of each color of RGB is changed and displayed according to a difference in image.
  • the poxel difference evaluation processing unit 40 selects a surface poxel. With this selection, the poxel difference evaluation processing unit 40 determines whether the selected poxel is a surface poxel. Then, if a surface poxel is selected by this determination, the poxel difference evaluation processing unit 40 determines the range of re-poxelization. That is, the part where the image does not match is determined as the range to be investigated in detail.
  • the poxel difference evaluation processing unit 40 re-poxels the determined range with the evaluation accuracy. That is, the surface poxels for which the range has been determined are further refined according to the evaluation accuracy. Then, the poxel difference evaluation processing unit 40 counts the difference poxels based on the three-dimensional surface converted into the fine poxels and the measurement surface / cross-sectional data. In other words, this corresponds to the number of non-matching poxels (the number of difference poxels) from each surface (zy, zx, and yx) of the three-dimensional shape of the reference poxel with respect to the fine poxels of the matching surface (reference poxels). Count).
  • the poxel difference evaluation processing section 40 determines the color depth as a display format based on the counted difference poxels, and colors the poxels. That is, the poxel difference evaluation processing unit 40 determines the display format (color) of each face of the reference poxel with respect to the fine poxels (reference poxels) of the matching surfaces.
  • the poxel difference evaluation processing unit 40 determines whether or not the poxels are colored. That is, the poxel difference evaluation processing unit 40 determines whether or not the three-dimensional images converted into the poxels have been compared.
  • FIG. 19 is a flowchart showing a second image comparison process in the three-dimensional image comparison process of FIG.
  • step S40 image comparison processing
  • Fig. 13 processing is performed according to the following flow.
  • this second image comparison processing is performed by the poxel difference evaluation processing unit 40 of the three-dimensional image comparison device.
  • the poxel on the surface adjacent to and coincident with the non-coincidence (difference) is described as the reference poxel, but in the following description, the surface adjacent to and coincident with the difference A poxel is selected and a range is determined, and the poxel when this range is re-poxelized with the evaluation accuracy is set as a reference poxel.
  • Poxel difference evaluation processing section 40 selects a surface poxel.
  • the selection of the surface poxels can be easily determined by, for example, that the inside of the surface poxels is on and the outside is off for a three-dimensional image that has been turned into a poxel.
  • Step S422 The poxel difference evaluation processing unit 40 determines whether the poxel selected in step S421 is a surface poxel. Here, the poxel difference evaluation processing unit 40 proceeds to step S 423 if a surface poxel is selected, and returns to step S 421 if a surface poxel is not selected.
  • Step S 4 2 3 Since the surface poxels have been selected in step S 4 22, the poxel difference evaluation processing section 40 determines the range of re-poxelization. That is, the part where the images do not match is determined as the range to be investigated in detail.
  • Step S424 The poxel difference evaluation processing unit 40 performs re-poxelization with evaluation accuracy on the range determined in step S424. In other words, the surface poxels for which the range has been determined are further refined into poxels according to the evaluation accuracy.
  • the poxel difference evaluation processing unit 40 counts the difference poxels again based on the three-dimensional surface converted into the poxels and the measured surface Z-section data. That is, this corresponds to the matching surface poxel (reference poxel). On the other hand, the number of non-matching poxels (depth of the number of difference poxels) is counted from each surface (zy plane, zx plane, yx plane) of the three-dimensional shape of the reference poxel.
  • the poxel difference evaluation processing unit 40 determines the color depth, which is the display format, based on the difference poxels counted in step S 4 25, and colorizes the poxels. Do. That is, the poxel difference evaluation processing unit 40 determines the display format (color) of each face of the poxel with respect to the matching face poxel (reference poxel).
  • Step S 4 27 The poxel difference evaluation processing section 40 determines whether or not the poxels are colored in step S 4 26. That is, the poxel difference evaluation processing unit 40 determines whether or not the three-dimensional images converted into the poxels have been compared. Here, the poxel difference evaluation processing unit 40 returns to step S 4 21 and repeats the same processing if the poxel is not colored. To step S40.
  • the comparison processing in the second embodiment is performed by replacing the pixels including the depth direction with three in order to reduce the computer resources (memory amount, disk area, and calculation time) for comparison.
  • the difference is displayed by the color of each surface of the dimensional shape. This reduces computational resources by a fraction.
  • the color that expresses the difference is determined by the depth proportional to the number of depth poxels. Furthermore, this density is averaged based on the maximum value of the total density.
  • the poxels on the surface that are adjacent and coincide with the non-coincidence part are described as the reference poxels, and the surface poxels that are adjacent and coincide with the difference part are selected and the range is determined.
  • the reference poxels the poxels obtained when this range is re-poxelized with the evaluation accuracy are used as the reference poxels, the poxels on the matching surface can be simply used as the reference poxels.
  • the third embodiment differs from the first embodiment in the form of comparison processing. That is, in the third embodiment, a three-dimensional image is compared on a per-poxel basis, and a predetermined display format corresponding to the ratio of the number of surface poxels and the volume in the entire difference part of the three-dimensional image is obtained. Allows This is to display the cell.
  • FIG. 20 is a diagram showing a difference evaluation (1) of the comparison processing of the three-dimensional image comparison program in the third embodiment.
  • the difference between the pixels generated in the first embodiment is only the part where the two images are different.
  • This different part is composed of one or more poxels.
  • the ratio between the number of surface poxels in these islands and the total number of poxels in the island is displayed as a color.
  • the color is determined as follows. In the following description, this island is simply referred to as a difference part.
  • a three-dimensional image A and an image B are compared using the poxels created in the first embodiment.
  • the shape of the difference part is expressed as follows. 'Volume: 1 2 Poxels
  • the ratio of the number of surface poxels Z to the total number of poxels was 2 5 12, indicating that the difference of the poxels for a certain difference portion was small.
  • the volume per surface area is set in colors as shown in Fig. 21.
  • FIG. 21 is a diagram showing colors (light colors) determined by the comparison processing of the three-dimensional image comparison program in the third embodiment.
  • the colors for displaying the different portions are described using RGB as an example.
  • the color of blue (B) is made variable to express the poxel difference of the entire difference part.
  • the RGB settings at the setting position ST11 are: 255 (R), 255 (G), 22 3 (B) is set.
  • This 2 23 (B) shows a pale color, which indicates that the poxel difference of a certain difference portion is light. If this is, for example, 96 (B), then the poxel differences across a given difference can be shown to be strong.
  • FIG. 22 shows the evaluation of the poxel difference when the poxel difference is deep.
  • FIG. 22 is a diagram showing a difference evaluation (2) of the comparison processing of the three-dimensional image comparison program in the third embodiment. Note that FIG. 22 describes a case where the pixel difference is deeper than that of FIG.
  • a three-dimensional image A and an image B are compared using the poxels created in the first embodiment.
  • the shape of the difference part is expressed as follows. • Volume: 4 3 Poxels
  • the ratio of the number of surface poxels Z to the total number of poxels was 3 2 Z 43, indicating that the difference of the poxels for a certain difference portion was deep.
  • the volume per surface area is set in colors as shown in Fig. 23.
  • FIG. 23 is a diagram showing colors (dark colors) determined by the comparison processing of the three-dimensional image comparison program in the third embodiment.
  • the colors for displaying the different portions are described using RGB as an example.
  • the color of blue (B) variable by making the color of blue (B) variable, the poxel difference of the entire difference part is expressed.
  • the RGB settings at the setting position ST12 are 255 (R), 255 (G), 96 (B) is set. This 96 (B) shows a dark color, indicating that the poxel difference in the entire difference part is deep.
  • FIG. 24 is a diagram showing the difference evaluation (3) of the comparison processing of the three-dimensional image comparison program in the third embodiment.
  • the separation (thick line) line in Fig. 24 indicates the separation position of the difference poxel compared with the original difference poxel by the fine poxel.
  • FIG. 24 shows a case where the poxel difference is evaluated in more detail with respect to FIG. The description of the difference evaluation method is omitted because the same processing as in FIGS. 20 and 22 is performed on more detailed poxels.
  • the functions in the above outline differ only in the function of the poxel difference evaluation processing unit 40 in the functional block diagram of FIG. Therefore, in the third embodiment, only new functions (including duplication) of the pixel difference evaluation processing unit 40 will be described.
  • the surface poxels that are adjacent and coincide with the non-coincidence part are selected, the range is determined, and the poxel obtained when the re-poxelization is performed with the evaluation accuracy is used as the reference poxel.
  • a poxel on the surface that is adjacent to and coincides with the non-coincidence portion may be used as the reference poxel.
  • a matching surface poxel may simply be used as the reference poxel.
  • a new function of the poxel difference evaluation processing unit 40 is to compare and evaluate the two images aligned by the image alignment processing unit 30 and display the evaluated two images in a predetermined display format.
  • the predetermined display format is, for example, that the depth of each color of RGB is changed and displayed according to a difference between images.
  • the poxel difference evaluation processing unit 40 selects a surface poxel. With this selection, the poxel difference evaluation processing unit 40 determines whether the selected poxel is a surface poxel. Then, since the surface poxels have been selected, the poxel difference evaluation processing unit 40 determines the range of the re-poxelization. That is, the part where the images do not match is determined as the range to be investigated in detail.
  • the poxel difference evaluation processing unit 40 re-poxels the determined range with the evaluation accuracy. That is, the surface poxels for which the range has been determined are further refined according to the evaluation accuracy. Then, the poxel difference evaluation processing unit 40 counts the number of difference poxels based on the three-dimensional surface converted into the poxels and the measured surface Z-section data. That is, the ratio between the number of surface poxels and the number of all poxels in the entire difference part of the three-dimensional difference image is obtained.
  • the poxel difference evaluation processing unit 40 determines the color depth, which is the display format, based on the obtained difference poxel ratio, and colors the poxels. In other words, the poxel difference evaluation processing unit 40 determines the display format (color) of the entire difference part composed of unmatched poxels.
  • the poxel difference evaluation processing unit 40 determines whether or not the poxels are colored. In other words, the poxel difference evaluation processing unit 40 executes the three-dimensional Determine if the images have been compared.
  • FIG. 25 is a flowchart showing a third image comparison process in the three-dimensional image comparison process of FIG.
  • step S40 image comparison processing
  • FIG. 25 processing is performed according to the following flow. Note that the third image comparison processing is performed by the poxel difference evaluation processing unit 40 of the three-dimensional image comparison device.
  • Poxel difference evaluation processing section 40 selects a surface poxel.
  • the selection of the surface poxels can be easily determined by, for example, that the inside of the surface poxels is on and the outside is off for a three-dimensional image that is turned into a poxel.
  • Step S 4 32 The poxel difference evaluation processing unit 40 determines whether the poxel selected in step S 431 is a surface poxel. Here, the poxel difference evaluation processing unit 40 proceeds to step S433 if a surface poxel has been selected, and returns to step S431 if a surface poxel has not been selected.
  • Step S 4 3 3 Since the surface poxels have been selected in step S 4 32, the poxel difference evaluation processing section 40 determines the range of re-poxelization. That is, the part where the images do not match is determined as the range to be investigated in detail.
  • Step S4334 The poxel difference evaluation processing unit 40 re-poxels the range determined in step S433 with evaluation accuracy. In other words, the surface poxels for which the range has been determined are further refined into poxels according to the evaluation accuracy.
  • the poxel difference evaluation processing unit 40 counts the number of difference poxels based on the three-dimensional surface converted into the poxels and the measured surface / cross-sectional data. That is, the ratio between the number of surface poxels in the entire difference portion of the three-dimensional difference image and the total number of poxels is obtained.
  • the poxel difference evaluation processing unit 40 determines the color depth of the display format based on the ratio of the difference poxels obtained in step S 4 3 5, and colors the poxels. I do. That is, the poxel difference evaluation processing unit 40 determines the display format (color) of the entire difference unit including the unmatched poxels.
  • Step S437 The poxel difference evaluation processing section 40 determines whether or not the poxels are colored in step S436. That is, the poxel difference evaluation processing unit 40 determines whether or not the three-dimensional images converted into the poxels have been compared. Here, the poxel difference evaluation processing unit 40 returns to step S 4 31 if the poxels are not colored, and repeats the same processing. If the poxels are colored, the processing ends and the processing ends. To step S40.
  • the comparison process in the third embodiment reduces the computer resources (memory capacity, disk area, and calculation time) for comparison, and instead of using the pixels including the depth direction, The difference is displayed in the color of the whole part.
  • the ratio of the number of surface poxels / total number of poxels is a maximum of 2 and a minimum of 0 (more precisely, a very small value close to 0). Therefore, the color depth is set so that 2 is the lightest color and 0 is the darkest color.
  • the size of the different part can be understood from the size of the entire different part when displayed in three dimensions.
  • computational resources are evaluated as follows.
  • this calculation requires calculation of tens of thousands of floating point operations to find one shortest distance. This is repeated for 100 billion points. Ultimately, the time required for this calculation is about 30 hours, assuming that 2 floating point operations can be performed simultaneously at 1 GHz CPU.
  • the method according to the embodiment of the present invention only divides the entire space into poxels with measurement accuracy for a volume including a plurality of surfaces and performs a logical operation for each poxel. 100 million floating point operations. Therefore, the same computer takes 25 seconds to calculate.
  • the method according to the first embodiment of the present invention can greatly reduce the processing time of the computer as compared with the conventional method.
  • the method according to the first embodiment of the present invention divides the entire space into poxels having measurement accuracy, so that comparison and calculation can be performed with the same accuracy as the conventional method.
  • Computer-readable recording media include magnetic recording media, optical discs, magneto-optical recording media, and semiconductor memories. Magnetic recording media include hard disks, flexible disks (FD), ZIP (a type of magnetic disk), and magnetic tape.
  • Optical discs include DVD (Digital Versatile Disc), DVD-RAM (DVD Random Access Memory), CD-ROM (Compact Disc Read Only Memory), CD-R (CD Recordable), and CD-RW (CD Rewritable) ) and so on.
  • Magneto-optical recording media include M ⁇ (Magneto Optical Disk).
  • Semiconductor memories include flash memory. When distributing computer programs, portable recording media such as DVDs and CD-ROMs on which the computer programs are recorded are sold.
  • the computer program can be stored in the storage device of the server, and the program can be transferred from the server to the client via the network.
  • two three-dimensional images converted into poxels are compared, and a surface based on the comparison is generated with further finer accuracy, and the comparison is performed. Since the display format is determined and the display format is used, it is possible to compare the three-dimensional images with high accuracy and high speed.
  • a display format corresponding to a ratio obtained by dividing the surface area of the entire difference part of the three-dimensional image by the volume is determined. Displaying in the display format enables 3D images to be compared with high accuracy and high speed.

Landscapes

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

Description

明 細 書
3次元イメージ比較プログラム、 3次元イメージ比較方法、 および 3次元ィメ一 ジ比較装置 技術分野
本発明は 3次元イメージを比較して差異を表示する 3次元イメージ比較プログ ラム、 3次元イメージ比較方法、 および 3次元イメージ比較装置に関し、 特に 3 次元ィメ一ジのポクセル差異を高精度で比較する 3次元ィメージ比較プログラム、 3次元イメージ比較方法、 および 3次元イメージ比較装置に関する。 背景技術
近年、 物体を 3次元イメージで表現する技術は、 さまざまな分野で利用されて いる。 この 3次元イメージには、 たとえば、 各種製造業での 3次元 C AD (Com puter Aided Design)による設計イメージや、 医療分野での 3次元超音波画像診 断装置による患部立体イメージ画像などがある。 また、 これらの 3次元イメージ を利用することにより、 たとえば、 以下のようなことが可能となる。
- 3次元 C ADによる設計イメージであれば、 デザインレビューや製品検査等 の際に設計担当者や検査担当者は、 製造された製品あるいは試作品と、 設計ィメ 一ジとを容易に確認 (製造検証) することができる。
• 3次元超音波画像診断装置による患部立体イメージ画像であれば、 診断の際 に医師は、 患部の異常や身体部位等の変形をより視覚的に確認 (把握) すること ができる。
このような目的で、 3次元イメージを効率的にかつ精度よく確認するためには、 3次元イメージを自動的に比較する技術が必要になる。 特に、 製造業では、 製品 品質の向上、 開発期間の短縮により、 高精度かつ高速な計測、 および、 比較が要 求されている。 比較されるイメージは、 C ADにより作成された 3次元自由曲面 からなる設計モデルと、 そのモデルから製造された製品、 あるいは、 部品を 3次 元形状計測器にて測定した点群情報、 もしくは、 その点群情報より作成された断 面ゃサ一フエス情報である。
従来の 3次元イメージ比較では、 ある表面の点ともう一方の表面上の点の距離 をもって、 その差異を評価していた。 これらの点は、 CT (Computed Tomogr aphy) スキャナや 3次元デジタイザにて計測したものと、 与えられた表面をあ る一定の規則に従って分割することによって得るものとがある。 表面と表面との 距離は、 各点から表面への垂線の距離となる。 この距離情報を表面の点にプロッ 卜し、 その値、 あるいは、 距離に比例した色の濃さ、 もしくは、 輝度により可視 ィ匕している。 そこで、 以下、 図 26を参照して具体的に説明する。
図 26は、 物体の 3次元イメージ (CADデータ: 3次元曲面) と 3次元点群 (計測データ) とを示す図である。
図 26のような場合、 3次元イメージの表面を表す方程式、
z=y3+3x2y · · · (1)
とある点 (x。, y。, z。) の距離は, { (x-x0) 2+ (y-y0) 2+ (z - Z 0) 2} 1/2の最小のものである。 この計算を、 点 (x。, y。, z 0) の近傍 無数にある表面の全てに対して行ない、 その中の最小のものをもって点とィメー ジの距離とする。
このようにして行われる従来の比較方法について、 概念図として図 27に示す。 図 27は、 3次元イメージでの従来の比較方法における概念図である。
従来の比較は、 図 27に示すように、 設計モデルである表面データ I ml 01 と、 写真等の表面データ I ml 02とを直接比較することはできない。 このため、 従来の比較では、 表面データ I ml 02の表面を点群に変換し、 点群データ P I ml 02を生成した後に、 表面データ I m 101と点群データ P I m 102との 比較を行っている。
しかしながら、 従来の比較方法では、 表面と点群の高精度比較であれば、 対象 となる点群は無数となる。 また、 実製品の CADモデルでは、 微少な面が多く、 比較すべき面数が膨大となる。 つまり、 実製品の CADモデルでは、 計算すべき 点と表面の組み合わせが無数となる。 その結果、 以下の問題が表出する。
•計算時間が膨大となる。
-ある点に最近接の表面の選択が困難となり、 正しくない値が与えられる場合が ある。
また、 実際の製造業において要求される計測点数は、 以下の表 (1 ) の通りで あり、 どの方法により比較した場合でも、 計測点数がネックとなり、 現実的に実 現が不可能となっている。
Figure imgf000005_0001
ただし、 ブラスティック製品では、 計測点数を上記部品分の表面サイズ (400 X 400mm) の 5倍と仮定している。 発明の開示
本発明はこのような点に鑑みてなされたものであり、 3次元イメージの比較を 高精度かつ高速に行うことができる 3次元イメージ比較ズログラム、 3次元ィメ ージ比較方法、 および 3次元イメージ比較装置を提供することを目的とする。 本発明の第 1の態様では上記課題を解決するために、 たとえば、 図 1のような 3次元ィメージ比較プログラムが提供される。 本発明の 3次元イメージ比較プロ グラムは、 3次元イメージを比較して差異を表示する場合に適用される。 また、 本発明の 3次元イメージ比較プログラムは、 コンピュータに以下の処理を実行さ せる。
予めポクセル化された第 1の 3次元ィメージと第 2の 3次元イメージとを比較 して、 3次元差異イメージを生成する (ステップ S 1 ) 。 また、 前記 3次元差異 イメージの表面に対し、 前記比較よりも細かい精度にて前記第 1の 3次元ィメー ジと前記第 2の 3次元イメージとの 3次元差異細密イメージを生成する (ステツ プ S 2 ) 。 さらに、 前記第 1の 3次元イメージの 3次元差異細密イメージと、 前 記第 2の 3次元ィメ一ジの 3次元差異細密イメージとを比較して、 前記 3次元差 異イメージの差異部の表示形式を決定する (ステップ S 3 ) 。 そして、 前記 3次 元差異ィメ一ジをポクセル比較結果として前記表示形式で表示する (ステップ S このような 3次元ィメ一ジ比較プログラムをコンピュータに実行させれば、 コ ンピュー夕により、 まず、 予めポクセル化された第 1の 3次元イメージと第 2の 3次元イメージとが比較されて、 3次元差異イメージが生成される。 また、 前記 3次元差異ィメ一ジの表面に対し、 前記比較よりも細かい精度にて前記第 1の 3 次元ィメージと前記第 2の 3次元イメージとの 3次元差異細密ィメ一ジが生成さ れる。 さらに、 前記 3次元差異イメージと前記 3次元差異細密イメージとに対し、 ポクセルの差異を表す表示形式が決定される。 そして、 前記 3次元差異イメージ がポクセル比較結果として前記表示形式で表示される。
また、 本発明の第 2の態様では上記課題を解決するために、 コンピュータに、 予めポクセル化された第 1の 3次元イメージと第 2の 3次元ィメージとを比較し て、 3次元差異イメージを生成し、 前記 3次元差異イメージの差異部を部分的に 抽出し、 一致している表面のポクセルに対して、 3次元形状の各面から垂直方向 に不一致ポクセルのポクセル数をカウントし、 前記カウントの数をもとにポクセ ルの差異を表す表示形式を決定し、 前記表示形式を決定した前記 3次元差異ィメ ージを、 前記表示形式により表示する、 処理を実行させることを特徴とする 3次 元イメージ比較プログラムが提供される。
このような 3次元イメージ比較プログラムによれば、 コンピュータにより、 ま ず、 予めポクセル化された第 1の 3次元ィメージと第 2の 3次元ィメ一ジとが比 較されて、 3次元差異イメージが生成される。 また、 前記 3次元差異イメージの 差異部が部分的に抽出され、 一致している表面のポクセルに対して、 3次元形状 の各面から垂直方向に不一致ポクセルのポクセル数がカウントされる。 さらに、 前記力ゥントの数をもとにポクセルの差異を表す表示形式が決定される。 そして、 前記表示形式を決定した前記 3次元差異イメージが、 前記表示形式により表示さ れる。
さらに、 本発明の第 3の態様では上記課題を解決するために、 コンピュータに、 予めポクセル化された第 1の 3次元イメージと第 2の 3次元イメージとを比較し て、 3次元差異イメージを生成し、 前記 3次元差異イメージの差異部における表 面ポクセル数と全体ポクセル数との比率を算出し、 前記 3次元差異イメージに対 し、 前記比率をもとにポクセルの差異を表す表示形式を決定し、 前記表示形式を 決定した前記 3次元差異イメージに対して、 前記比率に対応した所定の表示形式 により前記差異部を表示する、 処理を実行させることを特徴とする 3次元ィメー ジ比較プログラムが提供される。
このような 3次元イメージ比較プログラムによれば、 コンピュータにより、 ま ず、 予めポクセル化された第 1の 3次元イメージと第 2の 3次元イメージとを比 較して、 3次元差異イメージが生成される。 次に、 前記 3次元差異イメージの差 異部における表面ポクセル数と全体ポクセル数との比率が算出され、 前記 3次元 差異イメージに対し、 前記比率をもとにポクセルの差異を表す表示形式が決定さ れる。 そして、 前記表示形式を決定した前記 3次元差異イメージに対して、 前記 比率に対応した所定の表示形式により前記差異部が表示される。
本発明の上記および他の目的、 特徴および利点は本発明の例として好ましい実 施の形態を表す添付の図面と関連した以下の説明により明らかになるであろう。 図面の簡単な説明
図 1は、 第 1の実施の形態における原理構成図である。
図 2は、 第 1の実施の形態における 3次元イメージの比較方法の概念図である。 図 3は、 3次元ィメージ比較プログラムを実行可能なコンピュータのハードウ エア構成例を示す図である。
図 4は、 3次元イメージ比較装置の機能構成を示す機能プロック図である。 図 5は、 3次元イメージ比較プログラムの比較処理の概要を示す図である。 図 6は、 比較した 3次元イメージを、 さらに細かいポクセルで比較する手順を 示す図である。
図 7は、 ポクセル差異数により決定される色を設定する方法において、 共通ポ クセルに対する色決定を示す図である。
図 8は、 ポクセル差異数により決定される色を設定する方法において、 差異ポ クセル (基準イメージ) に対する色決定を示す図である。
図 9は、 ポクセル差異数により決定される色を設定する方法において、 差異ポ クセル (比較対象イメージ) に対する色決定を示す図である。
図 1 0は、 ポクセルの要素の番号と節点の関係を示す図である。 図 1 1は、 3次元イメージ上の節点おける一般的なデータ構造を示す図である。 図 1 2は、 ポクセルの節点におけるポクセル法のデータ構造を示す図である。 図 1 3は、 3次元イメージ比較処理の全体の流れを示すフローチャートである。 図 1 4は、 図 1 3の 3次元イメージ比較処理における比較情報入力処理を示す フローチャートである。
図 1 5は、 図 1 3の 3次元イメージ比較処理におけるポクセル処理を示すフロ 一チヤ一トである。
図 1 6は、 図 1 3の 3次元イメージ比較処理におけるイメージ位置合わせ処理 を示すフローチヤ一トである。
図 1 7は、 図 1 3の 3次元イメージ比較処理における第 1のイメージ比較処理 を示すフローチャートである。
図 1 8は、 第 2の実施の形態において、 3次元イメージ比較プログラムの比較 処理の概要を示す図である。
図 1 9は、 図 1 3の 3次元イメージ比較処理における第 2のイメージ比較処理 を示すフロ一チャートである。
図 2 0は、 第 3の実施の形態において、 3次元イメージ比較プログラムの比較 処理の差異評価 (1 ) を示す図である。
図 2 1は、 第 3の実施の形態において、 3次元イメージ比較プログラムの比較 処理により決定される色 (薄い色) を示す図である。
図 2 2は、 第 3の実施の形態において、 3次元イメージ比較プログラムの比較 処理の差異評価 (2 ) を示す図である。
図 2 3は、 第 3の実施の形態において、 3次元イメージ比較プログラムの比較 処理により決定される色 (濃い色) を示す図である。
図 2 4は、 第 3の実施の形態において、 3次元イメージ比較プログラムの比較 処理の差異評価 (3 ) を示す図である。
図 2 5は、 図 1 3の 3次元イメージ比較処理における第 3のイメージ比較処理 を示すフロ一チャートである。
図 2 6は、 物体の 3次元イメージ (C ADデ一夕: 3次元曲面) と 3次元点群 (計測データ) とを示す図である。 図 2 7は、 3次元イメージでの従来の比較方法における概念図である。 発明を実施するための最良の形態
本発明は、 例えば、 以下のような 3次元イメージに対し、 幾何学的形状を比較 する方法に適用できる。
•有形の物を製造するための C ADデータ
•コンピュータゲーム等の C G(Computer Graphics)画像
•有形物 (工業製品、 自然物、 生物およびその臓器等) を計測した結果
幾何学形状の比較は、 2つの 3次元イメージを共に均一の立方体、 または、 直 方体 (即ち、 ポクセル (Voxel) ) にて表現し、 個々のポクセル毎にブール演算 を施し、 2つのイメージの一致、 不一致を評価することにより行う。 ここで、 2 つのイメージは、 基準点、 あるいは、 特徴的な形状を元に位置合わせされる。 さ らに、 イメージの差異は、 一致しなかったポクセル数にて表される。
また、 C A D等の 3次元サ一フェスを比較する場合は、 一致しなかったポクセ ルの情報を当該のサ一フェス上に表示させる。
このような比較方法について、 以下に示す 3つの処理が考えられる。
• 3次元イメージをポクセル単位で比較し、 3次元イメージの差異を所定の表示 形式で表示する。 つまり、 3次元イメージの差異を点 (ポクセルの色) で表現す る。
·基準ポクセルの各面 (X方法、 y方向、 z方向) から深さ方向への差異ポクセ ル数を、 基準ポクセルの各面それぞれに対して、 所定の表示形式で表示する。 つ まり、 3次元イメージの差異を基準ポクセルの各面 (差異ポクセル数の深さに対 応した色) で表現する。
• 3次元イメージの差異を、 差異部全体の体積と表面積との比率をもとに所定の 表示形式で表示する。 つまり、 3次元イメージの差異を差異部全体 (比率に対応 した色) で表現する。
そこで、 以下の実施の形態では、 まず、 第 1の実施の形態として、 3次元ィメ ージの差異を点 (ポクセルの色) で表現することができる 3次元イメージ比較プ ログラム、 3次元イメージ比較方法、 および 3次元イメージ比較装置について説 明する。 次に、 第 2の実施の形態として、 3次元イメージの差異を基準ポクセル の各面 (差異ポクセル数の深さに対応した色) で表現することができる 3次元ィ メージ比較プログラム、 3次元イメージ比較方法、 および 3次元イメージ比較装 置について説明する。 そして、 第 3の実施の形態として、 3次元イメージの差異 を差異部全体 (比率に対応した色) で表現することができる 3次元イメージ比較 プログラム、 3次元イメージ比較方法、 および 3次元イメージ比較装置について 説明する。
以下、 本発明の実施の形態を図面を参照して説明する。
[第 1の実施の形態]
まず、 第 1の実施の形態について説明する。 第 1の実施の形態は、 3次元ィメ —ジをポクセル単位で比較し、 3次元イメージの差異を所定の表示形式で表示す るものである。
図 1は、 第 1の実施の形態における原理構成図である。 第 1の実施の形態に係 わる 3次元イメージ比較プログラムは、 3次元イメージの比較 Z差異表示をコン ピュー夕に実行させる。 コンピュータによる比較 Z差異表示は、 以下の手川頁で実 行される。 これは、 設計モデルをポクセル化した 3次元イメージ B 1と、 計測さ れたデータをポクセル化した 3次元イメージ B 2とをコンピュータにより比較し、 所定の表示形式により表示する手順を示している。
まず、 コンピュータは、 予めポクセル化された 3次元イメージ B 1と、 3次元 イメージ B 2とを比較して、 3次元差異イメージ V B 1を生成する (ステップ S 1 ) 。
また、 コンピュータは、 3次元差異イメージ V B 1の表面に対し、 ステップ S 1の比較よりも細かい精度にて 3次元イメージ B 1と 3次元イメージ B 2との 3 次元差異細密イメージ V B 2を生成する (ステップ S 2 ) 。
次に、 コンピュータは、 3次元イメージ B 1の 3次元差異細密イメージと 3次 元イメージ B 2の 3次元差異細密イメージとを比較して、 3次元差異イメージ V B 1の差異部の表示形式を決定する (ステップ S 3 ) 。 なお、 この表示形式とは、 赤、 緑、 および青の光 3原色を基本とした加法混色 (R G B) である。 また、 表 示形式は、 シアン、 マゼンダ、 および黄の補色 3原色を基本とした減法混色 (C MY) であってもよい。 ここで、 減法混色を使用する場合には、 シアン、 マゼン ダ、 および黄の補色 3原色の他に、 黒も追加することもできる。
そして、 コンピュータは、 表示形式を決定した 3次元差異イメージ V B 1を、 ボクセル比較結果として所定の表示形式でディスプレイ等に表示する (ステップ S 4 ) 。 なお、 ポクセル比較結果は、 ディスプレイ等への表示だけでなく、 プリ ン夕ゃプロッタ等への印刷も行うこともできる。
このように第 1の実施の形態では、 まず、 コンピュータにより、 予めポクセル 化された 3次元イメージ B 1と、 3次元イメージ B 2とが比較されて、 3次元差 異イメージ V B 1が生成される。
また、 コンピュータにより、 3次元差異イメージ V B 1の表面に対し、 ステツ プ S 1の比較よりも細かい精度にて 3次元イメージ B 1と 3次元イメージ B 2と の 3次元差異細密イメージ V B 2が生成される。
次に、 コンピュータにより、 3次元イメージ B 1の 3次元差異細密イメージと 3次元イメージ B 2の 3次元差異細密イメージとを比較して、 3次元差異ィメー ジ V B 1の差異部の表示形式が決定される。
そして、 コンピュータにより、 表示形式を決定した 3次元差異イメージ V B 1 が、 ポクセル比較結果として所定の表示形式でディスプレイ等に表示される。 これにより、 3次元イメージの比較において、 比較処理を高精度かつ高速に行 うことが可能となる。
次に、 第 1の実施の形態における 3次元イメージの比較方法の概念について、 図 2を参照して説明する。
図 2は、 第 1の実施の形態における 3次元ィメ一ジの比較方法の概念図である。 第 1の実施の形態における比較方法は、 図 2に示すように、 設計モデルである 表面データ I m l 1と、 他の表面デ一夕 I m l 2とを直接比較することはできな い。 このため、 第 1の実施の形態における比較方法では、 表面データ I m l 1, I m l 2それぞれ囲まれる体積をポクセルに変換し、 ポクセル B 1 1, B 1 2を それぞれ生成した後に、 ポクセル B 1 1とポクセル B 1 2との比較を行っている これにより、 2つの 3次元ィメ一ジの比較を容易に行うことができる。
以下に、 本発明における第 1の実施の形態について具体的に説明する。 図 3は、 3次元ィメ一ジ比較プログラムを実行可能なコンピュータのハードウ エア構成例を示す図である。 コンピュータ 1 0 0は、 C P U(Central Processing Unit) 101によって装置全体が制御されている。 C PU 101には、 バス 107を介して R AM(Random Access Memory) 102、 ハードディスクド ライブ (HDD:Hard Disk Drive) 103、 グラフィック処理装置 104、 入力 インタフェース 105、 および通信インタフェース 106が接続されている。
RAM 102には、 C P U 101に実行させる〇 S (Operating System)のプ ログラムやアプリケーションプログラムの少なくとも一部が一時的に格納される。 また、 RAMI 02には、 CPU101による処理に必要な各種データが格納さ れる。 HDD 103には、 OSや、 3次元イメージ比較プログラムなどのアプリ ケ一ションプログラムが格納される。
グラフィック処理装置 104には、 モニタ P 111が接続されている。 また、 グラフィック処理装置 104は、 CPU 101からの命令に従って、 画像をモニ 夕 P 11 1の画面に表示させる。 入力インタフェース 105には、 キ一ポ一ド P 112とマウス P 113とが接続されている。 また、 入力インタフェース 105 は、 キーボード P 112やマウス P 113から送られてくる信号を、 バス 107 を介して CPU 101に送信する。
通信インタフェース 106は、 ネットワーク 110に接続されている。 通信ィ ン夕フェース 106は、 ネットヮ一ク 110を介して、 他のコンピュータとの間 でデ一夕の送受信を行う。
以上のようなハ一ドウエア構成のコンピュータ 100に 3次元イメージ比較プ ログラムを実行させることで、 コンピュータ 100が 3次元イメージ比較装置と して機能する。 以下、 コンピュータが 3次元イメージ比較プログラムを実行する ことで実現される 3次元イメージ比較装置の処理機能について説明する。
図 4は、 3次元イメージ比較装置の機能構成を示す機能ブロック図である。 図 4によると、 本発明における 3次元イメージ比較装置は、 比較情報を入力す る比較情報入力部 10、 入力された 3次元イメージからポクセルを生成するポク セル処理部 20、 ポクセル化された 3次元ィメージを位置合わせするィメ一ジ位 置合わせ処理部 30、 ポクセルの差異を評価して表示するポクセル差異評価処理 部 4 0、 および、 表示するポクセルを重ね合わせるポクセル重ね合わせ処理部 5 0から構成される。 ここで、 比較情報入力部 1 0は、 自身の機能とは別に、 計測 精度入力部 1 1、 計測精度記憶部 1 2、 3次元サ一フェス入力部 1 3、 イメージ 記憶部 1 4、 点群処理切替部 1 5、 点群 (表面) 入力部 1 6、 点群 (表面) 記憶 部 1 7、 点群 (体積) 入力部 1 8、 および点群 (体積) 記憶部 1 9を有する。 ま た、 ポクセル処理部 2 0は、 自身の機能とは別に、 ポクセル発生部 2 1、 ポクセ ル発生部 2 2、 およびポクセル記憶部 2 3を有する。 以下、 これらの機能の詳細 を列記する。
比較情報入力部 1 0は、 ポクセル処理部 2 0、 3次元 C ADイメージ D B 1 1、 および 3次元計測イメージ D B 1 2に接続され、 比較情報である 3次元イメージ の C ADデータと計測データとを入力する。 ここで、 比較情報入力部 1 0は、 ュ 一ザによる 3次元イメージの入力指示を受け付ける。 また、 比較情報入力部 1 0 は、 受け付けた 3次元イメージのイメージ種別を判断する。 そして、 比較情報入 力部 1 0は、 イメージ記憶部 1 4、 点群 (表面) 記憶部 1 7、 あるいは点群 (体 積) 記憶部 1 9にて 3次元イメージが記憶されたかどうかを判定する。
計測精度入力部 1 1は、 計測精度記憶部 1 2に接続され、 ユーザによる計測精 度の入力を受け付ける。 なお、 計測精度とは、 ユーザが任意に決めることができ、 ポクセル化する 3次元イメージの最小単位である。 また、 計測精度には、 たとえ ば、 製造される製品の要求精度 (0.01mm等) を使用することが好ましい。
計測精度記憶部 1 2は、 計測精度入力部 1 1に接続され、 計測精度入力部 1 1 にて受け付けた計測精度を記憶する。
3次元サ一フェス入力部 1 3は、 イメージ記憶部 1 4と 3次元 C ADイメージ D B 1 1に接続され、 ユーザの指示した 3次元イメージのィメ一ジ種別が 3次元 サ一フェスであれば、 3次元 C ADイメージ D B 1 1から 3次元サ一フェスの入 力受け付けを行う。
イメージ記憶部 1 4は、 3次元サ一フェス入力部 1 3に接続され、 3次元サ一 フェス入力部 1 3にて受け付けた 3次元サ一フェスを、 内部の記憶装置に記憶す る。
点群処理切替部 1 5は、 点群 (表面) 入力部 1 6と点群 (体積) 入力部 1 8に 接続され、 点群処理の切り替えを行う。 ここで、 点群処理切替部 1 5は、 ユーザ の指示した 3次元イメージのイメージ種別が計測表面データ、 あるいは計測断面 データの場合に、 点群 (表面) 入力部 1 6と点群 (体積) 入力部 1 8との切り替 えを行う。
点群 (表面) 入力部 1 6は、 点群処理切替部 1 5、 点群 (表面) 記憶部 1 7、 および 3次元計測イメージ D B 1 2に接続され、 計測表面デ一夕の入力受け付け を行う。 ここで、 点群 (表面) 入力部 1 6は、 ユーザの指示した 3次元イメージ のイメージ種別が計測表面データであれば、 3次元計測イメージ D B 1 2から計 測表面デ一夕の入力受け付けを行う。
点群 (表面) 記憶部 1 7は、 点群 (表面) 入力部 1 6に接続され、 計測表面デ —夕を記憶する。 ここで、 点群 (表面) 記憶部 1 7は、 点群 (表面) 入力部 1 6 にて受け付けた計測表面データを、 内部の記憶装置に記憶する。
点群 (体積) 入力部 1 8は、 点群処理切替部 1 5、 点群 (体積) 記憶部 1 9、 および 3次元計測イメージ D B 1 2に接続され、 計測断面デ一夕の入力受け付け を行う。 ここで、 点群 (体積) 入力部 1 8は、 ユーザの指示した 3次元イメージ のイメージ種別が計測断面データであれば、 3次元計測イメージ D B 1 2から計 測断面データの入力受け付けを行う。
点群 (体積) 記憶部 1 9は、 点群 (体積) 入力部 1 8に接続され、 計測断面デ —夕を記憶する。 ここで、 点群 (体積) 記憶部 1 9は、 点群 (体積) 入力部 1 8 にて受け付けた計測断面データを、 内部の記憶装置に記憶する。
ポクセル処理部 2 0は、 比較情報入力部 1 0とイメージ位置合わせ処理部 3 0 とに接続され、 比較情報入力部 1 0にて入力された 3次元イメージからポクセル によるイメージを生成する。 ここで、 ポクセル処理部 2 0は、 比較情報入力部 1 0にて記録された比較情報を選択する。 また、 ポクセル処理部 2 0は、 比較情報 入力部 1 0にて記録された計測精度をもとに、 ポクセルサイズを決定する。 そし て、 ポクセル処理部 2 0は、 比較情報入力部 1 0にて記憶した 3次元イメージの イメージ種別を受け付けて、 次の処理を判断する。
この判断により、 ポクセル処理部 2 0は、 イメージ種別が 3次元サ一フェスで あれば、 比較情報入力部 1 0のイメージ記憶部 1 4から 3次元サーフェスの選 択-読み込みを行い、 この 3次元サ一フェスをポクセル発生部 2 1へ渡す。 また、 ポクセル処理部 2 0は、 イメージ種別が計測表面データであれば、 比較 情報入力部 1 0の点群 (表面) 記憶部 1 7から計測表面データの選択,読み込み を行い、 この計測表面データをポクセル発生部 2 2へ渡す。
さらに、 ポクセル処理部 2 0は、 イメージ種別が計測断面データであれば、 比 較情報入力部 1 0の点群 (体積) 記憶部 1 9から計測断面データの選択 ·読み込 みを行い、 この計測断面データをポクセル発生部 2 2へ渡す。
そして、 ポクセル処理部 2 0は、 ポクセル発生部 2 1、 ポクセル発生部 2 2に てポクセルが生成されたかどうかを判定する。
ポクセル発生部 2 1は、 渡された 3次元サ一フェスをもとに、 ポクセルによる イメージを生成する。
ポクセル発生部 2 2は、 渡された計測表面データをもとに、 ポクセルによるィ メージを生成する。 また、 ポクセル発生部 2 2は、 渡された計測断面データをも とに、 ポクセルによるイメージを生成する。
イメージ位置合わせ処理部 3 0は、 ポクセル処理部 2 0とポクセル差異評価処 理部 4 0に接続され、 ポクセル処理部 2 0にてポクセルに変換された 2つのィメ ージを位置合わせする。 ここで、 イメージ位置合わせ処理部 3 0は、 予め 2つの 3次元イメージに付加されている基準点を選択する。 なお、 この 2つの 3次元ィ メージとは、 C ADデータである 3次元サーフェスと、 計測データである計測表 面データあるいは計測断面データを指す。 また、 イメージ位置合わせ処理部 3 0 は、 選択した基準点からのオフセットを計算する。 さらに、 イメージ位置合わせ 処理部 3 0は、 計算したオフセットをもとに、 ポクセル変換した 3次元サ一フエ スを移動する。 そして、 イメージ位置合わせ処理部 3 0は、 移動した 3次元サ一 フェスを内部の記憶装置に記憶する。
ポクセル差異評価処理部 4 0は、 ィメ一ジ位置合わせ処理部 3 0とポクセル重 ね合わせ処理部 5 0とに接続され、 ィメ一ジ位置合わせ処理部 3 0にて位置合わ せされた 2つのイメージを比較して評価し、 評価した 2つのィメージを所定の表 示形式で表示する。 なお、 所定の表示形式とは、 たとえば、 R G Bのうちの各色 の濃さをイメージの差異により変更して表示することである。 ここで、 ポクセル 差異評価処理部 4 0は、 表面ポクセルを選択する。 この選択によりポクセル差異 評価処理部 4 0は、 選択されたポクセルが表面ポクセルかどうかを判断する。 そ して、 ポクセル差異評価処理部 4 0は、 この判断により表面ポクセルが選択され れば、 評価精度での再ポクセル化を行う。 すなわち、 評価精度に応じて、 表面ポ クセルをさらに細かくポクセル化する。 なお、 評価精度は、 ユーザにより予め決 めておくことができる。 また、 評価精度は、 計測精度の数倍から段階的に細かく することもできる。
また、 ポクセル差異評価処理部 4 0は、 ポクセルに変換された 3次元サ一フエ スと計測表面 Z断面データとをもとに、 差異ポクセルを比較して計数する。 なお、 この比較処理の詳細については後述する。
さらに、 ポクセル差異評価処理部 4 0は、 計数された差異ポクセルをもとに、 表示形式である色の濃さを決定し、 ポクセルの色付けを行う。
そして、 ポクセル差異評価処理部 4 0は、 ポクセルが色付けされたかどうかを 判定する。 すなわち、 ポクセル差異評価処理部 4 0は、 ポクセル化された 3次元 ィメージが比較されたかどうかを判定する。
ポクセル重ね合わせ処理部 5 0は、 ポクセル差異評価処理部 4 0に接続され、 ポクセル差異評価処理部 4 0にて表示された複数のイメージを重ね合わせる。
3次元 C ADイメージ D B 1 1は、 比較情報入力部 1 0に接続され、 C ADデ —夕である 3次元サーフェスを格納している。
3次元計測イメージ D B 1 2は、 比較情報入力部 1 0に接続され、 計測データ である計測表面デー夕と計測断面データとを格納している。
次に、 第 1の実施の形態における 3次元ィメ一ジ比較プログラムの比較処理に ついて、 図 5〜図 9を参照して具体的に説明する。
第 1の実施の形態は、 まず、 3次元形状をポクセルにより比較することで、 計 算時間の低減を図る。 例えば、 3次元の曲面に囲まれている C ADイメージの形 状では、 微細で均一な立方体 (ポクセル) にてイメージを分割し、 その内部のポ クセルをオンに、 外部をオフにすることによって 3次元ィメ一ジを表現する。 ま た、 X線 C Tスキャナにて計測されたイメージは、 その X線の密度が高いポクセ ルをオンに、 低いものをオフとする。 ここで、 ポクセルの比較処理の概要につい ては、 図 5を参照して説明する。
図 5は、 3次元イメージ比較プログラムの比較処理の概要を示す図である。 図 5のように、 比較処理は、 2つのイメージを 3次元イメージ B 21と 3次元 イメージ B 22のようにポクセル化して、 その重なっている部分 (3次元差異ィ メ一ジ VB 21) をブーリアン演算することで、 その差異を視覚的に表現するこ とができる。 これにより、 比較の演算量は従来法の 1Z数十となる。 なお、 この ブーリアン演算では、 たとえば、 3次元イメージ B 21のポクセルと、 3次元ィ メ一ジ B 22のポクセルとを排他的論理和 (XOR) することにより、 不一致ポ クセルを抜き出すことができる。
次に、 比較する 2つの 3次元イメージを評価精度より粗いポクセルにて表現し た後、 さらに、 この粗ポクセル毎に精細なポクセルを発生させ、 2つのイメージ の差異をその差異ポクセルの数にて表現する。 そして、 視覚的に差異を認識でき るように、 この数値を片方のイメージ上にて重ね合わせる。 この仕組みについて は、 図 6を参照して説明する。
図 6は、 比較した 3次元イメージを、 さらに細かいポクセルで比較する手順を 示す図である。
図 6のように、 まず、 粗いポクセル (3次元差異イメージ VB 21) を発生さ せる。 そして、 表面を含むポクセルを細密に分割 (3次元差異細密イメージ VB 22) し、 その差異ポクセルの数にて、 2つのイメージの相違を評価する。 なお、 図 6に示す例では、 2つのイメージの差異は " 8" となる。 これは、 たとえば、 粗いポクセル内において、 イメージ I 21, I 22の 3次元差異細密イメージ V B 22での細密ポクセルの差を示している。
次に、 差異数を元にポクセルの色の濃さを決定する。 結果として、 全体の差異 が可視的に理解可能となる。 たとえば、 図 6に示すように、 ポクセル Bx 1 1, Bx 12, Bx l 3では、 この部分の色 (例えば、 青色の濃さ) が決定している ので、 イメージ I 22が差異と判断できる。 また、 ポクセル Bx 16では、 この 部分の色 (例えば、 赤色の濃さ) が決定しているので、 イメージ I 21が差異で あることが分かる。 そして、 ポクセル Bx l 4, Bx l 5では、 この部分の色が 決定していないので、 差異なしと判断できる。 このように、 図 6に示す例では、 ポクセルをさらに相似比 1/4のポクセルで 細分化している。 つまり、 3次元では、 最大 64 (4X4X4) ポクセルの差異 となる。 この差異数により決定される色を設定する方法については、 図 7〜図 9 を参照して説明する。
図 7は、 ポクセル差異数により決定される色を設定する方法において、 共通ポ クセルに対する色決定を示す図である。
図 7によると、 ポクセルが共通である場合、 すなわち、 差異がない場合、 たと えば、 設定位置 ST 1における RGBの 3色とも、 輝度を最大値 (例えば、 25 5, 255, 255) にして設定する。
図 8は、 ポクセル差異数により決定される色を設定する方法において、 差異ポ クセル (基準イメージ) に対する色決定を示す図である。
図 8によると、 比較対象ィメ一ジに対して基準ィメージのポクセルが凹となつ ている場合、 たとえば、 設定位置 ST 2における緑 (G) の値を差異ポクセル数 に応じて減少させる。
図 9は、 ポクセル差異数により決定される色を設定する方法において、 差異ポ クセル (比較対象イメージ) に対する色決定を示す図である。
図 9によると、 基準イメージに対して比較対象イメージのポクセルが凹となつ ている場合、 たとえば、 設定位置 ST 3における赤 (R) の値を差異ポクセル数 に応じて減少させる。
以上、 図 7〜図 9に示すように、 色の値は、 それぞれのイメージのポクセル数 に比例して決定する。
次に、 3次元イメージ比較装置で使われるデータの構造について、 図 10〜図 12を参照して説明する。
図 10は、 ポクセルの要素の番号と節点の闋係を示す図である。
図 10によると、 3次元ィメ一ジから変換されたポクセルは、 たとえば、 8つ のポクセル B X 1〜B X 8とする。 この場合、 ポクセル B x lは、 節点 (articu lation point) (1, 2, 5, 4, 10, 11, 14, 13) から構成される。 なお、 節点の I噴番は、 上下の面それぞれにて、 反時計回りにカウントされる。 そ して、 ポクセルのデータ構造には、 このような各節点に対して、 座標情報が登録 されている。 そこで、 この座標情報の詳細について、 図 11と図 12を参照して 説明する。
図 11は、 3次元イメージの節点における一般的なデ一夕構造を示す図である。 なお、 この図 11は、 典型的な有限要素解析におけるデータ構造を示している。 図 11によると、 ポクセルのデ一夕構造 D 1は、 データの種別を示すデータ種 別エリア D 1 1、 各デ一夕種別におけるシリアル番号を示すデータ番号エリア D 12、 各データ種別におけるデータの位置情報を示す位置情報エリア D 13から 構成されている。
デ一夕種別エリア D 1 1は、 たとえば、 図 11に示すように "GR I D" 、 "CHEXA" の 2つからなる。 ここで、 "GR I D" は節点を示し、 "CHE
XA" は 3次元イメージを分割した要素を示している。
データ番号エリア D 12には、 たとえば、 図 11に示すようにデータ種別 "G
R I D" に対して、 1〜27のシリアル番号が付けられている。 このシリアル番 号 1〜27は、 図 10にて前述した節点番号である。 また、 データ番号エリア D 12には、 たとえば、 図 11に示すようにデータ種別 "CHEXA" に対して、
1〜 8のシリアル番号が付けられている。 このシリアル番号 1〜 8は、 図 10に て前述したポクセルの要素番号である。
位置情報エリア D13には、 たとえば、 図 11に示すようにデータ種別 "GR
I D" のデータ番号 1に対して、 "0. 000000, 0. 000000, 50. 00000" の各節点における座標情報が設定されている。 これは、 各デ一タ番 号 ( "GR I D" に対しては、 節点 I Dと称する) に続き、 x, y, zの座標位 置を浮動少数点にて表現している。 以下、 図のように設定される。 また、 位置情 報エリア D 13には、 たとえば、 図 11に示すようにデータ種別 "CHEXA" のデータ番号 1に対して、 "17, 21, 27, 22, 7, 10, 26, 1 " の各要素における節点 I Dが設定されている。 これは、 各データ番号 ( "CHE
XA" に対しては、 要素 IDと称する) に続き、 各要素を定義している。 以下、 図のように設定される。
図 12は、 ポクセルの節点におけるポクセル法のデータ構造を示す図である。 なお、 この図 12は、 ポクセル法にて 3次元を表現した場合のデータ構造を示し ている。
図 12によると、 ポクセルのデータ構造 D 2は、 予め定義された別々のレコ一 ドに対して、 要素の大きさと、 各要素の節点とをそれぞれ設定する。 なお、 要素 の大きさには、 X, y, zの辺の長さを浮動少数点にて表現する。 また、 各要素 の節点には、 各要素を定義する節点を設定する。 但し、 前述した要素 IDについ ては、 正規的に並んでいるため定義する必要はない。 ここで例えば、 図 12に示 すように、 レコード (RECORD 7) には、 指数表現 (浮動少数点表現) によ る要素の大きさ " 2. 50000000 e + 001, 2. 50000000 e + 001, 2. 50000004 e + 001" が設定される。 また、 レコード (R E CORD 9) には、 各要素の節点 " 1 , 2, 5, 4, 10, 11, 14, 1 3"〜 ", 14, 15, 18, 17, 23, 24, 27, 26" が設定される。 なお、 上記要素の大きさのうち、 "e" とは、 指数関数の底を示す。 たとえば、 上記の例では、 "e" が 10を示し、 "e + 001" が 10の 1乗を示している。 以上のような構成の 3次元ィメ一ジ比較装置により、 3次元イメージの比較処 理が行われる。 以下、 この 3次元イメージ比較処理の動作について、 具体的に説 明する。
図 13は、 3次元イメージ比較処理の全体の流れを示すフローチャートである。 この処理は、 3次元イメージ比較装置に対して、 電源投入、 プログラムの起動、 あるいは、 予め決められた任意のタイミングを契機に動作し、 CPU101にて 実行させる処理である。 以下、 図 13に示す処理をステップ番号に沿って説明す る。 なお、 本フ口一チャートにおける各機能の名称については、 図 4をもとに説 明する。
[ステツプ S 10 ] 比較情報入力部 10は、 比較情報である 3次元イメージ の CADデータと計測データとを、 入力する比較情報入力処理を行う。 なお、 比 較情報入力処理の詳細については、 図 14にて後述する。
[ステップ S 20] ポクセル処理部 20は、 ステップ S 10にて入力された 3次元イメージから、 ポクセルによるイメージを生成するポクセル処理を行う。 なお、 ポクセル処理の詳細については、 図 1 5にて後述する。
[ステップ S 30] イメージ位置合わせ処理部 30は、 ステップ S 20にて ポクセルに変換された 2つのイメージを、 位置合わせするイメージ位置合わせ処 理を行う。 なお、 イメージ位置合わせ処理の詳細については、 図 1 6にて後述す る。
[ステップ S 4 0 ] ポクセル差異評価処理部 4 0は、 ステップ S 3 0にて位 置合わせされた 2つのイメージを、 比較して評価するイメージ比較処理を行う。 なお、 イメージ比較処理の詳細については、 図 1 7にて後述する。
[ステップ S 5 0 ] ポクセル差異評価処理部 4 0は、 ステップ S 4 0にて比 較 ·評価された複数のイメージを、 所定の表示形式で表示する。 なお、 所定の表 示形式とは、 たとえば、 R G Bのうちの各色の濃さをイメージの差異により変更 して表示することである。
[ステップ S 6 0 ] ポクセル重ね合わせ処理部 5 0は、 ステップ S 5 0にて 表示された複数のイメージを重ね合わせる。 つまり、 比較結果を 3次元イメージ に重ね合わせる。 たとえば、 比較結果である 3次元差異イメージは、 ポクセルィ匕 された 2つのイメージをもとに、 排他的論理和により生成され差異部のみを可視 化することができる。 この場合、 ポクセル重ね合わせ処理部 5 0は、 ポクセル化 された 2つのイメージ、 3次元差異イメージ、 あるいは 3次元差異細密イメージ を重ね合わせることにより、 元のイメージと差異部を色分けして表示することが できる。
図 1 4は、 図 1 3の 3次元イメージ比較処理における比較情報入力処理を示す フローチャートである。 図 1 3のステップ S 1 0 (比較情報入力処理) が実行さ れると、 以下の流れに従って処理が行われる。 なお、 この比較情報入力処理は、 3次元イメージ比較装置の比較情報入力部 1 0により行われる。
[ステツプ S 1 0 1 ] 比較情報入力部 1 0の計測精度入力部 1 1は、 ユーザ による計測精度の入力を受け付ける。
[ステップ S 1 0 2 ] 比較情報入力部 1 0の計測精度記憶部 1 2は、 ステツ プ S 1 0 1にて受け付けた計測精度を記憶する。
[ステップ S 1 0 3 ] 比較情報入力部 1 0は、 ユーザによる 3次元イメージ の入力指示を受け付ける。
[ステップ S 1 0 4 ] 比較情報入力部 1 0は、 ステップ S 1 0 3にて受け付 けた 3次元イメージのイメージ種別を判断する。 ここで、 比較情報入力部 1 0は、 イメージ種別が 3次元サ一フェスであればステップ S 1 0 5へ進み、 イメージ種 別が計測表面データであればステップ S 1 0 7へ進み、 あるいは、 イメージ種別 が計測断面データであればステップ S 1 0 9へ進む。 なお、 イメージ種別が計測 表面データ、 あるいは計測断面デ一夕である場合には、 比較情報入力部 1 0の点 群処理切替部 1 5にて、 点群 (表面) 入力部 1 6と点群 (体積) 入力部 1 8との 切り替えが行われる。
[ステップ S 1 0 5 ] 比較情報入力部 1 0の 3次元サーフェス入力部 1 3は、 ステップ S 1 0 4にてイメージ種別が 3次元サ一フェスであるので、 3次元 C A Dイメージ D B 1 1から 3次元サーフェスの入力受け付けを行う。
[ステップ S 1 0 6 ] 比較情報入力部 1 0のイメージ記憶部 1 4は、 ステツ プ S 1 0 5にて受け付けた 3次元サ一フェスを、 内部の記憶装置に記憶する。
[ステップ S 1 0 7 ] 比較情報入力部 1 0の点群 (表面) 入力部 1 6は、 ス テツプ S 1 0 4にてイメージ種別が計測表面データであるので、 3次元計測ィメ —ジ D B 1 2から計測表面デ一夕の入力受け付けを行う。
[ステップ S 1 0 8 ] 比較情報入力部 1 0の点群 (表面) 記憶部 1 7は、 ス テツプ S 1 0 7にて受け付けた計測表面データを、 内部の記憶装置に記憶する。
[ステップ S 1 0 9 ] 比較情報入力部 1 0の点群 (体積) 入力部 1 8は、 ス テツプ S 1 0 4にてイメージ種別が計測断面データであるので、 3次元計測ィメ ージ D B 1 2から計測断面データの入力受け付けを行う。
[ステップ S 1 1 0 ] 比較情報入力部 1 0の点群 (体積) 記憶部 1 9は、 ス テツプ S 1 0 9にて受け付けた計測断面データを、 内部の記憶装置に記憶する。
[ステップ S 1 1 1 ] 比較情報入力部 1 0は、 ステップ S 1 0 6、 ステップ S 1 0 8、 あるいはステップ S 1 1 0にて 3次元イメージが記憶されたかどうか を判定する。 ここで、 比較情報入力部 1 0は、 3次元イメージが記録されていな ければ、 ステップ S 1 0 4へ戻り同様の処理を繰り返し、 3次元イメージが記録 されていれば、 本処理を終了して図 1 3のステップ S 1 0へ復帰する。
図 1 5は、 図 1 3の 3次元イメージ比較処理におけるポクセル処理を示すフロ 一チャートである。 図 1 3のステップ S 2 0 (ポクセル処理) が実行されると、 以下の流れに従って処理が行われる。 なお、 このポクセル処理は、 3次元ィメ一 ジ比較装置のポクセル処理部 2 0により行われる。
[ステップ S 2 0 1 ] ポクセル処理部 2 0は、 ステップ S 1 0にて記憶され た比較情報を選択する。
[ステップ S 2 0 2 ] ポクセル処理部 2 0は、 ステップ S 1 0にて記憶され た計測精度をもとに、 ポクセルサイズを決定する。
[ステップ S 2 0 3 ] ポクセル処理部 2 0は、 ステップ S 1 0にて記憶され た 3次元イメージの種別を受け付ける。
[ステップ S 2 0 4 ] ポクセル処理部 2 0は、 ステップ S 2 0 3にて受け付 けた 3次元イメージのイメージ種別を判断する。 ここで、 ポクセル処理部 2 0は、 イメージ種別が 3次元サ一フェスであればステップ S 2 0 5へ進み、 イメージ種 別が計測表面デ一夕であればステップ S 2 0 7へ進み、 あるいは、 イメージ種別 が計測断面データであればステップ S 2 0 9へ進む。
[ステップ S 2 0 5 ] ポクセル処理部 2 0は、 ステップ S 2 0 4にてィメ一 ジ種別が 3次元サ一フェスであるので、 比較情報入力部 1 0のイメージ記憶部 1 4から 3次元サ一フェスの選択 ·読み込みを行い、 この 3次元サーフェスをポク セル発生部 2 1へ渡す。
[ステップ S 2 0 6 ] ポクセル処理部 2 0のポクセル発生部 2 1は、 ステツ プ S 2 0 5にて渡された 3次元サ一フェスをもとに、 ポクセルによる 3次元ィメ ージを生成する。
[ステップ S 2 0 7 ] ポクセル処理部 2 0は、 ステップ S 2 0 4にてィメ一 ジ種別が計測表面デ一夕であるので、 比較情報入力部 1 0の点群 (表面) 記憶部 1 7から計測表面データの選択 ·読み込みを行い、 この計測表面データをポクセ ル発生部 2 2へ渡す。
[ステップ S 2 0 8 ] ポクセル処理部 2 0のポクセル発生部 2 2は、 ステツ プ S 2 0 7にて渡された計測表面デ一夕をもとに、 ポクセルによる 3次元ィメー ジを生成する。
[ステップ S 2 0 9 ] ポクセル処理部 2 0は、 ステップ S 2 0 4にてィメ一 ジ種別が計測断面データであるので、 比較情報入力部 1 0の点群 (体積) 記憶部 1 9から計測断面データの選択 '読み込みを行い、 この計測断面データをポクセ ル発生部 2 2へ渡す。
[ステップ S 2 1 0 ] ポクセル処理部 2 0のポクセル発生部 2 2は、 ステツ プ S 2 0 9にて渡された計測断面データをもとに、 ポクセルによる 3次元ィメ一 ジを生成する。
[ステップ S 2 1 1 ] ポクセル処理部 2 0は、 ステップ S 2 0 6、 ステップ S 2 0 8、 あるいはステップ S 2 1 0にて、 ポクセルによる 3次元イメージが生 成されたかどうかを判定する。 ここで、 ポクセル処理部 2 0は、 ポクセルが生成 されていなければ、 ステップ S 2 0 4へ戻り同様の処理を繰り返し、 ポクセルが 生成されていれば、 本処理を終了して図 1 3のステップ S 2 0へ復帰する。
図 1 6は、 図 1 3の 3次元イメージ比較処理におけるイメージ位置合わせ処理 を示すフローチャートである。 図 1 3のステップ S 3 0 (イメージ位置合わせ処 理) が実行されると、 以下の流れに従って処理が行われる。 なお、 このイメージ 位置合わせ処理は、 3次元イメージ比較装置のイメージ位置合わせ処理部 3 0に より行われる。
[ステップ S 3 0 1 ] イメージ位置合わせ処理部 3 0は、 予め 2つの 3次元 イメージに付加されている基準点を選択する。 なお、 この 2つの 3次元イメージ とは、 C ADデータである 3次元サーフェスと、 計測デ一夕である計測表面デー タあるいは計測断面デ一夕を指す。
[ステップ S 3 0 2 ] イメージ位置合わせ処理部 3 0は、 ステップ S 3 0 1 にて選択した基準点からのオフセットを計算する。
[ステップ S 3 0 3 ] イメージ位置合わせ処理部 3 0は、 ステップ S 3 0 2 にて計算されたオフセットをもとに、 ポクセル変換された 3次元イメージを移動 する。
[ステップ S 3 0 4 ] イメージ位置合わせ処理部 3 0は、 ステップ S 3 0 3 にて移動した 3次元ィメージを内部の記憶装置に記憶し、 本処理を終了して図 1 3のステップ S 3 0へ復帰する。
次に、 3次元イメージ比較処理について説明するが、 本発明の 3次元イメージ 比較処理は、 第 1〜第 3の実施の形態に対応して複数ある。 このため、 ここでは、 第 1の実施の形態における 3次元ィメ一ジ比較処理を、 第 1のィメージ比較処理 として説明する。
図 1 7は、 図 1 3の 3次元イメージ比較処理における第 1のイメージ比較処理 を示すフローチャートである。 図 1 3のステップ S 4 0 (イメージ比較処理) が 実行されると、 以下の流れに従って処理が行われる。 なお、 この第 1のイメージ 比較処理は、 3次元イメージ比較装置のポクセル差異評価処理部 4 0により行わ れる。 また、
[ステップ S 4 1 1 ] ポクセル差異評価処理部 4 0は、 表面ポクセルを選択 する。 なお、 表面ポクセルの選択は、 たとえば、 ポクセル化された 3次元ィメ一 ジに対して、 表面ポクセルの内側がオンであり、 外側がオフであることにより容 易に判断できる。
[ステップ S 4 1 2 ] ポクセル差異評価処理部 4 0は、 ステップ S 4 1 1に て選択されたポクセルが表面ポクセルかどうかを判断する。 ここで、 ポクセル差 異評価処理部 4 0は、 表面ポクセルが選択されていれば、 ステップ S 4 1 3へ進 み、 表面ポクセルが選択されていなければ、 ステップ S 4 1 1へ戻る。
[ステップ S 4 1 3 ] ポクセル差異評価処理部 4 0は、 ステップ S 4 1 2に て表面ポクセルが選択されたので、 評価精度での再ポクセル化を行う。 すなわち、 評価精度に応じて、 表面ポクセルをさらに細かくポクセル化する。
[ステップ S 4 1 4 ] ポクセル差異評価処理部 4 0は、 ポクセルに変換され た 3次元サーフェスと計測表面 Z断面デ一夕とをもとに、 差異ポクセルを比較し て計数する。
[ステップ S 4 1 5 ] ポクセル差異評価処理部 4 0は、 ステップ S 4 1 4に て計数された差異ポクセルをもとに、 表示形式である色の濃さを決定し、 ポクセ ルの色付けを行う。
[ステップ S 4 1 6 ] ポクセル差異評価処理部 4 0は、 ステップ S 4 1 5に てポクセルが色付けされたかどうかを判定する。 すなわち、 ポクセル差異評価処 理部 4 0は、 ポクセル化された 3次元イメージが比較されたかどうかを判定する。 ここで、 ポクセル差異評価処理部 4 0は、 ポクセルが色付けされていなければ、 ステップ S 4 1 1へ戻り同様の処理を繰り返し、 ポクセルが色付けされていれば、 本処理を終了して図 1 3のステップ S 4 0へ復帰する。
このように第 1の実施の形態では、 まず、 コンピュータにより、 予めポクセル 化された 3次元イメージ B 2 1と、 3次元イメージ B 2 2とが比較されて、 3次 元差異イメージ V B 2 1が生成される。
また、 コンピュータにより、 3次元差異イメージ V B 2 1の表面に対し、 前回 の比較よりも細かい精度にて 3次元イメージ B 2 1と 3次元イメージ B 2 2との 3次元差異細密イメージ V B 2 2が生成される。
次に、 コンピュータにより、 3次元イメージ B 2 1の 3次元差異細密イメージ と 3次元イメージ B 2 2の 3次元差異細密イメージとを比較して、 3次元差異ィ メージ V B 2 1の差異部の表示形式が決定される。 この決定の際には、 ポクセル の差異数に応じた色が決定される。
そして、 コンピュータにより、 表示形式を決定した 3次元差異イメージ V B 2 1が、 ポクセル比較結果として所定の表示形式でディスプレイ等に表示される。 これにより、 3次元イメージの比較において、 比較処理を高精度かつ高速に行 うことが可能となる。 また、 イメージの比較対象を 3次元形状のポクセルとする ことで、 無数の点群を比較する必要がなく、 メモリ資源の効率化及び演算回数の 低減を行うことが可能となる。 さらに、 3次元イメージの比較を表面のみとする ことで、 演算回数の大幅な低減が可能となる。
なお、 上記の説明では、 3次元イメージの表面ポクセル (粗ポクセル) に対し て、 3次元差異細密イメージを生成して比較する説明を行ったが、 3次元ィメ一 ジの表面ポクセルのうちの差異がある粗ポクセルのみ、 3次元差異細密イメージ を生成して比較することもできる。
[第 2の実施の形態]
次に、 第 2の実施の形態について説明する。 第 2の実施の形態は、 第 1の実施 の形態における比較処理の形態が異なる。 すなわち、 第 2の実施の形態は、 3次 元イメージをポクセル単位で比較し、 さらに 3次元イメージの差異については、 基準ポクセル各面からの差異ポクセル数の深さに対応した所定の表示形式により、 基準ポクセルの各面を表示するものである。
そこで、 第 2の実施の形態における 3次元イメージ比較プログラムの比較処理 について、 図 18および図 19を参照して具体的に説明する。
図 18は、 第 2の実施の形態において、 3次元イメージ比較プログラムの比較 処理の概要を示す図である。 第 2の実施の形態における比較処理は、 第 1の実施 の形態にて作成したポクセルを用い、 3次元イメージの一方におけるポクセルの 各面に対して、 差異を色にて表示する。 その色は、 以下の様に、 各面から垂直方 向 (奥行き) のポクセル数にて決定する。
図 18を参照して説明すると、 まず、 第 1の実施の形態にて作成したポクセル を用いて、 3次元イメージを比較する (ステップ S i b) 。 この比較された 3次 元差異イメージ VB 21の不一致部分を、 部分的に抽出する (ステップ S 2b) 。 次に、 ある一致している表面のポクセル (基準ポクセル BB 1) に対して、 3次 元形状の各面 (z y面、 z x面、 yx面) から垂直方向に不一致ポクセルのポク セル数 (差異ポクセル数の深さ) をカウントし、 表示形式を決定する (ステップ S 3b) 。 つまり、 A— A矢視断面に示す y z面と z X面とを例にすると、 yz 面では 3ポクセルの差異があり、 z X面では 2ポクセルの差異があると評価され る。 そして、 このカウント数に対応した所定の表示形式により、 一致している表 面のポクセル (基準ポクセル) の各面を表示する (ステップ S 4 b) 。 なお、 こ の表示形式は、 前述と同様の色 (RGBあるいは CMY) とすることができる。 このように本発明では、 コンピュータにより、 まず、 第 1の実施の形態にて作 成したポクセルを用いて、 3次元イメージが比較される。 また、 この比較された 3次元差異イメージ VB 21の不一致部分が、 部分的に抽出される。
次に、 コンピュータにより、 ある一致している表面のポクセル (基準ポクセル BB 1) に対して、 3次元形状の各面 (z y面、 z x面、 yx面) から垂直方向 に不一致ポクセルのポクセル数 (差異ポクセル数の深さ) がカウントされ、 表示 形式が決定される。
そして、 コンピュータにより、 このカウント数に対応した所定の表示形式をも とに、 一致している表面のポクセル (基準ポクセル) の各面が表示される。
上記した概要における機能は、 図 4の機能プロック図のうちのポクセル差異評 価処理部 40の機能が違うだけである。 したがって、 第 2の実施の形態では、 ポ クセル差異評価処理部 40の新たな機能 (重複含む) についてのみ説明する。 なお、 概要の説明では、 不一致部分 (差異部) に隣接し且つ一致している表面 のポクセルを基準ポクセルとして説明しているが、 以下の説明では、 差異部に隣 接し且つ一致している表面ポクセルを、 選択して範囲決定を行い、 この範囲を評 価精度にて再ポクセル化した際のポクセルを基準ポクセルとしている。
ポクセル差異評価処理部 4 0の新たな機能は、 ィメージ位置合わせ処理部 3 0 にて、 位置合わせされた 2つのイメージを比較するイメージ比較処理を行い、 比 較した 2つのイメージを所定の表示形式で表示する。 なお、 所定の表示形式とは、 たとえば、 R G Bのうちの各色の濃さをィメージの差異により変更して表示する ことである。 ここで、 ポクセル差異評価処理部 4 0は、 表面ポクセルを選択する。 この選択によりポクセル差異評価処理部 4 0は、 選択されたポクセルが表面ポク セルかどうか判断する。 そして、 ポクセル差異評価処理部 4 0は、 この判断によ り表面ポクセルが選択されれば、 再ポクセル化の範囲決定を行う。 すなわち、 ィ メージがー致しない部分を、 詳細に調査する範囲として決定する。
また、 ポクセル差異評価処理部 4 0は、 決定された範囲に対して、 評価精度で の再ポクセル化を行う。 すなわち、 評価精度に応じて、 範囲決定された表面ポク セルをさらに細かくポクセル化する。 そして、 ポクセル差異評価処理部 4 0は、 細密ポクセルに変換された 3次元サーフェスと計測表面/断面データとをもとに、 差異ポクセルを計数する。 すなわち、 これは、 一致している表面の細密ポクセル (基準ポクセル) に対して、 基準ポクセルの 3次元形状の各面 (z y面、 z x面、 y x面) から不一致ポクセルのポクセル数 (差異ポクセル数の深さ) をカウント する。
さらに、 ポクセル差異評価処理部 4 0は、 計数された差異ポクセルをもとに、 表示形式である色の濃さを決定し、 ポクセルの色付けを行う。 すなわち、 ポクセ ル差異評価処理部 4 0は、 一致している表面の細密ポクセル (基準ポクセル) に 対して、 基準ポクセルの各面それぞれの表示形式 (色) を決定する。
そして、 ポクセル差異評価処理部 4 0は、 ポクセルが色付けされたかどうかを 判定する。 すなわち、 ポクセル差異評価処理部 4 0は、 ポクセル化された 3次元 ィメージが比較されたかどうかを判定する。
次に、 第 2の実施の形態における 3次元イメージ比較プログラムの比較処理に ついて、 図 1 9を参照して具体的に説明する。 なお、 本発明の 3次元イメージ比 較処理は、 前述したように第 1〜第 3の実施の形態に対応して複数ある。 このた め、 ここでは、 第 2の実施の形態における 3次元ィメ一ジ比較処理を、 第 2のィ メージ比較処理として説明する。
図 1 9は、 図 1 3の 3次元イメージ比較処理における第 2のイメージ比較処理 を示すフローチャートである。 図 1 3のステップ S 4 0 (イメージ比較処理) が 実行されると、 以下の流れに従って処理が行われる。 なお、 この第 2のイメージ 比較処理は、 3次元イメージ比較装置のポクセル差異評価処理部 4 0により行わ れる。 なお、 概要の説明では、 不一致部分 (差異部) に隣接し且つ一致している 表面のポクセルを基準ポクセルとして説明しているが、 以下の説明では、 差異部 に隣接し且つ一致している表面ポクセルを、 選択して範囲決定を行い、 この範囲 を評価精度にて再ポクセル化した際のポクセルを基準ポクセルとしている。
[ステップ S 4 2 1 ] ポクセル差異評価処理部 4 0は、 表面ポクセルを選択 する。 なお、 表面ポクセルの選択は、 たとえば、 ポクセル化された 3次元ィメ一 ジに対して、 表面ポクセルの内側がオンであり、 外側がオフであることにより容 易に判断できる。
[ステップ S 4 2 2 ] ポクセル差異評価処理部 4 0は、 ステップ S 4 2 1に て選択されたポクセルが表面ポクセルかどうかを判断する。 ここで、 ポクセル差 異評価処理部 4 0は、 表面ポクセルが選択されていれば、 ステップ S 4 2 3へ進 み、 表面ポクセルが選択されていなければ、 ステップ S 4 2 1へ戻る。
[ステップ S 4 2 3 ] ポクセル差異評価処理部 4 0は、 ステップ S 4 2 2に て表面ポクセルが選択されたので、 再ポクセル化の範囲決定を行う。 すなわち、 イメージが一致しない部分を、 詳細に調査する範囲として決定する。
[ステップ S 4 2 4 ] ポクセル差異評価処理部 4 0は、 ステップ S 4 2 3に て決定された範囲に対して、 評価精度での再ポクセル化を行う。 すなわち、 評価 精度に応じて、 範囲決定された表面ポクセルをさらに細かくポクセル化する。
[ステップ S 4 2 5 ] ポクセル差異評価処理部 4 0は、 再度、 ポクセルに変 換された 3次元サ一フェスと計測表面 Z断面データとをもとに、 差異ポクセルを 計数する。 すなわち、 これは、 一致している表面のポクセル (基準ポクセル) に 対して、 基準ポクセルの 3次元形状の各面 (z y面、 z x面、 y x面) から不一 致ポクセルのポクセル数 (差異ポクセル数の深さ) をカウントする。
[ステップ S 4 2 6 ] ポクセル差異評価処理部 4 0は、 ステップ S 4 2 5に て計数された差異ポクセルをもとに、 表示形式である色の濃さを決定し、 ポクセ ルの色付けを行う。 すなわち、 ポクセル差異評価処理部 4 0は、 一致している表 面のポクセル (基準ポクセル) に対して、 ポクセルの各面それぞれの表示形式 (色) を決定する。
[ステップ S 4 2 7 ] ポクセル差異評価処理部 4 0は、 ステップ S 4 2 6に てポクセルが色付けされたかどうかを判定する。 すなわち、 ポクセル差異評価処 理部 4 0は、 ポクセル化された 3次元イメージが比較されたかどうかを判定する。 ここで、 ポクセル差異評価処理部 4 0は、 ポクセルが色付けされていなければ、 ステップ S 4 2 1へ戻り同様の処理を繰り返し、 ポクセルが色付けされていれば、 本処理を終了して図 1 3のステップ S 4 0へ復帰する。
このように、 第 2の実施の形態における比較処理は、 比較するための計算機資 源 (メモリ量、 ディスク領域、 計算時間) を軽減するため、 深さ方向を含むポク セルの替わりに、 3次元形状の各面の色にて差異を表示する。 これにより計算資 源は、 数分の 1となる。 なお、 差異を表現する色は、 奥行きのポクセル数に比例 する濃さにて決められる。 さらに、 この濃さは、 全体の濃さの最大値を基準に平 準化される。
なお、 上記説明では、 不一致部分 (差異部) に隣接し且つ一致している表面の ポクセルを基準ポクセルとしての説明と、 差異部に隣接し且つ一致している表面 ポクセルを、 選択して範囲決定を行い、 この範囲を評価精度にて再ポクセル化し た際のポクセルを基準ポクセルとしての説明とを行っているが、 単に、 一致して いる表面のポクセルを基準ポクセルとすることもできる。
[第 3の実施の形態]
次に、 第 3の実施の形態について説明する。 第 3の実施の形態は、 第 1の実施 の形態における比較処理の形態が異なる。 すなわち、 第 3の実施の形態は、 3次 元イメージをポクセル単位で比較し、 さらに 3次元イメージの差異部全体におけ る表面のポクセル数と、 その体積との比率に対応した所定の表示形式により、 ポ クセルを表示するものである。
そこで、 第 3の実施の形態における 3次元イメージ比較プログラムの比較処理 について、 図 2 0〜図 2 5を参照して具体的に説明する。
図 2 0は、 第 3の実施の形態において、 3次元イメージ比較プログラムの比較 処理の差異評価 (1 ) を示す図である。 第 1の実施の形態にて生成されたポクセ ルの差異部は、 2つのイメージが異なっている部分のみとなる。 この異なってい る部分 (島) は、 1つ、 あるいは, 複数のポクセルから構成されている。 そして、 第 3の実施の形態における比較処理は、 これらの島における表面のポクセル数と、 その島におけるポクセルの全数との比を色として表示する。 その色は、 以下の様 に決定する。 なお、 以下の説明では、 この島を単に差異部と表現する。
図 2 0を参照して説明すると、 第 3の実施の形態における比較処理では、 第 1 の実施の形態にて作成したポクセルを用いて、 3次元イメージであるイメージ A とイメージ Bとを比較する。 ここで、 差異部の形状は、 次のように表現される。 '体積: 1 2ポクセル
· Y方向の表面積: 2 5ポクセル面
この結果、 表面ポクセル数 Z全ポクセル数の比は 2 5 1 2となり、 ある差異 部に対するポクセルの差異が薄いことを示している。 また、 ある表面積当たりの 体積は、 図 2 1に示すような色で設定される。
図 2 1は、 第 3の実施の形態において、 3次元イメージ比較プログラムの比較 処理により決定される色 (薄い色) を示す図である。
図 2 1によると、 差異部を表示する色について、 R G Bを例にして説明してい る。 この例では、 青 (B ) の色を可変にすることで、 差異部全体のポクセル差異 を表現している。 図 2 1の例の場合、 ある表面に対するポクセルの差異が薄いこ と示すには、 たとえば、 設定位置 S T 1 1における R G B設定は、 2 5 5 (R) 、 2 5 5 (G) 、 2 2 3 (B ) と設定される。 この 2 2 3 (B) は、 薄い色を示し ており、 ある差異部全体のポクセル差異が薄いことを示している。 これが、 たと えば、 9 6 (B) であれば、 ある差異部全体におけるポクセル差異は、 濃いと示 すことができる。 以下、 このポクセル差異が濃い場合におけるポクセル差異評価 を図 2 2に示す。 図 2 2は、 第 3の実施の形態において、 3次元イメージ比較プログラムの比較 処理の差異評価 (2 ) を示す図である。 なお、 図 2 2は、 図 2 0に対してポクセ ル差異が濃い場合について説明するものである。
図 2 2を参照して説明すると、 第 3の実施の形態における比較処理では、 第 1 の実施の形態にて作成したポクセルを用いて、 3次元イメージであるイメージ A とイメージ Bとを比較する。 ここで、 差異部の形状は、 次のように表現される。 •体積: 4 3ポクセル
• Y方向の表面積: 3 2ポクセル面
この結果、 表面ポクセル数 Z全ポクセル数の比は 3 2 Z 4 3となり、 ある差異 部に対するポクセルの差異が濃いことを示している。 また、 ある表面積当たりの 体積は、 図 2 3に示すような色で設定される。
図 2 3は、 第 3の実施の形態において、 3次元イメージ比較プログラムの比較 処理により決定される色 (濃い色) を示す図である。
図 2 3によると、 差異部を表示する色について、 R G Bを例にして説明してい る。 この例では、 青 (B) の色を可変にすることで、 差異部全体のポクセル差異 を表現している。 図 2 3の例の場合、 ある表面に対するポクセルの差異が濃いこ と示すには、 たとえば、 設定位置 S T 1 2における R G B設定は、 2 5 5 (R) 、 2 5 5 (G) 、 9 6 (B) と設定される。 この 9 6 (B ) は、 濃い色を示してお り、 ある差異部全体のポクセル差異が濃いことを示している。
また、 差異が発生している部分を精細なポクセルにて、 再度差異を評価し、 差 異が連続的に発生している個所を同様に色付けする。 このとき、 色を設定する部 分を、 比較する 3次元サーフェスの部分サーフェスとして記憶する。 この差異評 価の例を図 2 4に図示する。
図 2 4は、 第 3の実施の形態において、 3次元イメージ比較プログラムの比較 処理の差異評価 (3 ) を示す図である。 ここで、 図 2 4の分離 (太線) ラインは、 元の差異ポクセルから、 精細なポクセルにて比較された差異ポクセルの分離位置 を示す。 なお、 図 2 4は、 図 2 0に対して、 ポクセル差異を更に細かく差異評価 した場合について示すものである。 また、 差異評価の方法の説明については、 更 に細かいポクセルに対して図 2 0と図 2 2と同様の処理を行うため省略する。 上記した概要における機能は、 図 4の機能ブロック図のうちのポクセル差異評 価処理部 4 0の機能が違うだけである。 したがって、 第 3の実施の形態では、 ポ クセル差異評価処理部 4 0の新たな機能 (重複含む) についてのみ説明する。 な お、 以下の説明では、 不一致部分 (差異部) に隣接し且つ一致している表面ポク セルを、 選択して範囲決定を行い、 評価精度にて再ポクセル化した際のポクセル を基準ポクセルとしている。 もちろん、 第 3の実施の形態では、 不一致部分 (差 異部) に隣接し且つ一致している表面のポクセルを基準ポクセルとしてもよい。 また、 単に、 一致している表面のポクセルを基準ポクセルとしてもよい。
ポクセル差異評価処理部 4 0の新たな機能は、 イメージ位置合わせ処理部 3 0 にて位置合わせされた 2つのイメージを比較して評価し、 評価した 2つのィメー ジを所定の表示形式で表示する。 なお、 所定の表示形式とは、 たとえば、 R G B のうちの各色の濃さをイメージの差異により変更して表示することである。 ここ で、 ポクセル差異評価処理部 4 0は、 表面ポクセルを選択する。 この選択により ポクセル差異評価処理部 4 0は、 選択されたポクセルが表面ポクセルかどうかを 判断する。 そして、 ポクセル差異評価処理部 4 0は、 表面ポクセルが選択された ので、 再ポクセル化の範囲決定を行う。 すなわち、 イメージが一致しない部分を、 詳細に調査する範囲として決定する。
また、 ポクセル差異評価処理部 4 0は、 決定された範囲に対して、 評価精度で の再ポクセル化を行う。 すなわち、 評価精度に応じて、 範囲決定された表面ポク セルをさらに細かくポクセル化する。 そして、 ポクセル差異評価処理部 4 0は、 ポクセルに変換された 3次元サーフェスと計測表面 Z断面データとをもとに、 差 異ポクセルを計数する。 すなわち、 3次元差異イメージの差異部全体における表 面ポクセル数と全体ポクセル数との比率を求める。
さらに、 ポクセル差異評価処理部 4 0は、 求められた差異ポクセルの比率をも とに、 表示形式である色の濃さを決定し、 ポクセルの色付けを行う。 すなわち、 ポクセル差異評価処理部 4 0は、 不一致ポクセルからなる差異部全体の表示形式 (色) を決定する。
そして、 ポクセル差異評価処理部 4 0は、 ポクセルが色付けされたかどうかを 判定する。 すなわち、 ポクセル差異評価処理部 4 0は、 ポクセル化された 3次元 ィメージが比較されたかどうかを判定する。
次に、 第 3の実施の形態における 3次元イメージ比較プログラムの比較処理に ついて、 図 2 5を参照して具体的に説明する。 なお、 本発明の 3次元イメージ比 較処理は、 前述したように第 1〜第 3の実施の形態に対応して複数ある。 このた め、 ここでは、 第 3の実施の形態における 3次元イメージ比較処理を、 第 3のィ メージ比較処理として説明する。
図 2 5は、 図 1 3の 3次元イメージ比較処理における第 3のイメージ比較処理 を示すフローチャートである。 図 1 3のステップ S 4 0 (イメージ比較処理) が 実行されると、 以下の流れに従って処理が行われる。 なお、 この第 3のイメージ 比較処理は、 3次元イメージ比較装置のポクセル差異評価処理部 4 0により行わ れる。
[ステップ S 4 3 1 ] ポクセル差異評価処理部 4 0は、 表面ポクセルを選択 する。 なお、 表面ポクセルの選択は、 たとえば、 ポクセル化された 3次元ィメー ジに対して、 表面ポクセルの内側がオンであり、 外側がオフであることにより容 易に判断できる。
[ステップ S 4 3 2 ] ポクセル差異評価処理部 4 0は、 ステップ S 4 3 1に て選択されたポクセルが表面ポクセルかどうか判断する。 ここで、 ポクセル差異 評価処理部 4 0は、 表面ポクセルが選択されていれば、 ステップ S 4 3 3へ進み、 表面ポクセルが選択されていなければ、 ステップ S 4 3 1へ戻る。
[ステップ S 4 3 3 ] ポクセル差異評価処理部 4 0は、 ステップ S 4 3 2に て表面ポクセルが選択されたので、 再ポクセル化の範囲決定を行う。 すなわち、 イメージが一致しない部分を、 詳細に調査する範囲として決定する。
[ステップ S 4 3 4 ] ポクセル差異評価処理部 4 0は、 ステップ S 4 3 3に て決定された範囲に対して、 評価精度での再ポクセル化を行う。 すなわち、 評価 精度に応じて、 範囲決定された表面ポクセルをさらに細かくポクセル化する。
[ステップ S 4 3 5 ] ポクセル差異評価処理部 4 0は、 ポクセルに変換され た 3次元サ一フェスと計測表面/断面データとをもとに、 差異ポクセルを計数す る。 すなわち、 3次元差異イメージの差異部全体における表面ポクセル数と全体 ポクセル数との比率を求める。 [ステップ S 4 3 6 ] ポクセル差異評価処理部 4 0は、 ステップ S 4 3 5に て求められた差異ポクセルの比率をもとに、 表示形式である色の濃さを決定し、 ポクセルの色付けを行う。 すなわち、 ポクセル差異評価処理部 4 0は、 不一致ポ クセルからなる差異部全体の表示形式 (色) を決定する。
[ステップ S 4 3 7 ] ポクセル差異評価処理部 4 0は、 ステップ S 4 3 6に てポクセルが色付けされたかどうかを判定する。 すなわち、 ポクセル差異評価処 理部 4 0は、 ポクセル化された 3次元イメージが比較されたかどうかを判定する。 ここで、 ポクセル差異評価処理部 4 0は、 ポクセルが色付けされていなければ、 ステップ S 4 3 1へ戻り同様の処理を繰り返し、 ポクセルが色付けされていれば、 本処理を終了して図 1 3のステップ S 4 0へ復帰する。
このように、 第 3の実施の形態における比較処理は、 比較するための計算機資 源 (メモリ量、 ディスク領域、 計算時間) を軽減するため、 深さ方向を含むポク セルの替わりに、 差異部全体の色にて差異を表示する。 これにより計算資源は、 数分の 1となる。 なお、 表面ポクセル数/全ポクセル数の比は、 最大 2, 最小 0 (正確には、 0に近い微小値) となる。 従って、 色の濃さは、 2を最も薄い色と し、 0にて最も濃い色となるよう様に設定する。 また、 差異部の大きさは、 3次 元表示した際の差異部全体の大きさで理解可能である。
[第 1の実施の形態における計算資源の評価結果例]
最後に、 第 1の実施の形態における計算資源の評価結果例について説明する。 ここで、 第 1の実施の形態における効果は、 以下の 3点となる。
•従来方式に比べて、 圧倒的に必要とする計算量が少なく、 結果を得る時間が短 い。
-詳細な比較を現実的な計算機資源 (メモリ、 ディスク量) にて実現可能となる。 •比較作業の効率化におけるメリットがある。
例えば、 計算資源について、 以下のように評価する。
まず、 従来技術における表 (1 ) にある計測点数を、 式 (1 ) に示す 3次元サ 一フェスの方程式を用いて最短距離を計算すると、 以下のような計算量となる。
•サーフェス上の最近節点の取得 数百 floating point operatio ns • この点から当該の計測点までの距離の計算 ' · '数十 floating point operatio ns
'近接すると思われるサ一フェスの数 1〜数十
したがって、 この計算は、 1つの最短距離を求めるために、 数万 floating poi nt operations の計算が必要となる。 これを 1 0 0億点に対して繰り返す。 結局、 この計算にかかる時間は、 1 G H zの C P Uにて 2 floating point operations が同時に行えるとして、 3 0時間程となる。
一方、 本発明の実施の形態による方法は、 複数のサーフェスを含む体積に対し て、 全空間を計測精度のポクセルに分割し、 ポクセル毎に論理演算を施すだけで あるため、 計算は 5 0 0億 floating point operationsである。 したがって、 同 じ計算機では、 計算にかかる時間が 2 5秒となる。 これにより、 本発明における 第 1の実施の形態による方法は、 従来方法に比べ計算機の処理時間を大幅に短縮 することが可能となった。 また、 本発明における第 1の実施の形態による方法は、 全空間を計測精度のポクセルに分割するので、 従来方法と同等の精度により比較 計算が可能となった。
以上説明した処理は、 コンピュータプログラムに記述し、 コンピュータで実行 することにより、 本発明の機能を実現することができる。 また、 コンピュータで 実行する際には、 コンピュータ内のハードディスク等にコンピュータプログラム を予め格納しておき、 メインメモリにロードして実行する。 なお、 コンピュータ プログラムは、 コンピュータで読み取り可能な記録媒体に記録しておくことがで きる。 コンピュータで読み取り可能な記録媒体としては、 磁気記録媒体、 光ディ スク、 光磁気記録媒体、 半導体メモリなどがある。 磁気記録媒体には、 ハードデ イスク、 フレキシブルディスク (FD) 、 Z I P (Zip:磁気ディスクの一種) 、 磁気テープなどがある。 光ディスクには、 D VD (Digital Versatile Disc) 、 D VD - R AM (DVD Random Access Memory) 、 C D - R OM (Compact Disc Read Only Memory) 、 C D— R (CD Recordable) 、 C D - RW (C D Rewritable) などがある。 光磁気記録媒体には、 M〇 (Magneto Optical Di sk) などがある。 半導体メモリには、 フラッシュ ·メモリ (Flash Memory) な どがある。 また、 コンピュータプログラムを流通させる場合には、 例えば各コンピュータ プログラムが記録された D V D、 C D— R OMなどの可搬型記録媒体が販売され る。 また、 コンピュータプログラムをサーバの記憶装置に格納しておき、 ネット ワークを介して、 サーバからクライアントへコンビュ一夕プログラムを転送する こともできる。
以上説明したように本発明の第 1の態様では、 ポクセル化された 2つの 3次元 イメージを比較し、 この比較による表面を更に細かい精度で生成して比較し、 こ の比較結果をもとに表示形式を決定して、 その表示形式で表示するようにしたの で、 3次元イメージの比較を高精度かつ高速に行うことが可能となる。
また、 本発明の第 2の態様では、 3次元イメージを比較した比較結果をもとに、 基準ポクセルの各面に対して、 基準ポクセル各面から垂直方向への差異ポクセル 数に応じた表示形式を決定して、 その表示形式で表示するようにしたので、 3次 元イメージの比較を高精度かつ高速に行うことが可能となる。
さらに、 本発明の第 3の態様では、 3次元イメージを比較した比較結果をもと に、 3次元ィメージの差異部全体における表面積を体積で除算した比率に応じた 表示形式を決定して、 その表示形式で表示するようにしたので、 3次元イメージ の比較を高精度かつ高速に行うことが可能となる。
上記については単に本発明の原理を示すものである。 さらに、 多数の変形、 変 更が当業者にとって可能であり、 本発明は上記に示し、 説明した正確な構成およ び応用例に限定されるものではなく、 対応するすべての変形例および均等物は、 添付の請求項およびその均等物による本発明の範囲とみなされる。

Claims

請 求 の 範 囲
1 . 3次元イメージを比較して差異を表示する 3次元イメージ比較プログラム において、
コンピュータに、
予めポクセル化された第 1の 3次元ィメ一ジと第 2の 3次元イメージとを比較 して、 3次元差異イメージを生成し、
前記 3次元差異イメージの表面に対し、 前記比較よりも細かい精度にて前記第 1の 3次元ィメージと前記第 2の 3次元ィメージとの 3次元差異細密ィメ一ジを 生成し、
前記第 1の 3次元イメージの 3次元差異細密ィメージと、 前記第 2の 3次元ィ メ一ジの 3次元差異細密イメージとを比較して、 前記 3次元差異イメージの差異 部の表示形式を決定し、
前記 3次元差異イメージをポクセル比較結果として前記表示形式で表示する、 処理を実行させることを特徴とする 3次元イメージ比較プログラム。
2 . 前記表示形式は、 色表示であることを特徴とする請求の範囲第 1項記載の 3次元イメージ比較プログラム。
3 . 前記色表示は、 赤、 緑、 および青の光 3原色を基本とした加法混色により 行われることを特徴とする請求の範囲第 2項記載の 3次元ィメージ比較プロダラ ム。
4 . 前記色表示は、 シアン、 マゼンダ、 および黄の補色 3原色を基本とした減 法混色により行われることを特徴とする請求の範囲第 2項記載の 3次元イメージ 比較プログラム。
5 . 前記減法混色には、 シアン、 マゼンダ、 および黄の補色 3原色の他に、 黒 も追加することを特徴とする請求の範囲第 4項記載の 3次元イメージ比較プログ ラム。
6 . 前記比較は、 第 1の 3次元イメージと第 2の 3次元イメージとの間におい て、 ポクセル情報のブーリアン演算により行われることを特徴とする請求の範囲 第 1項記載の 3次元イメージ比較プログラム。
7 . 前記第 1の 3次元イメージの 3次元差異細密イメージと、 前記第 2の 3次 元イメージの 3次元差異細密イメージとの比較は、 2つの 3次元差異細密ィメー ジの差をとることを特徴とする請求の範囲第 1項記載の 3次元イメージ比較プロ グラム。
8 . C ADデータからなる 3次元サ一フェス情報をポクセル化して、 前記第 1 の 3次元ィメージを作成することを特徴とする請求の範囲第 1項記載の 3次元ィ メ一ジ比較プログラム。
9 . 計測データからなる点群情報をポクセル化して、 前記第 2の 3次元ィメー ジを作成することを特徴とする請求の範囲第 1項記載の 3次元イメージ比較プロ クラム。
1 0 . 3次元イメージを比較して差異を表示する 3次元イメージ比較プロダラ ムにおいて、
コンピュータに、
予めポクセル化された第 1の 3次元イメージと第 2の 3次元イメージとを比較 して、 3次元差異イメージを生成し、
前記 3次元差異イメージの差異部を部分的に抽出し、 一致している表面のポク セルに対して、 3次元形状の各面から垂直方向に不一致ポクセルのポクセル数を カウントし、 前記力ゥントの数をもとにポクセルの差異を表す表示形式を決定し、 前記表示形式を決定した前記 3次元差異イメージを、 前記表示形式により表示 する、
処理を実行させることを特徴とする 3次元イメージ比較プログラム。
1 1 . 前記表示形式は、 色表示であることを特徴とする請求の範囲第 1 0項記 載の 3次元イメージ比較プログラム。
1 2 . 前記色表示は、 赤、 緑、 および青の光 3原色を基本とした加法混色によ り行われることを特徴とする請求の範囲第 1 1項記載の 3次元イメージ比較プロ グラム。
1 3 . 前記色表示は、 シアン、 マゼンダ、 および黄の補色 3原色を基本とした 減法混色により行われることを特徴とする請求の範囲第 1 1項記載の 3次元ィメ —ジ比較プログラム。
1 4. 前記減法混色には、 シアン、 マゼンダ、 および黄の補色 3原色の他に、 黒も追加することを特徴とする請求の範囲第 1 3項記載の 3次元イメージ比較プ ログラム。
1 5 . 前記表示は、 3次元形状の各面それぞれの前記カウントの数をもとに、 前記 3次元差異イメージのポクセル各面に対してそれぞれ行うことを特徴とする 請求の範囲第 1 0項記載の 3次元イメージ比較プログラム。
1 6 . C ADデータからなる 3次元サーフェス情報をポクセル化して、 前記第 1の 3次元ィメ一ジを作成することを特徴とする請求の範囲第 1 0項記載の 3次 元ィメ一ジ比較プログラム。
1 7 . 計測データからなる点群情報をポクセル化して、 前記第 2の 3次元ィメ ージを作成することを特徴とする請求の範囲第 1 0項記載の 3次元イメージ比較 プログラム。
1 8 . 3次元イメージを比較して差異を表示する 3次元イメージ比較プロダラ ムにおいて、
コンピュータに、
予めポクセル化された第 1の 3次元ィメージと第 2の 3次元イメージとを比較 して、 3次元差異イメージを生成し、
前記 3次元差異イメージの差異部における表面ポクセル数と全体ポクセル数と の比率を算出し、 前記 3次元差異イメージに対し、 前記比率をもとにポクセルの 差異を表す表示形式を決定し、
前記表示形式を決定した前記 3次元差異ィメージに対して、 前記比率に対応し た所定の表示形式により前記差異部を表示する、
処理を実行させることを特徴とする 3次元イメージ比較プログラム。
1 9 . 前記表示形式は、 色表示であることを特徴とする請求の範囲第 1 8項記 載の 3次元イメージ比較プログラム。
2 0 . 前記色表示は、 赤、 緑、 および青の光 3原色を基本とした加法混色によ り行われることを特徴とする請求の範囲第 1 9項記載の 3次元イメージ比較プロ グラム。
2 1 . 前記色表示は、 シアン、 マゼンダ、 および黄の補色 3原色を基本とした 減法混色により行われることを特徴とする請求の範囲第 1 9項記載の 3次元ィメ —ジ比較プログラム。
2 2 . 前記減法混色には、 シアン、 マゼンダ、 および黄の補色 3原色の他に、 黒も追加することを特徴とする請求の範囲第 2 1項記載の 3次元イメージ比較プ ログラム。
2 3 . C ADデータからなる 3次元サ一フェス情報をポクセルィヒして、 前記第 1の 3次元ィメージを作成することを特徴とする請求の範囲第 1 8項記載の 3次 元イメージ比較プログラム。
2 4 . 計測データからなる点群情報をポクセル化して、 前記第 2の 3次元ィメ ージを作成することを特徴とする請求の範囲第 1 8項記載の 3次元イメージ比較 プログラム。
2 5 . 3次元イメージを比較して差異を表示する 3次元イメージ比較方法にお いて、
予めポクセル化された第 1の 3次元イメージと第 2の 3次元イメージとを比較 して、 3次元差異イメージを生成し、
前記 3次元差異イメージの表面に対し、 前記比較よりも細かい精度にて前記第 1の 3次元イメージと前記第 2の 3次元イメージとの 3次元差異細密イメージを 生成し、
前記第 1の 3次元イメージの 3次元差異細密イメージと、 前記第 2の 3次元ィ メージの 3次元差異細密イメージとを比較して、 前記 3次元差異イメージの差異 部の表示形式を決定し、
前記 3次元差異イメージをポクセル比較結果として前記表示形式で表示する、 ことを特徴とする 3次元イメージ比較方法。
2 6 . 3次元イメージを比較して差異を表示する 3次元イメージ比較方法にお いて、
予めポクセルィ匕された第 1の 3次元イメージと第 2の 3次元イメージとを比較 して、 3次元差異イメージを生成し、
前記 3次元差異イメージの差異部を部分的に抽出し、 一致している表面のポク セルに対して、 3次元形状の各面から垂直方向に不一致ポクセルのポクセル数を カウントし、 前記カウントの数をもとにポクセルの差異を表す表示形式を決定し、 前記表示形式を決定した前記 3次元差異イメージを、 前記表示形式により表示 する、
ことを特徴とする 3次元ィメージ比較方法。
2 7 . 3次元イメージを比較して差異を表示する 3次元イメージ比較方法にお いて、
予めポクセルィ匕された第 1の 3次元イメージと第 2の 3次元イメージとを比較 して、 3次元差異イメージを生成し、
前記 3次元差異イメージの差異部における表面ポクセル数と全体ポクセル数と の比率を算出し、 前記 3次元差異イメージに対し、 前記比率をもとにポクセルの 差異を表す表示形式を決定し、
前記表示形式を決定した前記 3次元差異ィメージに対して、 前記比率に対応し た所定の表示形式により前記差異部を表示する、
ことを特徴とする 3次元イメージ比較方法。
2 8 . 3次元イメージを比較して差異を表示する 3次元イメージ比較装置にお いて、
予めポクセル化された第 1の 3次元イメージと第 2の 3次元イメージとを比較 して、 3次元差異イメージを生成する差異イメージ生成部と、
前記 3次元差異イメージの表面に対し、 前記比較よりも細かい精度にて前記第 1の 3次元ィメージと前記第 2の 3次元イメージとの 3次元差異細密イメージを 生成する差異細密イメージ生成部と、
前記第 1の 3次元イメージの 3次元差異細密ィメージと、 前記第 2の 3次元 イメージの 3次元差異細密イメージとを比較して、 前記 3次元差異ィメージの差 異部の表示形式を決定する表示形式決定部と、
前記 3次元差異ィメ一ジをポクセル比較結果として前記表示形式で表示する差 異表示部と、
を有することを特徴とする 3次元ィメージ比較装置。
2 9 . 3次元イメージを比較して差異を表示する 3次元イメージ比較装置にお いて、 予めポクセル化された第 1の 3次元ィメージと第 2の 3次元イメージとを比較 して、 3次元差異イメージを生成する差異イメージ生成部と、
前記 3次元差異イメージの差異部を部分的に抽出し、 一致している表面のポク セルに対して、 3次元形状の各面から垂直方向に不一致ポクセルのポクセル数を カウントし、 前記カウントの数をもとにポクセルの差異を表す表示形式を決定す る各面差異決定部と、
前記表示形式を決定した前記 3次元差異ィメ一ジを、 前記表示形式により表示 する各面差異表示部と、
を有することを特徴とする 3次元ィメ一ジ比較装置。
3 0 . 3次元イメージを比較して差異を表示する 3次元イメージ比較装置にお いて、
予めポクセル化された第 1の 3次元イメージと第 2の 3次元イメージとを比較 して、 3次元差異イメージを生成する差異イメージ生成部と、
前記 3次元差異イメージの差異部における表面ポクセル数と全体ポクセル数と の比率を算出し、 前記 3次元差異イメージに対し、 前記比率をもとにポクセルの 差異を表す表示形式を決定する差異比率決定部と、
前記表示形式を決定した前記 3次元差異イメージに対して、 前記比率に対応し た所定の表示形式により前記差異部を表示する差異比率表示部と、
を有することを特徴とする 3次元イメージ比較装置。
3 1 . 3次元イメージを比較して差異を表示するためのプログラムを記録した コンピュータ読み取り可能な記録媒体において、
前記コンピュータに、
予めポクセル化された第 1の 3次元イメージと第 2の 3次元イメージとを比較 して、 3次元差異イメージを生成し、
前記 3次元差異イメージの表面に対し、 前記比較よりも細かい精度にて前記第 1の 3次元ィメージと前記第 2の 3次元ィメージとの 3次元差異細密ィメ一ジを 生成し、
前記第 1の 3次元ィメ一ジの 3次元差異細密イメージと、 前記第 2の 3次元ィ メージの 3次元差異細密イメージとを比較して、 前記 3次元差異ィメージの差異 部の表示形式を決定し、
前記 3次元差異イメージをポクセル比較結果として前記表示形式で表示する、 処理を実行させることを特徴とするプログラムを記録したコンピュータ読み取 り可能な記録媒体。
3 2 . 3次元イメージを比較して差異を表示するためのプログラムを記録した コンピュー夕読み取り可能な記録媒体において、
前記コンピュータに、
予めポクセル化された第 1の 3次元ィメ一ジと第 2の 3次元イメージとを比較 して、 3次元差異イメージを生成し、
前記 3次元差異イメージの差異部を部分的に抽出し、 一致している表面のポク セルに対して、 3次元形状の各面から垂直方向に不一致ポクセルのポクセル数を カウントし、 前記力ゥントの数をもとにポクセルの差異を表す表示形式を決定し、 前記表示形式を決定した前記 3次元差異ィメ一ジを、 前記表示形式により表示 する、
処理を実行させることを特徴とするプログラムを記録したコンピュータ読み取 り可能な記録媒体。
3 3 . 3次元イメージを比較して差異を表示するためのプログラムを記録した コンピュータ読み取り可能な記録媒体において、
前記コンピュータに、
予めポクセル化された第 1の 3次元イメージと第 2の 3次元イメージとを比較 して、 3次元差異イメージを生成し、
前記 3次元差異イメージの差異部における表面ポクセル数と全体ポクセル数と の比率を算出し、 前記 3次元差異イメージに対し、 前記比率をもとにポクセルの 差異を表す表示形式を決定し、
前記表示形式を決定した前記 3次元差異イメージに対して、 前記比率に対応し た所定の表示形式により前記差異部を表示する、
処理を実行させることを特徵とするプログラムを記録したコンピュータ読み取 り可能な記録媒体。
PCT/JP2002/006621 2002-06-28 2002-06-28 3次元イメージ比較プログラム、3次元イメージ比較方法、および3次元イメージ比較装置 WO2004003850A1 (ja)

Priority Applications (6)

Application Number Priority Date Filing Date Title
PCT/JP2002/006621 WO2004003850A1 (ja) 2002-06-28 2002-06-28 3次元イメージ比較プログラム、3次元イメージ比較方法、および3次元イメージ比較装置
KR1020047021408A KR100639139B1 (ko) 2002-06-28 2002-06-28 3차원 이미지 비교 프로그램을 기록한 컴퓨터로 판독 가능한 기록 매체, 3차원 이미지 비교 방법 및 3차원 이미지 비교 장치
EP02736203A EP1519318A4 (en) 2002-06-28 2002-06-28 PROGRAM AND METHOD FOR COMPARING THREE-DIMENSIONAL IMAGE AND APPARATUS THEREFOR
JP2004517224A JP4001600B2 (ja) 2002-06-28 2002-06-28 3次元イメージ比較プログラム、3次元イメージ比較方法、および3次元イメージ比較装置
CNB028292359A CN100388317C (zh) 2002-06-28 2002-06-28 三维图象的比较程序、比较方法及比较装置
US10/989,464 US20050068317A1 (en) 2002-06-28 2004-11-17 Program, method, and device for comparing three-dimensional images in voxel form

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2002/006621 WO2004003850A1 (ja) 2002-06-28 2002-06-28 3次元イメージ比較プログラム、3次元イメージ比較方法、および3次元イメージ比較装置

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US10/989,464 Continuation US20050068317A1 (en) 2002-06-28 2004-11-17 Program, method, and device for comparing three-dimensional images in voxel form

Publications (1)

Publication Number Publication Date
WO2004003850A1 true WO2004003850A1 (ja) 2004-01-08

Family

ID=29808175

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2002/006621 WO2004003850A1 (ja) 2002-06-28 2002-06-28 3次元イメージ比較プログラム、3次元イメージ比較方法、および3次元イメージ比較装置

Country Status (6)

Country Link
US (1) US20050068317A1 (ja)
EP (1) EP1519318A4 (ja)
JP (1) JP4001600B2 (ja)
KR (1) KR100639139B1 (ja)
CN (1) CN100388317C (ja)
WO (1) WO2004003850A1 (ja)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100419794C (zh) * 2005-12-16 2008-09-17 北京中星微电子有限公司 一种图像问题的检测方法
JP2009104515A (ja) * 2007-10-25 2009-05-14 Fujitsu Ltd 差異強調プログラム、差異強調処理方法及び差異強調処理装置
JP2012103144A (ja) * 2010-11-11 2012-05-31 Univ Of Tokyo 3次元環境復元装置、3次元環境復元方法、及びロボット
JP2016103263A (ja) * 2015-08-31 2016-06-02 株式会社大林組 構造物状況把握支援装置、構造物状況把握支援方法及びプログラム
JP2018005862A (ja) * 2016-07-08 2018-01-11 富士通株式会社 ファセット化処理プログラム、ファセット抽出プログラム、ファセット化処理方法、ファセット抽出方法および情報処理装置
JP2018008403A (ja) * 2016-07-12 2018-01-18 学校法人慶應義塾 立体物製造装置、立体物製造方法及びプログラム
JP2018181056A (ja) * 2017-04-17 2018-11-15 富士通株式会社 差分検知プログラム、差分検知装置、差分検知方法
JP2021021672A (ja) * 2019-07-30 2021-02-18 日本電気通信システム株式会社 距離計測装置、システム、方法、及びプログラム

Families Citing this family (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1669056A (zh) * 2002-06-19 2005-09-14 西门子公司 3d数据组的平台交叉和针对数据的可视化
US7685024B2 (en) * 2005-02-24 2010-03-23 Dolphin Software Ltd. System and method for computerized ordering
CN101331380B (zh) * 2005-12-16 2011-08-03 株式会社Ihi 三维形状数据的存储/显示方法和装置以及三维形状的计测方法和装置
DE112006003363B4 (de) * 2005-12-16 2016-05-04 Ihi Corporation Verfahren und Vorrichtung zur Identifizierung der Eigenposition, und Verfahren und Vorrichtung zur Messung einer dreidimensionalen Gestalt
CN101331381B (zh) * 2005-12-16 2011-08-24 株式会社Ihi 三维形状数据的位置对准方法和装置
TWI409717B (zh) * 2009-06-22 2013-09-21 Chunghwa Picture Tubes Ltd 適用於電腦產品與影像顯示裝置的影像轉換方法
CN107481217A (zh) * 2009-06-24 2017-12-15 皇家飞利浦电子股份有限公司 对象内的植入设备的空间和形状表征
US20110307281A1 (en) * 2010-06-11 2011-12-15 Satterfield & Pontikes Construction, Inc. Model inventory manager
US8768029B2 (en) 2010-10-20 2014-07-01 Medtronic Navigation, Inc. Selected image acquisition technique to optimize patient model construction
US10373375B2 (en) * 2011-04-08 2019-08-06 Koninklijke Philips N.V. Image processing system and method using device rotation
GB201308866D0 (en) * 2013-05-16 2013-07-03 Siemens Medical Solutions System and methods for efficient assessment of lesion developemnt
US9715761B2 (en) 2013-07-08 2017-07-25 Vangogh Imaging, Inc. Real-time 3D computer vision processing engine for object recognition, reconstruction, and analysis
US9984498B2 (en) * 2013-07-17 2018-05-29 Microsoft Technology Licensing, Llc Sparse GPU voxelization for 3D surface reconstruction
US9710960B2 (en) 2014-12-04 2017-07-18 Vangogh Imaging, Inc. Closed-form 3D model generation of non-rigid complex objects from incomplete and noisy scans
US10007996B2 (en) * 2015-03-02 2018-06-26 Lawrence Livermore National Security, Llc System for detecting objects in streaming 3D images formed from data acquired with a medium penetrating sensor
AU2016250935B2 (en) * 2015-04-20 2020-12-03 Mars Bioimaging Limited Improving material identification using multi-energy CT image data
KR101644426B1 (ko) * 2015-09-30 2016-08-02 상명대학교서울산학협력단 변형3d모델에 대응되는 원본3d모델 인식방법
WO2017200527A1 (en) * 2016-05-16 2017-11-23 Hewlett-Packard Development Company, L.P. Generating a shape profile for a 3d object
US10380762B2 (en) 2016-10-07 2019-08-13 Vangogh Imaging, Inc. Real-time remote collaboration and virtual presence using simultaneous localization and mapping to construct a 3D model and update a scene based on sparse data
WO2018075628A1 (en) 2016-10-19 2018-04-26 Shapeways, Inc. Systems and methods for identifying three-dimensional printed objects
EP3355279A1 (en) * 2017-01-30 2018-08-01 3D Repo Ltd Method and computer programs for identifying differences between 3-dimensional scenes
US10839585B2 (en) 2018-01-05 2020-11-17 Vangogh Imaging, Inc. 4D hologram: real-time remote avatar creation and animation control
US11080540B2 (en) 2018-03-20 2021-08-03 Vangogh Imaging, Inc. 3D vision processing using an IP block
US10810783B2 (en) 2018-04-03 2020-10-20 Vangogh Imaging, Inc. Dynamic real-time texture alignment for 3D models
JP7084193B2 (ja) * 2018-04-10 2022-06-14 ザイオソフト株式会社 医用画像処理装置、医用画像処理方法、及び医用画像処理プログラム
US11170224B2 (en) 2018-05-25 2021-11-09 Vangogh Imaging, Inc. Keyframe-based object scanning and tracking
US10776949B2 (en) 2018-10-30 2020-09-15 Liberty Reach Inc. Machine vision-based method and system for measuring 3D pose of a part or subassembly of parts
US11232633B2 (en) 2019-05-06 2022-01-25 Vangogh Imaging, Inc. 3D object capture and object reconstruction using edge cloud computing resources
US11170552B2 (en) 2019-05-06 2021-11-09 Vangogh Imaging, Inc. Remote visualization of three-dimensional (3D) animation with synchronized voice in real-time
JP7376268B2 (ja) * 2019-07-22 2023-11-08 ファナック株式会社 三次元データ生成装置及びロボット制御システム
US11335063B2 (en) 2020-01-03 2022-05-17 Vangogh Imaging, Inc. Multiple maps for 3D object scanning and reconstruction
KR102238911B1 (ko) * 2020-04-23 2021-04-16 비아이밀리그램㈜ 건축물의 공정진행상태 산출 장치, 건축물의 공정진행상태 산출 방법 및 건축물의 공정진행상태 산출시스템

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5273040A (en) * 1991-11-14 1993-12-28 Picker International, Inc. Measurement of vetricle volumes with cardiac MRI
US5301271A (en) * 1991-07-11 1994-04-05 Matsushita Electric Industrial Co., Ltd. Image processing separately processing two-density-level image data and multi-level image data
JPH0935058A (ja) * 1995-07-17 1997-02-07 Nec Corp 画像認識方法
JPH09204532A (ja) * 1996-01-25 1997-08-05 Hitachi Ltd 画像認識方法および画像表示方法
EP1054384A2 (en) * 1999-05-20 2000-11-22 Mitsubishi Denki Kabushiki Kaisha Method and apparatus for translating and interfacing voxel memory addresses
US6175371B1 (en) * 1995-06-02 2001-01-16 Philippe Schoulz Process for transforming images into stereoscopic images, images and image series obtained by this process

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62219075A (ja) * 1986-03-20 1987-09-26 Hitachi Medical Corp 三次元画像の半透明表示方法
JPH04276513A (ja) * 1991-03-04 1992-10-01 Matsushita Electric Ind Co Ltd 形状測定処理方法
JPH0516435A (ja) * 1991-07-11 1993-01-26 Matsushita Electric Ind Co Ltd 画像形成装置
AU6341396A (en) * 1995-06-30 1997-02-05 Ultrapointe Corporation Method for characterizing defects on semiconductor wafers
JP3564531B2 (ja) * 1999-05-07 2004-09-15 株式会社モノリス 画像領域の追跡方法および装置
DE10000185A1 (de) * 2000-01-05 2001-07-12 Philips Corp Intellectual Pty Verfahren zur Darstellung des zeitlichen Verlaufs des Blutflusses in einem Untersuchungsobjekt

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5301271A (en) * 1991-07-11 1994-04-05 Matsushita Electric Industrial Co., Ltd. Image processing separately processing two-density-level image data and multi-level image data
US5273040A (en) * 1991-11-14 1993-12-28 Picker International, Inc. Measurement of vetricle volumes with cardiac MRI
US6175371B1 (en) * 1995-06-02 2001-01-16 Philippe Schoulz Process for transforming images into stereoscopic images, images and image series obtained by this process
JPH0935058A (ja) * 1995-07-17 1997-02-07 Nec Corp 画像認識方法
JPH09204532A (ja) * 1996-01-25 1997-08-05 Hitachi Ltd 画像認識方法および画像表示方法
EP1054384A2 (en) * 1999-05-20 2000-11-22 Mitsubishi Denki Kabushiki Kaisha Method and apparatus for translating and interfacing voxel memory addresses

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100419794C (zh) * 2005-12-16 2008-09-17 北京中星微电子有限公司 一种图像问题的检测方法
JP2009104515A (ja) * 2007-10-25 2009-05-14 Fujitsu Ltd 差異強調プログラム、差異強調処理方法及び差異強調処理装置
US8311320B2 (en) 2007-10-25 2012-11-13 Fujitsu Limited Computer readable recording medium storing difference emphasizing program, difference emphasizing method, and difference emphasizing apparatus
JP2012103144A (ja) * 2010-11-11 2012-05-31 Univ Of Tokyo 3次元環境復元装置、3次元環境復元方法、及びロボット
JP2016103263A (ja) * 2015-08-31 2016-06-02 株式会社大林組 構造物状況把握支援装置、構造物状況把握支援方法及びプログラム
JP2018005862A (ja) * 2016-07-08 2018-01-11 富士通株式会社 ファセット化処理プログラム、ファセット抽出プログラム、ファセット化処理方法、ファセット抽出方法および情報処理装置
JP2018008403A (ja) * 2016-07-12 2018-01-18 学校法人慶應義塾 立体物製造装置、立体物製造方法及びプログラム
JP2018181056A (ja) * 2017-04-17 2018-11-15 富士通株式会社 差分検知プログラム、差分検知装置、差分検知方法
JP2021021672A (ja) * 2019-07-30 2021-02-18 日本電気通信システム株式会社 距離計測装置、システム、方法、及びプログラム

Also Published As

Publication number Publication date
KR100639139B1 (ko) 2006-10-30
CN1630886A (zh) 2005-06-22
JP4001600B2 (ja) 2007-10-31
EP1519318A4 (en) 2008-11-19
EP1519318A1 (en) 2005-03-30
CN100388317C (zh) 2008-05-14
US20050068317A1 (en) 2005-03-31
KR20050010982A (ko) 2005-01-28
JPWO2004003850A1 (ja) 2005-10-27

Similar Documents

Publication Publication Date Title
WO2004003850A1 (ja) 3次元イメージ比較プログラム、3次元イメージ比較方法、および3次元イメージ比較装置
US8059120B2 (en) Image process apparatus for three-dimensional model
US7576737B2 (en) Image processing device and program
JPH1196374A (ja) 3次元モデリング装置、3次元モデリング方法および3次元モデリングプログラムを記録した媒体
CN102549622B (zh) 用于处理体图像数据的方法
KR20140142470A (ko) 나무 모델과 숲 모델 생성 방법 및 장치
CN110809788B (zh) 深度图像的融合方法、装置和计算机可读存储介质
US20200082618A1 (en) Isosurface generation method and visualization system
Besançon et al. An evaluation of visualization methods for population statistics based on choropleth maps
JP3611583B2 (ja) 特徴線の抽出を伴う三次元イメージ情報の処理装置
GB2399703A (en) Volumetric representation of a 3D object
CN111666788B (zh) 一种图像处理方法、装置及设备、存储介质
Rasoulzadeh et al. Strokes2Surface: Recovering Curve Networks From 4D Architectural Design Sketches
US7388584B2 (en) Method and program for determining insides and outsides of boundaries
JP4014140B2 (ja) 三次元モデリング装置及びその方法及びそのプログラム
EP3467764A1 (en) Image processing method and image processing apparatus
Balzer et al. Volumetric reconstruction applied to perceptual studies of size and weight
Xia et al. A novel approach for computing exact visual hull from silhouettes
Le Van et al. An effective RGB color selection for complex 3D object structure in scene graph systems
CN111462199A (zh) 基于gpu的快速散斑图像匹配方法
CN118334094B (zh) 一种基于三维点云的模型配准方法
CN117495693A (zh) 用于内窥镜的图像融合方法、系统、介质及电子设备
Tattersall et al. Reconstructing Creative Lego Models
CN117765115A (zh) 体数据可视化方法、系统和存储介质
WO2021024367A1 (ja) 形状データ処理装置、形状データ処理方法及び形状データ処理プログラム

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): CN JP KR US

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): AT BE CH CY DE DK ES FI FR GB GR IE IT LU MC NL PT SE TR

DFPE Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101)
121 Ep: the epo has been informed by wipo that ep was designated in this application
WWE Wipo information: entry into national phase

Ref document number: 2004517224

Country of ref document: JP

WWE Wipo information: entry into national phase

Ref document number: 2002736203

Country of ref document: EP

Ref document number: 10989464

Country of ref document: US

WWE Wipo information: entry into national phase

Ref document number: 1020047021408

Country of ref document: KR

Ref document number: 20028292359

Country of ref document: CN

WWP Wipo information: published in national office

Ref document number: 1020047021408

Country of ref document: KR

WWP Wipo information: published in national office

Ref document number: 2002736203

Country of ref document: EP