You can get model point positions as numpy array using arrayFromModelPoints and point scalars (if you have added point scalar array) using arrayFromModelPointData, but it would be much slower and more complicated to process the meshes using numpy.
Instead, you can extract non-zero voxels from your volume using vtkThreshold then interpolate values onto the isosurface using vtkPointInterpolator.