In [21]:

    

import numpy as np
import matplotlib.pyplot as plt
import scipy.ndimage as img
from skimage import io
%matplotlib inline


    

In [3]:

    




    

In [26]:

    

ampl=500;
a=np.zeros([256,256],dtype=np.uint16);
a[40:210,40:210]=ampl
a=a+100
fa=img.gaussian_filter(a,sigma=1.5)
nfa=fa+np.random.normal(0.0,ampl*0.05,size=a.shape)


    

In [37]:

    

plt.imshow(nfa.astype(np.uint16))
nfa=nfa.astype(np.uint16)
nfa.dtype


    

    
Out[37]:





dtype('uint16')






    

In [38]:

    

io.imsave('sharp.tif',a)
io.imsave('smooth.tif',fa)
io.imsave('noisy_smooth.tif',nfa)


    

    




/Users/kaestner/anaconda/lib/python3.6/site-packages/skimage/io/_io.py:132: UserWarning: sharp.tif is a low contrast image
  warn('%s is a low contrast image' % fname)
/Users/kaestner/anaconda/lib/python3.6/site-packages/skimage/io/_io.py:132: UserWarning: smooth.tif is a low contrast image
  warn('%s is a low contrast image' % fname)
/Users/kaestner/anaconda/lib/python3.6/site-packages/skimage/io/_io.py:132: UserWarning: noisy_smooth.tif is a low contrast image
  warn('%s is a low contrast image' % fname)

In [ ]: