Escribo una aplicación de procesamiento de imágenes que necesita hacer varias cosas y tiene que hacerlo lo más posible en tiempo real. La adquisición de los datos y su procesamiento se ejecuta en procesos separados (principalmente por motivos de rendimiento). Los datos en sí son bastante grandes (imágenes en escala de grises de 2MPix de […]
Estoy tratando de estimar la PSD de la variabilidad de la frecuencia cardíaca de una señal de ECG. Para probar mi código, extraje el intervalo RR de la base de datos de ECG de Fantasia . He extraído la señal que se puede acceder aquí . Para calcular el PSD, estoy usando el método de […]
Estoy tratando de entender la teoría del FTT y la convolución (correlación cruzada) y por esa razón he creado el siguiente código para entenderlo. El código es Matlab / Octave, sin embargo, también podría hacerlo en Python. En 1D: x = [5 6 8 2 5]; y = [6 -1 3 5 1]; x1 = […]
Miré en muchos ejemplos de scipy.fft y numpy.fft. Específicamente este ejemplo, el Análisis de Frecuencia de Scipy / Numpy FFT es muy similar a lo que quiero hacer. Por lo tanto, utilicé el mismo posicionamiento de subttwig y todo parece muy similar. Quiero importar datos de un archivo, que contiene una sola columna para que […]
Analicé los datos de sunspots.dat (abajo) usando fft, que es un ejemplo clásico en esta área. Obtuve resultados de FFT en partes reales e imaginarias. Luego intenté usar estos coeficientes (los primeros 20) para recrear los datos siguiendo la fórmula de la transformada de Fourier. Pensando partes reales corresponden a a_n e imaginarios a b_n, […]
Estoy tratando de controlar la funcionalidad fft de Python, y una de las cosas extrañas con las que me he topado es que el teorema de Parseval no parece aplicarse, ya que da una diferencia de aproximadamente 50 ahora, mientras que debería serlo. 0. import numpy as np import matplotlib.pyplot as plt import scipy.fftpack as […]
Obtuve este código de Python para el cálculo FFT de una señal de sonido: from math import * from cmath import exp, pi def fft(x): N = len(x) if N 0: print 20 * log10(sqr) # …in decibels else: print “-INF” Obtuve los siguientes resultados: SIN función de ventana (comentada): -20.3017238269 -16.9192604108 -12.5089302395 -8.97999530657 -5.96033201086 […]
Tengo varios puntos (coordenadas x, y, z) en un cuadro 3D con masas asociadas. Quiero dibujar un histogtwig de la densidad de masa que se encuentra en las esferas de un radio R dado. He escrito un código que, siempre que no haya cometido ningún error que pueda tener, funciona de la siguiente manera: Mi […]
Este es un seguimiento de una pregunta respondida que hice y que se puede encontrar aquí . Tengo varios puntos (coordenadas x, y, z) en un cuadro 3D con masas asociadas. Quiero dibujar un histogtwig de la densidad de masa que se encuentra en las esferas de un radio R dado. La idea es calcular […]
Tengo un lector de archivos que lee n bytes de un archivo y devuelve una cadena de caracteres que representan esos datos (binarios). Quiero leer n bytes en una gran cantidad de números y ejecutar una FFT en ella, pero tengo problemas para crear una matriz a partir de una cadena. Un par de líneas […]