mrcFETnormalizeBySD

• 報告者: kinoshita-eos
• 担当者: tacyas
• チケットの種類: バグ
• 状況: オープン [担当者決定済み]
• コンポーネント: base-v2.0
• マイルストーン: Small Toolsのドキュメント作成 (完了済み)
• 優先度: 5 - 中
• 重要度: 5 - 中
• 解決法: なし
• ファイル: なし

詳細 返信

【問題内容】
	FETファイルリストの(i, 0, 0)のデータ(θ=0 の幅、歪度、尖度)のみで正規化する動作になっています。
	FETファイル内の全てのθについて正規化する必要はないのでしょうか。
	(仕様である場合はドキュメントの方に追記します)

【ソースコード(一部)】
ファイル名: /src/Objects/DataManip/mrcImage/src/lmrcFETnormalizeBySD.c
(1) 平均、標準偏差の計算(list: 入力ファイルリスト、n: 特徴量の種類(=3)、count: ファイル数)
	関数名: calcSDofFET
	(146行目付近)
	while(NULL!=stringGetFromFile(s,"",list,stdout,1)){
		mrcFileRead(&cFile, s, "in lmrcFETnormalizeBySD", 0);
		for(i=0;i<n;i++){
			mrcPixelDataGet(&cFile,i,0,0,&data, mrcPixelRePart, mrcPixelHowNearest);
			ave[i] += data/(double)count;
		}
	}
	fseek(list, 0L, SEEK_SET);
	while(NULL!=stringGetFromFile(s,"",list,stdout,1)){
		mrcFileRead(&cFile, s, "in lmrcFETnormalizeBySD", 0);
		for(i=0;i<n;i++){
			mrcPixelDataGet(&cFile,i,0,0,&data, mrcPixelRePart, mrcPixelHowNearest);
			sd[i] += pow((data-ave[i]),2)/(double)count;
		}
	}
	for(i=0;i<n;i++)
		sd[i] = pow(sd[i],0.5);

(2) 正規化処理(list: 入力ファイルリスト、n: 特徴量の種類(=3))
	関数名: _normalizeBySDandDataWrite
	(119行目付近)
	while(NULL!=stringGetFromFile(s,"",list,stdout,1)){
		mrcFileRead(&cFile, s, "in normalizeBySDandDataWrite", 0);
		for(i=0;i<n;i++){
			mrcPixelDataGet(&cFile,i,0,0,&data, mrcPixelRePart, mrcPixelHowNearest);
			data -= ave[i];
			data /= sd[i];
			mrcPixelDataSet(&nord,i,1,0,data, mrcPixelRePart);
			if(info.flagOut) mrcPixelDataSet(info.Out,i,l,0,data, mrcPixelRePart);
		}
		mrcStatDataSet(&nord,0);
		strcat(s,".nord");
		mrcFileWrite(&nord,s,"in normalizeBySDandDataWrite", 0);
		l++;
	}
	if(info.flagOut) mrcStatDataSet(info.Out,0);

【補足】
	(1)では154, 162行目付近、(2)では123行目付近にて
	各ファイルの(i, 0, 0)のデータを得ていて、i毎に平均、標準偏差の計算、正規化を行っているようです。
	各ファイルはFETファイルですので、(i, 0, 0)のデータはθ=0でのシノグラムの(幅、歪度、尖度)に対応していると思います。

チケットの履歴 (3 件中 3 件表示)

2014-03-12 17:05 更新者: kinoshita-eos

• 新しいチケット "mrcFETnormalizeBySD" が作成されました

2014-05-01 04:00 更新者: tacyas

• コンポーネント が (未割り当て) から base-v2.0 に更新されました
• マイルストーン が (未割り当て) から Small Toolsのドキュメント作成 (完了済み) に更新されました

2014-09-01 09:01 更新者: None

• マイルストーン Small Toolsのドキュメント作成 (完了済み) は終了しました