asfortranarray中的運行時警告

我在中使用SPAMS來工作Dictionary Learning。asfortranarray中的運行時警告

這裏是我的代碼=>

import spams 
import numpy as np 
from PIL import Image 
import time 
import matplotlib.pyplot as plt 

img_file = 'gray-car.jpg' 
try: 
    img = Image.open(img_file) 
except: 
    print("Cannot load image %s : skipping test" %img_file) 
I = np.array(img)/255. 

print('Shape : ',I.shape) 

if I.ndim == 3: 
    A = np.asfortranarray(I.reshape((I.shape[0],I.shape[1] * I.shape[2]))) 
    rgb = True 
else: 
    A = np.asfortranarray(I) 
    rgb = False 

m = 8;n = 8; 
X = spams.im2col_sliding(A,m,n,rgb) 

X = X - np.tile(np.mean(X,0),(X.shape[0],1)) 
X = np.asfortranarray(X/np.tile(np.sqrt((X * X).sum(axis=0)),(X.shape[0],1)),dtype=np.float64) 
param = { 'K' : 100, # learns a dictionary with 100 elements 
      'lambda1' : 0.15, 'numThreads' : 4, 'batchsize' : 400, 
      'iter' : 1000} 

tic = time.time() 
D = spams.trainDL(X,**param) 
tac = time.time() 
t = tac - tic 
print('time of computation for Dictionary Learning: %f' %t) 

##param['approx'] = 0 
# save dictionnary as dict.png 
plt.imshow(D) 
plt.show() 
_objective(X,D,param,'dict')

這裏學習=>

noise.py:27: RuntimeWarning: invalid value encountered in true_divide 
X = np.asfortranarray(X/np.tile(np.sqrt((X * 
X).sum(axis=0)),(X.shape[0],1)),dtype=np.float64) 
/usr/local/lib/python3.4/dist-packages/numpy/core/fromnumeric.py:2699: 
VisibleDeprecationWarning: `rank` is deprecated; use the `ndim` 
attribute or function instead. To find the rank of a matrix see 
`numpy.linalg.matrix_rank`. VisibleDeprecationWarning)

，我這裏使用

輸入圖像=>

過程中，我得到的錯誤

來源

2017-03-08 Sudip Das

我想你在除數中有一個零（或nan/inf）。我認爲這與'np.asfortranarray'沒有任何關係。還請分享一些示例輸入，以便我們可以_actually_驗證代碼。 – MSeifert

@MSeifert我正在上傳我在那裏使用過的圖片 –

告訴我們關於'asfortran'，'X/np.tile（np.sqrt（（X * X）.sum（axis = 0））' 。形狀，dtype，是否運行等。 – hpaulj

讓我們猜想X莫名其妙地代表一個圖像，並有一個「白」的邊界：

In [127]: X=np.zeros((10,10)) 
In [128]: X[3:8,3:8]=1 
.... 
In [130]: np.sqrt(X*X).sum(axis=0) 
Out[130]: array([ 0., 0., 0., 5., 5., 5., 5., 5., 0., 0.])

平鋪師可以很容易地包含0，併產生運行時警告：

In [134]: X/np.tile(np.sqrt((X * X).sum(axis=0)),(X.shape[0],1)) 
/usr/local/bin/ipython3:1: RuntimeWarning: invalid value encountered in true_divide 
    #!/usr/bin/python3 
Out[134]: 
array([[  nan,  nan,  nan, 0.  , 0.  , 
     0.  , 0.  , 0.  ,  nan,  nan], 
     [  nan,  nan,  nan, 0.  , 0.  , 
     0.  , 0.  , 0.  ,  nan,  nan], 
      ...  , 
     0.  , 0.  , 0.  ,  nan,  nan]])

我不知道爲什麼它的使用asfortranarray()，但我沒有得到任何從那裏進一步的錯誤。

VisibleDeprecationWarning in python

是帖子問這個VisibleDeprecationWarning的一個例子。那裏有一個np.rank功能：

In [138]: np.rank(X) 
/usr/local/bin/ipython3:1: VisibleDeprecationWarning: `rank` is deprecated; use the `ndim` attribute or function instead. To find the rank of a matrix see `numpy.linalg.matrix_rank`. 
    #!/usr/bin/python3

一些舊的或模糊的一段代碼使用rank代替ndim。

我試圖怪怪asfortranarray，這是我以前沒見過的。但其目前的代碼只是：

return array(a, dtype, copy=False, order='F', ndmin=1)

不應有的不起眼的角落隱藏着一個rank電話。

我想我們需要了解更多關於X，尤其是它的形狀和dtype。也許它的步伐。 im2col_sliding表示這是一個由as_strided生成的滑動窗口數組。

來源

2017-03-08 22:05:25 hpaulj

asfortranarray中的運行時警告

回答

相關問題