offset_scale_cartesian_coordinates.cuh

前往此文件的文档。

/*
 * 版权所有 (c) 2023, NVIDIA CORPORATION。
 *
 * 根据 Apache 许可证 2.0 版（“许可证”）获得许可；
 * 除非符合许可证的规定，否则您不得使用此文件。
 * 您可以在以下位置获取许可证副本：
 *
 * https://apache.ac.cn/licenses/LICENSE-2.0
 *
 * 除非适用法律要求或书面同意，否则软件
 * 在许可证下分发是基于“原样”基础的，
 * 不附带任何明示或暗示的保证或条件。
 * 有关特定语言的权限和限制，请参阅许可证。
 *
 */
 
#pragma once
 
#include <cuproj/detail/utility/cuda.hpp>
#include <cuproj/ellipsoid.hpp>
#include <cuproj/operation/operation.cuh>
#include <cuproj/projection_parameters.hpp>
 
#include <thrust/iterator/transform_iterator.h>
 
#include <algorithm>
 
namespace cuproj {


template <typename Coordinate, typename T = typename Coordinate::value_type>
class offset_scale_cartesian_coordinates : operation<Coordinate> {
 public
 CUPROJ_HOST_DEVICE offset_scale_cartesian_coordinates(projection_parameters<T> const& params)
 : a_(params.ellipsoid_.a), ra_(T{1.0} / a_), x0_(params.x0), y0_(params.y0)
  {
  }

 CUPROJ_HOST_DEVICE Coordinate operator()(Coordinate const& coord, direction dir) const
  {
 if (dir == direction::FORWARD)
 return forward(coord);
 else
 return inverse(coord);
  }
 
 private
 CUPROJ_HOST_DEVICE Coordinate forward(Coordinate const& coord) const
  {
 return coord * a_ + Coordinate{x0_, y0_};
  };

 CUPROJ_HOST_DEVICE Coordinate inverse(Coordinate const& coord) const
  {
 return (coord - Coordinate{x0_, y0_}) * ra_;
  };
 
 T a_; // 椭球体半长轴
 T ra_; // 椭球体半长轴的倒数
 T x0_; // 投影原点 x
 T y0_; // 投影原点 y
};

 
} // namespace cuproj