Method run

plugins/matlab/filters/MatlabFilter.cpp:83–150

PointViewSet MatlabFilter::run(PointViewPtr view)
{
    log()->get(LogLevel::Debug) << "filters.matlab " << m_script <<
        " processing " << view->size() << " points." << std::endl;

    int logBufferSize(4096);
    std::unique_ptr<char[]> buf(new char[logBufferSize]);
    m_MatlabOutputBuffer.swap(buf);

    Engine* engine = mlang::Environment::get()->m_engine;
    engOutputBuffer(engine, m_MatlabOutputBuffer.get(), logBufferSize);

    Dimension::IdList dims;

    mxArray* matlabData = mlang::Script::setMatlabStruct(view, dims, m_pdalargs, m_tableMetadata, log());
    if (engPutVariable(engine, m_structName.c_str(), matlabData))
    {
        std::ostringstream oss;
        oss << "Could not push '" << m_structName << "' struct to Matlab";
        throwError(oss.str());
    }

    engEvalString(engine, m_script.m_source.c_str());

    std::string noise(m_MatlabOutputBuffer.get(), strlen(m_MatlabOutputBuffer.get()));
    log()->get(LogLevel::Debug) << "filters.matlab " << noise << std::endl;

    matlabData = engGetVariable(engine, m_structName.c_str());
    if (!matlabData)
        throwError("No 'PDAL' variable is available in Matlab scope!");

    PointViewSet viewSet;

    std::string logicalDimensionName = m_script.getLogicalMask(matlabData, log());
    if (logicalDimensionName.size())
    {
        PointViewPtr outview = view->makeNew();

        mxArray* f = mxGetField(matlabData, 0, logicalDimensionName.c_str());
        if (!f)
        {
            std::ostringstream oss;
            oss << "Unable to fetch mask dimension '" << logicalDimensionName << "'";
            throwError(oss.str());
        }

        mxLogical* logical = mxGetLogicals(f);
        if (!logical)
        {
            std::ostringstream oss;
            oss << "Unable to fetch logical mask for dimension '" << logicalDimensionName << "'";
            throwError(oss.str());
        }

        char *ok = (char *)logical;
        for (PointId idx = 0; idx < view->size(); ++idx)
            if (*ok++)
                outview->appendPoint(*view, idx);