Hello everyone, I have recently been attempting to modify some existing skimage code to display regionprops for a labeled image (e.g. area or eccentricity) I initially tried to translate a vectorized bit of old matlab code I had, but gave up on that and decided to alter the existing label2rgb skimage function I am attempting to change each label value to it's area property value similar to the label2rgb "avg" function. so I have: labels = a labeled image out = np.zeros_like(labels) #a blank array labels2 = np.unique(labels) #a vector of label vals out = np.zeros_like(labels) Props = regionprops(labels, ['Area']) bg_label=0 bg = (labels2 == bg_label) if bg.any(): labels2 = labels2[labels2 != bg_label] out[bg] = 0 for label in labels2: mask = (labels == label).nonzero() color = Props[label].area out[mask] = color but the "out" props image does not correspond to the correct area values? Can anyone help me with this? It also throws the following error: "list index out of range" It would certainly be useful to have a way to view the spatial distribution of label properties in this way - perhaps in a future skimage version?