Changeset 4636 for geovis/trunk/Renderer.cpp

Timestamp:

Sep 19, 2014 12:04:45 PM (9 years ago)

Author:

ldelgass

Message:

screen coords command now takes a list of coordinates

File:

• geovis/trunk/Renderer.cpp (modified) (3 diffs)

geovis/trunk/Renderer.cpp

@@ -12 +12 @@
 #include <cstdlib>
 
+#include <limits>
 #include <set>
 
@@ -986 +987 @@
 }
 
+bool Renderer::worldToScreen(std::vector<osg::Vec3d>& coords, bool invertY)
+{
+    if (!_viewer.valid()) {
+        ERROR("No viewer");
+        return false;
+    }
+    osg::Camera *cam = _viewer->getCamera();
+    osg::Matrixd MVP = cam->getViewMatrix() * cam->getProjectionMatrix();
+    const osg::Viewport *viewport = cam->getViewport();
+    for (std::vector<osg::Vec3d>::iterator itr = coords.begin();
+         itr != coords.end(); ++itr) {
+        // Get clip coords
+        osg::Vec4d pt;
+        pt = osg::Vec4d(*itr, 1.0) * MVP;
+        // Clip
+        if (pt.x() < -pt.w() ||
+            pt.x() > pt.w() ||
+            pt.y() < -pt.w() ||
+            pt.y() > pt.w() ||
+            pt.z() < -pt.w() ||
+            pt.z() > pt.w()) {
+            // Outside frustum
+            TRACE("invalid pt: %g,%g,%g,%g", pt.x(), pt.y(), pt.z(), pt.w());
+            itr->set(std::numeric_limits<double>::quiet_NaN(),
+                     std::numeric_limits<double>::quiet_NaN(),
+                     std::numeric_limits<double>::quiet_NaN());
+            continue;
+        }
+        TRACE("clip pt: %g,%g,%g,%g", pt.x(), pt.y(), pt.z(), pt.w());
+        // Perspective divide: now NDC
+        pt /= pt.w();
+#if 1
+        itr->x() = viewport->x() + viewport->width() * 0.5 + pt.x() * viewport->width() * 0.5;
+        itr->y() = viewport->y() + viewport->height() * 0.5 + pt.y() * viewport->height() * 0.5;
+        //double near = 0;
+        //double far = 1;
+        //itr->z() = (far + near) * 0.5 + (far - near) * 0.5 * pt.z();
+        itr->z() = 0.5 + 0.5 * pt.z();
+#else
+        *itr = osg::Vec3d(pt.x(), pt.y(), pt.z()) * viewport->computeWindowMatrix();
+#endif
+        if (invertY) {
+            itr->y() = viewport->height() - itr->y();
+        }
+        TRACE("screen: %g,%g,%g", itr->x(), itr->y(), itr->z());
+    }
+    return true;
+}
+
 /**
  * \brief Map screen mouse coordinates to map coordinates
@@ -1009 +1059 @@
     osg::Vec3d world;
     if (_mapNode->getTerrain()->getWorldCoordsUnderMouse(_viewer->asView(), mouseX, mouseY, world)) {
-        map.fromWorld(_mapNode->getMapSRS(), world);
+        map.fromWorld(getMapSRS(), world);
         return true;
     }