http://www.k-wave.org/forum/topic/transducer-design Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 00:17:36 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/transducer-design#post-348 Sat, 03 Mar 2012 00:50:39 +0000 Bradley Treeby 348@http://www.k-wave.org/forum/ <p>Hi Hussain,</p> <p>1. Do you mean the entire domain of interest is 5 x 5 x 5 cm, or there is a small optically absorbing sub-region that you are interested in? </p> <p>In 3D, <code>source.p0</code> is simply a 3D matrix which describes the initial pressure (in units of Pa) at each grid point within the domain. You can create this matrix and assign the grid point values (i.e., matrix elements) in the same way as you would create any other matrix in MATLAB. In your case, as you say, you can obtain this distribution from product of the fluence, absorption coefficient, and Gruneisen parameter at each grid point. </p> <p>The size of certain objects within the domain (as represented by a certain number of grid points or matrix elements) can be calculated using the values for <code>kgrid.dx</code>, <code>kgrid.dy</code>, and <code>kgrid.dz</code>.</p> <p>2. The <a href="http://www.k-wave.org/documentation/gaussianFilter.php">gaussianFilter</a> function is fairly straightforward to use. Simply choose a value for the center frequency and bandwidth (as a percentage of the centre frequency). For example:</p> <p><code>sensor_data_filtered = gaussianFilter(sensor_data, 1/kgrid.dt, center_freq, bandwidth);</code></p> <p>You can then apply this filter for different values of <code>center_freq</code> (in units of Hz) to see the effect on your signals.</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/transducer-design#post-344 Thu, 01 Mar 2012 16:36:28 +0000 Hussain 344@http://www.k-wave.org/forum/ <p>Dear Brad,</p> <p>Thanks a lot for your prompt reply. </p> <p>1. My source is based on initial value problem. Right now i'm using makeball function to create the source. I need to model the source replicating a tissue with dimensions 5*5*5cm. I have the value for initial pressure distribution which is the product of fluence and absorption co-efficient. How can i use the value for modelling in k-wave?</p> <p>2. As you pointed out, i tried using the height of the transducer and it's now working.. Thanks a lot!!!</p> <p>3. The center frequency of my detector transducer is say 1MHz. How do i use the gaussianFilter function to design the transducer. I would also like to see, the variations in the signal when the transducer frequency changes? </p> <p>Regards,</p> <p>Hussain </p> http://www.k-wave.org/forum/topic/transducer-design#post-339 Wed, 29 Feb 2012 05:46:54 +0000 Bradley Treeby 339@http://www.k-wave.org/forum/ <p>Hi Hussain,</p> <p>1. Quick question: Are you trying to model a photoacoustic source (initial value problem) or a time varying source? You can model arbitrary photoacoustic sources simply by assigning your source distribution to <code>source.p0</code> (which is a 3D matrix the same size as your grid). For time varying sources, you can similarly assign <code>source.p_mask</code> to be any distribution. You can then supply either one time series which is applied to all source points, or a separate one for each source point.</p> <p>2. Do you have dimensions for the height or width of your transducer in addition to the radius? This will also affect the effective focal length of the transducer you have created. To see what I mean in the case where the transducer is used a source, open up the m-file for the <a href="http://www.k-wave.org/documentation/example_tvsp_transducer_field_patterns.php">Simulating Transducer Field Patterns Example</a> and replace line 50 (the definition of the source mask) with</p> <p><code>height = 5;</code><br /> <code>source.p_mask = makeCircle(Nx, Ny, 61, Ny/2, 60, pi);</code><br /> <code>source.p_mask(height:end, :) = 0;</code></p> <p>Try changing the height value to see what happens.</p> <p>To simulate a particular frequency response, you will need to filter the pressure values returned by k-Wave. For example, you could try using <a href="http://www.k-wave.org/documentation/gaussianFilter.php">gaussianFilter</a>.</p> <p>As there seem to be a number of users interested in curved ultrasound transducers, we will try and include some additional functions and examples in the next k-Wave release.</p> <p>I hope that helps,</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/transducer-design#post-335 Tue, 28 Feb 2012 22:58:01 +0000 Hussain 335@http://www.k-wave.org/forum/ <p>Pls do find the attached image of my detector and object</p> <p><img src="http://img713.imageshack.us/img713/971/object3d.jpg" /> </p> http://www.k-wave.org/forum/topic/transducer-design#post-334 Tue, 28 Feb 2012 22:27:08 +0000 Hussain 334@http://www.k-wave.org/forum/ <p>Hello Bradley and Ben,</p> <p>I am indeed thankful to you for the latest version release. I'm using the focused detector example for my simulations. I have couple of questions </p> <p>1. In most of the examples discussed the source is either disk or sphere; My source depicts the model of tissue or blood vessel so is it possible to use a cubical object as source?</p> <p>2. I need to use a spherical focused transducer. I used program similar to the example "focused detector in 3D", the issue is i need to design a transducer of 1MHz frequency and the diameter of my transducer is lesser than the focal length. So after the simulations are finished my objects looks very close to the detector. I saw the example under the section ultrasound simulations, but the transducer is a linear one. Is there a possibility to design the transducer. Also how can i define the frequency of the transducer? </p> <p>Thanks a lot for your time to answer to the questions</p> <p>Regards,</p> <p>Hussain </p>