Retaining color when subsampling a triangulated surface: get indices from reducepatch?

You can post your problem related to MATLAB Projects here. We will try our best to help you out.
Joined:Sat Feb 18, 2017 7:12 am
Retaining color when subsampling a triangulated surface: get indices from reducepatch?

Post by Charles » Mon Mar 06, 2017 12:30 pm

I have a very densely tessellated surface which looks like this: White matter dense
This surface is too densely tessellated for me, so I subsample it to get a coarser surface. To do this, I used Matlab's reducepatch function. This works pretty well:

Code: Select all
White matter subsampled

Unfortunately, the coloring is based on a variable called sulcal_depth, which is defined for every vertex of my tessellated surface. So I need to retain sulcal depth information only from the vertices which remain after subsampling. Essentially, I need reducepatch to give me not just the subsampled version of the surface, but also the indices of vertex points that it retained. If I know the preserved indices, I can just index my sulcal_depth variable to get the new depth map.

Currently, I'm doing this as follows (this is also how I colored the subsampled version above):

Code: Select all

function indices = compute_reduced_indices(before, after)
%% Function to compute the indices of vertices preserved during an operation of
%  reducepatch. This allows you to use reducepatch to subsample a surface and
%  re-compute an original signal on the vertices for the new subsampled mesh

indices = zeros(length(after), 1);
for i = 1:length(after)
    dotprods = (before * after(i, :)') ./ sqrt(sum(before.^2, 2));
    [~, indices(i)] = max(dotprods);
But as you might imagine, this is pretty slow, because of the for loop over vertices. I don't have enough memory to vectorize the loop and compute the full dot product matrix in one go.

Is there a smart way to get reducepatch to give me indices, or an alternative approach (with or without reducepatch) that's faster?

Post Reply