How to export contour plot as shape file?
显示 更早的评论
Hi, I want to export result of plotting as shapefile. Before, I use several syntaxes, like shapewrite, but it does not work. Thank you.
[X,Y,Z] = peaks;
figure
A=contour(X,Y,Z,20);
shapewrite(A,'contour.shp')
采纳的回答
[X,Y,Z] = peaks; figure A=contour(X,Y,Z,20); % shapewrite(A,'contour.shp')
x = A(1,:) ;
y = A(2,:) ;
Data = struct([]) ; % initilaize structure
for i = 1:length(x)
Data(i).Geometry = 'Point' ;
Data(i).X=x(i) ; % x value
Data(i).Y =y(i) ; % y value
Data(i).Name = randseq(1) ; % some random attribute/ name
end
shapewrite(Data, 'myfile.shp')
5 个评论
Amra Rajuli
2018-4-19
KSSV
2018-4-19
Oh..... randseq is optional..it is not compulsory...you can skip that line completely.
Amra Rajuli
2018-4-19
Max Friedrich
2021-2-4
In your attached figure you are asking where a set of extra scattered points come from. It may come from the contour matrix A itself. For each new contour line extra information is stored. The code provided by KSSV is interpreting this information as coordinates. Following example (based on KSSV solution) solves that problem:
x_cont = A(1,:);
y_cont = A(2,:);
count_parts=1;
idx=1;
while idx<length(y_cont) % extracts each single contour
part(count_parts,:)=[idx+1 idx+y_cont(idx)];
idx=idx+y_cont(idx)+1;
count_parts=count_parts+1;
end
count_parts=count_parts-1;
Data = struct([]) ; % initilaize structure
for i = 1:count_parts
Data(i).Geometry = 'Line' ;
Data(i).X=x_cont(part(i,1):part(i,2)) ; % x value
Data(i).Y =y_cont(part(i,1):part(i,2)) ; % y value
Data(i).Z =x_cont(part(i,1)-1) ; % store the contour value
end
shapewrite(Data, 'myfile.shp')
Ozgur Kirca
2023-12-6
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Mapping Toolbox 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
