Artefacts in SART and OS_SART reconstruction for DBT Scanner
- Dominant language
- MATLAB
- Stars
- 803
- Forks
- 262
- Avg merge
- 9d 13h
- Merged PRs (30d)
- 8
Description
## Expected Behavior
I want to reconstruct a DBT data using 16 projection data into 48 images
## Actual Behavior
I tried FDK, OS_SART and SART but some king of duplicating structures are present in SART and OS_SART reconstructed images. Surprisingly these structures are considerably less in FDK but the contrast of output image is very poor. I have added all three output images below for your reference.



I used FDK to initialise SART and OS_SART.
Number of iterations : 8
Block size : 8
## Code to reproduce the problem (If applicable)
By examining the SART and OS_SART code, i found that for each iteration, a duplicate of these structures are added to the image. The line is given below
The **input res** and **res after one iteration** is shown in the image

```
if nesterov
% The nesterov update is quite similar to the normal update, it
% just uses this update, plus part of the last one.
ynesterov=res+ bsxfun(@times,1./V(:,:,jj),Atb(W(:,:,index_angles(:,jj)).*(proj(:,:,index_angles(:,jj))-Ax(res,geo,angles_reorder(:,jj),'gpuids',gpuids)),geo,angles_reorder(:,jj),'gpuids',gpuids));
res=(1-gamma)*ynesterov+gamma*ynesterov_prev;
else
res=res+lambda* bsxfun(@times,1./V(:,:,jj),Atb(W(:,:,index_angles(:,jj)).*(proj(:,:,index_angles(:,jj))-Ax(res,geo,angles_reorder(:,jj),'gpuids',gpuids)),geo,angles_reorder(:,jj),'gpuids',gpuids));
end
```
## Specifications
- MATLAB version: 2021b
- OS: windows 10
- CUDA version: 10.0
Contributor guide
Assessment
This issue has not been assessed yet.