Questions tagged [deconvolution]

in mathematics, the inverse operation of convolution signals. In general, the purpose of deconvolution is to find solutions of the convolution equation defined as: f * g = x. Where h is the recorded signal, and f is a signal that you want to recover, and we know that the first signal is obtained by convolution of the second with some known signal g. If the signal g is unknown, it has to be estimated (eg. statistical estimation).

Filter by
Sorted by
Tagged with
0 votes
2 answers
43 views

How to use deconvolution code with an impulse response to achieve the original signal

I am trying to write a code for my thesis to deconvolve a recording with an impulse response so that I can achieve the original audio signal. I have written a simple code so that I can implement this ...
user avatar
5 votes
2 answers
95 views

Can the deconvolution Wiener filter reduce noise without having a blurred image?

I am trying to denoise many several noises with several filters for a research i have, i found a deconvolution Wiener filter made by "mr.tranleanh" on Github, as you can see here . what I ...
user avatar
  • 53
1 vote
1 answer
77 views

Causal inverse of $h[n]=\delta[n]-\alpha\delta[n-1]$

Find the causal inverse of $$h[n]=\delta[n]-\alpha\delta[n-1]$$ we have $h[0]=1$ and $h[1]=-\alpha$ also $h[n]=0$ for $n>1$ From the formula $$ h_i[n]=\sum_{i=1}^n\frac{h[n]h_i[n-i]}{h[0]} $$ we ...
user avatar
  • 147
1 vote
0 answers
30 views

Deconvolution of signal with harmonic distortion

I came across these measurements done using a stepped sine sweep covering only some frequencies (~800 frequency points). I would like to get an impulse response from these measurements to convolve ...
user avatar
  • 11
4 votes
2 answers
148 views

Biexponential (double exponential) convolution of a function

Summary I am trying to run a convolution on some data that was originally calculated from a deconvolution (so the reverse). However I'm not getting the expected graph. Blue is expected, red is a ...
user avatar
  • 161
1 vote
0 answers
53 views

Deconvolution: how to build Discrete-Time Impulse Response matrix?

I am reading a paper about the Hunt problem: Consider the iput given by $$e(t) = \mathrm{e}^{-\left(\frac{t-400}{75}\right)^2} - \mathrm{e}^{-\left(\frac{t-600}{75}\right)^2},\quad 0\leq t \leq 1025$$ ...
user avatar
  • 111
1 vote
0 answers
138 views

Unexpected results of deconvolution with scipy.deconvolve

Below I have plotted the signal (Lifetime decay) I am trying to deconvolve from a known impulse response function (IRF), as well as the IRF itself. I'm using scipy.signal.deconvolve. Please note for ...
user avatar
0 votes
0 answers
109 views

How to handle zeros before FFT convolution / deconvolution?

I would like to calculate the input function (unknown) by deconvolution of the output and the "system response" signals. The output is a finite signal from a measure device so it presents ...
user avatar
0 votes
0 answers
23 views

Unexpected result of motion-blurred image deconvolution

In my case I have some images, captured by CMOS-camera (global shutter) during non accelerated motion (with fixed illumination and focus and known velocity and exposure time, so field of view travels ...
user avatar
  • 1
2 votes
3 answers
71 views

Approximating inverse of unstable difference of Gaussians filter

I am trying to invert a difference of Gaussians (DoG) filter. The inverse is not stable and so I am trying to find an approximation applied to a specific input. The DoG filter increases contrast at ...
user avatar
  • 21
0 votes
0 answers
81 views

How to smooth wiener deconvolution result in Python?

I'm wondering if it is possible to smooth the estimated response from a Wiener deconvolution in order to have a better representation of the original signal and to remove the side lobes. (Here an ...
user avatar
1 vote
2 answers
263 views

How to use deconvolution technique to find out impulse response?

I have been working to find out room for impulse response. I am using Logarithmic sweep sine wave as input say $x(n)$ and my recorded signal is $y(n)$. I know the room impulse response is ...
user avatar
0 votes
1 answer
45 views

What is $S$ in the Wiener filter exactly?

I am reading about the Wiener deconvolution in Wikipedia. In the expression of $G$ we have $$G=\frac{H^*S}{|H|^2S+N}$$ where $S = \mathbb{E}|X|^2$. Why do we have the $\mathbb{E}$ symbol? Isn't $X$ a ...
user avatar
5 votes
1 answer
62 views

Estimate the Image Using Multi Realizations of Its Convolution with a Known Filters Using Wiener Filter

Suppose we have a corrupted image $Y = H*X + \epsilon$ formed by taking an image $X$, convolving it with a point-spread function $H$, and adding gaussian noise $\epsilon$. Then we know that the Wiener ...
user avatar
2 votes
1 answer
802 views

Deconvolution with Python in real life

I have measured a signal which is convolved with the profile of the measuring apparatus. Now I want to remove this contribution to get the "real" signal. I am trying to do this with Python. ...
user avatar
3 votes
1 answer
71 views

Solve Efficiently the 1D $ L_1 $ Regularized Least Squares Problem (Denoising / Deblurring)

How to solve a 1D Least Squares with $ L_1 $ Regularization? I know gradient based method, I wonder how much faster / efficient I can get. Related to Solve Efficiently the 1D Total Variation ...
user avatar
  • 640
6 votes
1 answer
122 views

Solve Efficiently the 1D Total Variation Regularized Least Squares Problem (Denoising / Deblurring)

How to solve a 1D Least Squares with Total Variation Regularization? I know gradient based methods, I wonder how much faster / efficient I can get.
user avatar
  • 640
5 votes
1 answer
82 views

Converting Hadamard Product into Matrix Multiplication in Image Deconvolution with Total Variation (TV) Using ADMM

I would like to solve the following Image Deconvolution equation by ADMM. $$\mathbf { \min\frac{1}{2}\Vert{Cx-b}\Vert_2^2+\Vert w\circ (D x)\Vert_1 \tag 1}$$ Where, $x$ is a vector of unknown pixel ...
user avatar
4 votes
2 answers
943 views

Generate the Matrix Form of 1D Convolution Kernel

As a follow up to Generate the Matrix Form of 2D Convolution Kernel, could someone explain how to generate the matrix form of a 1D convolution kernel? How different convolutions shapes are handled? ...
user avatar
  • 640
6 votes
1 answer
465 views

Deblurring 1D Data Using Direct Inverse Filtering

In my assignment I have been given recorded temperature over a period of time (193 Values) and the Impulse Response (5 Values with n=0 corresponding to the middle value). data : data.csv | ...
user avatar
1 vote
0 answers
86 views

How to update point spread function of blind deconbolution by conjugate gradient?

There is an unblurred image $g$ and a blurred image $x$. Their relationship is expressed by the following formula using $psf$(point spread fucntion, size is $5×5$ kernel). $g = x \otimes psf\tag 1$ ...
user avatar
2 votes
2 answers
329 views

Deconvolution of system response in Python/Matlab

I had two sets of data, the output function of the system (time series with a length of 1292 entries) and the transfer function (similar to a gaussian with a length of 681 entries). I would like to ...
user avatar
4 votes
1 answer
113 views

How to Solve Blind Image Deblurring with Total Variation (TV) Prior Using ADMM?

As a continuation of the question How to Solve Non Blind Image Deblurring with Total Variation Prior Using ADMM? I would like to understand how could one solve the Blind Deblurring (Deconvolution) ...
user avatar
  • 415
6 votes
1 answer
149 views

How to Solve Non Blind Image Deblurring with Total Variation Prior Using ADMM?

How could one use the Total Variation frame work to solve the Deblurring problem? Specifically using the ADMM as a solver. One could assume the blurring operator is known, linear and shift invariant. ...
user avatar
  • 415
0 votes
2 answers
55 views

Deconvolution of shifted gaussian function in the frequency range

I have a signal defined as $$A(t)\cdot\exp\left(-i\omega_0t\right)$$ with $A$ the envelope function and $\omega_0$ the carrier frequency. I would like to transfer this signal into the fourier space ...
user avatar
  • 125
2 votes
1 answer
56 views

Differences Between Two $ {L}_{1} $ Norm Minimization Schemes

I was reading and working with L1 regularized least squares, where: $$ \arg \min_{\boldsymbol{x}} \frac{1}{2} {\left\| A \boldsymbol{x} - \boldsymbol{y} \right\|}_{2}^{2} + \lambda {\left\| \...
user avatar
  • 93
0 votes
1 answer
70 views

Signal deconvolution from exponential and linear tail

I would like to deconvolve the signal shown in the plot below from its linearly decaying tail and the exponential levelling. That means that after the deconvolutin, the tail should be at a constant ...
user avatar
  • 1
4 votes
1 answer
151 views

Frequency-domain deconvolution: "Direct" filtering vs "Wiener" filtering

Can someone help with clarifying the difference between two approaches to frequency domain "deconvolution: For the frequency domain problem: We want to find a filter $F(\omega)$ which will ...
user avatar
2 votes
0 answers
77 views

calculate or decompose a Fourier transform signal amplitudes with unknown weights on sources

migrated from math-se... I am trying to calculate , or approximate the solution of following Fourier-sine transform problem that can be expressed as a contributions of periodic sources $f_i(x)$ and ...
user avatar
  • 93
8 votes
4 answers
982 views

Deconvolution to Remove Gaussian Blur in 1D Signal (Wiener Filtering?)

I've got a set of biology data that I'm trying to denoise (effectively, a population statistic can only be measured convolved with a gaussian of known width) My problem is this: I can measure (f*g), ...
user avatar
0 votes
0 answers
28 views

Deconvolution of sidelobes in a point spread function?

It seems that most deconvolution algorithms mainly handle the main lobe of a point spread function (PSF) and assume that sidelobes can be safely neglected. For a direct algorithm trying to perform a ...
user avatar
  • 101
1 vote
0 answers
128 views

2D Deconvolution using a non-gaussian mask using C++

I am currently working on a project, where we record an electron beam profile using a target. The obtained image is a result of convolution of the actual beam profile and the aperture wherein the ...
user avatar
5 votes
2 answers
183 views

Wiener filtering/deconvolution for non-stationary noise

Consider a stationary discrete-time random process $x[k]$ which undergoes low pass filtering by a filter with impulse response $h[k]$ and is subject to additive, temporally uncorrelated noise $n[k]$ ...
user avatar
  • 355
2 votes
1 answer
249 views

What is the purpose of wrapping the negative times of a response function in discrete convolution?

I am trying to rationalize a figure given in the Numerical Recipes in C in the section of Fourier based convolution and deconvolution. The authors show the example of a continuous convolution with a ...
user avatar
  • 498
6 votes
1 answer
63 views

Stochastic Methods for Image Deconvolution Problem

If we convolve an image with a point spread function and from the resulting image to find the input image, can we use any stochastic approaches? I feel like we will not be able to. A single image ...
user avatar
  • 63
1 vote
0 answers
67 views

How to estimate the system characteristic function given experimental input and output

I have experimental signals $y_i(t)$ for $i = 1,\ldots,n$ that correspond to different excitation inputs to a system $x_i(t)$ for $i = 1,\ldots,n$. The goal is to find the system characteristic ...
user avatar
  • 141
7 votes
3 answers
387 views

Deconvolution of a 1D Time Domain Wave Signal Convolved with Series of Rect Signals

I have a synthesized signal (the bottom of the following figure), which is the convolution of the input signal (at the top) and the objective function (in the middle). The intention is to retrieve the ...
user avatar
  • 73
0 votes
0 answers
60 views

Deconvolution with noisy measurement of impulse response function

I observe a noisy 1D signal that has passed through some linear time-invariant filter: $$z(t) = (h\star x)(t) + n(t)$$ where $h(t)$ is the filter, and I am hoping to recover $x(t)$. $n(t)$ is some ...
user avatar
  • 269
2 votes
0 answers
60 views

Adapting Richardson Lucy (RL) Deconvolution for Shot Noise Limited Coherent Imaging

I am an experimental physicist who is collecting a series of coherent imaging of trapped gas. If you are familiar with phase contrast imaging, you may understand what I mean by coherent imaging. The ...
user avatar
  • 121
0 votes
0 answers
81 views

Richardson-Lucy iterative deconvolution is producing erroneous deconvoluted images

I am working on producing deconvoluted image from degraded image where the process of degradation is known to be linear and space invariant. I have used Richardson-Lucy iterative deconvolution method ...
user avatar
2 votes
2 answers
604 views

GnuRadio CC Decoder or How to use FEC Extended Decoder

I am currently trying to recover a satellite signal, encoded following CCSDS standards. This includes a convolutional code with rate 1/2 and constraint length 7. I am receiving the signal through an ...
user avatar
0 votes
1 answer
38 views

Recovering time domain signal from sum of independent and identically distributed signals [closed]

Background: I observe a sample of a variable z that is the sum of two independent and identically distributed variables x and <...
user avatar
  • 101
5 votes
1 answer
95 views

Removing Gaussian Noise from a Signal to Get Minimum Value

I have a signal that has a minimum value that I'm trying to read. The issue I'm having is that the signal is spread out by gaussian noise. I have the signal at a lot of timesteps (and expect the ...
user avatar
  • 153
-1 votes
3 answers
99 views

Recover an OFDM signal affected by AWGN

I have an OFDM signal affected by AWGN, if I know the power of the gaussian noise, can I recover the OFDM signal? Can I perform deconvolution to recover the original signal? thanks
user avatar
0 votes
1 answer
105 views

Obtaining an IR from a sine sweep without the dry sweep

I have a test vinyl record that includes an exponential sine sweep as one of the tracks. I have a recording of the sweep through my turntable/mixer, and now I want to calculate the system IR. The ...
user avatar
  • 1
-1 votes
1 answer
93 views

Receiver function, frequency domain deconvolution not giving logic results

I'm working on some code for receiver function method in seismology. For anyone one not into the topic, it's just a deconvolution of two time series (seismograms). This can be done in the time domain ...
user avatar
5 votes
3 answers
234 views

Deconvolving a 1d Signal Using a Lookup Table

assuming I measure a signal that has different PSFs per position in time. for example: ...
user avatar
  • 586
0 votes
0 answers
196 views

How to water level deconvolve a noisy signal if i have a zeros and poles file?

Thanks for your time and help! I am working with Apollo project passive seismic experiment (PSE) data, and I have a large set of seismic records (on digital counts) and the corresponding file of poles ...
user avatar
0 votes
1 answer
69 views

Convolution that outputs a unit impulse

Im thinking whether any convolutional operation can output a unit impulse, an example to further explain: where a convolution between system $h[n]$ and unknown system $g[n]$ would output $\delta[n]$. ...
user avatar
0 votes
0 answers
16 views

Whats the Correct Approach to Estimate the PSF of a Moving Detector

In lab, I did a bunch of scans using a radiation source and a detector. My source emits a gaussian beam (whose dimensions I know), and my detector is a uniform 7mmx7mm square. These are stationary and ...
user avatar