Adaptive Filter | ADPCM | ADSP | ADSP-2181 | Aliasing | AMR | Anti-Aliasing | ARMA | Autocorrelation | AutoCovariance | Beamforming | Bessel | Blackfin | Butterworth | C6713 | CCS | Chebyshev | CIC Filter | Circular Convolution | Code Composer Studio | Comb Filter | Compression | Convolution | Cross Correlation | DCT | Decimation | Deconvolution | Demodulation | DM642 | DSP Boards | DSP/BIOS | DTMF | Echo Cancellation | Equalization | Equalizer | ETSI | EZLITE (Ez-kit Lite) | FFT | FFTW | FIR Filter | Fixed Point | FSK | G.711 | G.723 | G.729 | Gaussian Noise | Goertzel | GPIO | Hilbert Transform | IFFT | IIR Filter | Interpolation | Invariance | JTAG | Kalman | Laplace Transform | Levinson | LPC | McBSP | MIPS | Modulation | MPEG | Multirate | Notch Filter | Nyquist | OFDM | Oversampling | Pink Noise | Pitch | PLL | Polyphase | QAM | QDMA | Quantization | Quantizer | Radar | Random Noise | Reed Solomon | Remez | Resampling | RTDX | Sampling | Sharc | TI C6711 | Undersampling | Viterbi | Wavelets | White Noise | Wiener Filter | Windowing | XDS510PP | Z Transform

Hi,
 I am new to this programming world, and i really need some serious help.
I have done a program in which i design a filter in frequency domain, then
i take a CT image, which is a 3d image, read each voxel and make a 1d
array. then i equal the number of samples for both the data set. so i have
equal samples. then i do convolution using fftw3 library. but what happens
is that i am doing circular convolution, I want help to do this as a
linear convolution. please help me with this. waiting for your replies.
any other data you need please let me know.

thank you.

______________________________
New DSP Code Snippets Section now Live.   Learn more about the reward program for contributors here.

usual process of getting the effect of linear convolution.

if x is data and h is filter

M = data length
N = filter length

then zero pad x with N-1 zeros
zero pad h with M-1 zeros

perform FFT of both x and h with length >= M+N-1

do multiplication in complex domain to get the effect of filter.

perform IFFT and get back y(n)

hope this helps,

regards
bharat pathak

Arithos Designs
www.Arithos.com

%--- DSP design consultancy and training company ------------








>Hi,
> I am new to this programming world, and i really need some serious
help.
>I have done a program in which i design a filter in frequency domain,
then
>i take a CT image, which is a 3d image, read each voxel and make a 1d
>array. then i equal the number of samples for both the data set. so i
have
>equal samples. then i do convolution using fftw3 library. but what
happens
>is that i am doing circular convolution, I want help to do this as a
>linear convolution. please help me with this. waiting for your replies.
>any other data you need please let me know.
>
>thank you.
>
>
>

______________________________
New DSP Code Snippets Section now Live.   Learn more about the reward program for contributors here.