CERN / CERN/TIGRE

SART algorithm stuck in "set_w" when the input image is large

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

Description

## Behavior

When I ran iterative methods like SART on a CT scan (512×512×463, can be downloaded [here](https://drive.google.com/file/d/1BGGxrHR984JRAF0AYlLAWksPivCU9hKw/view?usp=drive_link)), it was stuck in "set_w" (w=Ax(...)) when initializing IterativeReconAlg class ([link to the code](https://github.com/CERN/TIGRE/blob/4598aaa32a9aa3ae31aef45f5fe2f5daa01854ea/Python/tigre/algorithms/iterative_recon_alg.py#L222)). Do you know the possible reason and how to fix it?

Smaller images (like 256×256×256) work well. Running FDK on the 512×512×463 CT also works well.

## Code to reproduce the problem (If applicable)

```
image_pred = algs.sart(projections, geo, angles, niter=4)
```
where "projections, geo, angles" can be obtained from the CT pickle file I attached above.

## Specifications

- MATLAB/python version: python 3.9
- OS: Linux
- CUDA version: 11.3

(Data reference: Zha et. al, [NAF](https://github.com/Ruyi-Zha/naf_cbct))

Thanks!

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.