Scripts: volmin.py

# Compute the minimum density map value reading the data one plane at a time
# so it does not need to fit entirely in memory.

def map_min(v):
    d = v.data
    isz, jsz, ksz = d.size
    kmin = []
    for k in range(ksz):
        matrix = d.matrix((0,0,k), (isz,jsz,1))
        kmin.append(matrix.min())
    return min(kmin)

from VolumeViewer import volume_list
for v in volume_list():
    print('Map %s minimum %.5g' % (v.name, map_min(v)))