I have a solution -   the Computational Morphometry Toolkit can do this job. 
https://www.nitrc.org/projects/cmtk 

(For future reference, I previously tried pynrrd and np.count_nonzero , but that doesn't work)