Why is the function isinterior filling up memory?

1 次查看(过去 30 天)
Hi all,
I have matrices of x and y coordinates of size 751, 337 and a polygon with 3700 vertices. I am trying to extract the x,y coordinates of the matrix that are comprised within my polygon. I have reshaped my x and y matrices to two vectors xV and yV of size 751*337,1 and am trying to use the function isinterior:
isinterior(polyshape, xV, yV)
The code fills up both my RAM memory and HD and never completes. I am running that code on my MBP that has 16Gb of RAM and 50Gb of free space on my disk, and the activity monitor during the computation shows 44Gb of RAM used and fills up my HD until I have about 1Gb free.
Does anyone have any idea why that occurs and, alternatively, any idea of an alternative? Before 2017b I have commonly used inpolygon and inpoly (available on the file exchange) and although it has sometime be time consuming, I never had this problem of memory. I can't rely on them because my polygon has holes in it. And yes, I know my data is quite big!
Thanks a lot!
  1 个评论
Bruno Luong
Bruno Luong 2018-10-18
编辑:Bruno Luong 2018-10-18
If you have a hole, for workaround, just do
in = inpolygon(..., xextern, yextern) & ...
~inpolygon(..., xintern, yintern)

请先登录,再进行评论。

采纳的回答

Bruno Luong
Bruno Luong 2018-10-18
编辑:Bruno Luong 2018-10-18
Alternative? Easy just go back to INPOLYGON or even better: my FEX which is probably hard to beat in term of speed.
Just another day we speek about polyshape and I told that IMO TMW neglects the essential: speed and efficiency. Here again, it shows the very point.
theta=linspace(0,2*pi,3701);
theta(end)=[];
xv = cos(theta);
yv = sin(theta);
x = 4*rand(751,337)-2;
y = 4*rand(751,337)-2;
tic
in = inpolygon(x, y, xv, yv);
toc % Elapsed time is 7.574770 seconds..
tic
in = insidepoly(x, y, xv, yv); % My FEX
toc % Elapsed time is 0.033196 seconds.
close all
hold on
plot(xv,yv,'-');
plot(x(in),y(in),'.','markersize',0.1);
axis equal
  1 个评论
Seb Biass
Seb Biass 2018-10-19
Thanks for your rapid answer! As suggested, I went back to inpoly, everything was working in the blink of an eye. I was just really curious to see why the algorithm was filling up the memory.

请先登录,再进行评论。

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Loops and Conditional Statements 的更多信息

产品


版本

R2018a

Community Treasure Hunt

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

Start Hunting!

Translated by