tak0kadaの何でもノート

発声練習、生存確認用。

医学関連は 医学ノート へ

低体温療法のデータいじり(その2)

matplotlib 可視化

低体温療法続き。設定温度からの解離が旧システムと新システムで異なるか検定しようと考えたが、正しい手法を選択できているか怪しい。検定以前に背景の患者因子が揃っていないことに対する注意を払っていないのが大問題だが、そもそもデータが12人分しか手に入らなかった... boxplotで箱ひげ図を描いた後に普通のplot関数を使って生データを○で描きたしている。箱ひげ図にも色を付けている。zorderを指定するとどの要素が上書きされるか指定できる。

追記(2018/7/1): L1 Errorの計算の際に記録体温がNAになってしまっている時間の除外を忘れてしまっていることに気付く

続きを読む

C++コンパイラの名前修飾

C、C++

コンパイラは最終的にはバイナリの実行ファイルを生成するわけだが、少なくともCやC++の場合は分割コンパイルができるので複数のオブジェクトファイルをリンカで結びつけることになっている。リンカはC言語の形式でしか名前を処理できないので、コンパイラはC++の場合でもマングリングと呼ばれるC言語形式への書き換え処理を行ってくれている。

続きを読む

Gnome 3.28でデスクトップのアイコンが表示されなくなったらしい

linux

アップデートしたら見た目が変わっていた。Gnome 3.28 removes the option to display desktop icons. / Applications & Desktop Environments / Arch Linux Forumsによると、とりあえずCarlos Soriano / org.gnome.desktop-icons · GitLabからプラグインをインストールするか、nemoに乗り換えることを推奨されている模様。

Raspberry Pi Zero WHを買った

機械工作

モデルの比較を見て値段が安いので選んだ。正規代理店の方が全体的にamazonより安い感じだった。一緒にケースとブレッドボードとジャンパーワイヤを買って送料込み3,834円、LEDは単品では売っていなかった。Pi3 ModelBは性能が良さそうだけど高い。。GPIOという入出力ポートから外部へ接続するらしく、C++で使えるライブラリはWiringPiやpigpioがあるそうな。買っただけでまだ届いてもいないし届いたらLED点滅させるくらいはしないと。

ロボットアームの運動の行列表示など

機械工作 数値積分

<追記>進捗が出ておらず積みマイコンが1台生えたところで止まっているが、この件の方々は既成品を購入したとのこと。精密な動きが出来るものを作る能力があれば自分のおもちゃ制作に他人のお金を使えたということで自分の無能を大変後悔している。具体的にはハードのプロトタイプと動作イメージ確認のためのプログラムをさくっと作る能力もほしい。</追記>

ロボットアームを使って実験動物の手術をしてみたい人がいる。市販のものは少し大きいのと、自由度が足りないように見え、小さくて細かい動きが出来るものがあると便利だろう。

調べてみると自作している人はいるようで(リンク)、パルス波の周波数を変化させて角度を制御できるサーボモーターというものがあり、それで腕の角度を変化させるらしい。自作するとなるとモーター、アームの素材選びも必要だが、そもそもどう制御するかが分かっていないと手術をできるレベルのものは難しいだろう。以下はロボットアームの運動について調べものをしたメモ。

続きを読む

Cでのプログラミング作法

C、C++

リンク先も含めて安全なコードを書くための情報が多く参考になる。コード例も多い。ただノンプログラマがどこまでやるべきかという問題が...そもそも安全なプログラムが必要ならルールベースではなくてコンパイルの段階で弾かれるようになっていそうなRustを勉強したほうが良さそう。

qiita.com

dicomファイルをまとめて変換するスクリプト

dicom opencv 画像処理 matplotlib

病院で学術用途で匿名化された画像ファイルの出力を依頼すると、

/run/media/tak0kada/dvd/DICOM
|--00000000
|  |--00000000
|  |--00000001
|--00000001
...
|--_DCM_INF
|  |--00000000
|  |  |--_DCM_IDX.txt
|  |  |--overlay.txt
...

というディレクトリ構成のDVDをもらえる。dicomファイルは00000000、00000001と連番になっており、画像/動画データである。スクリプトでまとめて変換したい。pydicomだと難しそうだったがsitkとopencvで動画も出力できた。以下詳細。

続きを読む