void MStarguider::ProcessFrame ( const unsigned long  n, byte *  img, struct timeval *  tm  )

virtual

Reimplemented from PixClient.

Definition at line 1127 of file MStarguider.cc.

References FilterLed::DrawCircle(), MGImage::DrawColImg(), MGImage::DrawImg(), FilterLed::Execute(), fCaos, fDimData, fDisplay, fFindStarBox, fFindStarCut, fImage, FilterLed::FindStar(), fScreenshotName, fTPointMode, Ring::GetMag(), Led::GetMag(), MCaos::GetNumDetectedLEDs(), MCaos::GetNumDetectedRings(), Ring::GetPhi(), Ring::GetX(), Led::GetX(), Ring::GetY(), Led::GetY(), IDM_kCaosFilter, IDM_kFilter, IDM_kFindStar, IDM_kStretch, Interpolate(), FilterLed::MarkPoint(), MCaos::Run(), FilterLed::SetBox(), FilterLed::SetCut(), DimService::setData(), DimService::setQuality(), DimService::setTimestamp(), FilterLed::Stretch(), DimService::updateService(), and WritePNG().

{
    if (!Interpolate(n, img))
        return;

    if (fTPointMode==2)
    {
        fTPointMode=1;
        return;
    }

    byte cimg[768*576];
    memset(cimg, 0, 768*576);

    FilterLed f (img,  768, 576, 2.5); // 2.5
    FilterLed f2(cimg, 768, 576);      // former color 0xb0

    if (!fTPointMode && fScreenshotName.empty() && fDisplay->IsEntryChecked(IDM_kStretch))
        f.Stretch();

    // Visual Filter, whole FOV
    if (!fTPointMode && fDisplay->IsEntryChecked(IDM_kFilter))
    {
        vector<Led> leds;
        f.Execute(leds, 768/2, 576/2);
        for (auto it=leds.begin(); it!=leds.end(); it++)
            f.MarkPoint(*it);
    }

    // Find Center of Camera for Caos and Tpoints
    int numleds  = 0;
    int numrings = 0;
    Ring center(-1, -1);//(5, 5);

    if (fTPointMode || fDisplay->IsEntryChecked(IDM_kCaosFilter))
    {
        center   = fCaos->Run(img);
        numleds  = fCaos->GetNumDetectedLEDs();
        numrings = fCaos->GetNumDetectedRings();
    }

    //cout << "cx=" << center.GetX() << "   cy=" << center.GetY() << "   Nled=" << numleds << "   Nrings=" << numrings << endl;

    // Find Star at Center---for Tpoint Procedure
    Led star(-1, -1);
    if (center.GetX()>0 && center.GetY()>0)
    {
        if (fTPointMode || fDisplay->IsEntryChecked(IDM_kFindStar))
        {
            // Set search Paremeters (FIXME: Get them from user input!)
            f.SetCut(fFindStarCut+1.5);  // FindStar.CleaningLevel
            f.SetBox(fFindStarBox+12/*+80*/);  // FindStar.SizeBox

            // Try to find the star
            vector<Led> leds;
            f.FindStar(leds, (Int_t)center.GetX(), (Int_t)center.GetY(), true);

            // Check whether star found
            if (leds.size()>0)
            {
                //cout << "Found star @ " << flush;
                //leds[0].Print();
                f2.MarkPoint(leds[0].GetX(), leds[0].GetY(), 2<<2);
                star = leds[0];
            }
        }

        // DrawZoomImage(img);
        // DrawCosyImage(img);

        // Position corresponding to the camera center (53.2, 293.6)
        // Draw Circles around center of Camera
        if (fTPointMode || fDisplay->IsEntryChecked(IDM_kCaosFilter))
        {
            f2.DrawCircle(center, 0x0a);
            f2.DrawCircle(center,   7.0,
                          fDisplay->IsEntryChecked(IDM_kFindStar)?3:0xb0);
            //f2.DrawCircle(center, 115.0, 0x0a);
            //f2.DrawCircle(center, 230.0, 0x0a);
            //f2.DrawCircle(center, 245.0, 0x0a);
        }
    }

    if (fTPointMode ||
        fDisplay->IsEntryChecked(IDM_kCaosFilter) ||
        fDisplay->IsEntryChecked(IDM_kFindStar))
        fImage->DrawColImg(img, cimg);
    else
        fImage->DrawImg(img);

    if (fTPointMode && !fScreenshotName.empty())
    {
        WritePNG(fScreenshotName.c_str(), img, cimg);
        fScreenshotName = "";
        fTPointMode = 0;
        return;
    }


    if (star.GetX()<0 || star.GetY()<0 || fTPointMode==0)
        return;

    if (fTPointMode==1)
        fTPointMode=0;

    // Convert from Pixel to millimeter (1pix=2.6mm) [deg/pix / deg/mm = mm/pix]
    // Correct for abberation.

    //const float Dleds =     510; // 5.96344 deg
    //const float Dpix  = 2*237.58;

    // The DC reflector elongates light from off axis sources
    // This is a correction. It is 7% for MAGIC (1:1) and
    // less for FACT (1:1.4). This is an estimate from a
    // Orbit mode observation at 0.17deg distance to the
    // camera center [4.90m might also not be very accurate
    // depending on the position of the CCD camera]
    const double abberation  = 1.0638; //1.0713;
    const double sec_per_pix = 45.14;  //45.311;  (atan(510mm/2 / 4.90m) / 237.58)           // FACT LEDs

    const double conv = sec_per_pix/abberation;

    const double dx = star.GetX()-center.GetX();
    const double dy = star.GetY()-center.GetY();

    //const double dxx = - conv * (star.GetX()-center.GetX());
    //const double dyy =   conv * (star.GetY()-center.GetY());

    const double dphi = - center.GetPhi() * M_PI/180;

    // The sign is because the pixels are not counted in
    // both directions in the same direction as Zd/Az
    const double dxx = - conv * (dx*cos(dphi) - dy*sin(dphi));
    const double dyy =   conv * (dx*sin(dphi) + dy*cos(dphi));

    double arr[11] = {
        dxx, dyy,
        double(numleds), double(numrings),
        center.GetX(), center.GetY(), center.GetMag(),
        star.GetX(),   star.GetY(),   star.GetMag(),
        center.GetPhi()
    };

    fDimData.setData(arr, 11*sizeof(double));
    fDimData.setQuality(0);
    fDimData.setTimestamp(tm->tv_sec, tm->tv_usec/1000);
    fDimData.updateService();
}

Here is the call graph for this function: