You can use vtkImageHistogramStatistics with a stencil to get basic statistics (min, max, mean, stdev). You can compute a lot of features using ITK (SimpleITK or ITKPython), too. If you need even more then you can get the voxel data as a numpy array and use numpy, scipy, etc. packages.