三维图象的比较程序、比较方法及比较装置

摘要

本发明提供一种三维图象的比较程序、比较方法及比较装置，高精度且快速地进行三维图象的比较。首先，计算机比较将CAD数据体元化生成的三维图象(B1)和将计测数据体元化生成的三维图象(B2)，生成三维差异图象(VB1)。然后，对三维差异图象的表面进行高精度的再体元化，生成三维差异细密图象(VB2)。进而，基于三维差异图象(VB1)或三维差异细密图象(VB2)，用不同颜色显示原来的图象和差异部。

说明书

技术领域

本发明涉及一种比较三维图象并显示差异的三维图象比较程序、三维图象比较方法及三维图象比较装置，特别涉及高精度地比较三维图象的体元差异的三维图象比较程序、三维图象比较方法及三维图象比较装置。

背景技术

近年来，利用三维图象表现物体的技术在各种领域里得到应用。该三维图象例如有各种制造业的三维CAD(Computer Aided Design：计算机辅助设计)的设计图象，医疗领域中使用三维超声波图像诊断装置形成的患部立体图象图像等。并且，通过使用这些三维图象，例如可以进行以下事项。

●如果是三维CAD的设计图象，则在进行设计评审或产品检查等时，设计担当者和检查担当者可以容易地确认(制造验证)所制造的产品或试制品和设计图象。

●如果是使用三维超声波图像诊断装置形成的患部立体图象图像，在诊断时医生可以进一步在视觉上确认(把握)患部的异常或身体部位等的变形。

基于这些目的，为了高效且高精度地确认三维图象，需要自动比较三维图象的技术。特别是在制造业中，为了提高产品质量、缩短开发周期，要求进行高精度且快速的计测和比较。将要比较的图象是由利用CAD作成的三维自由曲面构成的设计模型、和利用三维形状计测仪器测定根据该模型制造的产品或部件得到的点组信息或根据该点组信息作成的剖面或界面(surface)信息。

在以往的三维图象比较中，根据某表面的点到另一个表面上的点的距离评价其差异。这些点有利用CT(computed tomography：计算层析成象)扫描仪或三维数字化仪器计测得到的点、或按照某种规律将所提供的表面分割得到的点。表面与表面的距离就是从各点到表面的垂线距离。把该距离信息反映到表面的点，根据其值或与距离成比例的颜色深度或亮度，实现可视化。因此，以下参照图26进行具体说明。

图26是表示物体的三维图象(CAD数据：三维曲面)和三维点组(计测数据)的图。

在图26所示情况下，表示三维图象的表面的方程式、

z＝y³+3x²y  ……(1)

和某点(x₀、y₀、z₀)的距离为{(x-x₀)²+(y-y₀)²+(z-z₀)²}^1/2的最小值。对位于点(x₀、y₀、z₀)附近的数量庞大的表面，全部进行该计算，以其中的最小值作为点和图象的距离。

关于这样进行的以往的比较方法，图27表示其概念图。

图27是三维图象的以往的比较方法的概念图。

以往的比较如图27所示，不能直接比较设计模型的表面数据Im101和照片等的表面数据Im102。因此，在以往的比较中，把表面数据Im102的表面变换为点组，生成点组数据PIm102，然后进行表面数据Im101和点组数据PIm102的比较。

但是，在以往的比较方法中，如果是表面和点组的高精度比较，则成为对象的点组数量庞大。并且，在实际产品的CAD模型中，微小的面多，必须比较的面数量庞大。即，在实际产品的CAD模型中，必须计算的点和表面的组合数量庞大。结果，表现出以下问题。

●计算时间很长。

●难以选择最接近某点的表面，有时提供不正确的值。

并且，实际制造业要求的计测点数如下表(1)所示，无论用哪种方法比较时，计测点数均成为瓶颈，现实中不可能实现。

  种类要求精度    部件尺寸    计测点数塑料制品0.01mm 400×400×20mm 100亿个点(仅部件表面)铸造制品0.1mm 400×400×350mm 500亿个点(含部件内部)

其中，在塑料制品中，把计测点数假定为上述部件的表面尺寸(400×400mm)的5倍。

发明内容

本发明就是鉴于上述情况而提出的，其目的在于，提供一种能够高精度且快速进行三维图象比较的三维图象比较程序、三维图象比较方法及三维图象比较装置。

为了达到上述目的，本发明的第一方式提供例如图1所示的三维图象比较程序。本发明的三维图象比较程序适用于比较三维图象并显示差异的场合。并且，本发明的三维图象比较程序使计算机执行以下处理。

比较预先被体元化的第1三维图象和第2三维图象，生成三维差异图象(步骤S1)。并且针对所述三维差异图象的表面，以比所述比较更精细的精度生成所述第1三维图象和所述第2三维图象的三维差异细密图象(步骤S2)。另外，比较所述第1三维图象的三维差异细密图象和所述第2三维图象的三维差异细密图象，决定所述三维差异图象的差异部的显示方式(步骤S3)。然后，把所述三维差异图象作为体元比较结果以所述显示方式进行显示(步骤S4)。

如果使计算机执行这种三维图象比较程序，则通过计算机首先比较预先被体元化的第1三维图象和第2三维图象，生成三维差异图象。并且，针对所述三维差异图象的表面，以比所述比较更精细的精度生成所述第1三维图象和所述第2三维图象的三维差异细密图象。另外，对所述三维差异图象和所述三维差异细密图象，决定表示体元差异的显示方式。然后，把所述三维差异图象作为体元比较结果以所述显示方式进行显示。

另外，为了达到上述目的，本发明的第二方式提供一种三维图象比较程序，其特征在于，使计算机执行以下处理：比较预先被体元化的第1三维图象和第2三维图象，生成三维差异图象，部分地抽取所述三维差异图象的差异部，对一致的表面的体元，从三维形状的各面计数在垂直方向不一致的体元的体元数，以所述计数的数为基础决定表示体元的差异的显示方式，将决定了所述显示方式的所述三维差异图象以所述显示方式进行显示。

根据这种三维图象比较程序，通过计算机首先比较预先被体元化的第1三维图象和第2三维图象，生成三维差异图象。并且，部分地抽取所述三维差异图象的差异部，对一致的表面体元，从三维形状的各面计数在垂直方向不一致的体元的体元数。进而，以所述计数的数为基础决定表示体元的差异的显示方式。然后，将决定了所述显示方式的所述三维差异图象以所述显示方式进行显示。

另外，为了达到上述目的，本发明的第三方式提供一种三维图象比较程序，其特征在于，使计算机执行以下处理：比较预先被体元化的第1三维图象和第2三维图象，生成三维差异图象，算出所述三维差异图象的差异部的表面体元数和整体体元数的比率，针对所述三维差异图象，以所述比率为基础决定表示体元的差异的显示方式，对决定了所述显示方式的所述三维差异图象，以与所述比率对应的规定的显示方式显示所述差异部。

根据这种三维图象比较程序，通过计算机首先比较预先被体元化的第1三维图象和第2三维图象，生成三维差异图象。然后，算出所述三维差异图象的差异部的表面体元数和整体体元数的比率，针对所述三维差异图象，以所述比率为基础决定表示体元的差异的显示方式。然后，对决定了所述显示方式的所述三维差异图象，以与所述比率对应的规定的显示方式显示所述差异部。

本发明的上述及其他目的、特征和优点，通过表示作为本发明的示例的优选实施方式的附图及以下相关说明将更加明确。

附图说明

图1是第1实施方式的原理结构图。

图2是第1实施方式的三维图象的比较方法的概念图。

图3是表示可以执行三维图象比较程序的计算机的硬件结构的示例图。

图4是表示三维图象比较装置的功能结构的功能方框图。

图5是表示三维图象比较程序的比较处理的概要图。

图6是表示利用更精细的体元比较进行了比较的三维图象的步骤的图。

图7是表示在设定由体元差异数决定颜色的方法中决定相同体元的颜色的图。

图8是表示在设定由体元差异数决定颜色的方法中决定差异体元(基准图象)颜色的图。

图9是表示在设定由体元差异数决定颜色的方法中决定差异体元(比较对象图象)的颜色的图。

图10是表示体元的要素的序号和节点的关系的图。

图11是表示三维图象上的节点的一般数据结构的图。

图12是表示体元的节点的体元法的数据结构的图。

图13是表示三维图象比较处理的整体流程的流程图。

图14是表示图13的三维图象比较处理的比较信息输入处理的流程图。

图15是表示图13的三维图象比较处理的体元处理的流程图。

图16是表示图13的三维图象比较处理的图象定位处理的流程图。

图17是表示图13的三维图象比较处理的第1图象比较处理的流程图。

图18是表示第2实施方式的三维图象比较程序的比较处理的概要图。

图19是表示图13的三维图象比较处理的第2图象比较处理的流程图。

图20是表示第3实施方式的三维图象比较程序的比较处理的差异评价(1)的图。

图21是表示第3实施方式的通过三维图象比较程序的比较处理决定的颜色(浅颜色)的图。

图22是表示第3实施方式的三维图象比较程序的比较处理的差异评价(2)的图。

图23是表示第3实施方式的通过三维图象比较程序的比较处理决定的颜色(深颜色)的图。

图24是表示第3实施方式的三维图象比较程序的比较处理的差异评价(3)的图。

图25是表示图13的三维图象比较处理的第3图象比较处理的流程图。

图26是表示物体的三维图象(CAD数据：三维曲面)和三维点组(计测数据)的图。

图27是表示三维图象的以往的比较方法的概念图。

具体实施方式

本发明例如可以适用于对以下三维图象进行几何形状比较的方法。

●用于制造有形物体的CAD数据

●计算机游戏等的CG(Computer Graphics：计算机图形)图像

●计测有形物体(工业产品、自然物、生物及其内脏器官等)的结果

几何学形状的比较是这样进行的：把两个三维图象均利用均匀的立方体或长方体(即体元(Voxel))表现，对各自的每个体元实施布尔运算，评价两个图象是否一致。此处，两个图象以基准点或特征形状为基础进行定位。进而，利用不一致的体元数表示图象的差异。

并且，在比较CAD等三维界面的情况下，在该界面上显示不一致的体元的信息。

关于这种比较方法可以考虑以下三种处理。

●以体元单位比较三维图象，以规定的显示方式显示三维图象的差异。即，利用点(体元的颜色)表现三维图象的差异。

●对基准体元的各面，利用规定的显示方式显示从基准体元的各面(x方向、y方向、z方向)起向深度方向的差异体元数。即，利用基准体元的各面(与差异体元数的深度对应的颜色)表现三维图象的差异。

●以差异部的整体体积和表面积的比率为基础，利用规定的显示方式显示三维图象的差异。即，利用整个差异部(与比率对应的颜色)表现三维图象的差异。

因此，在以下实施方式中，首先，作为第1实施方式，说明可以利用点(体元的颜色)表现三维图象的差异的三维图象比较程序、三维图象比较方法及三维图象比较装置。然后，作为第2实施方式，说明可以利用基准体元的各面(与差异体元数的深度对应的颜色)表现三维图象的差异的三维图象比较程序、三维图象比较方法及三维图象比较装置。并且，作为第3实施方式，说明可以利用整个差异部(与比率对应的颜色)表现三维图象的差异的三维图象比较程序、三维图象比较方法及三维图象比较装置。

以下，参照附图说明本发明的实施方式。

(第1实施方式)

首先，说明第1实施方式。第1实施方式以体元为单位比较三维图象，利用规定的显示方式显示三维图象的差异。

图1是第1实施方式的原理结构图。第1实施方式的三维图象比较程序使计算机执行三维图象的比较/差异显示。计算机进行的比较/差异显示是按以下步骤进行的。即，通过计算机比较将设计模型体元化的三维图象B1和将计测数据体元化的三维图象B2，并利用规定的显示方式进行显示的步骤。

首先，计算机比较预先被体元化的三维图象B1和三维图象B2，生成三维差异图象VB1(步骤S1)。

并且，计算机对三维差异图象VB1的表面，以比步骤S1的比较更精细的精度生成三维图象B1和三维图象B2的三维差异细密图象VB2(步骤S2)。

然后，计算机比较三维图象B1的三维差异细密图象和三维图象B2的三维差异细密图象，决定三维差异图象VB1的差异部的显示方式(步骤S3)。另外，所说的该显示方式是以红色、绿色、蓝色的光三原色为基础的加法混色(RGB)。并且，显示方式也可以是以青色、品红色、黄色的补色三原色为基础的减法混色(CMY)。此处，在使用减法混色的情况下，除了青色、品红色、黄色的补色三原色以外，还可以追加黑色。

并且，计算机把决定了显示方式的三维差异图象VB1作为体元比较结果，利用规定的显示方式显示在显示器等上(步骤S4)。另外，体元比较结果不仅可以显示在显示器等上，也可以在打印机或绘图仪等上进行打印。

这样，在第1实施方式中，首先，通过计算机比较预先被体元化的三维图象B1和三维图象B2，生成三维差异图象VB1。

并且，通过计算机对三维差异图象VB1的表面，以比步骤S1的比较更精细的精度，生成三维图象B1和三维图象B2的三维差异细密图象VB2。

然后，通过计算机比较三维图象B1的三维差异细密图象和三维图象B2的三维差异细密图象，决定三维差异图象VB1的差异部的显示方式。

并且，通过计算机把决定了显示方式的三维差异图象VB1作为体元比较结果，利用规定的显示方式显示在显示器上等。

由此，在三维图象的比较中，可以高精度且快速地进行比较处理。

下面，参照图2说明第1实施方式的三维图象的比较方法的概念。

图2是第1实施方式的三维图象的比较方法的概念图。

第1实施方式的比较方法如图2所示，作为设计模型的表面数据Im11和其他表面数据Im12不能直接进行比较。因此，在第1实施方式的比较方法中，把分别包围表面数据Im11、Im12的体积变换为体元，在分别生成体元B11、B12后，进行体元B11和体元B12的比较。由此，可以容易地进行两个三维图象的比较。

以下，具体说明本发明的第1实施方式。

图3是表示可以执行三维图象比较程序的计算机的硬件结构的示例图。计算机100通过CPU(Central Processing Unit：中央处理单元)101控制着整个装置。CPU101通过总线107连接着RAM(Random AccessMemory：随机存取存存储器)102、硬盘驱动器(HDD：Hard Disk Drive)103、图象处理装置104、输入接口105、和通信接口106。

在RAM102中临时存储使CPU101执行的OS(Operating System：操作系统)程序和应用程序的至少一部分。并且，在RAM102中存储CPU101进行处理所需的各种数据。在HDD103中存储OS、三维图象比较程序等应用程序。

图象处理装置104连接着监视器P111。并且，图象处理装置104按照来自CPU101的指令，使监视器P111的画面显示图像。输入接口105连接着键盘P112和鼠标P113。并且，输入接口105把从键盘P112、鼠标P113发送来的信号通过总线107发送给CPU101。

通信接口106连接着网络110。通信接口106通过网络110进行与其他计算机之间的数据收发。

通过使上述硬件结构的计算机100执行三维图象比较程序，计算机100发挥三维图象比较装置的作用。以下，说明计算机通过执行三维图象比较程序而实现的三维图象比较装置的处理功能。

图4是表示三维图象比较装置的功能结构的功能方框图。

根据图4，本发明的三维图象比较装置由以下部分构成：输入比较信息的比较信息输入部10；根据所输入的三维图象生成体元的体元处理部20；将被体元化了的三维图象定位的图象定位处理部30；评价体元的差异并显示的体元差异评价处理部40；以及将显示的体元叠合的体元叠合处理部50。此处，比较信息输入部10除了自身的功能之外，具有：计测精度输入部11、计测精度存储部12、三维界面输入部13、图象存储部14、点组处理切换部15、点组(表面)输入部16、点组(表面)存储部17、点组(体积)输入部18、和点组(体积)存储部19。并且，体元处理部20除了自身的功能之外，具有：体元产生部21、体元产生部22、和体元存储部23。以下，详细叙述这些功能。

比较信息输入部10连接着体元处理部20、三维CAD图象DB11、和三维计测图象DB12，输入作为比较信息的三维图象的CAD数据和计测数据。此处，比较信息输入部10接受使用者进行的三维图象的输入指示。并且，比较信息输入部10判断所接受的三维图象的图象种类。然后，比较信息输入部10判定是否在图象存储部14、点组(表面)存储部17和点组(体积)存储部19中存储有三维图象。

计测精度输入部11连接着计测精度存储部12，接受使用者输入的计测精度。另外，所说的计测精度可以由使用者任意决定，是体元化的三维图象的最小单位。并且，计测精度例如优选使用制造的产品的要求精度(0.01mm等)。

计测精度存储部12连接着计测精度输入部11，存储由计测精度输入部11接受的计测精度。

三维界面输入部13连接着图象存储部14和三维CAD图象DB11，如果使用者指示的三维图象的图象种类是三维界面，则从三维CAD图象DB11接受输入的三维界面。

图象存储部14连接着三维界面输入部13，把由三维界面输入部13接受的三维界面存储在内部的存储装置中。

点组处理切换部15连接着点组(表面)输入部16和点组(体积)输入部18，进行点组处理的切换。此处，点组处理切换部15在使用者指示的三维图象的图象种类是计测表面数据或计测剖面数据的情况下，进行点组(表面)输入部16和点组(体积)输入部18的切换。

点组(表面)输入部16连接着点组处理切换部15、点组(表面)存储部17和三维计测图象DB12，接受输入的计测表面数据。此处，如果使用者指示的三维图象的图象种类是计测表面数据，则点组(表面)输入部16从三维计测图象DB12接受输入的计测表面数据。

点组(表面)存储部17连接着点组(表面)输入部16，存储计测表面数据。此处，点组(表面)存储部17把由点组(表面)输入部16接受的计测表面数据存储在内部的存储装置中。

点组(体积)输入部18连接着点组处理切换部15、点组(体积)存储部19和三维计测图象DB12，接受输入的计测剖面数据。此处，如果使用者指示的三维图象的图象种类是计测剖面数据，则点组(体积)输入部18从三维计测图象DB12接受输入的计测剖面数据。

点组(体积)存储部19连接着点组(体积)输入部18，存储计测剖面数据。此处，点组(体积)存储部19把由点组(体积)输入部18接受的计测剖面数据存储在内部的存储装置中。

体元处理部20连接着比较信息输入部10和图象定位处理部30，根据由比较信息输入部10输入的三维图象生成利用体元形成的图象。此处，体元处理部20选择由比较信息输入部10记录的比较信息。并且，体元处理部20以由比较信息输入部10记录的计测精度为基础，决定体元尺寸。然后，体元处理部20接受由比较信息输入部10存储的三维图象的图象种类，对下一个处理进行判断。

根据该判断，如果图象种类是三维界面，则体元处理部20从比较信息输入部10的图象存储部14选择/读入三维界面，把该三维界面传送给体元产生部21。

另外，如果图象种类是计测表面数据，则体元处理部20从比较信息输入部10的点组(表面)存储部17选择/读入计测表面数据，把该计测表面数据传送给体元产生部22。

另外，如果图象种类是计测剖面数据，则体元处理部20从比较信息输入部10的点组(体积)存储部19选择/读入计测剖面数据，把该计测剖面数据传送给体元产生部22。

然后，体元处理部20判定是否已由体元产生部21、体元产生部22生成体元。

体元产生部21以传送来的三维界面为基础，生成利用体元形成的图象。

体元产生部22以传送来的计测表面数据为基础，生成利用体元形成的图象。并且，体元产生部22以传送来的计测剖面数据为基础，生成利用体元形成的图象。

图象定位处理部30连接着体元处理部20和体元差异评价处理部40，定位由体元处理部20变换为体元的两个图象。此处，图象定位处理部30选择预先附加在两个三维图象上的基准点。另外，所说的这两个三维图象是指作为CAD数据的三维界面和作为计测数据的计测表面数据或计测剖面数据。并且，图象定位处理部30根据所选择的基准点计算偏移。另外，图象定位处理部30以所计算的偏移为基础，移动体元变换后的三维界面。然后，图象定位处理部30把移动后的三维界面存储在内部的存储装置中。

体元差异评价处理部40连接着图象定位处理部30和体元叠合处理部50，比较并评价由图象定位处理部30定位的两个图象，利用规定的显示方式显示所评价的两个图象。另外，所说的规定的显示方式是指例如根据图象的差异变更RGB中各颜色的深度进行显示。此处，体元差异评价处理部40选择表面体元。通过该选择，体元差异评价处理部40判断所选择的体元是否是表面体元。并且，体元差异评价处理部40如果通过该判断选择了表面体元，则以评价精度进行再体元化。即，根据评价精度，使表面体元形成为更精细的体元。另外，评价精度可以由使用者预先决定。另外，评价精度也可以从计测精度的数倍起分阶段地使之精细化。

并且，体元差异评价处理部40以被变换为体元的三维界面和计测表面/剖面数据为基础，比较并计数差异体元。另外，关于该比较处理的详细情况将在后面叙述。

另外，体元差异评价处理部40以所计数的差异体元为基础，决定作为显示方式的颜色的深度，进行体元的着色。

并且，体元差异评价处理部40判定体元是否已被着色。即，体元差异评价处理部40判定被体元化的三维图象是否已被比较。

体元叠合处理部50连接着体元差异评价处理部40，叠合由体元差异评价处理部40显示的多个图象。

三维CAD图象DB11连接着比较信息输入部10，存储作为CAD数据的三维界面。

三维CAD图象DB12连接着比较信息输入部10，存储作为计测数据的计测表面数据和计测剖面数据。

下面，参照图5～图9具体说明第1实施方式的三维图象比较程序的比较处理。

第1实施方式首先利用体元比较三维形状，由此实现计算时间的减少。例如，在被三维曲面包围的CAD图象的形状中，利用细微、均匀的立方体(体元)分割图象，通过使其内部的体元开(on)、外部体元关(off)来表现三维图象。并且，利用X线CT扫描仪计测的图象使其X线密度高的体元开，使密度低的体元关。此处，参照图5说明体元的比较处理的概要。

图5是表示三维图象比较程序的比较处理的概要图。

如图5所示，比较处理把两个图象，象三维图象B21和三维图象B22那样体元化，对其重叠的部分(三维差异图象VB21)进行布尔运算，从而可以可视化地表现其差异。由此，比较的运算量为以往方法的几十分之一。另外，在该布尔运算中，例如，通过使三维图象B21的体元和三维图象B22的体元作求异运算(XOR)，可以挑出不一致的体元。

然后，利用比评价精度粗糙的体元表现将要比较的两个三维图象，之后对每个粗糙体元产生精细体元，用该差异体元数表现两个图象的差异。并且，为了可以可视化地识别差异，在一个图象上叠合该数值。关于该结构，参照图6进行说明。

图6是表示利用更精细的体元比较所比较的三维图象的步骤的图。

如图6所示，首先产生粗糙体元(三维差异图象VB21)。然后，细密地分割包括表面的体元(三维差异细密图象VB22)，利用该差异体元数评价两个图象的差异。另外，在图6的示例中，两个图象的差异为“8”。这表示例如在粗糙体元内，图象I21、I22的三维差异细密图象VB22的细密体元的差异。

然后，以差异数为基础决定体元的颜色深度。作为结果，能够可视地理解整体差异。例如图6所示，在体元Bx11、Bx12、Bx13中，已决定出该部分的颜色(例如蓝色的深度)，所以可以判断图象I22有差异。并且，在体元Bx16中，已决定出该部分的颜色(例如红色的深度)，所以判明图象I21有差异。并且，在体元Bx14、Bx15中，未决定出该部分的颜色，所以可以判断为没有差异。

这样，在图6的示例中，以相似比为1/4的体元再细分体元。即，三维时，最多有64(4×4×4)体元的差异。关于设定根据该差异数决定的颜色的方法，参照图7～图9进行说明。

图7是在设定由体元差异数决定的颜色的方法中表示决定相同体元的颜色的图。

根据图7，在体元相同的情况下，即没有差异的情况下，例如，将设定位置ST1的RGB三色的亮度均设定为最大值(例如，255，255，255)。

图8是在设定由体元差异数决定的颜色的方法中表示决定差异体元(基准图象)的颜色的图。

根据图8，在基准图象的体元相对比较对象图象为凹状的情况下，例如，根据差异体元数，减少设定位置ST2的绿色(G)的值。

图9是在设定由体元差异数决定的颜色的方法中表示决定差异体元(比较对象图象)的颜色的图。

根据图9，在比较对象图象的体元相对基准图象为凹状的情况下，例如，根据差异体元数，减少设定位置ST3的红色(R)的值。

以上，如图7～图9所示，颜色与各个图象的体元数成比例地决定。

下面，参照图10～图12说明在三维图象比较装置中使用的数据的结构。

图10是表示体元的要素的序号和节点的关系图。

根据图10，从三维图象被变换的体元例如是8个体元Bx1～Bx8。该情况下，体元Bx1由节点(articulation point)(1，2，5，4，10，11，14，13)构成。另外，节点的顺序在上下各面中逆时针计数。并且，在体元的数据结构中，登录有这些各节点的坐标信息。因此，参照图11和图12说明该坐标信息的详细情况。

图11是表示三维图象的节点的一般数据结构的图。另外，该图11表示典型的有限元分析中的数据结构。

根据图11，体元的数据结构D1由以下部分构成：表示数据的类别的数据类别区域D11；表示各数据类别的序列编号的数据编号区域D12；表示各数据类别的数据位置信息的位置信息区域D13。

数据类别区域D11例如如图11所示，由“GRID”、“CHEXA”这两部分构成。此处，“GRID”表示节点，“CHEXA”表示分割三维图象的要素。

在数据编号区域D12中，例如如图11所示，对数据类别“GRID”赋予了1～27的序列编号。该序列编号1～27是图10中的所述节点序号。并且，在数据编号区域D12中，例如如图11所示，对数据类别“CHEXA”赋予了1～8的序列编号。该序列编号1～8是图10中的所述体元的要素序号。

在位置信息区域D13中，例如如图11所示，对数据类别“GRID”的数据序号1设定了“0.000000，0.000000，50.00000”节点的坐标信息。这接着各数据序号(相对于“GRID”称为节点ID)，利用浮点数表现x、y、z的坐标位置。以后按照图示设定。并且，在位置信息区域D13中，例如如图11所示，对数据类别“CHEXA”的数据序号1设定了“17，21，27，22，7，10，26，14”的各要素的节点ID。这接着各数据序号(相对于“CHEXA”称为要素ID)，定义各要素。以后按照图示设定。

图12是表示有关体元的节点的体元法的数据结构的图。另外，该图12表示利用体元法进行三维表现时的数据结构。

根据图12，体元的数据结构D2对于预先定义的各个记录分别设定要素的大小和各要素的节点。另外，要素的大小中，利用浮点数表现x、y、z边的长度。并且，各要素的节点中设定定义各要素的节点。但是，关于前述的要素ID，由于是规则排列所以没必要定义。此处，例如图12所示，记录(RECORD7)中设定有使用指数表示(浮点数表示)的要素的大小“2.50000000e+001，2.50000000e+001，2.50000004e+001”。并且，记录(RECORD9)中设定有各要素的节点“1，2，5，4，10，11，14，13”～“14，15，18，17，23，24，27，26”。另外，上述要素的大小中的“e”表示指数函数的底。例如，在上述示例中，“e”表示10，“e+001”表示10的一次幂。

利用如上构成的三维图象比较装置进行三维图象的比较处理。以下，具体说明该三维图象比较处理的动作。

图13是表示三维图象比较处理的整体流程的流程图。该处理是使三维图象比较装置以电源接通、程序的起动、或者预先决定的任意定时为契机开始动作并通过CPU101执行的处理。以下，按照步骤序号说明图13所示处理。另外，关于本流程图中的各功能的名称以图4为基础进行说明。

(步骤S10)比较信息输入部10进行比较信息输入处理，即输入作为比较信息的三维图象的CAD数据和计测数据。另外，关于比较信息输入处理的详细情况将在后面的图14中说明。

(步骤S20)体元处理部20进行体元处理，即根据在步骤S10输入的三维图象生成使用体元形成的图象。另外，关于体元处理的详细情况将在后面的图15中说明。

(步骤S30)图象定位处理部30进行图象定位处理，即定位在步骤S20中被变换为体元的两个图象。另外，关于图象定位处理的详细情况将在后面的图16中说明。

(步骤S40)体元差异评价处理部40进行图象比较处理，即比较并评价在步骤S30中被定位的两个图象。另外，关于图象比较处理的详细情况将在后面的图17中说明。

(步骤S50)体元差异评价处理部40利用规定的显示方式显示在步骤S40中被比较/评价的多个图象。另外，所说的规定的显示方式，是指例如根据图象的差异变更RGB中的各颜色的深度并显示。

(步骤S60)体元叠合处理部50叠合在步骤S50中显示的多个图象。即，将比较结果叠合在三维图象上。例如，作为比较结果的三维差异图象以被体元化的两个图象为基础，可以仅将通过求异生成的差异部可视化。该情况时，体元叠合处理部50通过将被体元化的两个图象、三维差异图象、或三维差异细密图象叠合，可以用不同颜色显示原有图象和差异部。

图14是表示图13的三维图象比较处理中的比较信息输入处理的流程图。在执行图13的步骤S10(比较信息输入处理)时，按照以下流程进行处理。另外，该比较信息输入处理由三维图象比较装置的比较信息输入部10来进行。

(步骤S101)比较信息输入部10的计测精度输入部11接受由使用者进行的计测精度的输入。

(步骤S102)比较信息输入部10的计测精度存储部12存储在步骤S101中接受的计测精度。

(步骤S103)比较信息输入部10接受由使用者进行的三维图象的输入指示。

(步骤S104)比较信息输入部10判断在步骤S103中接受的三维图象的图象种类。此处，比较信息输入部10在图象种类为三维界面时转入步骤S105，在图象种类为计测表面数据时转入步骤S107，或者在图象种类为计测剖面数据时转入步骤S109。另外，在图象种类为计测表面数据或者计测剖面数据的情况下，由比较信息输入部10的点组处理切换部15进行点组(表面)输入部16和点组(体积)输入部18的切换。

(步骤S105)比较信息输入部10的三维界面输入部13由于在步骤S104中，图象种类是三维界面，所以从三维CAD图象DB11进行三维界面的输入接受。

(步骤S106)比较信息输入部10的图象存储部14把在步骤S105接受的三维界面存储在内部的存储装置中。

(步骤S107)比较信息输入部10的点组(表面)输入部16由于在步骤S104中，图象种类是计测表面数据，所以从三维计测图象DB12进行计测表面数据的输入接受。

(步骤S108)比较信息输入部10的点组(表面)存储部17把在步骤S107接受的计测表面数据存储在内部的存储装置中。

(步骤S109)比较信息输入部10的点组(体积)输入部18由于在步骤S104中，图象种类是计测剖面数据，所以从三维计测图象DB12进行计测剖面数据的输入接受。

(步骤S110)比较信息输入部10的点组(体积)存储部19把在步骤S109接受的计测剖面数据存储在内部的存储装置中。

(步骤S111)比较信息输入部10判定在步骤S106、步骤S108或步骤S110中是否已存储三维图象。此处，比较信息输入部10在没有记录三维图象时，返回步骤S104并重复进行相同的处理，在记录有三维图象时结束本处理，返回图13的步骤S10。

图15是表示图13的三维图象比较处理中的体元处理的流程图。在执行图13的步骤S20(体元处理)时，按照以下流程进行处理。另外，该体元处理由三维图象比较装置的体元处理部20进行。

(步骤S201)体元处理部20选择在步骤S10中存储的比较信息。

(步骤S202)体元处理部20以在步骤S10中存储的计测精度为基础，决定体元尺寸。

(步骤S203)体元处理部20接受在步骤S10中存储的三维图象的种类。

(步骤S204)体元处理部20判断在步骤S203接受的三维图象的图象种类。此处，体元处理部20在图象种类为三维界面时转入步骤S205，在图象种类为计测表面数据时转入步骤S207，或者在图象种类为计测剖面数据时转入步骤S209。

(步骤S205)体元处理部20由于在步骤S204判断图象种类是三维界面，所以从比较信息输入部10的图象存储部14选择/读入三维界面，把该三维界面传送给体元产生部21。

(步骤S206)体元处理部20的体元产生部21以在步骤S205中传送的三维界面为基础，生成使用体元形成的三维图象。

(步骤S207)体元处理部20由于在步骤S204中判断出图象种类是计测表面数据，所以从比较信息输入部10的点组(表面)存储部17选择/读入计测表面数据，把该计测表面数据传送给体元产生部22。

(步骤S208)体元处理部20的体元产生部22以在步骤S207中传送的计测表面数据为基础，生成使用体元形成的三维图象。

(步骤S209)体元处理部20由于在步骤S204中判断出图象种类是计测剖面数据，所以从比较信息输入部10的点组(体积)存储部19选择/读入计测剖面数据，把该计测剖面数据传送给体元产生部22。

(步骤S210)体元处理部20的体元产生部22以在步骤S209传送的计测剖面数据为基础，生成使用体元形成的三维图象。

(步骤S211)体元处理部20判定在步骤S206、步骤S208或步骤S210是否已生成使用体元形成的三维图象。此处，体元处理部20在没有生成体元时，返回步骤S204并重复进行相同的处理，在生成体元时结束本处理，返回图13的步骤S20。

图16是表示图13的三维图象比较处理中的图象定位处理的流程图。在执行图13的步骤S30(图象定位处理)时，按照以下流程进行处理。另外，该图象定位处理由三维图象比较装置的图象定位处理部30进行。

(步骤S301)图象定位处理部30选择预先附加给两个三维图象的基准点。另外，所说的这两个三维图象是指作为CAD数据的三维界面和作为计测数据的计测表面数据或计测剖面数据。

(步骤S302)图象定位处理部30根据在步骤S301中选择的基准点计算偏移。

(步骤S303)图象定位处理部30以在步骤S302中计算的偏移为基础，移动被变换为体元的三维图象。

(步骤S304)图象定位处理部30把在步骤S303中移动了的三维图象存储在内部的存储装置中，结束本处理并返回图13的步骤S30。

下面，说明三维图象比较处理，本发明对应于第1～第3实施方式存在多个三维图象比较处理。因此，此处把第1实施方式的三维图象比较处理作为第1图象比较处理进行说明。

图17是表示图13的三维图象比较处理中的第1图象比较处理的流程图。在执行图13的步骤S40(图象比较处理)时，按照以下流程进行处理。另外，该第1图象比较处理由三维图象比较装置的体元差异评价处理部40进行。

(步骤S411)体元差异评价处理部40选择表面体元。另外，表面体元的选择，例如对被体元化的三维图象，使表面体元的内侧体元为开、使外侧体元关，从而可以容易地判断。

(步骤S412)体元差异评价处理部40判断在步骤S411中选择的体元是否是表面体元。此处，体元差异评价处理部40在选择了表面体元时转入步骤S413，在没有选择表面体元时返回步骤S411。

(步骤S413)体元差异评价处理部40由于在步骤S412中选择了表面体元，所以以评价精度进行再体元化。即，根据评价精度，对表面体元进行更精细的体元化。

(步骤S414)体元差异评价处理部40以被变换为体元的三维界面和计测表面/剖面数据为基础，比较差异体元并计数。

(步骤S415)体元差异评价处理部40以在步骤S414中计数的差异体元为基础，决定作为显示方式的颜色的深度，进行体元的着色。

(步骤S416)体元差异评价处理部40判定在步骤S415中体元是否已被着色。即，体元差异评价处理部40判定被体元化的三维图象是否已被比较。此处，体元差异评价处理部40在体元未被着色时返回步骤S411并重复进行相同的处理，在体元已被着色时结束本处理，并返回图13的步骤S40。

这样，在第1实施方式中，首先，通过计算机比较预先被体元化的三维图象B21和三维图象B22，生成三维差异图象VB21。

并且，通过计算机对三维差异图象VB21的表面，以比前次比较更精细的精度生成三维图象B21和三维图象B22的三维差异细密图象VB22。

然后，通过计算机比较三维图象B21的三维差异细密图象和三维图象B22的三维差异细密图象，决定三维差异图象VB21的差异部的显示方式。在进行该决定时，确定与体元的差异数对应的颜色。

并且，通过计算机把决定了显示方式的三维差异图象VB21作为体元比较结果，以规定的显示方式显示在显示器等上。

这样，在三维图象的比较中，可以高精度且快速地进行比较处理。并且，通过把三维形状的体元作为图象的比较对象，无需比较数量庞大的点组，就可以实现存储器资源的高效化和运算次数的减少。另外，通过在进行三维图象的比较时仅作表面比较，可以大幅度减少运算次数。

另外，在上述说明中，进行了对三维图象的表面体元(粗糙体元)生成三维差异细密图象并比较的说明，也可以仅对三维图象的表面体元中有差异的粗糙体元，生成三维差异细密图象并比较。

(第2实施方式)

下面，说明第2实施方式。第2实施方式的比较处理形式与第1实施方式不同。即，第2实施方式以体元为单位比较三维图象，并且对三维图象的差异，以与距基准体元各面的差异体元数的深度对应的规定的显示方式，显示基准体元的各面。

因此，参照图18和图19具体说明第2实施方式的三维图象比较程序的比较处理。

图18是表示第2实施方式的三维图象比较程序的比较处理的概要图。第2实施方式的比较处理使用由第1实施方式作成的体元，在三维图象的体元的各面利用颜色显示差异。该颜色如下所述，由各面在垂直方向(深度方向)的体元数决定。

参照图18进行说明，首先，使用在第1实施方式中作成的体元，比较三维图象(步骤S1b)。部分地抽取该所比较的三维差异图象VB21的不一致部分(步骤S2b)。然后，对某个一致的表面体元(基准体元BB1)，从三维形状的各面(zy面、zx面、yx面)计数在垂直方向不一致的体元的体元数(差异体元数的深度)，决定显示方式(步骤S3b)。即，如果以在A-A向视剖面表示的yz面和zx面为例，则评价为在yz面有3个体元的差异，在zx面有2个体元的差异。进而，利用与该计数对应的规定的显示方式，显示一致的表面体元(基准体元)的各面(步骤S4b)。另外，该显示方式可以使用和前述相同的颜色(RGB或CMY)。

这样，在本发明中，通过计算机，首先使用在第1实施方式中作成的体元比较三维图象。并且，部分地抽取比较的三维差异图象VB21的不一致部分。

然后，通过计算机对某个一致的表面体元(基准体元BB1)，从三维形状的各面(zy面、zx面、yx面)计数垂直方向的不一致体元的体元数(差异体元数的深度)，决定显示方式。

并且，通过计算机以与该计数对应的规定的显示方式为基础，显示一致的表面的体元(基准体元)的各面。

上述的概要功能只有图4的功能方框图中的体元差异评价处理部40的功能不同。因此，在第2实施方式中，仅说明体元差异评价处理部40的新功能(包括重复)。

另外，在概要说明中，把与不一致部分(差异部)邻接并且一致的表面的体元作为基准体元进行了说明，但在以下说明中，选择与差异部邻接并且一致的表面体元进行范围确定，把以评价精度对该范围进行了再体元化时的体元作为基准体元。

体元差异评价处理部40的新功能是进行图象比较处理，即比较在图象定位处理部30中被定位的两个图象，利用规定的显示方式显示所比较的两个图象。另外，所说的规定的显示方式，是指例如根据图象的差异变更RGB中的各颜色的深度并显示。此处，体元差异评价处理部40选择表面体元。通过该选择，体元差异评价处理部40判断所选择的体元是否是表面体元。然后，体元差异评价处理部40在通过该判断选择了表面体元时，确定再体元化的范围。即，把图象不一致的部分确定为详细调查的范围。

并且，体元差异评价处理部40在所确定的范围以评价精度进行再体元化。即，根据评价精度将被确定范围的表面体元进行更精细的体元化。并且，体元差异评价处理部40以被变换为细密体元的三维界面和计测表面/剖面数据为基础，计数差异体元。即，对一致的表面的细密体元(基准体元)，从基准体元的三维形状的各面(zy面、zx面、yx面)计数不一致体元的体元数(差异体元数的深度)。

另外，体元差异评价处理部40以所计数的差异体元为基础，决定作为显示方式的颜色的深度，进行体元的着色。即，体元差异评价处理部40对一致的表面的细密体元(基准体元)，决定基准体元的各面的各自显示方式(颜色)。

并且，体元差异评价处理部40判定体元是否已被着色。即，体元差异评价处理部40判定被体元化的三维图象是否已被比较。

下面，参照图19具体说明第2实施方式的三维图象比较程序的比较处理。另外，如前面所述，对应于第1～第3实施方式有多个本发明的三维图象比较处理。因此，此处把第2实施方式的三维图象比较处理作为第2图象比较处理进行说明。

图19是表示图13的三维图象比较处理的第2图象比较处理的流程图。在执行图13的步骤S40(图象比较处理)时，按照以下流程进行处理。另外，该第2图象比较处理由三维图象比较装置的体元差异评价处理部40进行。并且，在概要说明中，把与不一致部分(差异部)邻接并且一致的表面的体元作为基准体元进行了说明，但在以下说明中，选择与差异部邻接并且一致的表面体元进行范围确定，把以评价精度对该范围进行了再体元化时的体元作为基准体元。

(步骤S421)体元差异评价处理部40选择表面体元。另外，表面体元的选择通过例如对被体元化的三维图象，使表面体元的内侧体元开、使外侧体元关，可以容易地进行判断。

(步骤S422)体元差异评价处理部40判断在步骤S421中选择的体元是否是表面体元。此处，体元差异评价处理部40在选择了表面体元时转入步骤S423，在没有选择表面体元时返回步骤S421。

(步骤S423)体元差异评价处理部40由于在步骤S422中选择了表面体元，所以确定再体元化的范围。即，把图象不一致的部分确定为详细调查的范围。

(步骤S424)体元差异评价处理部40对在步骤S423中确定的范围以评价精度进行再体元化。即，根据评价精度将被确定范围的表面体元进行更精细的体元化。

(步骤S425)体元差异评价处理部40再次以被变换为体元的三维界面和计测表面/剖面数据为基础，计数差异体元。即，对一致的表面的体元(基准体元)，从基准体元的三维形状的各面(zy面、zx面、yx面)计数不一致体元的体元数(差异体元数的深度)。

(步骤S426)体元差异评价处理部40以在步骤S425中计数的差异体元为基础，决定作为显示方式的颜色的深度，进行体元的着色。即，体元差异评价处理部40对一致的表面的体元(基准体元)，决定体元各面的各自显示方式(颜色)。

(步骤S427)体元差异评价处理部40判定在步骤S426中是否已将体元着色。即，体元差异评价处理部40判定被体元化的三维图象是否已被比较。此处，体元差异评价处理部40在体元未被着色时返回步骤S421并重复进行相同的处理，在体元已被着色时结束本处理，并返回图13的步骤S40。

这样，第2实施方式的比较处理为了减少用于比较的计算机资源(存储量、盘区域、计算时间)，取代包括深度方向的体元，利用三维形状的各面的颜色显示差异，。由此，计算资源只需原来的几分之一。另外，表现差异的颜色由与深度方向的体元数成比例的深度值决定。另外，该深度值以整体深度的最大值为基准进行了均等化。

另外，在上述说明中，进行了以与不一致部分(差异部)邻接并且一致的表面的体元作为基准体元的说明，以及选择与差异部邻接并且一致的表面体元进行范围确定，对把以评价精度对该范围进行了再体元化时的体元作为基准体元进行了说明，但是也可以简单地把一致的表面体元作为基准体元。

(第3实施方式)

下面，说明第3实施方式。第3实施方式的比较处理形式与第1实施方式不同。即，第3实施方式以体元单位比较三维图象，并且以与三维图象的差异部整体的表面体元数和其体积的比率对应的规定的显示方式显示体元。

因此，参照图20～图25具体说明第3实施方式的三维图象比较程序的比较处理。

图20是表示第3实施方式的三维图象比较程序的比较处理的差异评价(1)的图。在第1实施方式中生成的体元的差异部仅是两个图象不同的部分。该不同部分(岛)由一个或多个体元构成。而第3实施方式的比较处理把这些岛的表面体元数和该岛的体元总数之比作为颜色来显示。该颜色是如下决定的。另外，在以下说明中，简单地把该岛表述为差异部。

参照图20进行说明，在第3实施方式的比较处理中，使用在第1实施方式中生成的体元，比较作为三维图象的图象A和图象B。此处，差异部的形状表述如下。

●体积：12体元

●Y方向的表面积：25体元面

结果，表面体元数/体元总数的比为25/12，表示相对某差异部的体元差异较浅。并且，每一定体积的表面积利用图21所示的颜色设定。

图21是表示第3实施方式中通过三维图象比较程序的比较处理决定的颜色(浅颜色)的图。

根据图21，以RGB为例说明表示差异部的颜色。在该示例中，通过改变蓝色(B)的颜色，表现差异部整体的体元差异。图21的示例为了表示相对某表面的体元差异较浅，例如，设定位置ST11的RGB设定被设定为255(R)、255(G)、223(B)。该223(B)表示浅颜色，表示某差异部整体的体元差异较浅。这例如是96(B)时，某差异部整体的体元差异可以表示为深色。以下，图22表示的是该体元差异较深时的体元差异评价。

图22是表示第3实施方式的三维图象比较程序的比较处理的差异评价(2)的图。另外，图22说明相对于图20体元差异较深的情况。

参照图22进行说明，在第3实施方式的比较处理中，使用在第1实施方式中生成的体元，比较作为三维图象的图象A和图象B。此处，差异部的形状表述如下。

●体积：43体元

●Y方向的表面积：32体元面

结果，表面体元数/体元总数的比为32/43，表示相对某差异部的体元差异较深。并且，每一定体积的表面积利用图23所示的颜色设定。

图23是表示第3实施方式中通过三维图象比较程序的比较处理决定的颜色(深颜色)的图。

根据图23，以RGB为例说明表示差异部的颜色。在该示例中，通过改变蓝(B)的颜色，表现差异部整体的体元差异。图23的示例为了表示相对某表面的体元差异较深，例如，设定位置ST12的RGB设定被设定为255(R)、255(G)、96(B)。该96(B)表示深颜色，表示某差异部整体的体元差异较深。

并且，利用精细的体元对产生差异的部分再次评价差异，对连续产生差异的部位进行相同着色。此时，把设定颜色的部分，存储为将要比较的三维界面的部分界面。图24表示该差异评价的示例。

图24是表示第3实施方式中的三维图象比较程序的比较处理的差异评价(3)的图。此处，图24的分离(粗线)线表示在原有差异体元中利用精细体元比较的差异体元的分离位置。另外，图24表示相对于图20对体元差异进行更精细的差异评价的情况。并且，关于差异评价的方法说明，对更精细的体元进行和图20及图22相同的处理，所以省略。

上述的概要功能只有图4的功能方框图中的体元差异评价处理部40的功能不同。因此，在第3实施方式中，仅说明体元差异评价处理部40的新功能(包括重复)。另外，在以下说明中，选择与不一致部分(差异部)邻接并且一致的表面体元进行范围确定，把以评价精度进行了再体元化时的体元作为基准体元。当然，在第3实施方式中，也可以把与不一致部分(差异部)邻接并且一致的表面体元作为基准体元。并且，也可以简单地把一致的表面体元作为基准体元。

体元差异评价处理部40的新功能是：比较并评价由图象定位处理部30定位的两个图象，利用规定的显示方式显示所比较的两个图象。另外，所说的规定的显示方式是指例如根据图象的差异变更RGB中的各颜色的深度并显示。此处，体元差异评价处理部40选择表面体元。通过该选择，体元差异评价处理部40判断所选择的体元是否是表面体元。并且，体元差异评价处理部40由于选择了表面体元，确定进行再体元化的范围。即，把图象不一致的部分确定为详细调查的范围。

并且，体元差异评价处理部40对所确定的范围以评价精度进行再体元化。即，根据评价精度将被确定范围的表面体元进行更精细的体元化。并且，体元差异评价处理部40以被变换为体元的三维界面和计测表面/剖面数据为基础，计数差异体元。即，求出三维差异图象的差异部整体中的表面体元数和体元总数的比率。

另外，体元差异评价处理部40以所求出的差异体元的比率为基础，决定作为显示方式的颜色的深度，进行体元的着色。即，体元差异评价处理部40决定由不一致体元构成的差异部整体的显示方式(颜色)。

并且，体元差异评价处理部40判定体元是否已被着色。即，体元差异评价处理部40判定被体元化的三维图象是否已被比较。

下面，参照图25具体说明第3实施方式的三维图象比较程序的比较处理。另外，如前面所述，对应于第1～第3实施方式有多个本发明的三维图象比较处理。因此，此处把第3实施方式中的三维图象比较处理作为第3图象比较处理进行说明。

图25是表示图13的三维图象比较处理中的第3图象比较处理的流程图。在执行图13的步骤S40(图象比较处理)时，按照以下流程进行处理。另外，该第3图象比较处理由三维图象比较装置的体元差异评价处理部40进行。

(步骤S431)体元差异评价处理部40选择表面体元。另外，表面体元的选择例如通过对被体元化的三维图象使表面体元的内侧体元开、外侧体元关，可以容易地进行判断。

(步骤S432)体元差异评价处理部40判断在步骤S431中选择的体元是否是表面体元。此处，体元差异评价处理部40在选择了表面体元时转入步骤S433，在没有选择表面体元时返回步骤S431。

(步骤S433)体元差异评价处理部40由于在步骤S432中选择了表面体元，所以进行再体元化的范围确定。即，把图象不一致的部分确定为详细调查的范围。

(步骤S434)体元差异评价处理部40对在步骤S433中确定的范围以评价精度进行再体元化。即，根据评价精度，将被确定范围的表面体元进行更精细的体元化。

(步骤S435)体元差异评价处理部40以被变换为体元的三维界面和计测表面/剖面数据为基础，计数差异体元。即，求出三维差异图象的差异部整体中的表面体元数和体元总数的比率。

(步骤S436)体元差异评价处理部40以在步骤S435中求出的差异体元的比率为基础，决定作为显示方式的颜色的深度，进行体元的着色。即，体元差异评价处理部40决定由不一致体元构成的差异部整体的显示方式(颜色)。

(步骤S437)体元差异评价处理部40判定在步骤S436中是否已对体元进行了着色。即，体元差异评价处理部40判定被体元化的三维图象是否已被比较。此处，体元差异评价处理部40在体元未被着色时返回步骤S431并反复相同的处理，在体元已被着色时结束本处理，并返回图13的步骤S40。

这样，第3实施方式中的比较处理为了减少用于比较的计算机资源(存储量、盘区域、计算时间)，替代包括深度方向的体元，利用差异部整体的颜色显示差异，。由此，计算资源只需原来的几分之一。另外，表面体元数/体元总数之比最大为2、最小为0(准确讲是近似0的微小值)。因此，颜色的深度被设定成把2作为最浅的颜色，把0作为最深的颜色。并且，差异部的大小可以理解为三维显示时的差异部整体的大小。

(第1实施方式的计算资源的评价结果示例)

最后，说明第1实施方式的计算资源的评价结果示例。

此处，第1实施方式的效果有以下3点。

●与以往的方式相比，所需要的计算量绝对减少，获得结果的时间短。

●可以利用现有的计算机资源(存储器、盘量)实现精细的比较。

●具有比较运算高效化的优点。

例如，关于计算资源评价如下。

首先，如果对现有技术的表(1)中的计测点数，使用公式(1)所示的三维界面的方程式计算最短距离，则计算量如下。

●界面上的最近节点的获取…………数百个浮点运算

●从该最近节点到该计测点的距离的计算…………数十个浮点运算

●认为接近的界面的数量…………1～数十浮点运算

因此，该计算为了求出一个最短距离需要计算数万个浮点运算。对100亿个点重复该计算。结果，该计算所花费的时间如果可以利用1GHz的CPU同时进行2个浮点运算，将达到约30小时。

另一方面，本发明的实施方式的方法，对于包括多个界面的体积，把整个空间分割成计测精度的体元，对每个体元仅实施逻辑运算，所以计算是500亿浮点运算。因此，用相同计算机，计算所需时间为25秒。由此，本发明的第1实施方式的方法与以往的方法相比，可以大幅度缩短计算机的处理时间。并且，本发明的第1实施方式的方法把整个空间分割成计测精度的体元，所以能够利用和以往方法相同的精度进行比较计算。

以上说明的处理通过写成计算机程序并由计算机执行，可以实现本发明的功能。并且，在通过计算机执行时，把计算机程序预先存储在计算机内的硬盘等中，加载到主存储器中执行。另外，计算机程序可以预先记录在计算机可读取的记录介质上。作为计算机可读取的记录介质，有磁记录介质、光盘、光磁记录介质、半导体存储器等。磁记录介质有硬盘、软盘(FD)、ZIP(Zip：磁盘的一种)、磁带等。光盘有DVD(DigitalVersatile Disc：数字通用盘)、DVD-RAM(DVD Random Access Memory：随机存取存储器)、Cd-ROM(Compact Disc Read Only Memory：CD只读存储器)、CD-R(CD Recordable：可读CD)、CD-RW(CD-Rewritable：K可写CD)等。光磁记录介质有MO(Magneto Optical Disk：磁光盘)等。半导体存储器有闪存(Flash Memory)等。

并且，在使计算机程序流通的情况下，例如，可以销售记录有各种计算机程序的DVD、CD-ROM等可移动型记录介质。并且，可以把计算机程序存储在服务器的存储装置上，通过网络从服务器向客户发送计算机程序。

如上所述，在本发明的第1方式中，比较被体元化的两个三维图象，以更加精细的精度生成该比较的表面并进行比较，以该比较结果为基础决定显示方式，利用该显示方式进行显示，所以能够高精度且快速地进行三维图象的比较。

并且，在本发明的第2方式中，以比较三维图象的比较结果为基础，对基准体元的各面，决定与从基准体元各面向垂直方向的差异体元数对应的显示方式，利用该显示方式进行显示，所以能够高精度且快速地进行三维图象的比较。

并且，在本发明的第3方式中，以比较三维图象的比较结果为基础，决定与将三维图象的差异部整体的表面积除以体积所得的比率对应的显示方式，利用该显示方式进行显示，所以能够高精度且快速地进行三维图象的比较。

以上简单地说明了本发明的原理。另外，本行业人员可以进行多种变形、变更，本发明不限于如上所示说明的正确结构和应用示例，对应的所有变形示例及其等同物被视为处于附带的权利要求项及其等同物形成的本发明的范围之内。