<?xml version="1.0" encoding="UTF-8" ?><rdf:RDF 
  xmlns="http://purl.org/rss/1.0/"
  xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
  xmlns:atom="http://www.w3.org/2005/Atom"
  xmlns:dc="http://purl.org/dc/elements/1.1/"
  xml:lang="ja">
  <channel rdf:about="http://w.atwiki.jp/pianissimo/">
    <title>pianissimo @ ウィキ</title>
    <link>http://w.atwiki.jp/pianissimo/</link>
    <atom:link href="https://w.atwiki.jp/pianissimo/rss10.xml" rel="self" type="application/rss+xml" />
    <atom:link rel="hub" href="https://pubsubhubbub.appspot.com" />
    <description>pianissimo @ ウィキ</description>

    <dc:language>ja</dc:language>
    <dc:date>2010-08-27T11:54:26+09:00</dc:date>
    <utime>1282877666</utime>

    <items>
      <rdf:Seq>
                <rdf:li rdf:resource="https://w.atwiki.jp/pianissimo/pages/27.html" />
                <rdf:li rdf:resource="https://w.atwiki.jp/pianissimo/pages/1.html" />
                <rdf:li rdf:resource="https://w.atwiki.jp/pianissimo/pages/29.html" />
                <rdf:li rdf:resource="https://w.atwiki.jp/pianissimo/pages/28.html" />
                <rdf:li rdf:resource="https://w.atwiki.jp/pianissimo/pages/16.html" />
                <rdf:li rdf:resource="https://w.atwiki.jp/pianissimo/pages/12.html" />
                <rdf:li rdf:resource="https://w.atwiki.jp/pianissimo/pages/26.html" />
                <rdf:li rdf:resource="https://w.atwiki.jp/pianissimo/pages/13.html" />
                <rdf:li rdf:resource="https://w.atwiki.jp/pianissimo/pages/11.html" />
                <rdf:li rdf:resource="https://w.atwiki.jp/pianissimo/pages/10.html" />
              </rdf:Seq>
    </items>
	
		
    
  </channel>
    <item rdf:about="https://w.atwiki.jp/pianissimo/pages/27.html">
    <title>rmemo</title>
    <link>https://w.atwiki.jp/pianissimo/pages/27.html</link>
    <description>
      R memo

C言語とC#とR
単数の変数を扱うのであれば、ほとんど同じと思われる。違いはC言語とC#は型にうるさいのに対して、Rは型のことは気にしなくても欲しい結果が得られることだろう。
ベクトル、行列、テーブルを扱うときに、話が変わってくる。C言語でやろうとすると、forループが必須である。C#では、forでもできるが、個数を意識しなくてよいforeachや、私は使いこなせないがLINQを使うことによって、あたかもブロック単位に演算しているようなコードが書ける。Rではもっと気楽に、単数とベクトルと行列を同じ書式で記述できてしまう。

たとえばn個の実数(double)からなるベクトルaの成分をすべて2倍したものをbに代入するには、

C言語
 // 配列bはすでに確保されているものとして、、、
 for (i=0; i&lt;n; i++){
  b[i] = a[i] *2;
 }

C#
その1
 var b = new double[a.Length]; 
 for (inti i=0; i&lt;a.Length; i++){
  b[i] = a[i] * 2;
 }
 //これは単に 配列長 n が必要なくなったというだけ。

その2
 var b = from i in a select i * 2;
 //LINQを使った書き方。速度は遅くなるが、今後PLINQが使えるようになれば別かも。

その3
 //ベクトルクラスにして演算子*をオーバーロードすれば、
 var b=a*2;

R
 b = a*2    # b &lt;- a*2のほうが推奨される書き方らしいが・・・

ベクトル以上のデータを扱うときは、Rのほうが短いコードで済んでしまうであろう。プログラミングにかかる時間もRのほうが短いであろう（慣れれば）。しかし問題は実行時間である。Rはインタープリター言語なので、実行速度が遅い。プログラミングにかかる時間と実行時間をよく天秤にかける必要がある。1度実行すれば目的が達せられるプログラムであれば、仮にRではコーディングに1時間、実行に24時間、C#ではコーディングに3時間、実行に10分だとしても、Rでコードする人が多いかもしれない。しかし、これを100回実行する必要がある場合（たとえば患者100人分の画像データを処理するとか）、Rでは計算だけで100日かかるが、C#では数日で終わる。

Rはインストールしたばかりの初期状態でも極めて高機能であるが、驚くべきは付加できるパッケージの膨大さである。生物学関係だけでも圧倒されるほどたくさんある([[ref&gt;http://bioconductor.g.hatena.ne.jp/nakao_mitsuteru/20070618/bioc]])。自分の目的に合ったものを選び出し、正しい方法で使用するのは、相当な苦労がいるに違いない。

画像のフィルターは、convolve関数のtype=&quot;filter&quot;がポイント。
 a=rnorm(100)  #100個の正規乱数
 b=c(1/5,1/5,1/5,1/5,1/5)
 c=convolve(a,b,type=&quot;filter&quot;)
 plot(a,ylim=c(-3,3))
 par(new=T)
 plot(c,ylim=c(-3,3),pch=2)
なお、aの長さは100だが、cの長さは96になっている。

convolve関数の中身は、convolve [enter] で表示されるのでわかりやすいが、
要点は、
1. 被フィルター配列の「前」に、フィルター配列数-1個のゼロを加える
2. フィルター配列の「後」に、ゼロを加えて、上記の「ゼロ＋被フィルター配列」と同じ長さにする。
3. 被フィルター配列をfft、フィルター配列をConj(fft(XX)) して、両者をかけあわせる（普通に要素ごとのかけ算）。
4. その配列をfft(YY,inverse=T)する。
5. できた配列の前と後をトリムして、完成

plot関数の中で、ann=F,axes=F　と書いたのに軸のラベルや数値が表示された！
なぜだ？
調査の結果、Fという名前の変数に何かを代入してあったことが原因と判明した。
横着しないで、ann=FALSE, axes=FALSEと書くことで余計なバグフィックスの時間を節約できたはずだ。
次回から気をつけよう。

DOS窓のコマンドラインから、
R CMD BATCH test.r result.txt
とすると、test.rを実行して、result.txtに出力してくれる。    </description>
    <dc:date>2010-08-27T11:54:26+09:00</dc:date>
    <utime>1282877666</utime>
  </item>
    <item rdf:about="https://w.atwiki.jp/pianissimo/pages/1.html">
    <title>トップページ</title>
    <link>https://w.atwiki.jp/pianissimo/pages/1.html</link>
    <description>
      *このページは、核医学を勉強している大学院生のメモです。
Visitors &amp;counter() (Today &amp;counter(today), Yesterday &amp;counter(yesterday))

**標準脳
[[Talairach座標&gt;http://www.talairach.org/applet/]]は、
xは右向き（右に進むほど増加する）
yは前向き
zは上向きとなっており、
すべてECAT7と逆である。
一方、NEUROSTATは、hdrファイルの中身が
 !the right brain on the left   :=1
 !the anterior to the posterior :=1
 !the superior to the inferior  :=1
このようになっていれば、3軸ともECAT7と同じ向きになる。TalairachからNEUROSTATに変換するときは、3軸ともマイナスを掛ければよい。マイナス1ではなく、
 !slice thickness (mm/pixel) :=2.25
で指定される通常2.25にマイナスをつけてかければよい。

**C#
乱数について：[[メルセンヌツイスター Mersenne Twister C#&gt;http://takel.jp/mt/]]
C#のSystem.Randomは線形合同法というアルゴリズムで乱数を発生させるが、Numerical recipesによると、&quot;Newer use a generator principally based on a linear congruential generator&quot;だというのだ。[[良いアルゴリズム&gt;http://d.hatena.ne.jp/pashango_p/20090717/1247848900]]  [[俺式&gt;http://www.tatsuya-koyama.com/note/program001.html]]
[[c# リフレクションで静的メンバにアクセスする&gt;http://d.hatena.ne.jp/s-kita/20090705/1246806550]]
[[C# FastBitmap&gt;http://blogs.msdn.com/windowsmobile/archive/2008/04/15/faster-c.aspx]]
Dicom c# http://www.fujita-hu.ac.jp/~kmuto/webdas/CyberRad2007-T3-kmuto.pdf
[[プログラムと音楽&gt;http://thorshammer.blog95.fc2.com/blog-entry-253.html#more]]
.NETのツールバーやメニューがnon-active時に1回クリックで反応してくれない件について [[1&gt;http://www.atmarkit.co.jp/bbs/phpBB/viewtopic.php?topic=37772&amp;forum=7]] [[2&gt;http://www.atmarkit.co.jp/bbs/phpBB/viewtopic.php?topic=21432&amp;forum=7]]　Officeも同様の仕様に見えるのだが、決定的に違う点がある。.NETはマウスカーソルをそのメニューやボタンの上に持って行ったとき色が変化し、「いかにもクリックできる」印象を与えるのに対して、Officeではnon-active時にマウスカーソルを重ねても何も変化しない。.NETの設計ミスだと思う。
[[C#で画像をフーリエ変換&gt;http://csharpimage.blog60.fc2.com/blog-entry-14.html]]
[[C#で統計解析&gt;http://d.hatena.ne.jp/KrdLab/searchdiary?word=c%23&amp;.submit=%B8%A1%BA%F7&amp;type=detail]]


**R
[[私のメモ&gt;rmemo]]
[[RjpWiki&gt;http://www.okada.jp.org/RWiki/]]
[[Bioconductorの紹介&gt;http://bioconductor.g.hatena.ne.jp/nakao_mitsuteru/20070618/bioc]]
[[Rによる統計解析（青木繁伸先生）&gt;http://aoki2.si.gunma-u.ac.jp/]] [[dot_plot&gt;http://aoki2.si.gunma-u.ac.jp/R/dot_plot.html]]
[[R 筑波大学ミラー&gt;http://cran.md.tsukuba.ac.jp/]] &quot;Windows&quot;→&quot;base&quot;→最新版DL(2010/1/5現在2.10.1)
推奨の初期設定は、「全てインストール」「SDI」「CHMヘルプ」
[[R設定ファイル（オーム社）&gt;http://www.ohmsha.co.jp/data/link/978-4-274-06783-9/]]のRconsole、Rdevga、Rprofile.siteを、program files\R\以下から探して置換する。
[[エラーバー付きのバーグラフ&gt;http://hosho.ees.hokudai.ac.jp/~kato/unix/R.html#plotCI]]
post hoc [[1&gt;http://www.shiga-med.ac.jp/~koyama/stat/com-ph.html]] [[2&gt;http://www.ibaraki-kodomo.com/toukei/posthoc.html]] [[3:Rを用いるTukey-Kramer&gt;http://133.100.216.71/R_analysis/test_anova03.html]]
[[母回帰直線の係数の推定&gt;http://www.fmu.ac.jp/home/mathema/lec/excel-static/statistics-04.htm]] [[2&gt;http://aoki2.si.gunma-u.ac.jp/lecture/mb-arc/arc023/046.html]]


**統計関連website
[[SWOG&gt;http://www.swogstat.org/statoolsout.html]]　サンプルサイズ、alpha, powerを互いに計算してくれる southwest oncology group
[[StatTools&gt;http://amchang.net/StatTools/index.php]]　StatTools　統計のデザインだけではなく統計解析自体も充実している


**Maxima
[[maximaの使い方&gt;http://cosmo.phys.hirosaki-u.ac.jp/wiki.cgi/maxima]] [[同&gt;http://phys.hirosaki-u.ac.jp/wiki.cgi/maxima?page=FrontPage]]

[[RとMaximaの文法早見表&gt;RMax]]


**Scala
Windowsへのインストール：[[ここでDL&gt;http://www.scala-lang.org/downloads/]]して、zipを展開した後にbinフォルダーにパスを通せばインストール終了。
Macへのインストール：[[ここ&gt;http://d.hatena.ne.jp/tamakiii/20100504/1273000979]]
[[日経BPのScala記事&gt;http://itpro.nikkeibp.co.jp/article/COLUMN/20080613/307981/?ST=develop]]
[[JDK6 update 20&gt;http://java.sun.com/javase/downloads/index.jsp]]　Windows64bitならその対応版をDLすること。
[[Scala for C# programmers&gt;http://hestia.typepad.com/flatlander/2009/01/scala-for-c-programmers-part-1-mixins-and-traits.html]]

**ソフトウェアの使い方
[[Python入門&gt;http://www.pythonweb.jp/tutorial/index.html]]
[[medcon使い方をまとめた&gt;medcon]]


***MISC
[[wavelet解析&gt;http://laputa.cs.shinshu-u.ac.jp/~yizawa/InfSys1/]]
[[数学メモ&gt;math]]
Excelで満年齢を求めるとき、あるいは満月齢でも満日齢でも良いが、datedif()という関数があるので使おう([[ref&gt;http://www.relief.jp/itnote/archives/000423.php]])。この関数は関数一覧に出てこない([[ref&gt;http://www.relief.jp/itnote/archives/002258.php]])。
たとえば、A1セルに 2008/1/17  B1セルに 2010/2/8 が入っているとすると、=DATEDIF(A1,B1,&quot;Y&quot;)とすると満年齢が求められる。&quot;Y&quot;のかわりに、&quot;M&quot;,&quot;D&quot;とすると満月齢、満日齢が求められる。
Impact factorの検索のしかた：　googleで&quot;isi knowledge&quot;→追加情報→Journal of Citation Reports

Vistaでオフライン共有フォルダの検索ができなくて苦労したことがある。解決法は簡単だ。ドライブ名（アルファベット）を割り当てているとき、そこから入っていくとインデックスは作成されていないのだが、ふつうに\\コンピュータ名\\フォルダ名で入っていくと、ちゃんとインデックスは作成されている。
[[Pythagorean means  4種類の平均の幾何学的な意味について&gt;http://en.wikipedia.org/wiki/Pythagorean_means]]
[[WIndows7 TIPs&gt;http://www.microsoft.com/japan/athome/umall/win7/function/01.aspx#02]]
[[JPEG原理&gt;http://webs.lanset.com/crazy17/jp/lab/jpeg1.htm]]
[[クラシックスコア IMSLP&gt;http://imslp.org/wiki/Main_Page]]

[[JNM&gt;http://jnm.snmjournals.org/]] [[Radiology&gt;http://radiology.rsnajnls.org/]] [[NEJM&gt;http://content.nejm.org/]]
[[SNM今後の開催地&gt;http://interactive.snm.org/index.cfm?PageID=221&amp;RPID=221&amp;M=6&amp;Y=2010]]
June 5–9, 2010 Salt Lake City, Utah 済み
June 4–8, 2011 San Antonio, TX 
June 9–13, 2012 Miami, FL 
June 8–12, 2013 Vancouver, BC, Canada 
June 7–11, 2014 St. Louis, MO 
June 6–10, 2015 Baltimore, MD 
June 11–15, 2016 San Diego, California

*[[ウェブ・スクラップ&gt;ウェブ・スクラップ]]










**@wikiへようこそ
-ウィキはみんなで気軽にホームページ編集できるツールです。
-このページは自由に編集することができます。
-メールで送られてきたパスワードを用いてログインすることで、各種変更（サイト名、トップページ、メンバー管理、サイドページ、デザイン、ページ管理、等）することができます

**まずはこちらをご覧ください。
-[[@wikiの基本操作&gt;http://atwiki.jp/guide/category2.html]]
-[[用途別のオススメ機能紹介&gt;http://atwiki.jp/guide/category22.html]]
-[[@wikiの設定/管理&gt;http://atwiki.jp/guide/category6.html]]

**分からないことは？
-[[@wiki ご利用ガイド&gt;http://atwiki.jp/guide/]]
-[[よくある質問&gt;http://atwiki.jp/guide/category1.html]]
-[[無料で会員登録できるSNS内の@wiki助け合いコミュニティ&gt;http://sns.atfb.jp/view_community2.php?no=112]]
-[[@wiki更新情報&gt;http://www1.atwiki.jp/guide/pages/264.html]]
-[[@wikiへお問い合わせ&gt;http://atwiki.jp/guide/contact.html]]
等をご活用ください

**@wiki助け合いコミュニティの掲示板スレッド一覧
#atfb_bbs_list(112)

**その他お勧めサービスについて
-[[フォーラム型の無料掲示板は@bbをご利用ください&gt;&gt;http://atbb.jp/]]
-[[2ch型の無料掲示板は@chsをご利用ください&gt;&gt;http://atchs.jp/]]
-[[お絵かき掲示板は@paintをご利用ください&gt;&gt;http://atpaint.jp/]]
-[[その他の無料掲示板は@bbsをご利用ください&gt;&gt;http://atbbs.jp/]]
-[[無料ブログ作成は@WORDをご利用ください&gt;&gt;http://atword.jp/]]
-[[大容量１Ｇ、PHP/CGI、MySQL、FTPが使える無料ホームページは@PAGES&gt;&gt;http://atpages.jp/]]

**おすすめ機能
-[[気になるニュースをチェック&gt;http://atwiki.jp/guide/17_174_ja.html]]
-[[関連するブログ一覧を表示&gt;http://atwiki.jp/guide/17_161_ja.html]]

**その他にもいろいろな機能満載！！
-[[@wikiプラグイン&gt;http://atwiki.jp/guide/category17.html]]
-[[@wiki便利ツール&gt;http://atwiki.jp/guide/category32.html]]
-[[@wiki構文&gt;http://atwiki.jp/guide/category16.html]]
-[[@wikiプラグイン一覧&gt;http://www1.atwiki.jp/guide/pages/264.html]]

**バグ・不具合を見つけたら？ 要望がある場合は？
お手数ですが、メールでお問い合わせください。    </description>
    <dc:date>2010-08-13T08:59:25+09:00</dc:date>
    <utime>1281657565</utime>
  </item>
    <item rdf:about="https://w.atwiki.jp/pianissimo/pages/29.html">
    <title>RMax</title>
    <link>https://w.atwiki.jp/pianissimo/pages/29.html</link>
    <description>
      このページでは、RとMaximaの文法を比較します。
私はどちらのソフトもそれほど頻繁には使用しないのですが、たまに使用するときに文法を混乱して、毎回成書や解説websiteを当たっています。時間の無駄なので、頻繁に使用する文法だけでも早見表にまとめておこうという趣旨です。
参考文献
[[Maxima 入門ノート1.2.1 中川義行氏&gt;http://www.eonet.ne.jp/~kyo-ju/maxima.pdf]]
[[Maxima による数式処理 弘前大学 葛西 真寿氏&gt;http://phys.hirosaki-u.ac.jp/wiki.cgi/maxima?page=FrontPage]]

*代入と基本演算
||R|Maxima|
|べき乗|^|^|
|100÷7の剰余|100 %% 7|mod(100,7)|
|100÷7の整数商|100 %/% 7||
|日常使用関数|sqrt exp log sin cos tan abs|sqrt exp log sin cos tan abs|
|変数に代入|a=5|a:5|
|ベクトルを代入|a=c(1,2,3)|a:[1,2,3]|
|ベクトル&amp;bold(){a},&amp;bold(){b}の要素ごとの積|a*b|a*b|
|ベクトル&amp;bold(){a},&amp;bold(){b}の内積|a%*%b　または　sum(a*b)|a.b|
|ベクトル&amp;bold(){a},&amp;bold(){b}の外積|outer(a,b)|自前で定義が必要？|
|行列Aに代入&amp;ref(matrix.PNG)|A=matrix(c(1,2,3,4),2,byrow=T)|A:matrix([1,2],[3,4])|
|円周率を浮動小数点表示||float(%pi)|
|自然対数の底e||%e|
|虚数単位||%i|

*関数
||R|Maxima|
|関数を定義|f=function(x) x^2+1|f(x):=x^2+1|
|x^2をxで微分||diff(x^2,x)|
|x^2をxで2階微分||diff(x^2,x,2)|
|x^2をxで不定積分||integrate(x^2,x)|
|数値積分||romberg(x^2,x,0,1)|
|cos(x)を0～pi/2で積分||integrate(cos(x),x,0,%pi/2)|
|テイラー展開 exp(x)を3次まで||taylor(exp(x),x,0,3|
|式の簡単化||ratsimp(f(x))|
|因数分解||factor(f(x))|
|誤差関数&amp;ref(http://upload.wikimedia.org/math/4/d/4/4d4a6284c0034cad3b701f8799205132.png)||erf(x)|

*繰り返し演算
||R|Maxima|
|Σi^2 (i=1～10)|sum((1:10)^2)|sum(i^2,i,1,10)|

*Graphics
||R|Maxima|
|2Dグラフ sin(x)を-2pi～2piで|plot(sin,xlim=c(-2*pi,2*pi))|plot2d(sin(x),[x,-2*%pi,2*%pi])|
|3Dグラフ sin(x^2+y^2)|簡単な方法なし？|plot3d(cos(x^2+y^2),[x,-1,1],[y,-1,1])|
|3Dグラフ 球の媒介変数表示||plot3d([cos(s)*cos(t),cos(s)*sin(t),sin(s)],[s,-%pi/2,%pi/2],[t,0,2*%pi])|    </description>
    <dc:date>2010-04-30T19:23:14+09:00</dc:date>
    <utime>1272622994</utime>
  </item>
    <item rdf:about="https://w.atwiki.jp/pianissimo/pages/28.html">
    <title>math</title>
    <link>https://w.atwiki.jp/pianissimo/pages/28.html</link>
    <description>
      [[Gaussian Fitting&gt;http://nuclear.phys.tohoku.ac.jp/~ykoba/latex2html/gaussian-fitting/]]


問題
(x1, y1) (x2, y2) (x3, y3)の3点を通る放物線(2次関数)を求めよ。ただしx1≠x2, x1≠x3, x2≠x3は保証されているとする。

解答1
y=ax^2+bx+c
とおいて、ここに3点を代入すれば求められる。
かなり大変な計算になるが、
 a = {(y1-y2)/(x1-x2)-(y2-y3)/(x2-x3)}/(x1-x3)
 b = (y1-y2)/(x1-x2)-(x1+x2)a
 c = y1-ax1^2-bx1
となる。

解答2
3点を通る放物線は以下の式で書ける。[[ref&gt;http://detail.chiebukuro.yahoo.co.jp/qa/question_detail/q1113582701?fr=rcmd_chie_detail]]
y=a(x-x1)(x-x2)+b(x-x2)(x-x3)+c(x-x3)(x-x1)
ここへ3点を代入すると、簡単に、
 (x1,y1) → b = y1/{(x1-x2)(x1-x3)}
 (x2,y2) → c = y2/{(x2-x3)(x2-x1)}
 (x3,y3) → a = y3/{(x3-x1)(x3-x2)}
が求まる。

解答3
2次関数に限らず高次になっても応用しやすい最も一般的な解法は、行列を用いて連立方程式を解くことだろう。[[ref1&gt;http://cse.naro.affrc.go.jp/takezawa/r-tips/r/20.html]] [[ref2&gt;http://flashjp.com/column/gaussjordan.php]]    </description>
    <dc:date>2010-04-21T07:51:41+09:00</dc:date>
    <utime>1271803901</utime>
  </item>
    <item rdf:about="https://w.atwiki.jp/pianissimo/pages/16.html">
    <title>ウェブ・スクラップ</title>
    <link>https://w.atwiki.jp/pianissimo/pages/16.html</link>
    <description>
      - [[Vistaにおけるファイル同期&gt;http://www10.atwiki.jp/pianissimo/archive/20080511/59e8dbcce1c93f1356675b960e608fb2]] (2008-05-11 22:13:45)
- [[Vista usersが共有されるバグ&gt;http://www10.atwiki.jp/pianissimo/archive/20080512/7ccfd326e731052f17c312e8370dda01]] (2008-05-12 22:37:13)
- [[Excel 時間差を分だけで表示したいとき(75分など)...&gt;http://www10.atwiki.jp/pianissimo/archive/20080516/ecbd5d28874fad50cb58e5bdde6b43a6]] (2008-05-16 08:23:46)
- [[deinococcus radiodurans&gt;http://www10.atwiki.jp/pianissimo/archive/20080523/72e88529e3f5a17b4b175984332c670b]] (2008-05-23 13:19:07)
- [[ポリゴン内部にあることを判定する方法...&gt;http://www10.atwiki.jp/pianissimo/archive/20080709/7e9a088df172bf3639829d0466e81bc5]] (2008-07-09 05:27:37)
- [[Analyze format spec&gt;http://www10.atwiki.jp/pianissimo/archive/20080907/ade948c54682f9d2780b9010c50986b7]] (2008-09-07 13:03:09)
- [[dicom c#&gt;http://www10.atwiki.jp/pianissimo/archive/20090825/7be4258595ee8b0b3967a22292582627]] (2009-08-25 23:56:13)
- [[bisphosphonate&gt;http://www10.atwiki.jp/pianissimo/archive/20091031/cc2574444c504c6a3a09baedf93fefbf]] (2009-10-31 11:18:25)
- [[MDP&gt;http://www10.atwiki.jp/pianissimo/archive/20091031/51ad420d65c0287b01583504d263ab8d]] (2009-10-31 11:21:14)
- [[HMDP&gt;http://www10.atwiki.jp/pianissimo/archive/20091031/ad037ce3b51b6ae02f60afca76d7cf36]] (2009-10-31 11:21:28)
- [[テスト&gt;http://www10.atwiki.jp/pianissimo/archive/20100123/9a57a42091cfbdca725acb9cefab0527]] (2010-01-23 09:49:33)
- [[arial vs hervetica&gt;http://www10.atwiki.jp/pianissimo/archive/20100223/77dcaea356f2915faec88bbad8253589]] (2010-02-23 09:25:15)
- [[arial vs hervetica 2&gt;http://www10.atwiki.jp/pianissimo/archive/20100223/4a1b45b8483d4988e7c0a08a936f0cba]] (2010-02-23 09:27:42)
- [[filter&gt;http://www10.atwiki.jp/pianissimo/archive/20100309/9e13f816245455cd3cccc0d99a04e824]] (2010-03-09 12:51:43)
- [[parzan パルザン窓&gt;http://www10.atwiki.jp/pianissimo/archive/20100309/32d8e793252ca63fbbf5938c825c674f]] (2010-03-09 12:54:37)
#archive_log    </description>
    <dc:date>2010-03-09T12:54:37+09:00</dc:date>
    <utime>1268106877</utime>
  </item>
    <item rdf:about="https://w.atwiki.jp/pianissimo/pages/12.html">
    <title>テストページ</title>
    <link>https://w.atwiki.jp/pianissimo/pages/12.html</link>
    <description>
      
&lt;p&gt; &lt;/p&gt;
&lt;p&gt; &lt;/p&gt;
&lt;hr /&gt;&lt;p&gt; &lt;/p&gt;
    </description>
    <dc:date>2010-01-23T09:46:28+09:00</dc:date>
    <utime>1264207588</utime>
  </item>
    <item rdf:about="https://w.atwiki.jp/pianissimo/pages/26.html">
    <title>medcon</title>
    <link>https://w.atwiki.jp/pianissimo/pages/26.html</link>
    <description>
      *MedConの使い方

**インストール
[[MedCon website&gt;http://xmedcon.sourceforge.net/]]
このサイトで、左のコラムにあるDownloadを選択し、win32 versionをDLすればよい。
2009/12/30現在、xmedcon-0.10.5-1-win32.zipというファイルが取得される。
解凍して、setup.exeを実行する。
インストール後に表示されるreadme.txtの要点は、command lineで実行したいなら、\bin\ディレクトリにある実行ファイルにパスを通しておくこと。絶対パスでいうと、
C:\Program Files\XMedCon\bin
である。

**コマンドラインでの使用
そのパスで、medcon.exe -hを実行すると、下記のヘルプファイルが手に入る。

全文は最後に載せるとして、基本的な使い方を書いておきたい。

まず、もっとも単純な使い方は、ヘッダ内容の確認である。

例 medcon -f aaaaaa.v
aaaaaa.vというecat7ファイルを開いて、そのヘッダを表示する。
dicomヘッダも表示できるのでなかなかよい。



むしろ、MedConの主な使用目的は、ファイルのコンバートであろう。

例 medcon -c dicom -f *.v
ワイルドカードで指定した&quot;.v&quot;で終わるファイルをすべて、DICOM形式に変換する。
DICOMは1ボリューム1ファイルなので、いつもの感覚と違うが、スライス数は、(0054,0081) US[1] NumberOfSlices: 63 (2 bytes)のタグに入っているようだ。
いつもの1スライス1ファイルにするには、オプションで指定してやればよい。
例 medcon -c dicom -split3d -f *.v


出力形式には以下のものが選べる。

  &quot;ascii&quot; = Raw Ascii     (.asc)
  &quot;bin&quot;   = Raw Binary    (.bin)
  &quot;acr&quot;   = Acr/Nema      (.ima)
  &quot;anlz&quot;  = Analyze       (.hdr)+(.img)
  &quot;conc&quot;  = Concorde/uPET (.img.hdr)
  &quot;dicom&quot; = DICOM         (.dcm)
  &quot;ecat6&quot; = CTI ECAT 6    (.img)
  &quot;ecat7&quot; = CTI ECAT 7    (.v)
  &quot;gif&quot;   = Gif89a        (.gif)
  &quot;intf&quot;  = InterFile     (.h33)+(.i33)
  &quot;inw&quot;   = INW (RUG)     (.im)
  &quot;nifti&quot; = NIfTI         (.nii)
  &quot;png&quot;   = PNG           (.png)

pngを選択すると、全スライスがpng形式のモノクロ（黒→白）画像で出力される。



**ヘルプファイル全文

Usage:

  medcon [options] -f &lt;files&gt; ...

Flags:

  -f, --file, --files  file or list of files to handle

General: [-i -e -r -w] [-p -pa|-c &lt;format&gt; ...] [-o &lt;basename&gt;]

  -e, --extract            extract images from file
  -i, --interactive        read files after user input
  -o, --output-name        output name set from command-line
  -p, --print-values       print values of specified pixels
  -pa, --print-all-values  print all values without asking
  -r, --rename-file        rename file after user input
  -w, --overwrite-files    always overwrite files

  -c, --convert            give list of conversion &quot;&lt;format&gt;&quot; strings:

                &quot;ascii&quot; = Raw Ascii     (.asc)
                &quot;bin&quot;   = Raw Binary    (.bin)
                &quot;acr&quot;   = Acr/Nema      (.ima)
                &quot;anlz&quot;  = Analyze       (.hdr)+(.img)
                &quot;conc&quot;  = Concorde/uPET (.img.hdr)
                &quot;dicom&quot; = DICOM         (.dcm)
                &quot;ecat6&quot; = CTI ECAT 6    (.img)
                &quot;ecat7&quot; = CTI ECAT 7    (.v)
                &quot;gif&quot;   = Gif89a        (.gif)
                &quot;intf&quot;  = InterFile     (.h33)+(.i33)
                &quot;inw&quot;   = INW (RUG)     (.im)
                &quot;nifti&quot; = NIfTI         (.nii)
                &quot;png&quot;   = PNG           (.png)

Pixels: [-n] [-nf] [-qs|-qc|-q] [-b8|-b16[.12]] [-big|little]
        [-si=&lt;slope&gt;:&lt;intercept&gt;] [-cw=&lt;centre&gt;:&lt;width&gt;]

  -n, --negatives           enable negative pixels
  -nf, --norm-over-frames   normalize values over each frames
  -q, --quantitation        quantitation using all factors (-qc)
  -qs, --quantification     quantification (use one scale factor )
  -qc, --calibration        calibration    (use two scale factors)
  -b8, --unsigned-char      write unsigned char pixels  (8  bits)
  -b16, --signed-short      write signed short integers (16 bits)
  -b16.12                   write unsigned shorts, only 12 bits used
  -big, --big-endian        write files in big    endian
  -little, --little-endian  write files in little endian
  -si                       force slope/intercept rescaling
  -cw                       force specified contrast remapping

Fallback Read Format: [-fb-none|-fb-anlz|-fb-conc|-fb-ecat|fb-dicom]

  -fb-none, --without-fallback  fallback disabled
  -fb-anlz, --fallback-analyze  fallback on Analyze (SPM)
  -fb-conc, --fallback-concorde fallback on Concorde uPET
  -fb-ecat, --fallback-ecat     fallback on ECAT 6.4
  -fb-dicom, --fallback-dicom   fallback on DICOM 3.0

Slices Transform: [-fh -fv] [-rs -cs -cu] [-sqr | -sqr2] [-crop=&lt;X&gt;:&lt;Y&gt;:&lt;W&gt;:&lt;H&gt;]

                [-pad | -padtl | -padbr]
  -fh, --flip-horizontal        flip images horizontally (along x-axis)
  -fv, --flip-vertical          flip images vertically   (along y-axis)
  -sqr, --make-square           make square images (lagest dimension)
  -sqr2, --make-square-two      make square images (nearest power)
  -crop, --crop-images          crop image dimensions
  -rs, --reverse-slices         reverse slices sequence
  -cs, --cine-sorting           apply cine sorting
  -cu, --cine-undo              undo  cine sorting
  -pad, --pad-around            padding symmetrical around image
  -padtl, --pad-top-left        padding before first row and column
  -padbr, --pad-bottom-right    padding after last row and column (default)

Color Remap: [-24 | -8 [-g -dith -mh|-mr|-mi|-mc|-lut &lt;file&gt;]]

  -24, --true-color             color mode of 24 bits RGB triplets
  -8, --indexed-color           color mode of  8 bits indexed colormap
  -dith, --dither-color         apply dithering on color reduction

  -g, --make-gray               remap images to grayscale
  -mh, --map-hotmetal           load colormap hotmetal
  -mr, --map-rainbow            load colormap rainbow
  -mi, --map-inverted           load colormap gray inverted
  -mc, --map-combined           load colormap combined (gray/rainbow)
  -lut, --load-lut              load specified LUT colormap

Extras: [-alias -noprefix -preacq -preser -uin]
        [[-splits | -splitf] | [-stacks | -stackf]]

  -alias, --alias-naming        output name based  on patient/study id&#039;s
  -noprefix, --without-prefix   output name without default prefix
  -preacq, --prefix-acquisition use acquisition number as filename prefix
  -preser, --prefix-series      use series      number as filename prefix
  -uin, --use-institution-name  override default name of institution

  -split3d, --split-slices      split single image slices in separate files
  -split4d, --split-frames      split volume time frames  in separate files
  -stack3d, --stack-slices      stack single image slices into one 3D file
  -stack4d, --stack-frames      stack volume time  frames into one 4D file

Format Ecat/Matrix: [-byframe]

  -byframe, --sort-by-frame     sort ECAT images by frame (not anatomical)

Format Analyze: [-spm -optspm]

  -spm, --analyze-spm           use analyze format for SPM software

  -optspm, --options-spm        ask for SPM related options

Format DICOM:

 a) general: [-cw=&lt;center&gt;:&lt;width&gt;]
             [-contrast] [-gap] [-implicit] [-nometa]

  -contrast, --enable-contrast  enable support for contrast changes
  -gap, --spacing-true-gap      slice spacing is true gap or overlap
  -implicit, --write-implicit   output file as implicit little endian
  -nometa, --write-without-meta output file without (part 10) meta header
  -cw                           force window center/width contrast

 b) mosaic: [-mosaic | -fmosaic=&lt;W&gt;x&lt;H&gt;x&lt;N&gt; [-interl] [-mfixv]]

  -mosaic, --enable-mosaic      enable mosaic by &quot;detected&quot; stamps layout
  -fmosaic, --force-mosaic      force  mosaic by predefined stamps layout
  -mfixv, --mosaic-fix-voxel    rescale voxel sizes by mosaic factor
  -interl, --mosaic-interlaced  consider mosaic stamp slices as interlaced


Format Gif89a: [-optgif]

  -optgif, --options-gif        ask for GIF related options

Format InterFile: [-skip1 -nopath -one]

  -skip1, --skip-preview-slice  skip first preview slice
  -nopath, --ignore-path        ignore path in &#039;name of data file&#039; key
  -one, --single-file           write header and image to same file


Patient/Slice/Study: [-anon|-ident] [-vifi]

  -anon, --anonymous            make patient/study anonymous
  -ident, --identify            ask for patient/study information
  -vifi, --edit-fileinfo        edit internal entries (images/slice/origent)

Reslicing: [-tra|-sag|-cor]

  -tra, --tranverse             reslice images transverse
  -sag, --sagittal              reslice images sagittal
  -cor, --coronal               reslice images coronal

Debug/Mode: [-d -v -db -hackacr -ean]

  -d, --debug                   give debug information (printout FI)
  -s, --silent                  force silent mode, suppress all messages
  -v, --verbose                 run in verbose mode
  -db, --database               print database info
  -ean, --echo-alias-name       echo alias name on screen

  -hackacr, --hack-acrtags      try to locate and interpret acr tags in file    </description>
    <dc:date>2009-12-30T20:41:32+09:00</dc:date>
    <utime>1262173292</utime>
  </item>
    <item rdf:about="https://w.atwiki.jp/pianissimo/pages/13.html">
    <title>Sleipnir活用法2</title>
    <link>https://w.atwiki.jp/pianissimo/pages/13.html</link>
    <description>
          </description>
    <dc:date>2008-05-10T18:40:07+09:00</dc:date>
    <utime>1210412407</utime>
  </item>
    <item rdf:about="https://w.atwiki.jp/pianissimo/pages/11.html">
    <title>sleipnir活用法</title>
    <link>https://w.atwiki.jp/pianissimo/pages/11.html</link>
    <description>
          </description>
    <dc:date>2008-05-06T14:07:37+09:00</dc:date>
    <utime>1210050457</utime>
  </item>
    <item rdf:about="https://w.atwiki.jp/pianissimo/pages/10.html">
    <title>プラグイン/コメント</title>
    <link>https://w.atwiki.jp/pianissimo/pages/10.html</link>
    <description>
      * コメントプラグイン
@wikiのwikiモードでは
 #comment()
と入力することでコメントフォームを簡単に作成することができます。
詳しくはこちらをご覧ください。
＝＞http://atwiki.jp/guide/17_60_ja.html


-----
たとえば、#comment() と入力すると以下のように表示されます。

#comment    </description>
    <dc:date>2008-05-06T14:02:13+09:00</dc:date>
    <utime>1210050133</utime>
  </item>
  </rdf:RDF>
