http://www.k-wave.org/forum/topic/how-to-simulate-the-proton-beam-and-proton-induced-acoustic-emission Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 01:03:16 +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-simulate-the-proton-beam-and-proton-induced-acoustic-emission#post-7774 Wed, 26 Aug 2020 10:27:55 +0000 HustYjq 7774@http://www.k-wave.org/forum/ <p>Hi,zhuoxiaoyu<br /> Have you worked out this question?I have the same problem with you. </p> http://www.k-wave.org/forum/topic/how-to-simulate-the-proton-beam-and-proton-induced-acoustic-emission#post-6555 Thu, 16 Aug 2018 20:39:09 +0000 bencox 6555@http://www.k-wave.org/forum/ <p>Hi Hui,</p> <p>Each row of the matrix <code>source.p</code> should be a time series.<br /> Each row (each time series) is the source applied at one source point in <code>source.p_mask</code>.</p> <p>The source points are ordered columnwise, ie. the first time series (the first row of <code>source.p</code>) is the source applied at the source point with the smallest index according to Matlab's usual matrix numbering convention:<br /> <a href="https://uk.mathworks.com/help/matlab/math/matrix-indexing.html">https://uk.mathworks.com/help/matlab/math/matrix-indexing.html</a>, and so on. </p> <p>Because each row is associated with one grid point, it can be weighted to give the source at that point a different amplitude. The weighting will be determined beforehand, in your case using Geant4 or similar. </p> <p>Hope that's clear!<br /> Ben </p> http://www.k-wave.org/forum/topic/how-to-simulate-the-proton-beam-and-proton-induced-acoustic-emission#post-6530 Thu, 19 Jul 2018 07:19:37 +0000 zhuoxiaoyu 6530@http://www.k-wave.org/forum/ <p>I think I've probably figured out what (num_source_points, num_time_points) means.</p> <p>I modified my code, but I still had some problems.</p> <p>When I use </p> <p>x = 0:1e-7:1e-5;<br /> input_signal = gaussmf(x,[5e-7 8e-6]);<br /> source.p(1:sum(source.p_mask ~= 0), :) = source.p_mask(source.p_mask ~= 0) * input_signal;</p> <p>to define the source, the simulation breaks for "Subscripted assignment dimension mismatch".</p> <p>Could you tell me how to solve this and how to assign "mag_distr(i,j)"(Two dimensional matrix of the same size as the grids) to "source.p"?</p> <p>Thanks a lot! </p> http://www.k-wave.org/forum/topic/how-to-simulate-the-proton-beam-and-proton-induced-acoustic-emission#post-6529 Wed, 18 Jul 2018 08:53:46 +0000 zhuoxiaoyu 6529@http://www.k-wave.org/forum/ <p>Hi Brad.,</p> <p>When I use </p> <p>for i=1:Nx<br /> for j=1:Ny<br /> source.p(i,j)=mag_distr(i,j);<br /> end<br /> end</p> <p>to define the source, the simulation couldn't start for "Input signal must be a vector".<br /> ( mag_distr(i,j) represents pressure distribution about proton beam dose deposition</p> <p>Could you tell me how to define source.p as a matrix ? </p> http://www.k-wave.org/forum/topic/how-to-simulate-the-proton-beam-and-proton-induced-acoustic-emission#post-6278 Sun, 04 Feb 2018 10:36:11 +0000 Bradley Treeby 6278@http://www.k-wave.org/forum/ <p>Hi Hui,</p> <p>The source must always be a two-dimensional matrix with dimensions (num_source_points, num_time_points). The source points are indexed according to MATLAB's standard column-wise linear matrix index ordering.</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/how-to-simulate-the-proton-beam-and-proton-induced-acoustic-emission#post-6273 Mon, 29 Jan 2018 09:50:34 +0000 HuiLIN 6273@http://www.k-wave.org/forum/ <p>Dear Brad.,</p> <p> Could you speak more about how to describe a matrix source.p with different amplitude more detail? Is there appropriate example in k-Wave examples for ref?<br /> My source is an inhomogeneous spatial distribution with a periodic input_signal. When I use<br /> source.p(i,j,k) = source_strength(i,j,k).* input_signal;<br /> to define the source, the code always breaks for different dimension or “Index beyond the matrix dimension”.<br /> Thank you a lot in advance.</p> <p>Hui </p> http://www.k-wave.org/forum/topic/how-to-simulate-the-proton-beam-and-proton-induced-acoustic-emission#post-6245 Mon, 18 Dec 2017 10:22:44 +0000 Bradley Treeby 6245@http://www.k-wave.org/forum/ <p>Hi zhuoxiaoyu,</p> <p>I would suggest using <code>source.p_mask</code> to define the spatial distribution of the source, and then use something like a Gaussian shaped <code>source.p</code> with the appropriate width to model the energy deposition. You will need to define <code>source.p</code> as a matrix, with the amplitude scaled to the deposited dose at each spatial position. See <a href="http://onlinelibrary.wiley.com/doi/10.1118/1.4905047/abstract">this paper</a> for more details.</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/how-to-simulate-the-proton-beam-and-proton-induced-acoustic-emission#post-6244 Sat, 16 Dec 2017 03:41:05 +0000 zhuoxiaoyu 6244@http://www.k-wave.org/forum/ <p>I'm new in K-Wave, now I need to do a simulation about proton beam and proton-induced acoustic emission.</p> <p>In the simulation, The position of the Bragg peak will change with the energy of the proton beam. Now I already have a function about the dose distribution of the proton beam and Bragg, different dose values will produce different levels of acoustic emission. </p> <p>but I don't know how to simulate the acoustic emission images of the proton beam (including the Bragg peak) by using the function "kspaceFirstOrder3D".I've finished defining grids, medium, and sensors. however, I don't know how to use "source.p0"(or source.p_mask and source.p?) to express the proton beam and its acoustic emission. </p> <p>I want to get some help, thank you! </p>