WMSRasterSource.cpp

#include "WMSRasterSource.h"
 
#include "Foundation/Logging/Logger.h"
#include "Foundation/Collections/SmallVector.h"
 
namespace{
  using namespace Cogs::Core;
  Cogs::Logging::Log logger = Cogs::Logging::getLogger("WMSRasterSource");
 
  // Coordinate system WKIDs in WMS 1.3 where X,Y (Long,Lat) is switched to Y,X (Lat,Long)
    // The WKIDS are given in a list of ranges
  struct {
    unsigned a;
    unsigned b;
  } latLongCRSRanges[] = {
    { 4001, 4999 },
    {2044, 2045},   {2081, 2083},   {2085, 2086},   {2093, 2093},
    {2096, 2098},   {2105, 2132},   {2169, 2170},   {2176, 2180},
    {2193, 2193},   {2200, 2200},   {2206, 2212},   {2319, 2319},
    {2320, 2462},   {2523, 2549},   {2551, 2735},   {2738, 2758},
    {2935, 2941},   {2953, 2953},   {3006, 3030},   {3034, 3035},
    {3058, 3059},   {3068, 3068},   {3114, 3118},   {3126, 3138},
    {3300, 3301},   {3328, 3335},   {3346, 3346},   {3350, 3352},
    {3366, 3366},   {3416, 3416},   {20004, 20032}, {20064, 20092},
    {21413, 21423}, {21473, 21483}, {21896, 21899}, {22171, 22177},
    {22181, 22187}, {22191, 22197}, {25884, 25884}, {27205, 27232},
    {27391, 27398}, {27492, 27492}, {28402, 28432}, {28462, 28492},
    {30161, 30179}, {30800, 30800}, {31251, 31259}, {31275, 31279},
    {31281, 31290}, {31466, 31700}
  };
 
  bool isLatLong(TerrainProvider::WMSVersion version, const Cogs::Core::TerrainProvider::CoordSys& coordSys)
  {
    if (version != TerrainProvider::WMSVersion::v1_3_0) return false;
    if (coordSys.kind != TerrainProvider::CoordSys::Kind::EPSG) return false;
    for (auto& range : latLongCRSRanges) {
      if (range.a <= coordSys.id && coordSys.id <= range.b) return true;
    }
    return false;
  }
 
}
 
Cogs::Core::TerrainProvider::WMSRasterSource::WMSRasterSource(Context* context)
  : HTTPRasterSource(context)
{
}
 
Cogs::Core::TerrainProvider::WMSRasterSource::~WMSRasterSource()
{
}
 
bool Cogs::Core::TerrainProvider::WMSRasterSource::init(const WMSConfig& conf, std::unique_ptr<ICache>&& icache)
{
  layerDescriptions = conf.layerDescriptions;
  otherArguments = conf.otherArguments;
  version = conf.version;
  dataType = conf.dataType;
  serverIsLittleEndian = conf.serverIsLittleEndian;
  retryFailedRequests = conf.retryFailedRequests;
  return HTTPRasterSource::init(conf, std::move(icache));
}
 
void Cogs::Core::TerrainProvider::WMSRasterSource::getConfig(WMSConfig& conf) const
{
  HTTPRasterSource::getConfig(conf);
  conf.layerDescriptions = layerDescriptions;
  conf.otherArguments = otherArguments;
  conf.version = version;
  conf.dataType = dataType;
  conf.serverIsLittleEndian = serverIsLittleEndian;
  conf.retryFailedRequests = retryFailedRequests;
}
 
bool Cogs::Core::TerrainProvider::WMSRasterSource::createFetchTileUrl(std::string& url, TileId id)
{
  Collections::SmallVector<WMSLayerDescription*,4> layers;
 
  for (auto& item : layerDescriptions) {
    if (item.minLevel <= id.level && id.level <= item.maxLevel) {
      layers.push_back(&item);
    }
  }
  if (layers.size() == 0) {
    LOG_ERROR(logger, "No layer defined for level %u", id.level);
    return false;
  }
 
  double tileSize = std::exp2(-double(id.level));
 
  glm::dvec2 ta = tileSize * glm::dvec2(id.i, id.j);
  glm::dvec2 tb = tileSize * glm::dvec2(id.i + 1, id.j + 1);
 
  glm::dvec2 a = (glm::dvec2(1.0) - ta) * extent.min + ta * extent.max;
  glm::dvec2 b = (glm::dvec2(1.0) - tb) * extent.min + tb * extent.max;
 
  url.assign(Strings::getC(baseUrl));
  url.append("?VERSION=");
  switch (version) {
  case WMSVersion::Unknown:
    return false;
    break;
  case WMSVersion::v1_0_0:
    url.append("1.0.0");
    break;
  case WMSVersion::v1_1_0:
    url.append("1.1.0");
    break;
  case WMSVersion::v1_1_1:
    url.append("1.1.1");
    break;
  case WMSVersion::v1_3_0:
    url.append("1.3.0");
    break;
  default:
    assert(false && "Garbage in conf.wms.version");
    return false;
    break;
  }
  url.append("&REQUEST=GetMap");
  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));
 
  if (isLatLong(version, coordsys)) {
    url.append("&BBOX=");
    url.append(std::to_string(a.y)); url.append(","); url.append(std::to_string(a.x)); url.append(",");
    url.append(std::to_string(b.y)); url.append(","); url.append(std::to_string(b.x));
  }
  else {
    url.append("&BBOX=");
    url.append(std::to_string(a.x)); url.append(","); url.append(std::to_string(a.y)); url.append(",");
    url.append(std::to_string(b.x)); url.append(","); url.append(std::to_string(b.y));
  }
  url.append("&WIDTH=");  url.append(std::to_string(tiling.width));
  url.append("&HEIGHT="); url.append(std::to_string(tiling.height));
  url.append("&FORMAT=image/png");
  url.append("&LAYERS=");
  for (size_t i = 0, n=layers.size(); i < n; i++) {
    if (i != 0) url.append(",");
    url.append(Strings::getC(layers[i]->name));
  }
  url.append("&STYLES=");
  for (size_t i = 0, n = layers.size(); i < n; i++) {
    if (i != 0) url.append(",");
    url.append(Strings::getC(layers[i]->style));
  }
  return true;
}