Downloaded 54 times
![wall
hidden object
laser and detector
focus on this point
confocal sampling
histogram
timestamp (nanoseconds)
APJarvis [CC BY-SA 4.0]
lasers and detectors
illuminate and image
same points](https://image.slidesharecdn.com/2019siggraph-190803185738/75/Wave-Based-Non-Line-of-Sight-Imaging-Using-Fast-f-k-Migration-SIGGRAPH-2019-16-2048.jpg)
![NLOS image formation model:
measurements unknown volumetransport matrix
Backprojection [Velten 12, Buttafava 15]
Flops:
Memory:
Runtime: Approx. 10 min.
Iterative Inversion [Gupta 12, Wu 12, Heide 13]
Flops:
Memory:
per iter.
Runtime: > 1 hour
Problem: extremely large in practice
• for n=100, has 1 trillion elements
• for n=1000, sparse needs petabytes of memory
• even matrix-free is computationally intractable](https://image.slidesharecdn.com/2019siggraph-190803185738/75/Wave-Based-Non-Line-of-Sight-Imaging-Using-Fast-f-k-Migration-SIGGRAPH-2019-22-2048.jpg)
![NLOS image formation model:
measurements unknown volumetransport matrix
Backprojection [Velten 12, Buttafava 15]
Flops:
Memory:
Runtime: Approx. 10 min.
Iterative Inversion [Gupta 12, Wu 12, Heide 13]
Flops:
Memory:
per iter.
Runtime: > 1 hour](https://image.slidesharecdn.com/2019siggraph-190803185738/75/Wave-Based-Non-Line-of-Sight-Imaging-Using-Fast-f-k-Migration-SIGGRAPH-2019-23-2048.jpg)
![NLOS image formation model:
measurements unknown volumetransport matrix
Backprojection [Velten 12, Buttafava 15]
Flops:
Memory:
Runtime: Approx. 10 min.
Iterative Inversion [Gupta 12, Wu 12, Heide 13]
Flops:
Memory:
per iter.
Runtime: > 1 hour
3D Deconvolution (with Light-Cone Transform)
[O’Toole et al. 2018]
Flops:
Memory:
Runtime: < 1 second
measurements unknown volumeblur kernel
Confocal scanning and Light-Cone Transform:](https://image.slidesharecdn.com/2019siggraph-190803185738/75/Wave-Based-Non-Line-of-Sight-Imaging-Using-Fast-f-k-Migration-SIGGRAPH-2019-24-2048.jpg)
![NLOS image formation model:
measurements unknown volumetransport matrix
Backprojection [Velten 12, Buttafava 15]
Flops:
Memory:
Runtime: Approx. 10 min.
Iterative Inversion [Gupta 12, Wu 12, Heide 13]
Flops:
Memory:
per iter.
Runtime: > 1 hour
3D Deconvolution (with Light-Cone Transform)
[O’Toole et al. 2018]
Flops:
Memory:
Runtime: < 1 second
measurements unknown volumeblur kernel
Confocal scanning and Light-Cone Transform:](https://image.slidesharecdn.com/2019siggraph-190803185738/75/Wave-Based-Non-Line-of-Sight-Imaging-Using-Fast-f-k-Migration-SIGGRAPH-2019-25-2048.jpg)
![compute memory
f–k migration O(𝒏 𝟑
log 𝒏) O(𝒏 𝟑
)
occlusion NLOS [Heide 19] O(𝒏 𝟓) O(𝒏 𝟓)
LCT [O’Toole 18] O(𝒏 𝟑
log 𝒏) O(𝒏 𝟑
)
backprojection [Velten 12] O(𝒏 𝟓) O(𝒏 𝟑
)
discussion
• cost
• resolution
• runtime
• ~$80k-100k — same as other SPAD-based setups
• 5123 and 10243 hidden volume (2 m)](https://image.slidesharecdn.com/2019siggraph-190803185738/75/Wave-Based-Non-Line-of-Sight-Imaging-Using-Fast-f-k-Migration-SIGGRAPH-2019-51-2048.jpg)
![Fermat Paths [Xin et al. 2019]
Xin et al. 2019
Reconstructs surface geometry from
discontinuities in histogram measurements
- surface reconstruction
- works with specular objects
- single object per scene
- noise sensitivity (requires identifying
discontinuities)](https://image.slidesharecdn.com/2019siggraph-190803185738/75/Wave-Based-Non-Line-of-Sight-Imaging-Using-Fast-f-k-Migration-SIGGRAPH-2019-54-2048.jpg)
![wall
hidden object
laser and detector
focus on this point
confocal sampling
histogram
timestamp (nanoseconds)
APJarvis [CC BY-SA 4.0]
lasers and detectors
illuminate and image
same points](https://crownmelresort.com/image.slidesharecdn.com/2019siggraph-190803185738/75/Wave-Based-Non-Line-of-Sight-Imaging-Using-Fast-f-k-Migration-SIGGRAPH-2019-16-2048.jpg)
![NLOS image formation model:
measurements unknown volumetransport matrix
Backprojection [Velten 12, Buttafava 15]
Flops:
Memory:
Runtime: Approx. 10 min.
Iterative Inversion [Gupta 12, Wu 12, Heide 13]
Flops:
Memory:
per iter.
Runtime: > 1 hour
Problem: extremely large in practice
• for n=100, has 1 trillion elements
• for n=1000, sparse needs petabytes of memory
• even matrix-free is computationally intractable](https://crownmelresort.com/image.slidesharecdn.com/2019siggraph-190803185738/75/Wave-Based-Non-Line-of-Sight-Imaging-Using-Fast-f-k-Migration-SIGGRAPH-2019-22-2048.jpg)
![NLOS image formation model:
measurements unknown volumetransport matrix
Backprojection [Velten 12, Buttafava 15]
Flops:
Memory:
Runtime: Approx. 10 min.
Iterative Inversion [Gupta 12, Wu 12, Heide 13]
Flops:
Memory:
per iter.
Runtime: > 1 hour](https://crownmelresort.com/image.slidesharecdn.com/2019siggraph-190803185738/75/Wave-Based-Non-Line-of-Sight-Imaging-Using-Fast-f-k-Migration-SIGGRAPH-2019-23-2048.jpg)
![NLOS image formation model:
measurements unknown volumetransport matrix
Backprojection [Velten 12, Buttafava 15]
Flops:
Memory:
Runtime: Approx. 10 min.
Iterative Inversion [Gupta 12, Wu 12, Heide 13]
Flops:
Memory:
per iter.
Runtime: > 1 hour
3D Deconvolution (with Light-Cone Transform)
[O’Toole et al. 2018]
Flops:
Memory:
Runtime: < 1 second
measurements unknown volumeblur kernel
Confocal scanning and Light-Cone Transform:](https://crownmelresort.com/image.slidesharecdn.com/2019siggraph-190803185738/75/Wave-Based-Non-Line-of-Sight-Imaging-Using-Fast-f-k-Migration-SIGGRAPH-2019-24-2048.jpg)
![NLOS image formation model:
measurements unknown volumetransport matrix
Backprojection [Velten 12, Buttafava 15]
Flops:
Memory:
Runtime: Approx. 10 min.
Iterative Inversion [Gupta 12, Wu 12, Heide 13]
Flops:
Memory:
per iter.
Runtime: > 1 hour
3D Deconvolution (with Light-Cone Transform)
[O’Toole et al. 2018]
Flops:
Memory:
Runtime: < 1 second
measurements unknown volumeblur kernel
Confocal scanning and Light-Cone Transform:](https://crownmelresort.com/image.slidesharecdn.com/2019siggraph-190803185738/75/Wave-Based-Non-Line-of-Sight-Imaging-Using-Fast-f-k-Migration-SIGGRAPH-2019-25-2048.jpg)
![compute memory
f–k migration O(𝒏 𝟑
log 𝒏) O(𝒏 𝟑
)
occlusion NLOS [Heide 19] O(𝒏 𝟓) O(𝒏 𝟓)
LCT [O’Toole 18] O(𝒏 𝟑
log 𝒏) O(𝒏 𝟑
)
backprojection [Velten 12] O(𝒏 𝟓) O(𝒏 𝟑
)
discussion
• cost
• resolution
• runtime
• ~$80k-100k — same as other SPAD-based setups
• 5123 and 10243 hidden volume (2 m)](https://crownmelresort.com/image.slidesharecdn.com/2019siggraph-190803185738/75/Wave-Based-Non-Line-of-Sight-Imaging-Using-Fast-f-k-Migration-SIGGRAPH-2019-51-2048.jpg)
![Fermat Paths [Xin et al. 2019]
Xin et al. 2019
Reconstructs surface geometry from
discontinuities in histogram measurements
- surface reconstruction
- works with specular objects
- single object per scene
- noise sensitivity (requires identifying
discontinuities)](https://crownmelresort.com/image.slidesharecdn.com/2019siggraph-190803185738/75/Wave-Based-Non-Line-of-Sight-Imaging-Using-Fast-f-k-Migration-SIGGRAPH-2019-54-2048.jpg)
More Related Content
![[DL輪読会]医用画像解析におけるセグメンテーション](https://cdn.slidesharecdn.com/ss_thumbnails/20190301fujino4-190322072121-thumbnail.jpg?width=640&height=640&fit=bounds)
![[DL輪読会]Flow-based Deep Generative Models](https://cdn.slidesharecdn.com/ss_thumbnails/20190307-190328024744-thumbnail.jpg?width=640&height=640&fit=bounds)
![[DL輪読会]NeRF: Representing Scenes as Neural Radiance Fields for View Synthesis](https://cdn.slidesharecdn.com/ss_thumbnails/nerfdlseminar1-200327021512-thumbnail.jpg?width=640&height=640&fit=bounds)
![[DL輪読会]Neural Ordinary Differential Equations](https://cdn.slidesharecdn.com/ss_thumbnails/nnasode1-190111001755-thumbnail.jpg?width=640&height=640&fit=bounds)
![SSII2022 [TS3] コンテンツ制作を支援する機械学習技術〜 イラストレーションやデザインの基礎から最新鋭の技術まで 〜](https://cdn.slidesharecdn.com/ss_thumbnails/ts32022ssiiess-220607054523-e80be8dc-thumbnail.jpg?width=640&height=640&fit=bounds)
Similar to Wave-Based Non-Line-of-Sight Imaging Using Fast f–k Migration | SIGGRAPH 2019
Wave-Based Non-Line-of-Sight Imaging Using Fast f–k Migration | SIGGRAPH 2019
- 1. Wave-Based Non-Line-of-Sight Imaging UsingFast f–k Migration David B. Lindell, Gordon Wetzstein, Matthew O’Toole Stanford University Carnegie Mellon University 08/2019
- 2.
- 4.
- 5.
- 6. resolution: 128 x128 area: 2 m × 2 m
- 7.
- 8. x y z Dimensions: 2 mx 2 m x 1.5 m scene photo reconstruction
- 9. contributions • fast wave-basedimage formation model for NLOS imaging • complex surface reflectances, more robust to noise • new hardware prototype: room-sized scenes, interactive scanning
- 10. NLOS imaging forautonomous cars
- 11. LIDAR (light detectionand ranging) Velodyne VLS-128 NLOS imaging for autonomous cars
- 12. NLOS image formationand related work hardware prototype outline wave-based model
- 13. picosecond laser SPAD sensor * timestamp ofphoton event Time-to-Digital Converter (TDC) scene Num.Detections Time of Flight single-photon avalanche diodes (SPADs)
- 14.
- 15. wall hidden object laser detector histogram timestamp(nanoseconds) non-confocal sampling
- 16. wall hidden object laser anddetector focus on this point confocal sampling histogram timestamp (nanoseconds) APJarvis [CC BY-SA 4.0] lasers and detectors illuminate and image same points
- 17. wall hidden object confocal sampling samepath to the object and back
- 18. wall hidden object confocal sampling -simplified NLOS mathematical model - enables efficient NLOS reconstruction equivalent to one-way propagation at half-speed
- 19.
- 20.
- 21. sensor light source occluder imaging area 3Dmeasurements
- 22. NLOS image formationmodel: measurements unknown volumetransport matrix Backprojection [Velten 12, Buttafava 15] Flops: Memory: Runtime: Approx. 10 min. Iterative Inversion [Gupta 12, Wu 12, Heide 13] Flops: Memory: per iter. Runtime: > 1 hour Problem: extremely large in practice • for n=100, has 1 trillion elements • for n=1000, sparse needs petabytes of memory • even matrix-free is computationally intractable
- 23. NLOS image formationmodel: measurements unknown volumetransport matrix Backprojection [Velten 12, Buttafava 15] Flops: Memory: Runtime: Approx. 10 min. Iterative Inversion [Gupta 12, Wu 12, Heide 13] Flops: Memory: per iter. Runtime: > 1 hour
- 24. NLOS image formationmodel: measurements unknown volumetransport matrix Backprojection [Velten 12, Buttafava 15] Flops: Memory: Runtime: Approx. 10 min. Iterative Inversion [Gupta 12, Wu 12, Heide 13] Flops: Memory: per iter. Runtime: > 1 hour 3D Deconvolution (with Light-Cone Transform) [O’Toole et al. 2018] Flops: Memory: Runtime: < 1 second measurements unknown volumeblur kernel Confocal scanning and Light-Cone Transform:
- 25. NLOS image formationmodel: measurements unknown volumetransport matrix Backprojection [Velten 12, Buttafava 15] Flops: Memory: Runtime: Approx. 10 min. Iterative Inversion [Gupta 12, Wu 12, Heide 13] Flops: Memory: per iter. Runtime: > 1 hour 3D Deconvolution (with Light-Cone Transform) [O’Toole et al. 2018] Flops: Memory: Runtime: < 1 second measurements unknown volumeblur kernel Confocal scanning and Light-Cone Transform:
- 27. NLOS image formationand related work hardware prototype outline wave-based model
- 28. z x confocal measurements x wavefield wall (z= 0) hidden object t image formation model
- 29. z x x t general solution(time reversal) wall (z = 0) hidden object wavefield confocal measurements
- 30. general solution (timereversal) 1. approximate wave equation with finite differences 2. solve for previous timestep 3. repeatedly update at all grid cells finite-difference time-domain method
- 31. general solution (timereversal) 1. approximate wave equation with finite differences 2. solve for previous timestep 3. repeatedly update at all grid cells finite-difference time-domain method
- 32. z x x t frequency–wavenumber (f–k)Migration wall (z = 0) hidden object wavefield confocal measurements FLOPS: O(𝒏 𝟑 log 𝒏)
- 33.
- 34.
- 35. f-k Migration Measurements (z=0)Spectrum Hidden Volume (t=0)Interpolated Spectrum Resample
- 36. f-k Migration Express wavefieldas function of measurement spectrum (plane wave decomposition) wavefield Fourier transform of measurements Set t=0 to get migrated solution Almost an inverse Fourier Transform!
- 37. f-k Migration Set t=0to get migrated solution Almost an inverse Fourier Transform! Use dispersion relation1 to perform substitution of variables 1Georgi, Howard. The physics of waves. Englewood Cliffs, NJ: Prentice Hall, 1993.
- 38. Use dispersion relation1to perform substitution of variables
- 39. Use dispersion relation1to perform substitution of variables
- 40. Use dispersion relation1to perform substitution of variables The migrated solution is an inverse Fourier Transform! Resample
- 41. f-k Migration x y z Dimensions: 2x 2 m Exposure: 180 min Reconstruction time: ~1 min (CPU)
- 42.
- 43. NLOS image formationand related work hardware prototype outline wave-based model
- 44.
- 45.
- 46.
- 47. real-time scanning Framerate: 4Hz Resolution: 32 x 32 Dimensions: 2 m x 2 m x 2 m Reconstruction time: ~1 s per frame
- 48. noise sensitivity LCT f-k Migration 15 sec1 min 2 min 15 min x y z
- 49. noise comparison Filtered BackprojectionLCT f-k Migration 10 min. exposure
- 50. non-confocal measurements phaseretrieval without phase retrieval with phase retrieval non-planar scan surfaces non-planar ‘wall’
- 51. compute memory f–k migrationO(𝒏 𝟑 log 𝒏) O(𝒏 𝟑 ) occlusion NLOS [Heide 19] O(𝒏 𝟓) O(𝒏 𝟓) LCT [O’Toole 18] O(𝒏 𝟑 log 𝒏) O(𝒏 𝟑 ) backprojection [Velten 12] O(𝒏 𝟓) O(𝒏 𝟑 ) discussion • cost • resolution • runtime • ~$80k-100k — same as other SPAD-based setups • 5123 and 10243 hidden volume (2 m)
- 52. conclusion • confocal sampling+ wave-based models enables fast, robust NLOS imaging • fast scanning and outdoor reconstruction is possible with optimized hardware setup • moving towards practical NLOS imaging “in the wild”
- 53. lindell@stanford.edu davidlindell.com code & data: computationalimaging.org github.com/computational-imaging/nlos-fk Wave-BasedNon-Line-of-Sight Imaging Using Fast f–k Migration Gordon Wetzstein Stanford Matthew O’Toole CMU
- 54. Fermat Paths [Xinet al. 2019] Xin et al. 2019 Reconstructs surface geometry from discontinuities in histogram measurements - surface reconstruction - works with specular objects - single object per scene - noise sensitivity (requires identifying discontinuities)
Editor's Notes
- #2 I’m David Lindell, and I’ll be talking about a wave-based model for non-line-of-sight imaging
- #3 In non-line-of-sight imaging, we image objects that are hidden from direct line of sight. To do this we measure the time it takes for short pulses of light to bounce off a wall and scatter back from a hidden object.
- #4 These time delays give us the distance from the wall to the hidden object. But we can also use these types of measurements to recover the full 3D geometry.
- #5 Here’s an example scene that we captured outdoors using the side of a building. This is a room sized scene that we’re capturing under indirect sunlight
- #6 From the perspective of the imaging system, we can only see the wall that we’re illuminating and imaging.
- #7 From behind the scene, though, we can observe that the light scatters off of the wall and reflects back from the objects in the scene. This light eventually scatters back to our detector, where we capture the timing information.
- #8 The measurements capture an ultra high-speed view of light returning from the hidden objects and scattering against the wall.
- #9 From these measurements, we reconstruct the 3D geometry. The reconstruction we get is volumetric, with the intensity of each voxel corresponding to the brightness or albedo of the object.
- #10 The key contributions of our work are to introduce a fast wave-based image formation model for NLOS imaging which overcomes some of the limitations of previous approaches based on geometric optics. In particular, our approach can handle complex surface reflectances and is more robust to noise than previous approaches. We also built a new hardware prototype and demonstrate imaging room-sized scenes and NLOS imaging at interactive rates for the first time.
- #11 With these new capabilities we move towards enabling NLOS imaging in more practical applications like autonomous vehicle navigation.
- #12 The current LIDAR systems for these vehicles are using basically the same detectors and lasers that we use in our imaging system. However, they are limited to imaging LOS points. Being able to add the 3D information of objects around the corner and outside direct line of sight could help cars navigate more safely.
- #13 for the remainder of the talk, I’ll go into more detail about NLOS image formation and related work, explain our wave-based model, and then describe our hardware prototype and some additional results.
- #14 To measure the travel time of photons, we use a detector called a single-photon avalanche diode or SPAD We pair the SPAD with a picosecond laser that sends millions of short pulses of light into the scene every second. Each pulse interacts with the scene and sends some scattered light back to our detector. For each pulse emitted into the scene, the SPAD has a chance at detecting the arrival of a single photon. When a detection occurs, we get a timestamp of when the photon arrived with respect to when the laser pulse was emitted. These timestamps are accumulated and binned into a histogram showing how much light is arriving over picosecond time intervals.
- #15 In NLOS imaging, we align the detector and the laser so that we illuminate and image the same spot on the wall. The measurements capture the number of photon counts over time. We observe a large peak from the direct path to the wall (which I’ve positioned at time zero), and then some nanoseconds later a secondary peak from the hidden object.
- #16 We can also align the laser and detector to illuminate and image different points, which we call non-confocal sampling. In this case, the measured time delay gives us the round-trip distance from the laser point to the hidden object to the detector point.
- #17 Alternatively, confocal sampling is illuminating and imaging the same point. With this sampling geometry, the time delay in the measurements directly captures the distance between the wall and the hidden object. Confocal sampling is also nice because LOS LIDARs are designed to illuminate and image the same point.
- #18 Confocal sampling is equivalent to having each point on the hidden object emit light, and the resulting wavefront propagates at half the speed of light to the wall. This equivalent, but simplified model enables efficient reconstructions and existing LIDAR systems are designed for confocal sampling.
- #19 Confocal sampling is equivalent to having each point on the hidden object emit light, and the resulting wavefront propagates at half the speed of light to the wall. This equivalent, but simplified model enables efficient reconstructions and existing LIDAR systems are designed for confocal sampling.
- #20 Here’s an example measurement that we captured where we illuminate and image a spot on the wall. As I move back and forth with this exit sign, we can observe the time difference to the secondary peak.
- #21 While illuminating and imaging one point gives the distance to the hidden object
- #22 We can also scan a patch on the wall to capture a 3D measurement volume.
- #23 we use these measurements to reconstruct 3D geometry.
- #24 Using geometric optics, the measurements and hidden volume can be related through a transport matrix, A. Backprojection and direct inversion techniques have been proposed for reconstruction, but have computational requirements that scale as N^5. For example, for a low-resolution volume where n=100, A has 1 trillion elements. For n=1000 we would need petabytes of memory to store the matrix. But even matrix free implementations have intractable computational requirements. (replace “A” with latex font)
- #25 Thus these techniques are largely intractable for high-resolution reconstructions
- #26 To address this, another approach uses confocal scanning and the light-cone transform or LCT. This technique transforms the measurements to a domain where the reconstruction can be performed with deconvolution.
- #27 While this method is fast, it is limited to diffuse and retroreflective hidden objects.
- #28 Instead of using a geometric image formation model, we rely on a wave model based on the time-dependent wave equation. This model extends to varied surface reflectances and we can solve the reconstruction without sacrificing any computational speed. Wave-based models have been used in the computer graphics community for simulation, rendering, and sound. But wave based models have also been used for radar, sonar, and seismic imaging. This is incredibly useful because by recognizing the connections to these domains we can leverage decades of research for this new application of NLOS imaging.
- #29 One example of where this wave-based model is used…
- #30 is in seismic imaging, which turns out to have many similarities to NLOS imaging. Here, shockwaves caused by explosions on the surface cause waves to propagate downward, reflect off of hidden underground surfaces, and are picked up by receivers at the earth’s surface. The measurements are used to reconstruct the underground structures.
- #31 In optical non line of sight, the laser illumination is like the shockwave which is reflected back by the hidden object.
- #32 We can model the confocal measurements as the value of a wavefield over time that starts from the surface of the hidden object at time zero and propagates to the observed wall. I’ll show the whole wavefield over x,z, and time on the left, and then the captured measurements evaluated at the wall on the right.
- #33 Our goal is then to migrate the wavefield from z=0 back to t=0 to recover the surface of the hidden object. One way to do this is to time-reverse the measurements and solve the differential equation using finite difference techniques, but this is extremely expensive to compute. For example, this simple 2d simulation took 20 minutes to compute.
- #34 the finite difference technique consists of discretizing the measurement domain, and then approximating the differential wave equation with finite differences. If we do this, we can solve for the next time-reversed step in terms of our current and previous measurement steps. This corresponds to an update step that we repeatedly compute and apply until we reach t=0.
- #35 In general, though, this iterative method is slower than necessary if we only need to recover the value at time zero.
- #36 Another recent approach to solving this migration problem was called phasor fields. In this approach, the measurements are convolved with a virtual wave function– a wavelet for example– to simulate the scene response to a known wave. Then the resulting wavefield is propagated backwards using Rayleigh sommerfield diffraction or Fresnel propagation. The benefit here is that you can also propagate to arbitrary time values; the computational requirements depend on the propagation method and can be as fast as n^3 log n
- #37 Our approach adapts f-k migration which is used in radar and seismology to efficiently go from one boundary condition to another. In this case, we can directly migrate the measurements to the hidden surface. This approach also has the benefit of being fast and uncomplicated.
- #38 Now I’m going to show this for captured measurements in 3D. On the left is an image of the room-sized hidden scene with different surface reflectances. On the right are the captured measurements. The concentrated streaks at the beginning of the volume come from the glossy dragon, the bright spots are specular reflections from the discoball, and the broad elongated streaks at the most distant part correspond to diffuse reflections from the statue.
- #39 I’ll play back the measurements as an ultrafast video, and we can take a minute to watch the ripples of light that we observe at the wall from the hidden scene. (fade in music)
- #40 Here’s how the reconstruction works. We take the fourier transform of our measurements and resample. Then, an inverse fourier transform reconstructs the hidden volume. Let me focus on why this resampling step works.
- #41 We express the wavefield as a function of the measurement spectrum or as an integral over plane waves at different frequencies. To get the expression for the migrated wavefield, we set t=0. The resulting integral expression gives us the surface of the hidden volume but is expensive to compute numerically. However, we see that it is almost an inverse fourier transform except for the spatial and temporal frequency variables do not match
- #42 We can use the dispersion relation to massage this equation into the right form. This relation can be derived from the wave equation or found in physics textbooks and relates the spatial and temporal frequencies of electromagnetic waves.
- #43 The temporal frequency of the wave is given by one over the period
- #44 its spatial wavenumber is given by 2 pi over the wavelength
- #45 the dispersion relation describes how fast waves at one temporal frequency propagate in all spatial directions.
- #46 and so this equation describes exactly how to resample the hidden volume
- #47 this allows us to convert from temporal frequency f to spatial frequency k_z by rearranging the fourier coefficients
- #48 After the substitution of variables, we find that the migrated wavefield is given by an inverse fourier transform, and this gives us our reconstructed volume.
- #49 Here’s the reconstruction using f-k migration, we can even reconstruct the position of the disco ball, which sits on the edge and sticks out from the shelf.
- #50 These methods are really designed for diffuse reflectances so they contain artifacts in the reconstruction, f-k migration produces a cleaner result.
- #51 need to highlight key contributions, key differences, add outline (1. image formation model, 2. hardware, etc.)
- #52 We captured additional results with a new hardware prototype that is shown here.
- #53 The main part of the prototype consists of a high-power laser and spad which share an optical axis. Two scanning mirrors allow us to scan a wall or other scan surface to probe the hidden scene.
- #54 Here’s a video of the hardware prototype in action, the laser sends out 35 ps laser pulses 10 million times every second.
- #55 We also wanted to see how fast we could push the scanning of this system, and so I dressed up in a retroreflective suit to increase the amount of available signal. By scanning at 4 Hz, we can capture human pose and location at low resolution from around the corner.
- #56 Our empirical results show that f-k migration is robust at low exposure settings. The LCT result is more noisy, partly due to it not modeling the glossy surface. Also, f-k migration requires no parameter tuning.
- #57 The improved quality of f-k migration is also apparent here, especially compared to filtered backprojection for a 10 min. exposure time.
- #58 Since a confocal scanning arrangement can be somewhat restrictive and because not all scan surfaces are planar walls, we also derive in the paper computational corrections to apply f-k migration in these scenarios. Additionally, we show how phase retrieval techniques can be used to improve reconstruction results. (add concluding slide)
- #59 in terms of cost, our hardware prototype is expensive but about the same as any other SPAD-based setup. However, these costs will probably come down as more commercial lidar systems emerge. We reconstruct hidden volumes with a size of 512^3 voxels we also have data at 1024^3 resolution. However, at these resolutions for a 2x2x2 m volume we’re sampling at or above fundamental resolution limits for diffuse objects that we derived in our 2018 paper published in Nature. Higher resolutions can only be achieved with specular or retroreflective objects or priors that hallucinate details. handling occlusions could improve reconstruction quality, it comes at a cost for runtime and memory. Our method could also be implemented on a GPU for real-time reconstruction.
- #60 in terms of cost, our hardware prototype is expensive but about the same as any other SPAD-based setup. However, these costs will probably come down as more commercial lidar systems emerge. We reconstruct hidden volumes at 512^3 resolution and we also have data at 1024^3 resolution. However, at these resolutions we’re sampling at or above fundamental resolution limits for diffuse objects that we derived in our 2018 nature paper. Higher resolutions can only be achieved with specular or retroreflective objects or priors that hallucinate details. While handling occlusions could improve reconstruction quality, it comes at a cost for runtime and memory. Our method could also be implemented on a GPU for real-time reconstruction.
- #61 in conclusion, we’ve demonstrated NLOS imaging based on confocal scanning and wave-based models that is fast and more robust than other methods We’ve also demonstrated outdoor reconstruction and interactive scanning for the first time and take steps towards making NLOS imaging practical in everyday scenarios for applications like autonomous vehicle navigation
- #62 We’ve also made our code and data available on our project webpage.
- #63 At CVPR the best paper award was given to a paper that related the principle of fermat paths to nlos imaging. This works by identifying discontinuities in the histograms of photon counts and using these to recover a hidden object surface. This technique is really unique in being able to work on surfaces and surface normals rather than a volumetric representation. It also works with specular objects. However, in practice this technique has only been shown with a single object per scene, and it may be more difficult to have this work with many separate objects. Since the technique requires detecting discontinuities in the measurements, it may also be more sensitive to noise than our approach. It may also be possible to use both methods to further improve reconstruction quality.