(FEATURE) plot intersection lineation

Added a function to plot the intersection lineation for a folded foliation data
in the 3d viewer

Author: lachlangrose

LoopStructural/visualisation/model_visualisation.py

@@ -643,6 +643,24 @@ def add_data(self, feature, **kwargs):
         if interface.shape[0] > 0 and add_interface:
             self.add_points(self.model.rescale(interface[:,:3],inplace=False), name + "_interface_cp")
 
+    def add_intersection_lineation(self, feature, **kwargs):
+        name = feature.name
+        if 'name' in kwargs:
+            name = kwargs['name']
+            del kwargs['name']
+        intersection = feature.fold.foldframe.calculate_intersection_lineation(
+            feature.builder)
+        gpoints = feature.builder.interpolator.get_gradient_constraints()[:,:6]
+        npoints = feature.builder.interpolator.get_norm_constraints()[:,:6]
+        points = []
+        if gpoints.shape[0] > 0:
+            points.append(gpoints)
+        if npoints.shape[0] > 0:
+            points.append(npoints)
+        points = np.vstack(points)
+        if intersection.shape[0] > 0:
+            self.add_vector_data(self.model.rescale(points[:,:3],inplace=False), intersection, name + "_intersection")
+        
     def add_points(self, points, name, **kwargs):
         """