BF: Handle de-generate case of a surface which contains a single triangle

fsl/data/gifti.py

@@ -229,7 +229,7 @@ def loadGiftiMesh(filename):
                          'at least one pointset array'.format(filename))
 
     vertices = [ps.data for ps in pointsets]
-    indices  = triangles[0].data
+    indices  = np.atleast_2d(triangles[0].data)
 
     if len(vdata) == 0: vdata = None
     else:               vdata = prepareGiftiVertexData(vdata, filename)

fsl/data/mesh.py

@@ -687,7 +687,7 @@ def calcFaceNormals(vertices, indices):
     fnormals = np.cross((v1 - v0), (v2 - v0))
     fnormals = affine.normalise(fnormals)
 
-    return fnormals
+    return np.atleast_2d(fnormals)
 
 
 def calcVertexNormals(vertices, indices, fnormals):