Fill inside a polygon with a color

43 次查看（过去 30 天）
BN
BN 2020-5-3
评论： Patrick Laux 2021-7-13
Dear all,
I want to fill a polygon that I readed from shapefile with blue in order to visualzing it as lake. I searched alot and tried:
S2 = shaperead('lake.shp');
polygonwater1_y = S2(1).X;
polygonwater1_x = S2(1).Y;
fill (polygonwater1_y,polygonwater1_x,'b')
And:
patch (polygonwater1_y,polygonwater1_x,'b')
But as you can see below, unfortunately, it doesn't fill:
Do you know how can I do to fill it?
Thank you
  2 个评论
BN
BN 2020-5-3
Dear KSSV
I attached data here. There is my shapefile containing 3 lakes. I want to fill all of them by blue.
Thank you.

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

KSSV
KSSV 2020-5-3
编辑：KSSV 2020-5-3
It happned because, your region have multiple regions inside and they have NaN's.
shfile = "Lakes/lakes.shp" ;
S = shaperead(shfile) ;
N = length(S) ;
x = S(3).X ; y = S(3).Y ; % Picked the 3rd region
idx = find(isnan(x)) ; % find positions of NaNs
idx = [1 idx length(x)] ; % append first and last position
figure
hold on
for i = 1:length(idx)-1
pos = idx(i):idx(i+1) ; % gEt the required position
xi = x(pos) ; yi = y(pos) ; % GEt the corodinates
% Remove NaN's
xi(isnan(xi)) = [] ;
yi(isnan(yi)) = [] ;
patch(xi,yi,'b') ; % color the area
end
  4 个评论
显示 2更早的评论
KSSV
KSSV 2020-8-26
Because the user wanted the third region.

请先登录，再进行评论。

更多回答（1 个）

Steven Lord
Steven Lord 2020-8-26
Another way to plot and work with this data is to create a polyshape for each region.
S = shaperead('lakes.shp');
for k = 1:numel(S)
P(k) = polyshape(S(k).X, S(k).Y); % This will issue a warning each time
end
plot(P(3), 'FaceColor', 'r')
% Or plot all of them
figure
plot(P)
Some of the functions for working with polyshape objects may be of use to you. For example, computing the area of the polyshape agrees pretty well with the Shape_Area from the shape file.
[area(P(3)); S(3).Shape_Area; area(P(3))-S(3).Shape_Area]

类别

在 Help Center 和 File Exchange 中查找有关 Elementary Polygons 的更多信息

标签

产品


版本

R2020a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by