http://www.k-wave.org/forum/topic/reconstruction-from-experimental-data-for-a-circular-sensor Support for the k-Wave MATLAB toolbox en-US Tue, 12 May 2026 23:36:43 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/reconstruction-from-experimental-data-for-a-circular-sensor#post-539 Mon, 04 Jun 2012 07:35:25 +0000 Bradley Treeby 539@http://www.k-wave.org/forum/ <p>Hi Srijeeta,</p> <p>I would start with just using the Cartesian sensor mask and your measured data (without using interpCartData to remove the gaps), and see if you can reconstruct an image.</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/reconstruction-from-experimental-data-for-a-circular-sensor#post-535 Fri, 01 Jun 2012 11:09:04 +0000 Srijeeta 535@http://www.k-wave.org/forum/ <p>Hi Brad,</p> <p>I am a little confused but I think I have been using the sensor positions defined as a Cartesian sensor mask to run the time-reversal reconstruction. 'interpCartData' is used to interpolate the time-reversal reconstruction data to remove the gaps. Is that where I should modify the mask?</p> <p>Srijeeta </p> http://www.k-wave.org/forum/topic/reconstruction-from-experimental-data-for-a-circular-sensor#post-520 Tue, 22 May 2012 07:01:22 +0000 Bradley Treeby 520@http://www.k-wave.org/forum/ <p>Hi Srijeeta,</p> <p>I am sure we can get it working. Perhaps try first without using <code>interpCartData</code> and just your actual sensor positions defined as a Cartesian sensor mask to see if you can get an image. </p> <p>Brad. </p> http://www.k-wave.org/forum/topic/reconstruction-from-experimental-data-for-a-circular-sensor#post-519 Mon, 21 May 2012 15:04:43 +0000 Srijeeta 519@http://www.k-wave.org/forum/ <p>Hi Brad,</p> <p>The transducer used to record measurements has a centre frequency of 6.25 MHz with a bandwidth slightly higher than 80%. So I am using a grid that allows frequencies till 12-13 MHz. I have hard-coded the sensor positions and time step size into the code instead of using the functions to generate them, and also specified the sound speed. I am using 'interpCartData' to interpolate the Cartesian coordinates of the sensor positions. Is that something I should be worrying about? </p> <p>I am still hoping I can make the k-Wave toolbox work for my problem because it is difficult to find an alternative reconstruction algorithm for such high frequencies.</p> <p>Kind regards,<br /> Srijeeta </p> http://www.k-wave.org/forum/topic/reconstruction-from-experimental-data-for-a-circular-sensor#post-513 Thu, 17 May 2012 09:01:00 +0000 Bradley Treeby 513@http://www.k-wave.org/forum/ <p>Hi Srijeeta,</p> <p>I am not sure what else to suggest other than double checking that the grid size is sufficient to support the frequencies you are interested in, that the signals are filtered to only contain frequencies supported by the grid, and that the sensor positions, time step size, and sound speed match those used in your experiment. Probably best not to worry about attenuation compensation until you've got it working. Are you using <code>interpCartData</code>?</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/reconstruction-from-experimental-data-for-a-circular-sensor#post-504 Tue, 15 May 2012 14:46:48 +0000 Srijeeta 504@http://www.k-wave.org/forum/ <p>Hi Dr. Treeby,</p> <p>Yes, I do see the waves propagating from the sensor positions in circles but the final image does not show any build up at the inclusion locations. Thank you for suggesting 'filterTimeSeries', the measurements appear much less noisier now. However, the reconstruction still does not show much promise. I am not using attenuation compensation presently. Would it be better to do so?</p> <p>Kind regards,<br /> Srijeeta </p> http://www.k-wave.org/forum/topic/reconstruction-from-experimental-data-for-a-circular-sensor#post-498 Wed, 09 May 2012 03:22:59 +0000 Bradley Treeby 498@http://www.k-wave.org/forum/ <p>Hi Srijeeta, </p> <p>I'm not sure what else to suggest. When you run the reconstruction, do you see the waves propagating from the sensor positions into the domain? (Try setting <code>&#39;PlotScale&#39;</code> to <code>&#39;auto&#39;</code> if you're unsure about the pressure magnitudes.) You could also try filtering out the high frequency information in your signals that is not supported by the grid (for example, using <code>filterTimeSeries</code>). Are you using attenuation compensation?</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/reconstruction-from-experimental-data-for-a-circular-sensor#post-491 Wed, 02 May 2012 15:44:17 +0000 Srijeeta 491@http://www.k-wave.org/forum/ <p>Hi Dr. Treeby,</p> <p>The reconstructed image I get back after time reversal is solely noise. I believe I am going wrong somewhere.</p> <p>I used the experimental time points to create 'kgrid.t_array' instead of calling 'makeTime'. The data was acquired at a sampling frequency of 80 MHz, so I use dt = 1.25e-8 sec with 4096 time points. Based on the Courant stability criterion, dt = 1.5664e-8 sec for an area of 40.1x40.1 mm^2 with a grid size of 512x512 points.</p> <p>Thank you for clarifying the intent of the bandlimited sensors example. I suppose I could just concentrate on modifying 'example_pr_2D_TR_circular_sensor' for experimental reconstruction.</p> <p>I would be grateful if you could provide further insight that would help me obtain a reconstructed image from the noisy data.</p> <p>Srijeeta </p> http://www.k-wave.org/forum/topic/reconstruction-from-experimental-data-for-a-circular-sensor#post-486 Tue, 01 May 2012 23:38:01 +0000 Bradley Treeby 486@http://www.k-wave.org/forum/ <p>Hi Srijeeta,</p> <p>I am not sure what you mean when you say you don't get a reconstruction? Do you get an error or just something that doesn't look correct? Another thing to double check is that your array of time points (<code>kgrid.t_array</code>) is constructed to match your experimental acquisition settings. The pressure is in units of Pascals.</p> <p>You could change the bandlimited sensors example to use a Cartesian sensor mask. Note that this example just shows the effect of using a bandlimited sensor, it doesn't correct for this behaviour.</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/reconstruction-from-experimental-data-for-a-circular-sensor#post-485 Tue, 01 May 2012 20:39:23 +0000 Srijeeta 485@http://www.k-wave.org/forum/ <p>Hi Dr. Treeby,</p> <p>Thank you very much for the suggestions. The sensor locations were indeed overlapping and using a finer grid solved the problem. </p> <p>However, I am still unable to get a reconstruction using the experimental data which otherwise gives a reconstruction with the backprojection algorithm. Could you please tell me the unit for pressure used in the toolbox codes?</p> <p>Also, could I use Cartesian sensor locations in the example for reconstruction from bandlimited sensors?</p> <p>Kind regards,<br /> Srijeeta </p> http://www.k-wave.org/forum/topic/reconstruction-from-experimental-data-for-a-circular-sensor#post-482 Tue, 01 May 2012 01:26:36 +0000 Bradley Treeby 482@http://www.k-wave.org/forum/ <p>Hi Srijeeta,</p> <p>Without seeing your code, it's difficult to say exactly what might be going wrong. You could try double checking that your experimental data and sensor mask have the same number of points, and that they are indexed in the correct way (e.g., via <code>whos</code>) using the example file as reference. Another thing to check is whether your Cartesian sensor points are overlapping in discrete space. You could try plotting the sensor mask as k-Wave sees it using <code>imagesc</code> and <code>cart2grid</code>. If any of your sensor points are overlapping so they are not unique when mapped to the discrete grid, you'll need to use a finer grid spacing.</p> <p>Kind regards,</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/reconstruction-from-experimental-data-for-a-circular-sensor#post-479 Mon, 30 Apr 2012 07:14:03 +0000 Srijeeta 479@http://www.k-wave.org/forum/ <p>Hi Dr. Treeby and Dr. Cox,</p> <p>I am trying to get back reconstructed images from time-domain pressure data acquired using a curved array sensor. In order to do so, I have modified 'example_pr_2D_TR_circular_sensor.m' to use the experimental data and respective sensor locations. Could you please help me to sort out the following error?</p> <p>Error in ==&gt; kspaceFirstOrder_inputChecking at 227<br /> sensor.time_reversal_boundary_data(:, new_col_pos) = order_index;</p> <p>Error in ==&gt; kspaceFirstOrder2D at 427<br /> kspaceFirstOrder_inputChecking;</p> <p>The data in 'sensor.time_reversal_boundary_data' can be accessed as (sensor,time) and I have specified the sensor locations in Cartesian coordinates. Could I be going wrong with the sensor mask then?</p> <p>Thanks,<br /> Srijeeta </p>