CERN / CERN/TIGRE

'OS_ASD_POCS' reconstruction results in visible streaks

Open
#607 7 comments 0 reactions 0 assignees View on GitHub
Dominant language
MATLAB
Stars
803
Forks
262
Avg merge
9d 13h
Merged PRs (30d)
8

Description

Obviously see a clear streak in the reconfiguration result, I don't know where the parameter is set incorrectly

![image](https://github.com/user-attachments/assets/761611a2-1e41-4f27-9324-ad9fe397479b)
FDK
![image](https://github.com/user-attachments/assets/af04f980-29ed-407c-b910-3c34c2d20d2e)
OS_ASD_POCS
## Code to reproduce the problem (If applicable)
```
%% Demo 4: Simple Image reconstruction
%
%
% This demo will show how a simple image reconstruction can be performed,
% by using OS-SART and FDK
%
%--------------------------------------------------------------------------
%--------------------------------------------------------------------------
% This file is part of the TIGRE Toolbox
%
% Copyright (c) 2015, University of Bath and
% CERN-European Organization for Nuclear Research
% All rights reserved.
%
% License: Open Source under BSD.
% See the full license at
% https://github.com/CERN/TIGRE/blob/master/LICENSE
%
% Contact: tigre.toolbox@gmail.com
% Codes: https://github.com/CERN/TIGRE/
% Coded by: Ander Biguri
%--------------------------------------------------------------------------
%% Initialize

clear;
close all;
%% Geometry
geo=defaultGeometry('nVoxel',[128;128;128]);

%% Load data and generate projections
% define angles
angles=linspace(0,2*pi,100);
% Load thorax phantom data
head=headPhantom(geo.nVoxel);
% generate projections
projections=Ax(head,geo,angles,'interpolated');
% add noise
noise_projections=addCTnoise(projections);

%% Reconstruct image using OS-SART and FDK

% FDK
imgFDK=FDK(noise_projections,geo,angles);
% OS-SART

niter=50;
imgOSSART=OS_ASD_POCS(noise_projections,geo,angles,niter);

% Show the results
plotImg([imgFDK,imgOSSART],'Dim','Z');

```

## Specifications

- MATLAB version:2023b
- OS:
- CUDA version: v11.8

Contributor guide

Open the contributing guide

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.