00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00048
00049
00050
00051 #ifndef _PARZENDISTRIBUTIONFORCE_H
00052 #define _PARZENDISTRIBUTIONFORCE_H
00053
00054 #include "ForceGenerator.h"
00055 #include "Image.h"
00056 #include "Mask.h"
00057 #include "LevelSetFunction.h"
00058 #include "Histogram.h"
00059 #include <vector>
00060 #include <iostream>
00061
00062 namespace lsseg {
00063
00064
00084 class ParzenDistributionForce : public ForceGenerator
00085
00086 {
00087 public:
00100 ParzenDistributionForce(bool multireg_mode = false) :
00101 multi_region_mode_(multireg_mode),
00102 num_force_bins_(256),
00103 img_(0),
00104 mask_(0)
00105 {}
00106
00124 explicit ParzenDistributionForce(const Image<double>* img,
00125 Mask* mask = 0,
00126 bool multireg_mode = false) ;
00127
00128 virtual ~ParzenDistributionForce() {}
00129
00130 virtual void init(const Image<double>* img, const Mask* mask = 0);
00131
00132 virtual void update(const LevelSetFunction& phi);
00133 virtual double force2D(int x, int y) const;
00134 virtual double force3D(int x, int y, int z) const ;
00135 virtual double force(size_t ix) const;
00136 virtual void force(LevelSetFunction&, const Mask* mask = 0) const;
00137
00145 void saveChannelDistribution(int channel, std::ostream& os, bool inside) const;
00146
00156 void fixChannelDistributionTo(int channel, std::istream& is, bool inside);
00157
00166 void fixDistributionToUniform(int channel, bool inside);
00167
00174 void unfixChannelDistribution(int channel, bool inside);
00175
00177 const Image<double>* const baseImage() const {return img_;}
00178
00179 private:
00180
00184 const bool multi_region_mode_;
00185
00197 const int num_force_bins_;
00198
00199
00200
00202 const Image<double>* img_;
00203
00206 const Mask* mask_;
00207
00210 std::vector<std::pair<bool, bool> > is_fixed_;
00211
00215 std::vector<std::pair<Histogram, Histogram> > hist_;
00216
00220 std::vector<std::pair<double, double> > channel_ranges_;
00221
00229 std::vector<double> force_bin_factor_;
00230
00238 std::vector<std::vector<double> > precalc_force_;
00239
00240
00241 inline void report() const;
00242
00248 inline int getForceBin(double val, int channel) const;
00249
00262 void precompute_force(int channel,
00263 const Histogram& ihist,
00264 const Histogram& ohist,
00265 std::vector<double>& forcevec);
00266
00272 double getCenterBinValue(int channel,int ix);
00273
00285 inline void get_histogram(int channel,
00286 const LevelSetFunction& phi,
00287 Histogram& inside_hist,
00288 Histogram& outside_hist) const;
00289
00292 Histogram makeDefaultHistogram(int channel) const ;
00293 };
00294
00295 };
00296
00297 #endif // _PARZENDISTRIBUTIONFORCE_H
00298