基础摄影测量学_光束区域网平差法列加密区内第一次迭代的全部原始误差方程式

此为前两周课程实习的题目。

1.课程设计目的

按光束法区域网平差方法列出加密区内第一次迭代的全部原始误差方程式，掌握光束法区域网平差的基本原理。

2.数据说明

（1）影像资料：两条航线，每航线3张像片，共6张航片，保存在images文件夹；
（2）控制点资料：控制点刺点位置图片“控制点刺点图.pdf”和点地面坐标文件“control.txt”保存在文件夹control中。
（3）相机鉴定资料：
相机名：UCX30
像素大小：0.006mm
影像宽度（pixels）;11310
影像高度(pixels)：17310
数码影像边缘像元数：0
相机焦距（mm）：70.50
相机PPA主点X（mm):-0.00000
相机PPA主点Y（mm):0.00000

3.方案设计

1.像片编号
（1）像片编号：按航线从上到下，像片从左到右的顺序进行航带排列，并将像片编号标示在方便的位置。如下图所示：

2.点编号：将所有控制点和连接点按从左到右，从上到下的次序进行统一编号，并标在点的方便处。如下图所示：

3.像片连接点选取、转点
对比控制点刺点图标出控制点的位置，在ps cs6中用红色三角形表示，像点位置大概位于三角形中心；在各相片重叠区域中选取明显像点作为连接点，并在ps cs6用红圈表示，点位大约在红圈中心。在ps cs6中将所有控制点和连接点转到有重叠的临近像片上，并标出相应的点位。效果如下图所示：

4.坐标量测与内定向
逐片进行内定向，在ENVI软件中将光标逐移到该像片所有控制点和连接点上，读出像点的行列值，编写内定向程序获取像点的x,y坐标并记录。
将所量测的像点坐标变换到以像主点为原点的像平面坐标系中，该过程称为数字影像的内定向。其主要流程包括：（1）将行列值转化为像点坐标值；（2）原点的平移。

#include "pch.h"
#include <iostream>
#include <fstream>
#include<math.h>
using namespace std;
//输出m行、n列的矩阵a
void output(double *a, int m, int n)
{
	int i, j;
	for (i = 0; i < m; i++)
	{
		for (j = 0; j < n; j++)
		{
			cout << a[i*n + j] << "\t";
		}
		cout << "\n" << endl;
	}
}

int main()
{
	//读取量测的坐标
	ifstream f1("point.txt");
	double P[42][2];
	for (int i = 0; i < 42; i++)
	{
		for (int j = 0; j < 2; j++)
		{
			f1 >> P[i][j];
		}
	}
	f1.close();

//相机参数
double pixel = 0.006;//mm
double width = 11310;//pixel
double height = 17310;//pixel
double x_offset = (width*pixel) / 2;
double y_offset = -(height*pixel) / 2;
double P1[42][2] = { 0 };


//内定向
for (int i = 0; i < 42; i++)
{
	P1[i][0] = P[i][0] * pixel - x_offset;
	P1[i][1] = -P[i][1] * pixel - y_offset;

}
cout << "量测的像点坐标为：" << endl;
output(*P1, 42, 2);
system("pause");
return 0;

}
5.确定各像片外方位元素和连接点地面坐标的初值 
各张像片的外方位元素如下图所示；编写前方交会程序确定连接点坐标初始值，并将其记录。

利用立体像对两张像片的内方位元素、同名像点坐标和像对的外方位元素解算模型点坐标地面点坐标的工作就是前方交会。其工作流程如下方框图所示：

#include "pch.h"
#include <iostream>
#include <fstream>
#include