# Optimal Inverse Anscombe Transform
start_time = time.time()
print "Inverse Anscombe"
fimg = InverseAnscombe.OVST(fimg)
elapsed_time = time.time() - start_time
print "Elapsed time:", elapsed_time
outputstack = ImageStack(x, y, z )
for i in xrange(0, z):
# Get the slice at index i and assign array elements corresponding to it.
outputstack.setPixels(fimg[int(i*x*y):int((i+1)*x*y)], i+1)
print 'Preparing denoised image for display '
outputImp = ImagePlus("Output Image", outputstack)
print "OutputImage Stats:"
Stats = StackStatistics(outputImp)
print "mean:", Stats.mean, "minimum:", Stats.min, "maximum:", Stats.max
outputImp.setDisplayRange(Stats.min, Stats.max)
outputImp.show()
fs = FileSaver(outputImp)
