给定隐式圆锥曲线方程如何获取曲线上的点用于生成网格？

核心思路是先把带交叉项的隐式方程转换为无交叉项的标准型，采样后再逆变换回原坐标系，全程仅需数值计算，不依赖绘图库，输出的点集可直接用于后续仿射变换与网格生成。 #

前置步骤：判别圆锥曲线类型 #

首先用二次项判别式确认曲线类型，方便后续选择对应采样逻辑：

• 计算判别式 Δ = B² - 4AC
  • Δ < 0：椭圆（含圆形）
  • Δ = 0：抛物线
  • Δ > 0：双曲线

步骤1：消去交叉项，转换为标准型 #

原方程中的Bxy交叉项会增加采样难度，先通过旋转变换消去：

1. 计算主轴旋转角度：θ = 0.5 * arctan2(B, A - C)，若A=C可直接取θ = π/4
2. 旋转变换公式：将原坐标(x,y)映射为旋转后坐标系的(x', y')

   x = x' * cosθ - y' * sinθ
   y = x' * sinθ + y' * cosθ
   

3. 将上式代入原隐式方程，消去x'y'交叉项，得到新方程 A'x'² + C'y'² + D'x' + E'y' + F' = 0
4. 再做平移变换，消去一次项，将曲线中心移到旋转后坐标系的原点，得到标准型方程（如椭圆标准型x'²/a² + y'²/b² = 1）

步骤2：标准坐标系下采样 #

根据曲线类型选择对应采样逻辑，采样密度可根据你的精度要求自行调整：

• 椭圆/圆：用参数方程均匀采样，参数t取值范围[0, 2π]，步长越小精度越高：

  x' = a * cos(t)
  y' = b * sin(t)
  

• 抛物线：以抛物线开口方向的对应坐标为自变量，均匀取值后反算另一坐标即可
• 双曲线：根据开口方向选择对应分支的参数范围，采样后仅保留你需要的分支点集即可

步骤3：逆变换得到原平面坐标 #

将标准坐标系下采样得到的所有(x', y')点，先做平移逆变换、再做旋转逆变换，即可得到原xy平面下的圆锥曲线点集。

精度验证方法：将得到的(x,y)代入原隐式方程，计算Ax²+Bxy+Cy²+Dx+Ey+F的绝对值，只要小于预设阈值（如1e-6）即符合精度要求。

备选通用方案：无需类型判别 #

如果不想做坐标变换与类型判定，可直接用Marching Squares算法提取隐式等值线：对xy平面的矩形网格做遍历，提取Ax²+Bxy+Cy²+Dx+Ey+F = 0的等值点，该方法适配所有类型的隐式曲线，精度由网格步长决定。

内容的提问来源于stack exchange，提问作者Inyoung Kim 김인영