A 回答 (2件)
- 最新から表示
- 回答順に表示
No.2
- 回答日時:
>(∂/∂x)^2の行列表現とか・・
下にも書いたように、5×5の場合には、(∂/∂x)^2 つまりfxx、などを一意に定めることはできません。
なんで、(質問者さんが)なんらかの指針で、fxx、などを決めないといけません。
下の回答を、具体的に書き下してみますか。
f(x,y)で画素の値を表すことにすれば、(ただし-2≦x≦2、-2≦y≦2)
f(-2,-2) = f(0,0) - 2*fx - 2*fy + 2*fxx + 4*fxy + 2*fyy
f(-1,-2) = f(0,0) - fx - 2*fy + 1/2*fxx + 2*fxy + 2*fyy
f(0,-2) = f(0,0) - 2*fy + 2*fyy
f(1,-2) = f(0,0) + fx - 2*fy + 1/2*fxx - 2*fxy + 2*fyy
f(2,-2) = f(0,0) + 2*fx - 2*fy + 2*fxx - 4*fxy + 2*fyy
f(-2,-1) = f(0,0) - 2*fx - fy + 2*fxx + 2*fxy + 1/2*fyy
f(-1,-1) = f(0,0) - fx - fy + 1/2*fxx + fxy + 1/2*fyy
f(0,-1) = f(0,0) - fy + 1/2*fyy
f(1,-1) = f(0,0) + fx - fy + 1/2*fxx - fxy + 1/2*fyy
f(2,-1) = f(0,0) + 2*fx - fy + 2*fxx - 2*fxy + 1/2*fyy
…
って感じに25個の式ができます。
(ただし、1個は、f(0,0)=f(0,0)という恒等式になってしまうので、実質的には24個の式)
この24個の式を使って、5つの未知数、fx, fy, fxx, fxy, fxx を決めます。
未知数より式の数のほうが多いわけで、24個の式全てを完全に満たすことはできません。
なんで、ならかの指針で、適当に決める必要があります。
例えば、24本全ての式を同程度満たすべきとするなら、
各式の誤差(=左辺-右辺)の2乗を24本の式について総和をとったものが、最小になるように、fx, fy, fxx, fxy, fxxを決めればよいでしょう。
あるいは、下に書いたように、中心に近い画素の式のほうが、中心から遠い画素よりも重みを重くするべきだと考えたなら、
各式の誤差の2乗に、重みをかけたものの総和をとったものが、最小になるように、fx、fy, fxx, fxy, fxxを決めればよいでしょう。
いずれにせよ、そうすれば、fx, fy, fxx, fxy, fxx が25個の画素の値の関数として表せたことになります。
後は、それを、
(fxx*fy^2 +fyyfx^2 −2fxyfxfy)/(fx^2 + fy^2)
に代入すればよいです。
で、代入したものを、各画素の値について1次近似すれば、フィルタ行列になるでしょう。
No.1
- 回答日時:
5x5だと、一意には定まらないですが。
例えば、テイラー展開の式
f(x0+a,y0+b) = f(x0, y0) + fx(x0,y0)*a + fy(x0,y0)*b + 1/2fxx(x0,y0)*a^2 + fxy(x0,y0)*ab + 1/2fyy(x0,y0)*b^2
で、左辺を -2≦a≦2, -2≦b≦2で動かせば、24本の式ができます。
この24本の式から、fx, fy, fxx, fxy, fyyについて解きます。
ただし、未知数がfx, fy, fxx, fxy, fyyの5つで、式が24個あるので、通常の連立方程式では解けないので、
適当に回帰するなりして、とにかくfx, fy, fxx, fxy, fyyを決めます。(真ん中に近い画素の重みを周辺よりも高くした重み付き回帰をするのがいいと思います)
これで、fx, fy, fxx, fxy, fyyが、25個の画素の値で表せたわけで、
(fxx*fy^2 +fyyfx^2 −2fxyfxfy)/(fx^2 + fy^2)
に代入すればよい。
もし、行列であること(線形であること)が必須であるなら、さらに各画素値について1次近似すればよい。
この回答へのお礼
お礼日時:2017/01/26 22:27
ありがとうございます。
式が複雑で結果が出せておりません。
少し補足説明を加えましたので、見ていただければ幸いです。
よろしくお願いいたします。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- Excel(エクセル) Excel(エクセル)でフィルター抽出後、非表示の行を計算しないで、合計を算出する方法 【内容】 添 4 2023/01/30 17:17
- Excel(エクセル) Indirect関数について、Formulatextで抽出した数式を参照したい。 1 2022/12/15 11:16
- 物理学 なめらかな水平面の床の上に、質量 200 g の物体がある。床の面を xy 面とし、鉛直方向に z 1 2022/07/23 11:28
- Excel(エクセル) 特定文字(数字)で行挿入、挿入された行で合計したい 2 2023/03/13 14:30
- C言語・C++・C# このプログラミング誰か教えてくれませんか 1 2022/06/02 15:27
- 数学 あってますか? 1 2022/05/29 13:45
- 数学 多変数の極値、サドルの問題 f(x,y)=x^3-3x-y^2が極値かサドル、またはなにもないから調 3 2023/07/06 10:23
- Visual Basic(VBA) 別シートから年齢別の件数をカウントしたいの続き 5 2023/01/24 00:16
- Excel(エクセル) アウトラインの小計のやり方 1 2023/03/20 11:51
- その他(Microsoft Office) エクセルマクロ オートフィルターでで選択コピー 2 2022/04/18 11:05
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
オペアンプ、コンパレータ等を...
-
水管ボイラの循環ボイラと貫流...
-
DS漢字検定ソフト26問目でいつ...
-
低圧絶縁抵抗測定について教え...
-
導体の起電力に関する質問です...
-
福井大学 工学部に入学しました...
-
自動制御について質問です。 フ...
-
オペアンプ、コンパレータ等を...
-
オペアンプを使用したアナログ...
-
複素数のベクトル図の書き方を...
-
乾電池3本で動くファンがあり...
-
デジタル信号についてです。 線...
-
電検王さんの問題について
-
電子回路の記号
-
pn接合ダイオード およそ0.6Vか...
-
保線屋さんにとっては、基本中...
-
製図です、ハッチング狭すぎで...
-
デジタル信号処理について質問...
-
RFパワーコンバイナの片方のポ...
-
電気電子材料についてです。 c...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
ベクトル解析
-
関数fxのfは何の英単語の略です...
-
為替レート取得
-
市販の電卓で2⇔10⇔16進数ができ...
-
FXはどこで
-
多変数の極値、サドルの問題 f(...
-
関数電卓の使い方
-
fxである程度、稼いだらfxは辞...
-
中学生でもfxはしても良いので...
-
画像処理フィルターの行列表示
-
fxのぱや手法の詳細を探しています
-
大きく勝てている人
-
ドル円相場が158円を付けたGW前...
-
PayPayポイント運用 やめ時がわ...
-
小川ゆうり主催のLDFサロンにつ...
-
トラッキングトレードを儲けず...
-
FXで他人の資金を運用し利益還...
-
FXをやめたい。
-
FXなどで、月利10パーセント前...
-
エアトレードって何?
おすすめ情報
あ、どうも詳しい解説ありがとうございます。
線形1次におとしていけばよいことは何となくわかりました。
ただもし出来れば、計算結果の行列式を教えていただけると幸いです。
(厚かましくてすいません)
例えばラプラシアン(∂^2/∂x^2 + ∂^2/∂x^2)f(x,y)は3x3行列で表現すると
| 1 1 1 |
| 1 -8 1 |
| 1 1 1 |
とかだそうですその導出はわかったのですが、本件は式がもっと複雑なので、
行列まで落とせていません。
・・そもそもその前に(∂/∂x)^2の行列表現とか・・