| 1 | #undef EXPERT
|
|---|
| 2 | #define FACT
|
|---|
| 3 |
|
|---|
| 4 | #include "MStarguider.h"
|
|---|
| 5 |
|
|---|
| 6 | #ifdef HAVE_PNG
|
|---|
| 7 | #include <png.h>
|
|---|
| 8 | #endif
|
|---|
| 9 | #include <iostream>
|
|---|
| 10 | #include <fstream>
|
|---|
| 11 |
|
|---|
| 12 | #include <math.h>
|
|---|
| 13 |
|
|---|
| 14 | #include <TEnv.h>
|
|---|
| 15 | #include <TSystem.h>
|
|---|
| 16 |
|
|---|
| 17 | #include <TGMenu.h>
|
|---|
| 18 | #include <TGLabel.h>
|
|---|
| 19 | #include <TGButton.h>
|
|---|
| 20 | #include <TGSplitter.h> // TGHorizontal3DLine
|
|---|
| 21 | #include <TGTextEntry.h>
|
|---|
| 22 | #include <TGLayout.h>
|
|---|
| 23 |
|
|---|
| 24 | #include "MCaos.h"
|
|---|
| 25 | #include "MGImage.h"
|
|---|
| 26 | #include "Camera.h"
|
|---|
| 27 | #include "../src/DimSetup.h"
|
|---|
| 28 | #include "Led.h"
|
|---|
| 29 | //#include "Writer.h"
|
|---|
| 30 | #include "FilterLed.h"
|
|---|
| 31 | //#include "CaosFilter.h"
|
|---|
| 32 |
|
|---|
| 33 | using namespace std;
|
|---|
| 34 |
|
|---|
| 35 | enum {
|
|---|
| 36 | IDM_kFilter,
|
|---|
| 37 | IDM_kFindStar,
|
|---|
| 38 | IDM_kCaosFilter,
|
|---|
| 39 | IDM_kStarguider,
|
|---|
| 40 | IDM_kStretch,
|
|---|
| 41 | IDM_kInput,
|
|---|
| 42 | IDM_kSetup,
|
|---|
| 43 | IDM_kCut,
|
|---|
| 44 | IDM_kInterpol250,
|
|---|
| 45 | IDM_kInterpol125,
|
|---|
| 46 | IDM_kInterpol50,
|
|---|
| 47 | IDM_kInterpol25,
|
|---|
| 48 | IDM_kInterpol10,
|
|---|
| 49 | IDM_kInterpol5,
|
|---|
| 50 | IDM_kInterpol2,
|
|---|
| 51 | IDM_kInterpol1,
|
|---|
| 52 | //IDM_kCaosPrintRings,
|
|---|
| 53 | //IDM_kCaosPrintLeds,
|
|---|
| 54 | IDM_kCaosWriteStart,
|
|---|
| 55 | IDM_kCaosWriteStop,
|
|---|
| 56 | IDM_kResetHistograms,
|
|---|
| 57 | };
|
|---|
| 58 |
|
|---|
| 59 | Bool_t MStarguider::HandleTimer(TTimer *)
|
|---|
| 60 | {
|
|---|
| 61 | if (IsMapped())
|
|---|
| 62 | fImage->DoRedraw();
|
|---|
| 63 |
|
|---|
| 64 | return kTRUE;
|
|---|
| 65 | }
|
|---|
| 66 |
|
|---|
| 67 | #define kZOOM 96
|
|---|
| 68 | /*
|
|---|
| 69 | XY MStarguider::GetCoordinates() const
|
|---|
| 70 | {
|
|---|
| 71 | return fPZdAz->GetCoordinates();
|
|---|
| 72 | }
|
|---|
| 73 | */
|
|---|
| 74 | void MStarguider::InitGui(Int_t)
|
|---|
| 75 | {
|
|---|
| 76 | //fList = new MGList;
|
|---|
| 77 | //fList->SetOwner();
|
|---|
| 78 |
|
|---|
| 79 | /*
|
|---|
| 80 | const TGWindow *p=gClient->GetRoot();
|
|---|
| 81 |
|
|---|
| 82 | fChannel = new TGPopupMenu(p);
|
|---|
| 83 | fChannel->AddEntry("Starfield Camera", IDM_kChannel1);
|
|---|
| 84 | fChannel->AddEntry("TPoint Camera", IDM_kChannel2);
|
|---|
| 85 | //fChannel->AddEntry("Read from File", IDM_kChannel3);
|
|---|
| 86 | //if (channel<0)
|
|---|
| 87 | // fChannel->CheckEntry(IDM_kChannel3);
|
|---|
| 88 | //else
|
|---|
| 89 | fChannel->CheckEntry(channel==0?IDM_kChannel1:IDM_kChannel2);
|
|---|
| 90 | fChannel->Associate(this);
|
|---|
| 91 | fList->Add(fChannel);
|
|---|
| 92 |
|
|---|
| 93 | fFileType = new TGPopupMenu(p);
|
|---|
| 94 | fFileType->AddEntry("PP&M", IDM_kPPM);
|
|---|
| 95 | fFileType->AddEntry("&PNG", IDM_kPNG);
|
|---|
| 96 | fFileType->CheckEntry(IDM_kPNG);
|
|---|
| 97 | fFileType->Associate(this);
|
|---|
| 98 | fList->Add(fFileType);
|
|---|
| 99 |
|
|---|
| 100 | fWriteType = new TGPopupMenu(p);
|
|---|
| 101 | fWriteType->AddEntry("&Once", IDM_kOnce);
|
|---|
| 102 | fWriteType->AddEntry("&Continous", IDM_kContinous);
|
|---|
| 103 | fWriteType->CheckEntry(IDM_kOnce);
|
|---|
| 104 | fWriteType->Associate(this);
|
|---|
| 105 | fList->Add(fWriteType);
|
|---|
| 106 |
|
|---|
| 107 | fWriteRate = new TGPopupMenu(p);
|
|---|
| 108 | fWriteRate->AddEntry("25/s", IDM_kRate25ps);
|
|---|
| 109 | fWriteRate->AddEntry("5/s", IDM_kRate5ps);
|
|---|
| 110 | fWriteRate->AddEntry("1s", IDM_kRate1s);
|
|---|
| 111 | fWriteRate->AddEntry("5s", IDM_kRate5s);
|
|---|
| 112 | fWriteRate->AddEntry("30s", IDM_kRate30s);
|
|---|
| 113 | fWriteRate->AddEntry("1min", IDM_kRate1m);
|
|---|
| 114 | fWriteRate->AddEntry("5min", IDM_kRate5m);
|
|---|
| 115 | fWriteRate->CheckEntry(IDM_kRate1m);
|
|---|
| 116 | fWriteRate->Associate(this);
|
|---|
| 117 | fList->Add(fWriteRate);
|
|---|
| 118 |
|
|---|
| 119 | fWrtRate = 25*60;
|
|---|
| 120 |
|
|---|
| 121 | fLimMag = new TGPopupMenu(p);
|
|---|
| 122 | fLimMag->AddEntry("3", IDM_kLimMag3);
|
|---|
| 123 | fLimMag->AddEntry("4", IDM_kLimMag4);
|
|---|
| 124 | fLimMag->AddEntry("5", IDM_kLimMag5);
|
|---|
| 125 | fLimMag->AddEntry("6", IDM_kLimMag6);
|
|---|
| 126 | fLimMag->AddEntry("7", IDM_kLimMag7);
|
|---|
| 127 | fLimMag->AddEntry("8", IDM_kLimMag8);
|
|---|
| 128 | fLimMag->AddEntry("9", IDM_kLimMag9);
|
|---|
| 129 | fLimMag->CheckEntry(IDM_kLimMag9);
|
|---|
| 130 | fLimMag->Associate(this);
|
|---|
| 131 | fList->Add(fLimMag);
|
|---|
| 132 | */
|
|---|
| 133 | //fSao->SetLimitMag(9.0);
|
|---|
| 134 |
|
|---|
| 135 | const TGWindow *p=gClient->GetRoot();
|
|---|
| 136 |
|
|---|
| 137 | fInterpol = new TGPopupMenu(p);
|
|---|
| 138 | fInterpol->AddEntry("250", IDM_kInterpol250);
|
|---|
| 139 | fInterpol->AddEntry("125", IDM_kInterpol125);
|
|---|
| 140 | fInterpol->AddEntry("50", IDM_kInterpol50);
|
|---|
| 141 | fInterpol->AddEntry("25", IDM_kInterpol25);
|
|---|
| 142 | fInterpol->AddEntry("10", IDM_kInterpol10);
|
|---|
| 143 | fInterpol->AddEntry("5", IDM_kInterpol5);
|
|---|
| 144 | fInterpol->AddEntry("2", IDM_kInterpol2);
|
|---|
| 145 | fInterpol->AddEntry("Off", IDM_kInterpol1);
|
|---|
| 146 | fInterpol->Associate(this);
|
|---|
| 147 | //fList->Add(fInterpol);
|
|---|
| 148 |
|
|---|
| 149 | TString disp=gVirtualX->DisplayName();
|
|---|
| 150 | cout << "Display: " << disp << endl;
|
|---|
| 151 | if (disp.First(':')>=0)
|
|---|
| 152 | disp=disp(0, disp.First(':'));
|
|---|
| 153 |
|
|---|
| 154 | if (disp.IsNull() || disp==(TString)"localhost")
|
|---|
| 155 | {
|
|---|
| 156 | fInterpol->CheckEntry(IDM_kInterpol25);
|
|---|
| 157 | fIntRate = 25;
|
|---|
| 158 | }
|
|---|
| 159 | else
|
|---|
| 160 | {
|
|---|
| 161 | fInterpol->CheckEntry(IDM_kInterpol125);
|
|---|
| 162 | fIntRate = 125;
|
|---|
| 163 | }
|
|---|
| 164 |
|
|---|
| 165 | /*
|
|---|
| 166 | fCaosPrint = new TGPopupMenu(p);
|
|---|
| 167 | fCaosPrint->AddEntry("&Leds", IDM_kCaosPrintLeds);
|
|---|
| 168 | fCaosPrint->AddEntry("&Rings", IDM_kCaosPrintRings);
|
|---|
| 169 | fCaosPrint->Associate(this);
|
|---|
| 170 | fList->Add(fCaosPrint);
|
|---|
| 171 |
|
|---|
| 172 | fCaosWrite = new TGPopupMenu(p);
|
|---|
| 173 | fCaosWrite->AddEntry("&Start", IDM_kCaosWriteStart);
|
|---|
| 174 | fCaosWrite->AddEntry("Sto&p", IDM_kCaosWriteStop);
|
|---|
| 175 | fCaosWrite->DisableEntry(IDM_kCaosWriteStop);
|
|---|
| 176 | fCaosWrite->Associate(this);
|
|---|
| 177 | fList->Add(fCaosWrite);
|
|---|
| 178 |
|
|---|
| 179 | fCaosAnalyse = new TGPopupMenu(p);
|
|---|
| 180 | fCaosAnalyse->AddEntry("S&tart Analysis", IDM_kCaosAnalStart);
|
|---|
| 181 | fCaosAnalyse->AddEntry("St&op Analysis", IDM_kCaosAnalStop);
|
|---|
| 182 | fCaosAnalyse->DisableEntry(IDM_kCaosAnalStop);
|
|---|
| 183 | fCaosAnalyse->Associate(this);
|
|---|
| 184 | fList->Add(fCaosAnalyse);
|
|---|
| 185 | */
|
|---|
| 186 | fMenu = new TGMenuBar(this, 0, 0, kHorizontalFrame);
|
|---|
| 187 | fDisplay = fMenu->AddPopup("&Display");
|
|---|
| 188 | //fMode = fMenu->AddPopup("&Mode");
|
|---|
| 189 | //fWritePictures = fMenu->AddPopup("&WritePics");
|
|---|
| 190 | fSetup = fMenu->AddPopup("&Setup");
|
|---|
| 191 | //fOperations = fMenu->AddPopup("&Operations");
|
|---|
| 192 | fMenu->Resize(fMenu->GetDefaultSize());
|
|---|
| 193 | AddFrame(fMenu);
|
|---|
| 194 |
|
|---|
| 195 | //
|
|---|
| 196 | // Create Menu for MStarguider Display
|
|---|
| 197 | //
|
|---|
| 198 | //fDisplay = new MMGPopupMenu(p);
|
|---|
| 199 | fDisplay->AddEntry("&Filter", IDM_kFilter);
|
|---|
| 200 | fDisplay->AddEntry("Stretch", IDM_kStretch);
|
|---|
| 201 | fDisplay->AddSeparator();
|
|---|
| 202 | fDisplay->AddEntry("Find &Star", IDM_kFindStar);
|
|---|
| 203 | fDisplay->AddEntry("C&aos Filter", IDM_kCaosFilter);
|
|---|
| 204 | //fDisplay->AddSeparator();
|
|---|
| 205 | //if (channel>=0)
|
|---|
| 206 | // fDisplay->AddPopup("&Input", fChannel);
|
|---|
| 207 | // fDisplay->CheckEntry(IDM_kStretch);
|
|---|
| 208 | fDisplay->CheckEntry(IDM_kFindStar);
|
|---|
| 209 | fDisplay->CheckEntry(IDM_kCaosFilter);
|
|---|
| 210 | fDisplay->Associate(this);
|
|---|
| 211 |
|
|---|
| 212 | //fMode->AddEntry("Tpoint", IDM_kTpointMode);
|
|---|
| 213 | //fMode->Associate(this);
|
|---|
| 214 | /*
|
|---|
| 215 | fWritePictures->AddEntry("&Start", IDM_kStart);
|
|---|
| 216 | fWritePictures->AddEntry("Sto&p", IDM_kStop);
|
|---|
| 217 | fWritePictures->AddSeparator();
|
|---|
| 218 | //fWritePictures->AddPopup("File &Type", fFileType);
|
|---|
| 219 | fWritePictures->AddPopup("&Write Type", fWriteType);
|
|---|
| 220 | fWritePictures->AddPopup("Write &Rate", fWriteRate);
|
|---|
| 221 | fWritePictures->DisableEntry(IDM_kStop);
|
|---|
| 222 | fWritePictures->Associate(this);
|
|---|
| 223 | */
|
|---|
| 224 | fSetup->AddPopup("Disp. &Interpolation", fInterpol);
|
|---|
| 225 | fSetup->Associate(this);
|
|---|
| 226 |
|
|---|
| 227 | /*
|
|---|
| 228 | fCaOs = new TGPopupMenu(p);
|
|---|
| 229 | //fCaOs->AddPopup("&Write", fCaosWrite);
|
|---|
| 230 | fCaOs->AddPopup("&Print", fCaosPrint);
|
|---|
| 231 | //fCaOs->AddPopup("&Analyse", fCaosAnalyse);
|
|---|
| 232 | fCaOs->Associate(this);
|
|---|
| 233 | fList->Add(fCaOs);
|
|---|
| 234 | */
|
|---|
| 235 | //RA,Dec for catalog
|
|---|
| 236 | /*
|
|---|
| 237 | fCRaDec = new MGCoordinates(this, kETypeRaDec);
|
|---|
| 238 | fCRaDec->Move(4, fMenu->GetDefaultHeight()+584);
|
|---|
| 239 | AddFrame(fCRaDec);
|
|---|
| 240 |
|
|---|
| 241 | //telescope position
|
|---|
| 242 | fCZdAz = new MGCoordinates(this, kETypeZdAz, 2);
|
|---|
| 243 | fCZdAz->Move(240+12+28, fMenu->GetDefaultHeight()+597);
|
|---|
| 244 | AddFrame(fCZdAz);
|
|---|
| 245 |
|
|---|
| 246 | //starguider position
|
|---|
| 247 | fPZdAz = new MGCoordinates(this, kETypeZdAz, 2);
|
|---|
| 248 | fPZdAz->Move(240+12+28, fMenu->GetDefaultHeight()+640);
|
|---|
| 249 | AddFrame(fPZdAz);
|
|---|
| 250 |
|
|---|
| 251 | //mispointing
|
|---|
| 252 | fDZdAz = new MGCoordinates(this, kETypeZdAz, 2);
|
|---|
| 253 | fDZdAz->Move(240+12+28, fMenu->GetDefaultHeight()+683);
|
|---|
| 254 | AddFrame(fDZdAz);
|
|---|
| 255 |
|
|---|
| 256 | fSZdAz = new MGCoordinates(this, kETypeZdAz, 2);
|
|---|
| 257 | fSZdAz->Move(240+12+28, fMenu->GetDefaultHeight()+795);
|
|---|
| 258 | AddFrame(fSZdAz);
|
|---|
| 259 |
|
|---|
| 260 | fGNumStars = new MGNumStars(this, 235);
|
|---|
| 261 | fGNumStars->DrawText("Number of stars");
|
|---|
| 262 | fGNumStars->Move(278, fMenu->GetDefaultHeight()+713);
|
|---|
| 263 | fList->Add(fGNumStars);
|
|---|
| 264 |
|
|---|
| 265 | fTPoint = new TGTextButton(this, "TPoint");
|
|---|
| 266 | //fTPoint->Move(4, fMenu->GetDefaultHeight()+785);
|
|---|
| 267 | fTPoint->Move(170, fMenu->GetDefaultHeight()+785);
|
|---|
| 268 | fTPoint->AllowStayDown(kTRUE);
|
|---|
| 269 | AddFrame(fTPoint);
|
|---|
| 270 |
|
|---|
| 271 | fStargTPoint = new TGTextButton(this, "StargTPoint");
|
|---|
| 272 | fStargTPoint->Move(170, fMenu->GetDefaultHeight()+785);
|
|---|
| 273 | fStargTPoint->AllowStayDown(kTRUE);
|
|---|
| 274 | AddFrame(fStargTPoint);
|
|---|
| 275 |
|
|---|
| 276 | fFps = new TGLabel(this, "---fps");
|
|---|
| 277 | fFps->SetTextJustify(kTextRight);
|
|---|
| 278 | fFps->Move(650-495, fMenu->GetDefaultHeight()+714+23);
|
|---|
| 279 | AddFrame(fFps);
|
|---|
| 280 |
|
|---|
| 281 | fPosZoom = new TGLabel(this, "(----, ----) ----.--d/----.--d");
|
|---|
| 282 | fPosZoom->SetTextJustify(kTextLeft);
|
|---|
| 283 | fPosZoom->Move(4, fMenu->GetDefaultHeight()+765);
|
|---|
| 284 | AddFrame(fPosZoom);
|
|---|
| 285 |
|
|---|
| 286 | fSkyBright = new TGLabel(this, "Sky Brightness: --- ");
|
|---|
| 287 | fSkyBright->SetTextJustify(kTextLeft);
|
|---|
| 288 | fSkyBright->Move(4, fMenu->GetDefaultHeight()+785);
|
|---|
| 289 | AddFrame(fSkyBright);
|
|---|
| 290 |
|
|---|
| 291 | TGLabel *l = new TGLabel(this, "deg");
|
|---|
| 292 | l->SetTextJustify(kTextLeft);
|
|---|
| 293 | l->Move(606-412, fMenu->GetDefaultHeight()+669);
|
|---|
| 294 | AddFrame(l);
|
|---|
| 295 |
|
|---|
| 296 | l = new TGLabel(this, "arcsec/pix");
|
|---|
| 297 | l->SetTextJustify(kTextLeft);
|
|---|
| 298 | l->Move(606-412, fMenu->GetDefaultHeight()+692);
|
|---|
| 299 | AddFrame(l);
|
|---|
| 300 |
|
|---|
| 301 | l = new TGLabel(this, "sigma");
|
|---|
| 302 | l->SetTextJustify(kTextLeft);
|
|---|
| 303 | l->Move(606-412, fMenu->GetDefaultHeight()+715);
|
|---|
| 304 | AddFrame(l);
|
|---|
| 305 |
|
|---|
| 306 | fCZdAzText = new TGLabel(this, "Zd/Az telescope pointing at");
|
|---|
| 307 | fCZdAzText->SetTextJustify(kTextLeft);
|
|---|
| 308 | fCZdAzText->Move(240+12+20+7, fMenu->GetDefaultHeight()+584-5);
|
|---|
| 309 | AddFrame(fCZdAzText);
|
|---|
| 310 |
|
|---|
| 311 | fPZdAzText = new TGLabel(this, "Zd/Az starguider pointing at");
|
|---|
| 312 | fPZdAzText->SetTextJustify(kTextLeft);
|
|---|
| 313 | fPZdAzText->Move(240+12+20+7, fMenu->GetDefaultHeight()+630+20-5-23);
|
|---|
| 314 | AddFrame(fPZdAzText);
|
|---|
| 315 |
|
|---|
| 316 | fDZdAzText = new TGLabel(this, "Zd/Az mispointing");
|
|---|
| 317 | fDZdAzText->SetTextJustify(kTextLeft);
|
|---|
| 318 | fDZdAzText->Move(240+12+20+7, fMenu->GetDefaultHeight()+676+2*20-5-46);
|
|---|
| 319 | AddFrame(fDZdAzText);
|
|---|
| 320 | */
|
|---|
| 321 | // Set input box for rotation angle
|
|---|
| 322 | /*
|
|---|
| 323 | fAngle = new TGTextEntry(this, " ", IDM_kAngle);
|
|---|
| 324 | fAngle->SetAlignment(kTextCenterX);
|
|---|
| 325 | fAngle->Move(547-410, fMenu->GetDefaultHeight()+667);
|
|---|
| 326 | AddFrame(fAngle);
|
|---|
| 327 |
|
|---|
| 328 | //SetRotationAngle(-0.2);
|
|---|
| 329 |
|
|---|
| 330 | // Set input box for pixel size
|
|---|
| 331 | fPixSize = new TGTextEntry(this, " ", IDM_kPixSize);
|
|---|
| 332 | fPixSize->SetAlignment(kTextCenterX);
|
|---|
| 333 | fPixSize->Move(547-410, fMenu->GetDefaultHeight()+690);
|
|---|
| 334 | AddFrame(fPixSize);
|
|---|
| 335 | */
|
|---|
| 336 | //SetPixSize(48.9);
|
|---|
| 337 |
|
|---|
| 338 | // Set input box for cleaning cut
|
|---|
| 339 | //fCut = new TGTextEntry(this, " ", IDM_kCut);
|
|---|
| 340 | //fCut->SetAlignment(kTextCenterX);
|
|---|
| 341 | //fCut->Move(547-410, fMenu->GetDefaultHeight()+713);
|
|---|
| 342 | //AddFrame(fCut);
|
|---|
| 343 |
|
|---|
| 344 | //SetCut(3.0);
|
|---|
| 345 |
|
|---|
| 346 | // TGHorizontal3DLine *fLineSep = new TGHorizontal3DLine(this);
|
|---|
| 347 | // AddFrame(fLineSep, new TGLayoutHints (kLHintsNormal | kLHintsExpandX));
|
|---|
| 348 | // fList->Add(fLineSep);
|
|---|
| 349 |
|
|---|
| 350 | //
|
|---|
| 351 | // Create Image Display
|
|---|
| 352 | /*
|
|---|
| 353 | fZoomImage = new MGImage(this, kZOOM, kZOOM);
|
|---|
| 354 | // fZoomImage->Move(768-kZOOM-2, 700-kZOOM-2);
|
|---|
| 355 | fZoomImage->Move(4, 700-kZOOM-2+85);
|
|---|
| 356 | AddFrame(fZoomImage);
|
|---|
| 357 | */
|
|---|
| 358 | fImage = new MGImage(this, 768, 576);
|
|---|
| 359 | fImage->Move(0, fMenu->GetDefaultHeight());
|
|---|
| 360 | AddFrame(fImage);
|
|---|
| 361 |
|
|---|
| 362 | const Int_t w = 768;
|
|---|
| 363 | const Int_t h = 576;
|
|---|
| 364 | SetWMSizeHints(w, h, w, h, 1, 1); // set the smallest and biggest size of the Main frame
|
|---|
| 365 |
|
|---|
| 366 | //
|
|---|
| 367 | // Make everything visible
|
|---|
| 368 | //
|
|---|
| 369 | SetWindowName("TPoint Main Window");
|
|---|
| 370 | SetIconName("TPoint");
|
|---|
| 371 |
|
|---|
| 372 | MapSubwindows();
|
|---|
| 373 | //fTPoint->UnmapWindow();
|
|---|
| 374 | //fStargTPoint->UnmapWindow();
|
|---|
| 375 | //fGStarg->UnmapWindow();
|
|---|
| 376 | //fGNumStars->UnmapWindow();
|
|---|
| 377 | //fCRaDec->UnmapWindow();
|
|---|
| 378 | //fCZdAz->UnmapWindow();
|
|---|
| 379 | //fCZdAzText->UnmapWindow();
|
|---|
| 380 | //fPZdAz->UnmapWindow();
|
|---|
| 381 | //fPZdAzText->UnmapWindow();
|
|---|
| 382 | //fDZdAz->UnmapWindow();
|
|---|
| 383 | //fDZdAzText->UnmapWindow();
|
|---|
| 384 | //fSZdAz->UnmapWindow();
|
|---|
| 385 | //fSkyBright->UnmapWindow();
|
|---|
| 386 | MapWindow();
|
|---|
| 387 |
|
|---|
| 388 |
|
|---|
| 389 | //IconifyWindow();
|
|---|
| 390 |
|
|---|
| 391 | //------------------------------------------------------------
|
|---|
| 392 | // XY xy(3.819444, 24.05333);
|
|---|
| 393 | // fCRaDec->SetCoordinates(xy);
|
|---|
| 394 | // fRaDec->Set(xy.X()*360/24, xy.Y());
|
|---|
| 395 | //------------------------------------------------------------
|
|---|
| 396 | }
|
|---|
| 397 |
|
|---|
| 398 | MStarguider::MStarguider(Int_t channel) : TGMainFrame(gClient->GetRoot(), 768, 840),
|
|---|
| 399 | fDimData("TPOINT/DATA", "D:11", (void*)NULL, 0),
|
|---|
| 400 | fDimTPoint("TPOINT/EXECUTE", "", this),
|
|---|
| 401 | fDimScreenshot("TPOINT/SCREENSHOT", "B:1;C", this),
|
|---|
| 402 | fRadius(200), fFindStarCut(2.), fFindStarBox(30), fTPointMode(0)
|
|---|
| 403 | {
|
|---|
| 404 | cout << " #### FIXME: Make MCaos Thread safe!" << endl;
|
|---|
| 405 |
|
|---|
| 406 | // This means that all objects added with AddFrame are deleted
|
|---|
| 407 | // automatically, including all LayoutHints.
|
|---|
| 408 | SetCleanup();
|
|---|
| 409 |
|
|---|
| 410 | fCaos = new MCaos;
|
|---|
| 411 | fCaos->ReadResources("leds_fact.txt");
|
|---|
| 412 |
|
|---|
| 413 | InitGui(channel);
|
|---|
| 414 |
|
|---|
| 415 | fTimer=new TTimer(this, 1000/25); // 40ms
|
|---|
| 416 | fTimer->TurnOn();
|
|---|
| 417 |
|
|---|
| 418 | gVirtualX->GrabButton(fId, kButton2, 0, 0, 0, 0, kTRUE);
|
|---|
| 419 |
|
|---|
| 420 | fGetter = new Camera(*this, channel);
|
|---|
| 421 |
|
|---|
| 422 | DimClient::setNoDataCopy();
|
|---|
| 423 | DimServer::start("TPOINT");
|
|---|
| 424 | }
|
|---|
| 425 |
|
|---|
| 426 | MStarguider::~MStarguider()
|
|---|
| 427 | {
|
|---|
| 428 | DimServer::stop();
|
|---|
| 429 |
|
|---|
| 430 | fGetter->ExitLoop();
|
|---|
| 431 | delete fGetter;
|
|---|
| 432 |
|
|---|
| 433 | gVirtualX->GrabButton(fId, kButton2, 0, 0, 0, 0, kFALSE);
|
|---|
| 434 |
|
|---|
| 435 | fTimer->TurnOff();
|
|---|
| 436 | delete fTimer;
|
|---|
| 437 |
|
|---|
| 438 | delete fInterpol;
|
|---|
| 439 |
|
|---|
| 440 | //delete fList;
|
|---|
| 441 | delete fCaos;
|
|---|
| 442 |
|
|---|
| 443 | cout << "Camera Display destroyed." << endl;
|
|---|
| 444 | }
|
|---|
| 445 |
|
|---|
| 446 | /*
|
|---|
| 447 | void MStarguider::SetupEnv(TEnv &env)
|
|---|
| 448 | {
|
|---|
| 449 | fCaos->ReadEnv(env, "TPointLeds", kTRUE);
|
|---|
| 450 | //fStargCaos->ReadEnv(env, "StarguiderLeds", kTRUE);
|
|---|
| 451 |
|
|---|
| 452 | //SetRotationAngle(env.GetValue("Starguider.RotationAngle", fSao->GetRotationAngle()));
|
|---|
| 453 | //SetCut(env.GetValue("Starguider.CleaningLevel", atof(fCut->GetText())));
|
|---|
| 454 |
|
|---|
| 455 | //SetPixSize(env.GetValue("StarguiderLeds.ArcsecPerPixel", fSao->GetPixSize()));
|
|---|
| 456 |
|
|---|
| 457 | fRadius = env.GetValue("Leds.Radius", fRadius);
|
|---|
| 458 |
|
|---|
| 459 | fStarguiderW = env.GetValue("Starguider.Width", fStarguiderW);
|
|---|
| 460 | fStarguiderH = env.GetValue("Starguider.Height", fStarguiderH);
|
|---|
| 461 | fStarguiderX = env.GetValue("Starguider.X", fStarguiderX);
|
|---|
| 462 | fStarguiderY = env.GetValue("Starguider.Y", fStarguiderY);
|
|---|
| 463 |
|
|---|
| 464 | fSkyOffsetX = env.GetValue("Starguider.SkyOffsetX", fSkyOffsetX);
|
|---|
| 465 | fSkyOffsetY = env.GetValue("Starguider.SkyOffsetY", fSkyOffsetY);
|
|---|
| 466 |
|
|---|
| 467 | fFindStarBox = env.GetValue("FindStar.SizeBox", fFindStarBox);
|
|---|
| 468 | fFindStarCut = env.GetValue("FindStar.CleaningLevel", fFindStarCut);
|
|---|
| 469 | }
|
|---|
| 470 | */
|
|---|
| 471 |
|
|---|
| 472 | void MStarguider::Layout()
|
|---|
| 473 | {
|
|---|
| 474 | // Resize(GetDefaultSize());
|
|---|
| 475 | }
|
|---|
| 476 |
|
|---|
| 477 | void MStarguider::CloseWindow()
|
|---|
| 478 | {
|
|---|
| 479 | cout << "EventDisplay::CloseWindow: Exit Application Loop." << endl;
|
|---|
| 480 |
|
|---|
| 481 | //fClient.ExitLoop();
|
|---|
| 482 | // cout << "FIXME: ExitLoop not called!!!!!!" << endl;
|
|---|
| 483 | fGetter->ExitLoop();
|
|---|
| 484 | gSystem->ExitLoop();
|
|---|
| 485 | }
|
|---|
| 486 |
|
|---|
| 487 | void MStarguider::SwitchOff(TGPopupMenu *p, UInt_t id)
|
|---|
| 488 | {
|
|---|
| 489 | p->UnCheckEntry(id);
|
|---|
| 490 | p->DisableEntry(id);
|
|---|
| 491 | }
|
|---|
| 492 |
|
|---|
| 493 | /*
|
|---|
| 494 | void MStarguider::SetChannel()
|
|---|
| 495 | {
|
|---|
| 496 | if (fChannel->IsEntryChecked(IDM_kChannel3))
|
|---|
| 497 | {
|
|---|
| 498 | if (dynamic_cast<PngReader*>(fGetter)==0)
|
|---|
| 499 | {
|
|---|
| 500 | delete fGetter;
|
|---|
| 501 | fGetter=new PngReader(*this);
|
|---|
| 502 | }
|
|---|
| 503 | }
|
|---|
| 504 | else
|
|---|
| 505 | {
|
|---|
| 506 | const Int_t ch = fChannel->IsEntryChecked(IDM_kChannel1) ? 0 : 1;
|
|---|
| 507 | if (dynamic_cast<Camera*>(fGetter)==0)
|
|---|
| 508 | {
|
|---|
| 509 | delete fGetter;
|
|---|
| 510 | fGetter = new Camera(*this, ch);
|
|---|
| 511 | }
|
|---|
| 512 | else
|
|---|
| 513 | fGetter->SetChannel(ch);
|
|---|
| 514 | }
|
|---|
| 515 | }*/
|
|---|
| 516 |
|
|---|
| 517 | void MStarguider::Toggle(TGPopupMenu *p, UInt_t id)
|
|---|
| 518 | {
|
|---|
| 519 | if (p->IsEntryChecked(id))
|
|---|
| 520 | p->UnCheckEntry(id);
|
|---|
| 521 | else
|
|---|
| 522 | p->CheckEntry(id);
|
|---|
| 523 | }
|
|---|
| 524 |
|
|---|
| 525 |
|
|---|
| 526 | Bool_t MStarguider::ProcessMessage(Long_t msg, Long_t mp1, Long_t)
|
|---|
| 527 | {
|
|---|
| 528 | switch (GET_MSG(msg))
|
|---|
| 529 | {
|
|---|
| 530 | case kC_TEXTENTRY:
|
|---|
| 531 | if (GET_SUBMSG(msg)==kTE_ENTER)
|
|---|
| 532 | switch (mp1)
|
|---|
| 533 | {
|
|---|
| 534 | /*
|
|---|
| 535 | case IDM_kPixSize:
|
|---|
| 536 | {
|
|---|
| 537 | const Float_t pixsize = atof(fPixSize->GetText());
|
|---|
| 538 | gLog << all << "Pixel Size changed to " << pixsize << "\"/pix" << endl;
|
|---|
| 539 | fSao->SetPixSize(pixsize);
|
|---|
| 540 | return kTRUE;
|
|---|
| 541 | }
|
|---|
| 542 | case IDM_kAngle:
|
|---|
| 543 | {
|
|---|
| 544 | const Float_t angle = atof(fAngle->GetText());
|
|---|
| 545 | gLog << all << "Rotation Angle changed to " << angle << "deg" << endl;
|
|---|
| 546 | fSao->SetRotationAngle(angle);
|
|---|
| 547 | return kTRUE;
|
|---|
| 548 | }
|
|---|
| 549 | case IDM_kCut:
|
|---|
| 550 | {
|
|---|
| 551 | const Float_t cut = atof(fCut->GetText());
|
|---|
| 552 | gLog << all << "Starguider cleaning level changed to " << cut << " sigma." << endl;
|
|---|
| 553 | return kTRUE;
|
|---|
| 554 | }*/
|
|---|
| 555 | }
|
|---|
| 556 | return kTRUE;
|
|---|
| 557 |
|
|---|
| 558 | case kC_COMMAND:
|
|---|
| 559 | switch (GET_SUBMSG(msg))
|
|---|
| 560 | {
|
|---|
| 561 | case kCM_MENU:
|
|---|
| 562 | switch (mp1)
|
|---|
| 563 | {
|
|---|
| 564 | /*
|
|---|
| 565 | case IDM_kTpointMode:
|
|---|
| 566 | Toggle(fMode, IDM_kTpointMode);
|
|---|
| 567 |
|
|---|
| 568 | if (fMode->IsEntryChecked(IDM_kTpointMode))
|
|---|
| 569 | {
|
|---|
| 570 | //unchecking not needed items
|
|---|
| 571 | //general
|
|---|
| 572 | SwitchOff(fDisplay, IDM_kFilter);
|
|---|
| 573 | SwitchOff(fChannel, IDM_kChannel3);
|
|---|
| 574 |
|
|---|
| 575 | //from starguider
|
|---|
| 576 | //SwitchOff(fDisplay, IDM_kStargCaosFilter);
|
|---|
| 577 | //SwitchOff(fDisplay, IDM_kCatalog);
|
|---|
| 578 | //SwitchOff(fDisplay, IDM_kStarguider);
|
|---|
| 579 | //ToggleStarguider();
|
|---|
| 580 | //fMode->UnCheckEntry(IDM_kStarguiderMode);
|
|---|
| 581 | //SwitchOff(fOperations, IDM_kStargAnalysis);
|
|---|
| 582 | //ToggleStargAnalysis();
|
|---|
| 583 |
|
|---|
| 584 | //switch camera
|
|---|
| 585 | SwitchOff(fChannel, IDM_kChannel1);
|
|---|
| 586 | fChannel->CheckEntry(IDM_kChannel2);
|
|---|
| 587 |
|
|---|
| 588 | SetChannel();
|
|---|
| 589 |
|
|---|
| 590 | //checking needed items
|
|---|
| 591 | fDisplay->UnCheckEntry(IDM_kStretch);
|
|---|
| 592 | fDisplay->CheckEntry(IDM_kCaosFilter);
|
|---|
| 593 | ToggleCaosFilter();
|
|---|
| 594 | fDisplay->CheckEntry(IDM_kFindStar);
|
|---|
| 595 | fTPoint->MapWindow();
|
|---|
| 596 | }
|
|---|
| 597 | else
|
|---|
| 598 | {
|
|---|
| 599 | //enable
|
|---|
| 600 | //starguider items
|
|---|
| 601 | //fDisplay->EnableEntry(IDM_kStargCaosFilter);
|
|---|
| 602 | //fDisplay->EnableEntry(IDM_kCatalog);
|
|---|
| 603 | //fDisplay->EnableEntry(IDM_kStarguider);
|
|---|
| 604 | //fOperations->EnableEntry(IDM_kStargAnalysis);
|
|---|
| 605 |
|
|---|
| 606 | //general
|
|---|
| 607 | fDisplay->EnableEntry(IDM_kFilter);
|
|---|
| 608 | fChannel->EnableEntry(IDM_kChannel1);
|
|---|
| 609 | fChannel->EnableEntry(IDM_kChannel3);
|
|---|
| 610 |
|
|---|
| 611 | //tpoint
|
|---|
| 612 | fDisplay->UnCheckEntry(IDM_kCaosFilter);
|
|---|
| 613 | ToggleCaosFilter();
|
|---|
| 614 | fDisplay->UnCheckEntry(IDM_kFindStar);
|
|---|
| 615 | fTPoint->UnmapWindow();
|
|---|
| 616 | }
|
|---|
| 617 | return kTRUE;
|
|---|
| 618 | */
|
|---|
| 619 | case IDM_kFilter:
|
|---|
| 620 | Toggle(fDisplay, IDM_kFilter);
|
|---|
| 621 | return kTRUE;
|
|---|
| 622 |
|
|---|
| 623 | case IDM_kFindStar:
|
|---|
| 624 | Toggle(fDisplay, IDM_kFindStar);
|
|---|
| 625 | //ToggleFindStar();
|
|---|
| 626 | return kTRUE;
|
|---|
| 627 |
|
|---|
| 628 | case IDM_kStretch:
|
|---|
| 629 | Toggle(fDisplay, IDM_kStretch);
|
|---|
| 630 | return kTRUE;
|
|---|
| 631 |
|
|---|
| 632 | case IDM_kCaosFilter:
|
|---|
| 633 | Toggle(fDisplay, IDM_kCaosFilter);
|
|---|
| 634 | //ToggleCaosFilter();
|
|---|
| 635 | return kTRUE;
|
|---|
| 636 | /*
|
|---|
| 637 | case IDM_kCaosPrintLeds:
|
|---|
| 638 | case IDM_kCaosPrintRings:
|
|---|
| 639 | Toggle(fCaosPrint, mp1);
|
|---|
| 640 | return kTRUE;
|
|---|
| 641 |
|
|---|
| 642 | case IDM_kCaosAnalStart:
|
|---|
| 643 | fCaosAnalyse->DisableEntry(IDM_kCaosAnalStart);
|
|---|
| 644 | fCaosAnalyse->EnableEntry(IDM_kCaosAnalStop);
|
|---|
| 645 | //fCaos->InitHistograms();
|
|---|
| 646 | return kTRUE;
|
|---|
| 647 |
|
|---|
| 648 | case IDM_kCaosAnalStop:
|
|---|
| 649 | fCaosAnalyse->DisableEntry(IDM_kCaosAnalStop);
|
|---|
| 650 | fCaosAnalyse->EnableEntry(IDM_kCaosAnalStart);
|
|---|
| 651 | fCaos->ShowHistograms();
|
|---|
| 652 | fCaos->DeleteHistograms();
|
|---|
| 653 | return kTRUE;
|
|---|
| 654 |
|
|---|
| 655 | case IDM_kCaosWriteStart:
|
|---|
| 656 | fCaosWrite->DisableEntry(IDM_kCaosWriteStart);
|
|---|
| 657 | fCaosWrite->EnableEntry(IDM_kCaosWriteStop);
|
|---|
| 658 | fCaos->OpenFile();
|
|---|
| 659 | return kTRUE;
|
|---|
| 660 |
|
|---|
| 661 | case IDM_kCaosWriteStop:
|
|---|
| 662 | fCaosWrite->DisableEntry(IDM_kCaosWriteStop);
|
|---|
| 663 | fCaosWrite->EnableEntry(IDM_kCaosWriteStart);
|
|---|
| 664 | fCaos->CloseFile();
|
|---|
| 665 | return kTRUE;
|
|---|
| 666 |
|
|---|
| 667 | case IDM_kStart:
|
|---|
| 668 | fWritePictures->DisableEntry(IDM_kStart);
|
|---|
| 669 | fWritePictures->EnableEntry(IDM_kStop);
|
|---|
| 670 | return kTRUE;
|
|---|
| 671 |
|
|---|
| 672 | case IDM_kStop:
|
|---|
| 673 | fWritePictures->DisableEntry(IDM_kStop);
|
|---|
| 674 | fWritePictures->EnableEntry(IDM_kStart);
|
|---|
| 675 | return kTRUE;
|
|---|
| 676 |
|
|---|
| 677 | case IDM_kPNG:
|
|---|
| 678 | fFileType->CheckEntry(IDM_kPNG);
|
|---|
| 679 | fFileType->UnCheckEntry(IDM_kPPM);
|
|---|
| 680 | return kTRUE;
|
|---|
| 681 |
|
|---|
| 682 | case IDM_kPPM:
|
|---|
| 683 | fFileType->CheckEntry(IDM_kPPM);
|
|---|
| 684 | fFileType->UnCheckEntry(IDM_kPNG);
|
|---|
| 685 | return kTRUE;
|
|---|
| 686 |
|
|---|
| 687 | case IDM_kOnce:
|
|---|
| 688 | fWriteType->CheckEntry(IDM_kOnce);
|
|---|
| 689 | fWriteType->UnCheckEntry(IDM_kContinous);
|
|---|
| 690 | return kTRUE;
|
|---|
| 691 |
|
|---|
| 692 | case IDM_kContinous:
|
|---|
| 693 | fWriteType->CheckEntry(IDM_kContinous);
|
|---|
| 694 | fWriteType->UnCheckEntry(IDM_kOnce);
|
|---|
| 695 | return kTRUE;
|
|---|
| 696 |
|
|---|
| 697 | case IDM_kRate25ps:
|
|---|
| 698 | case IDM_kRate5ps:
|
|---|
| 699 | case IDM_kRate1s:
|
|---|
| 700 | case IDM_kRate5s:
|
|---|
| 701 | case IDM_kRate30s:
|
|---|
| 702 | case IDM_kRate1m:
|
|---|
| 703 | case IDM_kRate5m:
|
|---|
| 704 | for (int i=IDM_kRate25ps; i<=IDM_kRate5m; i++)
|
|---|
| 705 | if (mp1==i)
|
|---|
| 706 | fWriteRate->CheckEntry(i);
|
|---|
| 707 | else
|
|---|
| 708 | fWriteRate->UnCheckEntry(i);
|
|---|
| 709 | switch (mp1)
|
|---|
| 710 | {
|
|---|
| 711 | case IDM_kRate25ps:
|
|---|
| 712 | fWrtRate = 1;
|
|---|
| 713 | return kTRUE;
|
|---|
| 714 | case IDM_kRate5ps:
|
|---|
| 715 | fWrtRate = 5;
|
|---|
| 716 | return kTRUE;
|
|---|
| 717 | case IDM_kRate1s:
|
|---|
| 718 | fWrtRate = 25;
|
|---|
| 719 | return kTRUE;
|
|---|
| 720 | case IDM_kRate5s:
|
|---|
| 721 | fWrtRate = 5*25;
|
|---|
| 722 | return kTRUE;
|
|---|
| 723 | case IDM_kRate30s:
|
|---|
| 724 | fWrtRate = 30*25;
|
|---|
| 725 | return kTRUE;
|
|---|
| 726 | case IDM_kRate1m:
|
|---|
| 727 | fWrtRate = 60*25;
|
|---|
| 728 | return kTRUE;
|
|---|
| 729 | case IDM_kRate5m:
|
|---|
| 730 | fWrtRate = 5*60*25;
|
|---|
| 731 | return kTRUE;
|
|---|
| 732 | }
|
|---|
| 733 | return kTRUE;
|
|---|
| 734 |
|
|---|
| 735 | case IDM_kChannel1:
|
|---|
| 736 | case IDM_kChannel2:
|
|---|
| 737 | {
|
|---|
| 738 | const Int_t ch0 = fChannel->IsEntryChecked(IDM_kChannel1) ? 0 : 1;
|
|---|
| 739 | const Int_t ch1 = mp1==IDM_kChannel1 ? 0 : 1;
|
|---|
| 740 |
|
|---|
| 741 | if (ch0==ch1)
|
|---|
| 742 | return kTRUE;
|
|---|
| 743 |
|
|---|
| 744 | fChannel->CheckEntry (ch1==0?IDM_kChannel1:IDM_kChannel2);
|
|---|
| 745 | fChannel->UnCheckEntry(ch1==1?IDM_kChannel1:IDM_kChannel2);
|
|---|
| 746 |
|
|---|
| 747 | SetChannel();
|
|---|
| 748 | }
|
|---|
| 749 | return kTRUE;
|
|---|
| 750 | */
|
|---|
| 751 | case IDM_kInterpol250:
|
|---|
| 752 | case IDM_kInterpol125:
|
|---|
| 753 | case IDM_kInterpol50:
|
|---|
| 754 | case IDM_kInterpol25:
|
|---|
| 755 | case IDM_kInterpol10:
|
|---|
| 756 | case IDM_kInterpol5:
|
|---|
| 757 | case IDM_kInterpol2:
|
|---|
| 758 | case IDM_kInterpol1:
|
|---|
| 759 | for (int i=IDM_kInterpol250; i<=IDM_kInterpol1; i++)
|
|---|
| 760 | if (mp1==i)
|
|---|
| 761 | fInterpol->CheckEntry(i);
|
|---|
| 762 | else
|
|---|
| 763 | fInterpol->UnCheckEntry(i);
|
|---|
| 764 | switch (mp1)
|
|---|
| 765 | {
|
|---|
| 766 | case IDM_kInterpol1:
|
|---|
| 767 | fIntRate = 1;
|
|---|
| 768 | return kTRUE;
|
|---|
| 769 | case IDM_kInterpol2:
|
|---|
| 770 | fIntRate = 2;
|
|---|
| 771 | return kTRUE;
|
|---|
| 772 | case IDM_kInterpol5:
|
|---|
| 773 | fIntRate = 5;
|
|---|
| 774 | return kTRUE;
|
|---|
| 775 | case IDM_kInterpol10:
|
|---|
| 776 | fIntRate = 10;
|
|---|
| 777 | return kTRUE;
|
|---|
| 778 | case IDM_kInterpol25:
|
|---|
| 779 | fIntRate = 25;
|
|---|
| 780 | return kTRUE;
|
|---|
| 781 | case IDM_kInterpol50:
|
|---|
| 782 | fIntRate = 50;
|
|---|
| 783 | return kTRUE;
|
|---|
| 784 | case IDM_kInterpol125:
|
|---|
| 785 | fIntRate = 125;
|
|---|
| 786 | return kTRUE;
|
|---|
| 787 | case IDM_kInterpol250:
|
|---|
| 788 | fIntRate = 250;
|
|---|
| 789 | return kTRUE;
|
|---|
| 790 | }
|
|---|
| 791 | return kTRUE;
|
|---|
| 792 | }
|
|---|
| 793 | break;
|
|---|
| 794 | }
|
|---|
| 795 | break;
|
|---|
| 796 | }
|
|---|
| 797 |
|
|---|
| 798 | return kTRUE;
|
|---|
| 799 | }
|
|---|
| 800 | /*
|
|---|
| 801 | void MStarguider::DrawZoomImage(const byte *img)
|
|---|
| 802 | {
|
|---|
| 803 | byte zimg[kZOOM*kZOOM];
|
|---|
| 804 | for (int y=0; y<kZOOM; y++)
|
|---|
| 805 | for (int x=0; x<kZOOM; x++)
|
|---|
| 806 | zimg[x+y*kZOOM] = img[(fDx+(x-kZOOM/2)/2)+(fDy+(y-kZOOM/2)/2)*768];
|
|---|
| 807 |
|
|---|
| 808 | fZoomImage->DrawImg(zimg);
|
|---|
| 809 | }
|
|---|
| 810 |
|
|---|
| 811 | void MStarguider::DrawCosyImage(const byte *img)
|
|---|
| 812 | {
|
|---|
| 813 | if (!fCosy)
|
|---|
| 814 | return;
|
|---|
| 815 |
|
|---|
| 816 | byte simg[(768/2-1)*(576/2-1)];
|
|---|
| 817 | for (int y=0; y<576/2-1; y++)
|
|---|
| 818 | for (int x=0; x<768/2-1; x++)
|
|---|
| 819 | simg[x+y*(768/2-1)] = ((unsigned int)img[2*x+2*y*768]+img[2*x+2*y*768+1]+img[2*x+2*(y+1)*768]+img[2*x+2*(y+1)*768+1])/4;
|
|---|
| 820 |
|
|---|
| 821 | fCosy->GetWin()->GetImage()->DrawImg(simg);
|
|---|
| 822 | }
|
|---|
| 823 |
|
|---|
| 824 | Led MStarguider::FindStar(const FilterLed &f, const FilterLed &f2, const Ring ¢er, Int_t numleds, Int_t numrings)
|
|---|
| 825 | {
|
|---|
| 826 | // Get tracking coordinates
|
|---|
| 827 | //const XY xy = fCRaDec->GetCoordinates(); // [h, deg]
|
|---|
| 828 |
|
|---|
| 829 | if (center.GetX()<=0 && center.GetY()<=0)
|
|---|
| 830 | {
|
|---|
| 831 | cout << "Couldn't determine center of the camera." << endl;
|
|---|
| 832 | //if (fTPoint->IsDown() && fCosy && fCosy->GetDriveCom())
|
|---|
| 833 | // fCosy->GetDriveCom()->SendTPoint(false, 'T', fTPointStarMag, fTPointStarName, AltAz(), ZdAz(), xy, 0, 0, t, center, Led(), numleds, numrings); // Report
|
|---|
| 834 | return;
|
|---|
| 835 | }
|
|---|
| 836 |
|
|---|
| 837 | // Try to find the star
|
|---|
| 838 | Leds leds;
|
|---|
| 839 | f.FindStar(leds, (Int_t)center.GetX(), (Int_t)center.GetY(), true);
|
|---|
| 840 |
|
|---|
| 841 | // Check whether star found
|
|---|
| 842 | if (leds.size()==0)
|
|---|
| 843 | {
|
|---|
| 844 | cout << "No star found." << endl;
|
|---|
| 845 | //if (fTPoint->IsDown() && fCosy && fCosy->GetDriveCom())
|
|---|
| 846 | // fCosy->GetDriveCom()->SendTPoint(false, 'T', fTPointStarMag, fTPointStarName, AltAz(), ZdAz(), xy, 0, 0, t, center, Led(), numleds, numrings); // Report
|
|---|
| 847 | return;
|
|---|
| 848 | }
|
|---|
| 849 |
|
|---|
| 850 | cout << "Found star @ " << flush;
|
|---|
| 851 |
|
|---|
| 852 | const Led &star = leds.front();
|
|---|
| 853 |
|
|---|
| 854 | star.Print();
|
|---|
| 855 | f2.MarkPoint(star.GetX(), star.GetY(), 2<<2);
|
|---|
| 856 |
|
|---|
| 857 | return star;
|
|---|
| 858 |
|
|---|
| 859 | const RaDec rd(xy.X()*MAstro::HorToRad(), xy.Y()*TMath::DegToRad());
|
|---|
| 860 |
|
|---|
| 861 | // Initialize Star Catalog on the camera plane
|
|---|
| 862 | MGeomCamMagic geom;
|
|---|
| 863 | MAstroCamera ac;
|
|---|
| 864 | ac.SetGeom(geom);
|
|---|
| 865 | ac.SetRadiusFOV(3);
|
|---|
| 866 | //ac.SetObservatory(*fSao);
|
|---|
| 867 | ac.SetTime(t);
|
|---|
| 868 | ac.SetRaDec(rd.Ra(), rd.Dec());
|
|---|
| 869 |
|
|---|
| 870 | // Convert from Pixel to millimeter (1pix=2.6mm) [deg/pix / deg/mm = mm/pix]
|
|---|
| 871 | // Correct for abberation.
|
|---|
| 872 | const double sec_per_pix = 45.311;
|
|---|
| 873 | const double mm_to_deg = 0.011693;
|
|---|
| 874 | const Double_t conv = sec_per_pix/3600/mm_to_deg / 1.0713;
|
|---|
| 875 |
|
|---|
| 876 | // Adapt coordinate system (GUIs and humans are counting Y in different directions)
|
|---|
| 877 | const Double_t dx = (star->GetX()-center.GetX())*conv;
|
|---|
| 878 | const Double_t dy = (center.GetY()-star->GetY())*conv;
|
|---|
| 879 |
|
|---|
| 880 | // Convert offset from camera plane into local ccordinates
|
|---|
| 881 | Double_t dzd, daz;
|
|---|
| 882 | ac.GetDiffZdAz(dx, dy, dzd, daz);
|
|---|
| 883 |
|
|---|
| 884 | ZdAz zdaz(dzd,daz);
|
|---|
| 885 |
|
|---|
| 886 | // Check TPoint data set request
|
|---|
| 887 | if (!fMode->IsEntryChecked(IDM_kTpointMode) || !fTPoint->IsDown())
|
|---|
| 888 | return;
|
|---|
| 889 |
|
|---|
| 890 | // If no file open: open new file
|
|---|
| 891 | if (!fOutTp)
|
|---|
| 892 | {
|
|---|
| 893 | //
|
|---|
| 894 | // open tpoint file
|
|---|
| 895 | //
|
|---|
| 896 |
|
|---|
| 897 | const TString name = MCosy::GetFileName("tpoint", "tpoint", "txt");
|
|---|
| 898 | cout << "TPoint File ********* " << name << " ********** " << endl;
|
|---|
| 899 |
|
|---|
| 900 | fOutTp = new ofstream(name);
|
|---|
| 901 | *fOutTp << "Magic Model TPOINT data file" << endl;
|
|---|
| 902 | *fOutTp << ": ALTAZ" << endl;
|
|---|
| 903 | *fOutTp << "49 48 0 ";
|
|---|
| 904 | *fOutTp << t << endl;
|
|---|
| 905 | // temp(°C) pressure(mB) height(m) humidity(1) wavelength(microm) troplapserate(K/m
|
|---|
| 906 | }
|
|---|
| 907 |
|
|---|
| 908 | // Output Ra/Dec the drive system thinks that it is currently tracking
|
|---|
| 909 | //cout << "TPoint Star: " << xy.X() << "h " << xy.Y() << "°" << endl;
|
|---|
| 910 |
|
|---|
| 911 | // From the star position in the camera we calculate the Alt/Az
|
|---|
| 912 | // position we are currently tracking (real pointing position)
|
|---|
| 913 | fSao->SetMjd(t.GetMjd());
|
|---|
| 914 | AltAz za0 = fSao->CalcAltAz(rd)*kRad2Deg;
|
|---|
| 915 |
|
|---|
| 916 | //ZdAz za0 = fSao->GetZdAz();
|
|---|
| 917 | za0 -= AltAz(-dzd, daz);
|
|---|
| 918 | fAltAzOffsetFromTp = AltAz(-dzd, daz);
|
|---|
| 919 | fTimeFromTp=t;
|
|---|
| 920 |
|
|---|
| 921 | // From the Shaftencoders we get the current 'pointing' position
|
|---|
| 922 | // as it is seen by the drive system (system pointing position)
|
|---|
| 923 | const ZdAz za1 = fCosy->GetSePos()*360; // [deg]
|
|---|
| 924 |
|
|---|
| 925 |
|
|---|
| 926 | // Write real pointing position
|
|---|
| 927 | //cout << " Alt/Az: " << za0.Alt() << "° " << za0.Az() << "°" << endl;
|
|---|
| 928 | *fOutTp << setprecision(7) << za0.Az() << " " << za0.Alt() << " ";
|
|---|
| 929 |
|
|---|
| 930 | // Write system pointing position
|
|---|
| 931 | //cout << " SE-Pos: " << 90-za1.Zd() << "° " << za1.Az() << "°" << endl;
|
|---|
| 932 | *fOutTp << fmod(za1.Az()+360, 360) << " " << 90-za1.Zd();
|
|---|
| 933 |
|
|---|
| 934 | *fOutTp << " " << xy.X() << " " << xy.Y();
|
|---|
| 935 | *fOutTp << " " << -dzd << " " << -daz;
|
|---|
| 936 | *fOutTp << " " << setprecision(11) << t.GetMjd();
|
|---|
| 937 | *fOutTp << " " << setprecision(4) << center.GetMag();
|
|---|
| 938 | *fOutTp << " " << star->GetMag();
|
|---|
| 939 | *fOutTp << " " << center.GetX() << " " << center.GetY();
|
|---|
| 940 | *fOutTp << " " << star->GetX() << " " << star->GetY();
|
|---|
| 941 | *fOutTp << " " << numleds << " " << numrings;
|
|---|
| 942 | *fOutTp << " 0 0 0";
|
|---|
| 943 | *fOutTp << " " << fTPointStarMag << " " << fTPointStarName;
|
|---|
| 944 | *fOutTp << endl;
|
|---|
| 945 |
|
|---|
| 946 | gLog << all << "TPoint successfully taken." << endl;
|
|---|
| 947 |
|
|---|
| 948 | MLog &outrep = *fCosy->GetOutRep();
|
|---|
| 949 | if (outrep.Lock("MStarguider::FindStar"))
|
|---|
| 950 | {
|
|---|
| 951 | outrep << "FINDSTAR-REPORT 00 " << MTime(-1) << " " << setprecision(7);
|
|---|
| 952 | outrep << 90-za0.Alt() << " " << za0.Az() << " ";
|
|---|
| 953 | outrep << za1.Zd() << " " << za1.Az() << " ";
|
|---|
| 954 | outrep << xy.X() << " " << xy.Y() << " ";
|
|---|
| 955 | outrep << -dzd << " " << -daz << " ";
|
|---|
| 956 | outrep << star->GetX() << " " << star->GetY() << " ";
|
|---|
| 957 | outrep << center.GetX() << " " << center.GetY() << " ";
|
|---|
| 958 | outrep << dx/conv << " " << dy/conv << " " << star->GetMag();
|
|---|
| 959 | outrep << setprecision(11) << t.GetMjd() << endl;
|
|---|
| 960 | outrep.UnLock("MStarguider::FindStar");
|
|---|
| 961 | }
|
|---|
| 962 |
|
|---|
| 963 | // return zdaz;
|
|---|
| 964 |
|
|---|
| 965 | if (!fCosy)
|
|---|
| 966 | return;
|
|---|
| 967 |
|
|---|
| 968 | MDriveCom *com = fCosy->GetDriveCom();
|
|---|
| 969 | if (!com)
|
|---|
| 970 | return;
|
|---|
| 971 |
|
|---|
| 972 | com->SendTPoint(true, 'T', fTPointStarMag, fTPointStarName, za0, za1, xy, -dzd, -daz, t, center, *star, numleds, numrings); // Report
|
|---|
| 973 | }
|
|---|
| 974 | */
|
|---|
| 975 |
|
|---|
| 976 | void MStarguider::WritePNG(const char *name, const byte *gbuf, const byte *cbuf)
|
|---|
| 977 | {
|
|---|
| 978 | string fname(name);
|
|---|
| 979 | if (fname.size()<4 || fname.compare(fname.size()-4,4, ".png")==0)
|
|---|
| 980 | fname += ".png";
|
|---|
| 981 |
|
|---|
| 982 | cout << "Writing PNG '" << fname << "'" << endl;
|
|---|
| 983 |
|
|---|
| 984 | #ifdef HAVE_PNG
|
|---|
| 985 | //
|
|---|
| 986 | // open file
|
|---|
| 987 | //
|
|---|
| 988 | FILE *fd = fopen(fname.c_str(), "w");
|
|---|
| 989 | if (!fd)
|
|---|
| 990 | {
|
|---|
| 991 | cout << "Warning: Cannot open file for writing." << endl;
|
|---|
| 992 | return;
|
|---|
| 993 | }
|
|---|
| 994 |
|
|---|
| 995 | //
|
|---|
| 996 | // allocate memory
|
|---|
| 997 | //
|
|---|
| 998 | png_structp fPng = png_create_write_struct(PNG_LIBPNG_VER_STRING,
|
|---|
| 999 | NULL, NULL, NULL);
|
|---|
| 1000 |
|
|---|
| 1001 | if (!fPng)
|
|---|
| 1002 | {
|
|---|
| 1003 | cout << "Warning: Unable to create PNG structure" << endl;
|
|---|
| 1004 | fclose(fd);
|
|---|
| 1005 | return;
|
|---|
| 1006 | }
|
|---|
| 1007 |
|
|---|
| 1008 |
|
|---|
| 1009 | png_infop fInfo = png_create_info_struct(fPng);
|
|---|
| 1010 |
|
|---|
| 1011 | if (!fInfo)
|
|---|
| 1012 | {
|
|---|
| 1013 | cout << "Warning: Unable to create PNG info structure" << endl;
|
|---|
| 1014 | png_destroy_write_struct (&fPng, NULL);
|
|---|
| 1015 | fclose(fd);
|
|---|
| 1016 | return;
|
|---|
| 1017 | }
|
|---|
| 1018 |
|
|---|
| 1019 | fInfo->width = 768;
|
|---|
| 1020 | fInfo->height = 576;
|
|---|
| 1021 | fInfo->bit_depth = 8;
|
|---|
| 1022 | fInfo->color_type = PNG_COLOR_TYPE_RGB;
|
|---|
| 1023 | // fInfo->color_type = PNG_COLOR_TYPE_GRAY;
|
|---|
| 1024 |
|
|---|
| 1025 | //
|
|---|
| 1026 | // set jump-back point in case of errors
|
|---|
| 1027 | //
|
|---|
| 1028 | if (setjmp(fPng->jmpbuf))
|
|---|
| 1029 | {
|
|---|
| 1030 | cout << "longjmp Warning: PNG encounterd an error!" << endl;
|
|---|
| 1031 | png_destroy_write_struct (&fPng, &fInfo);
|
|---|
| 1032 | fclose(fd);
|
|---|
| 1033 | return;
|
|---|
| 1034 | }
|
|---|
| 1035 |
|
|---|
| 1036 | //
|
|---|
| 1037 | // connect file to PNG-Structure
|
|---|
| 1038 | //
|
|---|
| 1039 | png_init_io(fPng, fd);
|
|---|
| 1040 |
|
|---|
| 1041 | // png_set_compression_level (fPng, Z_BEST_COMPRESSION);
|
|---|
| 1042 |
|
|---|
| 1043 | //
|
|---|
| 1044 | // Write header
|
|---|
| 1045 | //
|
|---|
| 1046 | png_write_info(fPng, fInfo);
|
|---|
| 1047 |
|
|---|
| 1048 | png_byte buf[768*576*3];
|
|---|
| 1049 |
|
|---|
| 1050 | png_byte *d = buf;
|
|---|
| 1051 | const byte *g = gbuf;
|
|---|
| 1052 | const byte *c = cbuf;
|
|---|
| 1053 |
|
|---|
| 1054 | // d=destination, s1=source1, s2=source2, e=end
|
|---|
| 1055 | while (d<buf+768*576*3)
|
|---|
| 1056 | {
|
|---|
| 1057 | if (fScreenshotColor && *c)
|
|---|
| 1058 | {
|
|---|
| 1059 | *d++ = ((*c>>4)&0x3)*85;
|
|---|
| 1060 | *d++ = ((*c>>2)&0x3)*85;
|
|---|
| 1061 | *d++ = ((*c++ )&0x3)*85;
|
|---|
| 1062 | g++;
|
|---|
| 1063 | }
|
|---|
| 1064 | else
|
|---|
| 1065 | {
|
|---|
| 1066 | *d++ = *g;
|
|---|
| 1067 | *d++ = *g;
|
|---|
| 1068 | *d++ = *g++;
|
|---|
| 1069 | c++;
|
|---|
| 1070 | }
|
|---|
| 1071 | }
|
|---|
| 1072 |
|
|---|
| 1073 | //
|
|---|
| 1074 | // Write bitmap data
|
|---|
| 1075 | //
|
|---|
| 1076 | for (unsigned int y=0; y<768*576*3; y+=768*3)
|
|---|
| 1077 | png_write_row (fPng, buf+y);
|
|---|
| 1078 |
|
|---|
| 1079 | //
|
|---|
| 1080 | // Write footer
|
|---|
| 1081 | //
|
|---|
| 1082 | png_write_end (fPng, fInfo);
|
|---|
| 1083 |
|
|---|
| 1084 | //
|
|---|
| 1085 | // free memory
|
|---|
| 1086 | //
|
|---|
| 1087 | png_destroy_write_struct (&fPng, &fInfo);
|
|---|
| 1088 |
|
|---|
| 1089 | fclose(fd);
|
|---|
| 1090 | #else
|
|---|
| 1091 | cout << "Sorry, no png support compiled into tpoint." << endl;
|
|---|
| 1092 | #endif
|
|---|
| 1093 | }
|
|---|
| 1094 |
|
|---|
| 1095 | bool MStarguider::Interpolate(const unsigned long n, byte *img) const
|
|---|
| 1096 | {
|
|---|
| 1097 | if (fIntRate<=1)
|
|---|
| 1098 | return true;
|
|---|
| 1099 |
|
|---|
| 1100 | const int rate = fTPointMode>0 ? 5*25 : fIntRate;
|
|---|
| 1101 |
|
|---|
| 1102 | static unsigned short myimg[768*576];
|
|---|
| 1103 |
|
|---|
| 1104 | unsigned short *f = myimg;
|
|---|
| 1105 |
|
|---|
| 1106 | const byte *end = img+768*576;
|
|---|
| 1107 |
|
|---|
| 1108 | if (n%rate)
|
|---|
| 1109 | {
|
|---|
| 1110 | while (img<end)
|
|---|
| 1111 | *f++ += *img++;
|
|---|
| 1112 | return false;
|
|---|
| 1113 | }
|
|---|
| 1114 | else
|
|---|
| 1115 | {
|
|---|
| 1116 | while (img<end)
|
|---|
| 1117 | {
|
|---|
| 1118 | *img = (*img + *f)/rate;
|
|---|
| 1119 | ++img;
|
|---|
| 1120 | *f++ = 0;
|
|---|
| 1121 | }
|
|---|
| 1122 |
|
|---|
| 1123 | return true;
|
|---|
| 1124 | }
|
|---|
| 1125 | }
|
|---|
| 1126 |
|
|---|
| 1127 | void MStarguider::ProcessFrame(const unsigned long n, byte *img,
|
|---|
| 1128 | struct timeval *tm)
|
|---|
| 1129 | {
|
|---|
| 1130 | if (!Interpolate(n, img))
|
|---|
| 1131 | return;
|
|---|
| 1132 |
|
|---|
| 1133 | if (fTPointMode==2)
|
|---|
| 1134 | {
|
|---|
| 1135 | fTPointMode=1;
|
|---|
| 1136 | return;
|
|---|
| 1137 | }
|
|---|
| 1138 |
|
|---|
| 1139 | byte cimg[768*576];
|
|---|
| 1140 | memset(cimg, 0, 768*576);
|
|---|
| 1141 |
|
|---|
| 1142 | FilterLed f (img, 768, 576, 2.5); // 2.5
|
|---|
| 1143 | FilterLed f2(cimg, 768, 576); // former color 0xb0
|
|---|
| 1144 |
|
|---|
| 1145 | if (!fTPointMode && fScreenshotName.empty() && fDisplay->IsEntryChecked(IDM_kStretch))
|
|---|
| 1146 | f.Stretch();
|
|---|
| 1147 |
|
|---|
| 1148 | // Visual Filter, whole FOV
|
|---|
| 1149 | if (!fTPointMode && fDisplay->IsEntryChecked(IDM_kFilter))
|
|---|
| 1150 | {
|
|---|
| 1151 | vector<Led> leds;
|
|---|
| 1152 | f.Execute(leds, 768/2, 576/2);
|
|---|
| 1153 | for (auto it=leds.begin(); it!=leds.end(); it++)
|
|---|
| 1154 | f.MarkPoint(*it);
|
|---|
| 1155 | }
|
|---|
| 1156 |
|
|---|
| 1157 | // Find Center of Camera for Caos and Tpoints
|
|---|
| 1158 | int numleds = 0;
|
|---|
| 1159 | int numrings = 0;
|
|---|
| 1160 | Ring center(-1, -1);//(5, 5);
|
|---|
| 1161 |
|
|---|
| 1162 | if (fTPointMode || fDisplay->IsEntryChecked(IDM_kCaosFilter))
|
|---|
| 1163 | {
|
|---|
| 1164 | center = fCaos->Run(img);
|
|---|
| 1165 | numleds = fCaos->GetNumDetectedLEDs();
|
|---|
| 1166 | numrings = fCaos->GetNumDetectedRings();
|
|---|
| 1167 | }
|
|---|
| 1168 |
|
|---|
| 1169 | //cout << "cx=" << center.GetX() << " cy=" << center.GetY() << " Nled=" << numleds << " Nrings=" << numrings << endl;
|
|---|
| 1170 |
|
|---|
| 1171 | // Find Star at Center---for Tpoint Procedure
|
|---|
| 1172 | Led star(-1, -1);
|
|---|
| 1173 | if (center.GetX()>0 && center.GetY()>0)
|
|---|
| 1174 | {
|
|---|
| 1175 | if (fTPointMode || fDisplay->IsEntryChecked(IDM_kFindStar))
|
|---|
| 1176 | {
|
|---|
| 1177 | // Set search Paremeters (FIXME: Get them from user input!)
|
|---|
| 1178 | f.SetCut(fFindStarCut+1.5); // FindStar.CleaningLevel
|
|---|
| 1179 | f.SetBox(fFindStarBox+12/*+80*/); // FindStar.SizeBox
|
|---|
| 1180 |
|
|---|
| 1181 | // Try to find the star
|
|---|
| 1182 | vector<Led> leds;
|
|---|
| 1183 | f.FindStar(leds, (Int_t)center.GetX(), (Int_t)center.GetY(), true);
|
|---|
| 1184 |
|
|---|
| 1185 | // Check whether star found
|
|---|
| 1186 | if (leds.size()>0)
|
|---|
| 1187 | {
|
|---|
| 1188 | //cout << "Found star @ " << flush;
|
|---|
| 1189 | //leds[0].Print();
|
|---|
| 1190 | f2.MarkPoint(leds[0].GetX(), leds[0].GetY(), 2<<2);
|
|---|
| 1191 | star = leds[0];
|
|---|
| 1192 | }
|
|---|
| 1193 | }
|
|---|
| 1194 |
|
|---|
| 1195 | // DrawZoomImage(img);
|
|---|
| 1196 | // DrawCosyImage(img);
|
|---|
| 1197 |
|
|---|
| 1198 | // Position corresponding to the camera center (53.2, 293.6)
|
|---|
| 1199 | // Draw Circles around center of Camera
|
|---|
| 1200 | if (fTPointMode || fDisplay->IsEntryChecked(IDM_kCaosFilter))
|
|---|
| 1201 | {
|
|---|
| 1202 | f2.DrawCircle(center, 0x0a);
|
|---|
| 1203 | f2.DrawCircle(center, 7.0,
|
|---|
| 1204 | fDisplay->IsEntryChecked(IDM_kFindStar)?3:0xb0);
|
|---|
| 1205 | //f2.DrawCircle(center, 115.0, 0x0a);
|
|---|
| 1206 | //f2.DrawCircle(center, 230.0, 0x0a);
|
|---|
| 1207 | //f2.DrawCircle(center, 245.0, 0x0a);
|
|---|
| 1208 | }
|
|---|
| 1209 | }
|
|---|
| 1210 |
|
|---|
| 1211 | if (fTPointMode ||
|
|---|
| 1212 | fDisplay->IsEntryChecked(IDM_kCaosFilter) ||
|
|---|
| 1213 | fDisplay->IsEntryChecked(IDM_kFindStar))
|
|---|
| 1214 | fImage->DrawColImg(img, cimg);
|
|---|
| 1215 | else
|
|---|
| 1216 | fImage->DrawImg(img);
|
|---|
| 1217 |
|
|---|
| 1218 | if (fTPointMode && !fScreenshotName.empty())
|
|---|
| 1219 | {
|
|---|
| 1220 | WritePNG(fScreenshotName.c_str(), img, cimg);
|
|---|
| 1221 | fScreenshotName = "";
|
|---|
| 1222 | fTPointMode = 0;
|
|---|
| 1223 | return;
|
|---|
| 1224 | }
|
|---|
| 1225 |
|
|---|
| 1226 |
|
|---|
| 1227 | if (star.GetX()<0 || star.GetY()<0 || fTPointMode==0)
|
|---|
| 1228 | return;
|
|---|
| 1229 |
|
|---|
| 1230 | if (fTPointMode==1)
|
|---|
| 1231 | fTPointMode=0;
|
|---|
| 1232 |
|
|---|
| 1233 | // Convert from Pixel to millimeter (1pix=2.6mm) [deg/pix / deg/mm = mm/pix]
|
|---|
| 1234 | // Correct for abberation.
|
|---|
| 1235 |
|
|---|
| 1236 | //const float Dleds = 510; // 5.96344 deg
|
|---|
| 1237 | //const float Dpix = 2*237.58;
|
|---|
| 1238 |
|
|---|
| 1239 | // The DC reflector elongates light from off axis sources
|
|---|
| 1240 | // This is a correction. It is 7% for MAGIC (1:1) and
|
|---|
| 1241 | // less for FACT (1:1.4). This is an estimate from a
|
|---|
| 1242 | // Orbit mode observation at 0.17deg distance to the
|
|---|
| 1243 | // camera center [4.90m might also not be very accurate
|
|---|
| 1244 | // depending on the position of the CCD camera]
|
|---|
| 1245 | const double abberation = 1.0638; //1.0713;
|
|---|
| 1246 | const double sec_per_pix = 45.14; //45.311; (atan(510mm/2 / 4.90m) / 237.58) // FACT LEDs
|
|---|
| 1247 |
|
|---|
| 1248 | const double conv = sec_per_pix/abberation;
|
|---|
| 1249 |
|
|---|
| 1250 | const double dx = star.GetX()-center.GetX();
|
|---|
| 1251 | const double dy = star.GetY()-center.GetY();
|
|---|
| 1252 |
|
|---|
| 1253 | //const double dxx = - conv * (star.GetX()-center.GetX());
|
|---|
| 1254 | //const double dyy = conv * (star.GetY()-center.GetY());
|
|---|
| 1255 |
|
|---|
| 1256 | const double dphi = - center.GetPhi() * M_PI/180;
|
|---|
| 1257 |
|
|---|
| 1258 | // The sign is because the pixels are not counted in
|
|---|
| 1259 | // both directions in the same direction as Zd/Az
|
|---|
| 1260 | const double dxx = - conv * (dx*cos(dphi) - dy*sin(dphi));
|
|---|
| 1261 | const double dyy = conv * (dx*sin(dphi) + dy*cos(dphi));
|
|---|
| 1262 |
|
|---|
| 1263 | double arr[11] = {
|
|---|
| 1264 | dxx, dyy,
|
|---|
| 1265 | double(numleds), double(numrings),
|
|---|
| 1266 | center.GetX(), center.GetY(), center.GetMag(),
|
|---|
| 1267 | star.GetX(), star.GetY(), star.GetMag(),
|
|---|
| 1268 | center.GetPhi()
|
|---|
| 1269 | };
|
|---|
| 1270 |
|
|---|
| 1271 | fDimData.setData(arr, 11*sizeof(double));
|
|---|
| 1272 | fDimData.setQuality(0);
|
|---|
| 1273 | fDimData.setTimestamp(tm->tv_sec, tm->tv_usec/1000);
|
|---|
| 1274 | fDimData.updateService();
|
|---|
| 1275 | }
|
|---|
| 1276 |
|
|---|
| 1277 | /*
|
|---|
| 1278 | void MStarguider::UpdatePosZoom()
|
|---|
| 1279 | {
|
|---|
| 1280 | MString txt;
|
|---|
| 1281 | if (fDisplay->IsEntryChecked(IDM_kCatalog))
|
|---|
| 1282 | {
|
|---|
| 1283 | // FIXME: Necessary?
|
|---|
| 1284 | fSao->Now();
|
|---|
| 1285 | AltAz aa = fSao->CalcAltAzFromPix(fDx, fDy)*kRad2Deg;
|
|---|
| 1286 | if (aa.Az()<0)
|
|---|
| 1287 | aa.Az(aa.Az()+360);
|
|---|
| 1288 | txt.Form("(%d, %d) %.1fd/%.1fd", fDx, fDy, -aa.Alt(), aa.Az()-180);
|
|---|
| 1289 | }
|
|---|
| 1290 | else
|
|---|
| 1291 | txt.Form("(%d, %d)", fDx, fDy);
|
|---|
| 1292 | fPosZoom->SetText(txt);
|
|---|
| 1293 | }
|
|---|
| 1294 |
|
|---|
| 1295 | Bool_t MStarguider::HandleDoubleClick(Event_t *event)
|
|---|
| 1296 | {
|
|---|
| 1297 | const Int_t w = fImage->GetWidth();
|
|---|
| 1298 | const Int_t h = fImage->GetHeight();
|
|---|
| 1299 | const Int_t x = fImage->GetX();
|
|---|
| 1300 | const Int_t y = fImage->GetY();
|
|---|
| 1301 |
|
|---|
| 1302 | if (!(event->fX>x && event->fX<x+w && event->fY>y && event->fY<y+h))
|
|---|
| 1303 | return kTRUE;
|
|---|
| 1304 |
|
|---|
| 1305 | Int_t dx = event->fX-x;
|
|---|
| 1306 | Int_t dy = event->fY-y;
|
|---|
| 1307 |
|
|---|
| 1308 | if (dx<kZOOM/4) dx=kZOOM/4;
|
|---|
| 1309 | if (dy<kZOOM/4) dy=kZOOM/4;
|
|---|
| 1310 | if (dx>766-kZOOM/4) dx=766-kZOOM/4;
|
|---|
| 1311 | if (dy>574-kZOOM/4) dy=574-kZOOM/4;
|
|---|
| 1312 |
|
|---|
| 1313 | fDx = dx;
|
|---|
| 1314 | fDy = dy;
|
|---|
| 1315 |
|
|---|
| 1316 | //UpdatePosZoom();
|
|---|
| 1317 | return kTRUE;
|
|---|
| 1318 | }
|
|---|
| 1319 |
|
|---|
| 1320 | */
|
|---|
| 1321 |
|
|---|
| 1322 | /// Overwritten DimCommand::commandHandler
|
|---|
| 1323 | void MStarguider::commandHandler()
|
|---|
| 1324 | {
|
|---|
| 1325 | DimCommand *cmd = getCommand();
|
|---|
| 1326 | if (!cmd)
|
|---|
| 1327 | return;
|
|---|
| 1328 |
|
|---|
| 1329 | if (cmd==&fDimTPoint)
|
|---|
| 1330 | {
|
|---|
| 1331 | fTPointMode = 2;
|
|---|
| 1332 | fScreenshotName = "";
|
|---|
| 1333 | cout << "DimCommand[TPOINT]: " << cmd->itsSize << " " << string((char*)cmd->itsData, cmd->itsSize) << endl;
|
|---|
| 1334 | return;
|
|---|
| 1335 | }
|
|---|
| 1336 |
|
|---|
| 1337 | if (cmd==&fDimScreenshot && fTPointMode==0)
|
|---|
| 1338 | {
|
|---|
| 1339 | if (cmd->itsSize<2)
|
|---|
| 1340 | return;
|
|---|
| 1341 |
|
|---|
| 1342 | fTPointMode = 2;
|
|---|
| 1343 | fScreenshotColor = ((uint8_t*)cmd->itsData)[0];
|
|---|
| 1344 | fScreenshotName = string((char*)cmd->itsData+1, cmd->itsSize-1);
|
|---|
| 1345 | cout << "DimCommand[SCREENSHOT]: " << fScreenshotName << endl;
|
|---|
| 1346 | return;
|
|---|
| 1347 | }
|
|---|
| 1348 |
|
|---|
| 1349 | cout << "DimCommand[UNKNOWN]: " << cmd->itsSize << " " << string((char*)cmd->itsData, cmd->itsSize) << endl;
|
|---|
| 1350 | }
|
|---|