http://www.k-wave.org/forum/topic/making-an-animation Support for the k-Wave MATLAB toolbox en-US Tue, 12 May 2026 23:07:32 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/making-an-animation#post-6820 Wed, 03 Apr 2019 19:34:43 +0000 ahd123 6820@http://www.k-wave.org/forum/ <p>Just a heads up, the patch() command does not appear to be fully implemented in Octave, specifically transparency. As a result, running voxelPlot() does not work unless line 114<br /> <code>set(hp2, 'FaceAlpha', transparency);</code> is commented out. Also, the edge color does not seem to work. Output is not optimal, but is usable. </p> http://www.k-wave.org/forum/topic/making-an-animation#post-6810 Thu, 28 Mar 2019 21:00:02 +0000 ahd123 6810@http://www.k-wave.org/forum/ <p>Just started using k-wave with Octave. Unfortunately, Octave does not have the videoWriter functionality. It does have a video package which provides 'addframe', 'avifile','aviinfo',and 'aviread' functions. </p> <p>Does anyone know if these functions can be wrapped into a pseudofunction to provide videoWriter-like functionality that k-wave could use without having to rewrite the underlying K-wave movie functionality? </p> http://www.k-wave.org/forum/topic/making-an-animation#post-6030 Fri, 30 Jun 2017 02:41:22 +0000 Alexander 6030@http://www.k-wave.org/forum/ <p>Thank you Brad. Everything seems to be working great now. I look forward to using v1.2 soon. </p> http://www.k-wave.org/forum/topic/making-an-animation#post-6026 Thu, 29 Jun 2017 15:57:07 +0000 Bradley Treeby 6026@http://www.k-wave.org/forum/ <p>Hi eu7209,</p> <p>Yes, you can just set your sensor mask to cover the whole domain. See the "Defining A Sensor Mask By Opposing Corners Example". </p> <p>Brad. </p> http://www.k-wave.org/forum/topic/making-an-animation#post-6025 Thu, 29 Jun 2017 15:47:44 +0000 Alexander 6025@http://www.k-wave.org/forum/ <p>Thanks Brad,</p> <p>I'll re-post here once I test it out in octave. However, is there a way to capture the pressure field as it is? So for FOCUS I could get the pressure field as an array with properties (x,y,z,t) and just loop through t to plot the pressure field (as seen in the above code). Because from my understanding sensor_data is only used for image reconstruction, whereas I just what the pressure field as the simulation runs. </p> http://www.k-wave.org/forum/topic/making-an-animation#post-6021 Thu, 29 Jun 2017 15:01:45 +0000 Bradley Treeby 6021@http://www.k-wave.org/forum/ <p>Hi eu7209,</p> <p>We have actually switched to using <code>VideoWriter</code> for V1.2. The release should be out very soon. Unfortunately, I'm not sure about compatibility with Octave. </p> <p>Brad. </p> http://www.k-wave.org/forum/topic/making-an-animation#post-6014 Mon, 26 Jun 2017 00:30:12 +0000 Alexander 6014@http://www.k-wave.org/forum/ <p>Hi,</p> <p>I am currently using GNU/Octave with k-wave and wanted to make an animation of some of the simulations I am running. The default way to do this in k-wave isn't supported in octave (or newer versions of matlab) by passing ‘RecordMovie’ as true in kspaceFirstOrder2D(...). Therefore I am looking to capture the pressure field at each point in time and then plot it as a mesh plot. </p> <p>I was able to do something similar using FOCUS (<a href="http://www.egr.msu.edu/~fultras-web/index.php" rel="nofollow">http://www.egr.msu.edu/~fultras-web/index.php</a>) with the following code (using matlab)...</p> <pre><code>nt = size(p_tsd, 4); % The amount of time points in the array for it = 1:nt mesh(z*100, x*100, squeeze(p_tsd(:, :, :, it))) title([&#39;FNM TSD Example, t = &#39;, sprintf(&#39;%0.3f&#39;,(it/fs) * 1e6), &#39;\mu&#39;,&#39;s&#39;]) zlabel(&#39;pressure (Pa)&#39;) xlabel(&#39;z (cm)&#39;) ylabel(&#39;x (cm)&#39;) temp=axis(); temp(5)=-maxpressure; temp(6)=maxpressure; axis(temp); cb = colorbar(); caxis([-maxpressure maxpressure]); ylabel(cb, &#39;Pressure (Pa)&#39;); drawnow frame = getframe(1); im = frame2im(frame); [imind,cm] = rgb2ind(im,256); if it == 1 imwrite(imind,cm,&#39;x-z_transient_bone.gif&#39;,&#39;gif&#39;, &#39;Loopcount&#39;,inf); else imwrite(imind,cm,&#39;x-z_transient_bone.gif&#39;,&#39;gif&#39;,&#39;WriteMode&#39;,&#39;append&#39;,&#39;DelayTime&#39;, 0); end end</code></pre> <p>So would something similar be possible in k-wave? </p>