作业4--空间后方交会
【武汉大学-摄影测量学-单张相片解析】3.5.5单片空间后方交会
cos
0
s
in
0 0 1 0 0 0
1 0 0
X
YZ
R 1 R
R1
X
Y
Z
Xs
Ys Zs
0 R1 0
1
0 0 0
1X X s
0 0
Y Z
Ys Zs
武汉大学
摄影测量基础
偏导数－2-2
X
YZ
R 1
0 0
1
0 0 0
1 X
0 0
R
YZ
c1
X s
X)
f Z 2 (a1Z a3 X )
1
X
Z
(a1 f
f
Z
a3 )
1 Z
a1 f
a3 (x
x0 )
武汉大学
摄影测量基础
偏导数－1
x X s
1 Z
a1 f
a3(x x0 )
x Ys
1 Z
b1 f
b3(x x0 )
x Z s
1 Z
c1 f
c3(x x0 )
y X s
1 Z
c2 c3
0 0 0
a1 a1
a2 a3
bc11
a2 b2 c2
a3 b3 c3
X Y Z
0
aa23cc11
a1c2 a1c3
a1c2 a2c1 0
a3c2 a2c3
a1c3 a2c3
0
a3c1 a3c2
X
YZ
0 bb32
b3 0
b1
b2 b1
0
X
武汉大学
摄影测量基础
误差方程的建立
☺ 已知值 x0 , y0 , f , m, X, Y, Z ☺ 观测值 x, y
双向解析摄影测量
∵ ∴
X 1Y2 − X 2Y1 Y =− 2 X 2 Z1 − X 1Z 2 Z2
BxY1 Y = − 2 N2 X 2 Z1 − X 1Z 2 Z2
Bx (Y1Y2 + Z1Z 2 ) Y22 = −( Z 2 + ) N 2 X 2 Z1 − X 1Z 2 Z2
线性化方程
∂F ∂F ∂F ∂F ∂F ∆µ + ∆ν + ∆ϕ + ∆ω + ∆κ = 0 F =F + ∂µ ∂ν ∂ϕ ∂ω ∂κ
0
∂F = Bx ( X 2 Z1 − X1Z2 ) ∂µ
∂F ≈ B xY1 X 2 ∂φ
∂F = Bx ( X 1Y2 − X 2Y1 ) ∂ν
∂F ≈ − Bx X 2 Z1 ∂κ
X2 ∂ Y 2 0 Z 2 = cos ϕ cos ω ∂κ sin ω
− cos ϕ cos ω 0 − sin ϕ cos ω
− sin ω X 2 − Y2 Y ≈ X sin ϕ cos ω 2 2 0 Z2 0
Bz Z1 = 0 Z2
X 2 x2 Y2 = R y 2 Z − f 2
连续法解析相对定向原理 连续法解析相对定向原理
1 F = Bx X 1 X2
µ
Y1 Y2
ν
Z1 = 0 Z2
∂F ∂F ∂F ∂F ∂F ∆µ + ∆ν + ∆ϕ + ∆ω + ∆κ = 0 F=F + ∂µ ∂ν ∂ϕ ∂ω ∂κ
投影误差公式
问题 提出
5
6
RC30 ADS40（Leica/Helava） DMC
7
8
空中摄影就是 从空中一定高度上 摄取地面物体影像 的过程。
9
航摄像片的质量主要包括：
①构像质量: ②几何质量: ③表观质量:
影像的分解力，清晰度。 体现在影像的量测性能方面。 色调正常，反差适中。
航空摄影的飞行质量要求：
问题 提出
13
内 容 安 排
14
投影[Projection]
一个空间点按一定方式在一个平面上的构 像，叫做该空间点的投影。
物点、像点、投射线、像面（承影面）
平行投影[Parallel Projection ]
投射线互相平行的投影，叫做平行投影。
正射投影（垂直投影）
局部范围内，地形图可视为地面的垂直投影。
29
向量表达式
z y x S
SA λ Sa
-f
y
xA yA zA x y f
x
a
o
y
A
zA
xA
x f y f
xA zA yA zA
（1）
30
坐标变换
xA X T y A R Y z Z A
S
P yc
记为： a
f h f c c P0
倾斜误差 m0 a m a m0 yc0 hc xc,xc0 38 M
倾斜误差公式
a rc r0
rc sin sin f rc sin sin
2
2
i
yc , yc
0
r0 rc
y
c
摄影测量习题整理
习题——认识摄影测量1、下列关于摄影测量的发展阶段的说法中，不正确的是（B)（A）数字摄影测量处理的原始资料是数字影像或数字化影像（B）数字摄影测量和解析摄影测量使用的投影方式分别是数字投影和物理投影（C）模拟摄影测量和解析摄影测量处理的原始数据都是像片（均为数字投影）（D）数字摄影测量是以计算机视觉代替人的立体观测与识别，完成影像几何与物理信息的自动提取2、下列有关数字地面模型的相关概念的描述中，不正确的是（D)（A）数字地面模型DTM是地理信息系统的基础数据，可用于土地利用现状的分析、合理规划及洪水子险情预报等（B）数字地面模型DTM是地形表面形态等多种信息的一个数字表示（C）数字高程模型DEM只考虑DTM的地形分量，在计算机高级语言中，它就是一个三维数组或数学上的一个三维矩阵（D）与TIN相比，DEM能较好地顾及地貌特征点、线、面，表示复杂地形表面，其缺点是数据量较大，数据结构较复杂（与DEM相比，TIN能较好地顾及地貌特征点线、面，表示复杂地形表面，其缺点是数据量较大，数据结构较复杂。
）3、摄影测量的技术手段不包括（C）A模拟法B.解析法C实地法D数字法4、光学纠正仪是（A）时代的产品，其投影方式属于机械投影。
A模拟摄影测量B解析摄影测量C数字摄影测量D数字投影5、摄影机框架四边中点设有的框标记号，叫做光学框标。
（×）6、现有地图的数字化方法有手扶跟踪数字化和扫描数字化7、随着摄影测量技术的发展，摄影测量经历了（模拟摄影测量阶段）、（解析摄影测量阶段）与（数字摄影测量阶段 )三个发展阶段。
摄影测量的技术手段有（模拟法）、（解析法）、（数字法）。
8、摄影测量的基本问题，就是将中心投影的像片转换为正射投影的地形图9、数字摄影测量数字摄影测量是基于数字影像和摄影测量的基本原理，应用计算机技术、数字影像处理、影像匹配、模式识别等多学科的理论与方法，提取所摄对像以数字方式表达的几何与物理信息的摄影测量学的分支学科。
摄影测量学习题答案
摄影测量学习题一、名词解释:1、摄影测量学：摄影测量学是对研究的对象进行摄影,根据所获得的构像信息，从几何方面和物理方面加以分析研究，从而对所摄对象的本质提供各种资料的一门学科。
2、光圈号数×3、景深：光于点A，但在远景B和近景C之间这一段间隔内所有景物，在像面上仍可认为获得了清晰的构像。
此时，远景与近景之间的纵深距离称为景深4、超焦点距离：当物镜向无限远物体对光时,不仅远处的物体构象清晰，而且在离开物镜不小于某一距离H的所有物体,其构象都很清晰，这个距离H就称为超焦点距离或称为无限远起点。
5、视场×6、视场角×7、像场×8、像场角×9、反差系数×10、感光度×11、航向重叠（p):沿飞行方向上相邻像片所摄地面的重叠区.12、旁向重叠（q）：两相邻航带摄区之间的重叠。
13、摄影基线（B)：相邻像片摄影站(投影中心)之间的空间连线.14、透视旋转定律×15、内方位元素：确定物镜后节点和像片面相对位置的数据。
包括像主点在像片框标坐标系中的x0 坐标、y0和像片主距f。
16、外方位元素:确定摄影摄影机或像片的空间位置和姿态的参数，亦即摄影光束空间位置和姿态的数据。
17、航向倾角×18、旁向倾角×19、像片旋角×20、倾斜误差：因像片倾斜引起的像点位移.21、投影差：因地形起伏引起的像点位移。
22、摄影比例尺：构像比例尺：航摄像片上某一线段构像的长度与地面上相应线段水平距离之比.23、像片控制点:测定了地面坐标的像点24、像片判读×25、左右视差26、上下视差27、核点：基线延长线与左、右像片的交点k1、k2称为核点。
28、核线：核面与像片的交线称为核线29、核面：通过摄影基线S1S2与任一地面点A所作的平面W A，称为点A的核面。
30、投影基线：将摄影B缩小到若干分之一作为投影基线b.31、像片基线32、解析空中三角测量：是将建立的投影光束、单元模型或航带模型以至区域模型的数学模型，根据少量地面控制点，按最小二乘法原理进行平差计算，解求出各加密点的地面坐标。
摄影测量整理
第一章和第二章摄影测量学：是对研究的物体进行摄影，量测和解译所获得的影像，获取被摄物体的几何信息和物理信息的一门科学和技术。
分类：按用途分为地形摄影测量、非地形摄影测量按处理手段分为模拟摄影测量、解析摄影测量、数字摄影测量按摄影机所处的位置分为地面摄影测量、航空摄影测量、航天摄影测量发展：模拟摄影测量、解析摄影测量、数字摄影测量三个阶段。
光圈号越小，光圈就越大，能通过的光线就越多，聚焦在感光物上的亮度就越高（越亮），反之，使用的光圈号数越大，影像亮度就越小。
航空摄影的基本要求：一、摄影比例尺：摄影比例尺是指航摄设计中的像片比例尺。
像片比例尺是由摄影机的主距和摄影的高度来计算的。
摄影比例尺的确定取决于成图比例尺、摄影测量成图方法和成图精度，另外考虑经济性和摄影资料的可使用性。
摄影比例尺可以分为大、中、小三种比例尺，一般都应选择小的。
二、像片重叠度：航向重叠：沿航线飞行方向两相邻像片上的影像重叠。
（60%～65%）旁向重叠：两相邻行带像片之间的影像重叠。
（15%～30%）三度重叠：一般情况下在航线方向三张相邻像片必须要有公共重叠影像，称为三度重叠。
三、航带弯曲度：不超过3%四、像片旋偏角：相邻两像片的主点连线与像幅沿航带飞行方向的两框标连线之间的夹角称为像片的旋偏角。
一般要求像片旋偏角小于6°，个别最大不应大于8°，而且不能有连续三张超过6°的情况。
航摄像片的误差：1、底片变形2、航设计物镜畸变差3、大气折光差4、地球曲率影响三原色：红、绿、蓝（三种等量相加=白）互补色：黄和蓝、红和青、绿和品红色的形成：1、加色法：将三原色按照不同的比例进行光学混合，得到其他颜色2、减色法：从白色中减去主色而获得补色第三章单张航摄像片解析中心投影：投射线会聚于一点的投影。
投射线的会聚点S称为投影中心。
中心投影的构像过程也称透视变换。
中心投影的三种方式：航摄像片是地面的中心投影，而地形图则是地面的正摄投影。
后方交会MATLAB程序实习报告
《摄 影 测 量 学》单像空间后方交会实习报告班 级： XXXX姓 名： X X X学 号：XXXXXXXXXXXXX指导教师： X X X一、实习目的1、掌握空间后方交会的定义和实现算法；2、了解摄影测量平差的基本过程；3、熟练MATLAB 等程序编写。
二、实习原理利用至少三个已知地面控制点的坐标),,(A A A Z Y X A 、),,(B B B Z Y X B 、),,(C C C Z Y X C ，与其影像上对应的三个像点的影像坐标),(a a y x a 、),(b b y x b 、),(c c y x c ，根据共线方程，反求该像片的外方位元素κωϕ、、、、、S S S Z Y X 。
共线条件方程式：将共线式线性化并取一次小值项得：三、解算过程①获取已知数据。
包括影像比例尺1/m，平均摄影距离（航空摄影的航高）H，内方位元素x0、y0、f，控制点的空间坐标X、Y、Z。
②量测控制点的像点坐标并进行必要的影像坐标系统误差改正，得到像点坐标。
③确定未知数的初始值。
单像空间后方交会必须给出待定参数的初始值，在竖直航空摄影且地面控制点大体对称分布的情况下，Xs0和Ys0为均值，Zs0为航高，φ、ω、κ的初值都设为0。
或者κ的初值可在航迹图上找出或根据控制点坐标通过坐标正反变换求出。
④计算旋转矩阵R。
利用角元素近似值计算方向余弦值，组成R阵。
⑤逐点计算像点坐标的近似值。
利用未知数的近似值按共线条件式计算控制点像点坐标的近似值（x），(y)。
⑥逐点计算误差方程式的系数和常数项，组成误差方程式。
⑦计算法方程的系数矩阵ATA与常数项ATL，组成法方程式。
⑧解求外方位元素。
根据法方程，解求外方位元素改正数，并与相应的近似值求和，得到外方位元素新的近似值。
⑨检查计算是否收敛。
将所求得的外方位元素的改正数与规定的限差比较，通常对φ，ω，κ的改正数△φ，△ω，△κ给予限差，通常为0.000001弧度，当3个改正数均小于0.000001弧度时，迭代结束。
(完整版)摄影测量作业部分答案
1.什么是摄影测量学，摄影测量发展的三个阶段摄影测量学是对研究的物体进行摄影，测量和解译所获得的影像，获取被摄物体的几何信息和物理信息的一门科学。
三个阶段：模拟、解析、数字摄影测量。
2.根据成图的需要，规定了摄影比例尺后，如何选择航空摄影机与航摄高度？答：采用特宽角航摄机，航高值就小，采用常角或窄角航摄机，航高值就大。
（航高的大小将决定飞机实际升限和最小安全高度的限制，另外，测图的高程测定精度与航高有关（高程中误差与航高成正比））大比例尺单像测图，应选用常角或窄角航摄机，小比例尺立体测图应选用特宽角航摄机。
3.什么是航摄像片的内外方位元素，各有何作用？答：内方位元素包括三个参数，即摄影中心S到像片的垂距(主距)f及像主点o在像框标坐标系中的坐标00,y x，用其来恢复摄影光束。
确定摄影光束在摄影瞬间的空间位置和姿态的参数，称为外方位元素，一张的外方位元素包括六个参数，其中有三个是直线元素，用于描述摄影中心的空问坐标值；另外三个是角元素，用于表达像片面的空间姿态。
4.摄影测量中常用的坐标系有哪些，各有何作用？答：摄影测量中常用的坐标系有两大类。
一类是用于描述像点的位置，称为像方空间坐标系；另——类是用于描述地面点的位置．称为物方空间坐标系。
(1)像方空间坐标系①像平面坐标系像平面坐标系用以表示像点在像平面上的位置，通常采用右手坐标系，xy轴的选择按需要而定．在解析和数字摄影测量中，常根据框标来确定像平面坐标系，称为像框标坐标系。
②像空间坐标系为了便于进行空间坐标的变换，需要建立起描述像点在像空间位置的坐标系，即像空间坐标系。
以摄影中心S为坐标原点，xy,轴与像平面坐标系的xy轴平行，z轴与主光轴重合，形成像空间右手直角坐标系s-xyz③像空间辅助坐标系像点的像空间坐标可直接以像平面坐标求得，但这种坐标的待点是每张像片的像空间坐标系不统一，这给计算带来困难。
为此，需要建立一种相对统一的坐标系．称为像空间辅助坐标系，用s-XYZ表示。
摄影测量重点总结
摄影测量重点总结1、摄影测量中常用的坐标系有像平面直角坐标系、像空间直角坐标系、像空间辅助坐标系、地面摄影测量坐标系、地面测量坐标系。
2、解求单张像片的外方位元素最少需要3个平高地面控制点。
3、gps辅助空中三角测量的促进作用就是大量增加甚至全然免去地面控制点，缩为图周期，提升生产效率，降低生产成本。
4、两个空间直角坐标系间的坐标变换最少需要2个平高和1个高程地面控制点。
5、摄影测量的发展经历了模拟摄影测量、解析摄影测量和数字摄影测量三个阶段。
6、恢复立体像对左右像片的相互位置关系依据的是共面条件方程。
7、法方程消元的通式为=8、表示航摄像片的外方位角元素可以采用以y轴为主轴的?-ω-κ、以x轴为主轴的ω'-?'-κ'以z轴为主轴的a-a?k三种转角系统。
9、航摄像片是所覆盖地物的中心投影。
10、摄影测量加密按数学模型可以分成航带法、单一制模型法和光束法三种方法。
摄影测量加密按黄赤范围可以分成单模型法、航带法和区域网法三种方法。
11、从航摄像片上量测的像点坐标可能带有摄影材料变形、摄影机物镜畸变、大气折光误差和地球曲率误差四种系统误差。
12、要将地物点在摄影测量坐标系中的模型坐标转换到地面摄影测量坐标系，最少需要2个平高和1个高程地面控制点。
13、带状法方程系数矩阵的带宽是指法方程系数矩阵中主对角线元素起沿某一行到最远处的非零元素间所包含的未知数个数。
14、人眼观察两幅影像能产生立体视觉的基本条件是在不同摄站获取的具有一定重叠的两幅影像、观察时每只眼睛只能看一张像片、两幅影像的摄影比例尺尽量一致和两幅影像上相同地物的连线与眼基线尽量平行。
15、中心投影的共线条件方程抒发了摄影中心、像是点和对应地物点三点坐落于同一直线的几何关系，利用其解求单张像片6个外方位元素的方法称作单片空间后方交会，最少须要3个上恩地面控制点。
16、摄影测量中，为了恢复立体像对两张像片之间的相互位置关系，可以根据左右像片上的同名像点位于同一核面的几何条件，采用相对定向方法来实现，最少需要量测5对同名像点。
摄影测量考试试题及详细答案
摄影测量考试试题及详细答案1摄影测量学 2航向重叠3单像空间后⽅交会 4相对⾏⾼5像⽚纠正 6解析空中三⾓测量7透视平⾯旋转定律 8外⽅位元素9核⾯ 10绝对定向元素⼀、填空1摄影测量的基本问题，就是将_________转换为__________。
2物体的⾊是随着__________的光谱成分和物体对光谱成分固有不变的________、__________、和__________的能⼒⽽定的。
3⼈眼产⽣天然⽴体视觉的原因是由于_________的存在。
4相对定向完成的标志是__________。
5光束法区域⽹平差时，若像⽚按垂直于航带⽅向编号，则改化法⽅程系数阵带宽为_______，若按平⾏于航带⽅向编号，则带宽为_________。
三、简答题1两种常⽤的相对定向元素系统的特点及相对定向元素。
2倾斜位移的特性。
3单⾏带法相对定向后，为何要进⾏⽐例尺归化？为何进⾏？4独⽴模型法区域⽹平差基本思想。
5何谓正形变换？有何特点？四、论述题1空间后⽅交会的结算步骤。
2有三条航线，每条航线六像⽚组成⼀个区域，采⽤光束法区域⽹平差。
（1）写出整体平差的误差⽅程式的⼀般式。
（2）将像⽚进⾏合理编号，并计算带宽，存容量。
（3）请画出改化法⽅程系数阵结构简图。
A卷答案：⼀、1是对研究的对象进⾏摄影，根据所获得的构想信息，从⼏何⽅⾯和物理⽅⾯加以分析研究，从⽽对所摄影的对象本质提供各种资料的⼀门学科。
2供测图⽤的航测相⽚沿飞⾏⽅向上相邻像⽚的重叠。
3知道像⽚的⽅位元素，以及三个地⾯点坐标和量测出的相应像点的坐标，就可以根据共线⽅程求出六个外⽅位元素的⽅法。
4摄影瞬间航摄飞机相对于某⼀索取基准⾯的⾼度。
5将中⼼投影转换成正射投影时，经过投影变换来消除相⽚倾斜所引起的像点位移，使它相当于⽔平相⽚的构象，并符合所规定的⽐例尺的变换过程。
6是将建⽴的投影光束，单元模型或航带模型以及区域模型的数字模型，根据少数地⾯控制点，按最⼩⼆乘法原理进⾏平差计算，并求加密点地⾯坐标的⽅法。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实用 文档 作业报告 空间后方交会
专 业： 测绘工程 班 级： 2008级(1)班 姓 名： 闻亚 指导教师： 强
2010 年 4 月 16 日 实用
文档 1 作业任务 ------------------------------------------------------------------------------------ 3
2 作业思想 --------------------------------------------------------------------------------------- 3
3 作业条件及数据-------------------------------------------------------------------- 3
4 作业过程 --------------------------------------------------------------------------- 3
5 源程序 ----------------------------------------------------------------------------- 4
6 计算结果 --------------------------------------------------------------------------- 17
7心得体会与建议----------------------------------------------------------------------------- 17 实用
文档 1 作业任务 计算近似垂直摄影情况下后方交会解。即利用摄影测量空间后方交会的方法，获取相片的6个外方位元素。限差为0.1。
2作业思想 利用摄影测量空间后方交会的方法求解。该方法的基本思想是利用至少三个一直地面控制点的坐标A（XA，YA，ZA）、B（XB，YB，ZB）C（XC，YC，ZC），与其影像上对应的三个像点的影像坐标a（xa，ya）、b（xb，yb）、c（xc，yc），根据共线方程，反求该相片的外方位元素XS、YS、ZS、φ、ω、κ。
3作业条件及数据 已知摄影机主距f=153.24mm，四对点的像点坐标与相应的地面坐标列入下表： 表1 点号 像点坐标 地面坐标 x（mm） y(mm) X(m) Y(m) Z(m) 1 -86.15 -68.99 36589.41 25273.32 2195.17 2 -53.40 -82.21 37631.08 31324.51. 728.69 3 14.78 76.63 39100.97 24934.98 2386.50 4 10.46 64.43 40426.54 30319.81 757.31
4作业过程 4．1 获取已知数据 相片比例尺1/m=1：10000，方位元素f=153.24mm，x0，y0；获取控制点的地面测量坐标Xt、Yt、Zt。 4．2 量测控制点的像点坐标： 本次作业中为已知。见表1。 4．3 确定未知数的初始值： 在近似垂直摄影情况下，胶原素的初始值为0，即φ0 = ω0 = κ0=0；线元素中，ZS0=H=mf=1532.4m，XS0、YS0的取值可用四个控制点坐标的平均值，即：
XS0= 14 4ii1X=38437.00
YS0= 14 4ii1Y=89106.62 4．4 计算旋转矩阵R： 利用胶原素的近似值计算方向余弦值，组成R阵。 4．5 逐点计算像点坐标的近似值： 利用未知数的近似值按共线方程式计算控制点像点坐标的近似值（x）（y）。 4．6 组成误差方程： 实用 文档 逐点计算误差方程式的系数和常数项。 4．7 组成法方程式： 计算法方程的系数矩阵ATA与常数项ATL。 4．8 求解外方位元素： 根据法方程，由式X=（AtA）-1 ATL解求外方位元素改正数，并与相应的近似值求和，得到外方位元素新的近似值。 4．9 求解外方位元素： 将求得的外方位元素的改正数与规定的限差（0.1）比较，小于限差则计算终止，否则用新的近似值重复第4.4至4.8步骤的计算，知道满足要求为止。
5 源程序 #include #include #include
const double PRECISION=1e-5; typedef double DOUBLE[5]; int InputData(int &Num, DOUBLE *&Data,double &m,double &f); int Resection(const int &Num,const DOUBLE *&Data,const double &m,const double &f); int InverseMatrix(double *matrix,const int &row);
int main(int argc, char* argv[]) { DOUBLE *Data=NULL; int Num; double f(0),m(0); if(InputData(Num,Data,m,f)) { if (Data!=NULL) { delete []Data; } return 1; } if(Resection(Num,Data,m,f)) { if (Data!=NULL) { delete []Data; } return 1; } 实用 文档 if (Data!=NULL) { delete []Data; } printf("解算完毕...\n");
do{ printf("计算结果保存于\"结果.txt\"文件中\n" "请选择操作（输入P打开结果数据，R打开原始数据，其它退出程序）："); fflush(stdin); //刷新输入流 char order=getchar(); if ('P'==order || 'p'==order) { system("结果.txt"); } else if ('R'==order || 'r'==order) { system("data.txt"); } else break; system("cls"); }while(1); system("PAUSE");
return 0; }
/********************************************** *函数名：InputData *函数介绍：从文件（data.txt）中读取数据, *文件格式如下： *点数 m(未知写作0) * 方位元素(f x0 y0) *编号 x y X Y Z *实例： 4 0 153.24 0 0 1 -86.15 -68.99 36589.41 25273.32 2195.17 2 -53.40 82.21 37631.08 31324.51 728.69 3 -14.78 -76.63 39100.97 24934.98 2386.50 4 10.46 64.43 40426.54 30319.81 757.31
*参数：(in/out)Num(点数)， 实用 文档 *(in/out)Data(存放数据),m,f,x0,y0 *返回值：int ，0成功，1文件打开失败，2控制点个 *数不足，3文件格式错误 **********************************************/ int InputData(int &Num, DOUBLE *&Data,double &m,double &f) { double x0,y0; FILE *fp_input; if (!(fp_input=fopen("data.txt","r"))) { return 1; }
fscanf(fp_input,"%d%lf",&Num,&m); if (Num<4) { return 2; }
fscanf(fp_input,"%lf%lf%lf",&f,&x0,&y0); f/=1000; if (m<0 || f<0) { return 3; }
Data=new DOUBLE[Num]; double *temp= new double[Num-1]; double scale=0;
int i; for (i=0;i{ //读取数据，忽略编号 if(fscanf(fp_input,"%*d%lf%lf%lf%lf%lf", &Data[i][0],&Data[i][1],&Data[i][2], &Data[i][3],&Data[i][4])!=5) { return 3; } //单位换算成m Data[i][0]/=1000.0; Data[i][1]/=1000.0; }