[{"createTime":1735734952000,"id":1,"img":"hwy_ms_500_252.jpeg","link":"https://activity.huaweicloud.com/cps.html?fromacct=261f35b6-af54-4511-a2ca-910fa15905d1&utm_source=V1g3MDY4NTY=&utm_medium=cps&utm_campaign=201905","name":"华为云秒杀","status":9,"txt":"华为云38元秒杀","type":1,"updateTime":1735747411000,"userId":3},{"createTime":1736173885000,"id":2,"img":"txy_480_300.png","link":"https://cloud.tencent.com/act/cps/redirect?redirect=1077&cps_key=edb15096bfff75effaaa8c8bb66138bd&from=console","name":"腾讯云秒杀","status":9,"txt":"腾讯云限量秒杀","type":1,"updateTime":1736173885000,"userId":3},{"createTime":1736177492000,"id":3,"img":"aly_251_140.png","link":"https://www.aliyun.com/minisite/goods?userCode=pwp8kmv3","memo":"","name":"阿里云","status":9,"txt":"阿里云2折起","type":1,"updateTime":1736177492000,"userId":3},{"createTime":1735660800000,"id":4,"img":"vultr_560_300.png","link":"https://www.vultr.com/?ref=9603742-8H","name":"Vultr","status":9,"txt":"Vultr送$100","type":1,"updateTime":1735660800000,"userId":3},{"createTime":1735660800000,"id":5,"img":"jdy_663_320.jpg","link":"https://3.cn/2ay1-e5t","name":"京东云","status":9,"txt":"京东云特惠专区","type":1,"updateTime":1735660800000,"userId":3},{"createTime":1735660800000,"id":6,"img":"new_ads.png","link":"https://www.iodraw.com/ads","name":"发布广告","status":9,"txt":"发布广告","type":1,"updateTime":1735660800000,"userId":3},{"createTime":1735660800000,"id":7,"img":"yun_910_50.png","link":"https://activity.huaweicloud.com/discount_area_v5/index.html?fromacct=261f35b6-af54-4511-a2ca-910fa15905d1&utm_source=aXhpYW95YW5nOA===&utm_medium=cps&utm_campaign=201905","name":"底部","status":9,"txt":"高性能云服务器2折起","type":2,"updateTime":1735660800000,"userId":3}]
推导
三角形ABC仿射成为三角形DEF的变换矩阵M
猜测矩阵M=
[
[a1,b1,c1],
[a2,b2,c2]
]
仿射变换的数学联系
对于A点和D点
AX*a1+AY*b1+c1=DX
AX*a2+AY*b2+c1=DY
对于B点和E点
BX*a1+BY*b1+c1=EX
BX*a2+BY*b2+c2=EY
对于C点和F点
CX*a1+CY*b1+c1=FX
CX*a2+CY*b2+c2=FY
求解
对以上数学联系式均取第一条
AX*a1+AY*b1+c1=DX
BX*a1+BY*b1+c1=EX
CX*a1+CY*b1+c1=FX
可以看出这是一个三元一次方程组,我们可以借助scipy库linalg模块的solve函数进行求解,它的用法大家可以参考一下其他资料,这里不多赘述
对数学联系式均取第二条
AX*a2+AY*b2+c1=DY
BX*a2+BY*b2+c2=EY
CX*a2+CY*b2+c2=FY
也是一个三元一次方程组,同样可以借助scipy库linalg模块的solve函数进行求解
验证
根据以上各个点的坐标,分别用opencv的函数,和我们猜想的计算过程,观测对比两个结果是否一致
import numpy as np import cv2 as cv import math AX = 4 AY = 2 BX = 9 BY = 2 CX
= 4 CY = 7 DX = 13 DY = 8 EX = 13 EY = 6 FX = 11 FY = 8 ps1 = np.float32([ [AX,
AY], [BX, BY], [CX, CY] ]) ps2 = np.float32([ [DX, DY], [EX, EY], [FX, FY] ]) M
= cv.getAffineTransform(ps1, ps2) print(f"opencv_M=\n{M}") print() A =
np.array([ [AX, AY, 1], [BX, BY, 1], [CX, CY, 1] ]) b1 = np.array([ [DX], [EX],
[FX] ]) b2 = np.array([ [DY], [EY], [FY] ]) from scipy import linalg r1 =
linalg.solve(A, b1) r2 = linalg.solve(A, b2) r = np.array([ r1.T[0], r2.T[0] ])
print(f"self_M=\n{r}")
运行结果
可以看出,结果是一致的,函数的底层实现原理就是如此