adririquelme / DSE

Software to extract discontinuity sets from rock masses' 3D point clouds
GNU General Public License v2.0
30 stars 21 forks source link

Export Discontinuity info (Dip and Dip Angle) for further analysis (Dips or Stereonet) #20

Open hrvojelukacic opened 1 year ago

hrvojelukacic commented 1 year ago

HI,

can you please tell me how I can export discontinuity orientation data (dip and dip direction) of each discontinuity plane (cluster) for the purpose of further analysis in Dips or Stereonet?

Thank you in advance!

hrvojelukacic commented 1 year ago

The main purpose is to perform kinematic analysis using discontinuity poles.

adririquelme commented 1 year ago

Uncheck the box "Fix normal vector" to let each cluster have its own normal. Complete the analysis. Then, when you export the results find the file "XXX js-c-abcd.txt". For each cluster the plane equation is given by A, B and C. Then, for each cluster transform it into Dip and Dip Direction.

If you use MATLAB, have a look at this function: vnor2vbuz_v02.m

hrvojelukacic commented 1 year ago

Thank you for your answer. I have successfully done the first step, but I have a problem with transforming the ABCD equation to DD and D. Can you please help me with how to do it outside Matlab (because I am not fluent) or give me advice on how to import "XXX js-c-abcd.txt" into Matlab to run vnor2vbuz_v02.m. Thanks for your help.

adririquelme commented 1 year ago

First, load the file: A=load('filename');

[w,b]=vnor2vbuz_v02(A(:,4),A(:,5),A(:,6));

w: dip direction b: dip

This should work

hrvojelukacic commented 1 year ago

Thank you for your help. I tried, but I keep getting multiple errors. I am sending the code in a Txt file so maybe you can please have a look. Maybe there is a syntax error or something. DSE_test_code.txt

adririquelme commented 1 year ago

Analyse the point cloud and export the results.

% load the file in a matrix: A=load('filename'); % be sure that the DSE folder was added to the path, now: [w,b]=vnor2vbuz_v02(A(:,4),A(:,5),A(:,6)) % where w: dip direction % b: dip

hrvojelukacic commented 1 year ago

Thank you very much for your help !