在Matlab中光柵化Delaunay三角剖分

Question

我已經使用在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;