datab = numpy.empty_like(image_data)
numpy.multiply(255.0 / (display_max - display_min), image_data, out=datab)
Using numpy to efficiently convert 16-bit image data to 8 bit for display, with intensity scaling