平面度是对表面形状的一种度量，用于指示该表面上的所有点是否都在同一个平面上。平面度在几何尺寸和公差（GD&T）中用平行四边形表示，当两个表面必须装配在一起形成紧密密封时，平面度就特别有用。

使用平面度公差是为了确保某个即定表面位于两个假想的、完美的、平行的平面内。换句话说，公差范围位于被测表面上最高和最低可接受点之间。

测量平面度包括分析一个表面，找出它不完全平坦的程度。为此，首先要在表面上获取一些点，以便找到包括所有这些点的两个平行平面的包络。下一步是找到尽可能小的平面夹层，而不考虑方向（因为这些平面可以在空间自由移动）。最远的两个点之间的距离就是平面度。因此，这两个平面之间的空间越窄，这个平面就越平坦。

read_image (Image, 'E:/项目/蓝膜/3d平面度检测/深度/_2022_08_27_09_33_01.tiff')
dev_close_window ()
get_image_pointer1 (Image, Pointer, Type, Width, Height)
dev_open_window (0, 0, Width/8, Height/8, 'black', WindowHandle1)
dev_display (Image)
X_power:=0.18
Y_power:=0.18
Z_power:=	0.0037
Percentage1:=40
Percentage2:=60
Percentage11:=15
Percentage12:=85
dev_get_window (WindowHandle)
dev_set_draw ('margin')
gen_empty_obj (EmptyObject)
for Index1 := 0 to 17 by 1gen_rectangle2 (Rectangle1, 2290+228*Index1, 719, 0, 300, 60)concat_obj (EmptyObject, Rectangle1, EmptyObject)
endfor
for Index1 := 0 to 17 by 1gen_rectangle2 (Rectangle1, 2290+228*Index1, 1757, 0, 300, 60)concat_obj (EmptyObject, Rectangle1, EmptyObject)
endfor
stop ()
union1 (EmptyObject, RegionUnion)median_image (Image, ImageMedian, 'square', 17, 'cyclic')
Flatness_ALL (Image, RegionUnion, X_power, Y_power, Z_power, Percentage1, Percentage2, Percentage11, Percentage12, ObjectModel3D, flatness, result)
dev_display (Image)
dev_display (EmptyObject) 

tryget_image_size (Image, Width, Height)gen_image_surface_first_order (ImageX, 'real', X_power, 0, 0, 0, 0, Width, Height)gen_image_surface_first_order (ImageY, 'real', 0, Y_power, 0, 0, 0, Width, Height)convert_image_type (Image, ImageConverted, 'real')scale_image (ImageConverted, ImageZ, Z_power, 0)reduce_domain (ImageX, ROI_Base, ROI_BaseX)reduce_domain (ImageY, ROI_Base, ROI_BaseY)reduce_domain (ImageZ, ROI_Base, ROI_BaseZ)xyz_to_object_model_3d (ROI_BaseX ,ROI_BaseY,ROI_BaseZ, ObjectModel3D1)**********************去掉最大最小**************connection (ROI_Base, ConnectedRegions)count_obj (ConnectedRegions, Number)tuple_gen_const (0, 0, Point_X)tuple_gen_const (0, 0, Point_Y)tuple_gen_const (0, 0, Point_Z)****************每个矩形区域滤掉一部分*********************for Index := 0 to Number-1 by 1select_obj (ConnectedRegions, ObjectSelected,Index+1)dev_display (ObjectSelected)reduce_domain (ImageX, ObjectSelected, ROI_BaseX1)reduce_domain (ImageY, ObjectSelected, ROI_BaseY1)reduce_domain (ImageZ, ObjectSelected, ROI_BaseZ1)xyz_to_object_model_3d (ROI_BaseX1 ,ROI_BaseY1,ROI_BaseZ1, ObjectModel3D1s)get_object_model_3d_params (ObjectModel3D1s, 'num_points', NumPoints)get_object_model_3d_params (ObjectModel3D1s, 'point_coord_z' , DistanceDistribution)tuple_sort (DistanceDistribution, Sorted)Indexs1:=int(Percentage1*0.01*NumPoints)Indexs2:=int(Percentage2*0.01*NumPoints)select_points_object_model_3d (ObjectModel3D1s, 'point_coord_z', Sorted[Indexs1],  Sorted[Indexs2], ObjectModel3DThresholded)get_object_model_3d_params (ObjectModel3DThresholded, 'point_coord_x', BaseX)get_object_model_3d_params (ObjectModel3DThresholded, 'point_coord_y', BaseY)get_object_model_3d_params (ObjectModel3DThresholded, 'point_coord_z', BaseZ)if(Index == 0)Point_X:=BaseXPoint_Y:=BaseYPoint_Z:=BaseZelsePoint_X:=[Point_X,BaseX]Point_Y:=[Point_Y,BaseY]Point_Z:=[Point_Z,BaseZ]endifendforgen_object_model_3d_from_points (Point_X, Point_Y, Point_Z, ObjectModel3Ds)dev_get_window (WindowHandle)visualize_object_model_3d (WindowHandle, ObjectModel3Ds, [], [], ['lut', 'color_attrib'], ['color1', 'coord_z'], [], [], [], PoseOut)fit_primitives_object_model_3d (ObjectModel3Ds, ['primitive_type','fitting_algorithm'], ['plane','least_squares'], ObjectModel3DOut)visualize_object_model_3d (WindowHandle, [ObjectModel3Ds,ObjectModel3DOut], [], [], ['lut', 'color_attrib'], ['color1', 'coord_z'], [], [], [], PoseOut)get_object_model_3d_params (ObjectModel3DOut, 'primitive_type', GenParamValue)get_object_model_3d_params (ObjectModel3DOut, 'primitive_parameter', PlaneValue)get_object_model_3d_params (ObjectModel3DOut, 'point_coord_x', BaseX)get_object_model_3d_params (ObjectModel3DOut, 'point_coord_y', BaseY)get_object_model_3d_params (ObjectModel3DOut, 'point_coord_z', BaseZ)***************************************获取检测点数tuple_gen_const (0, 0, Points_X)tuple_gen_const (0, 0, Points_Y)tuple_gen_const (0, 0, Points_Z)****************每个矩形区域滤掉一部分*********************for Index := 0 to Number-1 by 1select_obj (ConnectedRegions, ObjectSelected,Index+1)dev_display (ObjectSelected)reduce_domain (ImageX, ObjectSelected, ROI_BaseX1)reduce_domain (ImageY, ObjectSelected, ROI_BaseY1)reduce_domain (ImageZ, ObjectSelected, ROI_BaseZ1)xyz_to_object_model_3d (ROI_BaseX1 ,ROI_BaseY1,ROI_BaseZ1, ObjectModel3D1s)get_object_model_3d_params (ObjectModel3D1s, 'num_points', NumPoints)get_object_model_3d_params (ObjectModel3D1s, 'point_coord_z' , DistanceDistribution)tuple_sort (DistanceDistribution, Sorted)Indexs1:=int(Percentage1*0.01*NumPoints)Indexs2:=int(Percentage2*0.01*NumPoints)select_points_object_model_3d (ObjectModel3D1s, 'point_coord_z', Sorted[Indexs1],  Sorted[Indexs2], ObjectModel3DThresholded)get_object_model_3d_params (ObjectModel3DThresholded, 'point_coord_z', BaseZ)tuple_median (BaseZ, TupleZ)area_center (ObjectSelected, Area, TupleX, TupleY)if(Index == 0)Points_X:=TupleXPoints_Y:=TupleYPoints_Z:=TupleZelsePoints_X:=[Points_X,TupleX]Points_Y:=[Points_Y,TupleY]Points_Z:=[Points_Z,TupleZ]endifendforfit_3d_plane_xyz (Points_X, Points_Y, Points_Z, Ox, Oy, Oz, Nx, Ny, Nz, MeanResidual)*平面度 *点到面的距离 计算base面的平面度d1:=(Nx*Points_X+Ny*Points_Y+Nz*Points_Z-PlaneValue[3])tuple_sort (d1, Sorted1)flatness:=max(d1)-min(d1)result:=1
catch (Exception)result:=-1
endtry
return ()