圆和椭圆外投影

1 圆外一点在缩放到圆上

请添加图片描述

圆方程:
x 2 + y 2 = 2 2 x^2+y^2 = 2^2 x2+y2=22

直线方程:
y = k x y = kx y=kx

圆外一点:
A ( 3 , 3 ) A(3,3) A(3,3)

求点B.

方法1-解方程

圆外一点 A ( 3 , 3 ) A(3,3) A(3,3),那么:直线 k = 1 k=1 k=1,直线方程: y = x y=x y=x

方程联立:
x 2 + x 2 = 4 x^2 + x^2 = 4 \\ x2+x2=4

x = 2 y = 2 x = \sqrt{2} \\ y = \sqrt{2} \\ x=2 y=2

方法2-向量
/**
 * 将圆外一点沿着圆心缩放到圆上,XOY平面
 * @param {number} circleRadius 圆半径
 * @param {Cesium.Cartesian3} point 园外点坐标
 */
function scaleToCircle(circleRadius=1, point=new Cesium.Cartesian3()) {
  const positionX = point.x;
  const positionY = point.y;

  // 将圆外一点带入圆方程
  const x2 = positionX * positionX;
  const y2 = positionY * positionY;

  // 圆外一点指向圆心方向向量的模长
  const squaredNorm = x2 + y2;
  const norm = Math.sqrt(squaredNorm);

  const ratio = circleRadius/norm;
  return Cesium.Cartesian3.multiplyByScalar(
    point,
    ratio,
    new Cesium.Cartesian3()
  );
}

const pointA = new Cesium.Cartesian3(3, 3, 0)
const pointB = scaleToCircle(2, cartesian1)
console.error(pointB);
2 椭圆外一点在椭圆上的投影

请添加图片描述

椭圆方程:
x 2 a 2 + y 2 b 2 = 1 \frac{x^2}{a^2} + \frac{y^2}{b^2} = 1 a2x2+b2y2=1

直线方程:
y = k x y = kx y=kx

圆外一点:
A ( 3 , 3 ) A(3,3) A(3,3)

求点B.

 /**
   *
   * @param {*} a
   * @param {*} b
   * @param {*} point
   * @returns
   */
  scaleToEllipsoid(a, b, point = new Cesium.Cartesian3()) {
    const positionX = point.x;
    const positionY = point.y;

    const x2 = positionX * positionX * (1 / (a * a));
    const y2 = positionY * positionY * (1 / (b * b));

    const squaredNorm = x2 + y2;
    const ratio = Math.sqrt(1/squaredNorm);

    return Cesium.Cartesian3.multiplyByScalar(
      point,
      ratio,
      new Cesium.Cartesian3()
    );
  }

证明:
椭圆外一点:
A ( c , d ) = n ⃗ A(c,d) = \vec{n} A(c,d)=n

椭圆方程:
x 2 a 2 + y 2 b 2 = 1 \frac{x^2}{a^2}+\frac{y^2}{b^2} = 1 a2x2+b2y2=1

step1:
将椭圆外一点带入椭圆方程:
c 2 a 2 + d 2 b 2 = S \frac{c^2}{a^2}+\frac{d^2}{b^2} = S a2c2+b2d2=S

step2:
方向向量为:
1 S n ⃗ = ( c c 2 a 2 + d 2 b 2 , d c 2 a 2 + d 2 b 2 , ) \frac{1}{\sqrt{S}} \vec{n} = ( \frac{c}{\sqrt{\frac{c^2}{a^2}+\frac{d^2}{b^2}}}, \frac{d}{\sqrt{\frac{c^2}{a^2}+\frac{d^2}{b^2}}}, ) S 1n =(a2c2+b2d2 c,a2c2+b2d2 d,)

step3:
将方向向量带入椭圆方程:
c c 2 a 2 + d 2 b 2 a 2 + d c 2 a 2 + d 2 b 2 b 2 = ? c 2 a 2 ( c 2 a 2 + d 2 b 2 ) + d 2 b 2 ( d 2 a 2 + d 2 b 2 ) = ? c 2 c 2 + a 2 d 2 b 2 + d 2 b 2 a 2 a 2 d 2 = ? b 2 c 2 b 2 c 2 + a 2 d 2 + a 2 d 2 b 2 c 2 + a 2 d 2 = 1 \begin{aligned} \frac{\frac{c}{\sqrt{\frac{c^2}{a^2}+\frac{d^2}{b^2}}}}{a^2} + \frac{\frac{d}{\sqrt{\frac{c^2}{a^2}+\frac{d^2}{b^2}}}}{b^2} &= ? \\ \frac{c^2}{a^2(\frac{c^2}{a^2}+\frac{d^2}{b^2})} + \frac{d^2}{b^2(\frac{d^2}{a^2}+\frac{d^2}{b^2})} &= ? \\ \frac{c^2}{c^2 + \frac{a^2d^2}{b^2}} + \frac{d^2}{\frac{b^2a^2}{a^2}d^2} &= ? \\ \frac{b^2c^2}{b^2c^2+a^2d^2} + \frac{a^2d^2}{b^2c^2+a^2d^2} &= 1 \end{aligned} a2a2c2+b2d2 c+b2a2c2+b2d2 da2(a2c2+b2d2)c2+b2(a2d2+b2d2)d2c2+b2a2d2c2+a2b2a2d2d2b2c2+a2d2b2c2+b2c2+a2d2a2d2=?=?=?=1

说明step2得到的方向向量是满足椭圆方程的。

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值