supercell-wx/scwx-qt/source/scwx/qt/map/radar_range_layer.cpp

#include <scwx/qt/map/radar_range_layer.hpp>
#include <scwx/qt/types/layer_types.hpp>
#include <scwx/qt/util/geographic_lib.hpp>
#include <scwx/util/logger.hpp>

#include <glm/glm.hpp>

namespace scwx::qt::map
{

static const std::string logPrefix_ = "scwx::qt::map::radar_range_layer";
static const auto        logger_    = scwx::util::Logger::Create(logPrefix_);

static std::shared_ptr<QMapLibre::Feature>
GetRangeCircle(float range, QMapLibre::Coordinate center);

void RadarRangeLayer::Add(std::shared_ptr<QMapLibre::Map> map,
                          float                           range,
                          QMapLibre::Coordinate           center,
                          const QString&                  before)
{
   static const QString layerId = QString::fromStdString(types::GetLayerName(
      types::LayerType::Data, types::DataLayer::RadarRange));

   logger_->debug("Add()");

   if (map->layerExists(layerId))
   {
      map->removeLayer(layerId);
   }
   if (map->sourceExists("rangeCircleSource"))
   {
      map->removeSource("rangeCircleSource");
   }

   std::shared_ptr<QMapLibre::Feature> rangeCircle =
      GetRangeCircle(range, center);

   map->addSource(
      "rangeCircleSource",
      {{"type", "geojson"}, {"data", QVariant::fromValue(*rangeCircle)}});
   map->addLayer(
      layerId, {{"type", "line"}, {"source", "rangeCircleSource"}}, before);
   map->setPaintProperty(layerId, "line-color", "rgba(128, 128, 128, 128)");
}

void RadarRangeLayer::Update(std::shared_ptr<QMapLibre::Map> map,
                             float                           range,
                             QMapLibre::Coordinate           center)
{
   std::shared_ptr<QMapLibre::Feature> rangeCircle =
      GetRangeCircle(range, center);

   map->updateSource("rangeCircleSource",
                     {{"data", QVariant::fromValue(*rangeCircle)}});
}

static std::shared_ptr<QMapLibre::Feature>
GetRangeCircle(float range, QMapLibre::Coordinate center)
{
   const GeographicLib::Geodesic& geodesic(
      util::GeographicLib::DefaultGeodesic());

   constexpr float angleDelta  = 0.5f;
   constexpr float angleDeltaH = angleDelta / 2.0f;

   float angle = -angleDeltaH;

   QMapLibre::Coordinates geometry;

   for (uint16_t azimuth = 0; azimuth <= 720; ++azimuth)
   {
      double latitude;
      double longitude;

      geodesic.Direct(center.first,
                      center.second,
                      angle,
                      range * 1000.0f,
                      latitude,
                      longitude);

      geometry.append({latitude, longitude});

      angle += angleDelta;
   }

   std::shared_ptr<QMapLibre::Feature> rangeCircle =
      std::make_shared<QMapLibre::Feature>(
         QMapLibre::Feature::LineStringType,
         std::initializer_list<QMapLibre::CoordinatesCollection> {
            std::initializer_list<QMapLibre::Coordinates> {geometry}});

   return rangeCircle;
}

} // namespace scwx::qt::map