python-opencv warpPerspective warpAffine 通過透視變化、仿射變化中的矩陣計算變換後的點

1、warpPerspective warpAffine的使用

詳見：https://blog.csdn.net/qq_27261889/article/details/80720359

2、原理分析

主要是得到一個3*3的矩陣M，在計算的時候，通過目标點去找原圖的點，進行采樣

官方文檔可見：https://docs.opencv.org/2.4/modules/imgproc/doc/geometric_transformations.html#warpperspective

通過Opencv自帶的warpPerspective warpAffine很容易對一副圖檔進行變換。但是為了可視化，你也許想知道某個點通過了這種變化後到了那兒。由原理分析可見，通過解上述方程即可。

3、代碼

def reCalculateBBS(self,BBS,M):
        M = np.transpose(M)
        for i in range(len(BBS)):
            k1=BBS[i]
            # x=((M[1][1]-M[2][1]*y0)*(M[2][2]*x0-M[0][2])-
            # (M[0][1]-M[2][1]*x0)*(M[2][2]*y0-M[1][2]))/((M[1][1]-M[2][1]*y0)*(M[0][0]-M[2][0]*x0)-
            # (M[0][1]-M[2][1]*x0)*(M[1][0]-M[2][0]*y0))

            # y=((M[1][0]-M[2][0]*y0)*(M[2][2]*x0-M[0][2])-
            # (M[0][0]-M[2][0]*x0)*(M[2][2]*y0-M[1][2]))/((M[1][0]-M[2][0]*y0)*(M[0][1]-M[2][1]*x0)-
            # (M[0][0]-M[2][0]*x0)*(M[1][1]-M[2][1]*y0))
            
            x = (k1[0]*M[0][0]+k1[1]*M[1][0]+M[2][0])/(k1[0]*M[0][2]+k1[1]*M[1][2]+M[2][2])
            y = (k1[0]*M[0][1]+k1[1]*M[1][1]+M[2][1])/(k1[0]*M[0][2]+k1[1]*M[1][2]+M[2][2])
            # print(k1[0],k1[1],"to",x,y)
            BBS[i]=(x,y)

        return BBS
    
    def reCalculateRotatedBBS(self,BBS,M):
        # M = np.transpose(M)
        for i in range(len(BBS)):
            x0,y0=BBS[i]
            x = M[0][0]*x0+M[0][1]*y0+M[0][2]
            y = M[1][0]*x0+M[1][1]*y0+M[1][2]
            BBS[i]=(x,y)

        return BBS