From 2391b30ef3a344cf8b8df398bece2cfc5d43732f Mon Sep 17 00:00:00 2001 From: Lachlan Grose Date: Mon, 29 Apr 2024 13:29:57 +1000 Subject: [PATCH] fix: put fault normal points on the trace --- .../modelling/features/builders/_fault_builder.py | 7 ++++++- 1 file changed, 6 insertions(+), 1 deletion(-) diff --git a/LoopStructural/modelling/features/builders/_fault_builder.py b/LoopStructural/modelling/features/builders/_fault_builder.py index 51f9fead..7f61bf73 100644 --- a/LoopStructural/modelling/features/builders/_fault_builder.py +++ b/LoopStructural/modelling/features/builders/_fault_builder.py @@ -192,7 +192,12 @@ def create_data_from_geometry( fault_frame_data["coord"] == 0, fault_frame_data["val"] == 0 ) fault_center = fault_frame_data.loc[trace_mask, ["X", "Y", "Z"]].mean(axis=0).to_numpy() - + dist = np.linalg.norm( + fault_center - fault_frame_data.loc[trace_mask, ["X", "Y", "Z"]].to_numpy(), axis=1 + ) + fault_center = fault_frame_data.loc[trace_mask, ["X", "Y", "Z"]].to_numpy()[ + np.argmin(dist), : + ] self.fault_normal_vector = fault_normal_vector self.fault_slip_vector = fault_slip_vector