フィットの実装

Question

こんにちは私は画像にpsfを収めようとしています。背景は、より低次の多項式によって近似されるべきである。私は定数を取る場合、それは正常に動作します：

def fitter(image, trueImage,psf,intensity): 
    p0 = [intensity]    
    p0.append(np.amin(trueImage)*4**2) 
    meritFun = lambda p: np.ravel(image-(p[0]*psf+p[1])) 
    p = least_squares(meritFun,p0,method='trf') 

は今、私は私の多項式のために、xとyのを定義する方法の問題があります。

#Does not work! 
def fitter(image, trueImage,psf,intensity): 
    p0 = [intensity]    
    p0.append(np.amin(trueImage)*4**2) 
    p0.append(1) 
    p0.append(1) #some clever initial guess 
    meritFun = lambda p: np.ravel(image-(p[0]*psf+p[1]+p[2]*x+p[3]*y)) 
    p = least_squares(meritFun,p0,method='trf') 

xとyはobviuosly指標私です、j私の画像配列は、どのように私のフィットルーチンを教えてください？

Answer 1

ポール装甲はコメントで述べたように、これを解決する一つの方法は、np.ogridを使用することです：あなたはnp.ogridの[ `のようなもののために

def fitter(image, trueImage,psf,intensity): 
    x = np.ogrid[:image.shape[0]] 
    y = np.ogrid[:image.shape[1]] 
    p0 = [intensity]    
    p0.append(np.amin(trueImage)*4**2) 
    p0.append(0) 
    p0.append(0) 
    meritFun = lambda p: np.ravel(image-(p[0]*psf+p[1]+p[2]*x+p[3]*y)) 
    p = least_squares(meritFun,p0,method='trf')