| 1 | #include "BasicGlCamera.h"
|
|---|
| 2 | #include "BasicGlCamera.h"
|
|---|
| 3 | #include <math.h>
|
|---|
| 4 | #include <fstream>
|
|---|
| 5 | #include <iostream>
|
|---|
| 6 | #include <string>
|
|---|
| 7 | #include <sstream>
|
|---|
| 8 | #include "src/tools.h"
|
|---|
| 9 |
|
|---|
| 10 | #include <algorithm>
|
|---|
| 11 | #include "src/Time.h"
|
|---|
| 12 |
|
|---|
| 13 | using namespace std;;
|
|---|
| 14 |
|
|---|
| 15 | //static variables
|
|---|
| 16 | PixelMap BasicGlCamera::fPixelMap;
|
|---|
| 17 | GLfloat BasicGlCamera::pixelsCoords[MAX_NUM_PIXELS][3];
|
|---|
| 18 | PixelsNeighbors BasicGlCamera::neighbors[MAX_NUM_PIXELS];
|
|---|
| 19 | int BasicGlCamera::hardwareMapping[NPIX];
|
|---|
| 20 | GLfloat BasicGlCamera::verticesList[NPIX*6][2];
|
|---|
| 21 | vector<edge> BasicGlCamera::patchesIndices[160];
|
|---|
| 22 | int BasicGlCamera::verticesIndices[NPIX][6];
|
|---|
| 23 | int BasicGlCamera::pixelsPatch[NPIX];
|
|---|
| 24 | int BasicGlCamera::softwareMapping[NPIX];
|
|---|
| 25 |
|
|---|
| 26 | //Coordinates of an hexagon of radius 1 and center 0
|
|---|
| 27 | GLfloat hexcoords[6][2] = {{-1./sqrt(3.), 1},
|
|---|
| 28 | { 1./sqrt(3.), 1},
|
|---|
| 29 | { 2./sqrt(3.), 0},
|
|---|
| 30 | { 1./sqrt(3.), -1},
|
|---|
| 31 | {-1./sqrt(3.), -1},
|
|---|
| 32 | {-2./sqrt(3.), 0}};
|
|---|
| 33 |
|
|---|
| 34 |
|
|---|
| 35 |
|
|---|
| 36 | BasicGlCamera::BasicGlCamera(QWidget* cParent)
|
|---|
| 37 | : QGLWidget(QGLFormat(QGL::DoubleBuffer |
|
|---|
| 38 | QGL::DepthBuffer |
|
|---|
| 39 | QGL::IndirectRendering),cParent)
|
|---|
| 40 | {
|
|---|
| 41 | QGL::setPreferredPaintEngine(QPaintEngine::OpenGL);
|
|---|
| 42 | fWhite = -1;
|
|---|
| 43 | fWhitePatch = -1;
|
|---|
| 44 | fMin = -1;
|
|---|
| 45 | fMax = -1;
|
|---|
| 46 | fScaleLimit = -0.5;
|
|---|
| 47 | fTextSize = 0;
|
|---|
| 48 | autoRefresh = false;
|
|---|
| 49 | logScale = false;
|
|---|
| 50 | cameraRotation = +90;
|
|---|
| 51 | fTextEnabled = true;
|
|---|
| 52 | unitsText = "";
|
|---|
| 53 | titleText = "";
|
|---|
| 54 | dataText = "";
|
|---|
| 55 | pixelContourColour[0] = 0.1f;
|
|---|
| 56 | pixelContourColour[1] = 0.1f;
|
|---|
| 57 | pixelContourColour[2] = 0.1f;
|
|---|
| 58 | patchesCoulour[0] = 0.1f;
|
|---|
| 59 | patchesCoulour[1] = 0.1f;
|
|---|
| 60 | patchesCoulour[2] = 0.1f;
|
|---|
| 61 | highlightedPatchesCoulour[0] = 0.6f;
|
|---|
| 62 | highlightedPatchesCoulour[1] = 0.6f;
|
|---|
| 63 | highlightedPatchesCoulour[2] = 0.6f;
|
|---|
| 64 | highlightedPixelsCoulour[0] = 0.8f;
|
|---|
| 65 | highlightedPixelsCoulour[1] = 0.8f;
|
|---|
| 66 | highlightedPixelsCoulour[2] = 0.8f;
|
|---|
| 67 |
|
|---|
| 68 | regularPalettePlease(true);
|
|---|
| 69 |
|
|---|
| 70 |
|
|---|
| 71 | hexRadius = 0.015f;
|
|---|
| 72 | hexTolerance = hexRadius/100.0f;
|
|---|
| 73 | viewSize = 1.0f;
|
|---|
| 74 | calculatePixelsCoords();
|
|---|
| 75 |
|
|---|
| 76 | /* ifstream fin2("MasterList-v3.txt");
|
|---|
| 77 | if (!fin2.is_open())
|
|---|
| 78 | {
|
|---|
| 79 | cout << "Error: file \"MasterList-v3\" missing. aborting." << endl;
|
|---|
| 80 | exit(-1);
|
|---|
| 81 | }
|
|---|
| 82 | int l = 0;
|
|---|
| 83 | string buf;
|
|---|
| 84 | while (getline(fin2, buf, '\n'))
|
|---|
| 85 | {
|
|---|
| 86 | buf = Tools::Trim(buf);
|
|---|
| 87 | if (buf[0]=='#')
|
|---|
| 88 | continue;
|
|---|
| 89 |
|
|---|
| 90 | unsigned int softid, hardid, dummy;
|
|---|
| 91 |
|
|---|
| 92 | stringstream str(buf);
|
|---|
| 93 |
|
|---|
| 94 | str >> softid;
|
|---|
| 95 | str >> dummy;
|
|---|
| 96 | str >> hardid;
|
|---|
| 97 |
|
|---|
| 98 | if (softid>=1440)
|
|---|
| 99 | continue;
|
|---|
| 100 |
|
|---|
| 101 | hardwareMapping[softid] = hardid;
|
|---|
| 102 | softwareMapping[hardid] = softid;
|
|---|
| 103 |
|
|---|
| 104 | l++;
|
|---|
| 105 | }*/
|
|---|
| 106 | buildVerticesList();
|
|---|
| 107 | /*
|
|---|
| 108 | ifstream fin1("Trigger-Patches.txt");
|
|---|
| 109 | if (!fin1.is_open())
|
|---|
| 110 | {
|
|---|
| 111 | cout << "Error: file \"Trigger-Patches.txt\" missing. Aborting." << endl;
|
|---|
| 112 | exit(-1);
|
|---|
| 113 | }
|
|---|
| 114 | l=0;
|
|---|
| 115 | while (getline(fin1, buf, '\n'))
|
|---|
| 116 | {
|
|---|
| 117 | buf = Tools::Trim(buf);
|
|---|
| 118 | if (buf[0]=='#')
|
|---|
| 119 | continue;
|
|---|
| 120 |
|
|---|
| 121 | stringstream str(buf);
|
|---|
| 122 | for (int i=0; i<9; i++)
|
|---|
| 123 | {
|
|---|
| 124 | unsigned int n;
|
|---|
| 125 | str >> n;
|
|---|
| 126 |
|
|---|
| 127 | if (n>=1440)
|
|---|
| 128 | continue;
|
|---|
| 129 |
|
|---|
| 130 | patches[l][i] = hardwareMapping[n];
|
|---|
| 131 | }
|
|---|
| 132 | l++;
|
|---|
| 133 | }
|
|---|
| 134 |
|
|---|
| 135 | //now construct the correspondance between pixels and patches
|
|---|
| 136 | for (int i=0;i<NTMARK;i++)
|
|---|
| 137 | for (int j=0;j<9;j++)
|
|---|
| 138 | pixelsPatch[softwareMapping[patches[i][j]]] = i;
|
|---|
| 139 |
|
|---|
| 140 | for (int i=0;i<1440;i++)
|
|---|
| 141 | updateNeighbors(i);
|
|---|
| 142 |
|
|---|
| 143 | buildPatchesIndices();
|
|---|
| 144 | */////////////////////////////////
|
|---|
| 145 | ss[0] = 0; ss[1] = 0.25f; ss[2] = 0.5f; ss[3] = 0.75f; ss[4] = 1.0f;
|
|---|
| 146 | rr[0] = 0.15; rr[1] = 0; rr[2] = 0; rr[3] = 1.0f; rr[4] = 0.85f;
|
|---|
| 147 | gg[0] = 0.15; gg[1] = 0; gg[2] = 1; gg[3] = 0; gg[4] = 0.85f;
|
|---|
| 148 | bb[0] = 0.15; bb[1] = 1; bb[2] = 0; bb[3] = 0; bb[4] = 0.85f;
|
|---|
| 149 |
|
|---|
| 150 | fPixelStride = 1;
|
|---|
| 151 | fcSlice = 0;
|
|---|
| 152 | fData.resize(1440);
|
|---|
| 153 | for (int i=0;i<NPIX;i++)
|
|---|
| 154 | fData[i] = (double)i;///1.44;//(double)(i)/(double)(ACTUAL_NUM_PIXELS);
|
|---|
| 155 |
|
|---|
| 156 | // setFont(QFont("Arial", 8));
|
|---|
| 157 | int buttonShift=0;
|
|---|
| 158 | scaleLabel = new QLabel("Scale", this);
|
|---|
| 159 | // buttonShift += scaleLabel->height();
|
|---|
| 160 |
|
|---|
| 161 | linearButton = new QRadioButton("Linear", this);
|
|---|
| 162 | linearButton->move(scaleLabel->width(), buttonShift);
|
|---|
| 163 | buttonShift += linearButton->height();
|
|---|
| 164 |
|
|---|
| 165 | logButton = new QRadioButton("Log", this);
|
|---|
| 166 | logButton->move(scaleLabel->width(), buttonShift);
|
|---|
| 167 | buttonShift += logButton->height()*1.1f;
|
|---|
| 168 |
|
|---|
| 169 | colorPaletteLabel = new QLabel("Colour\nPalette", this);
|
|---|
| 170 | colorPaletteLabel->move(0, buttonShift);
|
|---|
| 171 | // buttonShift += colorPaletteLabel->height();
|
|---|
| 172 |
|
|---|
| 173 | regularPaletteButton = new QRadioButton("Regular", this);
|
|---|
| 174 | regularPaletteButton->move(colorPaletteLabel->width(), buttonShift);
|
|---|
| 175 | buttonShift += regularPaletteButton->height();
|
|---|
| 176 |
|
|---|
| 177 | prettyPaletteButton = new QRadioButton("Pretty", this);
|
|---|
| 178 | prettyPaletteButton->move(colorPaletteLabel->width(), buttonShift);
|
|---|
| 179 | buttonShift += prettyPaletteButton->height();
|
|---|
| 180 |
|
|---|
| 181 | greyScalePaletteButton = new QRadioButton("Grey Scale", this);
|
|---|
| 182 | greyScalePaletteButton->move(colorPaletteLabel->width(), buttonShift);
|
|---|
| 183 | buttonShift += greyScalePaletteButton->height();
|
|---|
| 184 |
|
|---|
| 185 | glowingPaletteButton = new QRadioButton("Glowing", this);
|
|---|
| 186 | glowingPaletteButton->move(colorPaletteLabel->width(), buttonShift);
|
|---|
| 187 | buttonShift += glowingPaletteButton->height()*1.1f;
|
|---|
| 188 |
|
|---|
| 189 | rotationLabel = new QLabel("Camera\nRotation", this);
|
|---|
| 190 | rotationLabel->move(0, buttonShift);
|
|---|
| 191 | // buttonShift += rotationLabel->height();
|
|---|
| 192 |
|
|---|
| 193 | unsigned short utf16Array;
|
|---|
| 194 | utf16Array = 0x00b0;
|
|---|
| 195 | QString degreeSymbol(QString::fromUtf16(&utf16Array, 1));
|
|---|
| 196 | QString zerostr("0" + degreeSymbol);
|
|---|
| 197 | zeroRotationButton = new QRadioButton(zerostr, this);
|
|---|
| 198 | zeroRotationButton->move(rotationLabel->width(), buttonShift);
|
|---|
| 199 | buttonShift += zeroRotationButton->height();
|
|---|
| 200 | QString minus90str("+90" + degreeSymbol);
|
|---|
| 201 | minus90RotationButton = new QRadioButton(minus90str, this);
|
|---|
| 202 | minus90RotationButton->move(rotationLabel->width(), buttonShift);
|
|---|
| 203 | buttonShift += minus90RotationButton->height();
|
|---|
| 204 | QString plus90str("-90"+degreeSymbol);
|
|---|
| 205 | plus90Rotationbutton = new QRadioButton(plus90str, this);
|
|---|
| 206 | plus90Rotationbutton->move(rotationLabel->width(), buttonShift);
|
|---|
| 207 |
|
|---|
| 208 |
|
|---|
| 209 | scaleGroup = new QButtonGroup(this);
|
|---|
| 210 | colorGroup = new QButtonGroup(this);
|
|---|
| 211 | rotationGroup = new QButtonGroup(this);
|
|---|
| 212 | scaleGroup->addButton(linearButton);
|
|---|
| 213 | scaleGroup->addButton(logButton);
|
|---|
| 214 | colorGroup->addButton(regularPaletteButton);
|
|---|
| 215 | colorGroup->addButton(prettyPaletteButton);
|
|---|
| 216 | colorGroup->addButton(greyScalePaletteButton);
|
|---|
| 217 | colorGroup->addButton(glowingPaletteButton);
|
|---|
| 218 | rotationGroup->addButton(zeroRotationButton);
|
|---|
| 219 | rotationGroup->addButton(minus90RotationButton);
|
|---|
| 220 | rotationGroup->addButton(plus90Rotationbutton);
|
|---|
| 221 |
|
|---|
| 222 | linearButton->setChecked(true);
|
|---|
| 223 | regularPaletteButton->setChecked(true);
|
|---|
| 224 | // zeroRotationButton->setChecked(true);
|
|---|
| 225 | minus90RotationButton->setChecked(true);
|
|---|
| 226 | // linearButton->palette.setColor();
|
|---|
| 227 |
|
|---|
| 228 | linearButton->setAutoFillBackground(true);
|
|---|
| 229 | logButton->setAutoFillBackground(true);
|
|---|
| 230 | regularPaletteButton->setAutoFillBackground(true);
|
|---|
| 231 | prettyPaletteButton->setAutoFillBackground(true);
|
|---|
| 232 | greyScalePaletteButton->setAutoFillBackground(true);
|
|---|
| 233 | glowingPaletteButton->setAutoFillBackground(true);
|
|---|
| 234 | zeroRotationButton->setAutoFillBackground(true);
|
|---|
| 235 | minus90RotationButton->setAutoFillBackground(true);
|
|---|
| 236 | plus90Rotationbutton->setAutoFillBackground(true);
|
|---|
| 237 | scaleLabel->setAutoFillBackground(true);
|
|---|
| 238 | colorPaletteLabel->setAutoFillBackground(true);
|
|---|
| 239 | rotationLabel->setAutoFillBackground(true);
|
|---|
| 240 |
|
|---|
| 241 | linearButton->hide();
|
|---|
| 242 | logButton->hide();
|
|---|
| 243 | regularPaletteButton->hide();
|
|---|
| 244 | prettyPaletteButton->hide();
|
|---|
| 245 | greyScalePaletteButton->hide();
|
|---|
| 246 | glowingPaletteButton->hide();
|
|---|
| 247 | zeroRotationButton->hide();
|
|---|
| 248 | minus90RotationButton->hide();
|
|---|
| 249 | plus90Rotationbutton->hide();
|
|---|
| 250 | scaleLabel->hide();
|
|---|
| 251 | colorPaletteLabel->hide();
|
|---|
| 252 | rotationLabel->hide();
|
|---|
| 253 |
|
|---|
| 254 | connect(linearButton, SIGNAL(toggled(bool)),
|
|---|
| 255 | this, SLOT(linearScalePlease(bool)));
|
|---|
| 256 | connect(logButton, SIGNAL(toggled(bool)),
|
|---|
| 257 | this, SLOT(logScalePlease(bool)));
|
|---|
| 258 | connect(regularPaletteButton, SIGNAL(toggled(bool)),
|
|---|
| 259 | this, SLOT(regularPalettePlease(bool)));
|
|---|
| 260 | connect(prettyPaletteButton, SIGNAL(toggled(bool)),
|
|---|
| 261 | this, SLOT(prettyPalettePlease(bool)));
|
|---|
| 262 | connect(greyScalePaletteButton, SIGNAL(toggled(bool)),
|
|---|
| 263 | this, SLOT(greyScalePalettePlease(bool)));
|
|---|
| 264 | connect(glowingPaletteButton, SIGNAL(toggled(bool)),
|
|---|
| 265 | this, SLOT(glowingPalettePlease(bool)));
|
|---|
| 266 | connect(zeroRotationButton, SIGNAL(toggled(bool)),
|
|---|
| 267 | this, SLOT(zeroRotationPlease(bool)));
|
|---|
| 268 | connect(minus90RotationButton, SIGNAL(toggled(bool)),
|
|---|
| 269 | this, SLOT(plus90RotationPlease(bool)));
|
|---|
| 270 | connect(plus90Rotationbutton, SIGNAL(toggled(bool)),
|
|---|
| 271 | this, SLOT(minus90RotationPlease(bool)));
|
|---|
| 272 |
|
|---|
| 273 | connect(this, SIGNAL(signalUpdateCamera()),
|
|---|
| 274 | this, SLOT(timedUpdate()));
|
|---|
| 275 | }
|
|---|
| 276 | BasicGlCamera::~BasicGlCamera()
|
|---|
| 277 | {
|
|---|
| 278 | }
|
|---|
| 279 | void BasicGlCamera::assignPixelMap(const PixelMap& map)
|
|---|
| 280 | {
|
|---|
| 281 | fPixelMap = map;
|
|---|
| 282 |
|
|---|
| 283 | for (auto it=fPixelMap.begin(); it!=fPixelMap.end(); it++)
|
|---|
| 284 | {
|
|---|
| 285 | hardwareMapping[it->index] = it->hw();
|
|---|
| 286 | softwareMapping[it->hw()] = it->index;
|
|---|
| 287 | }
|
|---|
| 288 |
|
|---|
| 289 | //now construct the correspondance between pixels and patches
|
|---|
| 290 | for (int i=0;i<NTMARK;i++)
|
|---|
| 291 | for (int j=0;j<9;j++)
|
|---|
| 292 | pixelsPatch[softwareMapping[i*9+j]] = i;
|
|---|
| 293 |
|
|---|
| 294 | for (int i=0;i<1440;i++)
|
|---|
| 295 | {
|
|---|
| 296 | for (int j=0;j<6;j++)
|
|---|
| 297 | neighbors[i][j] = -1;
|
|---|
| 298 | updateNeighbors(i);
|
|---|
| 299 | }
|
|---|
| 300 | buildPatchesIndices();
|
|---|
| 301 |
|
|---|
| 302 | }
|
|---|
| 303 | void BasicGlCamera::enableText(bool on)
|
|---|
| 304 | {
|
|---|
| 305 | fTextEnabled = on;
|
|---|
| 306 | }
|
|---|
| 307 | void BasicGlCamera::setPatchColor(int id, float color[3])
|
|---|
| 308 | {
|
|---|
| 309 | for (int i=0;i<9;i++)
|
|---|
| 310 | for (int j=0;j<3;j++)
|
|---|
| 311 | pixelsColor[softwareMapping[id*9+i]][j] = color[j];
|
|---|
| 312 | }
|
|---|
| 313 | void BasicGlCamera::setUnits(const string& units)
|
|---|
| 314 | {
|
|---|
| 315 | unitsText = units;
|
|---|
| 316 | pixelColorUpToDate = false;
|
|---|
| 317 | if (isVisible() && autoRefresh)
|
|---|
| 318 | updateGL();
|
|---|
| 319 | }
|
|---|
| 320 | void BasicGlCamera::setTitle(const string& title)
|
|---|
| 321 | {
|
|---|
| 322 | titleText = title;
|
|---|
| 323 | pixelColorUpToDate = false;
|
|---|
| 324 | if (isVisible() && autoRefresh)
|
|---|
| 325 | updateGL();
|
|---|
| 326 | }
|
|---|
| 327 | void BasicGlCamera::SetWhite(int idx)
|
|---|
| 328 | {
|
|---|
| 329 | fWhite = idx;
|
|---|
| 330 | fWhitePatch = pixelsPatch[fWhite];
|
|---|
| 331 | if (isVisible() && autoRefresh)
|
|---|
| 332 | updateGL();
|
|---|
| 333 | // CalculatePatchColor();
|
|---|
| 334 | }
|
|---|
| 335 | void BasicGlCamera::SetMin(int64_t min)
|
|---|
| 336 | {
|
|---|
| 337 | // cout << "min: " << min << endl;
|
|---|
| 338 | fMin = min;
|
|---|
| 339 | pixelColorUpToDate = false;
|
|---|
| 340 | if (isVisible() && autoRefresh)
|
|---|
| 341 | updateGL();
|
|---|
| 342 | }
|
|---|
| 343 | void BasicGlCamera::setAutoscaleLowerLimit(float val)
|
|---|
| 344 | {
|
|---|
| 345 | fScaleLimit = val;
|
|---|
| 346 | if (isVisible() && autoRefresh)
|
|---|
| 347 | updateGL();
|
|---|
| 348 | }
|
|---|
| 349 |
|
|---|
| 350 | void BasicGlCamera::SetMax(int64_t max)
|
|---|
| 351 | {
|
|---|
| 352 | // cout << "max: " << max << endl;
|
|---|
| 353 | fMax = max;
|
|---|
| 354 | pixelColorUpToDate = false;
|
|---|
| 355 | if (isVisible() && autoRefresh)
|
|---|
| 356 | updateGL();
|
|---|
| 357 | }
|
|---|
| 358 | void BasicGlCamera::linearScalePlease(bool checked)
|
|---|
| 359 | {
|
|---|
| 360 | if (!checked) return;
|
|---|
| 361 | logScale = false;
|
|---|
| 362 | pixelColorUpToDate = false;
|
|---|
| 363 | emit colorPaletteHasChanged();
|
|---|
| 364 | if (isVisible() && autoRefresh)
|
|---|
| 365 | updateGL();
|
|---|
| 366 | }
|
|---|
| 367 | void BasicGlCamera::UpdateText()
|
|---|
| 368 | {
|
|---|
| 369 | ostringstream str;
|
|---|
| 370 | float min, max, median;
|
|---|
| 371 | int ii=0;
|
|---|
| 372 | for (;ii<ACTUAL_NUM_PIXELS;ii++)
|
|---|
| 373 | {
|
|---|
| 374 | if (finite(fData[ii]))
|
|---|
| 375 | {
|
|---|
| 376 | min = max = fData[ii];
|
|---|
| 377 | break;
|
|---|
| 378 | }
|
|---|
| 379 | }
|
|---|
| 380 | double mean = 0;
|
|---|
| 381 | double rms = 0;
|
|---|
| 382 | median = 0;
|
|---|
| 383 | if (ii==ACTUAL_NUM_PIXELS)
|
|---|
| 384 | {
|
|---|
| 385 | fmin = fmax = fmean = frms = fmedian = 0;
|
|---|
| 386 | return;
|
|---|
| 387 | }
|
|---|
| 388 |
|
|---|
| 389 | vector<double> medianVec;
|
|---|
| 390 | medianVec.resize(ACTUAL_NUM_PIXELS);
|
|---|
| 391 | auto it = medianVec.begin();
|
|---|
| 392 | int numSamples = 0;
|
|---|
| 393 | for (int i=0;i<ACTUAL_NUM_PIXELS;i++)
|
|---|
| 394 | {
|
|---|
| 395 | if (!finite(fData[i]))
|
|---|
| 396 | continue;
|
|---|
| 397 | if (fData[i] < min)
|
|---|
| 398 | min = fData[i];
|
|---|
| 399 | if (fData[i] > max)
|
|---|
| 400 | max = fData[i];
|
|---|
| 401 | mean += fData[i];
|
|---|
| 402 | rms += fData[i]*fData[i];
|
|---|
| 403 | //medianSet.insert(fData[i]);
|
|---|
| 404 | *it = fData[i];
|
|---|
| 405 | it++;
|
|---|
| 406 | numSamples++;
|
|---|
| 407 | }
|
|---|
| 408 |
|
|---|
| 409 | // vector<double> medianVec;
|
|---|
| 410 | // medianVec.resize(ACTUAL_NUM_PIXELS);
|
|---|
| 411 | // int iii=0;
|
|---|
| 412 | // for (auto it=medianVec.begin(); it != medianVec.end(); it++) {
|
|---|
| 413 | // *it = fData[iii];
|
|---|
| 414 | // iii++;
|
|---|
| 415 | // }
|
|---|
| 416 | sort(medianVec.begin(), medianVec.begin()+numSamples);
|
|---|
| 417 |
|
|---|
| 418 |
|
|---|
| 419 | mean /= numSamples;
|
|---|
| 420 | rms = sqrt((rms/numSamples) - (mean * mean));
|
|---|
| 421 |
|
|---|
| 422 | // multiset<double>::iterator it = medianSet.begin();
|
|---|
| 423 | auto jt = medianVec.begin();
|
|---|
| 424 | for (int i=0;i<(numSamples/2)-1;i++)
|
|---|
| 425 | {
|
|---|
| 426 | // it++;
|
|---|
| 427 | jt++;
|
|---|
| 428 | }
|
|---|
| 429 | median = *jt;
|
|---|
| 430 | // cout << *it << " " << *jt << endl;
|
|---|
| 431 | if (numSamples%2==0){
|
|---|
| 432 | jt++;
|
|---|
| 433 | median += *jt;
|
|---|
| 434 | median /= 2;}
|
|---|
| 435 |
|
|---|
| 436 | str << "Min: " << min << endl << " Max: " << max << " Mean: " << mean << " RMS: " << rms << " Median: " << median;
|
|---|
| 437 | str << " Units: " << unitsText;
|
|---|
| 438 | dataText = str.str();
|
|---|
| 439 |
|
|---|
| 440 | fmin = min;
|
|---|
| 441 | fmax = max;
|
|---|
| 442 | fmean = mean;
|
|---|
| 443 | frms = rms;
|
|---|
| 444 | fmedian = median;
|
|---|
| 445 | }
|
|---|
| 446 | void BasicGlCamera::DrawCameraText()
|
|---|
| 447 | {
|
|---|
| 448 | if (!fTextEnabled)
|
|---|
| 449 | return;
|
|---|
| 450 | glPushMatrix();
|
|---|
| 451 | glLoadIdentity();
|
|---|
| 452 | // cout << width() << " " << height() << endl;
|
|---|
| 453 | int textSize = (int)(height()*14/600);
|
|---|
| 454 | // setFont(QFont("Times", textSize));
|
|---|
| 455 | qglColor(QColor(255,223,127));
|
|---|
| 456 | float shiftx = 0.01f;//0.55f;
|
|---|
| 457 | float shifty = 0.01f;//0.65f;
|
|---|
| 458 | renderText(-shownSizex/2.f + shiftx, 0.f, 0.f, QString(dataText.c_str()));//-shownSizey/2.f + shifty, 0.f, QString(dataText.c_str()));
|
|---|
| 459 |
|
|---|
| 460 |
|
|---|
| 461 | // int textLength = titleText.size();
|
|---|
| 462 | renderText(-shownSizex/2.f + shiftx, shownSizey/2.f - textSize*pixelSize - shifty, 0.f, QString(titleText.c_str()));
|
|---|
| 463 |
|
|---|
| 464 | glPopMatrix();
|
|---|
| 465 |
|
|---|
| 466 | // textSize = (int)(600*14/600);
|
|---|
| 467 | // setFont(QFont("Times", textSize));
|
|---|
| 468 | }
|
|---|
| 469 | void BasicGlCamera::DrawScale()
|
|---|
| 470 | {
|
|---|
| 471 | glPushMatrix();
|
|---|
| 472 | glLoadIdentity();
|
|---|
| 473 | glPushAttrib(GL_POLYGON_BIT);
|
|---|
| 474 | glShadeModel(GL_SMOOTH);
|
|---|
| 475 | glBegin(GL_QUADS);
|
|---|
| 476 | float oneX = shownSizex/2.f - shownSizex/50.f;
|
|---|
| 477 | float twoX = shownSizex/2.f;
|
|---|
| 478 | float oneY = -shownSizey/2.f;
|
|---|
| 479 | float twoY = -shownSizey/4.f;
|
|---|
| 480 | float threeY = 0;
|
|---|
| 481 | float fourY = shownSizey/4.f;
|
|---|
| 482 | float fiveY = shownSizey/2.f;
|
|---|
| 483 | glColor3f(rr[0], gg[0], bb[0]);
|
|---|
| 484 | glVertex2f(oneX, oneY);
|
|---|
| 485 | glVertex2f(twoX, oneY);
|
|---|
| 486 | glColor3f(rr[1], gg[1], bb[1]);
|
|---|
| 487 | glVertex2f(twoX, twoY);
|
|---|
| 488 | glVertex2f(oneX, twoY);
|
|---|
| 489 |
|
|---|
| 490 | glVertex2f(oneX, twoY);
|
|---|
| 491 | glVertex2f(twoX, twoY);
|
|---|
| 492 | glColor3f(rr[2], gg[2], bb[2]);
|
|---|
| 493 | glVertex2f(twoX, threeY);
|
|---|
| 494 | glVertex2f(oneX, threeY);
|
|---|
| 495 |
|
|---|
| 496 | glVertex2f(oneX, threeY);
|
|---|
| 497 | glVertex2f(twoX, threeY);
|
|---|
| 498 | glColor3f(rr[3], gg[3], bb[3]);
|
|---|
| 499 | glVertex2f(twoX, fourY);
|
|---|
| 500 | glVertex2f(oneX, fourY);
|
|---|
| 501 |
|
|---|
| 502 | glVertex2f(oneX, fourY);
|
|---|
| 503 | glVertex2f(twoX, fourY);
|
|---|
| 504 | glColor3f(rr[4], gg[4], bb[4]);
|
|---|
| 505 | glVertex2f(twoX, fiveY);
|
|---|
| 506 | glVertex2f(oneX, fiveY);
|
|---|
| 507 | float zeroX = oneX - shownSizex/50.f;
|
|---|
| 508 | float zeroY = fiveY - shownSizey/50.f;
|
|---|
| 509 | glColor3fv(tooHighValueCoulour);
|
|---|
| 510 | glVertex2f(zeroX, fiveY);
|
|---|
| 511 | glVertex2f(oneX, fiveY);
|
|---|
| 512 | glVertex2f(oneX, zeroY);
|
|---|
| 513 | glVertex2f(zeroX, zeroY);
|
|---|
| 514 | glColor3fv(tooLowValueCoulour);
|
|---|
| 515 | glVertex2f(zeroX, -fiveY);
|
|---|
| 516 | glVertex2f(oneX, -fiveY);
|
|---|
| 517 | glVertex2f(oneX, -zeroY);
|
|---|
| 518 | glVertex2f(zeroX, -zeroY);
|
|---|
| 519 | glEnd();
|
|---|
| 520 | glTranslatef(0,0,0.1f);
|
|---|
| 521 |
|
|---|
| 522 | //draw linear/log tick marks
|
|---|
| 523 | glColor3f(0.f,0.f,0.f);
|
|---|
| 524 | glBegin(GL_LINES);
|
|---|
| 525 | float value;
|
|---|
| 526 | for (int i=1;i<10;i++)
|
|---|
| 527 | {
|
|---|
| 528 | if (logScale)
|
|---|
| 529 | value = log10(i);
|
|---|
| 530 | else
|
|---|
| 531 | value = (float)(i)/10.f;
|
|---|
| 532 | float yy = -shownSizey/2.f + value*shownSizey;
|
|---|
| 533 | glVertex2f(oneX, yy);
|
|---|
| 534 | glVertex2f(twoX, yy);
|
|---|
| 535 | }
|
|---|
| 536 | glEnd();
|
|---|
| 537 | glPopAttrib();
|
|---|
| 538 | glPopMatrix();
|
|---|
| 539 | }
|
|---|
| 540 | void BasicGlCamera::logScalePlease(bool checked)
|
|---|
| 541 | {
|
|---|
| 542 | if (!checked) return;
|
|---|
| 543 | logScale = true;
|
|---|
| 544 | pixelColorUpToDate = false;
|
|---|
| 545 | emit colorPaletteHasChanged();
|
|---|
| 546 | if (isVisible() && autoRefresh)
|
|---|
| 547 | updateGL();
|
|---|
| 548 | }
|
|---|
| 549 | void BasicGlCamera::regularPalettePlease(bool checked)
|
|---|
| 550 | {
|
|---|
| 551 | if (!checked) return;
|
|---|
| 552 | ss[0] = 0; ss[1] = 0.25f; ss[2] = 0.5f; ss[3] = 0.75f; ss[4] = 1.0f;
|
|---|
| 553 | rr[0] = 0.15; rr[1] = 0; rr[2] = 0; rr[3] = 1.0f; rr[4] = 0.85f;
|
|---|
| 554 | gg[0] = 0.15; gg[1] = 0; gg[2] = 1; gg[3] = 0; gg[4] = 0.85f;
|
|---|
| 555 | bb[0] = 0.15; bb[1] = 1; bb[2] = 0; bb[3] = 0; bb[4] = 0.85f;
|
|---|
| 556 | tooHighValueCoulour[0] = 1.f;
|
|---|
| 557 | tooHighValueCoulour[1] = 1.f;
|
|---|
| 558 | tooHighValueCoulour[2] = 0.f;
|
|---|
| 559 | tooLowValueCoulour[0] = 0.f;
|
|---|
| 560 | tooLowValueCoulour[1] = 1.f;
|
|---|
| 561 | tooLowValueCoulour[2] = 1.f;
|
|---|
| 562 | pixelColorUpToDate = false;
|
|---|
| 563 |
|
|---|
| 564 | emit colorPaletteHasChanged();
|
|---|
| 565 |
|
|---|
| 566 | if (isVisible() && autoRefresh)
|
|---|
| 567 | updateGL();
|
|---|
| 568 | }
|
|---|
| 569 | void BasicGlCamera::prettyPalettePlease(bool checked)
|
|---|
| 570 | {
|
|---|
| 571 | if (!checked) return;
|
|---|
| 572 | ss[0] = 0.f; ss[1] = 0.25f; ss[2] = 0.5f; ss[3] = 0.75f; ss[4] = 1.0f;
|
|---|
| 573 | rr[0] = 0.f; rr[1] = 0.35f; rr[2] = 0.85f; rr[3] = 1.0f; rr[4] = 1.f;
|
|---|
| 574 | gg[0] = 0.f; gg[1] = 0.10f; gg[2] = 0.20f; gg[3] = 0.73f; gg[4] = 1.f;
|
|---|
| 575 | bb[0] = 0.f; bb[1] = 0.03f; bb[2] = 0.06f; bb[3] = 0.00f; bb[4] = 1.f;
|
|---|
| 576 | tooHighValueCoulour[0] = 0.f;
|
|---|
| 577 | tooHighValueCoulour[1] = 1.f;
|
|---|
| 578 | tooHighValueCoulour[2] = 0.f;
|
|---|
| 579 | tooLowValueCoulour[0] = 0.f;
|
|---|
| 580 | tooLowValueCoulour[1] = 0.f;
|
|---|
| 581 | tooLowValueCoulour[2] = 1.f;
|
|---|
| 582 | pixelColorUpToDate = false;
|
|---|
| 583 |
|
|---|
| 584 | emit colorPaletteHasChanged();
|
|---|
| 585 |
|
|---|
| 586 | if (isVisible() && autoRefresh)
|
|---|
| 587 | updateGL();
|
|---|
| 588 | }
|
|---|
| 589 | void BasicGlCamera::greyScalePalettePlease(bool checked)
|
|---|
| 590 | {
|
|---|
| 591 | if (!checked) return;
|
|---|
| 592 | ss[0] = 0; ss[1] = 0.25f; ss[2] = 0.5f; ss[3] = 0.75f; ss[4] = 1.0f;
|
|---|
| 593 | rr[0] = 0; rr[1] = 0.25f; rr[2] = 0.5f; rr[3] = 0.75f; rr[4] = 1.0f;
|
|---|
| 594 | gg[0] = 0; gg[1] = 0.25f; gg[2] = 0.5f; gg[3] = 0.75f; gg[4] = 1.0f;
|
|---|
| 595 | bb[0] = 0; bb[1] = 0.25f; bb[2] = 0.5f; bb[3] = 0.75f; bb[4] = 1.0f;
|
|---|
| 596 | tooHighValueCoulour[0] = 0.f;
|
|---|
| 597 | tooHighValueCoulour[1] = 1.f;
|
|---|
| 598 | tooHighValueCoulour[2] = 0.f;
|
|---|
| 599 | tooLowValueCoulour[0] = 0.f;
|
|---|
| 600 | tooLowValueCoulour[1] = 0.f;
|
|---|
| 601 | tooLowValueCoulour[2] = 1.f;
|
|---|
| 602 | pixelColorUpToDate = false;
|
|---|
| 603 |
|
|---|
| 604 | emit colorPaletteHasChanged();
|
|---|
| 605 |
|
|---|
| 606 | if (isVisible() && autoRefresh)
|
|---|
| 607 | updateGL();
|
|---|
| 608 | }
|
|---|
| 609 | void BasicGlCamera::glowingPalettePlease(bool checked)
|
|---|
| 610 | {
|
|---|
| 611 | if (!checked) return;
|
|---|
| 612 | ss[0] = 0; ss[1] = 0.25f; ss[2] = 0.5f; ss[3] = 0.75f; ss[4] = 1.0f;
|
|---|
| 613 | rr[0] = 0.15; rr[1] = 0.5; rr[2] = 1.f; rr[3] = 0.0f; rr[4] = 1.f;
|
|---|
| 614 | gg[0] = 0.15; gg[1] = 0.5; gg[2] = 1.f; gg[3] = 0.5f; gg[4] = 0.5f;
|
|---|
| 615 | bb[0] = 0.15; bb[1] = 0.5; bb[2] = 1; bb[3] = 1.f; bb[4] = 0.f;
|
|---|
| 616 | tooHighValueCoulour[0] = 1.f;
|
|---|
| 617 | tooHighValueCoulour[1] = 0.f;
|
|---|
| 618 | tooHighValueCoulour[2] = 0.f;
|
|---|
| 619 | tooLowValueCoulour[0] = 0.f;
|
|---|
| 620 | tooLowValueCoulour[1] = 1.f;
|
|---|
| 621 | tooLowValueCoulour[2] = 0.f;
|
|---|
| 622 | pixelColorUpToDate = false;
|
|---|
| 623 |
|
|---|
| 624 | emit colorPaletteHasChanged();
|
|---|
| 625 |
|
|---|
| 626 | if (isVisible() && autoRefresh)
|
|---|
| 627 | updateGL();
|
|---|
| 628 | }
|
|---|
| 629 | void BasicGlCamera::SetAutoRefresh(bool on)
|
|---|
| 630 | {
|
|---|
| 631 | autoRefresh = on;
|
|---|
| 632 | if (isVisible() && autoRefresh)
|
|---|
| 633 | updateGL();
|
|---|
| 634 | }
|
|---|
| 635 | void BasicGlCamera::zeroRotationPlease(bool checked)
|
|---|
| 636 | {
|
|---|
| 637 | if (!checked) return;
|
|---|
| 638 | cameraRotation = 0;
|
|---|
| 639 | pixelColorUpToDate = false;
|
|---|
| 640 | if (isVisible() && autoRefresh)
|
|---|
| 641 | updateGL();
|
|---|
| 642 | }
|
|---|
| 643 | void BasicGlCamera::plus90RotationPlease(bool checked)
|
|---|
| 644 | {
|
|---|
| 645 | if (!checked) return;
|
|---|
| 646 | cameraRotation = 90;
|
|---|
| 647 | pixelColorUpToDate = false;
|
|---|
| 648 | if (isVisible() && autoRefresh)
|
|---|
| 649 | updateGL();
|
|---|
| 650 | }
|
|---|
| 651 | void BasicGlCamera::minus90RotationPlease(bool checked)
|
|---|
| 652 | {
|
|---|
| 653 | if (!checked) return;
|
|---|
| 654 | cameraRotation = -90;
|
|---|
| 655 | pixelColorUpToDate = false;
|
|---|
| 656 | if (isVisible() && autoRefresh)
|
|---|
| 657 | updateGL();
|
|---|
| 658 | }
|
|---|
| 659 |
|
|---|
| 660 | void BasicGlCamera::initializeGL()
|
|---|
| 661 | {
|
|---|
| 662 | qglClearColor(QColor(25,25,38));
|
|---|
| 663 | glShadeModel(GL_FLAT);
|
|---|
| 664 | glEnable(GL_DEPTH_TEST);
|
|---|
| 665 | glDisable(GL_CULL_FACE);
|
|---|
| 666 |
|
|---|
| 667 | // glEnable (GL_LINE_SMOOTH);
|
|---|
| 668 | // glEnable (GL_BLEND);
|
|---|
| 669 | // glBlendFunc (GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA);
|
|---|
| 670 | // glHint (GL_LINE_SMOOTH_HINT, GL_NICEST);
|
|---|
| 671 |
|
|---|
| 672 | }
|
|---|
| 673 | void BasicGlCamera::resizeGL(int cWidth, int cHeight)
|
|---|
| 674 | {
|
|---|
| 675 | glViewport(0, 0, cWidth, cHeight);
|
|---|
| 676 | glMatrixMode(GL_PROJECTION);
|
|---|
| 677 | glLoadIdentity();
|
|---|
| 678 | GLfloat windowRatio = (float)cWidth/(float)cHeight;
|
|---|
| 679 | if (windowRatio < 1)
|
|---|
| 680 | {
|
|---|
| 681 | windowRatio = 1.0f/windowRatio;
|
|---|
| 682 | gluOrtho2D(-viewSize, viewSize, -viewSize*windowRatio, viewSize*windowRatio);
|
|---|
| 683 | pixelSize = 2*viewSize/(float)cWidth;
|
|---|
| 684 | shownSizex = 2*viewSize;
|
|---|
| 685 | shownSizey = 2*viewSize*windowRatio;
|
|---|
| 686 | }
|
|---|
| 687 | else
|
|---|
| 688 | {
|
|---|
| 689 | gluOrtho2D(-viewSize*windowRatio, viewSize*windowRatio, -viewSize, viewSize);
|
|---|
| 690 | pixelSize = 2*viewSize/(float)cHeight;
|
|---|
| 691 | shownSizex = 2*viewSize*windowRatio;
|
|---|
| 692 | shownSizey = 2*viewSize;
|
|---|
| 693 | }
|
|---|
| 694 | glMatrixMode(GL_MODELVIEW);
|
|---|
| 695 |
|
|---|
| 696 | fTextSize = (int)(cWidth*12/600); //want a sized 12 font for a window of 600 pixels width
|
|---|
| 697 | setFont(QFont("Monospace", fTextSize));
|
|---|
| 698 | }
|
|---|
| 699 | void BasicGlCamera::paintGL()
|
|---|
| 700 | {
|
|---|
| 701 | glClear(GL_COLOR_BUFFER_BIT);
|
|---|
| 702 | glLoadIdentity();
|
|---|
| 703 |
|
|---|
| 704 | glTranslatef(0,-0.44,0);
|
|---|
| 705 | glScalef(1.5, 1.5, 1.5);
|
|---|
| 706 |
|
|---|
| 707 | drawCamera(true);
|
|---|
| 708 |
|
|---|
| 709 | drawPatches();
|
|---|
| 710 | }
|
|---|
| 711 | void BasicGlCamera::toggleInterfaceDisplay()
|
|---|
| 712 | {
|
|---|
| 713 | if (linearButton->isVisible())
|
|---|
| 714 | {
|
|---|
| 715 | linearButton->hide();
|
|---|
| 716 | logButton->hide();
|
|---|
| 717 | regularPaletteButton->hide();
|
|---|
| 718 | prettyPaletteButton->hide();
|
|---|
| 719 | greyScalePaletteButton->hide();
|
|---|
| 720 | glowingPaletteButton->hide();
|
|---|
| 721 | zeroRotationButton->hide();
|
|---|
| 722 | minus90RotationButton->hide();
|
|---|
| 723 | plus90Rotationbutton->hide();
|
|---|
| 724 | scaleLabel->hide();
|
|---|
| 725 | colorPaletteLabel->hide();
|
|---|
| 726 | rotationLabel->hide();
|
|---|
| 727 | }
|
|---|
| 728 | else
|
|---|
| 729 | {
|
|---|
| 730 | linearButton->show();
|
|---|
| 731 | logButton->show();
|
|---|
| 732 | regularPaletteButton->show();
|
|---|
| 733 | prettyPaletteButton->show();
|
|---|
| 734 | greyScalePaletteButton->show();
|
|---|
| 735 | glowingPaletteButton->show();
|
|---|
| 736 | zeroRotationButton->show();
|
|---|
| 737 | minus90RotationButton->show();
|
|---|
| 738 | plus90Rotationbutton->show();
|
|---|
| 739 | scaleLabel->show();
|
|---|
| 740 | colorPaletteLabel->show();
|
|---|
| 741 | rotationLabel->show();
|
|---|
| 742 | }
|
|---|
| 743 | }
|
|---|
| 744 |
|
|---|
| 745 | void BasicGlCamera::mousePressEvent(QMouseEvent *)
|
|---|
| 746 | {
|
|---|
| 747 |
|
|---|
| 748 | }
|
|---|
| 749 | void BasicGlCamera::mouseMoveEvent(QMouseEvent *)
|
|---|
| 750 | {
|
|---|
| 751 |
|
|---|
| 752 | }
|
|---|
| 753 | void BasicGlCamera::mouseDoubleClickEvent(QMouseEvent *)
|
|---|
| 754 | {
|
|---|
| 755 |
|
|---|
| 756 | }
|
|---|
| 757 | void BasicGlCamera::timedUpdate()
|
|---|
| 758 | {
|
|---|
| 759 | if (isVisible())
|
|---|
| 760 | updateGL();
|
|---|
| 761 | }
|
|---|
| 762 | void BasicGlCamera::updateCamera()
|
|---|
| 763 | {
|
|---|
| 764 | emit signalUpdateCamera();
|
|---|
| 765 | }
|
|---|
| 766 | void BasicGlCamera::drawCamera(bool alsoWire)
|
|---|
| 767 | {
|
|---|
| 768 | // cout << "Super PaintGL" << endl;
|
|---|
| 769 | glColor3f(0.5,0.5,0.5);
|
|---|
| 770 | glLineWidth(1.0);
|
|---|
| 771 | float color;
|
|---|
| 772 |
|
|---|
| 773 | for (int i=0;i<ACTUAL_NUM_PIXELS;i++)
|
|---|
| 774 | {
|
|---|
| 775 | color = float(fData[i*fPixelStride+fcSlice]);// + ]eventData[nRoi*i + whichSlice]+(VALUES_SPAN/2))/(float)(VALUES_SPAN-1);
|
|---|
| 776 | int index = 0;
|
|---|
| 777 | while (ss[index] < color)
|
|---|
| 778 | index++;
|
|---|
| 779 | index--;
|
|---|
| 780 | float weight0 = (color-ss[index]) / (ss[index+1]-ss[index]);
|
|---|
| 781 | float weight1 = 1.0f-weight0;
|
|---|
| 782 | pixelsColor[i][0] = weight1*rr[index] + weight0*rr[index+1];
|
|---|
| 783 | pixelsColor[i][1] = weight1*gg[index] + weight0*gg[index+1];
|
|---|
| 784 | pixelsColor[i][2] = weight1*bb[index] + weight0*bb[index+1];
|
|---|
| 785 | }
|
|---|
| 786 |
|
|---|
| 787 | for (int i=0;i<ACTUAL_NUM_PIXELS;i++)
|
|---|
| 788 | {
|
|---|
| 789 | // if (i == 690 ||
|
|---|
| 790 | // i == 70)
|
|---|
| 791 | // continue;
|
|---|
| 792 | glColor3fv(pixelsColor[i]);
|
|---|
| 793 | glLoadName(i);
|
|---|
| 794 |
|
|---|
| 795 | drawHexagon(i,true);
|
|---|
| 796 |
|
|---|
| 797 | }
|
|---|
| 798 | if (!alsoWire)
|
|---|
| 799 | return;
|
|---|
| 800 | glColor3f(0.0f,0.0f,0.0f);
|
|---|
| 801 | for (int i=0;i<ACTUAL_NUM_PIXELS;i++)
|
|---|
| 802 | {
|
|---|
| 803 | // if (i == 690 ||
|
|---|
| 804 | // i == 70)
|
|---|
| 805 | // continue;
|
|---|
| 806 | drawHexagon(i, false);
|
|---|
| 807 | }
|
|---|
| 808 | }
|
|---|
| 809 | void BasicGlCamera::drawPatches()
|
|---|
| 810 | {
|
|---|
| 811 | glLineWidth(2.0f);
|
|---|
| 812 | float backupRadius = hexRadius;
|
|---|
| 813 | hexRadius *= 0.95;
|
|---|
| 814 | glColor3f(0.5f, 0.5f, 0.3f);
|
|---|
| 815 | glBegin(GL_LINES);
|
|---|
| 816 | for (int i=0;i<NTMARK;i++)
|
|---|
| 817 | {
|
|---|
| 818 | for (unsigned int j=0;j<patchesIndices[i].size();j++)
|
|---|
| 819 | {
|
|---|
| 820 | glVertex2fv(verticesList[patchesIndices[i][j].first]);
|
|---|
| 821 | glVertex2fv(verticesList[patchesIndices[i][j].second]);
|
|---|
| 822 | }
|
|---|
| 823 | }
|
|---|
| 824 | glEnd();
|
|---|
| 825 | hexRadius = backupRadius;
|
|---|
| 826 | }
|
|---|
| 827 | int BasicGlCamera::PixelAtPosition(const QPoint &cPos)
|
|---|
| 828 | {
|
|---|
| 829 | const int MaxSize = 512;
|
|---|
| 830 | GLuint buffer[MaxSize];
|
|---|
| 831 | GLint viewport[4];
|
|---|
| 832 |
|
|---|
| 833 | makeCurrent();
|
|---|
| 834 |
|
|---|
| 835 | glGetIntegerv(GL_VIEWPORT, viewport);
|
|---|
| 836 | glSelectBuffer(MaxSize, buffer);
|
|---|
| 837 | glRenderMode(GL_SELECT);
|
|---|
| 838 |
|
|---|
| 839 | glInitNames();
|
|---|
| 840 | glPushName(0);
|
|---|
| 841 |
|
|---|
| 842 | glMatrixMode(GL_PROJECTION);
|
|---|
| 843 | glPushMatrix();
|
|---|
| 844 | glLoadIdentity();
|
|---|
| 845 | GLfloat windowRatio = GLfloat(width()) / GLfloat(height());
|
|---|
| 846 | gluPickMatrix(GLdouble(cPos.x()), GLdouble(viewport[3] - cPos.y()),
|
|---|
| 847 | 1.0, 1.0, viewport);
|
|---|
| 848 |
|
|---|
| 849 | if (windowRatio < 1)
|
|---|
| 850 | {
|
|---|
| 851 | windowRatio = 1.0f/windowRatio;
|
|---|
| 852 | gluOrtho2D(-viewSize, viewSize, -viewSize*windowRatio, viewSize*windowRatio);
|
|---|
| 853 | }
|
|---|
| 854 | else
|
|---|
| 855 | {
|
|---|
| 856 | gluOrtho2D(-viewSize*windowRatio, viewSize*windowRatio, -viewSize, viewSize);
|
|---|
| 857 | }
|
|---|
| 858 |
|
|---|
| 859 | glMatrixMode(GL_MODELVIEW);
|
|---|
| 860 | drawCamera(false);
|
|---|
| 861 | glMatrixMode(GL_PROJECTION);
|
|---|
| 862 | glPopMatrix();
|
|---|
| 863 |
|
|---|
| 864 | //for some reason that I cannot understand, the push/pop matrix doesn't do the trick here... bizarre
|
|---|
| 865 | //ok, so re-do the resizeGL thing.
|
|---|
| 866 | resizeGL(width(), height());
|
|---|
| 867 |
|
|---|
| 868 | if (!glRenderMode(GL_RENDER))
|
|---|
| 869 | return -1;
|
|---|
| 870 |
|
|---|
| 871 | return buffer[3];
|
|---|
| 872 | }
|
|---|
| 873 | void BasicGlCamera::drawHexagon(int index, bool solid)
|
|---|
| 874 | {
|
|---|
| 875 | /* float minX, maxX, minY, maxY;
|
|---|
| 876 | minX = minY = 1e10;
|
|---|
| 877 | maxX = maxY = -1e10;
|
|---|
| 878 | for (int i=0;i<1438;i++)
|
|---|
| 879 | {
|
|---|
| 880 | for (int j=0;j<6;j++)
|
|---|
| 881 | {
|
|---|
| 882 | if (verticesList[verticesIndices[i][j]][0] > maxX)
|
|---|
| 883 | maxX = verticesList[verticesIndices[i][j]][0];
|
|---|
| 884 | if (verticesList[verticesIndices[i][j]][0] < minX)
|
|---|
| 885 | minX = verticesList[verticesIndices[i][j]][0];
|
|---|
| 886 | if (verticesList[verticesIndices[i][j]][1] > maxY)
|
|---|
| 887 | maxY = verticesList[verticesIndices[i][j]][1];
|
|---|
| 888 | if (verticesList[verticesIndices[i][j]][1] < minY)
|
|---|
| 889 | minY = verticesList[verticesIndices[i][j]][1];
|
|---|
| 890 | }
|
|---|
| 891 | }
|
|---|
| 892 | cout << "Min, Max X: " << minX << " " << maxX << endl;
|
|---|
| 893 | cout << "Min, Max Y: " << minY << " " << maxY << endl;
|
|---|
| 894 | exit(0);*/
|
|---|
| 895 | if (solid)
|
|---|
| 896 | glBegin(GL_POLYGON);
|
|---|
| 897 | else
|
|---|
| 898 | glBegin(GL_LINE_LOOP);
|
|---|
| 899 |
|
|---|
| 900 | glVertex2fv(verticesList[verticesIndices[index][0]]);
|
|---|
| 901 | glVertex2fv(verticesList[verticesIndices[index][1]]);
|
|---|
| 902 | glVertex2fv(verticesList[verticesIndices[index][2]]);
|
|---|
| 903 | glVertex2fv(verticesList[verticesIndices[index][3]]);
|
|---|
| 904 | glVertex2fv(verticesList[verticesIndices[index][4]]);
|
|---|
| 905 | glVertex2fv(verticesList[verticesIndices[index][5]]);
|
|---|
| 906 | if (solid)
|
|---|
| 907 | glVertex2fv(verticesList[verticesIndices[index][0]]);
|
|---|
| 908 |
|
|---|
| 909 | glEnd();
|
|---|
| 910 | }
|
|---|
| 911 |
|
|---|
| 912 | void BasicGlCamera::updateNeighbors(int currentPixel)
|
|---|
| 913 | {
|
|---|
| 914 | float squaredDistance = 0;
|
|---|
| 915 | for (int i=0;i<currentPixel;i++)
|
|---|
| 916 | {
|
|---|
| 917 | squaredDistance = (pixelsCoords[i][0] - pixelsCoords[currentPixel][0])*
|
|---|
| 918 | (pixelsCoords[i][0] - pixelsCoords[currentPixel][0]) +
|
|---|
| 919 | (pixelsCoords[i][1] - pixelsCoords[currentPixel][1])*
|
|---|
| 920 | (pixelsCoords[i][1] - pixelsCoords[currentPixel][1]);
|
|---|
| 921 | if (squaredDistance < 4*hexRadius*hexRadius*(1.0f+hexTolerance))//neighbor !
|
|---|
| 922 | {//ok, but which one ?
|
|---|
| 923 | if (fabs(pixelsCoords[i][0] - pixelsCoords[currentPixel][0]) < hexTolerance &&
|
|---|
| 924 | pixelsCoords[i][1] < pixelsCoords[currentPixel][1]){//top
|
|---|
| 925 | neighbors[i][0] = currentPixel;
|
|---|
| 926 | neighbors[currentPixel][3] = i;
|
|---|
| 927 | continue;}
|
|---|
| 928 | if (fabs(pixelsCoords[i][0] - pixelsCoords[currentPixel][0]) < hexTolerance &&
|
|---|
| 929 | pixelsCoords[i][1] > pixelsCoords[currentPixel][1]){//bottom
|
|---|
| 930 | neighbors[i][3] = currentPixel;
|
|---|
| 931 | neighbors[currentPixel][0] = i;
|
|---|
| 932 | continue;}
|
|---|
| 933 | if (pixelsCoords[i][0] > pixelsCoords[currentPixel][0] &&
|
|---|
| 934 | pixelsCoords[i][1] > pixelsCoords[currentPixel][1]){//top right
|
|---|
| 935 | neighbors[i][4] = currentPixel;
|
|---|
| 936 | neighbors[currentPixel][1] = i;
|
|---|
| 937 | continue;}
|
|---|
| 938 | if (pixelsCoords[i][0] > pixelsCoords[currentPixel][0] &&
|
|---|
| 939 | pixelsCoords[i][1] < pixelsCoords[currentPixel][1]){//bottom right
|
|---|
| 940 | neighbors[i][5] = currentPixel;
|
|---|
| 941 | neighbors[currentPixel][2] = i;
|
|---|
| 942 | continue;}
|
|---|
| 943 | if (pixelsCoords[i][0] < pixelsCoords[currentPixel][0] &&
|
|---|
| 944 | pixelsCoords[i][1] > pixelsCoords[currentPixel][1]){//top left
|
|---|
| 945 | neighbors[i][2] = currentPixel;
|
|---|
| 946 | neighbors[currentPixel][5] = i;
|
|---|
| 947 | continue;}
|
|---|
| 948 | if (pixelsCoords[i][0] < pixelsCoords[currentPixel][0] &&
|
|---|
| 949 | pixelsCoords[i][1] < pixelsCoords[currentPixel][1]){//bottom left
|
|---|
| 950 | neighbors[i][1] = currentPixel;
|
|---|
| 951 | neighbors[currentPixel][4] = i;
|
|---|
| 952 | continue;}
|
|---|
| 953 | }
|
|---|
| 954 | }
|
|---|
| 955 | }
|
|---|
| 956 | void BasicGlCamera::skipPixels(int start, int howMany)
|
|---|
| 957 | {
|
|---|
| 958 | for (int i=start;i<MAX_NUM_PIXELS-howMany;i++)
|
|---|
| 959 | {
|
|---|
| 960 | pixelsCoords[i][0] = pixelsCoords[i+howMany][0];
|
|---|
| 961 | pixelsCoords[i][1] = pixelsCoords[i+howMany][1];
|
|---|
| 962 | }
|
|---|
| 963 | }
|
|---|
| 964 | void BasicGlCamera::calculatePixelsCoords()
|
|---|
| 965 | {
|
|---|
| 966 | if (pixelsCoords[0][1] >= 0.299999 && pixelsCoords[0][1] <= 0.300001)
|
|---|
| 967 | return;
|
|---|
| 968 | pixelsCoords[0][0] = 0;
|
|---|
| 969 | pixelsCoords[0][1] = 0.3;
|
|---|
| 970 | pixelsCoords[0][2] = 0;
|
|---|
| 971 | pixelsCoords[1][0] = 0;
|
|---|
| 972 | pixelsCoords[1][1] = 0.3+2*hexRadius;
|
|---|
| 973 | pixelsCoords[1][2] = 0;
|
|---|
| 974 | neighbors[0][0] = 1;
|
|---|
| 975 | neighbors[1][3] = 0;
|
|---|
| 976 | //from which side of the previous hexagon are we coming from ?
|
|---|
| 977 | int fromSide = 3;
|
|---|
| 978 | //to which side are we heading to ?
|
|---|
| 979 | int toSide = 0;
|
|---|
| 980 | for (int i=2;i<MAX_NUM_PIXELS;i++)
|
|---|
| 981 | {
|
|---|
| 982 | toSide = fromSide-1;
|
|---|
| 983 | if (toSide < 0)
|
|---|
| 984 | toSide =5;
|
|---|
| 985 | while (neighbors[i-1][toSide] >= 0)
|
|---|
| 986 | {
|
|---|
| 987 | toSide--;
|
|---|
| 988 | if (toSide < 0)
|
|---|
| 989 | toSide = 5;
|
|---|
| 990 | }
|
|---|
| 991 | fromSide = toSide + 3;
|
|---|
| 992 | if (fromSide > 5)
|
|---|
| 993 | fromSide -= 6;
|
|---|
| 994 | //ok. now we now in which direction we're heading
|
|---|
| 995 | pixelsCoords[i][0] = pixelsCoords[i-1][0];
|
|---|
| 996 | pixelsCoords[i][1] = pixelsCoords[i-1][1];
|
|---|
| 997 | pixelsCoords[i][2] = pixelsCoords[i-1][2];
|
|---|
| 998 | switch (toSide)
|
|---|
| 999 | {
|
|---|
| 1000 | case 0:
|
|---|
| 1001 | pixelsCoords[i][1] += 2*hexRadius;
|
|---|
| 1002 | break;
|
|---|
| 1003 | case 1:
|
|---|
| 1004 | pixelsCoords[i][0] += (2*hexRadius)*sin(M_PI/3.0);
|
|---|
| 1005 | pixelsCoords[i][1] += (2*hexRadius)*cos(M_PI/3.0);
|
|---|
| 1006 | break;
|
|---|
| 1007 | case 2:
|
|---|
| 1008 | pixelsCoords[i][0] += (2*hexRadius)*sin(M_PI/3.0);
|
|---|
| 1009 | pixelsCoords[i][1] -= (2*hexRadius)*cos(M_PI/3.0);
|
|---|
| 1010 | break;
|
|---|
| 1011 | case 3:
|
|---|
| 1012 | pixelsCoords[i][1] -= 2*hexRadius;
|
|---|
| 1013 | break;
|
|---|
| 1014 | case 4:
|
|---|
| 1015 | pixelsCoords[i][0] -= (2*hexRadius)*sin(M_PI/3.0);
|
|---|
| 1016 | pixelsCoords[i][1] -= (2*hexRadius)*cos(M_PI/3.0);
|
|---|
| 1017 | break;
|
|---|
| 1018 | case 5:
|
|---|
| 1019 | pixelsCoords[i][0] -= (2*hexRadius)*sin(M_PI/3.0);
|
|---|
| 1020 | pixelsCoords[i][1] += (2*hexRadius)*cos(M_PI/3.0);
|
|---|
| 1021 | break;
|
|---|
| 1022 | };
|
|---|
| 1023 |
|
|---|
| 1024 | updateNeighbors(i);
|
|---|
| 1025 | }
|
|---|
| 1026 | //Ok. So now we've circled around all the way to MAX_NUM_PIXELS
|
|---|
| 1027 | //do the required shifts so that it matches the fact camera up to ACTUAL_NUM_PIXELS pixels
|
|---|
| 1028 | //remember the location pixels 1438 and 1439, and re-assign them later on
|
|---|
| 1029 | GLfloat backupCoords[4];
|
|---|
| 1030 | skipPixels(1200, 1);
|
|---|
| 1031 | skipPixels(1218, 3);
|
|---|
| 1032 | skipPixels(1236, 1);
|
|---|
| 1033 | skipPixels(1256, 1);
|
|---|
| 1034 | skipPixels(1274, 3);
|
|---|
| 1035 | skipPixels(1292, 3);
|
|---|
| 1036 | skipPixels(1309, 6);
|
|---|
| 1037 | skipPixels(1323, 7);
|
|---|
| 1038 | skipPixels(1337, 6);
|
|---|
| 1039 | skipPixels(1354, 6);
|
|---|
| 1040 | skipPixels(1368, 7);
|
|---|
| 1041 | //la c'est dans 1390 qu'il y a 1439
|
|---|
| 1042 | backupCoords[0] = pixelsCoords[1390][0];
|
|---|
| 1043 | backupCoords[1] = pixelsCoords[1390][1];
|
|---|
| 1044 | skipPixels(1382, 9);
|
|---|
| 1045 | skipPixels(1394, 12);
|
|---|
| 1046 | skipPixels(1402, 15);
|
|---|
| 1047 | skipPixels(1410, 12);
|
|---|
| 1048 | //la c'est dans 1422 qu'il y a 1438
|
|---|
| 1049 | backupCoords[2] = pixelsCoords[1422][0];
|
|---|
| 1050 | backupCoords[3] = pixelsCoords[1422][1];
|
|---|
| 1051 | skipPixels(1422, 12);
|
|---|
| 1052 | skipPixels(1430, 15);
|
|---|
| 1053 |
|
|---|
| 1054 | pixelsCoords[1438][0] = backupCoords[2];
|
|---|
| 1055 | pixelsCoords[1438][1] = backupCoords[3];
|
|---|
| 1056 | pixelsCoords[1439][0] = backupCoords[0];
|
|---|
| 1057 | pixelsCoords[1439][1] = backupCoords[1];
|
|---|
| 1058 | }
|
|---|
| 1059 | void BasicGlCamera::buildVerticesList()
|
|---|
| 1060 | {
|
|---|
| 1061 | numVertices = 0;
|
|---|
| 1062 | GLfloat cVertex[2];
|
|---|
| 1063 | for (int i=0;i<NPIX;i++)
|
|---|
| 1064 | {
|
|---|
| 1065 | for (int j=0;j<6;j++)
|
|---|
| 1066 | {
|
|---|
| 1067 | for (int k=0;k<2;k++)
|
|---|
| 1068 | cVertex[k] = hexcoords[j][k]*hexRadius + pixelsCoords[i][k];
|
|---|
| 1069 |
|
|---|
| 1070 | cVertex[1] -= hexRadius;
|
|---|
| 1071 |
|
|---|
| 1072 | bool found = false;
|
|---|
| 1073 | for (int k=0;k<numVertices;k++)
|
|---|
| 1074 | {
|
|---|
| 1075 | if ((cVertex[0] - verticesList[k][0])*
|
|---|
| 1076 | (cVertex[0] - verticesList[k][0]) +
|
|---|
| 1077 | (cVertex[1] - verticesList[k][1])*
|
|---|
| 1078 | (cVertex[1] - verticesList[k][1]) < hexTolerance*hexTolerance)
|
|---|
| 1079 | {
|
|---|
| 1080 | found = true;
|
|---|
| 1081 | break;
|
|---|
| 1082 | }
|
|---|
| 1083 | }
|
|---|
| 1084 | if (!found)
|
|---|
| 1085 | {
|
|---|
| 1086 | for (int k=0;k<2;k++)
|
|---|
| 1087 | verticesList[numVertices][k] = cVertex[k];
|
|---|
| 1088 | numVertices++;
|
|---|
| 1089 | }
|
|---|
| 1090 | }
|
|---|
| 1091 | }
|
|---|
| 1092 | //cout << "numVertices: " << numVertices << endl;
|
|---|
| 1093 | for (int i=0;i<NPIX;i++)
|
|---|
| 1094 | {
|
|---|
| 1095 | for (int j=0;j<6;j++)
|
|---|
| 1096 | {
|
|---|
| 1097 | for (int k=0;k<2;k++)
|
|---|
| 1098 | cVertex[k] = hexcoords[j][k]*hexRadius + pixelsCoords[i][k];
|
|---|
| 1099 |
|
|---|
| 1100 | cVertex[1] -= hexRadius;
|
|---|
| 1101 |
|
|---|
| 1102 | for (int k=0;k<numVertices;k++)
|
|---|
| 1103 | {
|
|---|
| 1104 | if ((cVertex[0] - verticesList[k][0])*
|
|---|
| 1105 | (cVertex[0] - verticesList[k][0]) +
|
|---|
| 1106 | (cVertex[1] - verticesList[k][1])*
|
|---|
| 1107 | (cVertex[1] - verticesList[k][1]) < hexTolerance*hexTolerance)
|
|---|
| 1108 | {
|
|---|
| 1109 | verticesIndices[i][j] = k;
|
|---|
| 1110 | break;
|
|---|
| 1111 | }
|
|---|
| 1112 | }
|
|---|
| 1113 | }
|
|---|
| 1114 | }
|
|---|
| 1115 | }
|
|---|
| 1116 | void BasicGlCamera::buildPatchesIndices()
|
|---|
| 1117 | {
|
|---|
| 1118 | vector<edge>::iterator it;
|
|---|
| 1119 | bool erased = false;
|
|---|
| 1120 | // patchesIndices.resize(NTMARK);
|
|---|
| 1121 | for (int i=0;i<NTMARK;i++)//for all patches
|
|---|
| 1122 | {
|
|---|
| 1123 | patchesIndices[i].clear();
|
|---|
| 1124 | for (int j=0;j<9;j++)//for all cells of the current patch
|
|---|
| 1125 | {
|
|---|
| 1126 | if (softwareMapping[i*9+j] >= ACTUAL_NUM_PIXELS)
|
|---|
| 1127 | continue;
|
|---|
| 1128 | for (int k=0;k<6;k++)//for all sides of the current cell
|
|---|
| 1129 | {
|
|---|
| 1130 | int first = k-1;
|
|---|
| 1131 | int second = k;
|
|---|
| 1132 | if (first < 0)
|
|---|
| 1133 | first = 5;
|
|---|
| 1134 | erased = false;
|
|---|
| 1135 | for (it=(patchesIndices[i]).begin(); it != (patchesIndices[i]).end(); it++)//check if this side is here already or not
|
|---|
| 1136 | {
|
|---|
| 1137 | const int idx = i*9+j;
|
|---|
| 1138 |
|
|---|
| 1139 | if ((it->first == verticesIndices[softwareMapping[idx]][first] &&
|
|---|
| 1140 | it->second == verticesIndices[softwareMapping[idx]][second]) ||
|
|---|
| 1141 | (it->first == verticesIndices[softwareMapping[idx]][second] &&
|
|---|
| 1142 | it->second == verticesIndices[softwareMapping[idx]][first]))
|
|---|
| 1143 | {
|
|---|
| 1144 | patchesIndices[i].erase(it);
|
|---|
| 1145 | erased = true;
|
|---|
| 1146 | break;
|
|---|
| 1147 | }
|
|---|
| 1148 | }
|
|---|
| 1149 | if (!erased)
|
|---|
| 1150 | {
|
|---|
| 1151 | edge temp;
|
|---|
| 1152 | temp.first = verticesIndices[softwareMapping[i*9+j]][first];
|
|---|
| 1153 | temp.second = verticesIndices[softwareMapping[i*9+j]][second];
|
|---|
| 1154 | patchesIndices[i].push_back(temp);
|
|---|
| 1155 | }
|
|---|
| 1156 | }
|
|---|
| 1157 | }
|
|---|
| 1158 | }
|
|---|
| 1159 | // for (int i=0;i<NTMARK;i++)
|
|---|
| 1160 | // {
|
|---|
| 1161 | // cout << ".....................patch " << i << " size: " << patchesIndices[i].size() << endl;
|
|---|
| 1162 | // for (unsigned int j=0;j<patchesIndices[i].size();j++)
|
|---|
| 1163 | // {
|
|---|
| 1164 | // if (patchesIndices[i][j].first < 0 || patchesIndices[i][j].first > 3013)
|
|---|
| 1165 | // cout << patchesIndices[i][j].first << " and " << patchesIndices[i][j].second << " and " << j << endl;
|
|---|
| 1166 | // }
|
|---|
| 1167 | // }
|
|---|
| 1168 | }
|
|---|
| 1169 |
|
|---|