空间后方交会原理
空间后方交会的解算
空间后方交会的解算一. 空间后方交会的目的摄影测量主要利用摄影的方法获取地面的信息,主要是是点位信息,属性信息,因此要对此进行空间定位和建模,并首先确定模型的参数,这就是空间后方交会的目的,用以求出模型外方位元素。
二. 空间后方交会的原理空间后方交会的原理是共线方程。
共线方程是依据相似三角形原理给出的,其形式如下111333222333()()()()()()()()()()()()A S A S A S A S A S A S AS A S A S A S A S A S a X X b Y Y c Z Z x f a X X a Y Y a Z Z a X X b Y Y c Z Z y f a X X a Y Y a Z Z -+-+-=--+-+--+-+-=--+-+-上式成为中心投影的构线方程,我们可以根据几个已知点,来计算方程的参数,一般需要六个方程,或者要三个点,为提高精度,可存在多余观测,然后利用最小二乘求其最小二乘解。
将公式利用泰勒公式线性化,取至一次项,得到其系数矩阵A ;引入改正数(残差)V ,则可将其写成矩阵形式:V AX L =-其中111333222333[,]()()()()()()()()()()()()()()Tx y A S A S A S x A S A S A S A S A S A S y A S A S A S L l l a X X b Y Y c Z Z l x x x fa X X a Y Y a Z Z a X Xb Y Yc Z Z l y y y fa X X a Y Y a Z Z =-+-+-=-=+-+-+--+-+-=-=+-+-+- 则1()T T X A A A L -=X 为外方位元素的近似改正数,由于采用泰勒展开取至一次项,为减少误差,要将的出的值作为近似值进行迭代,知道小于规定的误差三. 空间后方交会解算过程1. 已知条件近似垂直摄影00253.24mmx y 0f ===2. 解算程序流程图MATLAB 程序format long;s1=xlsread('data.xls');%读取数据a1=s1(1:4,1:2);%影像坐标b1=s1(1:4,3:5);%地面摄影测量坐标a2=s1.*10^-3;%影像坐标单位转化j1=a2(1,:)-a2(2,:);j2=j1(1,1)^2+j1(1,2)^2;lengh_a1=sqrt(j2); %相片某一长度j1=b1(1,:)-b1(1,:);j2=j1(1,1)^2+j1(1,2)^2;lengh_b1=sqrt(j2); %地面对应的长度m=lengh_b1/lengh_a1;%求出比例尺n0=0;p0=0;q0=0;x0=mean(b1(:,1));y0=mean(b1(:,2));f=153.24*10^-3;z0=m*f;x001={x0,x0,x0,x0};X0=cell2mat(x001)';y001={y0,y0,y0,y0};Y0=cell2mat(y001)';z001={z0,z0,z0,z0};Z0=cell2mat(z001)';%初始化外方位元素的值aa1=cos(n0)*cos(q0)-sin(n0)*sin(p0)*sin(q0);aa2=-sin(q0)*cos(n0)-sin(n0)*sin(p0)*cos(q0);aa3=-sin(n0)*cos(p0);bb1=sin(q0)*cos(p0);bb2=cos(q0)*cos(p0);bb3=-sin(p0);cc1=sin(n0)*cos(q0)+sin(p0)*cos(n0)*sin(q0);cc2=-sin(n0)*sin(q0)+sin(p0)*cos(q0)*cos(n0);cc3=cos(n0)*cos(p0);%计算改正数XX1=aa1.*(b1(:,1)-X0)+bb1.*(b1(:,2)-Y0)+cc1.*(b1(:,3)-Z0); XX2=aa2.*(b1(:,1)-X0)+bb2.*(b1(:,2)-Y0)+cc2.*(b1(:,3)-Z0); XX3=aa3.*(b1(:,1)-X0)+bb3.*(b1(:,2)-Y0)+cc3.*(b1(:,3)-Z0); lx=a1(:,1)+f.*(XX1./XX3);ly=a1(:,2)+f.*(XX2./XX3);l={lx',ly'};L=cell2mat(l)';%方程系数A=[-3.969*10^-5 0 2.231*10^-5 -0.2 -0.04 -0.06899;0 -3.969*10^-5 1.787*10^-5 -0.04 -0.18 0.08615;-2.88*10^-5 0 1*10^-5 -0.17 0.03 0.08211;0 -2.88*10^-5 -1.54*10^-5 0.03 -0.2 0.0534;-4.14*10^-5 0 4*10^-6 -0.15 -7.4*10^-3 -0.07663;0 -4.14*10^-5 2.07*10^-5 -7.4*10^-3 -0.19 0.01478;-2.89*10^-5 0 -1.98*10^-6 -0.15 -4.4*10^-3 0.06443;0 -2.89*10^-5 -1.22*10^-5 -4.4*10^-3 -0.18 0.01046];%L=[-1.28 3.78 -3.02 -1.45 -4.25 4.98 -4.72 -0.385]'.*10^-2; %第一次迭代X=inv(A'*A)*A'*L;3.结果X=1492.41127406195-554.4015671761941425.68660973544-0.0383847815608609 0.00911624039769785 -0.105416434087641S=1492.41127406195-554.401567176194 1425.68660973544 38436.9616152184 27963.1641162404-0.105416434087641。
单像空间后方交会实习报告
单像空间后方交会实习报告一、实习目的单像空间后方交会是摄影测量中确定像片外方位元素的重要方法。
通过本次实习,旨在深入理解单像空间后方交会的基本原理和计算过程,熟练掌握相关软件的操作,提高对摄影测量数据处理的实践能力,并培养解决实际问题的思维和方法。
二、实习原理单像空间后方交会的目的是利用像片上的像点坐标以及相应的地面控制点坐标,通过数学模型求解像片的外方位元素(三个线元素 Xs、Ys、Zs 和三个角元素φ、ω、κ)。
其基本原理基于共线条件方程,即摄影中心、像点和相应的地面点位于同一条直线上。
共线条件方程可以表示为:\\begin{align}x x_0&= f\frac{a_1(X X_s) + b_1(Y Y_s) + c_1(Z Z_s)}{a_3(X X_s) + b_3(Y Y_s) + c_3(Z Z_s)}\\y y_0&= f\frac{a_2(X X_s) + b_2(Y Y_s) + c_2(Z Z_s)}{a_3(X X_s) + b_3(Y Y_s) + c_3(Z Z_s)}\end{align}\其中,\((x,y)\)为像点坐标,\((x_0,y_0)\)为主点坐标,\(f\)为摄影机焦距,\((X,Y,Z)\)为地面点的物方空间坐标,\((X_s,Y_s,Z_s)\)为摄影中心的物方空间坐标,\((a_1,b_1,c_1),(a_2,b_2,c_2),(a_3,b_3,c_3)\)为由角元素φ、ω、κ 构成的旋转矩阵的元素。
三、实习数据本次实习使用了一组航空像片,像片比例尺为 1:5000,焦距为152mm,像主点坐标为\((x_0,y_0)=(5000mm,5000mm)\)。
同时,提供了 6 个均匀分布在像片范围内的地面控制点的物方空间坐标和像点坐标。
四、实习步骤1、数据准备整理地面控制点的物方空间坐标和像点坐标,确保数据的准确性。
输入像片的基本参数,如像主点坐标、焦距等。
后方交会残差值误差范围
后方交会残差值误差范围后方交会是摄影测量中常用的一种方法,用于确定地面上各个点的空间坐标。
在实际应用中,由于各种误差的存在,后方交会的结果会产生一定的残差值误差。
误差范围的确定对于保证测量结果的准确性和可靠性非常重要。
本文将从后方交会的基本原理、误差来源、误差计算方法以及误差范围的确定等方面进行详细的分析和论述。
一、后方交会的基本原理后方交会是一种基于像对几何关系的摄影测量方法,通过对各个像点的位置测量和相对方位角的观测,计算出地面控制点的空间坐标。
其基本原理可以简述如下:1. 反投影原理:根据像点在像空间上的位置,利用摄影测量的几何关系反推出这些像点所对应的地面点在物空间上的位置。
反投影原理是后方交会的理论基础,也是误差产生的根源。
2. 控制点观测:确定一定数量的控制点,并测量其像点位置及相对方位角。
控制点的选择应满足精度要求和实际情况,通常采用地面测量或其他摄影测量方法进行。
3. 几何模型:根据反投影原理和控制点观测,建立几何模型,描述像空间与物空间之间的几何关系。
模型包括相机的内外参数、像点的位置和相对方位角等。
4. 误差方程:利用几何模型,建立误差方程,将测量值与真实值之间的误差表示出来。
误差方程是分析误差来源、计算误差范围的基础。
二、后方交会误差的来源后方交会的误差主要来自于以下几个方面:1. 相机内外参数的误差:相机的内外参数是后方交会的重要参数,包括焦距、主点位置、旋转矩阵、平移向量等。
由于摄影测量设备和仪器的制造和使用限制,这些参数会存在误差,从而影响后方交会的结果。
2. 观测误差:观测误差包括控制点的像点测量误差和方位角观测误差。
像点测量误差可以由像点测量精度来描述,方位角观测误差可以由方位角观测精度来描述。
观测误差是由测量设备、操作人员和环境等因素共同引起的。
3. 地面控制点的精度:后方交会的精度还受到地面控制点的精度限制。
如果地面控制点的精度较差,那么后方交会的精度也会受到影响。
单张像片的空间后方交会原理
单像空间后方交会原理非常简单,稍稍懂一点线性代数和微分学知识就很容易理解。
但它的意义很大,所有的基于视觉技术的定位(测量)都是从它扩展来的,没有它就没有后面的复杂算法的演进,以面向更复杂的场景以及需求。
这块内容一般来说就是在已知地面上若干点的地面坐标以后,反求该相应摄影光束的外参数,当用以作摄影机的标定时,还可借以同时求出摄影的内参数。
具体的说就是一个多元线性模型,根据不同的需求可以基于它进一步构建各种平差模型。
摄影测量学空间后方交会实验报告
摄影测量学实验报告实验一、单像空间后方交会学院:建测学院班级:测绘082姓名:肖澎学号: 15一.实验目的1.深入了解单像空间后方交会的计算过程;2.加强空间后方交会基本公式和误差方程式,法线方程式的记忆;3.通过上机调试程序加强动手能力的培养。
二.实验原理以单幅影像为基础,从该影像所覆盖地面范围内若干控制点和相应点的像坐标量测值出发,根据共线条件方程,求解该影像在航空摄影时刻的相片外方位元素。
三.实验内容1.程序图框图2.实验数据(1)已知航摄仪内方位元素f=153.24mm,Xo=Yo=0。
限差0.1秒(2)已知4对点的影像坐标和地面坐标:3.实验程序using System;using System.Collections.Generic;using System.Linq;using System.Text;namespace ConsoleApplication3{class Program{static void Main(){//输入比例尺,主距,参与平参点的个数Console.WriteLine("请输入比例尺分母m:\r");string m1 = Console.ReadLine();double m = (double)Convert.ToSingle(m1);Console.WriteLine("请输入主距f:\r");string f1 = Console.ReadLine();double f = (double)Convert.ToSingle(f1);Console.WriteLine("请输入参与平差控制点的个数n:\r");string n1 = Console.ReadLine();int n = (int)Convert.ToSingle(n1);//像点坐标的输入代码double[] arr1 = new double[2 * n];//1.像点x坐标的输入for (int i = 0; i < n; i++){Console.WriteLine("请输入已进行系统误差改正的像点坐标的x{0}值:\r", i+1);string u = Console.ReadLine();for (int j = 0; j < n; j += 2){arr1[j] = (double)Convert.ToSingle(u);}}//2.像点y坐标的输入for (int i = 0; i < n; i++){Console.WriteLine("请输入已进行系统误差改正的像点坐标的y{0}值:\r", i+1);string v = Console.ReadLine();for (int j = 1; j < n; j += 2){arr1[j] = (double)Convert.ToSingle(v);}}//控制点的坐标输入代码double[,] arr2 = new double[n, 3];//1.控制点X坐标的输入for (int j = 0; j < n; j++){Console.WriteLine("请输入控制点在地面摄影测量坐标系的坐标的X{0}值:\r", j+1);string u = Console.ReadLine();arr2[j , 0] = (double)Convert.ToSingle(u);}//2.控制点Y坐标的输入for (int k = 0; k < n; k++){Console.WriteLine("请输入控制点在地面摄影测量坐标系的坐标的Y{0}值:\r", k+1);string v = Console.ReadLine();arr2[k , 1] = (double)Convert.ToSingle(v);}//3.控制点Z坐标的输入for (int p =0; p < n; p++){Console.WriteLine("请输入控制点在地面摄影测量坐标系的坐标的Z{0}值:\r", p+1);string w = Console.ReadLine();arr2[p , 2] = (double)Convert.ToSingle(w);}//确定外方位元素的初始值//1.确定Xs的初始值:double Xs0 = 0;double sumx = 0;for (int j = 0; j < n; j++){double h = arr2[j, 0];sumx += h;}Xs0 = sumx / n;//2.确定Ys的初始值:double Ys0 = 0;double sumy = 0;for (int j = 0; j < n; j++){double h = arr2[j, 1];sumy += h;}Ys0 = sumy / n;//3.确定Zs的初始值:double Zs0 = 0;double sumz = 0;for (int j = 0; j <= n - 1; j++){double h = arr2[j, 2];sumz += h;}Zs0 = sumz / n;doubleΦ0 = 0;doubleΨ0 = 0;double K0 = 0;Console.WriteLine("Xs0,Ys0,Zs0,Φ0,Ψ0,K0的值分别是:{0},{1},{2},{3},{4},{5}", Xs0, Ys0, Zs0, 0, 0, 0);//用三个角元素的初始值按(3-4-5)计算各方向余弦值,组成旋转矩阵,此时的旋转矩阵为单位矩阵I:double[,] arr3 = new double[3, 3];for (int i = 0; i < 3; i++)arr3[i, i] = 1;}double a1 = arr3[0, 0]; double a2 = arr3[0, 1]; double a3 = arr3[0, 2];double b1 = arr3[1, 0]; double b2 = arr3[1, 1]; double b3 = arr3[1, 2];double c1 = arr3[2, 0]; double c2 = arr3[2, 1]; double c3 = arr3[2, 2];/*利用线元素的初始值和控制点的地面坐标,代入共线方程(3-5-2),* 逐点计算像点坐标的近似值*///1.定义存放像点近似值的数组double[] arr4 = new double[2 * n];//----------近似值矩阵//2.逐点像点坐标计算近似值//a.计算像点的x坐标近似值(x)for (int i = 0; i < 2 * n; i += 2){for (int j = 0; j < n; j++){arr4[i] = -f * (a1 * (arr2[j, 0] - Xs0) + b1 * (arr2[j, 1] - Ys0) + c1 * (arr2[j, 2] - Zs0)) / (a3 * (arr2[j, 0] - Xs0) + b3 * (arr2[j, 1] - Ys0) + c3 * (arr2[j, 2] - Zs0)); }}//b.计算像点的y坐标近似值(y)for (int i = 1; i < 2 * n; i += 2){for (int j = 0; j < n; j++){arr4[i] = -f * (a2 * (arr2[j, 0] - Xs0) + b2 * (arr2[j, 1] - Ys0) + c2 * (arr2[j, 2] - Zs0)) / (a3 * (arr2[j, 0] - Xs0) + b3 * (arr2[j, 1] - Ys0) + c3 * (arr2[j, 2] - Zs0)); }}//逐点计算误差方程式的系数和常数项,组成误差方程:double[,] arr5 = new double[2 * n, 6]; //------------系数矩阵(A)//1.计算dXs的系数for (int i = 0; i < 2 * n; i += 2){arr5[i, 0] = -1 / m; //-f/H == -1/m}//2.计算dYs的系数for (int i = 1; i < 2 * n; i += 2){arr5[i, 1] = -1 / m; //-f/H == -1/m}//3.a.计算误差方程式Vx中dZs的系数for (int i = 0; i < 2 * n; i += 2)arr5[i, 2] = -arr1[i] / m * f;}//3.b.计算误差方程式Vy中dZs的系数for (int i = 1; i < 2 * n; i += 2){arr5[i, 2] = -arr1[i] / m * f;}//4.a.计算误差方程式Vx中dΦ的系数for (int i = 0; i < 2 * n; i += 2){arr5[i, 3] = -f * (1 + arr1[i] * arr1[i] / f * f);}//4.a.计算误差方程式Vy中dΦ的系数for (int i = 1; i < 2 * n; i += 2){arr5[i, 3] = -arr1[i - 1] * arr1[i] / f;}//5.a.计算误差方程式Vx中dΨ的系数for (int i = 0; i < 2 * n; i += 2){arr5[i, 4] = -arr1[i] * arr1[i + 1] / f;}//5.b.计算误差方程式Vy中dΨ的系数for (int i = 1; i < 2 * n; i += 2){arr5[i, 4] = -f * (1 + arr1[i] * arr1[i] / f * f);}//6.a.计算误差方程式Vx中dk的系数for (int i = 0; i < 2 * n; i += 2){arr5[i, 5] = arr1[i + 1];}//6.b.计算误差方程式Vy中dk的系数for (int i = 1; i < 2 * n; i += 2){arr5[i, 5] = -arr1[i - 1];}//定义外方位元素组成的数组double[] arr6 = new double[6];//--------------------外方位元素改正数矩阵(X)//定义常数项元素组成的数组double[] arr7 = new double[2 * n];//-----------------常数矩阵(L)//计算lx的值for (int i = 0; i < 2 * n; i += 2)arr7[i] = arr1[i] - arr4[i]; //将近似值矩阵的元素代入}//计算ly的值for (int i = 1; i <= 2 * (n - 1); i += 2){arr7[i] = arr1[i] - arr4[i]; //将近似值矩阵的元素代入}/* 对于所有像点的坐标观测值,一般认为是等精度量测,所以权阵P为单位阵.所以X=(ATA)-1ATL *///1.计算ATdouble[,] arr5T = new double[6, 2 * n];for (int i = 0; i < 6; i++){for (int j = 0; j < 2 * n; j++){arr5T[i, j] = arr5[j, i];}}//A的转置与A的乘积,存放在arr5AA中double[,] arr5AA = new double[6, 6];for (int i = 0; i < 6; i++){for (int j = 0; j < 6; j++){arr5AA[i, j] = 0;for (int l = 0; l < 2 * n; l++){arr5AA[i, j] += arr5T[i, l] * arr5[l, j];}}}nijuzhen(arr5AA);//arr5AA经过求逆后变成原矩阵的逆矩阵//arr5AA * arr5T存在arr5AARATdouble[,] arr5AARAT = new double[6, 2 * n];for (int i = 0; i < 6; i++){for (int j = 0; j < 2 * n; j++){arr5AARAT[i, j] = 0;for (int p = 0; p < 6; p++){arr5AARAT[i, j] += arr5AA[i, p] * arr5T[p, j];}}}//计算arr5AARAT x L,存在arrX中double[] arrX = new double[6];for (int i = 0; i < 6; i++){for (int j = 0; j < 1; j++){arrX[i] = 0;for (int vv = 0; vv < 6; vv++){arrX[i] += arr5AARAT[i, vv] * arr7[vv];}}}//计算外方位元素值double Xs, Ys, Zs, Φ, Ψ, K;Xs = Xs0 + arrX[0];Ys = Ys0 + arrX[1];Zs = Zs0 + arrX[2];Φ = Φ0 + arrX[3];Ψ = Ψ0 + arrX[4];K = K0 + arrX[5];for (int i = 0; i <= 2; i++){Xs += arrX[0];Ys += arrX[1];Zs += arrX[2];Φ += arrX[3];Ψ += arrX[4];K += arrX[5];}Console.WriteLine("Xs,Ys,Zs,Φ,Ψ,K的值分别是:{0},{1},{2},{3},{4},{5}", Xs0, Ys0, Zs0, Φ, Ψ, K);Console.Read();}//求arr5AA的逆矩public static double[,] nijuzhen(double[,] a) {double[,] B = new double[6, 6];int i, j, k;int row = 0;int col = 0;double max, temp;int[] p = new int[6];for (i = 0; i < 6; i++){p[i] = i;B[i, i] = 1;}for (k = 0; k < 6; k++){//找主元max = 0; row = col = i;for (i = k; i < 6; i++){for (j = k; j < 6; j++){temp = Math.Abs(a[i, j]);if (max < temp){max = temp;row = i;col = j;}}}//交换行列,将主元调整到k行k列上if (row != k){for (j = 0; j < 6; j++){temp = a[row, j];a[row, j] = a[k, j];a[k, j] = temp;temp = B[row, j];B[row, j] = B[k, j];B[k, j] = temp;i = p[row]; p[row] = p[k]; p[k] = i; }if (col != k){for (i = 0; i < 6; i++){temp = a[i, col];a[i, col] = a[i, k];a[i, k] = temp;}}//处理for (j = k + 1; j < 6; j++){a[k, j] /= a[k, k];}for (j = 0; j < 6; j++){B[k, j] /= a[k, k];a[k, k] = 1;}for (j = k + 1; j < 6; j++){for (i = 0; j < k; i++){a[i, j] -= a[i, k] * a[k, j];}for (i = k + 1; i < 6; i++){a[i, j] -= a[i, k] * a[k, j];}}for (j = 0; j < 6; j++){for (i = 0; i < k; i++){B[i, j] -= a[i, k] * B[k, j];}for (i = k + 1; i < 6; i++){B[i, j] -= a[i, k] * B[k, j];}for (i = 0; i < 6; i++) {a[i, k] = 0;a[k, k] = 1;}}//恢复行列次序for (j = 0; j < 6; j++){for (i = 0; i < 6; i++) {a[p[i], j] = B[i, j]; }}for (i = 0; i < 6; i++){for (j = 0; j < 6; j++) {a[i, j] = a[i, j];}}return a;}4.实验结果四.实验总结此次实验让我深入了解单像空间后方交会的计算过程,加强了对空间后方交会基本公式和误差方程式,法线方程式的记忆。
空间后方—前方交会的原理
空间后方—前方交会的原理
空间后方交会和前方交会都是航天技术中的重要概念。
空间后方交会是通过地面控制点和像片上的像点来确定单幅影像的外方位元素的方法。
这种方法以单张像片为基础,求解每张像片的外方位元素,确定被摄物体与航摄像片的几何关系。
后方交会需要地面控制点,依靠单张像片就可求出其外方位元素。
而前方交会则是利用立体像对的同名像点求解地面点坐标的过程。
它利用同名像点坐标和像对的相对方位元素解算模型点坐标。
这种方法在确定了被摄物体与航摄像片的几何关系后,利用立体像对同名像点对对相交的原理确定物方空间坐标。
如需了解更多关于空间后方—前方交会原理的信息,建议查阅航天领域相关书籍或咨询专业人士。
空间后方—前方交会的原理
空间后方—前方交会的原理
以空间后方—前方交会的原理为题,我来为大家描述一下。
空间后方—前方交会是一种用于确定目标位置的方法,常用于航空、导航、测绘等领域。
它利用人眼的立体视觉和视差效应,通过观察目标在不同视角下的位置变化,来推断目标的实际位置。
这种方法可以较精确地确定目标的距离和方位,尤其适用于远距离观测。
在进行空间后方—前方交会时,我们首先需要选择两个观测点,它们之间的距离应足够远,以便产生明显的视差效应。
然后,我们分别在这两个观测点上观察目标,并记录下目标在两个观测点的位置。
接下来,我们需要测量观测点之间的距离,并确定观测点与目标之间的夹角。
这些数据将用于计算目标的实际位置。
通过对两个观测点的位置和距离进行几何分析,我们可以得到目标相对于观测点的位移向量。
然后,我们再将这个位移向量与观测点之间的夹角结合起来,就可以计算出目标相对于观测点的实际位置。
空间后方—前方交会的原理基于视差效应,即当我们观察远处的目标时,由于两只眼睛的视角不同,目标在两只眼睛中的位置也会有所不同。
通过比较这两个位置的差异,我们就可以推断出目标的实际位置。
总的来说,空间后方—前方交会是一种利用视差效应来确定目标位
置的方法。
它可以在远距离观测中提供较为准确的测量结果,具有广泛的应用前景。
五上、数字摄影测量学单片空间后方交会
总误差方程
法方程
V Ax L
x (AT A) 1 (AT L)
X s Ys V1 A1 l1 Z V2 A2 l2 s V , A , L , x , Vn An ln T T li xi ( xi ) yi ( yi ) , Vi v xi v yi a11 a12 a13 a14 a15 a16 Ai a21 a22 a23 a24 a25 a26
已知点必须多余点, 数据处理方法采用 最小二乘法!
这是所有测量的一个统一的基本原则! 摄影测量也不例外。
二、误差方程与法方程
已知值 x0 , y0 , f ,m, X, Y, Z 观测值 x , y 相应改正数 vx,vy 未知数 Xs, Ys, Zs, , , 泰勒级数展开
四、空间后方交会的精度
求解各未知数的精度可以通过法方程系数矩阵 求逆的方法,解出相应的权倒数 Qii
mi m0 Qii 按下式计算第i未知数的中误差:
式中,m0为单位权中误差,计算公式 为: m [VV ] 0 2n 6 ,其中n为控制点的点数。
空间后方交会用到的已知点越多,空间后方交会 的精度越高,此外空点的分布也空间后方交会计算 的精度。空间后方交会使用的控制点应当避免位于 一个圆柱面上,否则,会出现解不唯一的情况。
偏导数 1
x f X Z 2 ( Z X) X s Z X s X s f 2 ( a1Z a3 X ) Z 1 X (a1 f f a3 ) Z Z 1 (a1 f a3 x) Z
偏导数 2
x f X Z 2 ( Z X) Z
空间后方—前方交会的原理
空间后方—前方交会的原理空间后方—前方交会是一种导航技术,通过测量目标物体在不同观测点的角度,并利用三角测量原理计算目标物体的位置。
这种技术广泛应用于航空、航天、导航等领域,可以帮助人们准确地确定目标物体的位置和方向。
在空间后方—前方交会中,观测者需要站在不同的位置观测目标物体,并测量目标物体相对于观测点的角度。
观测者需要使用测角仪或其他测量工具来测量角度,并记录下相应的数据。
测量完所有观测点的角度后,观测者需要根据这些角度数据进行三角计算,以确定目标物体的位置。
三角计算是利用三角函数来计算角度和边长的关系,通过已知的角度和边长来计算未知的角度和边长。
观测者需要根据测量得到的角度数据和观测点之间的距离,使用三角函数计算目标物体的位置坐标。
空间后方—前方交会的原理简单明了,但在实际应用中需要考虑一些因素。
首先,观测者需要选择合适的观测点,观测点的位置应尽量避免遮挡物,以确保观测到目标物体的角度准确无误。
其次,观测者需要准确测量角度,并尽量避免误差的产生。
最后,观测者需要进行精确的三角计算,以确保计算出的目标物体位置准确无误。
空间后方—前方交会技术的应用非常广泛。
在航空领域,飞行员可以利用该技术确定飞机的位置和方向,以确保飞行安全。
在航天领域,航天员也可以利用该技术确定航天器的位置和方向,以实现精确的轨道控制。
此外,该技术还可以应用于导航系统中,帮助人们准确导航和定位。
空间后方—前方交会是一种通过测量目标物体在不同观测点的角度,并利用三角计算原理确定目标物体位置的导航技术。
该技术在航空、航天、导航等领域有着广泛的应用,可以帮助人们准确地确定目标物体的位置和方向。
通过合理的观测点选择、准确的角度测量和精确的三角计算,空间后方—前方交会技术可以为人们提供准确可靠的导航和定位服务。
单像空间后方交会原理
单像空间后方交会原理你知道单像空间后方交会吗?这可是摄影测量里一个超有趣的概念呢!咱们先来说说啥是单像空间后方交会。
想象一下,你拿着相机拍了一张照片,这张照片里有好多好多的景物。
那单像空间后方交会呢,就是通过这一张照片里的信息,去算出拍摄这张照片的时候,相机在空间里的位置和姿态。
比如说,照片里有一座山,还有一条河,还有几棵大树。
那咱们怎么通过这些东西来知道相机当时在哪,朝哪个方向呢?这就用到单像空间后方交会啦!这当中有几个关键的东西哦。
一个是控制点,就好像是我们的“小帮手”。
这些控制点是我们事先知道它们在空间里准确位置的点。
比如说,有个特别明显的大石头,我们知道它在地球上的坐标是多少。
然后呢,还有像片的内方位元素。
这就像是相机的“小秘密”,比如说相机的焦距啦等等。
那怎么通过这些来算出相机的位置和姿态呢?这就像是一个解谜的过程!咱们得先把照片上控制点的像点坐标找出来,这就像是在照片里给这些控制点“定位”。
然后呢,根据一些数学公式和算法,把这些坐标啊、内方位元素啊、控制点的空间坐标啊等等都放到一起,就像是把一堆拼图的碎片拼起来。
这个过程可不容易哦,得算好多好多的数学式子。
但是别担心,咱们聪明的科学家们早就想出了办法,有各种软件和工具能帮咱们完成这些复杂的计算。
你可能会想,这有啥用啊?用处可大啦!比如说,我们要做地图,要对一个地方进行测量,单像空间后方交会就能帮我们得到相机的位置和姿态,这样就能更准确地知道照片里的东西在实际空间里的位置啦。
而且哦,现在科技越来越发达,单像空间后方交会的精度也越来越高。
这就像是我们的眼睛越来越厉害,能看得更清楚,更准确!想象一下,如果没有单像空间后方交会,那我们看到的照片就只是一张好看的图片,没办法知道那么多背后的信息。
但是有了它,一张照片就像是一个装满了秘密的宝盒,我们可以一点点地解开,发现更多有趣的东西。
怎么样,是不是觉得单像空间后方交会很神奇很有趣呀?希望我讲得能让你明白这个有点复杂但又超级酷的原理!。
单像空间后方交会和双像解析空间后方-前方交会的算法程序实现
单像空间后方交会和双像解析空间后方-前方交会的算法程序实现遥感科学与技术摘要:如果已知每张像片的6个外方位元素,就能确定被摄物体与航摄像片的关系。
因此,利用单像空间后方交会的方法,可以迅速的算出每张像片的6个外方位元素。
而前方交会的计算,可以算出像片上点对应于地面点的三维坐标。
基于这两点,利用计算机强大的运算能力,可以代替人脑快速的完成复杂的计算过程。
关键词:后方交会,前方交会,外方位元素,C++编程0.引言:单张像片空间后方交会是摄影测量基本问题之一,是由若干控制点及其相应像点坐标求解摄站参数(X S,Y S,ZS,ψ、ω、κ)。
单像空间后方交会主要有三种方法:基于共线条件方程的平差解法、角锥法、基于直接线性变换的解法。
而本文将介绍第一种方法,基于共线条件方程反求象片的外方位元素。
而空间前方交会先以单张像片为单位进行空间后方交会,分别求出两张像片的外方位元素,再根据待定点的一对像点坐标,用空间前方交会的方法求解待定点的地面坐标。
可以说,这种求解地面点的坐标的方法是以单张像片空间后方交会为基础的,因此,单张像片空间后方交会成为解决这两个问题以及算法程序实现的关键。
1.单像空间后方交会的算法程序实现:(1)空间后方交会的基本原理:对于遥感影像,如何获取像片的外方位元素,一直是摄影测量工作者探讨的问题,其方法有:利用雷达(Radar)、全球定位系统(GPS)、惯性导航系统(I N S)以及星像摄影机来获取像片的外方位元素;也可以利用一定数量的地面控制点,根据共线方程,反求像片的外方位元素,这种方法称为单像空间后方交会(如图1所示)。
图中,地面坐标X i、Yi、Zi和对应的像点坐标x i、yi是已知的,外方位元素XS、Y S、ZS(摄站点坐标),ψ、ω、κ(像片姿态角)是待求的。
(2)空间后方交会数学模型:空间后方交会的数学模型是共线方程, 即中心投影的构像方程:式中X、Y、Z是地面某点在地面摄影测量坐标系中的坐标,x,y是该地面点在像片上的构像点的像片坐标,对于空间后方交会而言它们是已知的,还有主距f是已知的。
后方交会的原理
后方交会的原理后方交会是测量学中一种常用的方法,可以用来确定点的位置,尤其是在地图制图、土地测量、建筑设计等领域中。
本文将从定义、原理、步骤、误差分析等方面详细介绍后方交会的原理。
一、定义后方交会是一种测量学方法,通过已知控制点的坐标、方位角等信息,对未知点进行定位。
它是一种基于三角测量原理的方法,通过三角形相似原理计算出未知点的坐标。
后方交会可以用来测量平面坐标、高程、方位角等,适用于各种地形和地貌。
二、原理后方交会的原理基于三角测量原理,即以三个已知点为顶点的三角形,其内部所有角度和三个顶点的距离都是已知的。
通过测量三角形内角和三边长度,可以计算出未知点的坐标。
具体来说,后方交会的原理可以分为以下几个步骤:1.测量控制点的坐标和方位角首先需要测量控制点的坐标和方位角,以确定已知点的位置和方向。
控制点可以是已知点,也可以是通过前方交会计算出的点。
2.测量未知点与控制点之间的距离和方位角接下来需要测量未知点与控制点之间的距离和方位角。
距离可以通过测距仪、激光测距仪等设备进行测量,方位角可以通过经纬仪或罗盘进行测量。
3.计算三角形内角和三边长度通过测量三角形内角和三边长度,可以计算出未知点的坐标。
其中,内角可以通过三角函数计算,三边长度可以通过勾股定理计算。
4.检验后方交会的精度最后需要对后方交会的精度进行检验,可以通过误差分析等方法进行评估。
如果误差较大,可以通过增加控制点、提高测量精度等方式进行调整。
三、步骤后方交会的步骤可以总结为以下几个:1.确定控制点的位置和方向,测量其坐标和方位角。
2.测量未知点与控制点之间的距离和方位角。
3.计算三角形内角和三边长度,确定未知点的坐标。
4.检验后方交会的精度,评估误差。
四、误差分析后方交会的精度受到多种因素的影响,如测量仪器精度、控制点布设精度、天气环境等。
误差分析可以通过以下几种方法进行: 1.重复测量法重复测量法是通过多次测量同一点来评估测量精度的方法。
9-空间后方交会
a11 a12 Ai a21 a22
a13 a23
a14 a24
a15 a25
a16 a26
T
X dXS
dYS
dZS
d d d
li l x
ly
T
vi vx
vy
T
把所有像控点的误差方程列出后,构成总 误差方程,根据最小二乘间接评差原理可 列出法方程式:
五、空间后方交会实践
如何获取像片的六个外方位元素?
1)利用雷达;全球定位系统GPS;惯性导航系统.
2)空间后方交会:利用一定数量的地面控制点, 根据共线方程,反求像片的外方位元素。(已知 像片的内方位元素,至少三个地面点坐标并测出 相应的像点坐标)
计算要点:
1)计算的数学模型:共线方程按泰勒级数展开, 取一次项(线性化)。
2)在像片的四角选取四个或更多地面控制点,利 用最小二乘法平差计算。
计算步骤: 1)获取已知数据:比例尺1/m;H;内定向;控制点 地面坐标。
2)测量控制点的像点坐标。标刺,测量像框坐 标;像主点改正。 3)确定未知参数的初始值:竖直摄影时,角元 素初始值为零;线元素中,Zs0=H=mf; Xs0,Ys0取 控制点坐标的均值。
投 影 中 心 的 系 数
二、线性化-续
X Y Z
x xs 1 1 R R R y y s z zs
1 1
其中,R
R R R
1
1
1
把各偏导数代入整理得
f XX b2 Z Z ZZ x XX f sin XY f cos a 15 fsin ZZ ZZ x Yf a 16 Z fX f b1 YY f b2 XY y b 3 a 24 f b1 Z ZZ ZZ y XY f sin YY f cos a 25 fcos ZZ ZZ y X f a 26 Z x a 14 f Yf
空间后方交会基本原理
F y X F y S 0 d S X F Y y S 0 d S Y F Z y S 0 d S Z F y 0 d F y 0 d F y 0 d F y 0
因为
Fx xf
a1(XXS)b1(YYS)c1(ZZS) 0 a3(XXS)b3(YYS)c3(ZZS)
④在摄影过程中直接获取。
h
4
内
• 单像空间后方交会概述
容
• 共线方程的线性化(难点) • 利用共线条件方程解算像片的外方位元
素(重点)
安
排
[一]概述
1、什么叫单像空间后方交会 利用地面控制点及其在片像上的像点,确定一
张像片外方位元素的方法。
S (XS 、 YS 、 ZS)
c
Z
Y
b
a
C
A X 地面控B制点(Ground Control Point, GCP)
第十讲 空间后方交会
已学过的主要内容
绪论 摄影与航空摄影 第一部分 单张航摄像片的解析 第二部分 立体像对的基本知识 第三部分 作业基本理论
问题的引出
XS,YS,ZS
YXYXss((ZZZZs)s)ba1cx11xxb2acy22yyb3ca3f3ff
若 ai,bi,ci
f x,y ,Z
已知
X,Y
h
6
[一]概述
1、什么叫单像空间后方交会 利用地面控制点及其在像片上的像点,确定一
张像片外方位元素的方法。 2、单像空间后方交会的基本方法
角锥体法
S
cb a
ZT
C B
YT
A
D
XT
h
7
角锥体法
S
c
b
空间后方交会实验报告
空间后方交会实验报告1. 引言空间后方交会是一种常用的测量手段,用于确定目标在空间中的坐标位置。
本实验旨在通过空间后方交会实验,了解空间测量的基本原理和方法,并锻炼我们的观察、计算和推理能力。
2. 实验目的- 学习和掌握空间后方交会实验的基本原理和步骤;- 熟悉使用测量仪器和处理数据的方法;- 锻炼团队合作和问题解决能力。
3. 实验原理空间后方交会利用多个测量设备,通过测量目标物体在不同位置上的观测数据,计算出目标物体的空间坐标。
主要基于以下原理:- 多点定位原理:通过多次测量目标物体在不同位置上的观测数据,可以确定目标物体的坐标位置;- 观测误差校正原理:由于观测仪器的误差存在,需要对观测数据进行误差校正,以提高测量精度。
4. 实验步骤4.1 实验准备- 确定实验区域,搭建固定的测量基准点;- 使用全站仪对测量基准点进行校准,确保测量精度;- 准备至少三个可移动的晶体棱镜,用于固定在目标物体上;- 准备笔记本电脑和相应的测量软件,用于数据处理和计算。
4.2 实验操作1. 将晶体棱镜固定在目标物体上,并确定测量的起始位置;2. 使用全站仪测量目标物体在不同位置上的观测数据,包括水平角、垂直角和斜距;3. 将观测数据导入计算软件,并进行误差校正;4. 根据误差校正后的观测数据,计算目标物体的空间坐标;5. 重复步骤2至4,获得更多的观测数据,以提高计算精度;6. 对多次观测结果进行平均计算,得到最终的目标物体空间坐标;7. 将计算结果进行分析和评估,确定测量精度和可靠度。
5. 实验结果与讨论经过实验操作和数据处理,获得了目标物体的空间坐标。
通过对测量结果的分析和比较,可以得出以下结论:- 目标物体的空间坐标可以通过空间后方交会实验进行测量和计算;- 观测数据的误差校正对测量结果的精度和可靠度具有重要影响;- 重复观测和平均计算可以提高测量结果的准确性。
6. 实验总结空间后方交会实验是一种常用的测量手段,可以用于测量目标物体在空间中的坐标位置。
后方交会法原理
后方交会法原理在测量领域中,后方交会法是一种常用的方法,它可以用于确定已知控制点的位置以及未知点的坐标。
后方交会法是基于三角形相似原理的测量技术,它可以通过对已知控制点的距离和方位角度的测量,来计算出未知点的坐标。
本文将详细介绍后方交会法的原理和应用。
一、后方交会法原理后方交会法是基于三角形相似原理的测量技术,它是通过已知控制点的位置和方位角度的测量,来计算未知点的坐标。
具体来说,后方交会法的原理包括以下几个步骤:1.测量控制点的位置和方位角度在后方交会法中,需要测量已知控制点的位置和方位角度。
控制点是已知坐标的点,通常是在测量区域的边界或者是在地面上的显著点。
测量控制点位置的方法有很多种,包括全站仪、GPS等。
方位角度是指测量点相对于一个已知方向的角度,可以通过全站仪或者经纬仪等测量仪器来测量。
2.测量未知点与控制点之间的距离和角度在已知控制点的基础上,需要测量未知点与控制点之间的距离和角度。
距离可以通过测距仪等仪器来测量,角度可以通过全站仪等仪器来测量。
需要注意的是,测量时需要保证控制点与未知点之间的视线畅通,以确保测量的准确性。
3.计算未知点的坐标通过已知控制点的坐标、方位角度和未知点与控制点之间的距离和角度,可以利用三角形相似原理来计算未知点的坐标。
具体计算公式如下:X = X0 + L * sin(A + α)Y = Y0 + L * cos(A + α)其中,X0和Y0是已知控制点的坐标,A是控制点与未知点之间的方位角度,α是未知点与控制点之间的角度,L是未知点与控制点之间的距离。
通过以上公式,可以计算出未知点的坐标。
二、后方交会法的应用后方交会法在测量领域中应用广泛,可以用于确定已知控制点的位置以及未知点的坐标。
以下是后方交会法的一些应用场景:1.地形测量在地形测量中,后方交会法可以用于确定山顶、山脚、河流等地形特征点的坐标。
通过测量已知控制点的位置和方位角度,以及未知点与控制点之间的距离和角度,可以计算出未知点的坐标,从而确定地形特征点的位置。
空间后方交会
空间后方交会
科技名词定义
中文名称:空间后方交会
英文名称:space resection
定义:利用航摄像片上三个以上不在一条直线上的已知点按构像方程计算该像片外方位元
素的方法。
空间后方交会,是指利用航摄像片上三个以上不在一条直线上的控制点按构像方程计算该像片外方位元素的方法。
内容
仅在待定点上设站,向三个已知控制点观测两个水平夹角a、b,从而计算待定点的坐标,称为后方交会。
交会测量是加密控制点常用的方法,它可以在数个已知控制点上设站,分别向待定点观测方向或距离,也可以在待定点上设站向数个已知控制点观测方向或距离,而后计算待定点的坐标。
常用的交会测量方法有前方交会、后方交会、侧边交会和自由测站法。
后方交会法首先出现于测绘地形图工作中,测量上称为“三点题”,是用图解法作为加密图根点之用。
后来随着解析法、公式法的出现,在工程建设控制测量中也经常被采用。
比如隧道工程控制网往往由于隧道开工前测设完成,而洞口土石方施工完毕后,需补设洞口投点,以便控制隧道轴线,测设投点就要用到后方交会法;深水桥墩放样测量中的墩心定位也可以应用此法,还可用来测定施工控制导线的始终点等。
应用范围之广说明了此法的实用性很强。
其代表图形如下图所示。
图中三角形ABC是控制网中的一个三角形,P点即为后方交会点(需确定坐标的待定点),只要置棱镜于P点,用全圆测回法测定a,b,r 三个角值,即可应用解析法公式算出待定点P的坐标。
此法内外业工作量小,只要P点的点位精度符合施工放样要求或作为洞口投点的精度要求,就可以成为广大测绘科技工作者所乐意选用的方法之一。
空间后方交会原理
式2-6,2-7其中:X 和Y 是控制点的地面坐标;和y'是控制点在摄影测量坐标系下的坐标值。
由此可以解算出参数X Y a,b,T ,T ,据此确定参数的初始值:
κ式2-8,注意:κ 值的确定采用的是全圆反正切,这和初等数学中的反正切图像不一样。该函数是计算机系统中的一个数学函数,将最后所得值归算到-180°到+180°之间。
文章来自中国硕士论文网,更多工程硕士论文就在中国硕士论文网,如需转载请保留一个链接:
参考文献
张剑清等著.摄影测量学[M].武汉:武汉大学出版社
李德仁等著.摄影测量与遥感概论[M].北京:测绘出版社
张凤举,张华海等编著.控制测量学[M].北京:煤炭工业出版社出版社
计算方向余弦矩阵方向余弦矩阵的形式如下:
式方向余弦的各个元素的计算方法求解影像的外方位元素的参数根据控制点的个数和坐标值,确定法方程,以其中一点A 为例,解释如下:各个系数的值即是2-2 中的各个改正数的系数值。
如果有四个已知的控制点,那么矩阵B就是8× 6的矩阵。使用四个以上的控制点,可以解出Δ的值(Δ = (BT B)?1 (BTε ))。然后重复步骤和1.3 进行一次迭代运算,即可完成单像空间后方交会的工作。
空间后方交会名词解释
空间后方交会名词解释
空间后方交会,是指利用航摄像片上三个以上不在一条直线上的控制点按共线方程计算该像片外方位元素的方法。
是单幅影像解析过程中的一个步骤。
如果我们知道每幅影像的6个外方位元素,就能确定被摄物体与航摄影像的关系。
因此,如何获取影像的外方位元素,一直是摄影测量工作者所探讨的问题。
可采取的方法有:利用雷达、全球定位系统(GPS)、惯性导航系统(INS)以及星相摄影机来获取影像的外方位元素;也可利用影像覆盖范围内一定数量的控制点的空间坐标与影像坐标,根据共线条件方程反求该影像的外方位元素,这种方法称为单幅影像的空间后方交会。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
空间后方交会原理.txtゅ你不用一上线看见莪在线,就急着隐身,放心。
莪不会去缠你。
说好的不离不弃现在反而自己却做不到╮单幅航空影像空间后方交会程序设计
来源:
摘要:航空影像单像空间后方交会是利用影像所覆盖地面范围内若干控制点的已知地面坐标和相应点的像点坐标,根据共线方程,在最小二乘法的原则下,计算航空影像的外方位元素。
本文采用 编译平台和C#语言编写了单像空间后方交会的参数计算程序。
利用相关文献中的实验数据与结果,证明本程序运行良好,参数解算精确,可以满足工程项目要求。
关键词:摄影测量与遥感;外方位元素;共线方程;最小二乘原则;空间后方交会;
0 引言
确定影像或摄影光束在摄影瞬间的空间位置和姿态的参数,称为影像的外方位元素。
一幅影像的外方位元素包括6 个参数,其中3 个参数线参数,用于描述摄影中心S 相对于物方空间坐标系的位置(X,Y,Z);另外3 个是角元素,用于描述影像面在摄影瞬间的空中姿态。
角元素有三种不同的表达形式:(1)以Y 轴为主轴的φ-ω-κ系统(主轴是在旋转过程中,空间方向不变的一个固定轴):以Y 为主轴旋转φ角,然后绕X 轴旋转ω角,最后围绕Z 轴旋转κ角。
(2)以X 轴为主轴的φ?-ω?-κ?系统:以X 为主轴旋转ω?角,然后绕轴旋转φ?角,最后围绕Z 轴旋转κ?角。
(3)以Z 轴为主轴的A-α-κ系统:以Z 为主轴旋转A 角,然后绕Y 轴旋转α角,最后围绕Z 轴旋转κ角。
本文的角元素系统选择使用第一种表达方式,即φ-ω-κ[1]。
如果知道了每幅影像的6 个外方位元素,就能确定被摄物体与航摄影像的关系。
因此,如何获取影像的外方位元素,一直是摄影测量工作者所探讨的问题。
目前,从技术方面来说,解决这个问题主要有两种手段:一种是利用雷达、全球定位系统(GPS),惯性导航系统(INS)以及星相摄影机来获取影像的外方位元素;另一种是利用影像覆盖范围内一定数量的控制点的空间坐标与摄像坐标,根据共线条件方程,反求该影像的外方位元素,这种方法称为单幅影像的空间后方交会。
在无法直接得到影像外方位元素的情况下,只能通过单像空间后方交会的方式来获取影像的外方位元素[2]。
本文在最小二乘法的原则下,采用 编译平台和C#语言编写了单像空间后方交会的参数计算程序。
利用相关文献中的实验数据与结果,证明本程序运行良好,参数解算精确,可以满足工程项目要求。
1 航空影像单像空间后方交会的基本思想和流程设计
单像空间后方交会的基本思想是:以单幅航空影像为基础,从该影像所覆盖地面范围内若干控制点的已知地面坐标和相应点的像坐标量测值出发,根据共线条件方程,求解该影像在航空摄影时刻的外方位元素。
由于空间后方交会所采用的数学模型共线方程是非线性函数,为了方便外方位元素的求解,需要首先对共线方程进行线性化。
共线方程的非线性形式是:
式在共线方程线性化的基础之上,设计了航空影像单像空间后方交会算法的计算流程。
式即为共线方程的线性化形式。
由于采用了泰勒公式将非线性的共线方程化成了线性的共线方程形式,因此在程序设计的时候,需要采用迭代运算的方式使运算精度达到工程应用的要求。
设置外方位元素的初始值在航空影像近似水平的情况下,可以设置外方位元素的初始值:
ω使用控制点中的两个控制点,确定的航高值作为外方位元素L Z 的初始值。
航高值H 的确定流程如下:
式±式根据计算出的航高和控制点在摄影测量坐标系下的坐标值可以计算出, ,κ L L X Y 的
初始值。
控制点在摄影测量坐标系下的坐标值计算公式为:
式使用下列公式,在最小二乘法原则下确定, ,κ L L X Y 的初始值[3]:
式2-6,2-7其中:X 和Y 是控制点的地面坐标;和y'是控制点在摄影测量坐标系下的坐标值。
由此可以解算出参数X Y a,b,T ,T ,据此确定参数的初始值:
κ式2-8,注意:κ值的确定采用的是全圆反正切,这和初等数学中的反正切图像不一样。
该函数是计算机系统中的一个数学函数,将最后所得值归算到-180°到+180°之间。
计算方向余弦矩阵方向余弦矩阵的形式如下:
式方向余弦的各个元素的计算方法求解影像的外方位元素的参数根据控制点的个数和坐标值,确定法方程,以其中一点A 为例,解释如下:各个系数的值即是式2-2 中的各个改正数的系数值。
如果有四个已知的控制点,那么矩阵B就是8× 6的矩阵。
使用四个以上的控制点,可以解出Δ的值(Δ = (BT B)?1 (BTε ))。
然后重复步骤和1.3 进行一次迭代运算,即可完成单像空间后方交会的工作。
2 单像空间后方交会的程序设计
根据以上基本理论,采用 编译平台和C#语言编写了单像空间后方交会的参数计算程序。
程序界面设计如下:点击打开按钮后,选择控制点文件,即可在dataGridView 控件中可以自动显示出控制点文件中的内容。
然后点击确定按钮完成程序运算。
函数作用:计算摄影测量地面坐标系和大地地面坐标系之间的仿射变换参数。
输入:摄影测量地面坐标系坐标值,存放在泛形结构Csrc 中;大地地面坐标系坐标值,存放在泛形结构Cdst 中。
输出:摄影测量地面坐标系和大地地面坐标系之间的仿射变换参数。
函数名称:
函数作用:组成方向余弦矩阵,返回一个3*3 的数组。
输入:外方位元素的三个值:omt , fa , kama。
输出:方向余弦矩阵。
3 实验过程与结果
使用相关文献中提供的数据对本程序进行测试,测试结果表明程序设计合理,运算结果正确,可以应用于工程项目。
根据文献[5]P238 提供的相关测试数据,对程序进行检测。
4 结论本文介绍了航空影像外方位元素的概念和单幅航空影像后方交会的基本理论,在此基础之上,利用.NET 平台和C#语言设计了参数解算程序。
设计名称为的主函数,输入是控制点的影像坐标和空间坐标,输出单幅影像的六个外方位元素。
通过相关文献中提供的实验数据进行实验,结果表明本文开发的单像后方交会参数解算程序效果良好,运算正确,可以应用于工程项目。
文章来自中国硕士论文网,更多工程硕士论文就在中国硕士论文网,如需转载请保留一个链接:
参考文献
张剑清等著.摄影测量学[M].武汉:武汉大学出版社
李德仁等著.摄影测量与遥感概论[M].北京:测绘出版社
张凤举,张华海等编著.控制测量学[M].北京:煤炭工业出版社出版社
周礼著.C#和.NET3.0 第一步[M].北京:清华大学出版社.。