From a473baa76d6748b44438f8b7dabe450f241ad4fe Mon Sep 17 00:00:00 2001
From: Paul McCarthy <pauldmccarthy@gmail.com>
Date: Tue, 23 Feb 2021 19:29:29 +0000
Subject: [PATCH] BF: Handle de-generate case of a surface which contains a
single triangle
---
fsl/data/gifti.py | 2 +-
fsl/data/mesh.py | 2 +-
2 files changed, 2 insertions(+), 2 deletions(-)
diff --git a/fsl/data/gifti.py b/fsl/data/gifti.py
index d75c73a79..077b7c4e1 100644
--- a/fsl/data/gifti.py
+++ b/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)
diff --git a/fsl/data/mesh.py b/fsl/data/mesh.py
index b03a628e2..b31e101e4 100644
--- a/fsl/data/mesh.py
+++ b/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):
--
GitLab