f_hsploy.m

来自「MATLAB二维电阻抗断层成像算法!用于医学成像,里面包括有限元剖分正问题,及反」· M 代码 · 共 22 行

function fhsploy(rho,g,H)
%ploy the result by node by node 

%load hsresult2nd
%计算三角形单元的中心
NElement=max(size(H));

for i=1:NElement
    ele=H(i,:);
    Nx(i)=(g(ele(1),1)+g(ele(2),1)+g(ele(3),1))/3;
    Ny(i)=(g(ele(1),2)+g(ele(2),2)+g(ele(3),2))/3;
end
%Nx=Nx';Ny=Ny';
rho01=(rho-min(rho))/(max(rho)-min(rho));%归一化到〔0 1〕
%i=-15:30/250:15;  
i=min(g(:)):(max(g(:))-min(g(:)))/100:max(g(:));
[Xi,Yi]=meshgrid(i,i);  %网格化
Rho=griddata(Nx,Ny,rho01,Xi,Yi,'cubic');%网格化插值
Rho(find(isnan(Rho)))=1;
imagesc(Rho)
%colormap(gray(32))