http://www.k-wave.org/forum/topic/running-photoacoustic-simulation Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 18:31:42 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/running-photoacoustic-simulation#post-6422 Wed, 18 Apr 2018 09:23:59 +0000 Bradley Treeby 6422@http://www.k-wave.org/forum/ <p>Hi Pankaj,</p> <p>I don't understand your question. What is it you're stuck on?</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/running-photoacoustic-simulation#post-6409 Thu, 12 Apr 2018 10:14:02 +0000 pankaj warbal 6409@http://www.k-wave.org/forum/ <p>My leave algorithm is complimenting the image I have created on MS paint and taking it as input image for my reconstruction program though the image reconstruction is Ok. Please tell me how to correct this?<br /> Algorithm goes like this..<br /> clear all<br /> p0_magnitude = 1;<br /> p0 = p0_magnitude * loadImage('C:\Users\IIIT-A\Desktop\bin\k.png');<br /> % Grid<br /> Nx = 300;<br /> Ny = 300;<br /> dx = 500e-6;<br /> dy = 500e-6;<br /> kgrid = kWaveGrid(Nx, dx, Ny, dy);<br /> %PML<br /> PML_size = [20,20];<br /> PML_Alpha = [10,10];<br /> PML_inside = false;<br /> plot_PML = true;<br /> %Source<br /> source.p0 = resize(p0,[Nx, Ny]);%to resize the image<br /> medium.sound_speed = 1500*ones(Nx, Ny); % [m/s]<br /> medium.density = 1040; % [kg/m^3]<br /> %Sensor<br /> sensor_radius = 60e-3; % [m]<br /> num_sensor_points = 100;<br /> sensor.mask = makeCartCircle(sensor_radius, num_sensor_points);<br /> %input arguments<br /> input_args = {'PMLSize',PML_size,'PMLAlpha',PML_Alpha,'PMLInside',false,'PlotPML', true}<br /> %simulation<br /> sensor_data = kspaceFirstOrder2D(kgrid, medium, source, sensor, input_args{:}); </p>