| 1 | /*
|
|---|
| 2 | * Q3DCameraWidget.cc
|
|---|
| 3 | *
|
|---|
| 4 | * Created on: Aug 26, 2011
|
|---|
| 5 | * Author: lyard
|
|---|
| 6 | */
|
|---|
| 7 | #include "Q3DCameraWidget.h"
|
|---|
| 8 |
|
|---|
| 9 | #include <math.h>
|
|---|
| 10 | #include <sstream>
|
|---|
| 11 |
|
|---|
| 12 | #include <GL/glu.h>
|
|---|
| 13 |
|
|---|
| 14 | #include <QMouseEvent>
|
|---|
| 15 |
|
|---|
| 16 | Q3DCameraWidget::Q3DCameraWidget(QWidget* pparent) : BasicGlCamera(pparent),
|
|---|
| 17 | currentLoc()
|
|---|
| 18 | {
|
|---|
| 19 | _data.resize(432000);
|
|---|
| 20 | _colorR.resize(432000);
|
|---|
| 21 | _colorG.resize(432000);
|
|---|
| 22 | _colorB.resize(432000);
|
|---|
| 23 | _x.resize(432000);
|
|---|
| 24 | _y.resize(432000);
|
|---|
| 25 | _z.resize(432000);
|
|---|
| 26 | for (int i=0;i<432000;i++)
|
|---|
| 27 | {
|
|---|
| 28 | _data[i] = 0;
|
|---|
| 29 | _colorR[i] = 0;
|
|---|
| 30 | _colorG[i] = 0;
|
|---|
| 31 | _colorB[i] = 0;
|
|---|
| 32 | _x[i] = 0;
|
|---|
| 33 | _y[i] = 0;
|
|---|
| 34 | _z[i] = 0;
|
|---|
| 35 | }
|
|---|
| 36 | _warningWritten = false;
|
|---|
| 37 | }
|
|---|
| 38 | Q3DCameraWidget::~Q3DCameraWidget()
|
|---|
| 39 | {
|
|---|
| 40 |
|
|---|
| 41 | }
|
|---|
| 42 | void Q3DCameraWidget::timedUpdate()
|
|---|
| 43 | {
|
|---|
| 44 | updateGL();
|
|---|
| 45 | }
|
|---|
| 46 |
|
|---|
| 47 | int rotation =130;
|
|---|
| 48 | int rotationy = 30;
|
|---|
| 49 | float transZ = 0;
|
|---|
| 50 | void Q3DCameraWidget::mousePressEvent(QMouseEvent* cEvent)
|
|---|
| 51 | {
|
|---|
| 52 |
|
|---|
| 53 | if (cEvent->buttons() & Qt::LeftButton)
|
|---|
| 54 | {
|
|---|
| 55 | rotationy = -60 + (cEvent->pos().y()/(float)height())*120.f;
|
|---|
| 56 | rotation = 130 + (cEvent->pos().x()/(float)width())*180.f;
|
|---|
| 57 | }
|
|---|
| 58 | else if (cEvent->buttons() & Qt::RightButton)
|
|---|
| 59 | {
|
|---|
| 60 | if (cEvent->pos().y() > height()/2)
|
|---|
| 61 | transZ -= 0.5;
|
|---|
| 62 | else
|
|---|
| 63 | transZ += 0.5;
|
|---|
| 64 | }
|
|---|
| 65 | updateGL();
|
|---|
| 66 | }
|
|---|
| 67 | void Q3DCameraWidget::mouseDoubleClickEvent(QMouseEvent *cEvent)
|
|---|
| 68 | {
|
|---|
| 69 |
|
|---|
| 70 | }
|
|---|
| 71 | void Q3DCameraWidget::mouseMoveEvent(QMouseEvent *cEvent)
|
|---|
| 72 | {
|
|---|
| 73 | if (cEvent->buttons() & Qt::LeftButton) {
|
|---|
| 74 | mousePressEvent(cEvent);
|
|---|
| 75 | }
|
|---|
| 76 |
|
|---|
| 77 | }
|
|---|
| 78 | void Q3DCameraWidget::paintGL()
|
|---|
| 79 | {
|
|---|
| 80 | makeCurrent();
|
|---|
| 81 | glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
|
|---|
| 82 | glLoadIdentity();
|
|---|
| 83 |
|
|---|
| 84 | glTranslatef(-0.0,-0.0, -5);
|
|---|
| 85 | glTranslatef(0,0,(float)(transZ));
|
|---|
| 86 | glRotatef((float)rotationy,1.0,0.0,0.0);
|
|---|
| 87 | glRotatef((float)rotation, 0.0, 1.0, 0.0);
|
|---|
| 88 |
|
|---|
| 89 | glColor3f(1.0,0.0,0.0);
|
|---|
| 90 |
|
|---|
| 91 | glBegin(GL_TRIANGLES);
|
|---|
| 92 | for (int i=0;i<1439;i++)
|
|---|
| 93 | {
|
|---|
| 94 | for (int j=6;j<250;j++)
|
|---|
| 95 | {
|
|---|
| 96 | //get the 4 vertices that we need for drawing this patch
|
|---|
| 97 | glColor3f(_colorR[i*300+j],_colorG[i*300+j],_colorB[i*300+j]);
|
|---|
| 98 | glVertex3f(_x[i*300+j], _y[i*300+j], _z[i*300+j]);
|
|---|
| 99 | glColor3f(_colorR[i*300+j+1],_colorG[i*300+j+1],_colorB[i*300+j+1]);
|
|---|
| 100 | glVertex3f(_x[i*300+j+1], _y[i*300+j+1], _z[i*300+j+1]);
|
|---|
| 101 | glColor3f(_colorR[(i+1)*300+j],_colorG[(i+1)*300+j],_colorB[(i+1)*300+j]);
|
|---|
| 102 | glVertex3f(_x[(i+1)*300+j], _y[(i+1)*300+j], _z[(i+1)*300+j]);
|
|---|
| 103 |
|
|---|
| 104 | glColor3f(_colorR[i*300+j+1],_colorG[i*300+j+1],_colorB[i*300+j+1]);
|
|---|
| 105 | glVertex3f(_x[i*300+j+1], _y[i*300+j+1], _z[i*300+j+1]);
|
|---|
| 106 | glColor3f(_colorR[(i+1)*300+j+1],_colorG[(i+1)*300+j+1],_colorB[(i+1)*300+j+1]);
|
|---|
| 107 | glVertex3f(_x[(i+1)*300+j+1], _y[(i+1)*300+j+1], _z[(i+1)*300+j+1]);
|
|---|
| 108 | glColor3f(_colorR[(i+1)*300+j],_colorG[(i+1)*300+j],_colorB[(i+1)*300+j]);
|
|---|
| 109 | glVertex3f(_x[(i+1)*300+j], _y[(i+1)*300+j], _z[(i+1)*300+j]);
|
|---|
| 110 |
|
|---|
| 111 | }
|
|---|
| 112 | }
|
|---|
| 113 | glEnd();
|
|---|
| 114 |
|
|---|
| 115 | }
|
|---|
| 116 | void Q3DCameraWidget::calculateColorsAndPositions()
|
|---|
| 117 | {
|
|---|
| 118 | short min = 10000;
|
|---|
| 119 | short max = -10000;
|
|---|
| 120 | for (int k=0;k<1440;k++)
|
|---|
| 121 | for (int j=6;j<251;j++)
|
|---|
| 122 | {
|
|---|
| 123 | int i = k*300+j;
|
|---|
| 124 | if (_data[i] < min)
|
|---|
| 125 | min = _data[i];
|
|---|
| 126 | if (_data[i] > max)
|
|---|
| 127 | max = _data[i];
|
|---|
| 128 | }
|
|---|
| 129 | float span = max - min;
|
|---|
| 130 |
|
|---|
| 131 |
|
|---|
| 132 | //max should be at one, min at -1
|
|---|
| 133 |
|
|---|
| 134 | for (int i=0;i<1440;i++)
|
|---|
| 135 | {
|
|---|
| 136 | for (int j=6;j<251;j++)
|
|---|
| 137 | {
|
|---|
| 138 | _x[i*300+j] = -1 + (2.f*i)/1440.f;
|
|---|
| 139 | _y[i*300+j] = -0.5 + 1.0f*(_data[i*300+j] - min)/span;
|
|---|
| 140 | _z[i*300+j] = -1+(2.f*j)/300.f;
|
|---|
| 141 | float value = (_data[i*300 + j] - min)/span;
|
|---|
| 142 | if (value < 0.33)
|
|---|
| 143 | {
|
|---|
| 144 | _colorR[i*300+j] = 0;
|
|---|
| 145 | _colorG[i*300+j] = 0;
|
|---|
| 146 | _colorB[i*300+j] = value/0.33;
|
|---|
| 147 | }
|
|---|
| 148 | if (value >= 0.33 && value <= 0.66)
|
|---|
| 149 | {
|
|---|
| 150 | _colorR[i*300+j] = 0;
|
|---|
| 151 | _colorG[i*300+j] = (value-0.33)/0.33;
|
|---|
| 152 | _colorB[i*300+j] = 1 - ((value-0.33)/0.33);
|
|---|
| 153 | }
|
|---|
| 154 | if (value > 0.66)
|
|---|
| 155 | {
|
|---|
| 156 | _colorR[i*300+j] = (value-0.66)/0.33;
|
|---|
| 157 | _colorG[i*300+j] = 1 - ((value-0.66)/0.33);
|
|---|
| 158 | _colorB[i*300+j] = 0;
|
|---|
| 159 |
|
|---|
| 160 | }
|
|---|
| 161 | }
|
|---|
| 162 | }
|
|---|
| 163 |
|
|---|
| 164 |
|
|---|
| 165 |
|
|---|
| 166 | }
|
|---|
| 167 | void Q3DCameraWidget::setData(float* ddata)
|
|---|
| 168 | {
|
|---|
| 169 | if (!_warningWritten)
|
|---|
| 170 | {
|
|---|
| 171 | _warningWritten = true;
|
|---|
| 172 | cout << "Info : 3D plotter disabled. requires more work so that less than 300 slices per pixel can be loaded" << endl;
|
|---|
| 173 | cout << "Contact Etienne (etienne.lyard@unige.ch) for more information." << endl;
|
|---|
| 174 | }
|
|---|
| 175 | //disabled for now as I am working with 150 slices only
|
|---|
| 176 | /* for (int i=0;i<1440;i++)
|
|---|
| 177 | for (int j=0;j<300;j++)
|
|---|
| 178 | _data[i*300+j] = (short)(ddata[i*300 + j]);
|
|---|
| 179 | calculateColorsAndPositions();
|
|---|
| 180 | if (isVisible())
|
|---|
| 181 | updateGL();
|
|---|
| 182 | */
|
|---|
| 183 | }
|
|---|
| 184 | void Q3DCameraWidget::setData(short* ddata)
|
|---|
| 185 | {
|
|---|
| 186 | if (!_warningWritten)
|
|---|
| 187 | {
|
|---|
| 188 | _warningWritten = true;
|
|---|
| 189 | cout << "Info : 3D plotter disabled. requires more work so that less than 300 slices per pixel can be loaded" << endl;
|
|---|
| 190 | cout << "Contact Etienne (etienne.lyard@unige.ch) for more information." << endl;
|
|---|
| 191 | }
|
|---|
| 192 | /* for (int i=0;i<1440;i++)
|
|---|
| 193 | for (int j=0;j<300;j++)
|
|---|
| 194 | _data[i*300+j] = ddata[i* 300 + j];
|
|---|
| 195 | calculateColorsAndPositions();
|
|---|
| 196 | if (isVisible())
|
|---|
| 197 | updateGL();
|
|---|
| 198 | */
|
|---|
| 199 | }
|
|---|
| 200 | void Q3DCameraWidget::drawCameraBody()
|
|---|
| 201 | {
|
|---|
| 202 | glPolygonMode( GL_FRONT_AND_BACK, GL_LINE );
|
|---|
| 203 |
|
|---|
| 204 | GLfloat color[4] = {0.8f, 1.f, 1.f, 1.f};
|
|---|
| 205 | glMaterialfv(GL_FRONT_AND_BACK, GL_AMBIENT_AND_DIFFUSE, color);
|
|---|
| 206 | glMaterialf(GL_FRONT_AND_BACK, GL_SHININESS, 1.0f);
|
|---|
| 207 | gluCylinder( gluNewQuadric(),
|
|---|
| 208 | 0.62,
|
|---|
| 209 | 0.62,
|
|---|
| 210 | 1.83,
|
|---|
| 211 | 30,
|
|---|
| 212 | 2 );
|
|---|
| 213 | glPolygonMode( GL_FRONT_AND_BACK, GL_FILL );
|
|---|
| 214 |
|
|---|
| 215 |
|
|---|
| 216 |
|
|---|
| 217 | }
|
|---|
| 218 |
|
|---|
| 219 | void Q3DCameraWidget::initializeGL()
|
|---|
| 220 | {
|
|---|
| 221 | qglClearColor(QColor(25,25,38));
|
|---|
| 222 |
|
|---|
| 223 | glShadeModel(GL_SMOOTH);
|
|---|
| 224 | glEnable(GL_DEPTH_TEST);
|
|---|
| 225 | glDepthFunc(GL_LESS);
|
|---|
| 226 | glDisable(GL_LIGHTING);
|
|---|
| 227 | // glEnable(GL_LIGHTING);
|
|---|
| 228 | // glEnable(GL_LIGHT0);
|
|---|
| 229 | // glEnable(GL_AUTO_NORMAL);
|
|---|
| 230 | glDisable(GL_CULL_FACE);
|
|---|
| 231 | // glCullFace(GL_FRONT);
|
|---|
| 232 |
|
|---|
| 233 | glEnable(GL_POLYGON_SMOOTH);
|
|---|
| 234 | glEnable(GL_BLEND);
|
|---|
| 235 | glBlendFunc(GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA);
|
|---|
| 236 | glHint(GL_POLYGON_SMOOTH_HINT, GL_NICEST);
|
|---|
| 237 | }
|
|---|
| 238 | void Q3DCameraWidget::resizeGL(int cWidth, int cHeight)
|
|---|
| 239 | {
|
|---|
| 240 | glViewport(0,0,cWidth, cHeight);
|
|---|
| 241 | glMatrixMode(GL_PROJECTION);
|
|---|
| 242 | glLoadIdentity();
|
|---|
| 243 | GLfloat windowRatio = (float)cWidth/(float)cHeight;
|
|---|
| 244 | if (windowRatio < 1)
|
|---|
| 245 | {
|
|---|
| 246 | // windowRatio = 1.0f/windowRatio;
|
|---|
| 247 | gluPerspective(40.f, windowRatio, 1, 100);
|
|---|
| 248 | // gluOrtho2D(-viewSize, viewSize, -viewSize*windowRatio, viewSize*windowRatio);
|
|---|
| 249 | pixelSize = 2*viewSize/(float)cWidth;
|
|---|
| 250 | shownSizex = 2*viewSize;
|
|---|
| 251 | shownSizey = 2*viewSize*windowRatio;
|
|---|
| 252 | }
|
|---|
| 253 | else
|
|---|
| 254 | {
|
|---|
| 255 | gluPerspective(40.f, windowRatio,1, 8);
|
|---|
| 256 | // gluOrtho2D(-viewSize*windowRatio, viewSize*windowRatio, -viewSize, viewSize);
|
|---|
| 257 | pixelSize = 2*viewSize/(float)cHeight;
|
|---|
| 258 | shownSizex = 2*viewSize*windowRatio;
|
|---|
| 259 | shownSizey = 2*viewSize;
|
|---|
| 260 | }
|
|---|
| 261 | glMatrixMode(GL_MODELVIEW);
|
|---|
| 262 | }
|
|---|