# Resample to Image on point cloud?
**URL:** https://discourse.paraview.org/t/resample-to-image-on-point-cloud/10979
**Category:** ParaView Support
**Created:** [December 14, 2022, 9:51pm UTC](https://discourse.paraview.org/t/resample-to-image-on-point-cloud/10979 "2022-12-14T21:51:33Z")
**Posts on this page:** 9
**Page:** 1
### Author:  [@Devin\_Richard\_Bayly](https://discourse.paraview.org/u/Devin_Richard_Bayly)
#### Post date: [December 14, 2022, 9:51pm UTC](https://discourse.paraview.org/t/resample-to-image-on-point-cloud/10979/1 "2022-12-14T21:51:33Z")
I’m trying to use volumetric rendering with data that comes in as a collection of 24million 3d points. I searched through the documentation and it seems like [resample to image](https://docs.paraview.org/en/latest/UsersGuide/filteringData.html#resample-to-image) is a way to do it if I had a mesh that had attribute data on it, but I only have points.
Is there a method that essentially samples a set of points and turns it into voxel data where the data in a cell is the number of points within that voxel? Is this an option I’m missing in the resample to image filter?
Thanks!
---
### Author:  [@jfavre](https://discourse.paraview.org/u/jfavre)
#### Post date: [December 15, 2022, 8:58am UTC](https://discourse.paraview.org/t/resample-to-image-on-point-cloud/10979/2 "2022-12-15T08:58:36Z")
[https://vtk.org/doc/nightly/html/classvtkPointDensityFilter.html](https://vtk.org/doc/nightly/html/classvtkPointDensityFilter.html) is exactly what you need. The challenge is that it is not directly available in ParaView. Fortunately, you can easily integrate that into a ProgrammableFilter. Here is the code, tested with v5.11.
```auto
PointDensity = ProgrammableFilter(registrationName='PointDensity', Input=mergeBlocks1)
PointDensity.OutputDataSetType = 'vtkImageData'
PointDensity.Script = """from vtk.vtkFiltersPoints import vtkPointDensityFilter
bounds = inputs[0].GetBounds()
hBin = vtkPointDensityFilter()
hBin.SetInputData(inputs[0].VTKObject)
hBin.SetModelBounds(bounds)
hBin.SetSampleDimensions(self.dims)
hBin.Update()
output.ShallowCopy(hBin.GetOutput())
"""
PointDensity.RequestInformationScript = """
executive = self.GetExecutive ()
outInfo = executive.GetOutputInformation(0)
self.dims = [512]*3
outInfo.Set(executive.WHOLE_EXTENT(), 0, self.dims[0]-1 , 0, self.dims[1]-1 , 0, self.dims[2]-1)
"""
PointDensity.RequestUpdateExtentScript = ''
PointDensity.PythonPath = ''
PointDensity.CopyArrays = 0
```
Adding more options to the PointDensityFilter is left as an exercise for the reader.
N.B. This filter has been threaded with vtkSMPTools. It can execute pretty quickly.
---
### Author:  [@Devin\_Richard\_Bayly](https://discourse.paraview.org/u/Devin_Richard_Bayly)
#### Post date: [December 15, 2022, 4:45pm UTC](https://discourse.paraview.org/t/resample-to-image-on-point-cloud/10979/3 "2022-12-15T16:45:42Z")
Incredible! I’ll get right to checking this out, and I’m excited to learn more in the underlying vtk space. Thank you so much.
---
### Author:  [@Devin\_Richard\_Bayly](https://discourse.paraview.org/u/Devin_Richard_Bayly)
#### Post date: [December 15, 2022, 4:50pm UTC](https://discourse.paraview.org/t/resample-to-image-on-point-cloud/10979/4 "2022-12-15T16:50:52Z")
Oh this is totally a one stop shop for everything I wanted
---
### Author:  [@Kenneth\_Moreland](https://discourse.paraview.org/u/Kenneth_Moreland)
#### Post date: [December 15, 2022, 8:08pm UTC](https://discourse.paraview.org/t/resample-to-image-on-point-cloud/10979/5 "2022-12-15T20:08:44Z")
There is also the `Point Volume Interpolator`, which is available in ParaView now.
---
### Author:  [@Devin\_Richard\_Bayly](https://discourse.paraview.org/u/Devin_Richard_Bayly)
#### Post date: [December 15, 2022, 8:11pm UTC](https://discourse.paraview.org/t/resample-to-image-on-point-cloud/10979/6 "2022-12-15T20:11:36Z")
I gave the `Point Volume Interpolator` a try, but since my points didn’t come with attribute data on them it was greyed out?
---
### Author:  [@Joachim\_P](https://discourse.paraview.org/u/Joachim_P)
#### Post date: [December 20, 2022, 2:40am UTC](https://discourse.paraview.org/t/resample-to-image-on-point-cloud/10979/7 "2022-12-20T02:40:28Z")
As an extended exercise, instead of a Python Programmable filter, you can simply describe the filter in XML and load it as a Plugin.
The solution is attached 😉
[PointDensityFilter.xml](https://discourse.paraview.org/uploads/short-url/wxLFhdubpZZuIJW9AsKO1VBa5Vs.xml) (7.7 KB)
---
### Author:  [@jfavre](https://discourse.paraview.org/u/jfavre)
#### Post date: [December 20, 2022, 11:47am UTC](https://discourse.paraview.org/t/resample-to-image-on-point-cloud/10979/8 "2022-12-20T11:47:21Z")
Well done @Joachim_P 👏
---
### Author:  [@Devin\_Richard\_Bayly](https://discourse.paraview.org/u/Devin_Richard_Bayly)
#### Post date: [December 20, 2022, 3:55pm UTC](https://discourse.paraview.org/t/resample-to-image-on-point-cloud/10979/9 "2022-12-20T15:55:22Z")
So fun, I appreciate all the addon growth options from everyone!
Devin