<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<html><head><meta name="robots" content="noindex">
<meta http-equiv="Content-Type" content="text/html;charset=iso-8859-1">
<title>ITK: itkGaussianDensityFunction.h Source File</title>
<link href="DoxygenStyle.css" rel="stylesheet" type="text/css">
</head><body bgcolor="#ffffff">
<!-- Section customized for INSIGHT : Tue Jul 17 01:02:45 2001 -->
<center>
<a href="index.html" class="qindex">Main Page</a>
<a href="modules.html" class="qindex">Groups</a>
<a href="namespaces.html" class="qindex">Namespace List</a>
<a href="hierarchy.html" class="qindex">Class Hierarchy</a>
<a href="classes.html" class="qindex">Alphabetical List</a>
<a href="annotated.html" class="qindex">Compound List</a>
<a href="files.html" class="qindex">File
List</a>
<a href="namespacemembers.html" class="qindex">Namespace Members</a>
<a href="functions.html" class="qindex">Compound Members</a>
<a href="globals.html" class="qindex">File Members</a>
<a href="pages.html" class="qindex">Concepts</a></center>
<!-- Generated by Doxygen 1.5.9 -->
<div class="navpath"><a class="el" href="dir_02e23d841c198393a6e961e9a0678fff.html">Code</a> » <a class="el" href="dir_2c33022d19e8077162952ea9ed57656d.html">Numerics</a> » <a class="el" href="dir_e82fe55e9fab4124d6625d716e091a5f.html">Statistics</a>
</div>
<div class="contents">
<h1>itkGaussianDensityFunction.h</h1><a href="itkGaussianDensityFunction_8h.html">Go to the documentation of this file.</a><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 <span class="comment">/*=========================================================================</span>
<a name="l00002"></a>00002 <span class="comment"></span>
<a name="l00003"></a>00003 <span class="comment"> Program: Insight Segmentation & Registration Toolkit</span>
<a name="l00004"></a>00004 <span class="comment"> Module: $RCSfile: itkGaussianDensityFunction.h,v $</span>
<a name="l00005"></a>00005 <span class="comment"> Language: C++</span>
<a name="l00006"></a>00006 <span class="comment"> Date: $Date: 2005-07-26 15:54:55 $</span>
<a name="l00007"></a>00007 <span class="comment"> Version: $Revision: 1.15 $</span>
<a name="l00008"></a>00008 <span class="comment"></span>
<a name="l00009"></a>00009 <span class="comment"> Copyright (c) Insight Software Consortium. All rights reserved.</span>
<a name="l00010"></a>00010 <span class="comment"> See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.</span>
<a name="l00011"></a>00011 <span class="comment"></span>
<a name="l00012"></a>00012 <span class="comment"> This software is distributed WITHOUT ANY WARRANTY; without even </span>
<a name="l00013"></a>00013 <span class="comment"> the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR </span>
<a name="l00014"></a>00014 <span class="comment"> PURPOSE. See the above copyright notices for more information.</span>
<a name="l00015"></a>00015 <span class="comment"></span>
<a name="l00016"></a>00016 <span class="comment">=========================================================================*/</span>
<a name="l00017"></a>00017 <span class="preprocessor">#ifndef __itkGaussianDensityFunction_h</span>
<a name="l00018"></a>00018 <span class="preprocessor"></span><span class="preprocessor">#define __itkGaussianDensityFunction_h</span>
<a name="l00019"></a>00019 <span class="preprocessor"></span>
<a name="l00020"></a>00020 <span class="preprocessor">#include "<a class="code" href="itkArray_8h.html">itkArray.h</a>"</span>
<a name="l00021"></a>00021 <span class="preprocessor">#include "<a class="code" href="itkVariableSizeMatrix_8h.html">itkVariableSizeMatrix.h</a>"</span>
<a name="l00022"></a>00022 <span class="preprocessor">#include "vnl/algo/vnl_matrix_inverse.h"</span>
<a name="l00023"></a>00023 <span class="preprocessor">#include "vnl/algo/vnl_determinant.h"</span>
<a name="l00024"></a>00024 <span class="preprocessor">#include "vnl/vnl_math.h"</span>
<a name="l00025"></a>00025
<a name="l00026"></a>00026 <span class="preprocessor">#include "<a class="code" href="itkMatrix_8h.html">itkMatrix.h</a>"</span>
<a name="l00027"></a>00027 <span class="preprocessor">#include "<a class="code" href="itkDensityFunction_8h.html">itkDensityFunction.h</a>"</span>
<a name="l00028"></a>00028
<a name="l00029"></a>00029 <span class="keyword">namespace </span>itk{
<a name="l00030"></a>00030 <span class="keyword">namespace </span>Statistics{
<a name="l00031"></a>00031
<a name="l00054"></a>00054 <span class="keyword">template</span>< <span class="keyword">class</span> TMeasurementVector >
<a name="l00055"></a><a class="code" href="classitk_1_1Statistics_1_1GaussianDensityFunction.html">00055</a> <span class="keyword">class </span><a class="code" href="itkWin32Header_8h.html#7d4384322557205d75f0719c012dded6">ITK_EXPORT</a> <a class="code" href="classitk_1_1Statistics_1_1GaussianDensityFunction.html" title="GaussianDensityFunction class represents Gaussian Density Function.">GaussianDensityFunction</a> :
<a name="l00056"></a>00056 <span class="keyword">public</span> <a class="code" href="classitk_1_1Statistics_1_1DensityFunction.html" title="DensityFunction class defines common interfaces for density functions.">DensityFunction</a>< TMeasurementVector >
<a name="l00057"></a>00057 {
<a name="l00058"></a>00058 <span class="keyword">public</span>:
<a name="l00060"></a><a class="code" href="classitk_1_1Statistics_1_1GaussianDensityFunction.html#e5d2b18d6161d2d2644c370c1a58a33f">00060</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1Statistics_1_1GaussianDensityFunction.html" title="GaussianDensityFunction class represents Gaussian Density Function.">GaussianDensityFunction</a> <a class="code" href="classitk_1_1LightObject.html" title="Light weight base class for most itk classes.">Self</a>;
<a name="l00061"></a><a class="code" href="classitk_1_1Statistics_1_1GaussianDensityFunction.html#c91b2acb23e72a774531a2c4645f36c0">00061</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1Statistics_1_1DensityFunction.html" title="DensityFunction class defines common interfaces for density functions.">DensityFunction< TMeasurementVector ></a> <a class="code" href="classitk_1_1LightObject.html" title="Light weight base class for most itk classes.">Superclass</a> ;
<a name="l00062"></a><a class="code" href="classitk_1_1Statistics_1_1GaussianDensityFunction.html#8fc320bc6eafc0efd0b299bbcd5ef1f2">00062</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1SmartPointer.html">SmartPointer<Self></a> <a class="code" href="classitk_1_1SmartPointer.html">Pointer</a>;
<a name="l00063"></a><a class="code" href="classitk_1_1Statistics_1_1GaussianDensityFunction.html#c2b877f0192fa6b4597e36c2de6c9d94">00063</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1SmartPointer.html">SmartPointer<const Self></a> <a class="code" href="classitk_1_1SmartPointer.html">ConstPointer</a>;
<a name="l00064"></a>00064
<a name="l00066"></a>00066 <a class="code" href="itkMacro_8h.html#574c1f86661d409dd32c742307a0e96a">itkTypeMacro</a>(<a class="code" href="classitk_1_1Statistics_1_1GaussianDensityFunction.html" title="GaussianDensityFunction class represents Gaussian Density Function.">GaussianDensityFunction</a>, <a class="code" href="classitk_1_1Statistics_1_1DensityFunction.html" title="DensityFunction class defines common interfaces for density functions.">DensityFunction</a>);
<a name="l00067"></a><a class="code" href="classitk_1_1Statistics_1_1GaussianDensityFunction.html#40fff3ca1a8e0ef24ad000d98fe5614f">00067</a> <a class="code" href="itkMacro_8h.html#0fdc15685afa271857c57e7cd482f30c">itkNewMacro</a>(<a class="code" href="classitk_1_1LightObject.html" title="Light weight base class for most itk classes.">Self</a>);
<a name="l00069"></a>00069
<a name="l00071"></a>00071 <span class="keyword">typedef</span> TMeasurementVector MeasurementVectorType ;
<a name="l00072"></a>00072
<a name="l00074"></a><a class="code" href="classitk_1_1Statistics_1_1GaussianDensityFunction.html#674b67949b2b88f4f9f960dc4078fc16">00074</a> <span class="keyword">typedef</span> <span class="keyword">typename</span> Superclass::MeasurementVectorSizeType MeasurementVectorSizeType;
<a name="l00075"></a>00075
<a name="l00077"></a><a class="code" href="classitk_1_1Statistics_1_1GaussianDensityFunction.html#ebf9fbf8ffabf7a2b0d0ef10f13e0ef6">00077</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1Array.html">Array< double ></a> <a class="code" href="classitk_1_1Array.html">MeanType</a>;
<a name="l00078"></a>00078
<a name="l00080"></a><a class="code" href="classitk_1_1Statistics_1_1GaussianDensityFunction.html#0f52f6476added0740cd763375c8f532">00080</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1VariableSizeMatrix.html">VariableSizeMatrix< double ></a> <a class="code" href="classitk_1_1VariableSizeMatrix.html">CovarianceType</a>;
<a name="l00081"></a>00081
<a name="l00083"></a><a class="code" href="classitk_1_1Statistics_1_1GaussianDensityFunction.html#b37d83838d05d61fbea1e676edfa9608">00083</a> <span class="keywordtype">void</span> SetMean( <span class="keyword">const</span> <a class="code" href="classitk_1_1Array.html">MeanType</a> * mean )
<a name="l00084"></a>00084 {
<a name="l00085"></a>00085 <span class="keywordflow">if</span>( this->GetMeasurementVectorSize() )
<a name="l00086"></a>00086 {
<a name="l00087"></a>00087 <a class="code" href="classitk_1_1MeasurementVectorTraits.html#1f802b66e64e80af8b16b310a966ff65">MeasurementVectorTraits::Assert</a>(mean, this->GetMeasurementVectorSize(),
<a name="l00088"></a>00088 <span class="stringliteral">"GaussianDensityFunction::SetMean Size of measurement vectors in the sample must the same as the size of the mean."</span> );
<a name="l00089"></a>00089 }
<a name="l00090"></a>00090 <span class="keywordflow">else</span>
<a name="l00091"></a>00091 {
<a name="l00092"></a>00092 this->SetMeasurementVectorSize( mean-><a class="code" href="classitk_1_1Array.html#326c75c7d496c6471e35e9eff7f691a8">Size</a>() );
<a name="l00093"></a>00093 }
<a name="l00095"></a>00095
<a name="l00096"></a>00096 <span class="keywordflow">if</span> ( m_Mean != mean)
<a name="l00097"></a>00097 {
<a name="l00098"></a>00098 m_Mean = mean ;
<a name="l00099"></a>00099 this->Modified() ;
<a name="l00100"></a>00100 }
<a name="l00101"></a>00101 }
<a name="l00102"></a>00102
<a name="l00104"></a><a class="code" href="classitk_1_1Statistics_1_1GaussianDensityFunction.html#4199cceac3f4992ae7fd3927a4bea6f1">00104</a> <span class="keyword">const</span> <a class="code" href="classitk_1_1Array.html">MeanType</a> * GetMean()<span class="keyword"> const</span>
<a name="l00105"></a>00105 <span class="keyword"> </span>{ <span class="keywordflow">return</span> m_Mean ; }
<a name="l00106"></a>00106
<a name="l00110"></a>00110 <span class="keywordtype">void</span> SetCovariance(<span class="keyword">const</span> CovarianceType* cov);
<a name="l00111"></a>00111
<a name="l00113"></a>00113 <span class="keyword">const</span> CovarianceType* GetCovariance() <span class="keyword">const</span> ;
<a name="l00114"></a>00114
<a name="l00116"></a>00116 <span class="keywordtype">double</span> Evaluate(<span class="keyword">const</span> MeasurementVectorType &measurement) <span class="keyword">const</span> ;
<a name="l00117"></a>00117
<a name="l00118"></a>00118 <span class="keyword">protected</span>:
<a name="l00119"></a>00119 <a class="code" href="classitk_1_1Statistics_1_1GaussianDensityFunction.html" title="GaussianDensityFunction class represents Gaussian Density Function.">GaussianDensityFunction</a>(<span class="keywordtype">void</span>) ;
<a name="l00120"></a><a class="code" href="classitk_1_1Statistics_1_1GaussianDensityFunction.html#976306b462516e495d4a1cbd543482f4">00120</a> <span class="keyword">virtual</span> <a class="code" href="classitk_1_1Statistics_1_1GaussianDensityFunction.html#976306b462516e495d4a1cbd543482f4">~GaussianDensityFunction</a>(<span class="keywordtype">void</span>) {}
<a name="l00121"></a>00121 <span class="keywordtype">void</span> <a class="code" href="namespaceHardConnectedComponentImageFilter.html#920ae03048a5ad715d803ca2e92b52de">PrintSelf</a>(std::ostream& os, <a class="code" href="classitk_1_1Indent.html" title="Control indentation during Print() invocation.">Indent</a> indent) <span class="keyword">const</span>;
<a name="l00122"></a>00122
<a name="l00123"></a>00123 <span class="keyword">private</span>:
<a name="l00124"></a>00124 <span class="keyword">const</span> MeanType * m_Mean; <span class="comment">// mean</span>
<a name="l00125"></a>00125 <span class="keyword">const</span> CovarianceType * m_Covariance; <span class="comment">// covariance matrix</span>
<a name="l00126"></a>00126
<a name="l00127"></a>00127 <span class="comment">// inverse covariance matrix which is automatically calculated </span>
<a name="l00128"></a>00128 <span class="comment">// when covariace matirx is set. This speed up the GetProbability()</span>
<a name="l00129"></a>00129 CovarianceType m_InverseCovariance;
<a name="l00130"></a>00130
<a name="l00131"></a>00131 <span class="comment">// pre_factor which is automatically calculated </span>
<a name="l00132"></a>00132 <span class="comment">// when covariace matirx is set. This speeds up the GetProbability() </span>
<a name="l00133"></a>00133 <span class="keywordtype">double</span> m_PreFactor;
<a name="l00134"></a>00134
<a name="l00137"></a>00137 <span class="keywordtype">bool</span> m_IsCovarianceZero ;
<a name="l00138"></a>00138 };
<a name="l00139"></a>00139
<a name="l00140"></a>00140 } <span class="comment">// end of namespace Statistics</span>
<a name="l00141"></a>00141 } <span class="comment">// end namespace itk</span>
<a name="l00142"></a>00142
<a name="l00143"></a>00143 <span class="preprocessor">#ifndef ITK_MANUAL_INSTANTIATION</span>
<a name="l00144"></a>00144 <span class="preprocessor"></span><span class="preprocessor">#include "itkGaussianDensityFunction.txx"</span>
<a name="l00145"></a>00145 <span class="preprocessor">#endif</span>
<a name="l00146"></a>00146 <span class="preprocessor"></span>
<a name="l00147"></a>00147 <span class="preprocessor">#endif</span>
<a name="l00148"></a>00148 <span class="preprocessor"></span>
</pre></div></div>
<hr><address><small>
Generated at Thu May 7 22:47:48 2009 for ITK by <a href="http://www.stack.nl/~dimitri/doxygen/index.html"> <img
src="http://www.stack.nl/~dimitri/doxygen/doxygen.png" alt="doxygen"
align="middle" border=0 width=110 height=53>
</a> 1.5.9 written by <a href="mailto:dimitri@stack.nl">Dimitri van Heesch</a>,
© 1997-2000</small></address>
</body>
</html>
distrib
>
Mandriva
>
2010.0
>
i586
>
media
>
contrib-release
>
by-pkgid
>
d5e74628f0e673bb8680aebce32b2c04
>
files
>
15885
