How to interpolate empty gridcells in the 3D model?

I have a 3D model of size ss = 139 48 101 which has a lot inactive cell having values NaN and therefore appears empty space shown in the figure below. I do not want empty and want to interpolate the iactive cells by using nearby values. How I can do it - interpolation based on near by data values (instead of taking NaN as a single contant).

The data is attached here:

Answer 1

[Mathieu NOE]

this is you solution :

Inpaint over missing data in 1-D, 2-D, 3-D,... ND arrays - File Exchange - MATLAB Central (mathworks.com)

or this one

Inpainting nan elements in 3-d - File Exchange - MATLAB Central (mathworks.com)