Mean curvature in binary volume

Hi all,

I need to calculate local mean curvatures of a 3D surface. The surface is currently defined as a binary volume, in that I have a meshgrid of X,Y,Z and at every point if it's on the surface it's a 1, if not it's a 0. In other words, the surface is an isosurface at a constant value of 1.

The different options I've seen are using surfature, but that doesn't appear to work for isosurfaces like I have. It only appears to work if I have a 2D contour like "peaks" where the 3rd dimension is really just a scalar value on a 2D plane. Moreover, I've seen that you can calculate the integral of mean curvature by using this, but I need local mean curvature and not integral. Does anyone have any ideas?

Thanks