from pylab import * from scipy.ndimage import measurements L = 100 pv = [0.2,0.3,0.4,0.5,0.6,0.7] z = rand(L,L) for i in range(len(pv)): p = pv[i] m = z<p lw, num = measurements.label(m) area = measurements.sum(m, lw, index=arange(lw.max() + 1)) areaImg = area[lw] subplot(2,3,i+1) tit = ‘p=’+str(p) imshow(areaImg, origin=’lower’) title(tit) axis() I’m trying to run this code […]