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

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

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

胖射线追踪示意图

用动画显示胖射线追踪的过程。

MATLAB代码

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

clear 
close all
clc

xmin = 0;
xmax = 10;
ymin = 0;
ymax = 10;
nex = 10;
ney = 10;
dx = 1;
dy = 1;

x1 = 0;
y1 = 3;
x2 = 10;
y2 = 7;

clear xx yy;hold on
[xx,yy] = meshgrid(xmin:dx:xmax,ymin:dy:ymax);
plot(xx,yy,'g:',xx',yy','g:');
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*