Stack Exchange Network

Stack Exchange network consists of 175 Q&A communities including Stack Overflow, the largest, most trusted online community for developers to learn, share their knowledge, and build their careers.

Visit Stack Exchange

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).

5
votes
3answers
215 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, ...
0
votes
1answer
51 views

Blind Deconvolution

I have a signal X which is a sinusoid with asymmetrical clipping on the positive half of the waveform. The clipping characteristic is “soft”. I don’t know ahead of time either the clipping threshold ...
3
votes
0answers
52 views

The tail of scipy deconvolve

I have a dataset for the function $g(t) = \int_{0}^{\infty}f(t-\tau)h(\tau)d\tau$ I would like to deconvolve. My assumptions are The signal $f(t) = 0$ before the start of observation, that is for $t &...
-1
votes
1answer
32 views

Deconvolution of 2 vectors (1 know + 1 unknown)

I am currently trying to deconvolute 2 vectors (a & b) from 1 (c). Actually, I have access to the recorded data of (a) & (c) but not (b). All are signal vs time with signal totally random. I ...
1
vote
1answer
208 views

Deconvolution of a 1D Signal with Known Kernel (Square Wave)

I have a signal measured from a radiation detector in a narrow beam of radiation. The peaks I get are quasi-gaussian in shape, see attached picture. The signal is not a function of time, rather a ...
1
vote
3answers
176 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] ...
0
votes
0answers
97 views

extracting a room response [duplicate]

So I have a wav file containing a record of a sine sweep played on my pc and recorded by a mic. And I want to extract the room impulse of it, I've looked online but I couldn't find any clear answers. ...
0
votes
1answer
70 views

Wiener deconvolution in cpp giving back the same image

I was trying to deblur a noisy image using wiener deconvolution. I found this code which added noise to an image and removed it as well. Modifying this code only i tried to implement the exact formula ...
1
vote
0answers
139 views

scipy.signal.deconvolve returns an array of NaNs [closed]

I am trying to implement a deconvolution-based event detection algorithm in python, but scipy.signal.deconvolve doesn't seem to work in my case. Here is a basic example: ...
2
votes
2answers
111 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 ...
2
votes
0answers
96 views

Is it possible to “equalise” a signal by deconvolving the impulse response of the room in which it is to be played?

I am using a test sweep with a flat power spectrum and linear group delay (Optimized Aoshima's Time Stretched Pulse) to measure a room's frequency response. Having obtained the impulse response of the ...
1
vote
1answer
60 views

Help with deconvolving photon reads distributed across neighboring pixels

I'm building a radiation detector that collects photons in CCD pixels and we can relate the energy of the photon from the intensity of the pixel. To test the detector, I took the following spectrum of ...
1
vote
1answer
272 views

1D convolution and deconvolution using FFT

The task: there is some original signal, and there is some response function. I need to convolve them using FFT and then do deconvolution to restore original signal. The task graphical illustration ( ...
1
vote
1answer
94 views

Problems with IFFT not being symmetrical

I have two signals, a measurement and a reference which I have performed an FFT on. They have both been windowed with a Hanning window, and now I would like to deconvolve these to get the impulse ...
0
votes
2answers
66 views

finding the width of a point spread function

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 ...
1
vote
2answers
185 views

Invertibility of Room Impulse Response: Reproducing Research Paper

I have been trying to reproduce this paper¹. Few things which are unclear to me. The paper talks about finding whether a given Room Impulse Response(RIR) is invertible or not based on Nyquist plot. ...
2
votes
1answer
259 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 blur image $B$ is obtained as follows $$B = X * K$$ I want to solve the following general regularization problem $$...
4
votes
1answer
102 views

Underdetermined deconvolution of windowed output

Consider a discrete 'blurred' output $h[t]$ given by the convolution of filter $f[t]$ and signal $g[t]$. This question considers recovering $g[t]$ from a window (subset) of $h[t]$. This causes the ...
4
votes
2answers
129 views

Deconvolving known impulse response from a sampled noisy signal

I am interested in measuring a signal with significant energy content up to 1-2 kHz. However, I am only able to sample the signal after it has passed through a first-order lowpass filter of known ...
1
vote
2answers
95 views

Question related to Deconvolution and Fourier Transform

Currently I am completely new to Signal Processing and I am about to kick off a project related to the system identification. The general idea is given a pair of sent and received signals, I have to ...
1
vote
1answer
574 views

How to properly match bit growth of a FIR filter with a reference value

I am trying to implement a chain of CIC/FIR filters on an ZYNQ FPGA. Using the Xilinx FIR compiler works fine so far but I am unable to properly get all the math. At the moment I have 2 chains of ...
1
vote
1answer
217 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 <...
2
votes
0answers
179 views

NMR spectrum deconvolution on Matlab

I am doing an NMR experiment by applying an rf signal to a sample which is mounted near a magnet and in the end I get an FID (free induction decay) and do Fourier Transform to this FID . So we have a ...
0
votes
1answer
42 views

Is this a robust method for analysing correctness of a filtered signal? Is there a more robust method?

I have a recorded signal which contains several frequency components with stochastic phase and amplitude changes and want to extract the true signal with the other frequency components and noise ...
1
vote
1answer
128 views

What is the optimal Wiener-filter for this purpose?

Suppose I want to measure something with a sensor which in response to $u(t)$ produces an output signal in form of $v_{out}(t) = (1-e^{-t/ \tau})u(t)$, where $u(t)$ is the Heaviside-function and $\tau=...
0
votes
1answer
84 views

Is it possible to do deconvolution with two data sets that have different sampling rates?

I have some terahertz spectroscopy time series data, a reference set with 2048 data points taken every 0.0521 picoseconds, and the sample data set with 544 data points taken every 0.0781 picoseconds. ...
0
votes
2answers
897 views

What is the difference between image restoration and image reconstruction?

I am new to image processing. I don't know whether this is the right place to ask, but what is the difference between image restoration and image reconstruction?
0
votes
3answers
1k views

Convolution and Deconvolution in C

I understand the theoretical foundations of convolution, but now that I'm trying to program it I'm having some issues conceptually. Say I have two blocks of 64 audio samples each. I also have access ...
0
votes
0answers
234 views

Deconvolve using divsion of FFT of shifted signal in time domain

I have FFT of two signals. Y=120 , Y1=80 Hz. to convolve in time domain I can convolve them using their ffts. as ...
3
votes
1answer
714 views

How to implement the $\tt conv(x,h)$ in MATLAB without using loop?

I'm trying to calculate convolution of two given vectors in MATLAB without using loop, and of course without the function conv itself, but I can't remove the last ...
1
vote
3answers
686 views

Deconvolution and Polynomial division

I found this comment inside MATLAB's deconv.m function Deconvolution and polynomial division are the same operations as a digital filter's impulse response $B(z)/A(z)$. What does this statement ...
0
votes
2answers
198 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 ...
3
votes
3answers
1k views

Sweep signal + microphone-recorded signal => Impulse response

I played this signal A (a 20Hz to 20000Hz sinusoidal sweep in 10 seconds) with a studio monitor speaker in a big church, and I recorded the result B with good microphones. The result is very reverb-...
0
votes
1answer
339 views

python rewrite lfilter (iir) with for

It's a beginner question, but useful to users from python - signal.lfilter, I was using lfilter from Find reverse one pole ...
2
votes
2answers
372 views

Find reverse one pole lowpass filter

I need to find a filter that revert the one pole filter of the current signal, a function using Python (or MATLAB) scipy.signal.filtfilt or ...
0
votes
0answers
124 views

noise in audio deconvolution result

I'm testing some code to perform deconvolution of two audio signals to recover the impulse response. Presently, as part of testing, I am simply deconvovolving two identical signals. I have the ...
0
votes
0answers
57 views

deconvolution challenge

OK hive mind... i have two identical fixed-size signal buffers, each containing 4 channels of audio. I perform a 2d REAL to COMPLEX FFTW_ESTIMATE on both, resulting in two identical spectra of the ...
0
votes
1answer
46 views

Beginner level: deconvolving Gaussian noise to find initial PDF

noob question. An initial process has an unknown pdf, P(X), which is then subjected to additive Gaussian noise. Repeated sampling of the sum is performed so that $y=x+noise$ is repeatedly ...
3
votes
2answers
87 views

Modelling an existing DSP filter by analysing the differences between its inputs and outputs

A DSP filter produces always the same output. $$\boxed{\textrm{wav file}}{\longrightarrow}\boxed{\textrm{DSP}}{\longrightarrow}\boxed{\textrm{wav file(modified)}}$$ The DSP is a black box as I don't ...
1
vote
0answers
173 views

Can we use SVD to solve single channel deconvolution problem?

I have seen using singular value decomposition (SVD) to solve deconvolution problem for example truncated SVD (TSVD) . It appears there is also a connection between Tikhonov regularization and SVD. My ...
0
votes
2answers
496 views

What would be the outcome of the convolution of two images? [closed]

Lena Mask I have written a C# code which produces the following output. Is my output correct? Am I padding in the correct way?
0
votes
1answer
562 views

How does scipy.signal.deconvolve work?

The deconvolve(signal, divisor) function of scipy ...
1
vote
1answer
187 views

Signal enhancing algorithm (deconvolution??)

I need an algorithm that will detect the correct peak location in a "signal" (actually a series of data points). The signal is generated by a machine. The signal is nothing similar to an electrical ...
1
vote
0answers
67 views

Why is the term deconvolution used more for signals and not so much (or at all) for systems?

Wikipedia defines mathematical deconvolution here, and with the examples given and my experience, what I've read over the years is that deconvolution is used to determine an input signal provided a ...
1
vote
1answer
116 views

How to deconvolve matrix when a model of noise exists?

I have a matrix of N rows of time-series data. There is a specific noise contaminating measurement of the data that I have some information about. The noise in the data can be modeled as a poisson ...
4
votes
1answer
805 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} \...
1
vote
0answers
255 views

how to get an input test signal for convolution/ distortion method

i am trying to do some experiments with Total Harmonic Distortion. According to ...
1
vote
1answer
225 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 ...
0
votes
2answers
124 views

Why do we assume the matrix of impulse responses to be a Toeplitz matrix during deconvolution

$y(n)$ = output signal $x(n)$ = input signal $\mathbf H$ = system response as a toeplitz matrix $$\mathbf H = \begin{bmatrix}h(0)&&&\\h(1)&h(0)&&\\h(2)&h(1)&h(0)&\\...
1
vote
0answers
137 views

How to use deconv() instead of roots() on MATLAB to find roots to a polynomial [closed]

I recently read that for polynomials of degree 5 or more, when executed with the roots() command on MATLAB produces an error. The documentation said as follows: As a substitute, using the deconv() in-...