//计算定向 int num = 0; int NumDirected = 0;//定向点号 CAngle Z0;//定向坐标 NumDirected = find(pAngleObs, pAngleObs[i].pStation->strID, 0); T0 = Azi(pAngleObs[i], pAngleObs[i]);//近似方位角，现在是DMS /*double aa = pAngleObs[i].pStation->x; double bb = pAngleObs[i].pStation->y;*///检核所用 Z0 = Azi(pAngleObs[NumDirected], pAngleObs[NumDirected]);//为DMS //double cc = pAngleObs[NumDirected].pObs->x; //double dd = pAngleObs[NumDirected].pObs->y;检核所用 CAngle L0; L0(RAD) = -Z0(RAD) + T0(RAD); if (L0(DEG) < 0) { L0(DEG) += 360; } CAngle lij; lij(DEG) = pAngleObs[i].ObsAngle(DEG) - L0(DEG); L(i, 0) = lij(DEG) * 3600;//转为秒

这段代码用于计算定向角度。首先，定义了一些变量，包括num、NumDirected、Z0、T0、L0和lij。

通过调用find函数，根据观测站点的名称在角度观测值数组pAngleObs中查找对应的定向观测值的索引NumDirected。

然后，通过调用Azi函数，分别计算当前观测值和定向观测值的近似方位角T0和Z0，其中近似方位角的单位为DMS。

接下来，通过计算L0，即将-Z0（单位为RAD）与T0（单位为RAD）相加并取负，得到L0的值。如果L0的单位为DEG（度）小于0，则将其加上360，保证L0在0到360之间。

最后，通过计算lij，即当前观测值的观测角度（单位为DEG）减去L0的值，得到lij的值。将lij转换为秒，并存入矩阵L的第i行第0列。

需要注意的是，代码中的pAngleObs是角度观测值数组，pAngleObs[i].pStation表示当前观测值对应的测站指针，pAngleObs[i].ObsAngle表示当前观测值的观测角度，Azi函数用于计算近似方位角。另外，L是矩阵对象。