diff --git a/ukat/segmentation/shape_features.py b/ukat/segmentation/shape_features.py
index 004d6161..d49df8fd 100644
--- a/ukat/segmentation/shape_features.py
+++ b/ukat/segmentation/shape_features.py
@@ -65,7 +65,7 @@ def __init__(self, pixel_array, affine=None, zoom=None, kidneys=True,
         else:
             self.region_labels = np.arange(1, self.n_labels + 1)
 
-        self.features = self.get_features()
+        # self.features = self.get_features()
 
     def get_features(self):
         """
diff --git a/ukat/segmentation/tests/test_shape_features.py b/ukat/segmentation/tests/test_shape_features.py
index 74e31ef9..5ad96e13 100644
--- a/ukat/segmentation/tests/test_shape_features.py
+++ b/ukat/segmentation/tests/test_shape_features.py
@@ -41,7 +41,10 @@ def test_get_smoothed_mesh(self):
         npt.assert_allclose([vertex_stats["mean"], vertex_stats["std"],
                              vertex_stats["min"], vertex_stats["max"]],
                             expected, rtol=1e-6, atol=1e-4)
-    #
+
+    def test_init(self):
+        shape_features = ShapeFeatures(self.kidneys, self.affine)
+
     # def test_get_region_props(self):
     #     shape_features = ShapeFeatures(self.kidneys, self.affine)
     #     props_dict = shape_features._get_region_props(self.kidneys == 1)