Kinoshita: /* Case: sD=50000 */

2014-07-22T23:48:33Z

‎Case: sD=50000

Kinoshita: /* Case: sD=1000 */

2014-07-22T23:48:18Z

‎Case: sD=1000

Kinoshita: /* Case: sD=1000 */

2014-07-22T23:47:50Z

‎Case: sD=1000

Kinoshita at 00:55, 18 July 2014

2014-07-18T00:55:42Z

Kinoshita: Created page with "'''mrcImageBilateralFilter''' ismrcImageファイルを平滑化するための Eos's Command. It is a method for proposed by Tomasi and Manduchi, to smooth while ..."

2014-07-18T00:42:23Z

Created page with "'''mrcImageBilateralFilter''' ismrcImageファイルを平滑化するための Eos's Command. It is a method for proposed by Tomasi and Manduchi, to smooth while ..."

New page

'''mrcImageBilateralFilter''' is[[mrcImage]]ファイルを平滑化するための [[Eos]]'s [[Command]].
It is a method for proposed by Tomasi and Manduchi, to smooth while maintaining the edge. It is a smoothing of weighted by the intensity difference of the image. （intensity difference, distance）をガウス分布などにより重みを計算していることにより、強度差がある領域をぼかすことなく維持できる。

== Algorithm ==

It has two weights Wx(xi, xj) and WI(Ii,Ij), they are defined as such Lorenz function or Gaussian function of |xi-xj|, |Ii-Ij|. That is because the difference between the image becomes large near the edges、weight of Smoothing in density band that themselves belong is applied to them.

<math>
f_i =
</math>

fi = Sigma_j(W(i,j)*Ij)/ Sigma_j(Wj)
,where W(i,j)=W_x(xi,xj)*W_I(Ii,Ij)

------- linfo.mode : 0
W_x(xi,xj) = exp(-|xi-xj|^2/2 sigmaSpace^2 )
W_I(Ii,Ij) = exp(-|Ii-Ij|^2/2 sigmaDensity^2)
------- linfo.mode : 1
W_x(xi,xj) = 1/(1+|xi-xj|^2/sigmaSpace^2 )
W_I(Ii,Ij) = 1/(1+|Ii-Ij|^2/sigmaDensity^2)

== List of option ==
=== Main option ===
<table border="1">
<div align="left">
<tr>
<th>Option</th>
<th>Essential/Optional</th>
<th>Description</th>
<th>Default</th>
</tr>
<tr>
<td>-i</td>
<td>Essential</td>
<td>Input：[[mrcImage]]</td>
<td>NULL</td>
</tr>
<tr>
<td>-o</td>
<td>Optional</td>
<td>Output：[[mrcImage]]</td>
<td>stdout(The standard output by default）</td>
</tr>
<tr>
<td>-sD</td>
<td>Optional</td>
<td>SigmaDensity(float) [ImageAmplitude]</td>
<td>2.0</td>
</tr>
<tr>
<td>-kTD</td>
<td>Optional</td>
<td>Kernel Type of Density (int) </td>
<td>0</td>
</tr>
<tr>
<td>-sS</td>
<td>Optional</td>
<td>SigmaSpace(float) [pixel]</td>
<td>1.0</td>
</tr>
<tr>
<td>-sS3</td>
<td>Optional</td>
<td>SigmaSpace 3D (x:float, y:float, z:float) [pixel]</td>
<td>1.0 1.0 1.0 </td>
</tr>
<tr>
<td>-kTS</td>
<td>Optional</td>
<td>Kernel Type of Space (int) </td>
<td>0</td>
</tr>
<tr>
<td>-kS</td>
<td>Optional</td>
<td>Kernel Size (int) [pixel]</td>
<td>1.0</td>
</tr>
<tr>
<td>-kS3</td>
<td>Optional</td>
<td>Kernel Size 3D (x:int, y:int, z:int) [pixel]</td>
<td>1.0 1.0 1.0 </td>
</tr>
<tr>
<td>-c</td>
<td>Optional</td>
<td>ConfigurationFile</td>
<td>NULL</td>
</tr>
<tr>
<td>-m</td>
<td>Optional</td>
<td>Mode</td>
<td>0</td>
</tr>
<tr>
<td>-h</td>
<td>Optional</td>
<td>Help</td>
<td>　</td>
</tr>
</div>
</table>

=== -m details ===
<table border="1">
<tr>
<th>Value</th>
<th>Description</th>
</tr>
<tr>
<td>0</td>
<td></td>
</tr>
</table>

===-kTD details===
<table border="1">
<div align="left">
<tr>
<th>Value</th>
<th>Description</th>
</tr>
<tr>
<td>0</td>
<td>Gaussian function</td>
</tr>
<tr>
<td>1</td>
<td>Lorenz function</td>
</tr>
</div>
</table>

===-kTS details===
<table border="1">
<tr>
<th>Value</th>
<th>Description</th>
</tr>
<tr>
<td>0</td>
<td>Gaussian function</td>
</tr>
<tr>
<td>1</td>
<td>Lorenz function</td>
</tr>
</table>

== Execution example ==
===[[Media:Input-1VOM-N-2D.mrc | Input file]](2D) image===
<table>
<tr>
<td>[[File:Input-1VOM-N-2D.png]]</td>
<td><p align="left">Min<br>
Max<br>
Mean<br>
SD<br>
SE<br></p>
</td>
<td><p align="left">-18651.7 (10, 1, 0)<br>
52942.7 (24, 39, 0)<br>
7214.87<br>
10067.6<br>
125.845<br></p>
</td>
</tr>
</table>

==== Example of options only essential ====
<table>
<tr>
<td>[[File:Outdata-mrcImageBilateralFilter.png]]</td>
<td><p align="left">Min<br>
Max<br>
Mean<br>
SD<br>
SE<br></p>
</td>
<td><p align="left">-18651.7 (10, 1, 0)<br>
52942.7 (24, 39, 0)<br>
7214.87<br>
10067.6<br>
125.845<br></p>
</td>
</tr>
</table>

==== Option -sD====
=====Case: sD=50000=====
<table>
<tr>
<td>[[File:Outdata-sD-mrcImageBilateralFilter.png]]</td>
<td><p align="left">Min<br>
Max<br>
Mean<br>
SD<br>
SE<br></p>
</td>
<td><p align="left">-3963.33 (76, 66, 0)<br>
36164.2 (26, 29, 0)<br>
7213.95<br>
8271.46<br>
103.393<br></p>
</td>
</tr>
</table>

<div>In the following, Setting atsD=50000.</div>

==== Option -kTD====
=====Case: kTD=1=====
<table>
<tr>
<td>[[File:Outdata-kTD-mrcImageBilateralFilter.png]]</td>
<td><p align="left">Min<br>
Max<br>
Mean<br>
SD<br>
SE<br></p>
</td>
<td><p align="left">-4063.09 (56, 4, 0)<br>
36357.6 (26, 29, 0)<br>
7212.18<br>
8280.1<br>
103.501<br></p>
</td>
</tr>
</table>

==== Option -sS====
=====Case: sS=10=====
<table>
<tr>
<td>[[File:Outdata-sS-mrcImageBilateralFilter.png]]</td>
<td><p align="left">Min<br>
Max<br>
Mean<br>
SD<br>
SE<br></p>
</td>
<td><p align="left">-3196.78 (38, 0, 0)<br>
32196.3 (25, 39, 0)<br>
7210.81<br>
8039.59<br>
100.495<br></p>
</td>
</tr>
</table>

==== Option -kTS====
=====Case: kTS=1, sS=10=====
<table>
<tr>
<td>[[File:Outdata-kTS-mrcImageBilateralFilter.png]]</td>
<td><p align="left">Min<br>
Max<br>
Mean<br>
SD<br>
SE<br></p>
</td>
<td><p align="left">-3198.73 (38, 0, 0)<br>
32222.2 (25, 39, 0)<br>
7210.86<br>
8040.82<br>
100.51<br></p>
</td>
</tr>
</table>

==== Option -kS====
=====Case: kS=10=====
<table>
<tr>
<td>[[File:Outdata-kS-mrcImageBilateralFilter.png]]</td>
<td><p align="left">Min<br>
Max<br>
Mean<br>
SD<br>
SE<br></p>
</td>
<td><p align="left">-3871.71 (76, 66, 0)<br>
35969.2 (26, 29, 0)<br>
7213.84<br>
8256.51<br>
103.206<br></p>
</td>
</tr>
</table>

===[[Media:Input-1VOM-N.mrc | Input file]](3D) image===
<table>
<tr>
<td><p align="Center">[[File:Input-1VOM-N.png]]<br>
xy-plane<br></p>
</td>
<td><p align="Center">[[File:Input1-1VOM-N.png]]<br>
yz-plane<br></p>
</td>
<td><p align="left">Min<br>
Max<br>
Mean<br>
SD<br>
SE<br></p>
</td>
<td><p align="left">-2390.42 (77, 14, 69)<br>
4374.62 (17, 62, 46)<br>
72.1487<br>
577.82<br>
0.722275<br></p>
</td>
</tr>
</table>

==== Option -sD====
=====Case: sD=1000=====
<table>
<tr>
<td><p align="Center">[[File:Outdata1-sD-mrcImageBilateralFilter.png]]<br>
xy-plane<br></p>
</td>
<td><p align="Center">[[File:Outdata2-sD-mrcImageBilateralFilter.png]]<br>
yz-plane<br></p>
</td>
<td><p align="left">Min<br>
Max<br>
Mean<br>
SD<br>
SE<br></p>
</td>
<td><p align="left">-1272.88 (77, 14, 69)<br>
4131.48 (17, 62, 46)<br>
67.6465<br>
267.559<br>
0.334449<br></p>
</td>
</tr>
</table>

<div>In the following, Setting atsD=1000で実行する</div>

==== Option -sS3====
=====Case: sS3=(1, 0.7, 1)=====
<table>
<tr>
<td><p align="Center">[[File:Outdata-sS3-mrcImageBilateralFilter.png]]<br>
xy-plane<br></p>
</td>
<td><p align="Center">[[File:Outdata1-sS3-mrcImageBilateralFilter.png]]<br>
yz-plane<br></p>
</td>
<td><p align="left">Min<br>
Max<br>
Mean<br>
SD<br>
SE<br></p>
</td>
<td><p align="left">-1459.87 (77, 14, 69)<br>
4154.41 (17, 62, 46)<br>
68.4469<br>
284.532<br>
0.355665<br></p>
</td>
</tr>
</table>

==== Option -kS3====
=====Case: kS3=(2, 2, 1)=====
<table>
<tr>
<td><p align="Center">[[File:Outdata-kS3-mrcImageBilateralFilter.png]]<br>
xy-plane<br></p>
</td>
<td><p align="Center">[[File:Outdata1-kS3-mrcImageBilateralFilter.png]]<br>
yz-plane<br></p>
</td>
<td><p align="left">Min<br>
Max<br>
Mean<br>
SD<br>
SE<br></p>
</td>
<td><p align="left">-1966.48 (78, 4, 26)<br>
4354.27 (17, 62, 46)<br>
71.3418<br>
396.845<br>
0.496056<br></p>
</td>
</tr>
</table>

== References ==
* C. Tomasi and R. Manduchi, Bilateral Filtering for Gray and Color Images, Proceedings of the 1998 IEEE International Conference on Computer Vision, Bombay, India

← Older revision		Revision as of 23:48, 22 July 2014
Line 353:		Line 353:


−	<div>In the following, Setting ~~atsD~~=1000</div>	+	<div>In the following, Setting at sD=1000</div>

	==== Option -sS3====		==== Option -sS3====

@@ Line 215: / Line 215: @@
-<div>In the following, Setting atsD=50000.</div>
+<div>In the following, Setting at sD=50000.</div>
-−
 ==== Option -kTD====

@@ Line 353: / Line 353: @@
-<div>In the following, Setting atsD=1000で実行する</div>
+<div>In the following, Setting atsD=1000</div>
-−
 ==== Option -sS3====

@@ Line 1: / Line 1: @@
-'''mrcImageBilateralFilter''' is[[mrcImage]]ファイルを平滑化するための [[Eos]]'s [[Command]].
+'''mrcImageBilateralFilter''' is [[Eos]]'s [[Command]] that performs smoothing [[mrcImage]] file.
-It is a method for proposed by Tomasi and Manduchi, to smooth while maintaining the edge. It is a smoothing of weighted by the intensity difference of the image. （intensity difference, distance）をガウス分布などにより重みを計算していることにより、強度差がある領域をぼかすことなく維持できる。
+It is a method for proposed by Tomasi and Manduchi, to smooth while maintaining the edge. It is a smoothing of weighted by the intensity difference of the image. Since weights are calculated by Gauss distribution for the two deviations(intensity difference, distance), region in which intensity difference is maintained without blurred.
 == Algorithm ==

mrcImageBilateralFilter - Revision history

Kinoshita: /* Case: sD=50000 */

Kinoshita: /* Case: sD=1000 */

Kinoshita: /* Case: sD=1000 */

Kinoshita at 00:55, 18 July 2014

Kinoshita: Created page with "'''mrcImageBilateralFilter''' ismrcImageファイルを平滑化するための Eos's Command. It is a method for proposed by Tomasi and Manduchi, to smooth while ..."