0

问题

我有一个描述地质结构的数据集。每个结构都有一行具有两个属性 - 它的长度和方向（0-360 度）。

在这个数据集中，有两种类型的结构。

  • 类型 1：更少的数据点，但结构在物理上更大（大length，因此更重要）。
  • 类型 2：更多数据点，但结构在物理上更小（小length，因此不那么重要）。

我想创建一个玫瑰图来显示结构方向的分布。但是，我希望这个图也能结合它们所面临的方向来代表结构的重要性——考虑到长度。

是否可以在 MATLAB 中以某种方式按长度对其进行缩放，以便在结构很大时不会低估数量较少的子集？

例子

一个数据集可能包含：

  • 10座南北向结构，长50公里。
  • 100个东西向的结构，长0.5公里。

在这种情况下，根据绝对数字，东西方人口看起来比南北人口更重要。然而，实际上，对这一人口做出贡献的成员的长度要小得多，因此结构并不重要。

代码

这是我到目前为止的代码：

load('WG_rose_data.xy') 
azimuth = WG_rose_data(:,2); 
length = WG_rose_data(:,1); 
rose(azimuth,20);

哪里WG_rose_data.xy是具有 2 列的数据文件，其中包含地质结构的长度和方位角（方向）数据。

4

2 回答 2

2

对于数据中的每一行，您可以根据其长度值将其复制给定次数。因此，如果您有一个长度为 50 的结构，它计为 50 个数据点，而长度为 1 的结构仅计为 1 个数据点。当然，你必须四舍五入你的长度，因为你只能有整数行。

这可以通过矩阵中的示例数据来实现d

% Set up example data: 10 large vertical structures, 100 small ones perpendicular
d = [repmat([0, 50], 10, 1); repmat([90, .5], 100, 1)];
% For each row, duplicate the data in column 1, according to the length in column 2
d1 = [];
for ii = 1:size(d,1)
    % make d(ii,2) = length copies of d(ii,1) = orientation
    d1(end+1:end+ceil(d(ii,2))) = d(ii,1); 
end

输出玫瑰图：

玫瑰图

您可以微调如何复制数据以实现实际数据和长度权重的理想平衡。

于 2017-03-28T13:35:27.267 回答
0

感谢您提供的所有帮助。此代码是我的最终工作版本供参考：

clear all
close all

% Input dataset
original_data = load('WG_rose_data.xy');

d = [];
%reformat azimuth
d(:,1)= original_data(:,2);
%reformat length
d(:,2)= original_data(:,1);

% For each row, duplicate the data in column 1, according to the length in column 2

d1 = [];
for a = 1:size(d,1)
    d1(end+1:end+ceil(d(a,2))) = d(a,1); 
end

%create oposite directions for rose diagram
length_d1_azi = length(d1);
d1_op_azi=zeros(1,length_d1_azi);

for i = 1:length_d1_azi

    d1_op_azi(i)=d1(i)-180;
    if d1_op_azi(i) < 1;
        d1_op_azi(i) = 360 - (d1_op_azi(i)*-1);
    end    

end

%join calculated oposites to original input
new_length = length_d1_azi*2;
all=zeros(new_length,1);

for i = 1:length_d1_azi
    all(i)=d1(i);
end

for j = length_d1_azi+1:new_length;
    all(j)=d1_op_azi(j-length_d1_azi);
end

%convert input aray into radians to plot
d1_rad=degtorad(all);

rose(d1_rad,24)
set(gca,'View',[-90 90],'YDir','reverse');
于 2017-03-28T14:29:44.900 回答