Impala Documentation

template<class DstArrayT, class SrcArrayT>

void Impala::Core::Feature::WeibullNgbPnLoop< DstArrayT, SrcArrayT >::NextEl	(	int	x,
		int	y,
		SrcArithType	value
	)			`[inline]`

Processing one pixel.

Definition at line 123 of file WeibullNgbPnLoop.h.

00124     {
00125       SrcArithType dataPowerGamma;
00126       SrcArithType sumItem;
00127       SrcArithType data = value + 0.000001; 
00128       dataPowerGamma = pow(data, mGamma);
00129       int phase = mPhase%2;
00130         
00131       if(mVerbose)
00132         if (mGamma == 0.1 && phase == 0)
00133           std::cout << value << std::endl;
00134       //      if (mVerbose)
00135       //        std::cout << "      WeibullNgbPnLoop::next(" << x << "," << y << ","<< value << ") " << std::endl;
00136 
00137       if (mLast)
00138       {
00139         mData->SetValue(value, x*mNgbWidth + y, 0);
00140         return;
00141       } 
00142       switch (phase)
00143       {
00144       case 1 : 
00145         mSumDataPowerGamma += dataPowerGamma;
00146         mSumDataLog += dataPowerGamma*log(data);
00147         break;
00148         
00149       case 0 :
00150         sumItem = dataPowerGamma/mSumDataPowerGamma; 
00151         mF += log(mFrac*sumItem) - mFrac*(sumItem)*log(mFrac*sumItem);
00152         mFprime += (mSumDataPowerGamma/dataPowerGamma - mFrac*(log(mFrac*sumItem) + 1)) *
00153         (dataPowerGamma*(log(data)*mSumDataPowerGamma - mSumDataLog))/
00154           (mSumDataPowerGamma*mSumDataPowerGamma);
00155         break;
00156       }
00157         
00158     }