質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
87.20%
Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

解決済

Python3, curve_fitによる2次元画像データのガウシアンフィッティング

oknd1
oknd1

総合スコア0

Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

1回答

0評価

1クリップ

124閲覧

投稿2017/08/31 10:00

編集2022/01/12 10:58

###前提・実現したいこと
Python3とscipy.optimize.curve_fitを用いて二次元画像データからガウシアンフィッティングを行いたいのですが,

Result from function call is not a proper array of floats\.

のエラーメッセージの対処がわからず困っています。
対処法をお教えくださると幸いです。

![画像データ](a1fc50c2710e33efec0234dbe3084d37.png)
###該当のソースコード

Python

import numpy as np from PIL import Image as pil import matplotlib\.pyplot as plt from scipy\.optimize import curve_fit data = np\.array\(pil\.open\("figs/image\.png"\)\) data = data/np\.max\(data\) x,y=np\.meshgrid\(np\.linspace\(0,data\.shape\[1\],data\.shape\[1\]\),np\.linspace\(0,data\.shape\[0\],data\.shape\[0\]\)\) def twoDgaussian\(X,wx,wy,x0,y0\): x,y=X z=np\.exp\(-\(x-x0\)\*\*2/wx\*\*2-\(y-y0\)\*\*2/wy\*\*2\) return z initial=\(50,50,500,500\) popt,pcov=curve_fit\(twoDgaussian,\(x,y\),data,initial\)

###発生している問題・エラーメッセージ

------------------------------------------------------------------------- ValueError Traceback \(most recent call last\) ValueError: object too deep for desired array ------------------------------------------------------------------------- error Traceback \(most recent call last\) <ipython-input-133-8afc316156ee> in <module>\(\) ----> 1 popt,pcov=curve_fit\(twoDgaussian,\(x,y\),data,initial\) C:\\Users\\\*\*\*\\Anaconda3\\lib\\site-packages\\scipy\\optimize\\minpack\.py in curve_fit\(f, xdata, ydata, p0, sigma, absolute_sigma, check_finite, bounds, method, jac, \*\*kwargs\) 734 # Remove full_output from kwargs, otherwise we're passing it in twice\. 735 return_full = kwargs\.pop\('full_output', False\) --> 736 res = leastsq\(func, p0, Dfun=jac, full_output=1, \*\*kwargs\) 737 popt, pcov, infodict, errmsg, ier = res 738 cost = np\.sum\(infodict\['fvec'\] \*\* 2\) C:\\Users\\\*\*\*\\Anaconda3\\lib\\site-packages\\scipy\\optimize\\minpack\.py in leastsq\(func, x0, args, Dfun, full_output, col_deriv, ftol, xtol, gtol, maxfev, epsfcn, factor, diag\) 385 maxfev = 200\*\(n \+ 1\) 386 retval = _minpack\._lmdif\(func, x0, args, full_output, ftol, xtol, --> 387 gtol, maxfev, epsfcn, factor, diag\) 388 else: 389 if col_deriv: error: Result from function call is not a proper array of floats\.

良い質問の評価を上げる

以下のような質問は評価を上げましょう

  • 質問内容が明確
  • 自分も答えを知りたい
  • 質問者以外のユーザにも役立つ

評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

気になる質問をクリップする

クリップした質問は、後からいつでもマイページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

  • プログラミングに関係のない質問
  • やってほしいことだけを記載した丸投げの質問
  • 問題・課題が含まれていない質問
  • 意図的に内容が抹消された質問
  • 過去に投稿した質問と同じ内容の質問
  • 広告と受け取られるような投稿

評価を下げると、トップページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

まだ回答がついていません

会員登録して回答してみよう

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
87.20%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問

同じタグがついた質問を見る

Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。