Masuit.Tools/Masuit.Tools.Abstractions/Models/Sphere.cs

using System;

namespace Masuit.Tools.Models;

/// <summary>
/// 球体
/// </summary>
public class Sphere
{
    /// <summary>
    /// 地球
    /// </summary>
    public static Sphere Earth => new(6371.393);

    public Sphere(double radius)
    {
        Radius = radius;
    }

    /// <summary>
    /// 半径
    /// </summary>
    public double Radius { get; }

    /// <summary>
    /// 角度转弧度
    /// </summary>
    /// <param name="d">角度</param>
    /// <returns></returns>
    private static double Angle2Radian(double d)
    {
        return d / 180.0 * Math.PI;
    }

    /// <summary>
    /// 计算球体上两点的弧长
    /// </summary>
    /// <param name="lat1"></param>
    /// <param name="lng1"></param>
    /// <param name="lat2"></param>
    /// <param name="lng2"></param>
    /// <returns></returns>
    public double GetDistance(double lat1, double lng1, double lat2, double lng2)
    {
        var rad1 = Angle2Radian(lat1);
        var rad2 = Angle2Radian(lat2);
        return 2 * Math.Asin(Math.Sqrt(Math.Pow(Math.Sin((rad1 - rad2) / 2), 2) + Math.Cos(rad1) * Math.Cos(rad2) * Math.Pow(Math.Sin((Angle2Radian(lng1) - Angle2Radian(lng2)) / 2), 2))) * Radius;
    }

    /// <summary>
    /// 是否相交
    /// </summary>
    /// <param name="circle1"></param>
    /// <param name="circle2"></param>
    /// <returns></returns>
    public bool IsCrossWith(Circle circle1, Circle circle2)
    {
        var dis = GetDistance(circle1.Center.X, circle1.Center.Y, circle2.Center.X, circle2.Center.Y);
        return circle1.Radius - circle2.Radius < dis && dis < circle1.Radius + circle2.Radius;
    }

    /// <summary>
    /// 是否相切
    /// </summary>
    /// <param name="circle1"></param>
    /// <param name="circle2"></param>
    /// <returns></returns>
    public bool IsIntersectWith(Circle circle1, Circle circle2)
    {
        var dis = GetDistance(circle1.Center.X, circle1.Center.Y, circle2.Center.X, circle2.Center.Y);
        return Math.Abs(circle1.Radius - circle2.Radius - dis) < 1e-7 || Math.Abs(dis - (circle1.Radius + circle2.Radius)) < 1e-7;
    }

    /// <summary>
    /// 是否相离
    /// </summary>
    /// <param name="circle1"></param>
    /// <param name="circle2"></param>
    /// <returns></returns>
    public bool IsSeparateWith(Circle circle1, Circle circle2)
    {
        return !IsCrossWith(circle1, circle2) && !IsIntersectWith(circle1, circle2);
    }
}