在Matlab中光栅化Delaunay三角剖分

Question

我使用DelaunayTri()中的Matlab R2010b计算了Delaunay三角剖分。现在我想将三角形​​的子集（即12200个三角形）刻录成矩阵（Raster）。有没有快速/有效的方法呢？

我尝试使用poly2mask()和inpolygon()来选择三角形内的像素，但这相当慢。

Answer 1

我实施了这个快速解决方案：

qrypts=[xgridcoords, ygridccords]; %grid x and y coordinates
triids = pointLocation(dt, qrypts); %associate each grid point to its Delaunay triangle
Lia = ismember(triids,dtsubset); %keep subset of points contained in the desired triangles (dtsubset contains the indices of desired triangles)
IM=false(size(grid)); 
IM(Lia)=1;