Spaces:
Sleeping
Sleeping
| from torch import nn | |
| import torch | |
| import torch.nn.functional as F | |
| from src.facerender.sync_batchnorm import SynchronizedBatchNorm2d as BatchNorm2d | |
| from src.facerender.modules.util import KPHourglass, make_coordinate_grid, AntiAliasInterpolation2d, ResBottleneck | |
| class KPDetector(nn.Module): | |
| """ | |
| Detecting canonical keypoints. Return keypoint position and jacobian near each keypoint. | |
| """ | |
| def __init__(self, block_expansion, feature_channel, num_kp, image_channel, max_features, reshape_channel, reshape_depth, | |
| num_blocks, temperature, estimate_jacobian=False, scale_factor=1, single_jacobian_map=False): | |
| super(KPDetector, self).__init__() | |
| self.predictor = KPHourglass(block_expansion, in_features=image_channel, | |
| max_features=max_features, reshape_features=reshape_channel, reshape_depth=reshape_depth, num_blocks=num_blocks) | |
| # self.kp = nn.Conv3d(in_channels=self.predictor.out_filters, out_channels=num_kp, kernel_size=7, padding=3) | |
| self.kp = nn.Conv3d(in_channels=self.predictor.out_filters, out_channels=num_kp, kernel_size=3, padding=1) | |
| if estimate_jacobian: | |
| self.num_jacobian_maps = 1 if single_jacobian_map else num_kp | |
| # self.jacobian = nn.Conv3d(in_channels=self.predictor.out_filters, out_channels=9 * self.num_jacobian_maps, kernel_size=7, padding=3) | |
| self.jacobian = nn.Conv3d(in_channels=self.predictor.out_filters, out_channels=9 * self.num_jacobian_maps, kernel_size=3, padding=1) | |
| ''' | |
| initial as: | |
| [[1 0 0] | |
| [0 1 0] | |
| [0 0 1]] | |
| ''' | |
| self.jacobian.weight.data.zero_() | |
| self.jacobian.bias.data.copy_(torch.tensor([1, 0, 0, 0, 1, 0, 0, 0, 1] * self.num_jacobian_maps, dtype=torch.float)) | |
| else: | |
| self.jacobian = None | |
| self.temperature = temperature | |
| self.scale_factor = scale_factor | |
| if self.scale_factor != 1: | |
| self.down = AntiAliasInterpolation2d(image_channel, self.scale_factor) | |
| def gaussian2kp(self, heatmap): | |
| """ | |
| Extract the mean from a heatmap | |
| """ | |
| shape = heatmap.shape | |
| heatmap = heatmap.unsqueeze(-1) | |
| grid = make_coordinate_grid(shape[2:], heatmap.type()).unsqueeze_(0).unsqueeze_(0) | |
| value = (heatmap * grid).sum(dim=(2, 3, 4)) | |
| kp = {'value': value} | |
| return kp | |
| def forward(self, x): | |
| if self.scale_factor != 1: | |
| x = self.down(x) | |
| feature_map = self.predictor(x) | |
| prediction = self.kp(feature_map) | |
| final_shape = prediction.shape | |
| heatmap = prediction.view(final_shape[0], final_shape[1], -1) | |
| heatmap = F.softmax(heatmap / self.temperature, dim=2) | |
| heatmap = heatmap.view(*final_shape) | |
| out = self.gaussian2kp(heatmap) | |
| if self.jacobian is not None: | |
| jacobian_map = self.jacobian(feature_map) | |
| jacobian_map = jacobian_map.reshape(final_shape[0], self.num_jacobian_maps, 9, final_shape[2], | |
| final_shape[3], final_shape[4]) | |
| heatmap = heatmap.unsqueeze(2) | |
| jacobian = heatmap * jacobian_map | |
| jacobian = jacobian.view(final_shape[0], final_shape[1], 9, -1) | |
| jacobian = jacobian.sum(dim=-1) | |
| jacobian = jacobian.view(jacobian.shape[0], jacobian.shape[1], 3, 3) | |
| out['jacobian'] = jacobian | |
| return out | |
| class HEEstimator(nn.Module): | |
| """ | |
| Estimating head pose and expression. | |
| """ | |
| def __init__(self, block_expansion, feature_channel, num_kp, image_channel, max_features, num_bins=66, estimate_jacobian=True): | |
| super(HEEstimator, self).__init__() | |
| self.conv1 = nn.Conv2d(in_channels=image_channel, out_channels=block_expansion, kernel_size=7, padding=3, stride=2) | |
| self.norm1 = BatchNorm2d(block_expansion, affine=True) | |
| self.maxpool = nn.MaxPool2d(kernel_size=3, stride=2, padding=1) | |
| self.conv2 = nn.Conv2d(in_channels=block_expansion, out_channels=256, kernel_size=1) | |
| self.norm2 = BatchNorm2d(256, affine=True) | |
| self.block1 = nn.Sequential() | |
| for i in range(3): | |
| self.block1.add_module('b1_'+ str(i), ResBottleneck(in_features=256, stride=1)) | |
| self.conv3 = nn.Conv2d(in_channels=256, out_channels=512, kernel_size=1) | |
| self.norm3 = BatchNorm2d(512, affine=True) | |
| self.block2 = ResBottleneck(in_features=512, stride=2) | |
| self.block3 = nn.Sequential() | |
| for i in range(3): | |
| self.block3.add_module('b3_'+ str(i), ResBottleneck(in_features=512, stride=1)) | |
| self.conv4 = nn.Conv2d(in_channels=512, out_channels=1024, kernel_size=1) | |
| self.norm4 = BatchNorm2d(1024, affine=True) | |
| self.block4 = ResBottleneck(in_features=1024, stride=2) | |
| self.block5 = nn.Sequential() | |
| for i in range(5): | |
| self.block5.add_module('b5_'+ str(i), ResBottleneck(in_features=1024, stride=1)) | |
| self.conv5 = nn.Conv2d(in_channels=1024, out_channels=2048, kernel_size=1) | |
| self.norm5 = BatchNorm2d(2048, affine=True) | |
| self.block6 = ResBottleneck(in_features=2048, stride=2) | |
| self.block7 = nn.Sequential() | |
| for i in range(2): | |
| self.block7.add_module('b7_'+ str(i), ResBottleneck(in_features=2048, stride=1)) | |
| self.fc_roll = nn.Linear(2048, num_bins) | |
| self.fc_pitch = nn.Linear(2048, num_bins) | |
| self.fc_yaw = nn.Linear(2048, num_bins) | |
| self.fc_t = nn.Linear(2048, 3) | |
| self.fc_exp = nn.Linear(2048, 3*num_kp) | |
| def forward(self, x): | |
| out = self.conv1(x) | |
| out = self.norm1(out) | |
| out = F.relu(out) | |
| out = self.maxpool(out) | |
| out = self.conv2(out) | |
| out = self.norm2(out) | |
| out = F.relu(out) | |
| out = self.block1(out) | |
| out = self.conv3(out) | |
| out = self.norm3(out) | |
| out = F.relu(out) | |
| out = self.block2(out) | |
| out = self.block3(out) | |
| out = self.conv4(out) | |
| out = self.norm4(out) | |
| out = F.relu(out) | |
| out = self.block4(out) | |
| out = self.block5(out) | |
| out = self.conv5(out) | |
| out = self.norm5(out) | |
| out = F.relu(out) | |
| out = self.block6(out) | |
| out = self.block7(out) | |
| out = F.adaptive_avg_pool2d(out, 1) | |
| out = out.view(out.shape[0], -1) | |
| yaw = self.fc_roll(out) | |
| pitch = self.fc_pitch(out) | |
| roll = self.fc_yaw(out) | |
| t = self.fc_t(out) | |
| exp = self.fc_exp(out) | |
| return {'yaw': yaw, 'pitch': pitch, 'roll': roll, 't': t, 'exp': exp} | |