# Generate Rotational Maximum Intensity Projections

**URL:** https://discourse.vtk.org/t/generate-rotational-maximum-intensity-projections/3650
**Category:** Support
**Created:** [June 25, 2020, 3:22pm UTC](https://discourse.vtk.org/t/generate-rotational-maximum-intensity-projections/3650 "2020-06-25T15:22:00Z")
**Posts on this page:** 2
**Page:** 1

<div class="post-metadata">

### Author: ![divya\_chou](https://discourse.vtk.org/user_avatar/discourse.vtk.org/divya_chou/32/2023_2.png) [@divya\_chou](https://discourse.vtk.org/u/divya_chou)
#### Post date: [June 25, 2020, 3:22pm UTC](https://discourse.vtk.org/t/generate-rotational-maximum-intensity-projections/3650/1 "2020-06-25T15:22:00Z")

</div>

I have been redirected here from ITK discourse. I am very new to using ITK, and am working using SITK on Python.  
I have 3D PET scans. From each 3D matrix, I have generated 2 MIPs at 90 degrees - Sagittal and Coronal. This was easily achieved using the sitk.MaximumProjection function.

How can I generate Maximum Intensity Projection (MIP), at 10 degree rotations?  
So that every PET scan gives me 36 MIPs.

---

<div class="post-metadata">

### Author: ![lassoan](https://discourse.vtk.org/user_avatar/discourse.vtk.org/lassoan/32/50_2.png) [@lassoan](https://discourse.vtk.org/u/lassoan)
#### Post date: [June 25, 2020, 4:11pm UTC](https://discourse.vtk.org/t/generate-rotational-maximum-intensity-projections/3650/2 "2020-06-25T16:11:21Z")

</div>

You can use [VTK’s volume renderer](https://vtk.org/doc/nightly/html/classvtkGPUVolumeRayCastMapper.html) to generate maximum intensity projections. An alternative option is to use vtkImageReslice (slab mode), which is much slower but can create 16-bit or floating-point output.

If you find it easier to use a GUI then you can load your image to 3D Slicer and use Volume rendering module to generate MIP images. You can use Screen Capture module to save many renderings at different rotation angles to a list of files. You can automate all steps using Python scripting.
