Еще одна проблема увеличения масштаба изображения

У меня есть 2D-приложение, которое является средством просмотра изображений.
У меня есть возможность перемещаться по изображению в окне, и у меня есть базовый зум.
Теперь я хотел бы улучшить масштабирование, чтобы оно увеличивало масштаб до точки под мышью.
До сих пор я прочитал все возможное на Google, и я приблизился к чему-то работающему,
но все еще не могу найти что-то работающее, вот что я получил:

          glMatrixMode (GL_PROJECTION);
glLoadIdentity();

// w and h are respectivly the width  and height of the window
// newzoom_pos is the point in openGL coordinates  where the user requested the zoom
// zoomFactor is between [0.1,10.0]  , 0.1 -> 1.0 means downscale, 1.0 -> 10.0 means upscale
// transX and transY are used to pan around the image
float left = (0.f+transX -newzoom_pos.x())/zoomFactor  +newzoom_pos.x();
float right = (w+transX -newzoom_pos.x())/zoomFactor +newzoom_pos.x();
float bottom = (h-transY-newzoom_pos.y())/zoomFactor +newzoom_pos.y();
float top = (0.f-transY -newzoom_pos.y())/zoomFactor +newzoom_pos.y();

glOrtho(left, right, top, bottom, -1, 1);

glMatrixMode (GL_MODELVIEW);
glLoadIdentity ();

Он почти работает, но это не нормально, он не масштабируется точно до запрошенной точки, а когда изображение маленькое и еще не помещается полностью в средство просмотра (вокруг черное), то это грязно и не работает. ..

Кстати, вот код, когда пользователь увеличивает масштаб в обработчике события колеса мыши:

zoomFactor+=0.1;
QPoint oglpos = openGLpos(new_zoomPoint.x(), new_zoomPoint.y());
this->newzoom_pos = oglpos;

функция openGLpos выглядит следующим образом:

QPoint ViewerGL::openGLpos(int x,int y){
GLint viewport[4];
GLdouble modelview[16];
GLdouble projection[16];
GLfloat winX=0, winY=0, winZ=0;
GLdouble posX=0, posY=0, posZ=0;
glGetDoublev( GL_MODELVIEW_MATRIX, modelview );
glGetDoublev( GL_PROJECTION_MATRIX, projection );
glGetIntegerv( GL_VIEWPORT, viewport );
winX = (float)x;
winY = viewport[3]- y;
if(winY == 0) winY =1.f;
glReadPixels( x, winY, 1, 1, GL_DEPTH_COMPONENT, GL_FLOAT, &winZ );
GLint b=gluUnProject( winX, winY, winZ, modelview, projection, viewport, &posX, &posY, &posZ);
if(b==GLU_FALSE) cout << "failed unprojection" << endl;
return QPoint(posX,posY);
}

Кто-нибудь может мне помочь через это?