Files
2024-07-11 10:48:43 -06:00

81 lines
3.1 KiB
Python

# -*- coding: utf-8 -*-
# @Author: Peng Xiang
import torch
import torch.nn as nn
from .utils import MLP_Res, MLP_CONV
from .skip_transformer import SkipTransformer
from .attention import ResidualTransformerBlock
class SPD_crossattn(nn.Module):
def __init__(self, dim_feat=512, up_factor=2, i=0, radius=1, bounding=True, global_feat=True):
"""Snowflake Point Deconvolution"""
super().__init__()
self.i = i
self.up_factor = up_factor
self.bounding = bounding
self.radius = radius
self.global_feat = global_feat
self.ps_dim = 32 if global_feat else 64
self.mlp_1 = MLP_CONV(in_channel=3, layer_dims=[64, 128])
self.pcd_image_attn = ResidualTransformerBlock(
device=torch.device('cuda'),
dtype=torch.float32,
n_data=128,
width=128,
heads=8,
init_scale=1.0,
)
self.mlp_2 = MLP_CONV(in_channel=128 * 2 + dim_feat if self.global_feat else 128, layer_dims=[256, 128])
self.skip_transformer = SkipTransformer(in_channel=128, dim=64)
self.mlp_ps = MLP_CONV(in_channel=128, layer_dims=[64, self.ps_dim])
self.ps = nn.ConvTranspose1d(self.ps_dim, 128, up_factor, up_factor, bias=False) # point-wise splitting
self.up_sampler = nn.Upsample(scale_factor=up_factor)
self.mlp_delta_feature = MLP_Res(in_dim=256, hidden_dim=128, out_dim=128)
self.mlp_delta = MLP_CONV(in_channel=128, layer_dims=[64, 3])
def forward(self, pcd_prev, feat_global=None, K_prev=None):
"""
Args:
pcd_prev: Tensor, (B, 3, N_prev)
feat_global: Tensor, (B, dim_feat, 1)
K_prev: Tensor, (B, 128, N_prev)
Returns:
pcd_child: Tensor, up sampled point cloud, (B, 3, N_prev * up_factor)
K_curr: Tensor, displacement feature of current step, (B, 128, N_prev * up_factor)
"""
b, _, n_prev = pcd_prev.shape
feat_1 = self.mlp_1(pcd_prev)
# feat_1 = torch.cat([feat_1,
# torch.max(feat_1, 2, keepdim=True)[0].repeat((1, 1, feat_1.size(2))),
# feat_global.repeat(1, 1, feat_1.size(2))], 1) if self.global_feat else feat_1
feat_1 = torch.permute(feat_1, (0, 2, 1))
feat_global = torch.permute(feat_global, (0, 2, 1))
feat_1 = self.pcd_image_attn(feat_1, feat_global)
Q = torch.permute(feat_1, (0, 2, 1))
# Q = self.mlp_2(feat_1)
H = self.skip_transformer(pcd_prev, K_prev if K_prev is not None else Q, Q)
feat_child = self.mlp_ps(H)
feat_child = self.ps(feat_child) # (B, 128, N_prev * up_factor)
H_up = self.up_sampler(H)
K_curr = self.mlp_delta_feature(torch.cat([feat_child, H_up], 1))
delta = self.mlp_delta(torch.relu(K_curr))
if self.bounding:
delta = torch.tanh(delta) / self.radius**self.i # (B, 3, N_prev * up_factor)
pcd_child = self.up_sampler(pcd_prev)
pcd_child = pcd_child + delta
return pcd_child, K_curr