1000VARS

悲しいことに生後半年のパソコンが壊れたようです
ただ幸運なことに3年保障に入っていたので修理はだせる・・・はず
「はず」っていうのはサポート先があの悪名高きDell サポートだからなんですけどねー

ってなわけで今回は修理に出すまでの攻防を簡潔に書いていく
まず状況について「PCが何の前触れもなくフリーズするようになった。フリーズは1~120分以内に起こる」
意外と手に負えない壊れ方です。チェックディスク中にもフリーズ
教官にも見てもらったが「ハードが原因の可能性が高い」とのこと

さてDellサポートに電話をかける前に事前にチェックした情報を確認
・電話相手は中国人（日本人ではない）
・とりあえずOSのせいにされる
・その他いちゃもんを付けて修理に出させない
要点はこの３つ。これを念頭においてオペレーションを実行してみる

*以下の会話は文字数の問題で簡略化してあります
まず初めに対応してくれたのは日本人ではない男性
簡潔に状況を説明すると実に簡単な解等をくれた
「windowsのせいなので再インストールしろ」
HAHAHA…(怒)　とりあえず「チェックディスクをかけてが異常はなかった」と返す
「windowsの問題はそれでは引っかからない」との返答
そういわれると中級者にも満たない俺は返す言葉がない
とりあえず問題点を根本的なところに戻して
「windowsの再インストール中にフリーズすると予想される」って返答すると・・・
「windowsの再インストールをしないと修理にだせない」とのご回答
・・・・・・・・・・・・・・・はっ！？いやいやそれはおかしいだろｗｗｗ
金を払って保障に入ってるんだから、そこは引き取れよ！
とまぁ、派遣社員相手に自分の考えを述べても仕方ないので、そこで「わかりました」といって第一ラウンド終了

色々あって再インストール終了。だがこの後も色々文句をつけられる可能性があるので、担当者は毎回違うってのを利用して一気に乗り切ろうと考える

2度目の電話に出たのは日本人ではない女性
とりあえず相手に状況を尋ねられた際に状況と「再インストールなどの前の担当者にいわれた事はすべてしたが解決しなかった。だから修理に出す日程を教えてくれ」と一度にいってみた。YesNoで聞くんじゃなくて、自分に有利な具体的な質問をするのは会話の初歩だよね。どっかの偉い人も相手に選択肢は与えないっていってたし。
まぁ、相手がどう考えたかは知らないけど引取りの日程をすぐに教えてくれた。

パソコンが帰ってくる日は未定だけど私はフランスの忠犬の様にただただ待ち続ける日々なのだ。
感想：外人とのコミュニケーションが意外と難しい。微妙なニュアンスが通じない。そして完全なマニュアル対応もそれに拍車をかけている。あと今回一番ぇーってなったセリフが「windowsの再インストールの指導は有料です」・・・そうかいｗ

卒研もほぼ終了。あとは論文だけ・・・こんな論文見る奴いるのか（笑）

とまぁ、今回は卒研の〆でやったボリュームリンダリングについて話そう
ボリュームレンダリングとは簡単にいってしまえば３次元の物体を表示する方法
といっても１つ１つのvoxel(pixelの３次元Ver)に色を配色していくわけではなく、濃度や密度によってよりリアルに物体を再現する
主に医療用につくられた方法だね。必要なものは0~255の輝度の入ったデータ

今回はこの方法を簡単にやってくれるクラスの紹介をしよう（アルゴリズムが非常にやっかいで個人では面倒なのだ）

http://oshiro.bpe.es.osaka-u.ac.jp/people/staff/imura/products/VolumeRenderer

このクラスはGLEW と GLUT をビルドに必要とする。また実行には GLSL に対応したグラフィック カードが必要
実はボリュームレンダリングといえばVTKが一番有名。でも文章が英語なので読み解くのに時間がかかる
そこで日本語表記なところが一番の売り（？）のこれをお勧めする。

さてHPからrenderer.zipをインストールしたところから説明を始めよう。
まずGLUT を使えるようにするには,http://www.wakayama-u.ac.jp/~tokoi/opengl/libglut.html#2　に書いてあるの様するのがベスト
そしてGLEWも　http://glew.sourceforge.net/　から自分に合うものをダウンロードしよう。
GLEWのフォルダにはincludeとLibがあるがこれを下記に入れる
includeの中身はC:\Program Files (x86)\Microsoft SDKs\Windows\v7.0A\include
LibはC:\Program Files (x86)\Microsoft SDKs\Windows\v7.0A\Lib に移す（環境によっては6.0Aかも）

ではこのクラスの使い方について説明しよう
本記事を見ている方はとりあえずボリュームレンダリングをしたい、という人がほとんどと思う
だから自分が持っているデータ（３次元の輝度のデータ）をそのまま表示する方法を書こう
変更するのはMain.cppのInitVolumeRenderer( void ）
ここでdx,dy,dzにはその３次元空間の最大値
idataは輝度が入ったdx*dy*dzの大きさを持つ１次元の配列である
あらかじめファイルなどから読み込んでおこう
変更した例は
for ( int z = 0; z < dz; z++ ) {
for ( int y = 0; y < dy; y++ ) {
for ( int x = 0; x < dx; x++ ) {
data[ x + dx * y + dx * dy * z ]
= static_cast( idata[ x + dx * y + dx * dy * z ] );
}
}
}
volume_renderer->SetVolume( &data[ 0 ], sx, sy, sz );
これでビルドすれば綺麗なボリュームレンダリングが表示されるはず

ここで気になるのは色だが変更には/shader/density_to_color.fragの内容を変える
といっても場所ごとに色を変えるのにはどうすればいいかなどは研究中である
誰か分かった人がいればコメントでもしてくれると嬉しい

卒業前の最後の記事かもね（´・ω・｀）
来年以降も出来たらちびちび加筆するさ

先日までNIfTIファイルに関する研究を行っていたのだが、その研究がめでたくも完了（終了、打ち切りともいう）したので、ここにまとめてみようと思う。

まずNIfTIファイルとは一体何かを説明する。
簡単にいってしまえば、MRI(医療器具で体をスキャンする機械)から得られるデータフォーマットである。（つまり大抵の人には生涯縁のないフォーマット。ちなみにこのファイルの拡張子は.niiなんだが、これを開けるソフト自体がごく少ないほど珍しい。）
データの内容は、スキャンした体の部分を3次元データとして格納している。
今回の研究内容は、これをプログラムで読み込んで画像データとして抽出すること。
順序としては読み込むのが第一段階、画像に変換するのが第二段階といったところ。

それでは第一段階から説明していこう。
NIfTIファイルには、実にありがたいことに専用のライブラリが存在している。
その名も Niftilib 。NIfTIファイルの取り扱いはこれがないと出来ないといっても過言ではない。
というのは、NIfTIファイルの内部構造は公式ページにいってもほぼ非公開に近く、バイナリからプログラミングするのは現状では無理だからだ。
ライブラリの読み込み方とかは割愛させてもらおう。Linux専用？みたいな匂いがしたのでWindowsでやる場合は注意。
俺の場合はWindowsで実行したがインストールがうまくいかなかったので、ライブラリ内のヘッダー、Cファイルは直接プロジェクトにいれた。
ライブラリのヘッダー内にはNIfTIに関係する様々なことが書いてあるので、見分を広めたい人は是非。

実際のプログラミングに関して。
fslio.h nifti1.h nifti1_io.h znzlib.h　を入れてほしい（対応するソースも）。
かなりの関数があるが、実際に使う関数と変数は実に少ない（３次元データを呼び寄せるだけなら）
要所と抜き出すとたったこれだけ。
nifti_image * nim;
unsigned char * idata;
nim = nifti_image_read(inname, 1); //innameはファイル名
idata=(unsigned char *)nim->data;

idataには３次元データが格納されている。
またパラメータとしては３次元データの大きさはnim->nvoxに
x軸の大きさはnim->dim[1], y軸の大きさはnim->dim[2] ,z軸の大きさはnim->dim[3]にそれぞれ格納されている。
idataから( i , j , k )のデータを取り出す場合は idata[ i + j * nim->dim[1] + k *nim->dim[1] * nim->dim[2] ]という指定を。
idataは0 ~ 255の値がグレースケールとして入っている。

実に簡単な作業だが、ライブラリの説明は非常に分かりにくいのでご注意を。

第二段階はこれをBitMapに入れます。
ある一つの座標軸を固定して、残り２つの座標をそれぞれ入れることで３次元データの断面図を得ることができる。
またRGBの場合、idataから得られた値をそれぞれの各要素に入れる。( 例（N,N,N) )

以上で私の卒研発表を終わります。

感想:たかがこれだけのことに時間をかけたな・・・・