3dilize_anything / dust3r /inference.py
yansong1616's picture
Upload 90 files
56cd6b7 verified
# Copyright (C) 2024-present Naver Corporation. All rights reserved.
# Licensed under CC BY-NC-SA 4.0 (non-commercial use only).
#
# --------------------------------------------------------
# utilities needed for the inference
# --------------------------------------------------------
import tqdm
import torch
from dust3r.utils.device import to_cpu, collate_with_cat
from dust3r.model import AsymmetricCroCo3DStereo, inf # noqa: F401, needed when loading the model
from dust3r.utils.misc import invalid_to_nans
from dust3r.utils.geometry import depthmap_to_pts3d, geotrf
def load_model(model_path, device):
print('... loading model from', model_path)
ckpt = torch.load(model_path, map_location='cpu')
args = ckpt['args'].model.replace("ManyAR_PatchEmbed", "PatchEmbedDust3R")
if 'landscape_only' not in args:
args = args[:-1] + ', landscape_only=False)'
else:
args = args.replace(" ", "").replace('landscape_only=True', 'landscape_only=False')
assert "landscape_only=False" in args
print(f"instantiating : {args}")
net = eval(args)
print(net.load_state_dict(ckpt['model'], strict=False))
return net.to(device)
def _interleave_imgs(img1, img2):
res = {}
for key, value1 in img1.items():
value2 = img2[key]
if isinstance(value1, torch.Tensor):
value = torch.stack((value1, value2), dim=1).flatten(0, 1)
else:
value = [x for pair in zip(value1, value2) for x in pair]
res[key] = value
return res
def make_batch_symmetric(batch):
view1, view2 = batch
view1, view2 = (_interleave_imgs(view1, view2), _interleave_imgs(view2, view1))
return view1, view2
def loss_of_one_batch(batch, model, criterion, device, symmetrize_batch=False, use_amp=False, ret=None):
view1, view2 = batch # 输入模型的两张图片
for view in batch: # 将输入的图片放到GPU上
for name in 'img pts3d valid_mask camera_pose camera_intrinsics F_matrix corres'.split(): # pseudo_focal
if name not in view:
continue
view[name] = view[name].to(device, non_blocking=True) # 放到GPU上
if symmetrize_batch:
view1, view2 = make_batch_symmetric(batch)
with torch.cuda.amp.autocast(enabled=bool(use_amp)):
pred1, pred2 = model(view1, view2) # model:AsymmetricCroCo3DStereo
# loss is supposed to be symmetric
with torch.cuda.amp.autocast(enabled=False):# loss = None
loss = criterion(view1, view2, pred1, pred2) if criterion is not None else None
result = dict(view1=view1, view2=view2, pred1=pred1, pred2=pred2, loss=loss) #这里loss为None
return result[ret] if ret else result
@torch.no_grad()
def inference(pairs, model, device, batch_size=8):
print(f'>> Inference with model on {len(pairs)} image pairs') # 所有照片两两成一对
result = []
# first, check if all images have the same size
multiple_shapes = not (check_if_same_size(pairs))
if multiple_shapes: # force bs=1
batch_size = 1
for i in tqdm.trange(0, len(pairs), batch_size): # 将所有的pairs依次输入模型
res = loss_of_one_batch(collate_with_cat(pairs[i:i+batch_size]), model, None, device)
result.append(to_cpu(res))
result = collate_with_cat(result, lists=multiple_shapes) # view1、view2分别表示输入模型的两张图片
torch.cuda.empty_cache()
return result
def check_if_same_size(pairs):
shapes1 = [img1['img'].shape[-2:] for img1, img2 in pairs]
shapes2 = [img2['img'].shape[-2:] for img1, img2 in pairs]
return all(shapes1[0] == s for s in shapes1) and all(shapes2[0] == s for s in shapes2)
def get_pred_pts3d(gt, pred, use_pose=False):
if 'depth' in pred and 'pseudo_focal' in pred:
try:
pp = gt['camera_intrinsics'][..., :2, 2]
except KeyError:
pp = None
pts3d = depthmap_to_pts3d(**pred, pp=pp)
elif 'pts3d' in pred:
# pts3d from my camera
pts3d = pred['pts3d']
elif 'pts3d_in_other_view' in pred:
# pts3d from the other camera, already transformed
assert use_pose is True
return pred['pts3d_in_other_view'] # return!
if use_pose:
camera_pose = pred.get('camera_pose')
assert camera_pose is not None
pts3d = geotrf(camera_pose, pts3d)
return pts3d
def find_opt_scaling(gt_pts1, gt_pts2, pr_pts1, pr_pts2=None, fit_mode='weiszfeld_stop_grad', valid1=None, valid2=None):
assert gt_pts1.ndim == pr_pts1.ndim == 4
assert gt_pts1.shape == pr_pts1.shape
if gt_pts2 is not None:
assert gt_pts2.ndim == pr_pts2.ndim == 4
assert gt_pts2.shape == pr_pts2.shape
# concat the pointcloud
nan_gt_pts1 = invalid_to_nans(gt_pts1, valid1).flatten(1, 2)
nan_gt_pts2 = invalid_to_nans(gt_pts2, valid2).flatten(1, 2) if gt_pts2 is not None else None
pr_pts1 = invalid_to_nans(pr_pts1, valid1).flatten(1, 2)
pr_pts2 = invalid_to_nans(pr_pts2, valid2).flatten(1, 2) if pr_pts2 is not None else None
all_gt = torch.cat((nan_gt_pts1, nan_gt_pts2), dim=1) if gt_pts2 is not None else nan_gt_pts1
all_pr = torch.cat((pr_pts1, pr_pts2), dim=1) if pr_pts2 is not None else pr_pts1
dot_gt_pr = (all_pr * all_gt).sum(dim=-1)
dot_gt_gt = all_gt.square().sum(dim=-1)
if fit_mode.startswith('avg'):
# scaling = (all_pr / all_gt).view(B, -1).mean(dim=1)
scaling = dot_gt_pr.nanmean(dim=1) / dot_gt_gt.nanmean(dim=1)
elif fit_mode.startswith('median'):
scaling = (dot_gt_pr / dot_gt_gt).nanmedian(dim=1).values
elif fit_mode.startswith('weiszfeld'):
# init scaling with l2 closed form
scaling = dot_gt_pr.nanmean(dim=1) / dot_gt_gt.nanmean(dim=1)
# iterative re-weighted least-squares
for iter in range(10):
# re-weighting by inverse of distance
dis = (all_pr - scaling.view(-1, 1, 1) * all_gt).norm(dim=-1)
# print(dis.nanmean(-1))
w = dis.clip_(min=1e-8).reciprocal()
# update the scaling with the new weights
scaling = (w * dot_gt_pr).nanmean(dim=1) / (w * dot_gt_gt).nanmean(dim=1)
else:
raise ValueError(f'bad {fit_mode=}')
if fit_mode.endswith('stop_grad'):
scaling = scaling.detach()
scaling = scaling.clip(min=1e-3)
# assert scaling.isfinite().all(), bb()
return scaling