使用 scikit-image 去模糊图像

我正在尝试使用skimage.restoration.wiener http://scikit-image.org/docs/dev/api/skimage.restoration.html#skimage.restoration.wiener，但我总是得到一张带有一堆 1（或 -1）的图像，我做错了什么？原图来自滑铁卢大学 http://links.uwaterloo.ca/Repository.html.

import numpy as np
from scipy.misc import imread
from skimage import color, data, restoration
from scipy.signal import convolve2d as conv2

def main():
  image = imread("/Users/gsamaras/Downloads/boat.tif")
  psf = np.ones((5, 5)) / 25
  image = conv2(image, psf, 'same')
  image += 0.1 * image.std() * np.random.standard_normal(image.shape)

  deconvolved = restoration.wiener(image, psf, 0.00001)
  print deconvolved
  print image

if __name__ == "__main__":
    main()

Output:

[[ 1. -1.  1. ...,  1. -1. -1.]
 [-1. -1.  1. ..., -1.  1.  1.]
 [ 1.  1.  1. ...,  1.  1.  1.]
 ..., 
 [ 1.  1.  1. ...,  1. -1.  1.]
 [ 1.  1.  1. ..., -1.  1. -1.]
 [ 1.  1.  1. ..., -1.  1.  1.]]
[[  62.73526298   77.84202199   94.1563234  ...,   85.12442365
    69.80579057   48.74330501]
 [  74.79638704  101.6248559   143.09978769 ...,  100.07197414
    94.34431216   59.72199141]
 [  96.41589893  132.53865314  161.8286996  ...,  137.17602535
   117.72691238   80.38638741]
 ..., 
 [  82.87641732  122.23168689  146.14129645 ...,  102.01214025
    75.03217549   59.78417916]
 [  74.25240964  100.64285679  127.38475015 ...,   88.04694654
    66.34568789   46.72457454]
 [  42.53382524   79.48377311   88.65000364 ...,   50.84624022
    36.45044106   33.22771889]]

我尝试了几个值。我缺少什么？

到目前为止我最好的解决方案是：

import numpy as np
#import matplotlib.pyplot as plt
from scipy.misc import imfilter, imread
from skimage import color, data, restoration
from scipy.signal import convolve2d as conv2

def main():
  image = imread("/Users/gsamaras/Downloads/boat.tif")
  #plt.imshow(arr, cmap='gray')
  #plt.show()
  #blurred_arr = imfilter(arr, "blur")
  psf = np.ones((5, 5)) / 25
  image = conv2(image, psf, 'same')
  image += 0.1 * image.std() * np.random.standard_normal(image.shape)

  deconvolved = restoration.wiener(image, psf, 1, clip=False)
  #print deconvolved
  plt.imshow(deconvolved, cmap='gray')
  plt.show()
  #print image

if __name__ == "__main__":
    main()

中的值要小得多restoration.wiener()导致图像看起来像是在其上方放置了不透明的覆盖层（例如this http://1.bp.blogspot.com/-ljyhRfysqxo/UhcDG4gyrUI/AAAAAAABexk/cfX6P9tR2-Q/s1600/semedakia1.jpg）。另一方面，随着该值的增加，图像变得越来越模糊。接近 1 的值似乎效果最好，可以使图像去模糊。

值得注意的是，这个值越小（我的意思是balance http://scikit-image.org/docs/dev/api/skimage.restoration.html#skimage.restoration.wiener，图像尺寸越大。

_{PS - I am open to new answers.}