def errorBlend(self,bcs): # Arg is baryCoords of an image point. sumVec = (0.0, 0.0) for bc,errVec in zip(bcs,self.error): sumVec = geom.ptOffset(sumVec, errVec, -bc) # Negate the error. pass return sumVec
lineHead = m1.group(1) data = [float(f) for f in m1.groups()[1:]] # XXX Horrid Hack Warning - We need right-handed coordinates, # but the image Y coordinate goes down from 0 at the top. # Negate it internally. pixLoc = geom.ptBlend((data[0], -data[1]), (data[2], -data[3])) wcLoc = (data[4], data[5]) wAng = data[6] # Find the quadrant by looking at the center edges of the triangles. # We can actually blend linearly past the outer edges... for iTri in range(nPts-1): # Get the barycentric coords of the image point. bcs = triangles[iTri].baryCoords(pixLoc) if bcs[1] >= 0.0 and bcs[2] >= 0.0: # This is the triangle containing this image point. newLoc = geom.ptOffset(wcLoc, triangles[iTri].errorBlend(bcs)) ##print "pixLoc %s, iTri %d, bcs %s"%(pixLoc, iTri, bcs) ##print "triangles[%d] %s"%(iTri, triangles[iTri]) print "%s(%f,%f,%f)"%(lineHead, newLoc[0], newLoc[1], wAng) break pass pass gridLine = gridData.readline() pass
def errorBlend(self, bcs): # Arg is baryCoords of an image point. sumVec = (0.0, 0.0) for bc, errVec in zip(bcs, self.error): sumVec = geom.ptOffset(sumVec, errVec, -bc) # Negate the error. pass return sumVec