void FilterLed::Execute	(	std::vector< Led > &	leds,
		int	xc,
		int	yc,
		double &	bright
	)		const

Definition at line 344 of file FilterLed.cc.

References ClusterFinder::FindCluster().

Referenced by MStarguider::ProcessFrame(), MCaos::Run(), and SetCut().

{
    const int x0 = max(xc-fBoxX, 0);
    const int y0 = max(yc-fBoxY, 0);
    const int x1 = min(xc+fBoxX, fW);
    const int y1 = min(yc+fBoxY, fH);

    const int wx = x1-x0;
    const int hy = y1-y0;

    double sum = 0;
    double sq  = 0;

    for (int x=x0; x<x1; x++)
        for (int y=y0; y<y1; y++)
        {
            uint8_t &b = fImg[y*fW+x];

            // Skip saturating pixels
            if (b>0xf0)
                continue;

            sum += b;
            sq  += b*b;
        }

    sum /= wx*hy;
    sq  /= wx*hy;

    bright=sum;

    
    // 254 because b<=max and not b<max
    const double sdev = sqrt(sq-sum*sum);
    const uint8_t   max  = sum+fCut*sdev>254 ? 254 : (uint8_t)(sum+fCut*sdev);

    //
    // clean image from noise
    // (FIXME: A lookup table could accelerate things...
    //
    for (int x=x0; x<x1; x++)
        for (int y=y0; y<y1; y++)
        {
            uint8_t &b = fImg[y*fW+x];
            if (b<=max)
                b = 0;
        }

    ClusterFinder find(fImg, fW, fH);
    find.FindCluster(leds, x0, y0, x1, y1);
}

Here is the call graph for this function:

Here is the caller graph for this function: