BUG #1985: orthographic projection plot is empty

This is because proj4 sets points that are not visisble to infinity.
We set those points to 0 and hide them.

Author: danlipsa

Packages/vcs/vcs/vcs2vtk.py

@@ -4,6 +4,7 @@
 import numpy
 import json
 import os
+import math
 import meshfill
 from vtk.util import numpy_support as VN
 import cdms2
@@ -444,6 +445,20 @@ def genGrid(data1, data2, gm, deep=True, grid=None, geo=None, genVectors=False,
                                                data1.getAxis(-2))
         geo, geopts = project(pts, projection, getWrappedBounds(
             wc, [xm, xM, ym, yM], wrap))
+        # proj4 returns inf for points that are not visible. Set those to 0
+        # and hide them.
+        ghost = vg.AllocatePointGhostArray()
+        for i in range(vg.GetNumberOfPoints()):
+            point = geopts.GetPoint(i)
+            if (math.isinf(point[0]) or math.isinf(point[1])):
+                newPoint = list(point)
+                if (math.isinf(point[0])):
+                    newPoint[0] = 0
+                if (math.isinf(point[1])):
+                    newPoint[1] = 0
+                geopts.SetPoint(i, newPoint)
+                ghost.SetValue(i, vtk.vtkDataSetAttributes.HIDDENPOINT)
+
         # Sets the vertics into the grid
         vg.SetPoints(geopts)
     else:

testing/vcs/CMakeLists.txt

@@ -11,11 +11,16 @@ cdat_add_test(test_vcs_bad_png_path
   "${PYTHON_EXECUTABLE}"
   ${cdat_SOURCE_DIR}/testing/vcs/test_vcs_bad_png_path.py
 )
-cdat_add_test(test_vcs_boxfill_polar
-  "${PYTHON_EXECUTABLE}"
-  ${cdat_SOURCE_DIR}/testing/vcs/test_vcs_boxfill_polar.py
-  "${BASELINE_DIR}/test_vcs_boxfill_polar.png"
-  )
+
+foreach(projection polar mollweide lambert orthographic mercator polyconic robinson)
+  cdat_add_test(test_vcs_boxfill_${projection}
+    "${PYTHON_EXECUTABLE}"
+    ${cdat_SOURCE_DIR}/testing/vcs/test_vcs_boxfill_projection.py
+    "${BASELINE_DIR}/test_vcs_boxfill_${projection}.png"
+    ${projection}
+    )
+endforeach()
+
 cdat_add_test(test_vcs_create_get
   "${PYTHON_EXECUTABLE}"
   ${cdat_SOURCE_DIR}/testing/vcs/test_vcs_create_get.py

testing/vcs/test_vcs_boxfill_polar.py …sting/vcs/test_vcs_boxfill_projection.pytesting/vcs/test_vcs_boxfill_polar.py renamed to testing/vcs/test_vcs_boxfill_projection.py testing/vcs/test_vcs_boxfill_polar.py renamed to testing/vcs/test_vcs_boxfill_projection.py

@@ -1,16 +1,20 @@
 import os, sys, cdms2, vcs, testing.regression as regression
 
+baselineName = sys.argv[1]
+projection = sys.argv[2]
+
 
 f = cdms2.open(vcs.sample_data + "/clt.nc")
 a = f("clt")
 
 x = regression.init()
-p = x.getprojection("polar")
+p = x.getprojection(projection)
 b = x.createboxfill()
 b.projection = p
 x.plot(a(latitude=(90,-90)), b, bg=1)
 
-fileName = os.path.basename(__file__)
+fileName = os.path.basename(baselineName)
 fileName = os.path.splitext(fileName)[0]
 fileName += '.png'
-regression.run(x, fileName)
+
+regression.run(x, fileName)