墙面也能变镜子，只看影子就能还原视频，MIT新算法让摄像头无死角

边策 十三 发自 凹非寺
量子位 报道 | 公众号 QbitAI

你在看电影，墙上的影子也在动。如果只让你看到这样一段视频，你能猜出来屏幕上播放的是什么吗？

最近MIT人工智能实验室（CSAIL）开发出的算法可以做到：

而真实的视频是这样的：

算法还原的结果只是模糊了些，但已经能猜出视频的大致内容了。

有了这套算法，就可以通过观察视频中阴影和几何图形之间的相互作用，预测出光在场景中的传播方式，然后从观察到的阴影中估计隐藏的视频，甚至看出人的轮廓。

这种图像重建算法会有许多用途：自动驾驶汽车可以了解拐角处正在发生的事情，监控摄像头也可以发现在视线外的人。该论文已经被最近召开的NeurIPS 2019大会所收录。

根据影子，还原看不见的死角

在下面这个场景中，人摆弄玩具的镜头被隐藏，在人们正常视野范围内是无法看见的。

我们唯一能够捕捉到的就是打在墙上的影子。

MIT的这项研究就是仅仅利用这些影子，重新还原隐藏视频的原貌。

算法对场景中的光线传输做了预测。

下图左侧是通过算法估计出来的阴影，而右侧则是实际场景中的阴影。

根据光线传输的预测和估计，就可以重建隐藏物体的运动情况。

例如在下图中，隐藏在摄像头视野之外的人，双手不停的在做着运动。

而我们能够观察到的只是图中左侧单个物体中光影的变化。

就是利用这样简单的光影变化，便可以重构出如图中右侧的视频。

与隐藏视频相比，重构的视频已经可以大致再现双手运动的轮廓。

根据房间杂物乱七八糟的影子，同样也可以还原下图隐藏视频中人物走动的大致轮廓。

当然还有这样的。

以及这样的。

总体来说，MIT的这次研究，能够根据隐藏视频中的内容将光线传输分离出来，从而对它做一个大致的估计。

原理

图像的影子具有线性叠加的特性。如果依次点亮隐藏的屏幕上的两个像素，并墙上的影子图像求和，结果应该和一次同时点亮两个像素时得到的图像相同。

从数学上来看，无论是墙上的影子，还是屏幕上的画面，都是矩阵。而符合线性叠加的特性，等于是在这两个矩阵之间做线性变换。

我们不妨把这二者别看做两个矩阵Z和L，经过空间传输，画面L变成了影子Z，这就相当于做了一次矩阵乘法，T是空间传输矩阵。

Z=TL

问题是我们只看到了影子Z，对于T和L，我们一无所知。

这篇论文的第一作者Miika Aittala说：“这就像是我告诉你，我正在考虑两个秘密数字，它们的乘积为80。你能猜出它们是什么吗？也许是40和2，或是371.8和0.2152。”

对于这个问题也是类似，而且我们在每个像素上都会有一样的困扰。如何求出传输矩阵T成了问题的关键。

如果我们知道了光传输矩阵，那么求原图像L的操作就变成了最小化||Z-TL||2的最小二乘法问题。

作者通过DCT和PCA方法测量了T，然后通过求逆的方法恢复了原始图像。

因此知道了T，接下来恢复图像就好办了。

但是这篇文章要挑战更高的难度：如何在不知道T的情况下恢复图像。他们使用了去年一篇Deep Image Prior论文中的新的矩阵分解方法。

这篇文章曾经被CVPR 2018收录，在inpainting问题上收到了不错的效果。

过去也有一些矩阵分解方法，但是分解得到的矩阵一般都是低秩的，与图像差别很大，而且对初始值和优化的动力学都高度敏感，只能针对特定问题量身定制。

而作者使用的矩阵分解方法里，CNN随机初始化并“过拟合”，将两个噪声矢量映射到两个矩阵T和L，使它们的乘积与输入矩阵Z匹配。此过程将因式分解正则化为更接近于图像的结构。

结合上面的思想

Deep Image Prior

作者首先描述了一种基线方法，在知道Z和L的情况下求T，其实就是求||Z-TL||2最小值的最小二乘法问题。

使用了Deep Image Prior的方法，作者提出了一种在无法测得的光传输矩阵时，逆向求原图像的方法。逆向光传输矩阵的体系架构和数据流动如下图所示：

其中左下角是存储在U中的左奇异矢量的一个样本。L和Q是两个卷积神经网络，其余块是多维张量或矩阵，其尺寸显示在边缘。

L和Q生成各自矩阵的卷积神经网络中张量，然后在随后的网络操作中将结果重整为堆叠的矩阵表示形式，以便评估矩阵乘积。

传输矩阵T的分量可以表示为从输入视频的奇异值分解（SVD）获得的基本图像的线性组合。这样做计算效率高，又通过限制迭代和位于有效分解子空间中的解来指导优化。

通过输入Z预先计算的截断奇异值分解UΣVT带入到上面的网络中，计算出T，再将计算的TL与Z对比求得损失。

因此问题的核心就变成了用CNN得出Q，使得(UQ)L≈Z。

代码已开源

这么好玩的技术，代码当然开源啦~

GitHub链接如下：
https://github.com/prafull7/compmirrors

这份文档中列出了“矩阵分解”和“光线传输分解”的实现方法。

安装

git clone https://github.com/prafull7/compmirrors
cd compmirrors

配置要求：Python 3.7，以及还需要torch=1.0.1.post2、matplotlib、scipy、visdom等包。

矩阵分解

这步的实现代码在factorization_1d.py文件。可用如下代码运行：

python factorization_1d.py -T ./data/inputs_1d/lightfield.png -L ./data/inputs_1d/tracks_bg.png -o ./outdir_1d 

光线传输分解

一次性训练实现代码在factorization_light_transport.py文件。可用如下代码运行：

export FACTORIZE_DATA_DIR=/path/to/where/data/folders/
export FACTORIZE_OUT_DIR=/path/to/output/directory
python factorization_light_transport.py -d ./data/light_transport/ -f FOLDER_NAME -ds DATASET_NAME -s SEQUENCE_NAME -dev DEVICE_NUMBER

依旧是“鸡生蛋，还是蛋生鸡”问题

虽然这项技术能够重建被隐藏的内容，但是用户还是需要提前知道被隐藏的东西是存在的。

用数学的角度来打个比方。

A和B相乘得80，让你来猜A和B分别是哪两个数字。

可能是40和2，也可能是371.8和0.2152。

在重建工作中，每个像素都会遇到这样的问题——有多种选择。

要让计算机来做选择，那它就会做最简单的事情，得到的结果就是随机的图像。

因此需要规定的算法来做训练。

在此之前，量子位也曾报道过，通过墙壁漫反射的光影，来重建原始画面。

墙上的漫反射如下图所示：

算法还原的图像则是：

这个图像还原实验是，在房间中间随手放置了一个不明位置的遮挡物体，可以是一块不发光的板子，也可以是随手拽过来的一把椅子，阻挡一部分光线到达墙壁。

而这次是完全根据墙上的影子来做图像还原工作。

通过影子的变化可以大致了解房间里隐蔽区域发生了哪些运动。

研究人员同时也对接下来的工作做了展望：

未来，希望能够提高系统的总体分辨率，并最终在不受控制的环境中测试该技术。

传送门

博客：
https://news.mit.edu/2019/using-computers-view-unseen-computational-mirrors-mit-csail-1206

论文：
https://arxiv.org/pdf/1912.02314.pdf

代码：
https://github.com/prafull7/compmirrors

数据：
http://compmirrors.csail.mit.edu/data/dots-sequence.tar