IO Namespace Reference

Description

Read/write SANSDRaw files written by Nicos (*.001).

Functions
Datafield *	readNicosData (std::istream &input_stream)

Function Documentation

readNicosData()

Datafield * IO::readNicosData

(

std::istream &

input_stream

)

Definition at line 50 of file ReadWriteNicos.cpp.

{
    std::string line;
    int lineNumber = 0;
 
    unsigned int width = 0;
    unsigned int height = 0;
 
    // -- read dimensions
    bool inFileSection = false;
    bool fileSectionFound = false;
    bool typeFound = false;
    while (std::getline(input_stream, line)) {
        lineNumber++;
        line = BaseUtils::String::trimFront(line, " ");
 
        if (!inFileSection) {
            if (BaseUtils::String::startsWith(line, "%File")) {
                inFileSection = true;
                fileSectionFound = true;
                continue;
            }
            continue;
        }
 
        if (BaseUtils::String::startsWith(line, "%"))
            break; // next section
 
        if (BaseUtils::String::startsWith(line, "Type=")) {
            const auto parts = BaseUtils::String::split(line, "=");
            if (parts[1] != "SANSDRaw")
                throw std::runtime_error(
                    lineRelatedError("Unsupported file type '" + parts[1] + "'", lineNumber));
            typeFound = true;
        }
 
        if (BaseUtils::String::startsWith(line, "DataSizeX"))
            width = readAssignedPositiveIntValue(line, lineNumber);
        else if (BaseUtils::String::startsWith(line, "DataSizeY"))
            height = readAssignedPositiveIntValue(line, lineNumber);
        if (width != 0 && height != 0)
            break;
    }
 
    if (!fileSectionFound)
        throw std::runtime_error("Could not find 'File' section");
    if (!typeFound)
        throw std::runtime_error("File type not found");
    if (width == 0)
        throw std::runtime_error("Could not find 'DataSizeX' value");
    if (height == 0)
        throw std::runtime_error("Could not find 'DataSizeY' value");
 
    auto frame = new Frame(
        {new FixedBinAxis("x", width, 0.0, width),
         new FixedBinAxis("y", height, 0.0, height)});
    auto result = std::make_unique<Datafield>(frame);
 
    // -- read data
    bool inCountSection = false;
    bool countSectionFound = false;
    unsigned int dataRow = 0;
 
    while (std::getline(input_stream, line)) {
        lineNumber++;
        line = BaseUtils::String::trimFront(line, " ");
 
        if (!inCountSection) {
            if (BaseUtils::String::startsWith(line, "%Counts")) {
                inCountSection = true;
                countSectionFound = true;
                continue;
            }
            continue;
        }
 
        if (BaseUtils::String::startsWith(line, "%"))
            break; // next section
 
        // line is a data line
        line = BaseUtils::String::trim(line, " ");
        if (line.empty())
            continue;
 
        const auto valuesAsString = BaseUtils::String::split(line, ",");
        if (valuesAsString.size() != width)
            throw std::runtime_error(
                lineRelatedError("Number of found values (" + std::to_string(valuesAsString.size())
                                     + ") does not match DataSizeX (" + std::to_string(width) + ")",
                                 lineNumber));
 
        for (unsigned col = 0; col < width; ++col) {
            const size_t global_index = result->frame().toGlobalIndex(
                {col, height - 1 - dataRow}); // y-axis "0" is at bottom => invert y
                                              // to show first line at top of image
 
            int value = 0;
            if (!BaseUtils::String::to_int(valuesAsString[col], &value))
                throw std::runtime_error(lineRelatedError(
                    "Value '" + valuesAsString[col] + "' could not be converted to integer",
                    lineNumber));
 
            (*result)[global_index] = value;
        }
        dataRow++;
 
        if (dataRow == height)
            break;
    }
 
    if (!countSectionFound)
        throw std::runtime_error("Could not find 'Counts' section");
    if (dataRow != height)
        throw std::runtime_error("Number of found data rows (" + std::to_string(dataRow)
                                 + ") does not match DataSizeY (" + std::to_string(height) + ")");
 
    return result.release();
}

References BaseUtils::String::split(), BaseUtils::String::startsWith(), BaseUtils::String::to_int(), Frame::toGlobalIndex(), BaseUtils::String::trim(), and BaseUtils::String::trimFront().

Referenced by IOFactory::readDatafield().

Here is the call graph for this function: