MATLAB中的长度缩放方向数据

Question

问题

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

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

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

我想创建一个rose plot来展示结构的传播＆＃39;取向。但是，我希望这个图也能代表结构的重要性以及它们面临的方向 - 考虑到长度。

是否有可能以某种方式在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列，包含地质结构的长度和方位角（方向）数据。

Answer 1

对于数据中的每一行，您可以根据其长度值将其复制给定次数。因此，如果您的结构长度为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

输出玫瑰图：

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

Answer 2

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

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');