一种铁路线路及沿线区域形变信息的提取方法

1.一种铁路线路及沿线区域形变信息的提取方法，其特征在于，包括：获取铁路线路中线位置各里程对应的经纬度坐标P1；获取观测周期内铁路线路及沿线观测区域中的多个配置角反射器的形变基准点P2位置的GNSS连续观测值G1；应用天基微波遥感影像进行InSAR时序分析解算，获得所述观测周期内，观测区域内的地表目标体微波散射点的形变信息S1，S1中包括形变基准点P2位置的角反射点的微波遥感形变数据S11；根据所述观测区域内的地表目标体微波散射点的形变信息S1，以形变基准点P2位置的GNSS连续观测值G1测网平差结果G11作为优化条件，以同位置微波遥感形变数据S11与G11相对应时刻的形变差做为直接约束条件，拟合获得观测区域的形变曲面S2；根据拟合观测区域的形变曲面S2及铁路线路中线位置各里程对应的经纬度坐标P1，构建线路区域影响范围；提取所述观测区域面形变中，且在所述线路区域影响范围内的铁路线路形变S3；其中，拟合获得观测区域的形变曲面S2，包括：应用非线性最小二乘法拟合获得观测区域的形变曲面S2；获取微波散射点的形变点集Dsci＝x i ，y i ，d i ；其中i＝1,2，…m为观测区域微波散射点；x i ，y i ，d i 为微波散射点的经纬度坐标和形变探测值；将形变基准点P2位置的微波遥感形变数据S11与GNSS连续观测数据点中相应时间的形变数据差值集Dgi＝x j ，y j ，d j ，作为约束条件，其中j＝1,2,…,n为观测区域内形变基准点，x j ,y j ,d j 为形变基准点的经纬度坐标和形变观测值；以G1测网平差结果G11作为优化条件；拟合目标为观测区域的面形变Dz＝f(x，y)，目标函数满足d j ＝f(x j ，y j )，j＝1,2,…,n；根据形变基准点的分布情况对观测区域进行分段，寻找每段内的局部最优解并保持各段间的函数一次导数连续；在区段内满足公式(1)其中：i＝1,2…k,k<m的条件，且受约束条件d j ＝f(x j ，y j )，j＝1,2,…,n限制；进行局部寻优，获得观测区域内其他位置的形变信息f＝(x，y)，其中x,y不包括微波散射点和形变基准点位置；依次实施每段的拟合，获得观测区域的形变曲面S2。 2.根据权利要求1所述的一种铁路线路及沿线区域形变信息的提取方法，其特征在于，所述方法还包括：提取所述铁路线路区域影响范围外的铁路沿线区域地表结构形变S4。 3.根据权利要求1所述的一种铁路线路及沿线区域形变信息的提取方法，其特征在于，获取铁路线位中线位置各里程对应的经纬度坐标P1，包括：根据选取的形变基准点，将所述形变基准点的线路里程信息与线路位置建立映射关系；所述线路位置为空间经纬度坐标信息；将线路中线位置离散为点集，根据所述映射关系，获得线路里程离散点集的空间经纬度坐标。