def apply(self, x, y):
scale_matrix = np.array([[self.sxx, self.sxy], [self.sxy, self.syy]])
rot_matrix = rot2deg(self.rot_deg)
ha, dec = xy2hadec(x, y, 0, 0)
x1t, y1t = hadec2xy(ha, dec, self.dha, self.ddec)
xy = scale_matrix.dot(rot_matrix.dot(np.array([x1t, y1t])))
return xy[0], xy[1]
def apply_inverse(self, x, y):
det = self.sxx * self.syy - self.sxy**2
scale_matrix = np.array([[self.syy, -self.sxy], [-self.sxy, self.sxx]
]) / det
rot_matrix = rot2deg(-self.rot_deg)
ha, dec = xy2hadec(x, y, 0, 0)
x1t, y1t = hadec2xy(ha, dec, self.dha, self.ddec)
xy = rot_matrix.dot(scale_matrix.dot(np.array([x1t, y1t])))
return xy[0], xy[1]
def apply_inverse(self,x,y) :
det = self.sxx*self.syy - self.sxy**2
scale_matrix = np.array([[self.syy,-self.sxy],[-self.sxy,self.sxx]])/det
rot_matrix = rot2deg(-self.rot_deg)
xy=rot_matrix.dot(scale_matrix.dot(np.array([x-self.dx,y-self.dy])))
return xy[0],xy[1]
def apply(self,x,y) :
scale_matrix = np.array([[self.sxx,self.sxy],[self.sxy,self.syy]])
rot_matrix = rot2deg(self.rot_deg)
xy=scale_matrix.dot(rot_matrix.dot(np.array([x,y])))
return xy[0]+self.dx,xy[1]+self.dy