https://www.k-wave.org/forum/forum/feature-requests Support for the k-Wave MATLAB toolbox en-US Tue, 28 Jul 2026 11:11:34 +0000 http://bbpress.org/?v=1.0.2 q https://www.k-wave.org/forum/search.php https://www.k-wave.org/forum/topic/impose-dirichlet-boundary-conditions-in-k-wave#post-9181 Wed, 12 Feb 2025 21:01:57 +0000 HD Y 9181@https://www.k-wave.org/forum/ <p>Hi, is it possible to impose Dirichlet pressure condition in kspaceFirstOrder3D routine? </p> https://www.k-wave.org/forum/topic/please-add-focus-transducer#post-7419 Sun, 19 Apr 2020 13:38:15 +0000 lishine 7419@https://www.k-wave.org/forum/ <p>i often use points to make a bowl to form a focus transducer,and scan it. But it is so ....so hard! </p> https://www.k-wave.org/forum/topic/kwavearray-distributed_source_signal-and-alternative-source-assignment-methods#post-8327 Wed, 13 Oct 2021 19:20:37 +0000 vat 8327@https://www.k-wave.org/forum/ <p>Hello,</p> <p>I am trying to resolve some discrepancy/ambiguity about input signal assignment using kWaveArray objects to define my transducer/array geometry.</p> <p>Previously I use kWaveTransducer object passing input in m/s to be assigned to particle velocities as follows:</p> <p>% create the input signal using toneBurst<br /> input_signal = toneBurst(1/kgrid.dt, tone_burst_freq, tone_burst_cycles);</p> <p>% scale the source magnitude by the source strength in m/s<br /> input_signal = source_velocity.*input_signal;</p> <p>...</p> <p>transducer.input_signal = input_signal;</p> <p>% create the transducer using the defined settings<br /> transducer = kWaveTransducer(kgrid, transducer);</p> <p>Now a similar geometry defined kWaveArray has a different pressure output and I assume this relates partially to my use but perhaps to implementation? Here is essentially what I tried:</p> <p>karray = kWaveArray;</p> <p>% adding disc elements<br /> for col=1:N<br /> for row=1:N<br /> karray.addDiscElement(elem_pos, diameter, [100,0,0]);<br /> end<br /> end<br /> ...<br /> % assign binary mask<br /> source.p_mask = karray.getArrayBinaryMask(kgrid);<br /> ...<br /> source_sig = repmat(input_signal,count,1);%for assigning to number elements count </p> <p>% get distributed source signals (this automatically returns a weighted<br /> % source signal for each grid point that forms part of the source elements)<br /> source.p = karray.getDistributedSourceSignal(kgrid, source_sig);</p> <p>Based on documentation <a href="http://www.k-wave.org/documentation/example_us_defining_transducer.php" rel="nofollow">http://www.k-wave.org/documentation/example_us_defining_transducer.php</a><br /> I assume the particle velocities/grid points are assigned the input_signal in m/s on my attempt to use kWaveTransducer but here the weighting scheme reduces the values at each of the grid points?</p> <p>Is there a recommended approach that assures grid points for all elements of the kWaveArray are assigned specific particle velocity time traces I would like to pass to them in source_sig or some other structure? </p> <p>Thank you for the attention.</p> <p>Cheers,<br /> Vatche </p> https://www.k-wave.org/forum/topic/element-dependent-sources-in-kwavetransducer-and-kspacefirstorder3d#post-8144 Tue, 04 May 2021 18:46:41 +0000 kwrigby 8144@https://www.k-wave.org/forum/ <p>Although it is uncommon for commercial ultrasound imagers to have element-dependent transmit waveforms (I mean element-dependent waveform shapes, not just time-delays), there are some interesting things one can do with them. Any chance of adding support for this in kspaceFirstOrder3D? I'd try it myself, but I'm just barely Matlab literate.</p> <p>thanks,<br /> Wayne </p> https://www.k-wave.org/forum/topic/power-law-absorption-in-axisymmetric-simulations#post-7912 Tue, 03 Nov 2020 02:31:11 +0000 Karthik Nagabhushana 7912@https://www.k-wave.org/forum/ <p>Hi, </p> <p>I am currently trying to model imaging with a spherically focused transducer. Compared to the conventional 3D representation (kspaceFirstOrder3DC), I have found the axisymmetric representation (kspaceFirstOrderAS) to be extremely useful in terms of compute time and memory requirements for large grid sizes (&gt;1024^3). </p> <p>However, one constraint in the axisymmetric representation is that the full range of acoustic power law is not supported (or medium.alpha_power is set to 2). Is this a fundamental limitation of the axisymmetric approach? If not, do you have plans in the future to support the acoustic power law? </p> <p>Thanks in advance </p> <p>Regards<br /> Karthik </p> https://www.k-wave.org/forum/topic/nonlinear-waves-with-pstdelastic2d-or-pstdelastic3d#post-7438 Sat, 25 Apr 2020 23:42:19 +0000 Pezeril 7438@https://www.k-wave.org/forum/ <p>Hello,<br /> I love k-Wave. That's a great software.<br /> Still, I would like to simulate nonlinear acoustic waves in solids. To do so, I would need a version of pstdElastic where we could input a nonlinear parameter. Perhaps there is a trick to do so. Is it something that could happen in the next k-wave updates? Please let me know.<br /> Best.<br /> Thomas. </p> https://www.k-wave.org/forum/topic/2d-cylindrical-coords#post-132 Thu, 31 Mar 2011 01:32:33 +0000 Richard Stearns 132@https://www.k-wave.org/forum/ <p>I've begun to use your toolbox, and am very impressed with it. Thanks much for providing it as freeware. Most of my current interest is in modeling propagation and scattering from spherically focused acoustic beams. Have you considered allowing one to work in cylindrical coordinates, so that one could model focusing of a spherical beam with the speed of a 2-D simulation? That is, assume no phi dependence in a cylindrical system? I haven't tried it, but I suspect the cartesian 3-D modeling is prohibitively slow (I'm working in the 10MHz regime, with distances on the order of 2.5 cm). If you've done anything with a cylindrical coordinate system, I would be very interested.<br /> Thanks again,<br /> Rick Stearns </p> https://www.k-wave.org/forum/topic/sensor-aperture#post-6660 Sat, 24 Nov 2018 15:52:24 +0000 peiliu 6660@https://www.k-wave.org/forum/ <p>hi!<br /> I have used the k-wave for few weeks and have some problems.I want to do a simulation about the aperture of a single element transducer,so I define a transducer and a source use the 'makeBall' in the 3D.The code is shown blow:<br /> % create the computational grid<br /> Nx = 64; % number of grid points in the x direction<br /> Ny = 64; % number of grid points in the y direction<br /> Nz = 64; % number of grid points in the z direction<br /> dx = 100e-3/Nx; % grid point spacing in the x direction [m]<br /> dy = dx; % grid point spacing in the y direction [m]<br /> dz = dx; % grid point spacing in the z direction [m]<br /> kgrid = kWaveGrid(Nx, dx, Ny, dy, Nz, dz);</p> <p>% define the properties of the propagation medium<br /> medium.sound_speed = 1500; % [m/s]</p> <p>% create the time array<br /> kgrid.makeTime(medium.sound_speed);</p> <p>% create a concave sensor<br /> sphere_offset =10; % [grid points]<br /> diameter = Nx/2 + 1; % [grid points]<br /> bowl_pos = [1 + sphere_offset, Ny/2, Nz/2]; % [grid points]<br /> focus_pos = [Nx/2, Ny/2, Nz/2]; % [grid points]<br /> sensor.mask = makeBowl([Nx, Ny, Nz], bowl_pos, inf, diameter, focus_pos);</p> <p>% create initial pressure distribution using makeBall<br /> ball_magnitude = 5; % [Pa]<br /> ball_radius = 5; % [grid points]<br /> p0 = ball_magnitude * makeBall(Nx, Ny, Nz, 30, Ny/2, Nz/2, ball_radius);</p> <p>and I want to show the transducer and the source togather in the codes as follows:</p> <p>voxelPlot(sensor.mask+p0);<br /> view([130, 40]);</p> <p>my question is why the result is wrong.I execute the codes about the transducer and show it ,it's right,and I execute the codes about the source and show it,it's right as well.But when I try show them togather the transducer is can't be seen.I'm very confused,I need your help,thank you in advance! </p> https://www.k-wave.org/forum/topic/symmetrical-bc#post-6627 Thu, 01 Nov 2018 21:16:57 +0000 dtweten 6627@https://www.k-wave.org/forum/ <p>An option for a Symmetrical Boundary Condition would be a nice addition to k-Wave. This could significantly reduce the grid size for some 3D simulations (by half or one quarter). One of my current models is pushing the maximum grid size. The model is axisymmetric, so this BC option would allow a larger range\field for the simulation.</p> <p>Best Regards,<br /> Dennis </p> https://www.k-wave.org/forum/topic/voxel-type-for-defining-medium-properties-in-3d#post-6541 Tue, 31 Jul 2018 13:21:51 +0000 umitarabul 6541@https://www.k-wave.org/forum/ <p>Hi,</p> <p>I am using a large 3D grid 1024x1024x1024 and medium stores the same size matrix for each property. Instead, you could only store one matrix that stores the type of each voxel with an integer and a small property matrix that maps each voxel type to different property values.</p> <p>I have created an input.h5 with such medium that became over 30GB. I think if you make it as I suggested, it reduces the size significantly.</p> <p>Best, </p> https://www.k-wave.org/forum/topic/new-release-for-matlab-2016b#post-6033 Fri, 30 Jun 2017 21:41:19 +0000 katianne.alcantara 6033@https://www.k-wave.org/forum/ <p>Hi Bradley, </p> <p>First of all, thank you for your amazing work. Do you have any prediction to release a version that runs on Matlab 2016b?</p> <p>Best regards, </p> <p>Katianne. </p> https://www.k-wave.org/forum/topic/sensorrecord_interval#post-5778 Tue, 06 Dec 2016 22:20:52 +0000 nellens1 5778@https://www.k-wave.org/forum/ <p>Is it possible to sample the p or u field at a fixed timestep interval instead of saving everything after sensor.record_start_index, or only the final step (p_final, u_final)? </p> https://www.k-wave.org/forum/topic/heterogeneous-alpha_power#post-5274 Thu, 15 Oct 2015 07:53:24 +0000 Anthony 5274@https://www.k-wave.org/forum/ <p>Hi Brad, </p> <p>Would it be possible to set a heterogeneous <code>alpha_power</code> coefficient?</p> <p>I guess it is impossible (otherwise, you would have allowed it) but, out of curiosity, could you explain the fundamental reason? </p> <p>By the way, do you have a workaround for this when you deal with several different materials? Currently, I just split my simulation domain and run several simulations, one for each material...</p> <p>Best regards,<br /> Anthony </p> https://www.k-wave.org/forum/topic/simulating-ultrasound-speckle-noise#post-4655 Tue, 22 Jul 2014 13:06:23 +0000 f4rhad 4655@https://www.k-wave.org/forum/ <p>This is what I want to do:<br /> <a href="http://8pic.ir/images/uo48la69a8k8cuxlcn27.jpg" rel="nofollow">http://8pic.ir/images/uo48la69a8k8cuxlcn27.jpg</a></p> <p>but I do not want to add gaussian noise this way!<br /> I want to add the normal noise of ultrasound transducer(speckle) with simulation of sonography.<br /> is there any features helping me out?</p> <p>-BR<br /> Farhad </p> https://www.k-wave.org/forum/topic/organizing-code-in-packages#post-4369 Fri, 14 Mar 2014 14:42:22 +0000 DNF 4369@https://www.k-wave.org/forum/ <p>Hi.</p> <p>Great toolbox! But, have you considered putting k-Wave in a Matlab package, that is, organizing it in a '+'-folder? One small annoyance while using k-Wave, and frankly most Matlab code, is namespace clashes.</p> <p>I have a number of functions called 'smooth.m' or 'gaussian.m' and so on, and need to care that one does not shadow another at various times. It would be very convenient if i could write</p> <p><code>&gt;&gt; mat_sm = kwave.smooth(kgrid, mat);</code></p> <p>or </p> <p><code>&gt;&gt; import kwave.*</code><br /> <code>&gt;&gt; mat_sm = smooth(kgrid, mat);</code></p> <p>I believe that using packages is something that Mathworks is trying to encourage, as well. </p> https://www.k-wave.org/forum/topic/thermoelastic-deformation#post-4267 Tue, 14 Jan 2014 13:58:17 +0000 Seongjun Park 4267@https://www.k-wave.org/forum/ <p>Hi, </p> <p>This tool is really amazing!</p> <p>Do you have any plan to include the feature that simulates thermoelastic deformation of the sample when illuminated by a nanosecond laser pulse?<br /> The thermoelastic deformation occurs due to the nonuniform temperature distribution and can be described by thermoelastic wave equation.<br /> Thank you. </p> <p>Best regards,<br /> Seongjun Park </p> https://www.k-wave.org/forum/topic/shear-wave-propagation#post-1721 Mon, 25 Mar 2013 21:08:59 +0000 Petros 1721@https://www.k-wave.org/forum/ <p>I would like to see k-wave toolbox implementing Shear waves (elasto-dynamic equation). I believe you probably have already a working version (e.g., Firouzi et. al., 2012).<br /> Thank you for the great fast and easy to use open-source toolbox!<br /> Best regards,<br /> Petros </p> https://www.k-wave.org/forum/topic/real-directivity#post-937 Tue, 30 Oct 2012 17:18:10 +0000 styxed 937@https://www.k-wave.org/forum/ <p>Hello,</p> <p>Today I wanted to do simulations on single elements sensors with a specific directionality. I was happy to see that directivity was implemented (Sensor Element Directivity in 2D Example). Only to find out that it was just a big cheat. Sorry, maybe not the correct use of words. I was expecting something else.</p> <p>From the documentation of the example I understood that the 'gradient' and 'pressure' were something additional; which would be useful if you want a spatial averaging function. When I plotted the sensor_data of a single sensor element over time, I realized you increase the sensor element to a set of elements, perpendicular to the directivity angle, with a width of sensor.directivity_size to average over all of those elements. So for example the left element already senses the pressure wave sensor.directivity_size pixels earlier that when the wave is supposed to be on the element.</p> <p>What I expected is what I am now posting as feature request:<br /> Directivity in which a sensor output is given by the velocity of the medium in that point (differential of pressure), multiplied (inner product) by the directionality vector of the sensor. (So the inverse of creating a velocity source)</p> <p>Is that possible?</p> <p>Erik </p> https://www.k-wave.org/forum/topic/heterogeneous-absorption-values#post-149 Wed, 18 May 2011 19:05:51 +0000 DanR 149@https://www.k-wave.org/forum/ <p>Hi Brad,<br /> Heterogeneous velocity and density maps are great features that add to the flexibility of K-Wave. Is there any chance you can add absorption to the list? The configuration that we would like to simulate involves significant path lengths through breast tissue and water. Now I have to choose either breast or water (low) absorption for the entire volume, which doesn't represent the system as well as I would like.<br /> Thanks,<br /> Dan </p> https://www.k-wave.org/forum/topic/non-linear-wave-propagation#post-143 Wed, 11 May 2011 14:03:06 +0000 Steffen 143@https://www.k-wave.org/forum/ <p>Hi,<br /> is it possible using k-wave to simulate a high power focused ultrasound transducer? I mean in the focal region the local pressure amplitude will be high let's say 10^6 Pa so non linear effects are expected because the pressure is much higher than the atmospheric or surrounding DC pressure and absolute negative pressure is not expected. </p>