matrixを生成する。できなかったので、先生の助けなども借りながら考えてみた。研究インターンの面接の後、ぐずぐずしながらも、先日の続きであるIsomapの実装を行った。前回は失敗したけど、今回は無事成功。結果どうなるかは知っていたけど、非常にきれいな結果が得られました。Pythonlinalgを使う。laには固有値が、vには固有ベクトル(正確には列ごとに固有ベクトルが並んでいる行列)が代入される。言ってもらってやや頭が混乱してくる。1つ言えるのは自分がよりどころとする部分をはっきりさせないとね、ってことですね。最近はとりあえずPythonで実装してみて、よさげな結果が出そうならCで実装するという流れができつつある。今回の場合はPython依存の関数を使いまくっているので、Cで実装するときどうするか、ということが問題だなぁ。研究Numpyの使い方などを勉強しつつIsomapを実装してみましょう、としてみたけど失敗。固有値、固有ベクトルを求める辺りに問題がありそうななさそうな。Pythonそんなときはidentity()関数を使えばよい。同様のものにzeros()やones()がある。チュートリアルに従って、いろいろ試してみた。mat関数?を使えば行列を定義できる。mat('[1,2,3;4,5,6;7,8,9]')で3行3列の行列が定義できる。昨日の続きで、k近傍グラフとそれらの測地線距離を求めるプログラムを作成した。k近傍グラフが案外めんどくさかった。これから行列計算が連発してくるわけですが、いちいちプログラムを作るのはめんどくさいなぁというわけで、楽になるものを探し中。スイ...昨日の続きで、k近傍グラフとそれらの測地線距離を求めるプログラムを作成した。研究Isomapを実装するために、とりあえず適当なデータが欲しいところ。論文に書いてあったスイスロールがいいののかな。プログラムを考えてみるか。3次元だとx-z軸かy-z軸方向に回転しつつ半径を小さくしながら移動して、移動している位置から適当な...Isomapを実装するために、とりあえず適当なデータが欲しいところ。何か面白い文献はないかなぁといろいろと探していたら面白いものを見つけた。微妙にてこずりながらも完成。こんな感じでいいはず。大きさによるけどデータは500点くらい使えばすかすか感はなくなる。100点くらいだと寂しい感じだった。次はこのデータで近傍グラフを作成します。k近傍グラフでk=3,4,5辺りで...スイスロール作りました2008-10-23論文カキカキ研究OpenCVで遊んでばかりもいられないので、論文をカキカキしている。日本語の投稿論文を書き終え、研究室内レビュー*1もパスしたのでようやく一息、といいたいところだったのに、今度は国際学会の原稿を送らないといけないというワナ。日本語の論文についてはほぼ完成しつつある。後は、明日最終の見直しと最大の関門である先生のチェックをすればおk。論文は細かいところに気を使うのがめんどくさすぎる。投稿論文を書いている。