matlab怎么绘制山峰,matlab画山区地形图

clear

clc

fid = fopen('A_data.txt');

data = textscan(fid,'%s');

fclose(fid);

datas = data{1};

zdata = zeros(1,19);

xdata = [];

ydata = [];

for k = 1:size(data{1},1)

if k<39

if mod(k,2)==0

xydata(k/2) = datas(k);

else

temparray = datas(k);

s = size(temparray{1},2);

tempdata = temparray{1};

temp = 0;

for m=1:s

if tempdata(m)>='0'&&tempdata(m)<='9'

cover_data = str2double(tempdata(m));

temp = temp*10+cover_data;

end

end

zdata(ceil(k/2)) = temp;  % 提出高度值

end

end

end

% xydata

% zdata

s=size(xydata,2);

xydatas=cell(1,s);

for n=1:s

temp = xydata{n};

s1 = size(temp,2);

lo = find(temp=='('|temp==','|temp==')');

temp(lo)=' ';

xydatas{1,n} = str2num(temp);

end

for n=1:s

ss = size(xydatas{1,n},2);

for nums=1:ss

if mod(nums,2)==1

xdata{n}(1,ceil(nums/2)) = xydatas{1,n}(nums);

else

ydata{n}(1,nums/2) = xydatas{1,n}(nums);

end

end

end

for n=1:s

sss = size(xdata{n},2);

zdatas{n}(1:sss) = zdata(n);

end

% 山峰坐标值

% xdata

% ydata

% zdatas

x=[];

y=[];

z=[];

for n=1:s

x = [x xdata{n}];

y = [y ydata{n}];

z = [z zdatas{n}];

end

% x

% y

% z

plot3(x,y,z,'g','LineWidth',5)

grid on

hold on

scatter3(x,y,z)