http://www.k-wave.org/forum/topic/how-to-control-the-source-frequency-in-photoacoustic-image-reconstruction Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 01:07:29 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/how-to-control-the-source-frequency-in-photoacoustic-image-reconstruction#post-6451 Sun, 13 May 2018 19:30:50 +0000 Bradley Treeby 6451@http://www.k-wave.org/forum/ <p>Hi Hui,</p> <p>In photoacoustics, it is generally assumed that both stress confinement and thermal confinement are satisfied. This means that the wave propagation can be modelled as an initial value problem. In this case, the intrinsic bandwidth of the source will depend only on the source geometry. In an experiment, the measured bandwidth will also depend on acoustic absorption (which acts to proportionally reduce high-frequencies), and the bandwidth of the detector.</p> <p>In simulation, the bandwidth is also limited by resolution of the grid, where the maximum frequency is given by <code>c/(2*dx)</code>. To avoid oscillations in the simulated pressure field due to the bandlimited interpolant, k-Wave spatially smooths the initial pressure distribution (see the <a href="http://www.k-wave.org/documentation/example_na_source_smoothing.php">Source Smoothing Example</a>). Thus, numerically, the bandwidth will depend on:<br /> (1) the shape of the source,<br /> (2) the grid resolution,<br /> (3) the smoothing applied by k-Wave (which is related to 1), and<br /> (4) any acoustic absorption in the simulation.</p> <p>If stress confinement is not satisfied, i.e., if the acoustic wave starts to travel some appreciable distance while the optical energy is being deposited, then the bandwidth will also depend on the time profile of the laser energy. This can be modelled in k-Wave using a time-varying source, e.g., in the shape of a Gaussian. If you want your simulations to be consistent, the spatial smoothing can also be applied when using a time-varying source. This is explained in the <a href="http://www.k-wave.org/documentation/example_na_filtering.php">Filtering A Delta Function Input Signal Example</a>. </p> <p>The situation in diagnostic ultrasound is quite different as you mention, and depends on the bandwidth of the transducer and driving conditions, etc.</p> <p>Hope that helps,</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/how-to-control-the-source-frequency-in-photoacoustic-image-reconstruction#post-6429 Thu, 19 Apr 2018 11:35:13 +0000 HuiLIN 6429@http://www.k-wave.org/forum/ <p>Dear Bradley,</p> <p> I have studied the examples in Photoacoustic Image Reconstruction and Diagnostic Ultrasound Simulation. I am wondering the differences between them?<br /> 1) There are clear source frequency (i.e. tone_burst_freq) and pulse loop count ( tone_burst_cycles) in US examples, however there are no these in PR examples.<br /> 2) There are FFT and TR Image Reconstruction processes in PR examples, but there are no these in US examples.</p> <p> I am wondering whether the first point means the laser (T ~ ns) in Photoacoustic is a continuous beam, whereas the microwave (T ~ us) in Ultrasound is pulse beam? Could the user change source frequency in Photoacoustic example? And how?<br /> The second question is, could we use FFT and TR Image Reconstruction processes in US examples? It seems that US image is pulse and dynamic, and it cannot use the smooth process required by Image Reconstruction. How to solve this problem?</p> <p> Excuse my stupid questions if there are some common senses. I am novice and aim to link the pulse beam with the image reconstruction.<br /> Thank your time a lot. </p> <p> Hui </p>