Pude hacer esto con una transformación afín (gracias a this question). Después de la transformación afín, el triángulo de destino se dibuja en una máscara y luego se pega a la imagen de destino. Esto es lo que ocurrió:
import Image
import ImageDraw
import numpy
def transformblit(src_tri, dst_tri, src_img, dst_img):
((x11,x12), (x21,x22), (x31,x32)) = src_tri
((y11,y12), (y21,y22), (y31,y32)) = dst_tri
M = numpy.array([
[y11, y12, 1, 0, 0, 0],
[y21, y22, 1, 0, 0, 0],
[y31, y32, 1, 0, 0, 0],
[0, 0, 0, y11, y12, 1],
[0, 0, 0, y21, y22, 1],
[0, 0, 0, y31, y32, 1]
])
y = numpy.array([x11, x21, x31, x12, x22, x32])
A = numpy.linalg.solve(M, y)
src_copy = src_img.copy()
srcdraw = ImageDraw.Draw(src_copy)
srcdraw.polygon(src_tri)
src_copy.show()
transformed = src_img.transform(dst_img.size, Image.AFFINE, A)
mask = Image.new('1', dst_img.size)
maskdraw = ImageDraw.Draw(mask)
maskdraw.polygon(dst_tri, fill=255)
dstdraw = ImageDraw.Draw(dst_img)
dstdraw.polygon(dst_tri, fill=(255,255,255))
dst_img.show()
dst_img.paste(transformed, mask=mask)
dst_img.show()
im100 = Image.open('test100.jpg')
im250 = Image.open('test250.jpg')
tri1 = [(10,10), (20,20), (10,20)]
tri2 = [(35,30), (75,19), (50,90)]
transformblit(tri1, tri2, im100, im250)
El origen de la imagen de 100x100 se ve así (triángulo sobrepuesto en blanco):
El destino de imagen de 250x250 se ve así (región triangular rellena con blanco):
y luego, después de la transformación y pegar, el des tino imagen se ve así:
¿No deseaba calcular la transformación a mano? – carlosdc
Sí, pero no parecía que hubiera una manera de hacer eso con PIL.No fue demasiado problema ya que obtuve la fórmula de la publicación referenciada y numpy resuelve el sistema de ecuaciones para mí. Lo que realmente quería evitar era copiar píxel por píxel. – jterrace