projection.cuh

前往此文件的文档。

/*
 * Copyright (c) 2023-2024, NVIDIA CORPORATION.
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
 * You may obtain a copy of the License at
 *
 * https://apache.ac.cn/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 */
 
#pragma once
 
#include <cuproj/detail/pipeline.cuh>
#include <cuproj/ellipsoid.hpp>
#include <cuproj/operation/operation.cuh>
#include <cuproj/projection_parameters.hpp>
 
#include <rmm/cuda_stream_view.hpp>
#include <rmm/exec_policy.hpp>
 
#include <thrust/device_vector.h>
#include <thrust/transform.h>
 
#include <iterator>
#include <type_traits>
 
namespace cuproj {


template <typename Coordinate>
using device_projection = typename detail::pipeline<Coordinate>;

template <typename Coordinate, typename T = typename Coordinate::value_type>
class projection {
 public
 projection(std::vector<operation_type> const& operations,
 projection_parameters<T> const& params,
 direction dir = direction::FORWARD)
 : params_(params), constructed_direction_(dir)
  {
 setup(operations);
  }

 device_projection<Coordinate> get_device_projection(direction dir) const
  {
 dir = (constructed_direction_ == direction::FORWARD) ? dir : reverse(dir);
 return device_projection<Coordinate>{
 params_, operations_.data().get(), operations_.size(), dir};
  }

 template <class InputCoordIter, class OutputCoordIter>
 void transform(InputCoordIter first,
 InputCoordIter last,
 OutputCoordIter result,
 direction dir,
 rmm::cuda_stream_view stream = rmm::cuda_stream_default) const
  {
 thrust::transform(rmm::exec_policy(stream), first, last, result, get_device_projection(dir));
  }
 
 private
 void setup(std::vector<operation_type> const& operations)
  {
 std::for_each(operations.begin(), operations.end(), [&](auto const& op) {
 switch (op) {
 case operation_type::TRANSVERSE_MERCATOR: {
 auto op = transverse_mercator<Coordinate>{params_};
 params_ = op.setup(params_);
 break;
        }
 // TODO: some ops don't have setup. Should we make them all have setup?
 default: break;
      }
    });
 
 operations_.resize(operations.size());
 thrust::copy(operations.begin(), operations.end(), operations_.begin());
  }
 
 thrust::device_vector<operation_type> operations_;
 projection_parameters<T> params_;
 direction constructed_direction_{direction::FORWARD};
};

 
} // namespace cuproj