【学习记录】射线补偿原理与实现：从物理模型到MATLAB代码

【学习记录】射线补偿原理与实现从物理模型到MATLAB代码在经颅超声聚焦中颅骨会引入严重的相位畸变和幅度衰减。射线补偿Ray Compensation是一种快速预补偿方法通过沿直线积分介质参数计算每个通道的额外延时和衰减然后反向补偿。本文从原理、数学公式、离散实现、补偿逻辑到完整实现链条系统梳理射线补偿的全部知识点并对照实际脚本说明其局限。 目录原理一句话物理模型数学公式离散实现梯形积分补偿逻辑完整实现链条对应脚本局限总结一、原理一句话把声波从阵元到焦点的路径近似成一条直线沿这条线积分介质的慢度和衰减得到每通道的额外延时和额外衰减然后反向补偿。二、物理模型2.1 声波走直线真实声波在骨-水界面会折射路径弯曲。射线法忽略折射假设波沿直线从阵元走到焦点。⚠️这是一个近似但计算极快适合快速预检。2.2 沿直线积分对每个阵元画一条到焦点的直线。直线上每一点的介质参数声速、衰减从 k-Wave 的介质网格里插值得到。沿这条线积分慢度积分→ 总传播时间衰减积分→ 总衰减2.3 减去纯水参考得到“额外延时”和“额外衰减”这两个差值才是颅骨引入的畸变需要补偿。三、数学公式3.1 额外延时总传播时间纯水参考时间额外延时3.2 额外衰减幂律衰减模型沿直线积分额外衰减相对幅度四、离散实现梯形积分连续积分在代码里用梯形法离散% 1. 把直线分成 nIntervals 段nIntervalsceil(L/sampleStepMm);% 采样间隔 网格步长的一半tlinspace(0,1,nIntervals1);% 2. 计算每个采样点的坐标pointsMmelementMmt*(focusMm-elementMm);% 3. 插值查介质参数soundSpeedinterpolants.soundSpeed(pointsMm);alphaCoeffinterpolants.alphaCoeff(pointsMm);% 4. 梯形权重weightsones(nSamples,1);weights(1)0.5;weights(end)0.5;% 5. 积分dsMmL/nIntervals;travelTimeSec(dsMm*1e-3)*sum(weights.*(1./soundSpeed));alphaIntegral(dsMm/10)*sum(weights.*alphaCoeff);梯形法把积分区间分成小段每段用梯形面积近似再求和。五、补偿逻辑5.1 延时补偿extraDelayUs是绝对延时——穿骨越厚越大。补偿refmax(extraDelayUs);compDelayref-extraDelayUs;物理直觉穿骨最厚的通道延时最大 → 补偿量 0它已经最晚穿骨最薄的通道延时最小 → 补偿量最大它需要等这样所有通道同时到达焦点。5.2 幅度补偿relativeAmplitudeScale是接收幅度比。补偿gainrelativeAmplitudeScale.^(-gamma);% gamma 控制强度gaingain/mean(gain);% 归一化gain(1-blend)blend*gain;% 与均匀混合gainmax(minGain,min(maxGain,gain));% 限幅gammagamma值含义gamma 0不补幅度gamma 0.5半补ray_compensatedgamma 1全补⚠️限幅防止把噪声放大。六、完整实现链条1. 读探头坐标 颅骨介质 ↓ 2. 对每个有效阵元 a. 画直线到焦点 b. 沿直线采样介质参数 c. 梯形积分得时间和衰减 d. 减去纯水参考 ↓ 3. 输出 focus_channel_attenuation.csv extra_delay_us, relative_amplitude_scale ↓ 4. 补偿 compDelay max(extraDelay) - extraDelay gain relativeAmplitudeScale .^ (-gamma) ↓ 5. 构造发射信号 source.p(k,:) gain(k) * circshift(pulse, compDelay(k)/dt) ↓ 6. 跑 k-Wave记录焦平面算指标七、对应脚本脚本职责simulate_skull_attenuation画直线、积分、输出补偿表skull_focus_ray_phase_compensated_export只补延时幅度均匀skull_focus_ray_compensated_export延时 轻度幅度skull_focus_simulation_export构造发射信号、跑仿真、算指标八、局限假设真实后果直线传播骨-水界面强折射相位残差几十度无多次散射颅骨内反射、模式转换能量预测偏差点源阵元有限尺寸 2.5 mm高估焦点锐度静态介质颅骨不均匀个体差异⚠️折射是最大误差源。这就是为什么射线补偿只能到 20–30%而全波时反能到 50–80%。九、总结射线补偿 沿直线积分介质参数得额外延时和衰减反向补偿。维度内容原理直线近似 数值积分实现griddedInterpolant查参数 梯形法积分补偿延时取max - 当前幅度取scale^(-gamma)局限忽略折射精度有限定位快速预检不是最终方案一句话总结射线补偿是经颅超声聚焦的快速预补偿方法沿直线积分介质慢度和衰减得到每通道的额外延时和衰减反向补偿后让所有通道同时到达焦点。但它忽略了折射这一最大误差源精度有限只适合快速预检。