#ifndef _OTB_TRISKELE_MOI_ATTRIBUTES_TPP
#define _OTB_TRISKELE_MOI_ATTRIBUTES_TPP

inline
MoIAttributes::MoIAttributes (const Tree &tree, const AreaAttributes &areaAttributes, const XYAttributes &xyAttributes)
  : CompAttributeC<double> (tree) {
  compute (areaAttributes, xyAttributes);
}

inline
MoIAttributes::~MoIAttributes () {
}

template<typename PixelT>
inline void
MoIAttributes::cut (vector<vector<PixelT> > &allBands, const AttributeProfiles<PixelT> &attributeProfiles,
		   const vector<double> &thresholds) const {
  DEF_LOG ("MoIAttributes::cut", "thresholds:" << thresholds.size ());
  double maxValue = CompAttribute<double>::getMaxValue ();
  cerr << "moi max value:" << maxValue << endl;
  CompAttributeC<double>::cut (allBands, attributeProfiles, 0,
			       CompAttributeC<double>::getScaledThresholds (thresholds, maxValue));
}

inline void
MoIAttributes::compute (const AreaAttributes &areaAttributes, const XYAttributes &xyAttributes) {
  const AverageXY *xy = xyAttributes.getValues ();
  const DimImg *areas = areaAttributes.getValues ();
  computeSameCompLevel ([this, &xy, &areas] (const DimImg &parentId) {
      double tmpMoI = 0.;
      const double xa = xy[parentId].x;
      const double ya = xy[parentId].y;
      tree.forEachChildTI (parentId, [this, &tmpMoI, &xa, &ya, &xy, &areas, &parentId] (const bool &isLeaf, const DimImg &childId) {
	  if (isLeaf)
	    return;
	  const double dx = xy[childId].x - xa;
	  const double dy = xy[childId].y - ya;
	  tmpMoI += (dx*dx+dy*dy)*areas[childId];
	});
      const double card = areas[parentId];
      CompAttribute<double>::values[parentId] = tmpMoI/(card*card);
    });
}

#endif // _OTB_TRISKELE_MOI_ATTRIBUTES_TPP