Python is a high-level interpreted general programming language. Python has popular numerical and scientific libraries/packages, most notably numpy and scipy. The scipy.signal package is a powerful signal processing software collection. This tag covers the use of numpy, scipy, and other Python ...

learn more… | top users | synonyms

0
votes
0answers
10 views

Implementing capons spectral estimator in python

I am trying to implement capons high resolution frequency-wavenumber spectral estimator (mvdr) for use with data from an isotropic array. My results were unstable so I use diagonal loading. I am a bit ...
1
vote
0answers
8 views

(python+opencv) haar cascade training using cloud computing

I will train my own cascade, using python and opencv. I want it to be as precise as possible, so as I read it is common that training takes one or two weeks to compile. I want to use cloud computing ...
0
votes
1answer
23 views

Apply FIR filters to Extract Stego from WAV data

I appreciate anyone that takes a moment to help me with this problem. I've been banging my head against the keyboard for a while, searching forums and DSP tutorials and I can't figure this problem ...
0
votes
1answer
23 views

Confused with convolutions in scipy

I tried to write my own circular convolution function in python using the fact that for two signals $f$ and $g$ we have $$ \widehat{(f * g)} = f \cdot g $$ So I tried this ...
0
votes
1answer
42 views

Removing periodic noise from image using Fourier transform

I'm trying to get rid of some periodic flicker noise through post-processing of the recorded images. The reason for these artifacts is that the electronic rolling shutter of the camera reads each line ...
0
votes
1answer
35 views

Histogram matching (specification) in Python

I'm trying to implement an algorithm in which I first pad each row of the image with a fixed amount of new pixels in a certain range, apply Gaussian smoothing to the row cumulative histograms in ...
3
votes
1answer
106 views

Is it possible to hack camera calibration without having access to the camera?

Many times we have a stream of video to process without access to camera. Having access to camera matrix would be beneficial for various processing techniques. Is it possible to hack camera ...
0
votes
1answer
57 views

OpenCV Motion Analysis: Optical Flow vs Motion History

OpenCV Motion Analysis Documentation lists multiple techniques for motion tracking: Sparse Optical Flow | calcOpticalFlowPyrLK, Dense Optical Flow | calcOpticalFlowFarneback, and Motion History ...
0
votes
0answers
10 views

choosing how many DWT Decomposition levels for 1D

I am building application decomposition of 1D signal sampling frequency $f_s=256$, number of samples $N=512$ - how many levels if the above signal passed through LPF of $70 Hz$ $f_s=256$, number ...
1
vote
1answer
55 views

Sub pixel shift causes weird results

I am trying to shift an image on a subpixel level, with the method described here: ...
0
votes
1answer
41 views

The Nyquist-Theorm behind the fft

I saw a code in Python in the scipy.org (http://docs.scipy.org/doc/scipy-dev/reference/tutorial/fftpack.html) and found this: First question: Is it right that T = 1/800 <==> fs = 800 Hz fulfill ...
0
votes
1answer
25 views

Record and playback has echo/feedback

I didn't mess around with the code here, it was straight from their website, I just want to see if it would work. However, whenever I run it, it has feedback in it, sometimes I can't hear anything at ...
0
votes
1answer
43 views

Why do we need to normalize a WAV file before calculating the FFT?

I am trying to calculate the FFT of a signal stored in a WAV file using SciPy. I found a solution here, but it seems like we need to perform this step before the FFT: ...
0
votes
1answer
24 views

One representative time series signal from multiple time series data

I came across this basic problem in DSP which is not my area owf expertise but I consider this problem to be a common problem which others might be addressing. I have multiple time series data and I ...
0
votes
0answers
28 views

Feature extraction - DWT

I am new to DSP, I know a quite bit about Python. I have started to code to classify audio (wav files) ...
0
votes
0answers
21 views

Template matching results when the template is not on the original image

For testing purposes I'm building a bot that clicks on buttons in the screen. I have a folder with all the possible buttons that needs to be pressed and the application being tested only has one ...
2
votes
1answer
158 views

Applying filter in scipy.signal: Use lfilter or filtfilt?

I saw in a SO thread a suggestion to use filtfilt which performs backwards/forwards filtering instead of lfilter. What is the ...
0
votes
2answers
94 views

Butterworth filter in python

I'm trying to use a Butterworth filter in Python as described in this thread with these functions: ...
0
votes
1answer
30 views

Smooth X,Y,Z data with discontinuous, artificial jumps

I have some data $(X,Y,Z)$, which is a set of measurements $Z$ over a $2D$ space $X$,$Y$. The $Z$ data on this space is continuous except for some discontinuous jumps in certain domains of $X$,$Y$ ...
0
votes
0answers
52 views

Python implementation of multidimensional power spectral density with Welch method

I have done my best to write Welch method implementation for python for multidimensional time series and still in the case of one dimensional time series I am getting inconsistent response compared to ...
2
votes
0answers
81 views

Extracting the phase of a jittering sine wave using python

The background: I conducted an experiment where I rotated an object (with three-fold symmetry) many times and recorded the projection of the object. I should be able to tomographically reconstruct ...
0
votes
0answers
105 views

In a signal there is a peak, I want to find the start of the peak and the end of the peak

I want to the start time and end time for the peak shown below. The start and end are shown by the marked circles in the fig As there are number of such plots and want to know start and end time for ...
1
vote
0answers
81 views

GCC-PHAT cross correlation

I am trying to implement GCC-PHAT in python. The approach is similar to the following two links: link1 and link2 It seems the only difference between GCC-PHAT and normal cross-correlation using FFT ...
0
votes
3answers
64 views

float image to integer type image

I have a multispectral satellite image and I would like to "transform" each of the image bands into one with integer values. For example to stretch? the image between 0 and 255. Not sure how to ...
1
vote
1answer
64 views

correct for shift from butterworth filter?

I have implemented a python butterworth filter but it seems to shift the data a lot, the x location of my data is very important to me for my analysis. Is there a way to correct for this shift ...
0
votes
2answers
72 views

How to change the start point of resampling process?

I'm trying to resample a digital communication signal in SciPy or Matlab. Suppose the signal was sampled at 3 samples per second, and I want to resample it to 2 samples per second, the thing is I also ...
0
votes
0answers
84 views

Using Python with scipy.signal firwin, how do I know how many taps/coefficients should be used?

I need to alter an FIR filter which I know very little about. I have a signal which need to be altered, and looking at a Python script used for this task, I see that it uses Scipy's signal library, ...
0
votes
1answer
56 views

Advice on filtering a CW signal for processing

This question is related to my other recent questions. In the last question, I learned that I should filter about 4 * WPM (words per minute) around the carrier frequency if I want to recover Morse ...
0
votes
1answer
60 views

How to decode Morse code (CW)

I have two related questions on this topic. First, about how wide should the filter be around the carrier frequency? What is used to determine this? Second, in Python, I must use 'numtaps' when ...
0
votes
1answer
242 views

python; audio; detecting silence in audio signal

I have real-time audio application written in python, and I have question about detecting silence. What's the most efficient way to distinguish silence from non-silence, in an input audio signal? By ...
0
votes
2answers
44 views

Color artifacts in Fourier transformed image

I am working with images and Fourier transforms. I am trying to understand what might be causing some artifacts in my output image. I am starting with a 512x512, RGB image of Lenna. I FFT the ...
1
vote
1answer
86 views

Simulation and theory of first-order low pass filter

I am trying to teach myself some signal processing and control theory. I'm currently working through THIS tutorial. I've started an ipython notebook to keep my notes in. Unfortunately, something is ...
-1
votes
2answers
76 views

What is wrong with this method of detecting the presence of a signal?

I am trying to decode Morse code from an audio (wav) file. Here is the technique I am using to detect the presence of the signal: I have a window of size (sampling_freq/tone), where tone is the ...
0
votes
2answers
47 views

Amplifying part of a signal

If I want to amplify a certain frequency in a digital signal, can I just take the FFT, increase the real part of the value for the frequency I want (leaving the imaginary part alone), and then take ...
1
vote
2answers
103 views

Filtering a Morse code (CW, OOK) signal

I am building a morse code decoder in Python. I am using SciPy and Numpy for FFT, filter, etc. It is coming along ok, but the signal is noisy. I am extracting the dominant frequency from the FFT of ...
0
votes
0answers
49 views

Gettiing a rectangular waveform from OOK in Python

I have a wave (wav) sound file with on-off keying. I can read the signal value at any point in time (granularity of the sample rate) from an array. I want to create another array, the same size as ...
4
votes
2answers
258 views

Image Processing: Flatten a 3D Ball

There are many descriptions of how to turn a 2D image into a 3D one, however I want to do the opposite, in particular to a ball. As an example, consider the following ball: If only the (relatively) ...
1
vote
1answer
85 views

Filtering method that minimizes ringing in time domain, implementable in Python

I am acquiring data through an oscilloscope, which I am then filtering. I am doing all of this in Python. I need a bandpass filter with stop bands from 0 to .4MHz, and 4MHz and up. My pass band needs ...
3
votes
1answer
108 views

Confusion about result of FIR All-pass filter design

I try to design a FIR All-pass filter with random phase in the frequency domain. I am a bit confused by my result and am not sure if the reason is a programming error or a misconception about ...
1
vote
1answer
155 views

Transfer function determination from input and output data

I have some input and output data that I believe adequately includes excitation of the important dynamics of a system. I know it is at most a 4th-order transfer function. How can I identify the ...
1
vote
0answers
92 views

Signal leveling using scipy

Apologies if I'm using incorrect terms, and if this question's already been answered here. I have a data set which consists of a high DC signal modulated by AC signal (taken from custom-built ...
7
votes
3answers
340 views

Image processing coding

I was wondering which is the best language for image processing? I know MATLAB has a good library and user community for it. I havent checked it out, its there for python as well. What is the most ...
1
vote
0answers
120 views

How can I specify the scale parameter of the wavelet in PyWavelets (pywt)?

When computing wavelet transforms with PyWavelets (Discrete Wavelet Transform in Python), e.g. ...
0
votes
0answers
82 views

PSD Estimation: Welch's method with median averaging

I would like to use median averaging when making a PSD estimate via Welch's method. I notice that this is available as an option in R welchPSD, but not in SciPy or Matlab (in which the default is mean ...
0
votes
2answers
168 views

How can an FFT plot not be symmetrical?

From discussions all over the internet, I understand that in a 512 sample point fft, only the first 256 contain relevant information provided that my signal is real. The other half contains complex ...
0
votes
2answers
164 views

How to interpret results of Fourier transform by using Python

Connected with this topic , I would like to have some clarifications about the Fourier Transform and the Python tools that make it. Having my data set whose plot is showed in the above link, my array ...
1
vote
1answer
302 views

How to find the inflection point in a noisy curve?

I have the the noisy curve defined by numpy 2D array: As you can see, it has the first flat segment, then rise, peak and decay phases. I need to find the starting point of the rise phase, marked ...
0
votes
1answer
2k views

Plot real-time FFT using matplotlib

I'm trying to plot fft in python. I use pyalsaaudio for capturing audio in PCM (S16_LE) format. I use the ion() and draw() functions in matplotlib to have the fft plotted in real time. This is the ...
0
votes
0answers
13 views

Reference Pattern in Background

I want to apply a nonlinear matched filter to an image to find the reference pattern in the background. The image is just a white box in the foreground with a black background. I've included the code ...
1
vote
0answers
59 views

Lowpass then Inverse filtering in python

I wrote a simple lowpass filter in python to run against lena. Now I'd like to add Gaussian noise to the lowpass filtered data and then run an inverse filter against the lowpass and try to get the ...