发明名称 一种旋转摄影测量系统多片整体平差标定方法
摘要 一种旋转摄影测量系统多片整体平差标定方法,利用单个测站获取多张影像的外方位元素和旋转平台的水平、垂直角,进行旋转摄影测量系统标定。从而在其他测站已知第一张影像外方位元素的前提下,可根据定标结果和旋转平台的水平、垂直旋角,高精度地反推其他影像的外方位元素。本发明具有如下优点:利用单个测站获取多张影像的外方位元素和旋转平台的水平、垂直旋角,通过多片整体平差解算,该标定方法能够实现旋转摄影测量系统高精度的标定;在其他测站已知第一张影像外方位元素的前提下,仅需提供其他影像成像时旋转平台的水平、垂直旋角,便可高精度的自动解算该影像的外方位元素。
申请公布号 CN102519484A 申请公布日期 2012.06.27
申请号 CN201110385815.2 申请日期 2011.11.29
申请人 武汉大学 发明人 张永军;胡堃;柯涛;黄山
分类号 G01C25/00(2006.01)I 主分类号 G01C25/00(2006.01)I
代理机构 武汉科皓知识产权代理事务所(特殊普通合伙) 42222 代理人 薛玲
主权项 1.一种旋转摄影测量系统多片整体平差标定方法,其特征在于,包括以下步骤:步骤1,导入第一个测站获取的影像数目<img file="171246DEST_PATH_IMAGE001.GIF" wi="20" he="20" />、各张影像获取时的外方位角元素<img file="2011103858152100001DEST_PATH_IMAGE002.GIF" wi="74" he="26" />和外方位线元素<img file="150703DEST_PATH_IMAGE003.GIF" wi="90" he="26" />、旋转平台绕垂直旋转轴的水平角<img file="2011103858152100001DEST_PATH_IMAGE004.GIF" wi="28" he="26" />和旋转平台绕水平旋转轴的垂直角<img file="930440DEST_PATH_IMAGE005.GIF" wi="22" he="26" />,其中下标1标识测站号<img file="2011103858152100001DEST_PATH_IMAGE006.GIF" wi="9" he="18" />=1,下标<img file="385037DEST_PATH_IMAGE007.GIF" wi="14" he="21" />表示影像号,<img file="2011103858152100001DEST_PATH_IMAGE008.GIF" wi="78" he="22" />,影像数目<img file="261726DEST_PATH_IMAGE001.GIF" wi="20" he="20" />大于等于3;由外方位角元素<img file="84188DEST_PATH_IMAGE002.GIF" wi="74" he="26" />根据<img file="351222DEST_PATH_IMAGE009.GIF" wi="65" he="18" />转角系统构建旋转矩阵<img file="2011103858152100001DEST_PATH_IMAGE010.GIF" wi="25" he="26" />,由外方位线元素<img file="612439DEST_PATH_IMAGE003.GIF" wi="90" he="26" />构建线元素的坐标矩阵<img file="281317DEST_PATH_IMAGE011.GIF" wi="144" he="28" />;步骤2,构建相机相对于物方坐标系的旋转模型,旋转模型公式如下<img file="2011103858152100001DEST_PATH_IMAGE012.GIF" wi="94" he="25" />其中,<img file="336998DEST_PATH_IMAGE013.GIF" wi="17" he="18" />为物方坐标系下相机外方位角元素在<img file="91327DEST_PATH_IMAGE009.GIF" wi="65" he="18" />转角系统下旋转矩阵;<img file="2011103858152100001DEST_PATH_IMAGE014.GIF" wi="27" he="25" />为物方坐标系与标准摄影测量坐标系间的旋转矩阵;<img file="890656DEST_PATH_IMAGE015.GIF" wi="30" he="25" />为标准摄影坐标系与旋转坐标系间的旋转矩阵;<img file="2011103858152100001DEST_PATH_IMAGE016.GIF" wi="21" he="25" />为旋转坐标系与像空间坐标系间的旋转矩阵;根据第一个测站获取的前三张影像分别对应的旋转模型<img file="414041DEST_PATH_IMAGE017.GIF" wi="115" he="27" />,解算第一个测站对应的<img file="906203DEST_PATH_IMAGE009.GIF" wi="65" he="18" />转角系统下旋转矩阵<img file="2011103858152100001DEST_PATH_IMAGE018.GIF" wi="31" he="25" />的转角<img file="147828DEST_PATH_IMAGE019.GIF" wi="94" he="25" />和旋转矩阵<img file="485269DEST_PATH_IMAGE016.GIF" wi="21" he="25" />的转角<img file="DEST_PATH_IMAGE020.GIF" wi="62" he="25" />,并以所得解算结果作为未知数的初值,其中,下标1标识测站号<img file="128740DEST_PATH_IMAGE006.GIF" wi="9" he="18" />=1,下标<img file="526223DEST_PATH_IMAGE007.GIF" wi="14" he="21" />表示影像号,<img file="255144DEST_PATH_IMAGE021.GIF" wi="62" he="22" />,<img file="333959DEST_PATH_IMAGE010.GIF" wi="25" he="26" />为步骤1构建的旋转矩阵;步骤3,将旋转模型转化为误差方程,以转角<img file="897183DEST_PATH_IMAGE019.GIF" wi="94" he="25" />和转角<img file="668830DEST_PATH_IMAGE020.GIF" wi="62" he="25" />为未知数,对第一个测站获取的所有影像逐张进行线性化处理,根据未知数的初值按照最小二乘原理构建法方程,整体平差求解法方程,得到未知数改正数;步骤4,若未知数改正数最大值小于预设的未知数改正数阈值或者迭代次数超过预设的迭代次数阈值,则执行步骤5;否则,以当前未知数改正数为未知数的初值,返回迭代执行步骤3;步骤5,输出最后一次迭代中获取的未知数改正数,作为转角<img file="619468DEST_PATH_IMAGE019.GIF" wi="94" he="25" />和转角<img file="564291DEST_PATH_IMAGE020.GIF" wi="62" he="25" />的标定结果,根据转角<img file="916774DEST_PATH_IMAGE019.GIF" wi="94" he="25" />的标定结果重新组建旋转矩阵<img file="859323DEST_PATH_IMAGE018.GIF" wi="31" he="25" />,跟据转角<img file="359574DEST_PATH_IMAGE020.GIF" wi="62" he="25" />的标定结果重新组建旋转矩阵<img file="780191DEST_PATH_IMAGE016.GIF" wi="21" he="25" />;步骤6,构建相机相对于物方坐标系的偏心模型,偏心模型公式如下<img file="DEST_PATH_IMAGE022.GIF" wi="121" he="25" />其中<img file="315078DEST_PATH_IMAGE023.GIF" wi="109" he="26" />,为物方坐标系下相机外方位线元素的坐标矩阵;<img file="DEST_PATH_IMAGE024.GIF" wi="18" he="25" />为像空间坐标系原点在旋转坐标系下的坐标;<img file="428527DEST_PATH_IMAGE025.GIF" wi="121" he="26" />,为旋转中心在物方坐标系下的坐标矩阵;根据相机相对于物方坐标系的偏心模型<img file="DEST_PATH_IMAGE026.GIF" wi="145" he="26" />,根据第一个测站获取的影像逐张构建误差方程和法方程,按照最小二乘原理,整体平差解算像空间坐标系原点在旋转坐标系下的坐标<img file="416075DEST_PATH_IMAGE024.GIF" wi="18" he="25" />和旋转中心在物方坐标系下的坐标<img file="640383DEST_PATH_IMAGE027.GIF" wi="24" he="25" />,输出坐标<img file="764197DEST_PATH_IMAGE024.GIF" wi="18" he="25" />的解算结果,其中,下标1标识测站号<img file="48547DEST_PATH_IMAGE006.GIF" wi="9" he="18" />=1,下标<img file="461074DEST_PATH_IMAGE007.GIF" wi="14" he="21" />表示影像号,<img file="551390DEST_PATH_IMAGE008.GIF" wi="78" he="22" />,<img file="DEST_PATH_IMAGE028.GIF" wi="22" he="26" />为步骤1构建的坐标矩阵,<img file="732973DEST_PATH_IMAGE018.GIF" wi="31" he="25" />为步骤5重新组建的旋转矩阵;步骤7,根据步骤5重新组建的旋转矩阵<img file="984962DEST_PATH_IMAGE016.GIF" wi="21" he="25" />和步骤6所得像空间坐标系原点在旋转坐标系下的坐标<img file="884785DEST_PATH_IMAGE024.GIF" wi="18" he="25" />的解算结果、已知的其余测站第一张影像获取时的外方位角元素<img file="450896DEST_PATH_IMAGE029.GIF" wi="68" he="25" />和外方位线元素<img file="DEST_PATH_IMAGE030.GIF" wi="82" he="25" />、旋转平台绕垂直旋转轴的水平角<img file="835389DEST_PATH_IMAGE031.GIF" wi="26" he="25" />和绕水平旋转轴的垂直角<img file="DEST_PATH_IMAGE032.GIF" wi="20" he="25" />,通过旋转模型和偏心模型重新解算在其余测站条件下,物方坐标系与标准摄影测量坐标系间的旋转矩阵<img file="258280DEST_PATH_IMAGE033.GIF" wi="30" he="25" />和旋转中心在物方坐标系下的坐标<img file="DEST_PATH_IMAGE034.GIF" wi="22" he="25" />;其中下标<img file="645399DEST_PATH_IMAGE006.GIF" wi="9" he="18" />表示测站号,<img file="77517DEST_PATH_IMAGE035.GIF" wi="64" he="22" />,下标1表示影像号<img file="968113DEST_PATH_IMAGE007.GIF" wi="14" he="21" />=1;步骤8,根据已知的其余测站第一张以外其他影像获取时旋转平台绕垂直旋转轴的水平角<img file="DEST_PATH_IMAGE036.GIF" wi="25" he="26" />和绕水平旋转轴的垂直角<img file="827484DEST_PATH_IMAGE037.GIF" wi="20" he="26" />,通过旋转模型和偏心模型解算其余测站第一张以外其他影像获取时的外方位角元素<img file="DEST_PATH_IMAGE038.GIF" wi="65" he="26" />和外方位线元素<img file="498637DEST_PATH_IMAGE039.GIF" wi="80" he="26" />,其中下标<img file="672130DEST_PATH_IMAGE006.GIF" wi="9" he="18" />表示测站号,<img file="745128DEST_PATH_IMAGE035.GIF" wi="64" he="22" />,下标<img file="713084DEST_PATH_IMAGE007.GIF" wi="14" he="21" />表示影像号,<img file="DEST_PATH_IMAGE040.GIF" wi="66" he="22" />。
地址 430072 湖北省武汉市武昌区珞珈山武汉大学