Python腳本返回一個數組而不是單個值

我有一個使用numpy的Python腳本，它應該在返回單個值之前拍攝圖像並執行一些計算。當我單獨執行每條線時，它按預期工作。當我將它放在一個.py腳本中並從命令行或Canopy內運行時，它將返回一個數組。Python腳本返回一個數組而不是單個值

我已經修改了代碼略微不要求通常的圖像輸入，結果是一樣的：

import numpy as np 

# Instead of loading an image, generate a test case (w or wo structured noise) 
roi = np.random.poisson(38,(256,256)); 
blob = np.random.poisson(5,(128,128)); 
roi[64:192,64:192] = roi[64:192,64:192]+blob; 
# Load the other variables if necessary (i.e., no DICOM to load) 
[xDim,yDim] = [512,512]; 
roiLength = xDim/2; 
pix = 1.18958; 

# Declare memory for the FFTs 
sizeFFT = xDim; 
NPS2D = np.zeros((sizeFFT,sizeFFT)); # declare memory for fft results 
fftslice = np.zeros((sizeFFT,sizeFFT)); 

# Set the dimension of the ROI and pull the pixel size. This will be 
# used for the scaling factor in the 2D NPS. 
deltaX = pix; 
deltaY = pix; 
scaleFactor = (deltaX/roiLength)*(deltaY/roiLength); 

# Calculate the NPS 
roiMean = np.mean(roi); 
fftslice = np.fft.fft2((roi-roiMean),s=[sizeFFT,sizeFFT]); 
NPS2D = scaleFactor*np.fft.fftshift(np.multiply(fftslice,np.conj(fftslice))); 
NPS2D = NPS2D.real; 

# Subtract the white noise from the NPS to get the structured NPS 
stNPS = NPS2D - roiMean*deltaX*deltaY; 

# Calculate SNI 
SNI=sum(stNPS)/sum(NPS2D); 

# Display the result 
print SNI;

如果我執行每一行是0.107213670449（或類似的，因爲它是再生結果隨機數組）。如果我使用python foo.py從命令行運行腳本，或單擊Canopy中的播放按鈕，結果是一個512長度的數組[4.64940089e-03 ... -4.59789051e-02 -7.15113682e-02]，我已經手動刪除了509個條目。

有什麼想法？我錯過了明顯的東西嗎？

來源

2015-04-06 Nick Bevins

使用內置sum功能比使用numpy.sum或數組的sum方法不同。

對>一維數組，python的sum會給一個非常不同的結果：

In [1]: import numpy as np 

In [2]: x = np.arange(100).reshape(10, 10) 

In [3]: x 
Out[3]: 
array([[ 0, 1, 2, 3, 4, 5, 6, 7, 8, 9], 
     [10, 11, 12, 13, 14, 15, 16, 17, 18, 19], 
     [20, 21, 22, 23, 24, 25, 26, 27, 28, 29], 
     [30, 31, 32, 33, 34, 35, 36, 37, 38, 39], 
     [40, 41, 42, 43, 44, 45, 46, 47, 48, 49], 
     [50, 51, 52, 53, 54, 55, 56, 57, 58, 59], 
     [60, 61, 62, 63, 64, 65, 66, 67, 68, 69], 
     [70, 71, 72, 73, 74, 75, 76, 77, 78, 79], 
     [80, 81, 82, 83, 84, 85, 86, 87, 88, 89], 
     [90, 91, 92, 93, 94, 95, 96, 97, 98, 99]]) 

In [4]: sum(x) 
Out[4]: array([450, 460, 470, 480, 490, 500, 510, 520, 530, 540]) 

In [5]: x.sum() 
Out[5]: 4950 

In [6]: np.sum(x) 
Out[6]: 4950

這是因爲Python的總和基本上總結for循環在對象上。

循環> 1d數組將沿第一個軸返回切片。例如。

In [7]: for item in x: 
    ...:  print item 
    ...: 
[0 1 2 3 4 5 6 7 8 9] 
[10 11 12 13 14 15 16 17 18 19] 
[20 21 22 23 24 25 26 27 28 29] 
[30 31 32 33 34 35 36 37 38 39] 
[40 41 42 43 44 45 46 47 48 49] 
[50 51 52 53 54 55 56 57 58 59] 
[60 61 62 63 64 65 66 67 68 69] 
[70 71 72 73 74 75 76 77 78 79] 
[80 81 82 83 84 85 86 87 88 89] 
[90 91 92 93 94 95 96 97 98 99]

在這種情況下，Python的sum有效地給你列的款項（即row1 + row2 + row3 ...）

來源

2015-04-06 13:56:16

感謝注意到，。我從MATLAB修改了這段代碼，由於使用逐行執行方法在python中工作，我沒有添加np。這兩種執行方法之間'sum'的工作方式有什麼不同？ – 2015-04-06 14:02:27

SNI=sum(stNPS)/sum(NPS2D)

以默認python方式跨列進行求和。所以，你會得到長度爲512

數組而是嘗試從sum numpy的

SNI=stNPS.sum()/NPS2D.sum()

來源

2015-04-06 13:56:08 Zero

Python腳本返回一個數組而不是單個值

回答

相關問題