# 2D Projectors¶

Projectors determine the (implicit) weight matrix of the geometries. All CPU reconstruction algorithms require a projector.

Each projector is only suitable for specific projection geometries.

## Creation¶

```
id = astra_create_projector(...);
```

Create a projector object. Projectors determine the (implicit) weight matrix of the geometries. All CPU reconstruction algorithms require a projector.

This script is a wrapper around astra_mex_projector(‘create’, ...) with a more convenient interface. See below for specifics.

The allocated object must be freed after use with a call to

```
astra_mex_projector('delete', id)
```

## Parallel beam¶

The projectors in this section can be used with the parallel projection geometry.

**line**

```
proj_id = astra_create_projector('line', proj_geom, vol_geom);
```

The weight of a ray/pixel pair is given by the length of the intersection of the pixel and the ray, considered as a zero-thickness line.

**strip**

```
proj_id = astra_create_projector('strip', proj_geom, vol_geom);
```

The weight of a ray/pixel pair is given by the area of the intersection of the pixel and the ray, considered as a strip with the same width as a detector pixel.

**linear**

```
proj_id = astra_create_projector('linear', proj_geom, vol_geom);
```

A ray is traced through successive columns or rows (depending on which are most orthogonal to the ray). The contribution of this column/row to this ray is then given by linearly interpolating between the two nearest volume pixels of the intersection of the ray and the column/row.

This is also known as the Joseph kernel, or a slice-interpolated kernel.

## Fan beam¶

**line_fanflat**

```
proj_id = astra_create_projector('line_fanflat', proj_geom, vol_geom);
```

The weight of a ray/pixel pair is given by the length of the intersection of the pixel and the ray, considered as a zero-thickness line. This projector can be used with the fanflat and fanflat_vec geometries.

**strip_fanflat**

```
proj_id = astra_create_projector('strip_fanflat', proj_geom, vol_geom);
```

The weight of a ray/pixel pair is given by the area of the intersection of the pixel and the ray. The ray is considered as a 2D cone from the source to the full width of the detector pixel. The projector can only be used with the fanflat geometry.

NB: This mathematical model does not properly take into account the fan beam magnification effect.

## Miscellaneous¶

**sparse_matrix**

```
proj_id = astra_create_projector('sparse_matrix', proj_geom, vol_geom);
```

This projector uses a sparse matrix projection geometry. See the documentation for that geometry for details.

**cuda**

```
proj_id = astra_create_projector('cuda', proj_geom, vol_geom);
```

This projector does not directly specify a weight matrix, but instead is intended to let algorithms use GPU/CUDA code. It can be used with parallel, fanflat and fanflat_vec projection geometries.

NB: This functionality has not yet been implemented everywhere. This projector is currently only usable in astra_create_sino and astra_create_backprojection.

## astra_mex_projector¶

astra_mex_projector is used to manage projector objects.

It is a wrapper around the MEX file astra_mex_projector_c.

astra_mex_projector contains the following commands.

- create
- matrix
- volume_geometry
- projection_geometry
- is_cuda
- delete
- clear
- info

**create**

```
id = astra_mex_projector('create', cfg);
```

Create a projector from a config object. This is called internally by the script [[astra_create_projector]], which is the recommended way to create most projectors.

**matrix**

```
matrix_id = astra_mex_projector('matrix', projector_id);
```

Create an explicit sparse matrix for the weight matrix encoded by this projector.

This is only implemented for 2D CPU projectors.

The returned matrix_id can be further manipulated with [astra_mex_matrix]. In particular, it can be retrieved as a Matlab sparse matrix with

```
S = astra_mex_matrix('get', matrix_id);
```

It has to be freed after use with

```
astra_mex_matrix('delete', matrix_id);
```

NB: Such a matrix can be very large for large geometries.

**volume_geometry**

```
vol_geom = astra_mex_projector('volume_geometry', id);
```

Get the volume geometry attached to the given projector object.

NB: This is not fully implemented yet and the return value may not accurately represent the geometry.

**projection_geometry**

```
proj_geom = astra_mex_projector('projection_geometry', id);
```

Get the projection geometry attached to the given projector object.

NB: This is not fully implemented yet and the return value may not accurately represent the geometry.

**is_cuda**

```
c = astra_mex_projector('is_cuda', id);
```

Return if the the projector is a CUDA projector.

**delete**

```
astra_mex_projector('delete', id)
astra_mex_projector('delete', id1, id2, ...)
```

Free a single or multiple projector(s).

**clear**

```
astra_mex_projector('clear')
```

Free all projectors.

**info**

```
astra_mex_projector('info')
```

Print basic information about all allocated projector objects.