// mahotas.features.haralick bricks itself when provided a dtype larger than uint8 (version 1.4.3)
pixel_data = skimage.exposure.rescale_intensity(pixel_data,out_range=(0,gray_levels-1)).astype(numpy.uint8)
features = numpy.empty((n_directions, 13, len(unique_labels)))
for index, label in enumerate(unique_labels):
label_data = numpy.zeros_like(pixel_data)
After Change
pixel_data = skimage.exposure.rescale_intensity(pixel_data,out_range=(0,gray_levels-1)).astype(numpy.uint8)
props = skimage.measure.regionprops(labels, pixel_data)
per_label = [self.run_mahotas(prop, scale, n_directions) for prop in props]
features = dask.compute(per_label, scheduler="threads")features = numpy.array(features)[0].transpose(1,2,0)
for direction, direction_features in enumerate(features):
for feature_name, feature in zip(F_HARALICK, direction_features):
statistics += self.record_measurement(