椭圆及其标准方程
椭圆是一种特殊的几何图形，可以定义为平面上到两个固定点（焦点）的距离之和为常数的点的轨迹。在笛卡尔坐标系中，当椭圆的中心位于 (h, k) 且其主轴与坐标轴平行时，其标准方程为：
(x - h)² / a² + (y - k)² / b² = 1
其中：
- (x, y) 是椭圆上任意一点的坐标。
- (h, k) 是椭圆的中心坐标。
- a 是水平半轴的长度（如果主轴平行于X轴）。
- b 是垂直半轴的长度（如果主轴平行于Y轴）。
需要注意的是，a 和 b 分别代表半长轴和半短轴的长度，它们分别对应着在X轴和Y轴方向上的“半径”。如果 a > b，则水平方向是长轴；如果 b > a，则垂直方向是长轴。
判断点是否在椭圆内部的原理
要判断一个点 (pt_x, pt_y) 是否在椭圆的内部或边界上，我们只需将该点的坐标代入椭圆的标准方程的左侧，然后将计算结果与 1 进行比较。
- 如果 (pt_x - h)² / a² + (pt_y - k)² / b² < 1，则点在椭圆的内部。
- 如果 (pt_x - h)² / a² + (pt_y - k)² / b² = 1，则点在椭圆的边界上。
- 如果 (pt_x - h)² / a² + (pt_y - k)² / b² > 1，则点在椭圆的外部。
综合来看，如果点在椭圆内部或边界上，则满足条件： (pt_x - h)² / a² + (pt_y - k)² / b² <= 1
示例：判断点是否在特定椭圆内
假设我们有一个椭圆，其属性如下：
- 中心 center = [10, 10]
- 半径 radius = [3, 4]，其中 radius[0] 是垂直半径（y轴方向），radius[1] 是水平半径（x轴方向）。
根据这些属性，我们可以确定：
- h = 10
- k = 10
- 水平半轴 a = 4 (对应 x 坐标)
- 垂直半轴 b = 3 (对应 y 坐标)
因此，该椭圆的方程为： (x - 10)² / 4² + (y - 10)² / 3² = 1 即 (x - 10)² / 16 + (y - 10)² / 9 = 1
要判断一个点 pt = [pt_x, pt_y] 是否在此椭圆内，我们需要计算： (pt_x - 10)² / 16 + (pt_y - 10)² / 9 并将结果与 1 进行比较。
JavaScript 实现示例
以下是一个JavaScript函数，用于实现上述判断逻辑：
/**
* 判断一个点是否在给定椭圆的内部或边界上。
*
* @param {Array<number>} pt 点的坐标，例如 [pt_x, pt_y]。
* @param {Array<number>} center 椭圆中心的坐标，例如 [h, k]。
* @param {Array<number>} radii 椭圆的半轴长度，例如 [verticalRadius, horizontalRadius]。
* 注意：radii[0] 对应 y 轴方向的半轴，radii[1] 对应 x 轴方向的半轴。
* @returns {boolean} 如果点在椭圆内部或边界上，则返回 true；否则返回 false。
*/
function isPointInEllipse(pt, center, radii) {
const pt_x = pt[0];
const pt_y = pt[1];
const h = center[0]; // 椭圆中心x坐标
const k = center[1]; // 椭圆中心y坐标
// 根据问题描述，radii[0] 是垂直半径 (y轴方向)，radii[1] 是水平半径 (x轴方向)
const horizontalSemiAxis = radii[1]; // 对应方程中的 'a'
const verticalSemiAxis = radii[0]; // 对应方程中的 'b'
// 避免除以零的情况
if (horizontalSemiAxis === 0 || verticalSemiAxis === 0) {
// 如果任何一个半轴为零，则椭圆退化为一条线段或一个点。
// 在这种情况下，可以根据具体业务需求定义行为。
// 这里简单地认为点不在“内部”。
return false;
}
// 计算椭圆方程左侧的值
const value = Math.pow((pt_x - h), 2) / Math.pow(horizontalSemiAxis, 2) +
Math.pow((pt_y - k), 2) / Math.pow(verticalSemiAxis, 2);
// 判断值是否小于等于 1
return value <= 1;
}
// 示例使用：
const center = [10, 10];
const radii = [3, 4]; // verticalRadius = 3, horizontalRadius = 4
// 测试点
const pointInside = [12, 11]; // 应该在内部
const pointOnBoundary = [14, 10]; // 应该在边界上 (x=14, y=10) => (4^2/4^2 + 0/3^2) = 1
const pointOutside = [15, 10]; // 应该在外部
console.log(`点 ${pointInside} 是否在椭圆内: ${isPointInEllipse(pointInside, center, radii)}`); // 预期: true
console.log(`点 ${pointOnBoundary} 是否在椭圆内: ${isPointInEllipse(pointOnBoundary, center, radii)}`); // 预期: true (在边界上)
console.log(`点 ${pointOutside} 是否在椭圆内: ${isPointInEllipse(pointOutside, center, radii)}`); // 预期: false
// 更多测试点
console.log(`点 [10, 10] (中心) 是否在椭圆内: ${isPointInEllipse([10, 10], center, radii)}`); // 预期: true
console.log(`点 [10, 13] (上边界) 是否在椭圆内: ${isPointInEllipse([10, 13], center, radii)}`); // 预期: true
console.log(`点 [10, 14] (上外部) 是否在椭圆内: ${isPointInEllipse([10, 14], center, radii)}`); // 预期: false注意事项与总结
- 半轴对应关系：务必明确 a 和 b 哪个对应水平半轴，哪个对应垂直半轴。在提供的例子中，radius = [3, 4] 明确指出 3 是垂直半径，4 是水平半径，因此在方程中 x 坐标差的平方除以 4²，y 坐标差的平方除以 3²。
- 边界情况：当计算结果恰好等于 1 时，表示点位于椭圆的边界上。根据需求，可能需要区分“严格内部”和“内部或边界”。本教程的实现 value <= 1 包含了边界点。
- 退化椭圆：如果任何一个半轴长度为零，椭圆将退化为一条线段或一个点。在实际应用中，需要考虑如何处理这些特殊情况，例如在代码中添加对 radii 值为零的检查。
- 旋转椭圆：本教程讨论的是主轴与坐标轴平行的椭圆。对于经过旋转的椭圆，其方程会更加复杂，通常需要使用旋转矩阵或更通用的二次曲线方程来处理。
通过掌握椭圆的标准方程和点包含的原理，我们可以高效且准确地判断任意点与非旋转椭圆的相对位置，这在图形学、游戏开发、地理信息系统等领域都有广泛的应用。
