http://www.k-wave.org/forum/topic/frequency-content-signal-processing-in-fft-reconstruction Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 06:48:01 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/frequency-content-signal-processing-in-fft-reconstruction#post-4837 Tue, 11 Nov 2014 19:24:47 +0000 Bradley Treeby 4837@http://www.k-wave.org/forum/ <p>Hi Bill,</p> <p>What you are seeing is the combination of two classical problems when using linear piezoelectric arrays for photoacoustic imaging. </p> <p>The first is the limited aperture problem. This arises because the finite size of the imaging aperture means you will not detect the acoustic waves travelling in all directions from your photoacoustic source. This will result in arc-like artefacts in the reconstructed image, which you can see in the <a href="http://www.k-wave.org/documentation/example_pr_2D_fft_line_sensor.php">2D FFT Reconstruction For A Line Sensor Example</a>. A more detailed explanation of the limited view problem is given <a href="http://scitation.aip.org/content/aapm/journal/medphys/31/4/10.1118/1.1644531">here</a>.</p> <p>The second problem is the use of band-limited detectors. As photoacoustic waves are inherently broadband, you will lose some high and low frequency information when using narrowband piezoelectric detectors. This results in positive and negative oscillations in the reconstructed image, which you can see in the <a href="http://www.k-wave.org/documentation/example_pr_2D_tr_bandlimited_sensors.php">Image Reconstruction With Bandlimited Sensors Example</a>. Deconvolution <em>might</em> help, but you can't correct for frequencies that simply haven't been recorded. You'll also need to be careful you don't amplify noise.</p> <p>Regarding increasing the resolution of your image, you could use nearest neighbour interpolation to spatially upsample your data before calling <code>kspaceLineRecon</code>.</p> <p>Hope that helps,</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/frequency-content-signal-processing-in-fft-reconstruction#post-4831 Thu, 06 Nov 2014 20:30:44 +0000 wcvogt 4831@http://www.k-wave.org/forum/ <p>Hello,</p> <p>I am trying to use the FFT method kspaceLineRecon to reconstruct images from a 128-element linear array with 4-11 MHz bandwidth, sampled at 31.25 MHz. However, absorptive objects in reconstructed images appear not as expected point targets but as multiple bands or arcs with oscillating intensity between positive and negative values, especially along the axial (time-associated) direction.</p> <p>My two hypotheses as to where this artifact comes from are:<br /> 1) Insufficient grid density. Is there a way to adjust the grid density for the FFT methods? The grid seems hard-coded based on the dataset size. </p> <p>2) The oscillations appear similar to the temporal impulse response of our array as measured by directly illuminating it. Is impulse deconvolution critical for PAT? What signal processing is necessary before inputting data to the FFT method? </p> <p>Thanks,</p> <p>-Bill Vogt, Ph.D<br /> Postdoctoral Fellow<br /> US Food and Drug Administration<br /> Silver Spring MD </p>