• 締切済み

Octaveで表面粗さRaの計算

Octaveを勉強中です。 今， X=[位置座標，高さデータ] のデータから，表面粗さRaを計算するmファイルを作成したいと思っています。 表面粗さRaを求めるには 　・測定データ→位相補償フィルタ処理→うねり曲線/粗さ曲線分離 　　→粗さ曲線から算術平均粗さ（Ra）を算出 するといった手順になり， 位相補償フィルタ処理はカットオフλcとした場合， 　　S(x)= e －π(x/(αλc))^2／(αλc) 　 で表される重み関数ガウシアンフィルタで畳み込み積分するれば良い， http://www.aichi-inst.jp/other/up_docs/no100_02.pdf というところまで調べたのですが，具体的にOctaveの言語での記述ができず困っています。 畳み込み積分のコマンドとしてconvがあるようなのですが，どのように使えばよいか 教えてください。

みんなの回答

  • tac351115
  • ベストアンサー率66% (109/164)
回答No.1

つぎの資料のｐ１９以降辺りを参考にしてください。 http://members2.jcom.home.ne.jp/keiko0/octavefourier.pdf

PeppekepeE
質問者

お礼

ご返答ありがとうございます。紹介していただいた資料が大変参考になりました。がんばってプログラム作ってみます。

  1. カテゴリ
  2. 学問・教育
  3. 数学・算数

関連するQ&A

  • 畳み込み積分と畳み込み演算

    こんにちは。 今、ディジタル信号処理についてレポートを作成しています。 その中で、畳み込みについて書こうと調べているんですが、 手元にある参考書などをみてもよくわからないんです。 まずひとつ、 Σx(t)y(t-i) この式で表されるものを、畳み込み演算と言うのだと思っていました。 しかし過去の質問を見て、これは実は畳み込み和と言うのではないかなと思ったのですが、畳み込み和でいいのでしょうか？ そして、畳み込み積分と畳み込み和は、一体何が違うのかがよくわかりません。 畳み込み積分と畳み込み和は、具体的に何が違うのか？ どのような時は畳み込み積分で、どのようなとき畳み込み和なのか？ また、畳み込み和→畳み込み積分(orその逆）のように、 式の変形で導き出すことはできるのでしょうか？ 畳み込み和は、行列の形に表せること、 これをシステムとして実現したものがFIRフィルタであることはわかります。 畳み込み積分については、 ∫h(t-τ)x(τ)dτ＝h(t)＊x(t) という式であることぐらいしかわかりません； できるだけ詳しく教えていただけると助かります。 よろしくお願いします。

  • ＪＩＳ　表面粗さ　カットオフ値について

    ＪＩＳの表面粗さのカットオフ値の説明文「位相補償形高域フィルタの利得が５０％になる周波数に対応する波長」とあるんですがよく分かりません。もっとわかりやすい表現はないでしょうか？教えてください。

  • Octaveで画像が読めなくて困ってます。

    はじめまして！ 家で夏休みの課題でMATLABの演習をやってるいるのですが、 MATLABは学校でしか使えません。だから、無料配布のOctaveを使っています。 いま、Octaveで二次元FIRフィルタで画像を処理したいのですが、imread('ファイル名')で画像が読み込めません。エラーのメッセージはこんな感じです。 error: could not read file: C:/PROGRA~1/GNUOCT~1.36/tmp/oct-610bfb74.0 error: evaluating if command near line 178, column 4 error: called from `imread' in file `/usr/share/octave/2.1.36/m/image/imread.m' error: near line 2 of file `/octave_files/octave13_3.m 処理したい画像ファイルはoctave_filesに入れてやっていいるのですが、 違うところに入れないといけないのか。 あるいはたんにプログラムに問題ありなのか。 わかる方がいればよろしくお願いします。

  • 表面粗さの評価法について

    表面粗さの評価法について質問させてください。 現在使用している表面粗さ測定器が大変古いため、サーフコーダではJIS B0601:1982でRaなどが出てきます。 しかし、最新のJIS0601:2001で評価したいと思っております。 どのようなフィルタがあるか教えてください。 測定器から出てきた基の数値データからエクセルなどを使用して、JIS 2001のRz,Rzjis，Raなどを求めることはできるのでしょうか。 当方、研究を始めたばかりの素人ですので、上手な質問ができませんがよろしくお願いします。

  • ガウシアンフィルタのCプログラム

    画像を平滑化する手法にガウシアンフィルタというものがあります。 現在、このガウシアンフィルタをCで作成しようとしていますが、いまいち分かりません。 自分なりに調べてみたところ、平滑化の移動平均フィルタやメディアンフィルタに関しては様々な書籍やサイトがあるのですが、ガウシアンフィルタに関してはあまり見つかりませんでした。 行いたい処理は、簡単に次の通りです。 画像ファイルの読み込み ↓ フィルタ処理 ↓ 処理後の画像をファイル出力 ガウシアンフィルタに関して詳しく書いてある書籍やサイトがあればお教え願います。できれば、ソースプログラムも書いてあると助かります。 当方、参考書籍として、「C言語で学ぶ実践画像処理-井上他」という本を使用しています。

  • 表面粗さ測定データにおけるクルトシス値（Ｒｋｕ）…

    表面粗さ測定データにおけるクルトシス値（Ｒｋｕ）について 表題の件につきまして質問させてください。 　2つの同じプラスチック樹脂成形品の内1枚を強アルカリ溶剤に60分程度浸漬処理し、未処理のもう1枚と状態を比較するべく表面粗さ測定をそれぞれ行いました。 そこで得られた結果は、算術平均粗さRa等のデータにそれほど大きな違いは見られなかったのですが、粗さ曲線のクルトシス(Rku)というデータに違いが見られました。 このRKuについては今まで特に、注目した事が無かったデータでもあり、調べたところ、「表面の凹凸の山の尖り度を表し、この値が大きい程山が尖っている。」という意味までは何となく分かったのですが、そもそものこのクルトシス値の数値の計算式等がまったく理解出来ません。 もし誰か知っていらっしゃる方がいたら分かりやすい解説をしていただけないでしょうか？ 宜しくお願いします。

  • 表面あらさとうねりの測定

    初心者なので、どなたか考え方を教えてください。 粗さとうねり（輪郭精度)を,それぞれ指定されてい ケースがある場合、昔はRmax表記していたので 「表面あらさ＝うねり」と判断出来ていたと思いま すが、最近はRaやRzなど粗さ曲線？を用いていると 思います。 例えば、R3±0.05の形状で、面粗さはRa1.6、輪郭精度0.003 などと指定された場合、それぞれどのような測定器で、評価 の仕方をするものなのでしょうか？ 形状測定器で形状と粗さ、うねり（輪郭精度）を測定出来る ものなのでしょうか？ 粗さのパラメータを指定した場合、その粗さ成分を除したもの をうねり（輪郭精度）をしてデータを抽出出来るものなのでし ょうか？

  • 畳み込み積分と畳み込み和の違いについて

    畳み込みについて勉強し始めたばかりなのですがどうしてもわからないとことがあったため投稿させていただきます。 私は以下のurlから畳み込み和について考えていました。 http://www.ic.is.tohoku.ac.jp/~swk/lecture/yaruodsp/conv.html ここでは、単位インパルスδ(t)から応答関数h(t)を考え入力信号x(t)を離散化し、応答関数をそれぞれの入力の大きさで倍にしてその値を重ね合わせで出力信号y(t)を計算していると私は理解しています。 ここで、畳み込み積分を考えますと畳み込み積分は y(t)=∫h(τ)x(t-τ)dτ となりますので、応答関数h(t)と入力信号x(t)に関してその時間に応じた面積を計算しているんだろうなーという風に考えています。 以下のurlにある入力信号と出力信号が矩形波の図からそのことを考えました。 http://ja.wikipedia.org/wiki/%E7%95%B3%E3%81%BF%E8%BE%BC%E3%81%BF そこで疑問が浮かびました。 上で示す畳み込み和による考えならば wikiに書いてあるような矩形波と矩形波の畳み込みを考えるときに 単純に出力を足し合わせていくと最大値が1ではなくもっと大きな値になるんではないでしょうか？？ おそらく自分が理解できていないところがまだあるんだと思います。 馬鹿な質問かもしれませんが、何日か悩んでもよくわかりません。 よろしくお願いします。

  • 畳み込み積分の交換律の証明

    今、学部レベルのデジタル信号処理の教科書を読み直しているのですが、畳み込み積分の記述で疑問に思ったことがあります。 畳み込み積分では交換律が成り立ちますが、 f(t) * g(t) = g(t) * f(t) ... (1) これをどのように証明したらよいか分かりません。 ちなみにWikipediaではこれについての解説が「積分演算に由来する性質」とあるのみで参考になりませんでした。 (1)の両辺を展開すると \int f(x)・g(t-x) dx = \int g(x)・f(t-x) dx と、なります、多分、これの両辺が等しいことを証明すればいいと思うのですが、そこから先に進めません。(僕のデスクトップでは2バイト文字の積分記号が出てこないので\intと書きました) どうぞよろしくお願いいたします。

  • Octaveのprintコマンドについて

    質問を拝見いただき、ありがとうございます。 掲題の件につきまして、質問があります。 　先に私のoctave環境について記載いたします。 [環境] 　octave　version　3.4.3 　mingw　i686コンパイル 　gnuplotでのplot 　setenv（’GNUTERM’,’windows’） 　CPU　dual2core　メモリ4GB [質問内容] 　plotした画像を出力するために 　printコマンドを使用しております。 　 　print　-djpg　filename.jpg 　その際、データ点数が1000点数などのような小規模での 　print出力では数秒でjpgファイルを出力しますが、 　データ点数が数万のような規模をplotする場合、 　数分かかります。 　print関数で呼び出しているgcwin32c.exeが、CPU使用率50%を 　示し、数分かかるという処理の遅さのためよい解決方法を 　探しております。 　JPGSAVEや、ImageMagickといったOctaveに組み込めそうなのは 　見つけましたが、浅学のため組み込み方もわかりません。 　上記の解決方法について、上記のようなライブラリ？ソフト？を 　　Octave　pacageで組み込める（pkg listののち、pkg load）まで 　できる方法などを、詳細にお教えいただけますと幸いです。 　よろしくお願いします。

注目のQ&A

カテゴリ

専門家に質問してみよう

職業から探して質問する

あなたにピッタリな商品が見つかる！ OKWAVEセレクト

質問する