WCSRasterSource.cpp

#include "WCSRasterSource.h"
 
#include "Foundation/Logging/Logger.h"
 
namespace {
  using namespace Cogs::Core;
  Cogs::Logging::Log logger = Cogs::Logging::getLogger("WMSRasterSource");
 
}
 
Cogs::Core::TerrainProvider::WCSRasterSource::WCSRasterSource(Context* context)
  : HTTPRasterSource(context)
{
}
 
Cogs::Core::TerrainProvider::WCSRasterSource::~WCSRasterSource()
{
}
 
bool Cogs::Core::TerrainProvider::WCSRasterSource::init(const WCSConfig& conf, std::unique_ptr<ICache>&& icache)
{
  layer = conf.layer;
  if (layer == NoString) {
    LOG_ERROR(logger, "No layer specified");
    return false;
  }
 
  format = conf.format;
  if (format == NoString) {
    LOG_ERROR(logger, "No format specified");
    return false;
  }
  switch (Strings::get(format).hash()) {  // FIXME: case-insensitive?
  case Cogs::hash("GeoTIFF"):   // http ://www.remotesensing.org/geotiff/geotiff.html
  case Cogs::hash("GEOTIFF"):
  //case Cogs::hash("HDF-EOS"):   //
  //case Cogs::hash("DTED"):      // http://www.nima.mil/publications/specs/printed/89020A/89020A_DTED.pdf
  //case Cogs::hash("NITF"):      // http://www.ismc.nima.mil/ntb/baseline/1999.html
  //case Cogs::hash("GML"):       // http://www.opengis.org/techno/documents/02-023r4.pdf
  //case Cogs::hash("image/png"):
    break;
  default:
    LOG_ERROR(logger, "Unrecognized format '%s'", Strings::getC(format));
    return false;
  }
 
  interpolation = conf.interpolation;
  interpolationFieldName = conf.interpolationFieldName;
  return HTTPRasterSource::init(conf, std::move(icache));
}
 
void Cogs::Core::TerrainProvider::WCSRasterSource::getConfig(WCSConfig& conf) const
{
  HTTPRasterSource::getConfig(conf);
  conf.layer = layer;
  conf.format = format;
  conf.interpolation = interpolation;
  conf.interpolationFieldName = interpolationFieldName;
}
 
 
bool Cogs::Core::TerrainProvider::WCSRasterSource::createFetchTileUrl(std::string& url, TileId id)
{
  double divisor = std::exp2(static_cast<double>(id.level));
  double xres = tiling.size.x / divisor;
  double yres = tiling.size.y / divisor;
  double west = extent.min.x + (xres * id.i);
  double east = extent.min.x + (xres * (id.i + 1));
  double south = extent.min.y + (yres * id.j);
  double north = extent.min.y + (yres * (id.j + 1));
 
  url.assign(Strings::getC(baseUrl));
  //url.append("SERVICE=WCS&VERSION=1.0.0&REQUEST=GetCoverage");
  url.append("&COVERAGE=");
  url.append(Strings::getC(layer));
  url.append("&FORMAT=");
  url.append(Strings::getC(format));  
  url.append("&CRS=");
  switch (coordsys.kind) {
  case CoordSys::Kind::EPSG:
    url.append("EPSG:");
    break;
  default:
    assert(false && "Garbage in conf.coordSys.kind");
    return false;
  }
  url.append(std::to_string(coordsys.id));
  url.append("&WIDTH=");
  url.append(std::to_string(tiling.width));
  url.append("&HEIGHT=");
  url.append(std::to_string(tiling.height));
  url.append("&BBOX=");
  url.append(std::to_string(west)); url.append(",");
  url.append(std::to_string(south)); url.append(",");
  url.append(std::to_string(east)); url.append(",");
  url.append(std::to_string(north));
  if (interpolation != NoString && interpolationFieldName != NoString) {
    url.append("&");
    url.append(Strings::getC(interpolationFieldName));
    url.append("=");
    url.append(Strings::getC(interpolation));
  }
  LOG_DEBUG(logger, "%s", url.c_str());
  return true;
}