コード例 #1
0
 def LoadData(self, fileName):
     planar = triPackage.get_data(fileName)
     vertices = []
     segments = []
     segmentsMark = []
     for vertex in planar['vertices']:
         vertices.append((vertex[0], vertex[1]))
     for segment in planar['segments']:
         segments.append((segment[0], segment[1]))
         segmentsMark.append(((segment[0], segment[1]), 1))
     self.ruper = Ruper.Ruper(vertices, segments, segmentsMark)
コード例 #2
0
 def LoadData(self, fileName):
     planar = triPackage.get_data(fileName)
     vertices = []
     segments = []
     segmentsMark = []
     for vertex in planar['vertices']:
         vertices.append((vertex[0], vertex[1]))
     for segment in planar['segments']:
         segments.append((segment[0], segment[1]))
         segmentsMark.append(((segment[0], segment[1]), 1))
     self.ruper = Ruper.Ruper(vertices, segments, segmentsMark)
コード例 #3
0
def tri(x, y):
    # plot the plan
    box = triangle.get_data('box')
    #
    # ax1 = plt.subplot(121, aspect='equal')
    # triangle.plot.plot(ax1, **box)

    t = triangle.triangulate(box, 'pc')

    ax2 = plt.subplot(111)  # , sharex=ax1, sharey=ax1)
    plot.plot(ax2, **t)
    plt.plot(x, y, 'yo-')

    plt.show()
コード例 #4
0
ファイル: layout.py プロジェクト: dlarson314/crease-solver
def foo2():
    # From http://dzhelil.info/triangle/delaunay.html
    # (not my code)
    import triangle
    import triangle.plot as plot

    face = triangle.get_data('face')
    print face

    ax1 = mpl.subplot(121, aspect='equal')
    plot.plot(ax1, **face)

    t = triangle.triangulate(face, 'p')

    ax2 = mpl.subplot(122, sharex=ax1, sharey=ax1)
    triangle.plot.plot(ax2, **t)

    mpl.show()
コード例 #5
0

import triangle
import triangle.plot as plot
import matplotlib.pyplot as plt

box = triangle.get_data('box')

ax1 = plt.subplot(121, aspect='equal')
triangle.plot.plot(ax1, **box)

t = triangle.triangulate(box, 'pc')

ax2 = plt.subplot(122, sharex=ax1, sharey=ax1)
plot.plot(ax2, **t)

plt.show()
コード例 #6
0
ファイル: quality2.py プロジェクト: weaselp/triangle
import matplotlib.pyplot as plt

import triangle as tr

plt.figure(figsize=(8, 7))

la = tr.get_data('la')
ax1 = plt.subplot(311)
tr.plot(ax1, **la)

t = tr.triangulate(la, 'pq')
ax2 = plt.subplot(312, sharex=ax1, sharey=ax1)
tr.plot(ax2, **t)

t = tr.triangulate(la, 'pqa')
ax2 = plt.subplot(313, sharex=ax1, sharey=ax1)
tr.plot(ax2, **t)

plt.show()
コード例 #7
0
import triangle
import triangle.plot as plot
import matplotlib.pyplot as plt

ax = plt.axes()
A = triangle.get_data('A')
t = triangle.triangulate(A, 'p')
plot.plot(ax, **t)
plt.show()
コード例 #8
0
ファイル: quality1.py プロジェクト: PolarNick239/triangle
import matplotlib.pyplot as plt
import triangle
import triangle.plot

spiral = triangle.get_data("spiral")

t = triangle.triangulate(spiral, "a.2")
ax1 = plt.subplot(111, aspect="equal")
triangle.plot.plot(ax1, **t)

plt.show()
コード例 #9
0
ファイル: quality2.py プロジェクト: PolarNick239/triangle
import triangle
import triangle.plot
import matplotlib.pyplot as plt

la = triangle.get_data("la")

ax1 = plt.subplot(311, aspect="equal")
triangle.plot.plot(ax1, **la)

t = triangle.triangulate(la, "pq")
ax2 = plt.subplot(312, sharex=ax1, sharey=ax1)
triangle.plot.plot(ax2, **t)

t = triangle.triangulate(la, "pqa")
ax2 = plt.subplot(313, sharex=ax1, sharey=ax1)
triangle.plot.plot(ax2, **t)

plt.show()
コード例 #10
0
import matplotlib.pyplot as plt

import triangle as tr

box = tr.get_data('box')
t = tr.triangulate(box, 'pc')

tr.compare(plt, box, t)
plt.show()
コード例 #11
0
ファイル: Ruper.py プロジェクト: fagan2888/Rupers_Algorithm
        elif os.operation == 'remove':
            print 'remove'
        else:
            print os.operation
    
    def Show(self):
        planar = {'vertices' : np.array(self.vertices), 'segments' : np.array(self.segments.keys())}
        tri = {'vertices' : np.array(self.vertices), 'triangles' : np.array(self.triangles.keys())}
        plot(plt.axes(), **planar)
        plot(plt.axes(), **tri)
        plt.plot(self.vertices[-1][0], self.vertices[-1][1], 'bo')
        plt.show()
        plt.clf()


if __name__ == '__main__':
    planar = triPackage.get_data('happy')
    vertices = []
    segments = []
    segmentsMark = []
    for vertex in planar['vertices']:
        vertices.append((vertex[0], vertex[1]))
    for segment in planar['segments']:
        segments.append((segment[0], segment[1]))
        segmentsMark.append(((segment[0], segment[1]), 1))
    ruper = Ruper(vertices, segments, segmentsMark)
    while ruper.stage != 4:
        os = ruper.NextStep()
    ruper.Show()
    
コード例 #12
0
        while len(self.queueT) > 0:
            tri = self.queueT.popleft()
            if self.IsSkinny(tri):
                self.InsertCircleCenter(self.vertices[tri[0]],
                                        self.vertices[tri[1]],
                                        self.vertices[tri[2]])

    def Start(self):
        self.Triangulate()
        self.Show()
        self.InitializeSegmentQueue()
        self.InitializeTriangleQueue()
        self.EliminateSegment()
        self.Show()
        self.RemoveOutside()
        self.Show()
        self.InitializeTriangleQueue()
        self.EliminateAngle()

if __name__ == '__main__':
    planar = triangle.get_data('layers')
    planar['segments_type'] = {}
    for segment in planar['segments']:
        a, b = sorted((segment[0], segment[1]))
        planar['segments_type'][(a, b)] = 1
    ruper = Ruper(planar)
    now = time.time()
    ruper.Start()
    print time.time() - now
    ruper.Show()
コード例 #13
0
import triangle
import triangle.plot
import matplotlib.pyplot as plt

box1 = triangle.get_data('bbox.1')
box2 = triangle.triangulate(box1, 'rpaa.5')
triangle.plot.plot(plt.axes(), **box2)
plt.show()
コード例 #14
0
ファイル: delaunay1.py プロジェクト: BiaoLiu/LBI.Solar3Dcity
import triangle
import triangle.plot as plot
import matplotlib.pyplot as plt

spiral = triangle.get_data('spiral')
ax1 = plt.subplot(121, aspect='equal')
plot.plot(ax1, **spiral)

t = triangle.triangulate(spiral)

ax2 = plt.subplot(122, sharex=ax1, sharey=ax1)
plot.plot(ax2, **t)

plt.show()
コード例 #15
0
ファイル: Refinement.py プロジェクト: caow13/Rupers_Algorithm
    def EliminateAngle(self):
        while len(self.queueT) > 0:
            tri = self.queueT.popleft()
            if self.IsSkinny(tri):
                self.InsertCircleCenter(self.vertices[tri[0]], self.vertices[tri[1]], self.vertices[tri[2]])

    def Start(self):
        self.Triangulate()
        self.Show()
        self.InitializeSegmentQueue()
        self.InitializeTriangleQueue()
        self.EliminateSegment()
        self.Show()
        self.RemoveOutside()
        self.Show()
        self.InitializeTriangleQueue()
        self.EliminateAngle()

if __name__ == '__main__':
    planar = triangle.get_data('layers')
    planar['segments_type'] = {}
    for segment in planar['segments']:
        a, b = sorted((segment[0], segment[1]))
        planar['segments_type'][(a, b)] = 1
    ruper = Ruper(planar)
    now = time.time()
    ruper.Start()
    print time.time() - now
    ruper.Show()
    
コード例 #16
0
ファイル: bndries.py プロジェクト: weaselp/triangle
import matplotlib.pyplot as plt

import triangle as tr


def bndry(dict):
    return {k: dict[k] for k in ('vertices', 'segments')}


face = tr.get_data('face.1')
A = tr.triangulate(face, 'rc')
B = tr.triangulate(face, 'rpc')

tr.compare(plt, bndry(A), bndry(B))

plt.show()
コード例 #17
0
ファイル: quality2.py プロジェクト: Geodels/tribad
import triangle
import triangle.plot
import matplotlib.pyplot as plt

la = triangle.get_data('la')

ax1 = plt.subplot(311, aspect='equal')
triangle.plot.plot(ax1, **la)

t = triangle.triangulate(la, 'pq')
ax2 = plt.subplot(312, sharex=ax1, sharey=ax1)
triangle.plot.plot(ax2, **t)

t = triangle.triangulate(la, 'pqa')
ax2 = plt.subplot(313, sharex=ax1, sharey=ax1)
triangle.plot.plot(ax2, **t)

plt.show()
コード例 #18
0
import matplotlib.pyplot as plt

import triangle as tr

box1 = tr.get_data('bbox.1')
box2 = tr.triangulate(box1, 'rpa')
tr.plot(plt.axes(), **box2)
plt.show()
コード例 #19
0
ファイル: convex_hull.py プロジェクト: Geodels/tribad
import triangle
import triangle.plot as plot
import matplotlib.pyplot as plt

dots = triangle.get_data('dots')
pts = dots['vertices']
segs = triangle.convex_hull(pts)

plot.plot(plt.axes(), vertices=pts, segments=segs)

plt.show()
コード例 #20
0
ファイル: holes_cavities.py プロジェクト: weaselp/triangle
import matplotlib.pyplot as plt

import triangle as tr

face = tr.get_data('face')
t = tr.triangulate(face, 'p')

tr.compare(plt, face, t)
plt.show()
コード例 #21
0
import matplotlib.pyplot as plt

import triangle as tr

spiral = tr.get_data('spiral')
t = tr.triangulate(spiral)

tr.compare(plt, spiral, t)

plt.show()
コード例 #22
0
import triangle
import triangle.plot as plot
import matplotlib.pyplot as plt

A = triangle.get_data('A')
t = triangle.triangulate(A, 'pq0D')
plot.plot(plt.axes(), **t)
plt.show()
コード例 #23
0
ファイル: voronoi.py プロジェクト: BiaoLiu/LBI.Solar3Dcity
import triangle
import triangle.plot
import matplotlib.pyplot as plt

pts = triangle.get_data('dots')['vertices']

ax1 = plt.subplot(121, aspect='equal')
triangle.plot.plot(ax1, vertices=pts)
lim = ax1.axis()

points, edges, ray_origin, ray_direct = triangle.voronoi(pts)
d = dict(vertices=points, edges=edges, ray_origins=ray_origin, ray_directions=ray_direct)
ax2 = plt.subplot(122, sharex=ax1, sharey=ax1)
triangle.plot.plot(ax2, **d)
ax2.axis(lim)

plt.show()
コード例 #24
0
import triangle
import triangle.plot as plot
import matplotlib.pyplot as plt

spiral = triangle.get_data('spiral')

ax1 = plt.subplot(221, aspect='equal')
triangle.plot.plot(ax1, vertices=spiral['vertices'])

a = triangle.triangulate(spiral)
ax2 = plt.subplot(222, sharex=ax1, sharey=ax1)
plot.plot(ax2, **a)

b = triangle.triangulate(spiral, 'q')
ax3 = plt.subplot(223, sharex=ax1, sharey=ax1)
plot.plot(ax3, **b)

c = triangle.triangulate(spiral, 'q32.5')
ax4 = plt.subplot(224, sharex=ax1, sharey=ax1)
plot.plot(ax4, **c)

plt.show()
コード例 #25
0
ファイル: braille2.py プロジェクト: bttaylor/TactileMaps
def stl(string,fname,coarseness,pos,angle,label_cnt):
	global printSize
	import triangle
	import triangle.plot
	import matplotlib.pyplot as plt
	import numpy as np
	from shapely.geometry import box, LineString
	#buildBasePolyFile(string,pos,fname,coarseness,angle)


	print '/Users/bttaylor/Documents/TactileMap/PolyFiles/b_' + fname + '_frame'
	frames = triangle.get_data('/Users/bttaylor/Documents/TactileMap/PolyFiles/b_' + fname + '_frame');
	framesT = triangle.triangulate(frames,'p');
	for tri in framesT['triangles']:
		pt1 = Point(framesT['vertices'][tri[0]][0],framesT['vertices'][tri[0]][1])
		pt2 = Point(framesT['vertices'][tri[1]][0],framesT['vertices'][tri[1]][1])
		pt3 = Point(framesT['vertices'][tri[2]][0],framesT['vertices'][tri[2]][1])
		tri_v1 = [pt1.x,pt1.y,base_height]
		tri_v2 = [pt2.x,pt2.y,base_height]
		tri_v3 = [pt3.x,pt3.y,base_height]
		elevation2stl.printTriangle(tri_v1,tri_v2,tri_v3,fname)

	for i in range(1,label_cnt + 1):
		print '/Users/bttaylor/Documents/TactileMap/PolyFiles/b_' + fname + '_' + str(i)
		nondot = triangle.get_data('/Users/bttaylor/Documents/TactileMap/PolyFiles/b_' + fname + '_' + str(i));
		non = triangle.triangulate(nondot,'p');
		for tri in non['triangles']:
			pt1 = Point(non['vertices'][tri[0]][0],non['vertices'][tri[0]][1])
			pt2 = Point(non['vertices'][tri[1]][0],non['vertices'][tri[1]][1])
			pt3 = Point(non['vertices'][tri[2]][0],non['vertices'][tri[2]][1])
			tri_v1 = [pt1.x,pt1.y,base_height]
			tri_v2 = [pt2.x,pt2.y,base_height]
			tri_v3 = [pt3.x,pt3.y,base_height]
			elevation2stl.printTriangle(tri_v1,tri_v2,tri_v3,fname)



	# nondot = triangle.get_data('/Users/bttaylor/Documents/TactileMap/PolyFiles/b_' + fname);
	# #print nondot
	# edges = box(0,0,printSize,printSize);
	# edge_pts = []

	# non = triangle.triangulate(nondot,'p');
	# for tri in non['triangles']:
	# 	pt1 = Point(non['vertices'][tri[0]][0],non['vertices'][tri[0]][1])
	# 	pt2 = Point(non['vertices'][tri[1]][0],non['vertices'][tri[1]][1])
	# 	pt3 = Point(non['vertices'][tri[2]][0],non['vertices'][tri[2]][1])
	# 	tri_v1 = [pt1.x,pt1.y,base_height]
	# 	tri_v2 = [pt2.x,pt2.y,base_height]
	# 	tri_v3 = [pt3.x,pt3.y,base_height]
	# 	elevation2stl.printTriangle(tri_v1,tri_v2,tri_v3,fname)
	# for seg in non['segments']:
	# 	line = LineString(non['vertices'][seg])
	# 	if line.touches(edges):
	# 		XYpt1 = Point(non['vertices'][seg[0]])
	# 		XYpt2 = Point(non['vertices'][seg[1]])
	# 		edge_pts.append(XYpt1.coords[0]);
	# 		edge_pts.append(XYpt2.coords[0]);
	# #print edge_pts
	# edge_pts.append([0,0])
	# edge_pts.append([0,printSize])
	# edge_pts.append([printSize,0])
	# edge_pts.append([printSize,printSize])
	# edge_pt_array = np.array(edge_pts)
	# base = dict(vertices=edge_pt_array);
	# baseT = triangle.triangulate(base);
	# #triangle.plot.compare(plt,base,baseT);


	# for tri in baseT['triangles']:
	# 	pt1 = Point(baseT['vertices'][tri[0]][0],baseT['vertices'][tri[0]][1])
	# 	pt2 = Point(baseT['vertices'][tri[1]][0],baseT['vertices'][tri[1]][1])
	# 	pt3 = Point(baseT['vertices'][tri[2]][0],baseT['vertices'][tri[2]][1])
	# 	tri_v1 = [pt1.x,pt1.y,floor]
	# 	tri_v2 = [pt2.x,pt2.y,floor]
	# 	tri_v3 = [pt3.x,pt3.y,floor]

	# 	elevation2stl.printTriangle(tri_v1,tri_v2,tri_v3,fname)

	#print edges
	pt1 = [0, 0, floor]
	pt2 = [0, 0, base_height]
	pt3 = [0, printSize, floor]
	pt4 = [0, printSize, base_height]
	elevation2stl.printTriangle(pt1,pt2,pt4,fname)
	elevation2stl.printTriangle(pt1,pt3,pt4,fname)
	pt3 = [printSize, 0, floor]
	pt4 = [printSize, 0, base_height]
	elevation2stl.printTriangle(pt1,pt2,pt4,fname)
	elevation2stl.printTriangle(pt1,pt3,pt4,fname)
	pt1 = [printSize, printSize, floor]
	pt2 = [printSize, printSize, base_height]
	elevation2stl.printTriangle(pt1,pt2,pt4,fname)
	elevation2stl.printTriangle(pt1,pt3,pt4,fname)
	pt3 = [0, printSize, floor]
	pt4 = [0, printSize, base_height]
	elevation2stl.printTriangle(pt1,pt2,pt4,fname)
	elevation2stl.printTriangle(pt1,pt3,pt4,fname)
	#print bottom
	pt1 = [0, 0, floor]
	pt2 = [0, printSize, floor]
	pt3 = [printSize, printSize, floor]
	pt4 = [printSize, 0, floor]
	elevation2stl.printTriangle(pt1,pt3,pt2,fname)
	elevation2stl.printTriangle(pt1,pt4,pt3,fname)
コード例 #26
0
ファイル: Refinement.py プロジェクト: caow13/Rupers_Algorithm
            if rmBool == True:
                rmKeys.append(tuple(sorted((tri[0], tri[1], tri[2]))))
        for rm in rmKeys:
            self.delaunay_triangles.pop(rm)


                    

    def Start(self):
        self.Triangulate()
        self.EliminateSegment()
        self.RemoveOutside()
        self.Show()
#        self.RemoveOutside()
#        self.delaunay = self.Triangulate(rem = True)
#        self.EliminateAngle()


if __name__ == '__main__':
    planar = triangle.get_data('A') 
    planar['segments_type'] = {}
    for segment in planar['segments']:
        a, b = sorted((segment[0], segment[1]))
        planar['segments_type'][(a, b)] = 1
    ruper = Ruper(planar)
#    now = time.time()
    ruper.Start()
#    print time.time() - now
#    ruper.Show()
    
コード例 #27
0
ファイル: voronoi1.py プロジェクト: Geodels/tribad
import triangle
import triangle.plot
import matplotlib.pyplot as plt

pts = triangle.get_data('diamond_02_00009')['vertices']

ax1 = plt.subplot(121, aspect='equal')
triangle.plot.plot(ax1, vertices=pts)
lim = ax1.axis()

d = triangle.get_data('diamond_02_00009.1.v')
ax2 = plt.subplot(122, sharex=ax1, sharey=ax1)
triangle.plot.plot(ax2, **d)
ax2.axis(lim)

plt.show()
コード例 #28
0
ファイル: delaunay.py プロジェクト: weaselp/triangle
import matplotlib.pyplot as plt

import triangle as tr

d0 = tr.get_data('dots')
d1 = tr.triangulate(d0)
tr.compare(plt, d0, d1)
plt.show()
コード例 #29
0
ファイル: PSLG.py プロジェクト: PolarNick239/triangle
import triangle
from triangle.plot import plot
import matplotlib.pyplot as plt

A = triangle.get_data("A")
ax = plt.axes()
plot(ax, **A)
plt.show()
コード例 #30
0
import matplotlib.pyplot as plt

import triangle as tr

pts = tr.get_data('dots')['vertices']

A = dict(vertices=pts)

points, edges, ray_origin, ray_direct = tr.voronoi(pts)
B = dict(vertices=points,
         edges=edges,
         ray_origins=ray_origin,
         ray_directions=ray_direct)
tr.compare(plt, A, B)
plt.show()
コード例 #31
0
import matplotlib.pyplot as plt

import triangle as tr

face = tr.get_data('corridor')
t = tr.triangulate(face, 'p')

tr.compare(plt, face, t)
plt.show()
コード例 #32
0
import matplotlib.pyplot as plt

import triangle as tr

box = tr.get_data('double_hex3')
print(box)
t = tr.triangulate(box, 'p')
v = t['vertices'].tolist()
# print(box)
tList = t['triangles'].tolist()
tCoords = []  ## coordinates of the triangles: [[[0,0], [0,1], [1,0]], [...]]
for tL in tList:
    tCoords.append([v[tL[0]], v[tL[1]], v[tL[2]]])

print(tCoords)

tr.compare(plt, box, t)
plt.show()
コード例 #33
0
import triangle
import triangle.plot as plot
import matplotlib.pyplot as plt

box = triangle.get_data('box')

ax1 = plt.subplot(121, aspect='equal')
triangle.plot.plot(ax1, **box)

t = triangle.triangulate(box, 'pc')

ax2 = plt.subplot(122, sharex=ax1, sharey=ax1)
plot.plot(ax2, **t)

plt.show()
コード例 #34
0
import triangle
import triangle.plot as plot
import matplotlib.pyplot as plt


face = triangle.get_data('face')
ax1 = plt.subplot(121, aspect='equal')
plot.plot(ax1, **face)

t = triangle.triangulate(face, 'pq10')

ax2 = plt.subplot(122, sharex=ax1, sharey=ax1)
plot.plot(ax2, **t)

plt.show()
コード例 #35
0
import matplotlib.pyplot as plt

import triangle as tr

A = tr.get_data('A')
t = tr.triangulate(A, 'pq30')
tr.plot(plt.axes(), **t)
plt.show()
コード例 #36
0
        rm = [0] * len(delaunay['triangles'])
        for ind, triangle in enumerate(delaunay['triangles']):
            a, b, c = self.vertices[triangle[0]], self.vertices[
                triangle[1]], self.vertices[triangle[2]]
            u = (a + b + c) / 3.0
            for delta in [-1e-9, 0, 1e-9]:
                v = np.array([1e5, u[1] + delta])
                count = self.CrossCount(u, v)
                if count != None and count % 2 == 0:
                    rm[ind] = 1
        ind = map(lambda x: x[0],
                  filter(lambda x: x[1] == 0, list(enumerate(rm))))
        delaunay['triangles'] = delaunay['triangles'][ind]
        return delaunay

    def Start(self):
        self.delaunay = self.Triangulate()
        self.EliminateSegment()
        self.delaunay = self.Triangulate(rem=True)
        self.EliminateAngle()


if __name__ == '__main__':
    planar = triangle.get_data('A')
    planar['segments_type'] = {}
    for segment in planar['segments']:
        planar['segments_type'][(segment[0], segment[1])] = 1
    ruper = Ruper(planar)
    ruper.Start()
    ruper.Show()
コード例 #37
0
import triangle
import triangle.plot
import matplotlib.pyplot as plt

pts = triangle.get_data('dots')['vertices']

ax1 = plt.subplot(121, aspect='equal')
triangle.plot.plot(ax1, vertices=pts)
lim = ax1.axis()

points, edges, ray_origin, ray_direct = triangle.voronoi(pts)
d = dict(vertices=points,
         edges=edges,
         ray_origins=ray_origin,
         ray_directions=ray_direct)
ax2 = plt.subplot(122, sharex=ax1, sharey=ax1)
triangle.plot.plot(ax2, **d)
ax2.axis(lim)

plt.show()