doc/salome/examples/min_distance.py

   1 # Minimal Distance
   2
   3 import geompy
   4
   5 # Create two curves with three closest points
   6 Vertex_1 = geompy.MakeVertex(0, 0, 0)
   7 Vertex_2 = geompy.MakeVertex(0, 70, 0)
   8 Vertex_3 = geompy.MakeVertex(30, 100, 0)
   9 Vertex_4 = geompy.MakeVertex(60, 70, 0)
  10 Vertex_5 = geompy.MakeVertex(90, 100, 0)
  11 Vertex_6 = geompy.MakeVertex(120, 70, 0)
  12 Vertex_7 = geompy.MakeVertex(120, 0, 0)
  13 Vertex_8 = geompy.MakeVertex(90, -30, 0)
  14 Vertex_9 = geompy.MakeVertex(60, 0, 0)
  15 Vertex_10 = geompy.MakeVertex(30, -30, 0)
  16
  17 geompy.addToStudy( Vertex_1, 'Vertex_1' )
  18 geompy.addToStudy( Vertex_2, 'Vertex_2' )
  19 geompy.addToStudy( Vertex_3, 'Vertex_3' )
  20 geompy.addToStudy( Vertex_4, 'Vertex_4' )
  21 geompy.addToStudy( Vertex_5, 'Vertex_5' )
  22 geompy.addToStudy( Vertex_6, 'Vertex_6' )
  23 geompy.addToStudy( Vertex_7, 'Vertex_7' )
  24 geompy.addToStudy( Vertex_8, 'Vertex_8' )
  25 geompy.addToStudy( Vertex_9, 'Vertex_9' )
  26 geompy.addToStudy( Vertex_10, 'Vertex_10' )
  27
  28 Curve_a = geompy.MakeInterpol([Vertex_2, Vertex_3, Vertex_4, Vertex_5, Vertex_6], False, True)
  29 Curve_b = geompy.MakeInterpol([Vertex_1, Vertex_7, Vertex_8, Vertex_9, Vertex_10], False, True)
  30
  31 geompy.addToStudy( Curve_a, 'Curve_a' )
  32 geompy.addToStudy( Curve_b, 'Curve_b' )
  33
  34 # Get all closest points
  35 [nbSols, listCoords] = geompy.ClosestPoints(Curve_a, Curve_b)
  36
  37 for i in range(nbSols):
  38   v1 = geompy.MakeVertex(listCoords[i*6 + 0], listCoords[i*6 + 1], listCoords[i*6 + 2])
  39   v2 = geompy.MakeVertex(listCoords[i*6 + 3], listCoords[i*6 + 4], listCoords[i*6 + 5])
  40
  41   geompy.addToStudy(v1, 'MinDist_%d_Curve_a'%(i+1))
  42   geompy.addToStudy(v2, 'MinDist_%d_Curve_b'%(i+1))
  43
  44 # Get minimum distance
  45 print "Minimal distance between Curve_a and Curve_b is", geompy.MinDistance(Curve_a, Curve_b)
  46
  47 # Get minimum distance with components along axes
  48 [aDist, DX, DY, DZ] = geompy.MinDistanceComponents(Curve_a, Curve_b)
  49 print "Minimal distance between Curve_a and Curve_b is (", DX, ",", DY, ",", DZ, ")"