k-Wave

1. cyl
   

   Member
   Joined: Mar '24
   Posts: 4

   Hello, I am currently researching the use of a single element transducer to generate ultrasound penetration through the skull, and then testing the sound field beneath the skull to obtain the attenuation coefficient of the skull. Here is my code. Currently, I have encountered an issue: during simulation, it will display‘transducer=Transducer (kgrid, Transducer)；there is an error Cannot use a value of type kWaveGrid as an index’ , which I have researched a lot of information but still cannot solve. I hope to receive your help. Thank you. Additionally, if there are any issues with my code, I hope you can point them out. Thank you very much again

   Posted 5 months ago #

2. cyl
   

   Member
   Joined: Mar '24
   Posts: 4

   %%
   % 定义网格
   Nx = 300; % number of grid points in the x direction
   Ny = 300; %number of grid points in the y direction
   Nz = 300; % number of grid points in the z direction
   dx =0.1e-3; % grid point spacing in the x direction [m]
   dy =0.2e-3; % grid point spacing in the y direction [m]
   dz =0.1e-3; % grid point spacing in the z direction [m]
   kgrid = kWaveGrid(Nx, dx, Ny, dy, Nz, dz);
   %%
   % define the properties of the propagation medium
   medium.sound_speed = 1540; % [m/s]
   medium.density = 1000; % [kg/m^3]
   medium.sound_speed = medium.sound_speed * ones(Nx, Ny, Nz); % [m/s]
   medium.sound_speed(:,75:80,:) = 3500 ;%填颅骨的宽度和声速
   medium.density = medium.density * ones(Nx, Ny, Nz); % [kg/m^3]
   medium.density(:,75:80,:) = 2000 ;%填颅骨的宽度和密度
   medium.alpha_coeff=0.05*ones(Nx,Ny,Nz); %[dB/(MHz^y cm)]
   medium.alpha_coeff(:,75:80,:) = 2.7; %颅骨衰减系数
   %%
   % tempo e step de simulação
   kgrid.makeTime(medium.sound_speed);
   %%
   % target ping signal
   source_strength = 1e3; % [Pa]
   tone_burst_freq = 1e6; % [Hz]
   tone_burst_cycles = 1;
   input_signal = toneBurst(1/kgrid.dt, tone_burst_freq, tone_burst_cycles);
   %%
   %Setting up ultrasonic transducers
   Transducer.number_elements=1; %阵元个数
   Transducer.element_width = 5e-3; %阵元宽度
   Transducer.element_length = 12 ;%阵元高度
   Transducer.element_spaceing = 0 ;%阵元间距
   Transducer.radius = inf;%换能器曲率半径
   %计算换能器总宽度
   transducer_width = Transducer.number_elements*Transducer.element_width;
   transducer_x_pos = Nx/23;%grid points
   transducer_y_pos = 25;%grid points
   transducer_z_pos = 170;%grid points
   Transducer.positoin = round([transducer_x_pos,transducer_y_pos,transducer_z_pos]);
   % 用于波束合成延迟叠加的参数
   Transducer.sound_speed = 1500; % 声速 [m/s]
   Transducer.focus_distance = 20e-3; % 焦距 [m]
   Transducer.elevation_focus_distance = 19e-3; % elevation plane 仰角平面焦距 [m]
   Transducer.steering_angle = 0; % steering angle 偏转角 [degrees]
   % 设置变迹 apodization
   Transducer.transmit_apodization = 'Rectangular';
   Transducer.receive_apodization = 'Rectangular';
   % 设置当前激活的换能器阵元 define the transducer elements that are currently active
   Transducer.active_elements = 1;
   %设置用于驱动换能器的输入信号
   Transducer.input_signal = input_signal;
   %创建换能器
   transducer = Transducer(kgrid, Transducer);
   %%
   %创建具有检测传感器
   sensor.mask = zeros(Nx,Ny,Nz);
   sensor.mask(Nx/2,275,170)=1;
   %run the simulation
   [sensor_data] = kspaceFirstOrder3D(kgrid , medium , Transducer, sensor,input_args{:});
   %计算输入信号和每个传感器位置记录信号的幅度
   [f_input, as_input] = spect([input_signal,zeros(1,2*length(input_signal))],1/kgrid.dt);
   [~,as_1]=spect(sensor_data(1,:),1/kgrid.dt);
   [sensor_data] = kspaceFirstOrder3D(kgrid, medium, transducer, sensor, input_args{:});

   Posted 5 months ago #

RSS feed for this topic

Reply

You must log in to post.