Surface bounded by a polygon

Question

Sarah Andrade 2017-9-26

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/358418-surface-bounded-by-a-polygon

评论： Daniel Faller 2021-8-7

I have a polygon and a surface that covers only part of this polygon (this is in 2-D x and y axis). Is there any way to delimit the surface by the polygon? I want only the part that is inside the polygon.

Thanks.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Chad Greene 2017-9-26

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/358418-surface-bounded-by-a-polygon#answer_283223

Without any sample data or illustrations it's mighty difficult to know exactly what you mean, but it sounds like inpolygon is the function you're looking for. Just feed it the x,y locations of your surface (generated by meshgrid) and it will tell you which corresponding surface locations are inside the polygon.

4 个评论
显示 2更早的评论隐藏 2更早的评论

Chad Greene 2017-9-26

在 MATLAB Online 中打开

Hi Sarah,

Thanks for providing a minimal working example. That helps tremendously. Before we go any further, I want to make a fundamental change to your code. For this kind of stuff we often think of longitudes as being like x values and latitudes as similar to y values. However, the use of variable names X and Y can imply that the spatial dimensions have been projected into meters, so it gets confusing to mix lat, long, X, and Y all on the same map with the same units. If I were faced with your problem I'd start by plotting the raw scattered data in unprojected coordinates, but letting long and lat be the horizontal and vertical axes of a plot, respectively (note, this is a change from your method):

figure
scatter(long,lat,60,variavel,'filled'); 
hold on
plot(long2,lat2)
xlabel 'longitude'
ylabel 'latitude'

I've made a habit of typically using lowercase variable names (lat, long, x, y) for scattered data, and uppwercase (LAT, LONG, X, Y) for gridded data. That's a bit more intuitive for me than letting gridded longitudes become Z and latitudes X. So I'd rewrite your gridding section like this:

% No need to change variable names, just change their cases: 
[LONG,LAT] = meshgrid(floor(min(long)):.01:ceil(max(long)),...
   floor(min(lat)):0.01:ceil(max(lat)));
VARIAVEL = griddata(long,lat,variavel,LONG,LAT,'cubic');
s = surface(LONG,LAT,VARIAVEL); 
shading interp

Now use inpolygon to find which grid cells are inside the triangle:

% These points are inside the polygon: 
in = inpolygon(LONG,LAT,long2,lat2); 
% Set points not inside (~in) the polygon to NaN: 
VARIAVEL(~in) = NaN; 
% delete the surface we just plotted: 
delete(s) 
% Plot the new surface: 
s = surface(LONG,LAT,VARIAVEL); 
shading interp