Анализ попадания точек внутрь полигона
Функция Inpolygon используется для анализа того, попадают ли заданные точки внутрь полигона:
- IN=inpolygon(X,Y.xv.yv) – возвращает матрицу IN того же размера, что X и Y. Каждый элемент матрицы IN принимает одно из значений – 1, 0.5 или 0 – в зависимости от того, находится ли точка с координатами (X(p,q),Y(p,q)) внутри полигона, вершины которого определяются векторами xv и yv:
- IN(p,q) = 1 – если то.чка (X(p.q),Y(p,q)) лежит внутри полигона;
- IN(p,q) = 0.5 – если точка (X(p,q),Y(p,q)) лежит на границе полигона;
- IN(p.q) = 0 – если точка (X(p.q),Y(p,q)) лежит вне полигона. Пример:
>> L = linspace(0.2*pi,8); >> yv = sin(L)`; >> xv= cos(L)`; >> x= randn(100.1); у = randn(100.1); >> IN = inpolygon(x,y,xv,yv); >> plot(xv,yv,'k',x(IN),y(IN),'r*',x(~IN),y(~IN),'b*')Построенные в этом примере массив точек и полигон представлены на рис. 17.4.
Рис. 17.4. Пример применения функции inpolygon
Точки, попавшие внутрь полигона, обозначены символом звездочки, а точки вне полигона обозначены кружками.