asfortranarrayでの実行時の警告

Question

SPAMSをPythonで使用してDictionary Learningに取り組んでいます。ここで

は私のコードです=>ここで

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) 

私はここで使用している

入力画像=>

中に取得していますエラー

Answer 1

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が含まれており、実行時の警告作り出すことができる：それはasfortranarray()を使っている理由を私は知らないが、私はいずれかを取得しない

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]]) 

をそれ以降のエラー。

VisibleDeprecationWarning in python

このVisibleDeprecationWarningについて尋ねポストの一例です。

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 

コードのいくつかの古いまたはあいまいな部分ではなくndimのrankを使用している：np.rank機能がTheresのです。

私は前に使用したことがないasfortranarrayを責めて誘惑されています。しかし、現在のコードは次のとおりです。

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

rankコールを隠してあいまいなコーナーを持つべきではありません。

私たちはXについてもっと知る必要があると思います。それは形とdtypeです。多分その進歩。 im2col_slidingは、as_stridedで生成されたスライディングウィンドウ配列であることを示しています。