저는 공액 그라디언트 방법을 사용하여 MATLAB의 도움으로 이미지의 사이 노 그램을 제공하는 SPECT 이미지를 재구성해야하는 프로젝트에서 작업하고 있습니다. 선형 방정식 Ax = b를 풀어야한다는 것을 깨달았습니다. 여기서 b는 알려진 사이 노 그램 값을 포함하는 N 차원 벡터이고 벡터는 재구성 될 알 수없는 이미지를 포함하는 M 차원 벡터입니다. 시스템 벡터입니다. 내가 가지고있는 사이 노 그램으로 b를 쉽게 만들 수 있지만, 나는 어떻게하면 A를 얻는 지 정말로 모르겠다. 내가 원한 이미지는 128X128이고 사이 노 그램은 128X32 (즉, 32 개의 투영법)이다.공액 그라디언트 2D 이미지 재구성 Matlab
아무도 내 이미지를 재구성하기 위해 matlabs PCG 기능의 도움을 사용할 수있는 방법에 대한 아이디어가 있습니까?
투영 행렬 A의 n 번째 열은 n 번째 픽셀의 라돈 변환이며, 즉 n 번째 요소는 1이고 다른 모든 요소는 0입니다. 이것은 matlab 라돈 함수와 함께 할 수 있지만, 나는 매우 효율적이라고 생각한다. 왜냐하면 A는 희소 할 것이고 루프는 꽤 오랜 시간이 걸릴 것이라고 생각한다. 비록 128x128x128x32가 여전히 관리 가능해야하고 (희박한) 행렬이 .mat 파일이 지오메트리가 변경되지 않는 한. 당신이
N=128; % 128x128 pixels
p=128; % Number of rays per projection
alpha=(0:32)*180/32; % List of projection angles in degrees
A=paralleltomo(N, alpha, p);
나는 희망이 도움을 호출하여 어디서
당신은 예를 http://www2.imm.dtu.dk/~pch/AIRtools/ 를 들어, 다른 MATLAB 스크립트를 시도 할 수 있습니다.