当前位置：网站首页>Matlab based medical imaging technology filtering backprojection simulation, including direct backprojection, S-L filtering, R-L filtering, LeWitt filtering

Matlab based medical imaging technology filtering backprojection simulation, including direct backprojection, S-L filtering, R-L filtering, LeWitt filtering

2022-07-27 00:26:00 【I love c programming】

Catalog

1. Algorithm description

Medical imaging technology filter back projection contain R-L wave filtering , R-S wave filtering ,Lewitt wave filtering The reconstructed image is clear . Simulation data of human brain phantom and actual CT Image data filtering back projection reconstruction experiment , Compare direct back projection with R——L and S——L Reconstruction result of filtering method .

2. Partial procedure

    function produce(filter)
        for m=1:180/delta;
            creat_nProj(m);% Zero compensation , For convolution 
            
            c_h1=conv(nProj,h1);% Convolution 
            c_h2=conv(nProj,h2);
            c_h3=conv(nProj,h3);
            c=0.5*(N-1)*(1-cos(m*delta)-sin(m*delta));
            for i=1:N;
                for j=1:N;
                    L(i,j)=c+(i-1)*cos(m*delta)+(j-1)*sin(m*delta); % Insert 
                    n = fix(L(i,j));
                    cL=L(i,j)-n;
                    if strcmp(filter,'ori')
                       if (n>0)&(n<255)
                            rProj1(i,j)=rProj1(i,j)+(1-cL)*nProj(n+N)+cL*nProj(n+1+N);
                        elseif n==255
                            rProj1(i,j) =rProj1(i,j)+nProj(n+N);
                        elseif n==0
                            rProj1(i,j) =rProj1(i,j)+nProj(n+N+1);
                       end
                    elseif strcmp(filter,'sl')
                        if (n>0)&(n<255)
                            rProj2(i,j)=rProj2(i,j)+(1 - cL)*c_h1(n+N+N-1)+cL*c_h1(n+N+N);
                        elseif n==255
                            rProj2(i,j) =rProj2(i,j)+c_h1(n+N+N-1);
                        elseif n==0
                            rProj2(i,j)=rProj2(i,j)+c_h1(n+N+N);
                       end
                    elseif strcmp(filter,'rl')
                       if (n>0)&(n<255)
                            rProj3(i,j)=rProj3(i,j)+(1 - cL)*c_h2(n+N+N-1)+cL*c_h2(n+N+N);
                        elseif n==255
                            rProj3(i,j)=rProj3(i,j)+c_h2(n+N+N-1);
                        elseif n==0
                            rProj3(i,j) =rProj3(i,j)+c_h2(n+N+N);
                       end
                    elseif strcmp(filter,'lew')
                        if (n>0)&(n<255)
                            rProj4(i,j)=rProj4(i,j)+(1 - cL)*c_h3(n+N+N-1)+cL*c_h3(n+N+N);
                        elseif n==255
                            rProj4(i,j) =rProj4(i,j)+c_h3(n+N+N-1);
                        elseif n==0
                            rProj4(i,j) =rProj4(i,j)+c_h3(n+N+N);
                        end
                    end
                end
            end
        end
    end