function v0=loadPixels(segm, patientID)
    fs=string(segm.Alias_ID);
    for i=1:size(fs,1)
        fs(i)=erase(fs(i)," ");
    end

    regionIds=segm.Region_Id(fs==patientID);
    nr=size(regionIds,1);
    vtemp=zeros(nr,3);
    %rotate to match convention below
    vtemp(:,3)=segm.X(fs==patientID);
    vtemp(:,2)=segm.Y(fs==patientID);
    vtemp(:,1)=segm.Z(fs==patientID);

    v0=vtemp;
    for j=1:nr
        v0(regionIds(j)+1,:)=vtemp(j,:);
    end
end