作者 | 李秋键

责编 | 李雪敬

头图 | CSDN 下载自视觉中国

出品 | AI科技大本营（ID：rgznai100）

近段时间，一个让蒙娜丽莎图像动起来的项目火遍了朋友圈。而今天我们就将实现让图片中的人物随着视频人物一起产生动作。

其中通过在静止图像中动画对象产生视频有无数的应用跨越的领域兴趣，包括电影制作、摄影和电子商务。更准确地说，是图像动画指将提取的视频外观结合起来自动合成视频的任务一种源图像与运动模式派生的视频。

近年来，深度生成模型作为一种有效的图像动画技术出现了视频重定向。特别是，可生成的对抗网络(GANS)和变分自动编码器(VAES)已被用于在视频中人类受试者之间转换面部表情或运动模式。

根据论文FirstOrder Motion Model for Image Animation可知，在姿态迁移的大任务当中，Monkey-Net首先尝试了通过自监督范式预测关键点来表征姿态信息，测试阶段估计驱动视频的姿态关键点完成迁移工作。在此基础上，FOMM使用了相邻关键点的局部仿射变换来模拟物体运动，还额外考虑了遮挡的部分，遮挡的部分可以使用image inpainting生成。

而今天我们就将借助论文所分享的源代码，构建模型创建自己需要的人物运动。具体流程如下。

实验前的准备

首先我们使用的python版本是3.6.5所用到的模块如下：

• imageio模块用来控制图像的输入输出等。

• Matplotlib模块用来绘图。

• numpy模块用来处理矩阵运算。

• Pillow库用来加载数据处理。

• pytorch模块用来创建模型和模型训练等。

• 完整模块需求参见requirements.txt文件。

模型的加载和调用

通过定义命令行参数来达到加载模型，图片等目的。

（1）首先是训练模型的读取，包括模型加载方式：

def load_checkpoints(config_path, checkpoint_path, cpu=False):with open(config_path) as f:config = yaml.load(f)generator = OcclusionAwareGenerator(**config['model_params']['generator_params'],**config['model_params']['common_params'])if not cpu:generator.cuda()kp_detector = KPDetector(**config['model_params']['kp_detector_params'],**config['model_params']['common_params'])if not cpu:kp_detector.cuda()if cpu:checkpoint = torch.load(checkpoint_path, map_location=torch.device('cpu'))else:checkpoint = torch.load(checkpoint_path)generator.load_state_dict(checkpoint['generator'])kp_detector.load_state_dict(checkpoint['kp_detector'])if not cpu:generator = DataParallelWithCallback(generator)kp_detector = DataParallelWithCallback(kp_detector)generator.eval()kp_detector.eval()return generator, kp_detector

（2）然后是利用模型创建产生的虚拟图像，找到最佳的脸部特征：

def make_animation(source_image, driving_video, generator, kp_detector, relative=True, adapt_movement_scale=True, cpu=False):with torch.no_grad():predictions = []source = torch.tensor(source_image[np.newaxis].astype(np.float32)).permute(0, 3, 1, 2)if not cpu:source = source.cuda()driving = torch.tensor(np.array(driving_video)[np.newaxis].astype(np.float32)).permute(0, 4, 1, 2, 3)kp_source = kp_detector(source)kp_driving_initial = kp_detector(driving[:, :, 0])for frame_idx in tqdm(range(driving.shape[2])):driving_frame = driving[:, :, frame_idx]if not cpu:driving_frame = driving_frame.cuda()kp_driving = kp_detector(driving_frame)kp_norm = normalize_kp(kp_source=kp_source, kp_driving=kp_driving,kp_driving_initial=kp_driving_initial, use_relative_movement=relative,use_relative_jacobian=relative, adapt_movement_scale=adapt_movement_scale)out = generator(source, kp_source=kp_source, kp_driving=kp_norm)           predictions.append(np.transpose(out['prediction'].data.cpu().numpy(), [0, 2, 3, 1])[0])return predictions
def find_best_frame(source, driving, cpu=False):import face_alignmentdef normalize_kp(kp):kp = kp - kp.mean(axis=0, keepdims=True)area = ConvexHull(kp[:, :2]).volumearea = np.sqrt(area)kp[:, :2] = kp[:, :2] / areareturn kpfa = face_alignment.FaceAlignment(face_alignment.LandmarksType._2D, flip_input=True,device='cpu' if cpu else 'cuda')kp_source = fa.get_landmarks(255 * source)[0]kp_source = normalize_kp(kp_source)norm  = float('inf')frame_num = 0for i, image in tqdm(enumerate(driving)):kp_driving = fa.get_landmarks(255 * image)[0]kp_driving = normalize_kp(kp_driving)new_norm = (np.abs(kp_source - kp_driving) ** 2).sum()if new_norm < norm:norm = new_normframe_num = i
return frame_num

（3）接着定义命令行调用参数加载图片、视频等方式：

parser = ArgumentParser()parser.add_argument("--config", required=True, help="path to config")parser.add_argument("--checkpoint", default='vox-cpk.pth.tar', help="path to checkpoint to restore")parser.add_argument("--source_image", default='sup-mat/source.png', help="path to source image")parser.add_argument("--driving_video", default='sup-mat/source.png', help="path to driving video")parser.add_argument("--result_video", default='result.mp4', help="path to output")parser.add_argument("--relative", dest="relative", action="store_true", help="use relative or absolute keypoint coordinates")parser.add_argument("--adapt_scale", dest="adapt_scale", action="store_true", help="adapt movement scale based on convex hull of keypoints")parser.add_argument("--find_best_frame", dest="find_best_frame", action="store_true", help="Generate from the frame that is the most alligned with source. (Only for faces, requires face_aligment lib)")parser.add_argument("--best_frame", dest="best_frame", type=int, default=None,  help="Set frame to start from.")parser.add_argument("--cpu", dest="cpu", action="store_true", help="cpu mode.")parser.set_defaults(relative=False)parser.set_defaults(adapt_scale=False)opt = parser.parse_args()source_image = imageio.imread(opt.source_image)reader = imageio.get_reader(opt.driving_video)fps = reader.get_meta_data()['fps']driving_video = []try:for im in reader:driving_video.append(im)except RuntimeError:passreader.close()source_image = resize(source_image, (256, 256))[..., :3]driving_video = [resize(frame, (256, 256))[..., :3] for frame in driving_video]generator, kp_detector = load_checkpoints(config_path=opt.config, checkpoint_path=opt.checkpoint, cpu=opt.cpu)if opt.find_best_frame or opt.best_frame is not None:i = opt.best_frame if opt.best_frame is not None else find_best_frame(source_image, driving_video, cpu=opt.cpu)print ("Best frame: " + str(i))driving_forward = driving_video[i:]driving_backward = driving_video[:(i+1)][::-1]predictions_forward = make_animation(source_image, driving_forward, generator, kp_detector, relative=opt.relative, adapt_movement_scale=opt.adapt_scale, cpu=opt.cpu)predictions_backward = make_animation(source_image, driving_backward, generator, kp_detector, relative=opt.relative, adapt_movement_scale=opt.adapt_scale, cpu=opt.cpu)predictions = predictions_backward[::-1] + predictions_forward[1:]else:predictions = make_animation(source_image, driving_video, generator, kp_detector, relative=opt.relative, adapt_movement_scale=opt.adapt_scale, cpu=opt.cpu)
imageio.mimsave(opt.result_video, [img_as_ubyte(frame) for frame in predictions], fps=fps)

模型的搭建

整个模型训练过程是图像重建的过程，输入是源图像和驱动图像，输出是保留源图像物体信息的含有驱动图像姿态的新图像，其中输入的两张图像来源于同一个视频，即同一个物体信息，那么整个训练过程就是驱动图像的重建过程。大体上来说分成两个模块，一个是motion estimation module，另一个是imagegeneration module。

（1）其中通过定义VGG19模型建立网络层作为perceptual损失。

其中手动输入数据进行预测需要设置更多的GUI按钮，其中代码如下：

class Vgg19(torch.nn.Module):"""Vgg19 network for perceptual loss. See Sec 3.3."""def __init__(self, requires_grad=False):super(Vgg19, self).__init__()vgg_pretrained_features = models.vgg19(pretrained=True).featuresself.slice1 = torch.nn.Sequential()self.slice2 = torch.nn.Sequential()self.slice3 = torch.nn.Sequential()self.slice4 = torch.nn.Sequential()self.slice5 = torch.nn.Sequential()for x in range(2):self.slice1.add_module(str(x), vgg_pretrained_features[x])for x in range(2, 7):self.slice2.add_module(str(x), vgg_pretrained_features[x])for x in range(7, 12):self.slice3.add_module(str(x), vgg_pretrained_features[x])for x in range(12, 21):self.slice4.add_module(str(x), vgg_pretrained_features[x])for x in range(21, 30):self.slice5.add_module(str(x), vgg_pretrained_features[x])self.mean = torch.nn.Parameter(data=torch.Tensor(np.array([0.485, 0.456, 0.406]).reshape((1, 3, 1, 1))),requires_grad=False)self.std = torch.nn.Parameter(data=torch.Tensor(np.array([0.229, 0.224, 0.225]).reshape((1, 3, 1, 1))),requires_grad=False)if not requires_grad:for param in self.parameters():param.requires_grad = Falsedef forward(self, X):X = (X - self.mean) / self.stdh_relu1 = self.slice1(X)h_relu2 = self.slice2(h_relu1)h_relu3 = self.slice3(h_relu2)h_relu4 = self.slice4(h_relu3)h_relu5 = self.slice5(h_relu4)out = [h_relu1, h_relu2, h_relu3, h_relu4, h_relu5]return out

（2）创建图像金字塔计算金字塔感知损失：

class ImagePyramide(torch.nn.Module):"""Create image pyramide for computing pyramide perceptual loss. See Sec 3.3"""def __init__(self, scales, num_channels):super(ImagePyramide, self).__init__()downs = {}for scale in scales:downs[str(scale).replace('.', '-')] = AntiAliasInterpolation2d(num_channels, scale)self.downs = nn.ModuleDict(downs)def forward(self, x):out_dict = {}for scale, down_module in self.downs.items():out_dict['prediction_' + str(scale).replace('-', '.')] = down_module(x)return out_dict

（3）等方差约束的随机tps变换

class Transform:"""Random tps transformation for equivariance constraints. See Sec 3.3"""def __init__(self, bs, **kwargs):noise = torch.normal(mean=0, std=kwargs['sigma_affine'] * torch.ones([bs, 2, 3]))self.theta = noise + torch.eye(2, 3).view(1, 2, 3)self.bs = bsif ('sigma_tps' in kwargs) and ('points_tps' in kwargs):self.tps = Trueself.control_points = make_coordinate_grid((kwargs['points_tps'], kwargs['points_tps']), type=noise.type())self.control_points = self.control_points.unsqueeze(0)self.control_params = torch.normal(mean=0,std=kwargs['sigma_tps'] * torch.ones([bs, 1, kwargs['points_tps'] ** 2]))else:self.tps = Falsedef transform_frame(self, frame):grid = make_coordinate_grid(frame.shape[2:], type=frame.type()).unsqueeze(0)grid = grid.view(1, frame.shape[2] * frame.shape[3], 2)grid = self.warp_coordinates(grid).view(self.bs, frame.shape[2], frame.shape[3], 2)return F.grid_sample(frame, grid, padding_mode="reflection")def warp_coordinates(self, coordinates):theta = self.theta.type(coordinates.type())theta = theta.unsqueeze(1)transformed = torch.matmul(theta[:, :, :, :2], coordinates.unsqueeze(-1)) + theta[:, :, :, 2:]transformed = transformed.squeeze(-1)if self.tps:control_points = self.control_points.type(coordinates.type())control_params = self.control_params.type(coordinates.type())distances = coordinates.view(coordinates.shape[0], -1, 1, 2) - control_points.view(1, 1, -1, 2)distances = torch.abs(distances).sum(-1)result = distances ** 2result = result * torch.log(distances + 1e-6)result = result * control_paramsresult = result.sum(dim=2).view(self.bs, coordinates.shape[1], 1)transformed = transformed + resultreturn transformeddef jacobian(self, coordinates):new_coordinates = self.warp_coordinates(coordinates)grad_x = grad(new_coordinates[..., 0].sum(), coordinates, create_graph=True)grad_y = grad(new_coordinates[..., 1].sum(), coordinates, create_graph=True)jacobian = torch.cat([grad_x[0].unsqueeze(-2), grad_y[0].unsqueeze(-2)], dim=-2)return jacobian    

（4）生成器的定义：生成器，给定的源图像和和关键点尝试转换图像根据运动轨迹引起要点。部分代码如下：

class OcclusionAwareGenerator(nn.Module):def __init__(self, num_channels, num_kp, block_expansion, max_features, num_down_blocks,num_bottleneck_blocks, estimate_occlusion_map=False, dense_motion_params=None, estimate_jacobian=False):super(OcclusionAwareGenerator, self).__init__()if dense_motion_params is not None:self.dense_motion_network = DenseMotionNetwork(num_kp=num_kp, num_channels=num_channels,estimate_occlusion_map=estimate_occlusion_map,**dense_motion_params)else:self.dense_motion_network = Noneself.first = SameBlock2d(num_channels, block_expansion, kernel_size=(7, 7), padding=(3, 3))down_blocks = []for i in range(num_down_blocks):in_features = min(max_features, block_expansion * (2 ** i))out_features = min(max_features, block_expansion * (2 ** (i + 1)))down_blocks.append(DownBlock2d(in_features, out_features, kernel_size=(3, 3), padding=(1, 1)))self.down_blocks = nn.ModuleList(down_blocks)up_blocks = []for i in range(num_down_blocks):in_features = min(max_features, block_expansion * (2 ** (num_down_blocks - i)))out_features = min(max_features, block_expansion * (2 ** (num_down_blocks - i - 1)))up_blocks.append(UpBlock2d(in_features, out_features, kernel_size=(3, 3), padding=(1, 1)))self.up_blocks = nn.ModuleList(up_blocks)self.bottleneck = torch.nn.Sequential()in_features = min(max_features, block_expansion * (2 ** num_down_blocks))for i in range(num_bottleneck_blocks):self.bottleneck.add_module('r' + str(i), ResBlock2d(in_features, kernel_size=(3, 3), padding=(1, 1)))self.final = nn.Conv2d(block_expansion, num_channels, kernel_size=(7, 7), padding=(3, 3))self.estimate_occlusion_map = estimate_occlusion_mapself.num_channels = num_channels

（5）判别器类似于Pix2PixGenerator。

def __init__(self, num_channels=3, block_expansion=64, num_blocks=4, max_features=512,sn=False, use_kp=False, num_kp=10, kp_variance=0.01, **kwargs):super(Discriminator, self).__init__()down_blocks = []for i in range(num_blocks):down_blocks.append(DownBlock2d(num_channels + num_kp * use_kp if i == 0 else min(max_features, block_expansion * (2 ** i)),min(max_features, block_expansion * (2 ** (i + 1))),norm=(i != 0), kernel_size=4, pool=(i != num_blocks - 1), sn=sn))self.down_blocks = nn.ModuleList(down_blocks)self.conv = nn.Conv2d(self.down_blocks[-1].conv.out_channels, out_channels=1, kernel_size=1)if sn:self.conv = nn.utils.spectral_norm(self.conv)self.use_kp = use_kpself.kp_variance = kp_variancedef forward(self, x, kp=None):feature_maps = []out = xif self.use_kp:heatmap = kp2gaussian(kp, x.shape[2:], self.kp_variance)out = torch.cat([out, heatmap], dim=1)for down_block in self.down_blocks:feature_maps.append(down_block(out))out = feature_maps[-1]prediction_map = self.conv(out)return feature_maps, prediction_map

最终通过以下代码调用模型训练“python demo.py--config config/vox-adv-256.yaml --driving_video path/to/driving/1.mp4--source_image path/to/source/7.jpg --checkpointpath/to/checkpoint/vox-adv-cpk.pth.tar --relative --adapt_scale”

效果如下：

完整代码：

https://pan.baidu.com/s/1nPE13oI1qOerN0ANQSH92g

提取码：e4kx

作者简介：

李秋键，CSDN博客专家，CSDN达人课作者。硕士在读于中国矿业大学，开发有taptap竞赛获奖等等。

更多精彩推荐
☞谷歌软件工程师薪资百万，大厂薪资有多高？
☞CSDN 创始人蒋涛：选择长沙作“大本营”，打造开发者中心城市
☞杜甫在线演唱《奇迹再现》、兵马俑真人还原……用AI技术打破次元壁的大谷来参加腾讯全球数字生态大会啦！
☞开放源码，华为鸿蒙HarmonyOS 2.0来了
☞20张图，带你搞懂高并发中的线程与线程池！
☞跨链，该怎么跨？

点分享点点赞点在看