4

我有这个代码工作：

typedef model::point<double, 2, cs::spherical_equatorial<degree> > degree_point;

degree_point FlindersSE(-37.0, 144.0);

和这个：

quantity<plane_angle> Flinders = 0.375 * radians; //this works 0.375 radians

但我想做度分和秒，然后再转换成弧度。

我花了一天时间试图了解升压系统是如何工作的——例子有点少，所以我想知道是否有人可以展示一个简单的例子？

在此先感谢 8+)

编辑

//quantity<degree_base_unit> FlindersSDeg2.value(-37.0);
//quantity< angle::arcminute_base_unit> FlindersSMin = 57.0;
//quantity< angle::arcsecond_base_unit> FlindersSSec = 3.72030;

我想我需要更好地理解声明的工作原理。:)

编辑2：

非常感谢 - 也许我花了整个时间寻找使用 boost 的方法，但设施不存在！我认为这可能是因为我在这里找到了这个过时的代码http://www.boost.org/doc/libs/1_47_0/libs/geometry/doc/doxy/doxygen_input/sourcecode/doxygen_1.cpp

void example_dms()
{
/*
Extension, other coordinate system:
// Construction with degree/minute/seconds
boost::geometry::dms<boost::geometry::east> d1(4, 53, 32.5);

// Explicit conversion to double.
std::cout << d1.as_value() << std::endl;

// Conversion to string, with optional strings
std::cout << d1.get_dms(" deg ", " min ", " sec") << std::endl;

// Combination with latitude/longitude and cardinal directions
{
    using namespace boost::geometry;
    point_ll<double, boost::geometry::cs::geographic<boost::geometry::degree> >        canberra(
        latitude<>(dms<south>(35, 18, 27)),
        longitude<>(dms<east>(149, 7, 27.9)));
    std::cout << canberra << std::endl;
}
*/
}
4

3 回答 3

3

以下是我对升压单位和角度使用的一些转换方法：

double ToDegrees(const Angle & angle)
{
    return static_cast<boost::units::quantity<boost::units::degree::plane_angle>>(angle).value();
}

double ToRadians(const Angle & angle)
{
    return static_cast<boost::units::quantity<boost::units::si::plane_angle>>(angle).value();
}

这些由类型安全工厂补充：

Angle Degrees(double angleInDegrees)
{
    return angleInDegrees * boost::units::degree::degrees;
}

Angle Radians(double angleInRadians)
{
    return Angle(angleInRadians * boost::units::si::radians);
}

要捕获度数、分钟数、秒数，请将上述度数加倍替换为如下所示的转换结构：

struct DMS
{
    DMS(double value)
    {
        degrees = std::floor(value);
        double rem = (value-degrees) * 60;
        minutes = std::floor(rem);
        seconds = (rem-minutes) * 60;
    }

    operator double() const
    {
        return degrees + minutes/60 + seconds/3600;
    }

    double degrees;
    double minutes;
    double seconds;
};
于 2013-05-09T15:57:33.263 回答
0

您的第一个类型安全工厂没有返回角度。

也许你可以试试：

return Angle( angleInDegrees * boost::units::degree::degrees);
于 2017-05-24T14:08:05.163 回答
0

并确保您的类型正确，因此也可以进行正确的数字转换（即 PI/180）。使用这些函数进行测试：

///  theta is radian units  (si)
template<class Y>
bool isRad(const boost::units::quantity<si::plane_angle, Y>& theta)
{
    return true;
}
///  theta in other angular units 
template<class System, class Y>
bool isRad(const boost::units::quantity<boost::units::unit<boost::units::plane_angle_dimension, 
System>, Y>& theta)
{
    return false;
}
于 2017-05-24T17:26:24.123 回答