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
17 votes
13 answers
36k views

Deconvolution of 1D Signals Blurred by a Gaussian Kernel

I have convolved a random signal with a a Gaussian and added noise (Poisson noise in this case) to generate a noisy signal. Now I would like to deconvolve this noisy signal to extract the original ...
  • 179
16 votes
2 answers
4k views

How to Deduce a Linear System's Impulse Response from a Set of Input and Output Signals?

I want to know how to solve those types of problems.. is it by inspection ? Consider the linear system below. When the inputs to the system $x_1[n]$, $x_2[n]$ and $x_3[n]$, the responses of the ...
  • 169
14 votes
2 answers
2k views

Deconvolution - Richardson Lucy vs. Wiener Filter

I am studying some deconvolution techniques, In order to remove motion blur, like: Richardson-Lucy Wiener Are there any pros / cons of using one versus another? For example which are the pros / cons ...
  • 305
12 votes
2 answers
415 views

Deconvolution Question on Article "Deriving Intrinsic Images from Image Sequences" by Yair Weiss

there are n derivative filters: $f_i$, and denote $f_i^r$ as $f_i$'s reverse filter such that $$f_i(x,y)=f_i^r(-x, -y)$$ $r_i, f_i$ given, to find $r$ from the equations: $$f_i * r = r_i, (1 \leq i \...
12 votes
1 answer
3k views

Tutorial on 2nd generation wavelets (with lifting)?

For some denoising and deconvolution experiments, I'd like to apply a 2nd generation wavelet transform (using lifting steps) to images. I know that there are several implementations available, but ...
  • 3,887
11 votes
4 answers
4k views

Deconvolution by Convolution

This is now a second time I am attempting to ask this very important but simple question here. What I want to know is can you do deconvolution by convolving a signal. It is often stated that, for ...
  • 141
9 votes
3 answers
1k views

Deconvolution of Synthetic 1D Signals - How To?

I convolved a square wave with a Gaussian wave using linear convolution. Can I get the original square wave back by deconvolving my output with the Gaussian function? I took the FFT of both signals, ...
9 votes
1 answer
5k views

Estimate the Filter Coefficients of 1D Filtration (Convolution)

I have an output signal $y$ which is an input signal $x$ convolved $\star$ with an impulse response function $h$ with some added noise $n$ : $$y(t) = h(t) \star x(t) + n(t)$$ I know the input signal $...
  • 442
9 votes
1 answer
5k views

1D Deconvolution with Gaussian Kernel (MATLAB)

Suppose that I know the output and the transfer functions of a system and I would like to calculate the input function using deconvolution. To get a grasp of the idea I have created a simple ...
9 votes
4 answers
1k 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), ...
8 votes
2 answers
6k views

Using the Inverse Filter to Correct a Spatially Convolved Image (Deconvolution)

As part of a homework assignment, we are implementing the Inverse Filter. Degrade an image then recover with an Inverse Filter. I convolve the image in the spatial domain with a 5x5 box filter. I FFT ...
8 votes
3 answers
522 views

Deconvolution (Linear Convolution) with an Under Determined System of Equations?

If I have a measured signal $\mathbf{y}$ which is the result of the true signal $\mathbf{x}$ convolved with another function (linear and not circular convolution), I always seem to get an ...
  • 259
8 votes
2 answers
3k views

How Is the Formula for the Wiener Deconvolution Derived?

Wikipedia - Wiener Deconvolution shows this formula: $$ \ G(f) = \frac{H^*(f)S(f)}{ |H(f)|^2 S(f) + N(f) } $$ Could someone explain the derivation? Specifically, where does the squaring ($|H(f)|^2$) ...
  • 81
8 votes
1 answer
870 views

How to Use the DFT (FFT) to Solve a Least Squares Regularization Problem (Inverse Problem)?

Let $X$ and $K$ be an image and a Point Spread Function (PSF), respectively. The blurred image $B$ is obtained as follows $$B = X * K$$ I want to solve the following general regularization problem $$\...
8 votes
1 answer
2k views

Use MATLAB to Restore a Signal from a Given Degraded Signal Using Tikhonov Regularization

Anyone could share how to develop an application in MATLAB to restore the signal from a given degraded signal using Tikhonov regularization i.e restoring the signal $f$ via solving $$ \min || g - f ...
7 votes
3 answers
1k 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 ...
  • 73
7 votes
2 answers
3k views

What's the Difference between Convolution Kernel and Point Spread Function in the Context of Image Convolution?

When you use the deconvolution method to make the blurry image sharper, you will have to estimate the Point Spread Function. Is there a difference between this PSF and an image kernel? Second ...
  • 71
7 votes
2 answers
838 views

What Are Less Computationally Demanding Alternatives to the Viterbi Decoder?

What are less computationally demanding alternatives to the Viterbi Decoder? Ideally what I would like is a list of the most commonly used approximate methods, along with brief pros and cons.
  • 293
7 votes
1 answer
375 views

Why Sparse Priors Like Total Variation Opts to Concentrate Derivatives at a Small Number of Pixels?

When performing image deconvolution (deblurring), people often make use of priors to get rid of the illness of the problem. One very common prior is total variation, a sparse prior. Compared to ...
  • 557
7 votes
2 answers
3k views

Differences Using Maximum Likelihood or Maximum a Posteriori for Deconvolution / Deblur

Are there any differences if you use Maximum Likelihood or Maximum a Posteriori to estimate the Point Spread Function for image deconvolution?
  • 131
7 votes
1 answer
160 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.
  • 752
7 votes
2 answers
338 views

What Is the Correct Way to Apply Richarson Lucy Deconvolution to Luminance Data?

My question concerns the Richardson–Lucy deconvolution algorithm, which is described in Richardson's original paper - William Hadley Richardson - Bayesian Based Iterative Method of Image Restoration. ...
  • 71
7 votes
2 answers
700 views

Finding the Best Gaussian Smoothing Kernel to Minimize the Discrepancy Between Two Images

Suppose we have two grayscale images, $A$ and $B$. $A$ and $B$ very strongly resemble each other, such that the mean of the absolute difference $\lvert A - B\rvert$ is fairly small. Suppose further ...
  • 173
7 votes
2 answers
457 views

Using reference objects to estimate the point spread function?

I have a well-defined object and a clear image matrix of it. In subsequent frames the object moves, causing motion blur. I want to use the object as a reference to "guide" the deconvolution and ...
  • 806
6 votes
3 answers
237 views

Estimating the Signal by Deconvolution with a Prior on the Filter Coefficients and the Signal Samples

Assume I have signal $y[n]$ which is a result of convolution between channel $h[n]$ and signal $x[n]$. which means: $$y[n] = h[n] \ast x[n]$$ where $\ast$ is the convolution operation The signal $...
  • 640
6 votes
2 answers
353 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 ...
  • 63
6 votes
2 answers
195 views

Deconvolution of an Image Acquired by a Square Uniform Detector

So, I acquired some images by scanning a radiation source with a square detector like in the following gif. Where the dashed grid represents reality, the 3x3 square my detector, and the 4x4 my ...
6 votes
1 answer
4k views

Relationship between Discrete Deconvolution and Toeplitz Matrices

I have 2 vectors, $a$ & $c$, both of length M. I know they are related by $a*b=c$. My goal is to recover $b$. Obviously $b=$deconv$(c,a)$. I am only interested in the first M elements of the ...
6 votes
1 answer
182 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. ...
  • 415
6 votes
4 answers
323 views

Estimate Filter Coefficients from the Result of Linear Convolution with a Known Signal

If I have samples of input say x(1:500) and it passes through FIR filter with 9 taps and some unknown coefficients. The output y(1:508) is also known. The aim is to estimate the filter response in ...
  • 61
6 votes
1 answer
502 views

Deconvolution Using Response to an Heavy Side

I'm measuring a "charge" signal in function of time from an amplifier. Here is a measured signal (x-axis is the time in some arbitrary units, y-axis is the charge in ADU): I would like to get the "...
6 votes
1 answer
70 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 ...
6 votes
1 answer
408 views

What Are the Types of Deconvolution?

I am totally new to image processing and wanted to ask you if you could confirm what I understood. It is about deconvolution: From what I read we find 2 main types of deconvolution: 1. Analytical <...
  • 123
6 votes
2 answers
485 views

How to Create a Convolution Matrix with a Variable Condition Number (CN)

I want to know the performance of a deconvolution algorithm with different CN, so I'm convolving my signal with different convolution matrices(different CNs) and then applying the deconvolution ...
6 votes
1 answer
111 views

Why Does Image Deconvolution Still Work with Image without Sharp Edges?

Why does image deconvolution still work with image without sharp edges? Take this image for example:
  • 131
6 votes
1 answer
489 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 ...
6 votes
1 answer
488 views

How to Select Point Spread Function Empirically for Image Deconvolution?

When the captured image is blurring, one way of obtaining a clear image is via image deconvolution technique. In order to perform deconvolution successfully, usually we need to pay attention to the ...
  • 557
6 votes
1 answer
66 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 ...
  • 63
6 votes
1 answer
845 views

Deblurring algorithm to precede thresholding - speed over accuracy

I'm writing an app that recognizes Sudoku puzzles from a camera input. I'd like to remove camera blur from the images to improve recognition. Here is an example image: Since I'm processing a ...
  • 113
5 votes
1 answer
3k views

Compensating Loudspeaker frequency response in an audio signal

I have been working on a project in which I was required to work on the audio signals recorded from the loudspeaker kept in front of a filter. So, to simply explain it: $$\boxed{\rm LoudSpeaker} \...
  • 51
5 votes
4 answers
1k views

Can every type of linear filter be modelled by a convolution?

I have an input time series going through a filter that creates another time series as output. If I assume in first approximation that my filter is linear, does it necessarily mean that I can model ...
5 votes
2 answers
362 views

Inverse Problem / Deconvolution with Pink Noise

Hi I dived somewhat into deconvolution of systems which can be described as: $s(t) = o(t) * h(t) + n(t)$ where $s$ is my measured 1D time resolved signal, $o$ is the original signal $h$ is the ...
  • 53
5 votes
3 answers
744 views

Recover Filter Coefficients from Filtered Noise

I have a digital signal which may be represented as a pulse noise source filtered with an FIR (finite impulse response) filter. Suppose that the noise consists of discrete pulses (nonzero samples ...
  • 200
5 votes
3 answers
242 views

Deconvolving a 1d Signal Using a Lookup Table

assuming I measure a signal that has different PSFs per position in time. for example: ...
  • 616
5 votes
2 answers
326 views

Finding the Width of a Point Spread Function (PSF)

I have a binary image that has been convolved with the following PSF: $$\mathrm{PSF}(x,y)=\cosh^{-2}\left(\frac{\sqrt{x^2+y^2}}{w}\right)$$ where $w$ is unknown. The image binary image generally ...
  • 153
5 votes
2 answers
1k 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? ...
  • 752
5 votes
2 answers
419 views

The Gradient / Derivative of Least Squares of 2D Image Convolution

Given the objective function: $$ \frac{1}{2} {\left\| h \ast x - y \right\|}_{2}^{2} $$ Where $ h $ is the 2D convolution kernel and $ x $ is the 2D convolution image and $ y $ is a given 2D image. ...
  • 42.4k
5 votes
1 answer
107 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 ...
  • 153
5 votes
3 answers
2k views

Deconvolution Using Complex Division in The Frequency Domain

Consider these two signals: a = [1 1 0 0 0 0 0 0] b = [1 0 1 0 0 0 0 0] their convolution is c = a * b = [1 1 1 1 0 0 0 0] ...
5 votes
1 answer
96 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 ...