CN102147249B - 基于直线特征的星载光学线阵影像精确纠正处理方法 - Google Patents

Abstract

本发明涉及一种基于直线特征的星载光学线阵影像精确纠正处理方法，不使用控制点，直接利用遥感影像与地面的直线特征，对星载光学线阵影像进行精确纠正，具体步骤如下：1)建立模型；2)计算模型初始参数；3)计算模型精确参数；4)纠正星载光学遥感影像。与现有技术相比，本发明具有不用选取控制点实现星载遥感影像的精确几何纠正等优点。

Description

基于直线特征的星载光学线阵影像精确纠正处理方法

技术领域

本发明涉及摄影测量与遥感，尤其是涉及一种基于直线特征的星载光学线阵影像精确纠正处理方法。

背景技术

经典的遥感影像几何纠正方法，都是选取足够的对应控制点，依据相应的几何纠正模型求解成像的模型参数，进而实现遥感影像的精确几何纠正。目前的高分辨率遥感影像主要有两种成像模式：框幅式成像(包括面阵CCD影像)和线性阵列成像，高分辨率星载光学影像多为线性阵列成像。使用足够的对应控制点可以实现两种成像影像的精确几何纠正。但是，控制点需要精确选择明显的目标点、角点、道路交叉点等，通常选取困难；特别是在地物贫乏地区没有足够的特征信息时，控制点无法选取，而少量的线特征就成为遥感影像纠正的重要依据。目前，对于框幅式光学影像的线特征定向、空中三角测量、影像配准纠正方法有一些研究成果，基于线特征的星载线阵影像精确纠正文献很少，其难点主要在于：(1)线性阵列影像的每行象元的成像参数随时间变化，每条直线上象元的成像参数不同；(2)一条直线的不同象元具有不同外方位元素，使建立的直线方程与摄影测量共线模型方程产生相关性，无法正确解算出模型参数。

发明内容

本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于直线特征的星载光学线阵影像精确纠正处理方法。

本发明的目的可以通过以下技术方案来实现：

一种基于直线特征的星载光学线阵影像精确纠正处理方法，其特征在于，不使用控制点，直接利用遥感影像与地面的直线特征，对星载光学线阵影像进行精确纠正，具体步骤如下：

1)建立模型；

2)计算模型初始参数；

3)计算模型精确参数；

4)纠正星载光学遥感影像。

所述的步骤1)建立模型如下：

对于星载线阵光学影像，其直线特征上的不同点具有不同的外方位参数，基于星载遥感影像的几何特点，利用精确仿射变换模型建立像点与物点之间的几何关系，线性化后的计算模型和条件方程误差方程式为：

$\left\{\begin{array}{c}{V}_x={\mathrm{XdL}}_1+{\mathrm{YdL}}_2+{\mathrm{ZdL}}_3+{\mathrm{dL}}_4-(x-x_i)-l_x\\ V_y=\frac{1}{D}(\mathrm{Xd}{L}_5+{\mathrm{YdL}}_6+{\mathrm{ZdL}}_7+{\mathrm{dL}}_8+\frac{y}{(f-y\tan \mathrm{\ω}){\cos }^2\mathrm{\ω}}\·(\frac{Z_i\sin \mathrm{\ω}}{D}-y)\mathrm{d\ω})-(y-\frac{A}{D_i})-l_y\end{array}\right.---\left(1\right)$

其中，

l_x＝(y-y_a)×(x_b-x_a)÷(y_b-y_a)+x_a-x

l_y＝(x-x_a)×(y_b-y_a)÷(x_b-x_a)+y_a-y

$D=\frac{f-Z_i/\cos \mathrm{\ω}}{f-y\tan \mathrm{\ω}}$

(x_a，y_a)和(x_b，y_b)为确定影像上对应直线的两个端点a和b，l_x、l_b为像方直线任意点在像方x轴和y轴的投影，ω为传感器侧视角，(x，y)为影像上得到的像点坐标，(X，Y，Z)为地图上获取的高斯坐标。

所述的步骤2)计算模型初始参数如下：

21)依据像元的行坐标，从影像的星历数据中获取该行成像时刻的视方向，并利用中心行影像记录扫描时间、每行扫描间隔以及像点所在行，确定像点的姿态记录时刻，使用拉格朗日插值运算，由成像姿态数据计算所选像点成像时刻的成像位置、速度、成像姿态角：

位置与速度由下列公式计算

$\stackrel{\→}{P}\left(t\right)=\underset{j=1}{\overset{8}{\mathrm{\Σ}}}\frac{\stackrel{\→}{P}\left(t_j\right)\×\underset{i\≠j}{\overset{8}{\underset{i=1}{\mathrm{\Π}}}}(t-t_i)}{\underset{i\≠j}{\overset{8}{\underset{i=1}{\mathrm{\Π}}}}(t_j-t_i)}---\left(2\right)$

$\stackrel{\→}{V}\left(t\right)=\underset{j=1}{\overset{8}{\mathrm{\Σ}}}\frac{\stackrel{\→}{V}\left(t_j\right)\×\underset{i\≠j}{\overset{8}{\underset{i=1}{\mathrm{\Π}}}}(t-t_i)}{\underset{i\≠j}{\overset{8}{\underset{i=1}{\mathrm{\Π}}}}(t_j-t_i)}---\left(3\right)$

式中：

为卫星位置坐标，为卫星速度坐标，t_i为相应位置与速度的宇宙时间；

成像姿态角可以使用线性内插得到某一观测时间t对应影像某行的姿态值[a_p(t)，a_γ(t)，a_y(t)]，成像姿态角内插计算公式：

$a_p\left(t\right)=a_p\left(t_i\right)+(a_p\left(t_{i+1}\right)-a_p\left(t_i\right))\×\frac{t-t_i}{t_{i+1}-t_i}---\left(4\right)$

式中：a_p(t)和a_p(t_i)分别为t和t_i时刻绕俯仰轴的转角，a_γ(t)为t时刻绕滚动轴的转角，a_y(t)为t时刻绕偏轨轴的转角，a_γ(t)和a_y(t)内插计算同(4)式；

22)将轨道坐标系统的视方向转换到地心坐标系下的视方向，再在地心坐标系下求解视方向与地球椭球的交点：

计算像点在轨道坐标系统中的视方向

${\stackrel{\→}{u}}_2=\frac{{\stackrel{\→}{u}}_2^{\′}}{\left|\right|{\stackrel{\→}{u}}_2^{\′}\left|\right|}---\left(5\right)$

式中：

${\stackrel{\→}{u}}_2^{\′}=M_p\·M_r\·M_y\·{\stackrel{\→}{u}}_1$

$M_p=\left[\begin{array}{ccc}1& 0& 0\\ 0& \cos \left(a_p\left(t\right)\right)& \sin \left(a_p\left(t\right)\right)\\ 0& -\sin \left(a_p\left(t\right)\right)& \cos \left(a_p\left(t\right)\right)\end{array}\right]$

$M_r=\left[\begin{array}{ccc}\cos \left(a_r\left(t\right)\right)& 0& -\sin \left(a_r\left(t\right)\right)\\ 0& 1& 0\\ \sin \left(a_r\left(t\right)\right)& 0& \cos \left(a_r\left(t\right)\right)\end{array}\right]$

$M_y=\left[\begin{array}{ccc}\cos \left(a_y\left(t\right)\right)& -\sin \left(a_y\left(t\right)\right)& 0\\ \sin \left(a_y\left(t\right)\right)& \cos \left(a_y\left(t\right)\right)& 0\\ 0& 0& 1\end{array}\right]$

将轨道坐标系统的视方向转换到地心坐标系下的视方向

转换模型如下：

${\stackrel{\→}{u}}_3=\left[\begin{array}{ccc}{\left(a_p\left(t\right)\right)}_p& {\left(a_{\mathrm{\γ}}\left(t\right)\right)}_p& {\left(a_y\left(t\right)\right)}_p\\ {\left(a_p\left(t\right)\right)}_{\mathrm{\γ}}& {\left(a_{\mathrm{\γ}}\left(t\right)\right)}_{\mathrm{\γ}}& {\left(a_y\left(t\right)\right)}_{\mathrm{\γ}}\\ {\left(a_p\left(t\right)\right)}_y& {\left(a_{\mathrm{\γ}}\left(t\right)\right)}_y& {\left(a_y\left(t\right)\right)}_y\end{array}\right]\·{\stackrel{\→}{u}}_2---\left(6\right)$

式中：a_p(t)，a_γ(t)和a_y(t)分别为俯仰轴、滚动轴和偏轨轴。

23)将该交点坐标转换为高斯平面直角坐标，即得到像点对应的地面点空间坐标：

在地心坐标系下求解视方向与地球椭球的交点，从卫星所处位置P(t)得到视方向

可以计算在地球椭球体上方h高处视方向与椭球体上的交点(X，Y，Z)，公式如下：

$\left\{\begin{array}{c}X=X_P+\mathrm{\μ}\×{\left(u_3\right)}_p\\ Y=Y_P+\mathrm{\μ}\×{\left(u_3\right)}_{\mathrm{\γ}}\\ Z=Z_P+\mathrm{\μ}\×{\left(u_3\right)}_y\end{array}\right.---\left(7\right)$

式中：(X_p，Y_p，Z_p)为成像时刻卫星的空间位置，由式(2)求得；μ为需要求解的系数。由于椭球体上的点(X，Y，Z)满足椭球方程：

$\frac{X^2+Y^2}{A^2}+\frac{Z^2}{B^2}=1$

式中：A＝a+h，B＝b+h，a，b，h分别为椭球长半轴、短半轴和地面高程，没有地面高程模型数据时h＝0，将式(7)代入椭球方程并整理可得：

$[\frac{{\left(u_3\right)}_p^2+{\left(u_3\right)}_{\mathrm{\γ}}^2}{A^2}+\frac{{\left(u_3\right)}_y^2}{B^2}]\×{\mathrm{\μ}}^2+2\×[\frac{X_P{\left(u_3\right)}_P+Y_P{\left(u_3\right)}_{\mathrm{\γ}}}{A^2}+\frac{Z_P{\left(u_3\right)}_y}{B^2}]\×\mathrm{\μ}+[\frac{X_P^2+Y_P^2}{A^2}+\frac{Z_P^2}{B^2}]=1$

求解该一元二次方程，得到两个不同解(μ₁，μ₂)，取最小解(μ_min)为式(7)的系数，并代入式(7)可以求得地面点，最后将该坐标转换为高斯平面直角坐标，即得到像点对应的地面点空间坐标。

24)将求得足够的像点及对应地面点空间坐标代入式(1)，使式(1)中的l_x＝0、l_y＝0即为控制点参数解算模型，求得模型参数的初始值。

所述的步骤3)计算模型精确参数如下：

通过在像方和物方选取的同名直线，以及每对直线上任意选取的点坐标对，由式(1)建立误差方程和直线方程，由于任意选取的点坐标对不是同名点，通过直线方程的迭代逼近寻找同名点，每个点坐标对依据直线的方向仅列出一个方程，当直线与水平方向夹角大于45°时取x方向的误差方程，当直线与水平方向夹角小于45°取y方向的误差方程，误差方程式是依据直线的方向建立的，选择x方向4条和y方向5条以上直线作为控制线，对x式和y式分开迭代求解，使得v_x和v_y小于一定的限差或者dL₁…dL₈小于一定的限差时迭代收敛，求得纠正模型的精确参数。

所述的步骤4)纠正星载光学遥感影像如下：

41)由式(1)求解模型参数；

42)使用间接纠正法对影像进行几何纠正，其中几何纠正模型为：

$\left\{\begin{array}{c}x=L_{1}X+L_{2}Y+L_{3}Z+L_4\\ y=\frac{(L_{5}X+L_{6}Y+L_{7}Z+L_8)\·f\cos \mathrm{\ω}}{f\cos \mathrm{\ω}-Z_i+(L_{5}X+L_{6}Y+L_{7}Z+L_8)\sin \mathrm{\ω}}\end{array}\right.---\left(8\right)$

所述的点坐标对为像点和地面点。

所述的步骤42)中的间接纠正法为通过地面点坐标及模型参数，求得地面点对应的像点位置；经灰度内插获取地面点的灰度值，逐点计算完成影像纠正。

与现有技术相比，本发明具有精确度高，通过在遥感影像和地形图(栅格或矢量地图)上选取足够的同名直线，结合遥感影像的卫星星历参数，实现星载遥感影像的精确几何纠正。

具体实施方式

下面结合附表和具体实施例子对本发明进行详细说明。

实施例

一种基于直线特征的星载光学线阵影像精确纠正处理方法，不使用控制点，直接利用遥感影像与地面的直线特征，对星载光学线阵影像进行精确纠正，具体步骤如下：

1)建立模型；

2)计算模型初始参数；

3)计算模型精确参数；

4)纠正星载光学遥感影像。

以下结合一个实施例子对具体实现方法进行说明，即要对一幅高分辨率星载光学遥感影像(SPOT5影像)实现基于线特征的几何纠正，依据所使用的遥感数据SPOT5影像，经过基于线特征的几何纠正，其精度可以满足1∶5万地形图的要求。

实现过程如下：

(1)数据准备。使用SPOT5原始1A级影像，并有星历参数*.DIM文件；使用不小于1∶5万比例尺的地形图，5万比例尺地形图扫描时分辨率要大于400dpi，以保证地形图坐标获取的精度。直线是由选取的起点和终点坐标进行描述，影像上为像素坐标，地形图上为高斯坐标(X，Y，Z)。为了有效地获取高斯坐标，实际工作中开发了处理软件，以便通过鼠标位置直接获取地形图上直线的起点和终点平面坐标，高程值由人工输入。

(2)直线特征采集。分别在遥感影像和地形图上选取对应地物的直线特征，记录对应直线起点和终点坐标。

(3)计算模型初始参数。模型初始参数计算使用控制点坐标，可以选取每条直线的起点或终点坐标，结合星历表参数，直接求得像点对应的地面点空间坐标，由于精度不高，所以求解的模型参数用作模型精确参数计算的初始参数使用。表1为像点到地面点的计算精度，因此，模型参数获取不准确。

(4)计算模型精确参数。依据公式(1)，使用每组直线坐标列一个误差方程式，误差方程式的确定，可以按照直线的斜率确定x或y方向误差方程式，步骤(3)计算的模型参数作为初值，然后对式(1)进行迭代计算，迭代收敛后可以得到精确的模型参数。

表1像点到地面点计算精度

控制点号	地面X坐标误差(米)	地面Y坐标误差(米)	点位误差(米)
				1	20.72678	22.51614	30.60353
2	44.79316	16.81269	47.84448
				3	50.24334	2.144567	50.28909
4	20.34163	24.41845	31.78117
				5	35.97434	10.56796	37.49447
6	52.12521	0.623645	52.12894
				7	51.5903	6.687276	52.02191
8	54.00086	9.293314	54.7947
				9	46.99083	15.49383	49.47926
10	37.43913	11.8349	39.26517
				11	14.75742	15.24587	21.21835
12	-6.33322	13.14872	14.59447
				13	-1.84421	10.16476	10.33071
14	34.27994	11.02705	36.00986
				平均误差(米)	32.50611	12.14137	37.70401

(5)纠正星载光学遥感影像。求得精确模型参数后，将参数代入公式(8)，可以建立像点与地面点之间的严密几何关系。通常使用间接法几何纠正，在正射影像范围和像元地面分辨率确定后，可以逐点计算(X_ij，Y_ij，Z_ij)(ij为纠正影像点的行列号)对应的像点位置，对像点位置进行灰度内插，可以获得纠正的正射影像。表2为公式(8)代入精确模型参数后，计算地面点对应像点坐标与实际量测像点坐标误差。

表2基于线特征的线阵影像纠正计算结果

点号	x方向像元误差(像元)	y方向像元误差(像元)	像点位置误差(像元)
				1	-6.21979	1.038673	6.305923
2	3.584152	1.351295	3.830424
				3	8.260362	-0.61607	8.283304
4	1.520927	0.624349	1.64409
				5	-3.28601	-1.77276	3.733704
6	7.204223	-0.37207	7.213824
				7	4.71162	-0.46773	4.734779
8	8.651549	0.91448	8.699746
				9	5.142953	1.203605	5.281915
10	3.809867	-0.34838	3.825762
				11	-1.55956	-0.98447	1.844289
12	1.249474	-2.12605	2.466022
				13	2.161462	-3.02807	3.720363
14	3.055856	-0.69801	3.134562
				平均误差	2.866577	1.034735	3.047612

Claims (3)

1.一种基于直线特征的星载光学线阵影像精确纠正处理方法，其特征在于，直接利用遥感影像与地面的直线特征，对星载光学线阵影像进行精确纠正，具体步骤如下：

1）建立模型；

2）计算模型初始参数；

3）计算模型精确参数；

4）纠正星载光学遥感影像；

所述的步骤1）建立模型如下：

对于星载光学线阵影像，其直线特征上的不同点具有不同的外方位参数，基于星载光学遥感影像的几何特点，利用精确仿射变换模型建立像点与物点之间的几何关系，线性化后的计算模型和条件方程误差方程式为：

$\left\{\begin{array}{c}{V}_x={\mathrm{XdL}}_1+{\mathrm{YdL}}_2+{\mathrm{ZdL}}_3+{\mathrm{dL}}_4-(x-x_i)-l_x\\ V_y=\frac{1}{D}({\mathrm{XdL}}_5+{\mathrm{YdL}}_6+{\mathrm{ZdL}}_7+{\mathrm{dL}}_8+\frac{y}{(f-y\tan \mathrm{\ω}){\cos }^2\mathrm{\ω}}\·(\frac{Z_i\sin \mathrm{\ω}}{D}-y)\mathrm{d\ω})-(y-\frac{A}{D_i})-l_y\end{array}\right.---\left(1\right)$

其中，

l_x＝(y-y_a)×(x_b-x_a)÷(y_b-y_a)+x_a-x

l_y＝(x-x_a)×(y_b-y_a)÷(x_b-x_a)+y_a-y

$D=\frac{f-Z_i/\cos \mathrm{\ω}}{f-y\tan \mathrm{\ω}}$

(x_a,y_a)和(x_b,y_b)为确定影像上对应直线的两个端点a和b，l_x、l_b为像方直线任意点在像方x轴和y轴的投影，ω为传感器侧视角，(x,y)为影像上得到的像点坐标，(X,Y,Z)为地图上获取的高斯坐标；

所述的步骤2）计算模型初始参数如下：

21）依据像元的行坐标，从影像的星历数据中获取该行成像时刻的视方向，并利用中心行影像记录扫描时间、每行扫描间隔以及像点所在行，确定像点的姿态记录时刻，使用拉格朗日插值运算，由成像姿态数据计算所选像点成像时刻的成像位置、速度、成像姿态角：

位置与速度由下列公式计算

$\stackrel{\→}{P}\left(t\right)=\underset{j=1}{\overset{8}{\mathrm{\Σ}}}\frac{\stackrel{\→}{P}\left(t_j\right)\×\underset{\underset{i\≠j}{i=1}}{\overset{8}{\mathrm{\Π}}}(t-t_i)}{\underset{\underset{i\≠j}{i=1}}{\overset{8}{\mathrm{\Π}}}(t_j-t_i)}---\left(2\right)$

$\stackrel{\→}{V}\left(t\right)=\underset{j=1}{\overset{8}{\mathrm{\Σ}}}\frac{\stackrel{\→}{V}\left(t_j\right)\×\underset{\underset{i\≠j}{i=1}}{\overset{8}{\mathrm{\Π}}}(t-t_i)}{\underset{\underset{i\≠j}{i=1}}{\overset{8}{\mathrm{\Π}}}(t_j-t_i)}---\left(3\right)$

式中：

为卫星位置坐标，

为卫星速度坐标，t_i为相应位置与速度的宇宙时间；

成像姿态角使用线性内插得到某一观测时间t对应影像某行的姿态值[a_p(t),a_γ(t),a_y(t)]，成像姿态角内插计算公式：

$a_p\left(t\right)=a_p\left(t_i\right)+(a_p\left(t_{i+1}\right)-a_p\left(t_i\right))\×\frac{t-t_i}{t_{i+1}-t_i}---\left(4\right)$

式中：a_p(t)和a_p(t_i)分别为t和t_i时刻绕俯仰轴的转角，a_γ(t)为t时刻绕滚动轴的转角,a_y(t)为t时刻绕偏轨轴的转角，a_γ(t)和a_y(t)内插计算同（4）式；

22）将轨道坐标系统的视方向转换到地心坐标系下的视方向，再在地心坐标系下求解视方向与地球椭球的交点：

计算像点在轨道坐标系统中的视方向

${\stackrel{\→}{u}}_2=\frac{{\stackrel{\→}{u}}_2^{\′}}{\left|\right|{\stackrel{\→}{u}}_2^{\′}\left|\right|}---\left(5\right)$

式中：

${\stackrel{\→}{u}}_2^{\′}=M_p\·M_r\·M_y\·{\stackrel{\→}{u}}_1$

$M_p=\left[\begin{array}{ccc}1& 0& 0\\ 0& \cos \left(a_p\left(t\right)\right)& \sin \left(a_p\left(t\right)\right)\\ 0& -\sin \left(a_p\left(t\right)\right)& \cos \left(a_p\left(t\right)\right)\end{array}\right]$

$M_r=\left[\begin{array}{ccc}\cos \left(a_r\left(t\right)\right)& 0& -\sin \left(a_r\left(t\right)\right)\\ 0& 1& 0\\ \sin \left(a_r\left(t\right)\right)& 0& \cos \left(a_r\left(t\right)\right)\end{array}\right]$

$M_y=\left[\begin{array}{ccc}\cos \left(a_y\left(t\right)\right)& -\sin \left(a_y\left(t\right)\right)& 0\\ \sin \left(a_y\left(t\right)\right)& \cos \left(a_y\left(t\right)\right)& 0\\ 0& 0& 1\end{array}\right]$

将轨道坐标系统的视方向转换到地心坐标系下的视方向

转换模型如下：

${\stackrel{\→}{u}}_3=\left[\begin{array}{ccc}{\left(a_p\left(t\right)\right)}_p& {\left(a_{\mathrm{\γ}}\left(t\right)\right)}_p& {\left(a_y\left(t\right)\right)}_p\\ {\left(a_p\left(t\right)\right)}_{\mathrm{\γ}}& {\left(a_{\mathrm{\γ}}\left(t\right)\right)}_{\mathrm{\γ}}& {\left(a_y\left(t\right)\right)}_{\mathrm{\γ}}\\ {\left(a_p\left(t\right)\right)}_y& {\left(a_{\mathrm{\γ}}\left(t\right)\right)}_y& {\left(a_y\left(t\right)\right)}_y\end{array}\right]\·{\stackrel{\→}{u}}_2---\left(6\right)$

23）将该交点坐标转换为高斯平面直角坐标，即得到像点对应的地面点空间坐标：

在地心坐标系下求解视方向与地球椭球的交点，从卫星所处位置

得到视方向

计算在地球椭球体上方h高处视方向与椭球体上的交点(X,Y,Z)，公式如下：

$\left\{\begin{array}{c}X=X_P+\mathrm{\μ}\×{\left(u_3\right)}_p\\ Y=Y_P+\mathrm{\μ}\×{\left(u_3\right)}_{\mathrm{\γ}}\\ Z=Z_P+\mathrm{\μ}\×{\left(u_3\right)}_y\end{array}\right.---\left(7\right)$

式中：(X_p,Y_p,Z_p)为成像时刻卫星的空间位置，由式（2）求得；μ为需要求解的系数；由于椭球体上的点(X,Y,Z)满足椭球方程：

$\frac{X^2+Y^2}{A^2}+\frac{Z^2}{B^2}=1$

式中：A＝a+h，B＝b+h，a,b,h分别为椭球长半轴、短半轴和地面高程，没有地面高程模型数据时h＝0，将式（7）代入椭球方程并整理可得：

$[\frac{{\left(u_3\right)}_p^2+{\left(u_3\right)}_{\mathrm{\γ}}^2}{A^2}+\frac{{\left(u_3\right)}_y^2}{B^2}]\×{\mathrm{\μ}}^2+2\×[\frac{X_P{\left(u_3\right)}_p+Y_P{\left(u_3\right)}_{\mathrm{\γ}}}{A^2}+\frac{Z_p{\left(u_3\right)}_y}{B^2}]\×\mathrm{\μ}+[\frac{X_P^2+Y_P^2}{A^2}+\frac{Z_P^2}{B^2}]=1$

求解该一元二次方程，得到两个不同解(μ₁,μ₂)，取最小解(μ_min)为式（7）的系数，并代入式（7）求得地面点，最后将该坐标转换为高斯平面直角坐标，即得到像点对应的地面点空间坐标；

24）将求得足够的像点及对应地面点空间坐标代入模型（1），使式（1）中的l_x＝0、l_y＝0即为控制点参数解算模型，求得模型参数的初始值；

所述的步骤3）计算模型精确参数如下：

通过在像方和物方选取的同名直线，以及每对直线上任意选取的点坐标对，由式（1）建立误差方程和直线方程，由于任意选取的点坐标对不是同名点，通过直线方程的迭代逼近寻找同名点，每个点坐标对依据直线的方向仅列出一个方程，当直线与水平方向夹角大于45°时取x方向的误差方程，当直线与水平方向夹角小于45°取y方向的误差方程，误差方程式是依据直线的方向建立的，选择x方向4条和y方向5条以上直线作为控制线，对x式和y式分开迭代求解，使得v_x和v_y小于一定的限差或者dL₁…dL₈小于一定的限差时迭代收敛，求得纠正模型的精确参数；

所述的步骤4）纠正星载光学遥感影像如下：

41）由式（1）求解模型参数；

42）使用间接纠正法对影像进行几何纠正，其中几何纠正模型为：

$\left\{\begin{array}{c}x=L_{1}X+L_{2}Y+L_{3}Z+L_4\\ y=\frac{(L_{5}X+L_{6}Y+L_{7}Z+L_8)\·f\cos \mathrm{\ω}}{f\cos \mathrm{\ω}-Z_i+(L_{5}X+L_{6}Y+L_{7}Z+L_8)\sin \mathrm{\ω}}\end{array}\right.---\left(8\right)$

2.根据权利要求1所述的一种基于直线特征的星载光学线阵影像精确纠正处理方法，其特征在于，所述的点坐标对为像点和地面点。

3.根据权利要求1所述的一种基于直线特征的星载光学线阵影像精确纠正处理方法，其特征在于，所述的步骤42）中的间接纠正法为通过地面点坐标及模型参数，求得地面点对应的像点位置；经灰度内插获取地面点的灰度值，逐点计算完成影像纠正。