simtk.org

Code: Select all

osimc3d = osimC3D(fileName, 1);
markers = osimc3d.getTable_markers();
forces = osimc3d.getTable_forces();

Code: Select all

function placeNaNs(self)
% Since OpenSim can not deal well with zeros beeing in the data
% instead of NaNs, this function can be used to replace them.

import org.opensim.modeling.*

hasNaNs = zeros(self.getNumTrajectories(), 1);
for iRow = 0 : self.markers.getNumRows()-1
    dataRow = self.markers.getRowAtIndex(iRow);
    rowChanged = 0;
    for iMarker = 0 : self.getNumTrajectories()-1
        data = dataRow.get(iMarker);
        if data.isNumericallyEqual(0) % Put NaNs if at [0, 0, 0]
            hasNaNs(iMarker+1) = 1;
            rowChanged = 1;
            data.setToNaN();
            dataRow.set(iMarker, data);
        end
    end
    if rowChanged
        self.markers.setRowAtIndex(iRow, dataRow);
    end
end

% Tell user what we did
for iMarker = find(hasNaNs)-1
    fprintf('Zeros were replaced by NaNs for marker %s. OpenSim might not like NaNs but can deal with it better than with zeros. \n', self.markers.getColumnLabel(iMarker));
end
end