IDL编程学习之三点定圆、四点定球

15195775117

3点定圆，4点定球，IDL有专门的函数
IDL> X = [1.0, 2.0, 3.0]
IDL> Y = [1.0, 2.0, 1.0]
IDL> CIR_3PNT, X, Y, R, X0, Y0
结果： r=1.0；x0=2.0；y0=1.0
4点定球的函数：
IDL> X = [1, 2, 1, 1]
IDL> Y = [1, 1, 0, 0]
IDL> Z = [0, 2, 3, 1]
IDL> SPH_4PNT, X, Y, Z, Xc, Yc, Zc, R
IDL> print,Xc, Yc, Zc, R
    -0.500000      2.00000      2.00000      2.69258