基于Matlab的跨孔CT胖射线追踪算法（四）

基于Matlab的跨孔CT胖射线追踪算法（四）

CT技术是一种无损的工程物探检测技术，因其方法简单、分辨率高、理论上更可靠、结果更直观，被广泛的应用于各种工程。胖射线追踪是CT技术的一种正演算法，本文的代码和示意图，供学习借鉴。¹

胖射线追踪

MATLAB代码

逐条射线追踪实现胖射线追踪。

clear 
close all
clc

xmin = 0;
xmax = 5;
ymin = 0;
ymax = 5;
nex = 5;
ney = 5;
dx = 1;
dy = 1;
    
x1 = 0;
y1 = 0.85;
x2 = 5;
y2 = 4.2;

clear xx yy;hold on
[xx,yy] = meshgrid(xmin:dx:xmax,ymin:dy:ymax);
plot(xx,yy,'k:',xx',yy','k:');
line([x1,x2],[y1,y2],'color','r');
scatter(x1,y1,'r*');scatter(x2,y2,'rs');
syms x y
a1 = y2-y1;
b1 = x1-x2;
c1 = x2*y1-x1*y2;
f1 = @(x,y) a1*x+b1*y+c1;

lamda = 2; 
c = (1/2)*sqrt((y2-y1)^2+(x2-x1)^2);
b = sqrt(lamda*(sqrt((y2-y1)^2+(x2-x1)^2))/8);
a = sqrt(c^2 + b^2);
theta = atan((y2-y1)/(x2-x1));
x0 = (x1+x2)/2;
y0 = (y1+y2)/2;
f2 = @(x,y) (a^2-c^2*cos(theta)^2)*(x-x0).^2 + ...
     (a^2-c^2*sin(theta)^2)*(y-y0).