faceswap / face_swapper.py
coxapi's picture
Upload 31 files
ba9144c verified
import time
import torch
import onnx
import cv2
import onnxruntime
import numpy as np
from tqdm import tqdm
import torch.nn as nn
from onnx import numpy_helper
from skimage import transform as trans
import torchvision.transforms.functional as F
import torch.nn.functional as F
from utils import mask_crop, laplacian_blending
arcface_dst = np.array(
[[38.2946, 51.6963], [73.5318, 51.5014], [56.0252, 71.7366],
[41.5493, 92.3655], [70.7299, 92.2041]],
dtype=np.float32)
def estimate_norm(lmk, image_size=112, mode='arcface'):
assert lmk.shape == (5, 2)
assert image_size % 112 == 0 or image_size % 128 == 0
if image_size % 112 == 0:
ratio = float(image_size) / 112.0
diff_x = 0
else:
ratio = float(image_size) / 128.0
diff_x = 8.0 * ratio
dst = arcface_dst * ratio
dst[:, 0] += diff_x
tform = trans.SimilarityTransform()
tform.estimate(lmk, dst)
M = tform.params[0:2, :]
return M
def norm_crop2(img, landmark, image_size=112, mode='arcface'):
M = estimate_norm(landmark, image_size, mode)
warped = cv2.warpAffine(img, M, (image_size, image_size), borderValue=0.0)
return warped, M
class Inswapper():
def __init__(self, model_file=None, batch_size=32, providers=['CPUExecutionProvider']):
self.model_file = model_file
self.batch_size = batch_size
model = onnx.load(self.model_file)
graph = model.graph
self.emap = numpy_helper.to_array(graph.initializer[-1])
self.session_options = onnxruntime.SessionOptions()
self.session = onnxruntime.InferenceSession(self.model_file, sess_options=self.session_options, providers=providers)
def forward(self, imgs, latents):
preds = []
for img, latent in zip(imgs, latents):
img = img / 255
pred = self.session.run(['output'], {'target': img, 'source': latent})[0]
preds.append(pred)
def get(self, imgs, target_faces, source_faces):
imgs = list(imgs)
preds = [None] * len(imgs)
matrs = [None] * len(imgs)
for idx, (img, target_face, source_face) in enumerate(zip(imgs, target_faces, source_faces)):
matrix, blob, latent = self.prepare_data(img, target_face, source_face)
pred = self.session.run(['output'], {'target': blob, 'source': latent})[0]
pred = pred.transpose((0, 2, 3, 1))[0]
pred = np.clip(255 * pred, 0, 255).astype(np.uint8)[:, :, ::-1]
preds[idx] = pred
matrs[idx] = matrix
return (preds, matrs)
def prepare_data(self, img, target_face, source_face):
if isinstance(img, str):
img = cv2.imread(img)
aligned_img, matrix = norm_crop2(img, target_face.kps, 128)
blob = cv2.dnn.blobFromImage(aligned_img, 1.0 / 255, (128, 128), (0., 0., 0.), swapRB=True)
latent = source_face.normed_embedding.reshape((1, -1))
latent = np.dot(latent, self.emap)
latent /= np.linalg.norm(latent)
return (matrix, blob, latent)
def batch_forward(self, img_list, target_f_list, source_f_list):
num_samples = len(img_list)
num_batches = (num_samples + self.batch_size - 1) // self.batch_size
for i in tqdm(range(num_batches), desc="Generating face"):
start_idx = i * self.batch_size
end_idx = min((i + 1) * self.batch_size, num_samples)
batch_img = img_list[start_idx:end_idx]
batch_target_f = target_f_list[start_idx:end_idx]
batch_source_f = source_f_list[start_idx:end_idx]
batch_pred, batch_matr = self.get(batch_img, batch_target_f, batch_source_f)
yield batch_pred, batch_matr
def paste_to_whole(foreground, background, matrix, mask=None, crop_mask=(0,0,0,0), blur_amount=0.1, erode_amount = 0.15, blend_method='linear'):
inv_matrix = cv2.invertAffineTransform(matrix)
fg_shape = foreground.shape[:2]
bg_shape = (background.shape[1], background.shape[0])
foreground = cv2.warpAffine(foreground, inv_matrix, bg_shape, borderValue=0.0)
if mask is None:
mask = np.full(fg_shape, 1., dtype=np.float32)
mask = mask_crop(mask, crop_mask)
mask = cv2.warpAffine(mask, inv_matrix, bg_shape, borderValue=0.0)
else:
assert fg_shape == mask.shape[:2], "foreground & mask shape mismatch!"
mask = mask_crop(mask, crop_mask).astype('float32')
mask = cv2.warpAffine(mask, inv_matrix, (background.shape[1], background.shape[0]), borderValue=0.0)
_mask = mask.copy()
_mask[_mask > 0.05] = 1.
non_zero_points = cv2.findNonZero(_mask)
_, _, w, h = cv2.boundingRect(non_zero_points)
mask_size = int(np.sqrt(w * h))
if erode_amount > 0:
kernel_size = max(int(mask_size * erode_amount), 1)
structuring_element = cv2.getStructuringElement(cv2.MORPH_RECT, (kernel_size, kernel_size))
mask = cv2.erode(mask, structuring_element)
if blur_amount > 0:
kernel_size = max(int(mask_size * blur_amount), 3)
if kernel_size % 2 == 0:
kernel_size += 1
mask = cv2.GaussianBlur(mask, (kernel_size, kernel_size), 0)
mask = np.tile(np.expand_dims(mask, axis=-1), (1, 1, 3))
if blend_method == 'laplacian':
composite_image = laplacian_blending(foreground, background, mask.clip(0,1), num_levels=4)
else:
composite_image = mask * foreground + (1 - mask) * background
return composite_image.astype("uint8").clip(0, 255)