Fast Transient Imaging
A Python framework revolving around reconstruction of impulse responses from AMCW lidar measurements.
GetPisarenkoEstimate.m File Reference

Functions

fliplr(Polynomialroots ()
 

Variables

 function [Phase, Weight]
 
 Polynomial =Eigenvectors(:,iMinEigenvalue)
 
Compute the weights via a least squares system nMoment =size(TrigonometricMoments,1)
 
 Vandermonde =repmat(DiracDeltaLocations.',nMoment,1).^repmat((1:nMoment)',1,nMoment)
 
 Weight =real(Vandermonde\TrigonometricMoments)
 
Complete the output Angle =angle(DiracDeltaLocations)
 
 Phase =Angle+((Angle<0.0).*(2.0*pi))
 

Function Documentation

fliplr(Polynomial' roots ( )
virtual

Variable Documentation

function[Phase, Weight]
Initial value:
=GetPisarenkoEstimate(TrigonometricMoments)
%GetPisarenkoEstimate Given a vector of trigonometric moments without (!) a zeroth
%moment this function computes a distribution consisting of
%size(TrigonometricMoments,1) Dirac-delta distributions which describes this
%measurement correctly.
%\return Phase holds phases in radians (zero to two pi) and Weight holds
% corresponding weights. Both arrays have the same length as the given
% vector. Each entry describes one Dirac-delta.
% Construct the Toeplitz matrix and extract the eigenvector for the smallest
% eigenvalue
Toeplitz=toeplitz([0.0;TrigonometricMoments]')
Toeplitz
Definition: GetMaximumEntropySpectralEstimate.m:9
def GetPisarenkoEstimate(TrigonometricMoments)
Given a vector of trigonometric moments without (!) a zeroth moment this function computes a distribu...
Definition: Reconstruction.py:46
Phase
Definition: GetPisarenkoEstimate.m:23
function[Phase, Weight]
Definition: GetPisarenkoEstimate.m:1
Weight
Definition: GetPisarenkoEstimate.m:20
Polynomial =Eigenvectors(:,iMinEigenvalue)
Compute the weights via a least squares system nMoment =size(TrigonometricMoments,1)
Vandermonde =repmat(DiracDeltaLocations.',nMoment,1).^repmat((1:nMoment)',1,nMoment)
Weight =real(Vandermonde\TrigonometricMoments)
Complete the output Angle =angle(DiracDeltaLocations)
Phase =Angle+((Angle<0.0).*(2.0*pi))