A(x1,y1,z1), B(x2,y2,z2), и dA(dx, dy, dz)
!!A и B вычисляются переводом сферических координат в декартову!!
старое расстояние R1=Корень(квадрат(x1-x2)+квадрат(y1-y2)+квадрат(z1-z2))
новое расстояние R2=Корень(квадрат(x1-x2+dx)+квадрат(y1-y2+dy)+квадрат(z1-z2+dz))
разница расстояний dR = R2-R1;
прочитайте про перевод из сферических координат в декартову тут, и будет Вам счастье
По другому никак
UPD:
Есть вектор A(ϴa,φa), Вектор B(ϴb,φb), и вектор A2(ϴa+dϴ,φa+dφ)
они же в декартовой системе A(Xa,Ya,Za) B(Xb,Yb,Zb),
Где Xi=sin(ϴi)*cos(φi), Yi=sin(ϴi)*sin(φi), Zi=cos(ϴi)
Найдем угол между ними αab=ArcCos(ab'/(a'*b'))
где ab'=XaXb+YaYb+Za*Zb, a'=Корень(Xa^2+Ya^2+Za^2), b'=Корень(Xb^2+Yb^2+Zb^2)
и длинну дуги Lab=αab*R
т.е. примерно так:
public double GetSphereSectorLen(double ϴ1,double φ1,double ϴ2,double φ2, double R)
{
double X1=Math.Sin(ϴ1)*Math.Cos(φ1);
double Y1=Math.Sin(ϴ1)*Math.Sin(φ1);
double Z1=Math.Cos(ϴ1);
double X2=Math.Sin(ϴ2)*Math.Cos(φ2);
double Y2=Math.Sin(ϴ2)*Math.Sin(φ2);
double Z2=Math.Cos(ϴ2);
double ab_=X1*X2+Y1*Y2+Z1*Z2;
double a_=Math.Sqrt(Math.Pow(X1,2)+Math.Pow(Y1,2)+Math.Pow(Z1,2));
double b_=Math.Sqrt(Math.Pow(X2,2)+Math.Pow(Y2,2)+Math.Pow(Z2,2));
double αab=Math.Acos(ab_/(a_*b_));
return αab*R;
}
double LAB=GetSphereSectorLen(ϴa,φa,ϴb,φb);
double LA2B=GetSphereSectorLen(ϴa+dϴ,φa+dφ,ϴb,φb);
double dL = LAB-LA2B;