求两个相交圆的交点的公式

最新推荐文章于 2023-10-31 11:19:14 发布
最新推荐文章于 2023-10-31 11:19:14 发布
阅读量2.7w 收藏 18
点赞数 5
分类专栏： 数学算法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_24378389/article/details/53736852
版权

半径为R的圆心为A，坐标为(x,y)  半径为S的圆心为B，坐标为(a,b)  两圆交点为C,D   AB与CD的交点为E，坐标为(X0,Y0)  过C点垂线与过E点水平线交点为F   令L为AB长度，K1为线AB的斜率，K2为线CD的斜率 

 则L=√[(a-x)²+(b-y)²]   K1=(b-y)/(a-x)   K2=-1/K1   CE²=R²-AE²    CE²=S²-EB²=S²-(AB-AE)²=S²-(L-AE)²=S²-L²-AE²+2L*AE    故AE=(R²-S²+L²)/2L    AE/L=(R²-S²+L²)/2L²      X0=x+(a-x)AE/L=x+(a-x)(R²-S²+L²)/(2L²)     Y0=y+K1(X0-x)      R2=CE²=R²-(X0-x)²-(Y0-y)²     R2=CF²+EF²=(K2EF)²+EF²=(1+K2²)EF²      故EF=√[R2/(1+K2²)]  


所以，C,D坐标计算公式为  Xc=X0-EF=X0-√[R2/(1+K2²)]     Yc=Y0+K2(Xc-X0) 

                                           Xd=X0+EF=X0+√[R2/(1+K2²)]   Yd=Y0+K2(Xd-X0)

TA-0-1
关注
  • 5
    点赞
  • 踩
  • 18
    收藏
    觉得还不错? 一键收藏
  • 8
    评论
专栏目录
两圆交点 坐标算法（2元2次方程）
naifen00的专栏
12-23 1728
var ax = -1;var ay = 0;var bx = 1;var by = 0;var ab = Math.sqrt(Math.pow(ax - bx, 2) + Math.pow(ay - by, 2));var ac = 1;var bc = 1; //验证数据//var ax = 4088.368;//var ay = 1992.79;
计算任意两个圆的交点
weixin_30809173的博客
12-24 1402
这里是算法的数学思想： http://mathworld.wolfram.com/Circle-CircleIntersection.html 其实懂了之后发现也很简单，中文概况是这样： 1 计算两个圆心之间的距离L2 如果 L>r1+r2 无交点 3 如果 l = r1 + r2 则一个交点， 用圆心的坐标可以计算出具体点4 假定 两个交点的连线和圆心连线的交点是D，D到第一个圆...
二维几何计算-两个圆的交点
最新发布
pipi_xia的博客
10-31 105
几何计算
求解两圆相交的交点坐标
热门推荐
JBXZ的专栏
12-11 1万+
判断并求出两个圆的交点（平面几何）
Akiyama_sou的博客
08-19 3555
参考Knut Andreas Lie的Unstructured PEBI-grids Adapting to Geological Features in Subsurface Reservoirs这篇文章，用C++实现了判断并求圆交点的函数。 ——式（1） 对于圆C1，可以看作是圆C0偏移距离d得到。根据勾股定理可以得到上述2个式子。其中h⊥d，a是C0到垂足的距离。 如果d==0或R0==R1，说明圆C0和C1的交点为空，或者是重合的。 式（1）转......
求两个圆交点的坐标-解析表达式推导及编程
JPC27的博客
07-13 2735
最开始用增量搜索法计算，这是一种暴力求解方法，不精确，不稳定。我做学生时，面对这么复杂的式子，内心是崩溃的，是完全没信心计算出来的。2+Dx+Ey+F=0，（1）-（2）求出x=a*y+b （3）算例1在CAD中验证（有舍入误差）
计算两个圆相交的交点坐标值
12-15
本范例通过扩展，实现计算地图上绘制的两个圆的相交交点的坐标值，并将其展现在地图上。
两圆交点 （公式）
一条咸鱼
01-25 9541
http://blog.csdn.net/qq_24378389/article/details/53736852 Xc=X0-EF=X0-√[R2/(1+K2²)]     Yc=Y0+K2(Xc-X0)  Xd=X0+EF=X0+√[R2/(1+K2²)]   Yd=Y0+K2(Xd-X0)
计算几何---求解两圆相交的交点坐标
mrcrack的博客
03-07 8252
摘自https://blog.csdn.net/zx3517288/article/details/53326420 方法一仔细检查了，公式没啥问题，但在实践使用时发现当y1，y2相等时，k2是求不出来的，特殊情况没考虑，建议把最后求y的式子用圆的方程联立 问：大佬们，方法一的CF=EF*K2怎么得到的呀，看不懂啊 答： 1.CF是垂直x轴，EF与x轴平行。∠FEC 为斜边与x轴正方向的夹...
已知两圆圆心坐标及半径求两圆交点 (C语言|参数方程求解)
科技浪人
09-10 7544
已知两圆圆心坐标及半径求两圆交点 (C语言|参数方程求解)在一个二维平面上给定两个圆的圆心横纵坐标、半径共6个参数, 求交点. 这个问题无非是解二元二次方程组.普通二元二次方程联立消元求解的困难在于, 中间过程里的系数会变得非常复杂, 从而导致容易出错—因为公式毕竟还是要人来推导,人的出错率比计算机要高得多得多—改用圆的参数方程求解, 可以在显著地减轻这个负担. 现在谈谈这问题的求解过程. 选
java求两圆相交坐标
09-28
使用java求两个圆相交的坐标，下载前请参见博客：http://blog.csdn.net/u013780605/article/details/52673223，文章中有具体的过程，该资源为工程文件，可直接编译运行。
lineEllipse:获取直线与椭圆/圆的交点-matlab开发
05-29
给定长轴（a），次轴（b）， 椭圆的中心 (O(h,k)) 和位于直线上的两个点 (A(x1,y1), B(x2,y2))。 这个函数给出两点交点 C1 和 C2。 如果直线和椭圆不相交，则相交点返回 NaN。 a = randperm(10,1); % 一些随机主轴...
直线的交点坐标与距离公式.doc
05-19
教学重点是根据直线的方程判断两直线的位置关系和已知两相交直线求交点。教学难点是对方程组系数的分类讨论与两直线位置关系对应情况的理解。 知识点四：直线的交点坐标与二元一次方程组的关系 通过学习，学生将...
平面直线求交算法 C++
09-05
平面直线求交算法，C++类封装
生成两直线的交点
10-25
两直线相交，公式较为简洁，注释详细，易懂
两圆交点
bad
06-19 968
计算方法 圆心距 d = (x2-x1)
计算几何之圆与圆的交点
龙进的博客
02-17 3202
计算圆与圆的交点，需要用到余弦定理 步骤如下： 求出两个圆的圆心距d 求出向量c2.c-c1.c与c1.c到某交点的向量夹角a 求出向量c2.c-c1.c与x轴的夹角t 那么，两个交点就分别是以c1.c为起点，大小为c1.r，角度为t+a、t-a的两个向量 题目：CGL_7_E AC代码： #include <iostream> #include <cstdio> #include <math.h> using namespace std; #def..
求两圆相交的交点的方法
陈鹏万里
07-25 1万+
最近研究三角定位算法，研究了下怎样求两圆的交点，现给出求解方法：
求两圆交点坐标
Sooowhat
06-11 4161
大概也可以说成是解二元二次方程组.是从...网站上copy的:The following note describes how to find the intersection point(s) between two circles on a plane, the following notation is used. The aim is to find the two point
Java 已知两个坐标经纬度以及两个坐标方向（应该是方向角度数）求两点相交得交点经纬度
05-24
这个问题需要用到一些基本的地理学知识和公式，以下是一个简单的 Java 实现： ```java public class IntersectionPoint { // 地球半径，单位为米 private static final double EARTH_RADIUS = 6378137; // 计算两点间的距离，返回单位为米 private static double distance(double lat1, double lng1, double lat2, double lng2) { double radLat1 = Math.toRadians(lat1); double radLat2 = Math.toRadians(lat2); double a = radLat1 - radLat2; double b = Math.toRadians(lng1) - Math.toRadians(lng2); double s = 2 * Math.asin(Math.sqrt(Math.pow(Math.sin(a / 2), 2) + Math.cos(radLat1) * Math.cos(radLat2) * Math.pow(Math.sin(b / 2), 2))); return s * EARTH_RADIUS; } // 计算方向角，返回值为角度数 private static double direction(double lat1, double lng1, double lat2, double lng2) { double dx = Math.cos(Math.toRadians(lat2)) * Math.sin(Math.toRadians(lng2 - lng1)); double dy = Math.cos(Math.toRadians(lat1)) * Math.sin(Math.toRadians(lat2)) - Math.sin(Math.toRadians(lat1)) * Math.cos(Math.toRadians(lat2)) * Math.cos(Math.toRadians(lng2 - lng1)); double angle = Math.atan2(dx, dy); return (Math.toDegrees(angle) + 360) % 360; } // 计算交点坐标，返回值为经纬度数组 private static double[] intersection(double lat1, double lng1, double dir1, double lat2, double lng2, double dir2) { double d = distance(lat1, lng1, lat2, lng2); double radDir1 = Math.toRadians(dir1); double radDir2 = Math.toRadians(dir2); double radLat1 = Math.toRadians(lat1); double radLat2 = Math.toRadians(lat2); double radLng1 = Math.toRadians(lng1); double radLng2 = Math.toRadians(lng2); double x = Math.sin(radDir1) * Math.cos(radDir2) * Math.sin(radLat2 - radLat1); double y = Math.sin(radDir1) * Math.cos(radDir2) * Math.cos(radLat2 - radLat1) * Math.sin(radLng2 - radLng1); double z = Math.cos(radDir1) * Math.sin(radDir2) * Math.sin(radLat2 - radLat1) * Math.cos(radLng2 - radLng1); double a = Math.pow(x, 2) + Math.pow(y, 2) + Math.pow(z, 2); double b = 2 * EARTH_RADIUS * (x * (radLng2 - radLng1) + y * (radLat2 - radLat1) + z * (Math.PI / 2 - radLat2)); double c = Math.pow(EARTH_RADIUS, 2) * Math.pow(d, 2) - Math.pow(x * EARTH_RADIUS * (radLng2 - radLng1), 2) - Math.pow(y * EARTH_RADIUS * (radLat2 - radLat1), 2) - Math.pow(z * EARTH_RADIUS * (Math.PI / 2 - radLat2), 2); double delta = Math.pow(b, 2) - 4 * a * c; if (delta < 0) { return null; } double t = (-b + Math.sqrt(delta)) / (2 * a); double lat = Math.asin(Math.sin(radLat1) * Math.cos(t) + Math.cos(radLat1) * Math.sin(t) * Math.cos(radDir1)); double lng = radLng1 + Math.atan2(Math.sin(t) * Math.sin(radDir1) * Math.sin(radLat1), Math.cos(t) - Math.sin(radLat1) * Math.sin(lat)); return new double[]{Math.toDegrees(lat), Math.toDegrees(lng)}; } public static void main(String[] args) { double lat1 = 39.9042; double lng1 = 116.4074; double dir1 = 30; double lat2 = 31.2304; double lng2 = 121.4737; double dir2 = 220; double[] intersection = intersection(lat1, lng1, dir1, lat2, lng2, dir2); if (intersection == null) { System.out.println("No intersection point."); } else { System.out.printf("Intersection point: (%f, %f)\n", intersection[0], intersection[1]); } } } ``` 以上代码假设输入的经纬度为 WGS84 坐标系下的经纬度，计算结果也是 WGS84 坐标系下的经纬度。如果需要用其他坐标系，需要进行转换。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 8
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值