ITRF瞬时历元空间直角坐标转为CGCS2000坐标再思考

     以前也发过类似的文章，关于ITRF2020当前历元坐标转为CGCS2000坐标：https://blog.csdn.net/yisdong/article/details/135668976?spm=1001.2014.3001.5502，那篇文章初步实现了坐标的转换，不过计算公式不是很严密，水平有限，才做到2～5cm准确度，如有高见，请指正。最近根据生产实践和同行的讨论中发现一些新的计算方法和软件，来给大家分享下。

1. 已知点

       1.1   利用北京房山IGS观测站数据计算的空间直角坐标（ITRF2014框架，历元2021.0），bjfs   -2148744.5882   4426641.1511   4044655.7902

       1.2   IGS站已知点 空间直角坐标（ITR97框架，历元2000.0）

       2.计算方法

这里涉及3个历元，我们现在测一个点框架是ITRF2014，历元是测量时间约为2021（历元1），ITRF2014->ITRF1997的参数对应历元是2010(历元2)，CGCS2000对应的历元是2000（历元3）。  

   规范说的不是太清楚，转换过程目前有两种观点：
  观点一：①：ITRF2014->ITRF1997参数由2010(历元2)归化到2021（历元1），进行框架转换。
    ②：2021（历元1）到2000（历元3）的板块运动改正。
   观点二：    ①：ITRF2014->ITRF1997参数由2021(历元1)归化到2010（历元2），进行框架和历元转换。    ②：20