// 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)
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):