Method LoadRasterSource

src/extractor/raster_source.cpp:80–123

Load raster source into memory

// Load raster source into memory
int RasterContainer::LoadRasterSource(const std::string &path_string,
                                      double xmin,
                                      double xmax,
                                      double ymin,
                                      double ymax,
                                      std::size_t nrows,
                                      std::size_t ncols)
{
    const auto _xmin = static_cast<std::int32_t>(util::toFixed(util::FloatLongitude{xmin}));
    const auto _xmax = static_cast<std::int32_t>(util::toFixed(util::FloatLongitude{xmax}));
    const auto _ymin = static_cast<std::int32_t>(util::toFixed(util::FloatLatitude{ymin}));
    const auto _ymax = static_cast<std::int32_t>(util::toFixed(util::FloatLatitude{ymax}));

    const auto itr = RasterCache::getInstance().getLoadedSourcePaths().find(path_string);
    if (itr != RasterCache::getInstance().getLoadedSourcePaths().end())
    {
        util::Log() << "[source loader] Already loaded source '" << path_string << "' at source_id "
                    << itr->second;
        return itr->second;
    }

    int source_id = static_cast<int>(RasterCache::getInstance().getLoadedSources().size());

    util::Log() << "[source loader] Loading from " << path_string << "  ... ";
    TIMER_START(loading_source);

    std::filesystem::path filepath(path_string);
    if (!std::filesystem::exists(filepath))
    {
        throw util::RuntimeError(
            path_string, ErrorCode::FileOpenError, SOURCE_REF, "File not found");
    }

    RasterGrid rasterData{filepath, ncols, nrows};

    RasterSource source{std::move(rasterData), ncols, nrows, _xmin, _xmax, _ymin, _ymax};
    TIMER_STOP(loading_source);
    RasterCache::getInstance().getLoadedSourcePaths().emplace(path_string, source_id);
    RasterCache::getInstance().getLoadedSources().push_back(std::move(source));

    util::Log() << "[source loader] ok, after " << TIMER_SEC(loading_source) << "s";

    return source_id;
}

// External function for looking up nearest data point from a specified source
RasterDatum RasterContainer::GetRasterDataFromSource(unsigned int source_id, double lon, double lat)