Hi Moustafa,
From the provided code snippet, it is evident that you are trying to plot the data and based on the varying values of the data, you want to determine the type of fault.
You can implement the functionalities as outlined below:
1. You can use the MATLAB 'plot' function to represent the Earth's surface, adjusting the 'LineWidth' parameter to define the line thickness as follows:
plot(x_range, zeros(size(x_range, 2)), 'Color', 'k', 'LineWidth',5);
2. To extend the faults to Earth’s surface, you can apply the ‘polyfit’ function along with the ‘roots’ function inside the for loop to compute the x-intercept of the line, then plot a line up to the x-intercept value, extending the fault line to the Earth's surface:
poly_res = polyfit(x(idx-1:idx), z(idx-1:idx), 1);
x_intercept = roots(poly_res);
plot([x(idx-1), x_intercept], [z(idx-1), 0], 'r-');
3. To detect the fault type, you can use the 'atan2d' function to avoid NaN values that occur with 'atand' function. Based on the absolute values, you can classify the fault type as follows:
dip_angles = abs(atan2d(diff(z), diff(x)));
normal_faults = dip_angles~=0 & dip_angles <= 90;
reverse_faults = dip_angles~=0 & dip_angles > 90;
scatter(x(normal_faults), z(normal_faults), 'g', 'filled');
scatter(x(reverse_faults), z(reverse_faults), 'm', 'filled');
For more information, refer to the following MathWorks Documentation:
- ‘polyfit’ function: https://www.mathworks.com/help/matlab/ref/polyfit.html
- ‘roots’ function: https://www.mathworks.com/help/matlab/ref/roots.html
- ‘atan2d’ function: https://www.mathworks.com/help/matlab/ref/double.atan2d.html
I hope this helps you!
