def testHalfSquareOverSquare(self): """ Sovrapposizione di un mezzo quadrato su un quadrato """ sq = Path([P(0,0), P(10,0), P(17,0), P(80,0), P(100,0), P(100,100), P(0,100), P(0,0)]) hsq = Path([P(100,100), P(0,100), P(0,0)]) * rot(random.random()*2*math.pi) * xlate(P(random.random()*self.sf, random.random()*self.sf)) mr = pathOverPath(hsq, sq, 10) self.assertAlmostEquals(mr.sse/mr.n, 0.0) A = X(mr.A) for p1 in hsq: # controllo della distanza punto-path self.assertAlmostEquals(sq.project(p1*A).dist, 0, tollerance)
def testHalfCircleOverHalfCircle(self): """ La funzione deve sovrapporre due semicerchi identici, dopo che uno di questi e' stato traslato e ruotato """ O1 = P(0,0) O2 = P(100,100) c1 = Path( circle(O1, 100, minpts= 100, maxerr=100)[0:50] ) c2 = Path( c1 ) * rot(random.random()*2*math.pi) * xlate(O2) mr = pathOverPath(c2, c1, 10) self.assertAlmostEquals(mr.sse/mr.n, 0.0) A = X(mr.A) for p1, p2 in zip(c1, c2): # controllo della distanza punto-path self.assertAlmostEquals(c1.project(p2*A).dist, 0, tollerance)
def testQuarterCircleOverHalfCircle(self): """ La funzione deve sovrapporre un quarto-di-cerchio a un semicerchio, dopo che il quarto-di-cerchio e' stato traslato e ruotato """ O1 = P(0,0) O2 = P(100,100) c1 = Path( circle(O1, 100, minpts= 100, maxerr=100)[0:50] ) c2 = Path( c1[0:25] ) * rot(random.random()*2*math.pi) * xlate(O2) mr = pathOverPath(c2, c1, 10) self.assertAlmostEquals(mr.sse/mr.n, 0.0) A = X(mr.A) for p2 in c2: # controllo della distanza punto-path self.assertAlmostEquals(c1.project(p2*A).dist, 0, tollerance)
def testReverseQuarterCircleOverHalfCircle(self): """ La funzione deve sovrapporre un quarto-di-cerchio a un semicerchio, dopo che il quarto-di-cerchio e' stato traslato e ruotato. L'ordine dei punti del quarto di cerchio viene invertito (senso antiorario). Si nota che l'errore di accostamento dei due path sale a causa di una non piu' perfetta corrispondenza nella scelta dei punti (dato che i path "cominciano" da due punti differenti) """ O1 = P(0,0) O2 = P(100,100) c1 = Path( circle(O1, 100, minpts= 100, maxerr=100)[0:50] ) tmp = c1[25:50] tmp.reverse() c2 = Path( tmp ) * rot(random.random()*2*math.pi) * xlate(O2) # Per passare questo test e' necessario alzare la risoluzione a 1mm mr = pathOverPath(c2, c1, 1) self.assertAlmostEquals(mr.sse/mr.n, 0.0) A = X(mr.A) for p2 in c2: # controllo della distanza punto-path # Per passare questo test e' necessario alzare la tolleranza self.assertAlmostEquals(c1.project(p2*A).dist,0, tollerance)