实验 ERDAS视窗的基本操作 遥感图像的几何校正 遥感图像合成

实验 ERDAS视窗的基本操作 遥感图像的几何校正 遥感图像合成
实验 ERDAS视窗的基本操作 遥感图像的几何校正 遥感图像合成

实验 ERDAS视窗的基本操作遥感图像的几何校正遥感图像合成

一、实验目的:初步了解目前主流的感遥图象处理软件ERDAS的主要功能模块,在此基础上,掌握视窗操作模块的功能和操作技能,为遥感图像的几何校正等后续实习奠定基础。

二、实验内容:视窗功能介绍;文件菜单操作;实用菜单操作;显示菜单操作;矢量和删格菜单操作等。

视窗操作是ERDAS软件操作的基础, ERDAS所有模块都涉及到视窗操作。要求掌握视窗的基本功能,熟练掌握图像显示操作和矢量菜单操作,从而为深入理解和学习ERDAS软件打好基础。

1、视窗功能简介

二维视窗是显示删格图像、矢量图形、注记文件、AOI等数据层的主要窗口。通过实际操作,掌握视窗菜单的主要功能、视窗工具功能。

重点掌握ERDAS图表面板菜单条;ERDAS图表面板工具条;掌握视窗菜单功能和视窗工具功能等基本操作。

2、图像显示操作(Display an Image)

第一步:启动程序(Start Program)

视窗菜单条:File→open→ RasterLayer→Select Layer To Add对话框。

第二步:确定文件(Determine File)

在Select Layer To Add对话框中有File和Raster Option两个选择项,其中File就是用于确定图像文件的,具体内容和操作实例如表。

表1-1 图像文件确定参数

参数项含义实例

Look in 确定文件目录examples

File name 确定文件名xs_truecolor_sub.img

File of type 确定文件类型IMAGINE Image(*.img)

------

Recent 选择近期操作过的文

Go to 改变文件路径-------

第三步:设置参数(Raster option)

图1-1 参数设置

第四步:打开图像(Open Raster Layer)

4、实用菜单操作

了解光标查询功能;量测功能;数据叠加功能;文件信息操作;三维图像操作等。

5、显示菜单操作

掌握文件显示顺序(图1-2);显示比例;显示变换操作等。

6、矢量菜单操作

矢量菜单操作功能是ERDAS软件将遥感与地理信息系统相结合的一个体现。

图1-2 图层显示顺序

三、实验感想:初步了解了感遥图象处理软件ERDAS的主要功能模块,并在此基础上,基本掌握视窗操作模块的功能和操作技能,为遥感图像的几何校正等后续实习奠定基础。

实验二、遥感图像的几何校正

一、实验目的:通过实习操作,掌握遥感图像几何校正的基本方法和步骤,深刻理解图像几何校正的意义。

二、实验内容:ERDAS软件中图像预处理模块下的图像几何校正。

几何校正就是将图像数据投影到平面上,使其符合地图投影系统的过程。而将地图投影系统赋予图像数据的过程,称为地理参考(Geo-referencing)。由于所有地图投影系统都遵循一定的地图坐标系统,因此几何校正的过程包含了地理参考过程。

1、图像几何校正的途径

ERDAS图标面板工具条:点击DataPrep图标,→Image Geometric Correction →打开Set Geo-Correction Input File对话框(图2-1)。

ERDAS图标面板菜单条:Main→Data Preparation→Image Geometric Correction→打开Set Geo-Correction Input File对话框(图2-1)。

图2-1 Set Geo-Correction Input File对话框在Set Geo-Correction Input File对话框(图1)中,需要确定校正图像,有两种选择情况:

其一:首先确定来自视窗(FromViewer),然后选择显示图像视窗。

其二:首先确定来自文件(From Image File),然后选择输入图像。

2、图像几何校正的计算模型(Geometric Correction Model)

ERDAS提供的图像几何校正模型有7种,具体功能如下:

表2-1 几何校正计算模型与功能

模型功能

Affine 图像仿射变换(不做投影变换)

Polynomial 多项式变换(同时作投影变换)

Reproject 投影变换(转换调用多项式变换)

Rubber Sheeting 非线性变换、非均匀变换

Camera 航空影像正射校正

Landsat Lantsat卫星图像正射校正

Spot Spot卫星图像正射校正

3、图像校正的具体过程

第一步:显示图像文件(Display Image Files)

首先,在ERDAS图标面板中点击Viewer图表两次,打开两个视窗(Viewer1/Viewer2),并将两个视窗平铺放置,操作过程如下:

ERDAS图表面板菜单条:Session→Title Viewers

然后,在Viewer1中打开需要校正的Lantsat图像:tmAtlanta,img

在Viewer2中打开作为地理参考的校正过的SPOT图像:panAtlanta,img 第二步:启动几何校正模块(Geometric Correction Tool)

Viewer1菜单条:Raster→ Geometric Correction

→打开Set Geometric Model对话框(2)

→选择多项式几何校正模型:Polynomial→OK

→同时打开Geo Correction Tools对话框(3)和Polynomial Model Properties对话框(4)。

在Polynomial Model Properties对话框中,定义多项式模型参数以及投影参数:

→定义多项式次方(Polynomial Order):2

→定义投影参数:(PROJECTION):略

→Apply→Close

→打开GCP Tool Referense Setup 对话框(5)

图2-2 Set Geometric Model对话框

图2-3 Geo Correction Tools对话框

图2-4 Polynomial Properties对话框

图2-5 GCP Tool Referense Setup 对话框

第三步:启动控制点工具(Start GCP Tools)

图2-6 Viewer Selection Instructions

首先,在GCP Tool Referense Setup对话框(图5)中选择采点模式:→选择视窗采点模式:Existing Viewer→OK

→打开Viewer Selection Instructions指示器(图2-6)

→在显示作为地理参考图像panAtlanta,img的Viewer2中点击左键

→打开reference Map Information 提示框(图2-7);→OK

→此时,整个屏幕将自动变化为如图7所示的状态,表明控制点工具被启动,进入控制点采点状态。

图2-7 reference Map Information 提示框

图2-8 控制点采点

第四步:采集地面控制点(Ground Control Point)

GCP的具体采集过程:

在图像几何校正过程中,采集控制点是一项非常重要和繁重的工作,具体过程如下:

1、在GCP工具对话框中,点击Select GCP图表,进入GCP选择状态;

2、在GCP数据表中,将输入GCP的颜色设置为比较明显的黄色。

3、在Viewer1中移动关联方框位置,寻找明显的地物特征点,作为输入

GCP。

4、在GCP工具对话框中,点击Create GCP图标,并在Viewer3中点击左

键定点,GCP数据表将记录一个输入GCP,包括其编号、标识码、X坐标和Y坐标。

5、在GCP对话框中,点击Select GCP图标,重新进入GCP选择状态。

6、在GCP数据表中,将参考GCP的颜色设置为比较明显的红色,

7、在Viewer2中,移动关联方框位置,寻找对应的地物特征点,作为参

考GCP。

8、在GCP工具对话框中,点击Create GCP图标,并在Viewer4中点击左

肩顶巅,系统将自动将参考点的坐标(X、Y)显示在GCP数据表中。

9、在GCP对话框中,点击SelectGCP图标,重新进入GCP选择状态,并将光标移回到Viewer1中,准备采集另一个输入控制点。

10、不断重复1-9,采集若干控制点GCP,直到满足所选定的几何模型为止,尔后,没采集一个InputGCP,系统就自动产生一个Ref. GCP,通过移动Ref. GCP可以优化校正模型。

第五步:采集地面检查点(Ground Check Point)

以上采集的 GCP的类型均为控制点,用于控制计算,建立转换模型及多项式方程,。下面所要采集的GCP类型是检查点。

第六步:计算转换模型(Compute Transformation)

在控制点采集过程中,一般是设置为自动转换计算模型。所以随着控制点采集过程的完成,转换模型就自动计算生成。

在Geo-Correction Tools对话框中,点击Display Model Properties 图表,可以查阅模型。

第七步:图像重采样(Resample the Image)

重采样过程就是依据未校正图像的像元值,计算生成一幅校正图像的过程。原图像中所有删格数据层都要进行重采样。

ERDAS IMAGE 提供了三种最常用的重采样方法。

图像重采样的过程:

首先,在Geo-Correction Tools对话框中选择Image Resample 图标。

然后,在Image Resample对话框中,定义重采样参数;

→输出图像文件明(OutputFile):rectify.img

→选择重采样方法(Resample Method):Nearest Neighbor

→定义输出图像范围:

→定义输出像元的大小:

→设置输出统计中忽略零值:

→定义重新计算输出缺省值:

第八步:保存几何校正模式(Save rectification Model)

在Geo-Correction Tools对话框中点击Exit按钮,推出几何校正过程,按照系统提示,选择保存图像几何校正模式,并定义模式文件,以便下一次直接利用。

第九步:检验校正结果(Verify rectification Result)

基本方法:同时在两个视窗中打开两幅图像,一幅是矫正以后的图像,一幅是当时的参考图像,通过视窗地理连接功能,及查询光标功能进行目视定性检验。

三、实验感想:通过本次实验操作操作,初步掌握了遥感图像几何校正的基本方法和步骤,深刻理解了图像几何校正的意义。

实验三、遥感图像合成

一、实验目的:理解计算机图像合成的基本原理以及过程,达到能熟练地对遥感图像进行合成的目的。

二、实验基本要求:ERDAS遥感图像合成。

三、设备及器材配备:计算机、交换器。计算机配置为:内存256M,128M 独立显卡,CPU为奔腾4处理器。

四、实验项目与内容提要:ERDAS图标面板菜单条:点击Interpreter图标→Utilities→Layers Stack→Layer Selection and Stacking 对话框。

在Layer Selection and Stacking 对话框中,依次选择并加载(add)单

波段图像:→输出单波段文件(imput file:*.img):tm1.img→add

→输出单波段文件(imput file:*.img):tm2.img→add

→输出单波段文件(imput file:*.img):tm3.img→add

……

→输出多波段文件(output file:*.img):dff.img

→输出类型(output data type):unsigned 8 bit

→波段组合选择(output option):union

→输出统计忽略零值:ignore zero in stats

Ok(关闭 layer selection and stacking)对话框,执行波段组合。

五、实验感想:通过本次试验,主要掌握了计算机图像合成的基本原理、过程,以及遥感图像合成的步骤过程,初步了解了遥感图像合成的方法。

遥感图像的几何校正(配准)

遥感图像的几何校正(配准) 1.实验目的与任务: (1)了解几何校正的原理; (2)学习使用ENVI软件进行几何校正; 2.实验设备与数据: 设备:遥感图像处理系统ENVI 数据:TM数据 3 几何校正的过程: 注意:几何校正一种是影像对影像,一种是影像对地图,下面介绍的是影像对影像的配 准或几何校正。 1.打开参考影像(base)和待校正影像:分别打开,即在display#1,display#2中打开;2.在主菜单上选择map->Registration->select GCPs:image to image 3.出现窗口Image to Image Registration,分别在两边选中DISPLAY 1(左),和DISPLAY 2(右)。BASE图像指参考图像而warp则指待校正影像。选择OK! 4.现在就可以加点了:将两边的影像十字线焦点对准到自己认为是同一地物的地方, 就可以选择ADD POINT添加点了。(PS:看不清出别忘记放大)如果要放弃该点选择 右下脚的delete last point,或者点show point弹出image to image gcp list窗口,从中选择 你要删除的点,也可以进行其他很多操作,自己慢慢研究,呵呵。选好4个点后就可以 预测:把十字叉放在参考影像某个地物,点选predict则待校正影像就会自动跳转到与参 考影像相对应的位置,而后再进行适当的调整并选点。 5.选点结束后,首先把点保存了:ground control points->file->save gcp as ASCII.. 当然你没有选完点也可以保存,下次就直接启用就可以:ground control points->file->restore gcps from ASCII... 6.接下来就是进行校正了:在ground control points.对话框中选择: options->warp file(as image to map) 在出现的imput warp image中选中你要校正的影像,点ok进入registration parameters 对话框: 首先点change proj按钮,选择坐标系 然后更改象素的大小,如果本身就是你所需要大小则不用改了 最后选择重采样方法(resampling),一般都是选择双线性的(bilinear),最后的最后选择保存路径就OK了

ENVI中的几何校正

几何校正 1.遥感图像产生几何畸变的原因 地物目标发出的电磁波被卫星上所载传感器接收,这些电磁波上记录和传达了地物目标的信息,这是遥感图像成像的过程也是它的内在规律。在这个过程中图像的几何畸变也随即产生了,其中原因很多,主要表现在以下几个方面: 1. 1卫星位置和运动状态变化的影响 卫星围绕地球按椭圆轨道运动,引起卫星航高和飞行速度的变化,导致图像对应产生偏离与在卫星前进方向上的位置错动。另外,运动过程中卫星的偏航、翻滚和俯仰变化也能引起图像的畸变。 以上误差总的来说,都是因为传感器相对于地物的位置、姿态和运动速度变化产生的,属于外部误差。此外,由于传感器本身原因产生的误差,即内部误差,这类误差一般很小,通常人们不作考虑。 1. 2地球自转的影响 大多数卫星都是在轨道运行的降段接收图像,即当地球自西向东自转时,卫星自北向南运动。这种相对运动的结果会使卫星的星下位置产生偏离,从而使所成图像产生畸变。 1. 3地球表面曲率的影响 地球表面是不规则的曲面,这使卫星影像成像时像点发生移动,像元对应于地面的宽度不等。特别是当传感器扫描角度较大时,影响更加突出。 1. 4地形起伏的影响 当地形存在起伏时,使原来要反映的理想的地面点被垂直在其上的实际某高点所代替,引起图像上像点也产生相应的偏离。 1. 5大气折射的影响 由于大气圈的密度是不均匀分布的,从下向上越来越小,使得整个大气圈的折射率不断变化,当地物发出的电磁波穿越大气圈时,经折射后的传播路径不再是直线而是一条曲线,从而导致传感器接收的像点发生位移。 2.进行几何校正并保证精度的必要性 遥感图像几何校正的精确与否直接关系到应用遥感信息反应地表地物的地理位置和面积的精确度,关系到从图像上获取的信息准确与否,因此在选择控制点上要十分小心,尽可能提高其精度,并且要对校正结果进行反复的分析比较,必要时还要进行多次校正。几何校正让图像上地物对应的像元出现在它应该在的地方,再通过辐射校正、影像增强等遥感图像处理技术,还图像以“本来面目”。然后通过对图像的识别、分类、解译处理实现地面空间上各类资源信息的空间分析研究,使遥感技术投入到实际生产应用中。 3.几何精校正 遥感影像图的几何校正目前有3种方案,即系统校正、利用控制点校正以及混合校正。遥感数据接收后,首先由接收部门进行校正,这种校正叫系统校正(又叫几何粗校正) ,即把遥感传感器的校准数据、传感器的位置、卫星姿态等测量值代入理论校正公式进行几何畸变校正;而用户拿到这种产品后,由于使用目的不同或投影及比例尺不同,仍旧需要做进一步的几何校正,这就需要对其进行几何精校正即利用地面控制点GCP ( GroundContr ol2Point,遥感图像上易于识别,并可精确定位的点)对因其他因素引起的遥感图像几何畸变进行纠正。混合校正则是由一般地面站提供的遥感CCT已经完成了第一阶段的几何粗校正,用户所要完成的仅仅是对图像做进一步的几何精校正。 几何精校正就是利用地面控制点GCP对各种因素引起的遥感图像几何畸变进行校正。从数学上说,其原理是通过一组GCP建立原始的畸变图像空间与校正空间的坐标变换关系,

遥感实验报告

遥感原理与应用 实验报告 姓名:学号:学院:专业: 年月日 实验一: erdas视窗的认识实验 一、实验目的 初步了解目前主流的遥感图象处理软件erdas的主要功能模块,在此基础上,掌握几个视窗操作模块的功能和操作技能,为遥感图像的几何校正等后续实习奠定基础。 二、实验步骤 打开imagine 视窗 启动数据预处理模块 启动图像解译模块 启动图像分类模块 imagine视窗 1.数据预处理(data dataprep) 2.图像解译(image interpreter) 主成份变换 色彩变换 3.图像分类(image classification) 非监督分类 4. 空间建模(spatial modeler) 模型制作工具 三、实验小结 通过本次试验初步了解遥感图象处理软件erdas的主要功能模块,在此基础上,基本掌握了几个视窗操作模块的功能和用途。为后续的实验奠定了基础。 实验二遥感图像的几何校正 掌握遥感图像的纠正过程 二、实验原理 校正遥感图像成像过程中所造成的各种几何畸变称为几何校正。几何校正就是将图像数据投影到平面上,使其符合地图投影系统的过程。而将地图投影系统赋予图像数据的过程,称为地理参考(geo-referencing)。由于所有地图投影系统都遵循一定的地图坐标系统,因此几何校正的过程包含了地理参考过程。 几何校正包括几何粗校正和几何精校正。地面接收站在提供给用户资料前,已按常规处理方案与图像同时接收到的有关运行姿态、传感器性能指标、大气状态、太阳高度角对该幅图像几何畸变进行了几何粗校正。利用地面控制点进行的几何校正称为几何精校正。一般地面站提供的遥感图像数据都经过几何粗校正,因此这里主要进行一种通用的精校正方法的实验。该方法包括两个步骤:第一步是构建一个模拟几何畸变的数学模型,以建立原始畸变图像空间与标准图像空间的某种对应关系,实现不同图像空间中像元位置的变换;第二步是利用这种对应关系把原始畸变图像空间中全部像素变换到标准图像空间中的对应位置上,完成标准图像空间中每一像元亮度值的计算。 三、实验内容 根据实验的数据,对两张图片进行几何纠正 四、实验流程

几何校正操作步骤(精)

几何校正操作步骤 实验目的: 通过实习操作,掌握遥感图像几何校正的基本方法和步骤,深刻理解遥感图像几何校正的意义。 实验内容: ERDAS软件中图像预处理模块下的图像几何校正。几何校正就是将图像数据投影到平面上,使其符合地图投影系统的过程。而将地图投影系统赋予图像数据的过程,称为地里参考(Geo-referencing)。由于所有地图投影系统都遵循一定的地图坐标系统,因此几何校正的过程包含了地理参考过程。 1、图像几何校正的途径 ERDAS图标面板工具条:点击DataPrep图标,→Image Geometric Correction →打开Set Geo-Correction Input File对话框(图2-1)。 ERDAS图标面板菜单条:Main→Data Preparation→Image Geometric Correction→打开Set Geo-Correction Input File对话框(图2-1)。 在Set Geo-Correction Input File对话框(图1)中,需要确定校正图像,有两种选择情况: 其一:首先确定来自视窗(FromViewer),然后选择显示图像视窗。 其二:首先确定来自文件(From Image File),然后选择输入图像。 2、图像几何校正的计算模型(Geometric Correction Model) ERDAS提供的图像几何校正模型有7种,具体功能如下:

3、图像校正的具体过程 第一步:显示图像文件(Display Image Files) 首先,在ERDAS图标面板中点击Viewer图表两次,打开两个视窗(Viewer1/Viewer2),并将两个视窗平铺放置,操作过程如下: ERDAS图表面板菜单条:Session→Title Viewers 然后,在Viewer1中打开需要校正的Lantsat图像:xiamen,img 在Viewer2中打开作为地理参考的校正过的(图象或)矢量图层:xmdis3.shp 第二步:启动几何校正模块(Geometric Correction Tool)Viewer1菜单条:Raster→Geometric Correction →打开Set Geometric Model对话框(2-2) →选择多项式几何校正模型:Polynomial→OK →同时打开Geo Correction Tools对话框(2-3)和Polynomial Model Properties对话框(4)。 在Polynomial Model Properties对话框中,定义多项式模型参数以及投影参数:→定义多项式次方(Polynomial Order)(图2-4):2 →定义投影参数:(PROJECTION):略 →Apply→Close →打开GCP Tool Referense Setup 对话框(2-5)

实验三 遥感图像的几何校正

实验法三遥感图像的几何校正 一实验目的 通过实验操作,掌握遥感图像几何校正的基本方法和步骤,深刻理解遥感图像几何校正的意义。 二实验内容 ERDAS软件中图像预处理模块下的图像几何校正。 几何校正就是将图像数据投影到平面上,使其符合地图投影系统的过程。而将地图投影系统赋予图像数据的过程,称为地理参考(Geo-referencing)。由于所有地图投影系统都遵循一定的地图坐标系统,因此几何校正的过程包含了地理参考过程。 1、图像几何校正的途径 ERDAS图标面板工具条:点击DataPrep图标,→Image Geometric Correction →打开Set Geo-Correction Input File对话框(图1)。 ERDAS图标面板菜单条:Main→Data Preparation→Image Geometric Correction→打开Set Geo-Correction Input File对话框(图1)。 图1 Set Geo-Correction Input File对话框 在Set Geo-Correction Input File对话框(图1)中,需要确定校正图像,有两种选择情况: 其一:首先确定来自视窗(From Viewer),然后选择显示图像视窗。 其二:首先确定来自文件(From Image File),然后选择输入图像。 2、图像几何校正的计算模型(Geometric Correction Model) ERDAS提供的图像几何校正模型有7种,具体功能如下: 表1 几何校正计算模型与功能 模型功能 Affine 图像仿射变换(不做投影变换) Polynomial 多项式变换(同时作投影变换) Reproject 投影变换(转换调用多项式变换) Rubber Sheeting 非线性变换、非均匀变换 Camera 航空影像正射校正 Landsat Lantsat卫星图像正射校正 Spot Spot卫星图像正射校正 其中,多项式变换(Polynomial)在卫星图像校正过程中应用较多,在调用多项式模型时,需要确定多项式的次方数(Order),通常整景图像选择3次方。次方数与所需要的最

遥感实验报告

重庆交通大学 学生实验报告 实验课程名称遥感原理与应用 开课实验室测量与空间信息处理实验室 学院 2013 年级测绘工程专业 1班学生姓名刘文洋 学号 631301040126 开课时间 2015 至 2016 学年第 1 学期

目录 实验一 ENVI 视窗的基本操作 (2) 实验二遥感图像的几何校正 (4) 实验三遥感图像的增强处理 (8) 实验四遥感图像的变换 (12) 实验五遥感信息的融合 (15) 实验六遥感图像分类 --- 监督分类 (17) 实验七遥感图像分类 --- 非监督分类 (19) 实验八遥感图像分类后处理 (22)

实验一ENVI 视窗的基本操作 一、实验目的 初步了解目前主流的遥感图象处理软件 ENVI 的主要功能模块,在此基础上,掌握视窗操作模块的功能和操作技能,为遥感图像的几何校正等后续实习奠定基础。 二、实验内容 视窗功能介绍;文件菜单操作;显示数据;裁剪数据;合并波段 三、实验步骤 1、首先打开ENVI4.7软件,看见的只有菜单栏,如图所示: 2、打开每个下拉菜单浏览其下拉栏中都有哪些功能,比如:我们如果需要打开遥感文件,则可以选择File下的打开功能open image file,打开遥感图像如下图:

裁剪数据打开basic tools的resize data功能,如果需要对图像进行一系列处理,可以利用Transform,Classification等功能进行操作,在后续实验中我们也会用到其中的一些功能进行图像的一系列操作,到时候在详细叙述。 3、再熟悉了ENVI4.7的一些基本知识后我们可以简单地操作下,比如对一组数据分别用Gray Scale和Load RGB导入,看看两幅图的区别以及各自的优缺点。 四、实验结果分析 在这次的实验中,我们简单的熟悉了下ENVI4.7的一些功能,发现它是可以对遥感图像进行图像几何纠正,直方图均衡,监督分类,非监督分类等一系列操作,为我们后续利用软件对遥感图像处理打下了基础。

ERDAS遥感图像的几何校正

遥感图像的几何校正 实验目的:通过实习操作,掌握遥感图像几何校正的基本方法和步骤,深刻理解遥感图像几何校正的意义。 实验内容:ERDAS软件中图像预处理模块下的图像几何校正。 几何校正就是将图像数据投影到平面上,使其符合地图投影系统的过程。而将地图投影系统赋予图像数据的过程,称为地理参考(Geo-referencing)。由于所有地图投影系统都遵循一定的地图坐标系统,因此几何校正的过程包含了地理参考过程。 1、图像几何校正的途径 ERDAS图标面板工具条:点击DataPrep图标,→Image Geometric Correction →打开Set Geo-Correction Input File对话框(图1)。 ERDAS图标面板菜单条:Main→Data Preparation→Image Geometric Correction→打开Set Geo-Correction Input File对话框(图1)。 图1 Set Geo-Correction Input File对话框 在Set Geo-Correction Input File对话框(图1)中,需要确定校正图像,有两种选择情况: 其一:首先确定来自视窗(FromViewer),然后选择显示图像视窗。 其二:首先确定来自文件(From Image File),然后选择输入图像。 2、图像几何校正的计算模型(Geometric Correction Model) ERDAS提供的图像几何校正模型有7种,具体功能如下:

3、图像校正的具体过程 第一步:显示图像文件(Display Image Files) 首先,在ERDAS图标面板中点击Viewer图表两次,打开两个视窗(Viewer1/Viewer2),并将两个视窗平铺放置,操作如下:ERDAS图表面板菜单条:Session→Title Viewers 然后,在Viewer1中打开需要校正的Lantsat图像:tmatlanta.img 在Viewer2中打开作为地理参考的校正过的SPOT图像:panatlanta.img 第二步:启动几何校正模块(Geometric Correction Tool) Viewer1菜单条:Raster→Geometric Correction →打开Set Geometric Model对话框,如图2

实验二遥感图像的几何校正与镶嵌实验报告

实验二遥感图像的几何校正与镶嵌实验报告 实验目的: 通过本实验熟练操作遥感图像处理的专业软件进行基础图像处理,包括图像几何校正、镶嵌等。 实验内容: 1、熟悉图像几何校正、镶嵌的基本原理; 2、学习图像几何校正具体操作; 3、学习图像镶嵌正具体操作。 本实验的图像几何校正是通过“像图配准”的方式获取地面控制点的方里网坐标的,并对传统的从纸质地形图上量算坐标的方法进行改进,利用Auto CAD或Photoshop等软件从扫描后的电子地形图上直接量算坐标。 实验步骤: 第一步、熟悉图像几何校正、镶嵌的基本原理 第二步、图像几何校正 运行PCI,选择GCPWorks模块,在Source of GCPs选择User Entered Coordinates(用户输入投影坐标系统),点击Accept后,弹出校正模块:

选择第一项加载需要校正的图像(由实验一方法导出的125-42.pix )->点击Default->Load & Close->得到下图:

选择第二项,选择Other确定投影系统: 注意输入6度带的中央经度与向东平移500公里(500000米):

点击Earth Model确定地球模型: 点击Accept:

选择第三项采集地面控制点。在采集地面控制点之前,利用Photoshop软件打开扫描后的电子地形图。 分别在遥感图像和地形图中找到一个同名点,如下图(可以用放大遥感图)。 然后在地形图中量算出该点的大地坐标,精确到米,X坐标为6位(要去掉2位6度带的带号),Y坐标7位(运用测出)。再将坐标输入到GCP编辑窗口中,并点击Accept as GCP接受为一个控制点。

数字图像处理-畸变校正

数字图像处理

图像畸变及校正 1 图像畸变介绍 从数字图像处理的观点来考察畸变校正, 实际上是一个图像恢复的过程, 是对一幅退化了的图像进行恢复。在图像处理中,图像质量的改善和校正技术,也就是图像复原,当初是在处理从人造卫星发送回来的劣质图像的过程中发展、完善的。目前,图像畸变校正的应用领域越来越广,几乎所有涉及应用扫描和成像的领域都需要畸变校正。图像在生成和传送的过程中,很可能会产生畸变,如:偏色、模糊、几何失真、几何倾斜等等。前几种失真主要是体现在显示器上,而后一种失真则多与图像集角度有关。不正确的显影,打印、扫描,抓拍受反射光线的影响等方式,都会使图像产生偏色现像。模糊、几何畸变主要是在仪器采集图片过程中产生,大多是因机器故障或操作不当影响导致,如在医学成像方面。而几何空间失真广泛存在于各种实际工程应用中,尤其是在遥感、遥测等领域。 2 畸变产生的原因 在图像的获取或显示过程中往往会产生各种失真(畸变):几何形状失真、灰度失真、颜色失真。引起图像失真的原因有:成像系统的象差、畸变、带宽有限、拍摄姿态、扫描非线性、相对运动等;传感器件自身非均匀性导致响应不一致、传感器件工作状态、非均匀光照条件或点光源照明等;显示器件光电特性不一致;图像畸变的存在影响视觉效果,也是影响图像检测系统的形状检测和几何尺寸测量精度的重要因素之一。 3 图像畸变校正过程所用到的重要工具 灰度直方图是关于灰度级分布的函数,是对图象中灰度级分布的统计。灰度直方图是将数字图象中的所有像素,按照灰度值的大小,统计其所出现的频度。

通常,灰度直方图的横坐标表示灰度值,纵坐标为想像素个数。直方图上的一个点的含义是,图像存在的等于某个灰度值的像素个数的多少。这样通过灰度直方图就可以对图像的某些整体效果进行描述。从数学上讲,图像的灰度直方图是图像各灰度值统计特征与图像灰度值出现的频率。从图形上来讲,它是一个一维曲线,表征了图像的最基本的统计特征。 作为表征图像特征的信息而在图像处理中起着重要的作用。由于直方图反映了图像的灰度分布状况,所以从对图像的观察与分析,到对图像处理结果的评价,灰度直方图都可以说是最简单、最有效的工具。 4 图像颜色畸变校正介绍 图像颜色畸变现象可以是由摄像器材导致,也可以是由于真实环境本身就偏色导致,还有的是由于图像放置过久氧化、老化导致。无论其产生的原因如何,其校正方法都是类似的。 如果用Matlab显示颜色畸变的图像RGB基色直方图,发现相对正常图像,颜色畸变的图像的直方图的三种基色的直方图中至少有一个直方图的像素明显集中集中在一处,或则集中在0处或则集中在255处,而另一部分有空缺,或则集中在中间而两边空,因此通过调整该直方图的像素点的像素值在区间[0,255]上的分布来解决图像颜色畸变问题。如果直方图中像素集中在0一边则说明该基色偏暗,如果集中在255处则说明该基色偏亮。下图是一有颜色畸变的图像的基色B 的直方图。

多时相遥感图像配准实验报告

Harbin Institute of Technology 多时相遥感图像配准 实验报告 课程名称:遥感信息处理导论 院系:电子与信息工程学院 姓名: 学号: 授课教师: 哈尔滨工业大学

1. 实验目的 对于通过遥感平台获取的遥感数字影像,由于扫描过程中受地球曲率、地球自转、平台姿态、扫描方式等因素的影响,所获取的图像往往会产生不同程度的几何误差,这些误差如若不经处理,对图像的后续应用如分类、目标检测等会产生很大影响。本实验的目的就是利用多项式映射和重采样等方法对两幅不同时间采集的遥感图像进行几何校正处理,从而方便从配准图像中找到不同时相内发生变化的地物目标信息。 2. 实验原理 图像的配准过程主要可以分为三个过程:控制点对选取(自动或手动)、坐标映射函数拟合、映射后像素重采样。 1) 控制点对选取 对于图像匹配过程,首先要通过控制点对建立两幅图像间的坐标联系。控制点的选取有可分为自动和手动两种方式。对于自动选取,可以采用SIFT 等算法实现;对于手动选取,需要实验人员通过目测观察的方法找到不同时相的遥感图像中的相同目标,这些控制点可以是建筑物的顶点、道路等的交叉点等或其他较容易分辨的位置。 选出控制点后将这些点的坐标信息记录下来,值得注意的是所选控制点的数量、分布情况以及精度会直接影响配准结果的精度和质量。 2) 图像几何校正 图像的配准实质上就是通过图像的几何校正,将产生几何失真的图像转换为标准的数字图像。采用的方法就是通过一定的映射函数将原图像的像素坐标转换为标准图像中的坐标 (,)u f x y = (,)v g x y = 其中(,)x y 表示原图像中的像素坐标值,(,)u v 表示参考图像中的像素坐标值。 映射函数的选择可以有多种形式,一般较为简单常用的是多项式函数,以二次多项式函数为例,映射函数的具体形式为 22 01234522 012345u a a x a y a xy a x a y v b b x b y b xy b x b y ?=+++++?=+++++? 其中,05, ,a a ,05,,b b 分别表示二次多项式横纵坐标映射函数中对应项的系 数,通过选定控制点,可以将这些系数项求解出来。 假设在两幅图像中总共选择了n 对控制点,根据n 对控制点的坐标 (,)~(,)i i i i x y u v 可以得到如下两个方程组,

ERDAS 遥感图像的几何校正

学号: 遥感数字图像处理软件 实验报告 (2011~2012学年第二学期) 学院:地环学院 班级: 09地科2 姓名: 指导老师:

实验二:遥感图像的几何校正 一.实验平台:ERDAS IMAGINE 9.1 二.实验目的 通过实验操作,掌握遥感图像几何校正的基本方法和步骤,能够理解遥感图像的几何校正的意义。 三.实验内容 资源卫星数据的校正;图像的放射变换;航片的正射校正;图像的镶嵌;图像投影变换等内容。 四.实验步骤 1.资源卫星数据landsat 的校正 1.1打开图像文件 在ERDAS 图标面板中打开两个视窗(Viewer#1 和Viewer#2);在ERDAS 图标面板菜单条点击Session,点击 Tile Viewers 并将Viewer#1 和Viewer#2 平铺放置;在Viewer#1 中打开需要校正的Landsat TM 图像:tmAtlanta.img (图1),在Viewer#2 中打开作为地理参考的校正过的SPOT 图像:panAtlanta.img(图2)。 图 1 图 2 1.2 启动几何校正模块 (1)在Viewer#1 视窗菜单条中点击Raster (2)点击Geometric Correction,打开 Set Geometric Model 对话框(图3)。 图 3 图 4 图 5 (3)在Set Geometric Model 对话框中选择多项式几何校正模型:Polynomial

(4)同时打开Geo Correction Tools(图4)和Polynomial Model Properties(图5)对话框. (5)在Polynomial Model Properties 对话框中定义多项式模 型参数及投影参数:定义多项式次方: 2. (6)点击Apply (7)点击Close ,打开GCP Tool Reference Setup 对话框(图6) 图 6 1.3 启动控制点工具 (1)在GCP Tool Reference Setup 对话框中选择视窗采点模式Existing Viewer。 (2)点击OK,打开Viewer Selection Instructions 指示器。 (3)在显示作为地理参考图像panAtlanta.img 的Viewer#2 中点击左键,打开Reference Map Information 提示框(图7),显示参考图像的投影信息。 图 7 (4)点击OK,整个屏幕自动变为两个主视窗、两个放大窗口(Viewers #3 和#4)、两个关联方框(分别位于两个视窗中,指示放大视窗与主视窗的关系)、控制点工具对话框、几何校正工具等(图8)。表面控制点工具被启动,进入控制点采集状态。 图 8

ENVI实验报告

一、实验目的 ENVI是一套功能齐全的遥感图像处理系统,是处理、分析并显示多光谱数据、高光谱数据 和雷达数据的高级工具。此次实习主要是学习一些关于ENVI的基本操作,如:图像预处理,影像分析,图像增强,几何校正,监督分类以及专题制图等步骤。 二、实验数据 ELEVATION_SOURCE = "GLS2000" EPHEMERIS_TYPE = "DEFINITIVE" SPACECRAFT_ID = "Landsat7" SENSOR_ID = "ETM+" SENSOR_MODE = "SAM" ACQUISITION_DATE = 2000-09-24 WRS_PATH = 144 三、实验内容 1.遥感数据下载 2.数据格式转化及多波段的合成 3.图像边框的裁剪 4.图像的镶嵌 6.图像几何校正 7.裁剪 8.监督分类 9.专题出图 四、实验步骤 1.遥感数据下载 输入地名及日期,搜索符合自己通途的数据,查看并下载

2.数据格式转化及多波段的合成 由于下载图像的格式tif格式,需要转成img格式。并将单波段融合成为多波段,这一步在erdas中完成。 tif转img 单波段融合成多波段 3.图像边框的裁剪 通过主图像窗口-tools-region?of?interest-roi?tool,然后自己在图像上定义出一块感兴趣区(实际中难于确定严格准确的裁剪界限)。 4.图像的镶嵌 Map→mosaicking→georeferenced 启动镶嵌模板 6.图像几何校正 设置参数,参数在原始数据中的txt中可以查看 启动校正模块,并利用google地球进行校正 7.裁剪 分别打开需要裁剪的图像和所在区域的shapefile文件 将矢量文件转成ROI 利用ROI裁剪图像 选择主菜单→basic tools→subset data via ROIS,选择裁剪图像 在spatial subset via ROI parameters中,选择由矢量生成的

实验二 遥感图像的几何校正

实验二、遥感图像的几何校正 实验目的:通过实习操作,掌握遥感图像几何校正的基本方法和步骤,深刻理解遥感图像几何校正的意义。 实验内容:ERDAS软件中图像预处理模块下的图像几何校正。 几何校正就是将图像数据投影到平面上,使其符合地图投影系统的过程。而将地图投影系统赋予图像数据的过程,称为地里参考(Geo-referencing)。由于所有地图投影系统都遵循一定的地图坐标系统,因此几何校正的过程包含了地理参考过程。 1、图像几何校正的途径 ERDAS图标面板工具条:点击DataPrep图标,→Image Geometric Correction →打开Set Geo-Correction Input File对话框(图2-1)。 ERDAS图标面板菜单条:Main→Data Preparation→Image Geometric Correction→打开Set Geo-Correction Input File对话框(图2-1)。 图2-1 Set Geo-Correction Input File对话框 在Set Geo-Correction Input File对话框中,选择输入图像,确定校正图像。 2、图像几何校正的计算模型(Geometric Correction Model) ERDAS提供的图像主要几何校正模型,具体功能如下:

表2-1 几何校正计算模型与功能 模型功能 Affine 图像仿射变换(不做投影变换) Polynomial 多项式变换(同时作投影变换) Reproject 投影变换(转换调用多项式变换) Rubber Sheeting 非线性变换、非均匀变换 Camera 航空影像正射校正 Landsat Lantsat卫星图像正射校正 Spot Spot卫星图像正射校正 3、图像校正的具体过程 第一步:显示图像文件(Display Image Files) 首先,在ERDAS图标面板中点击Viewer图表两次,打开两个视窗(Viewer1/Viewer2),并将两个视窗平铺放置,操作过程如下: 在Viewer1中打开需要校正的图像(或通过图2-1已打开):tmAtlanta.img 在Viewer2中打开作为地理参考的校正过的图像:panAtlanta.img 第二步:启动几何校正模块(Geometric Correction Tool) Viewer1菜单条:Raster→ Geometric Correction →打开Set Geometric Model对话框(2-2) →选择多项式几何校正模型:Polynomial→OK →同时打开Geo Correction Tools对话框(2-3)和Polynomial Model Properties对话框(2-4)。 在Polynomial Model Properties对话框中,定义多项式模型参数以及投影参数: →定义多项式次方(Polynomial Order):2 →定义投影参数:(PROJECTION):略 →Apply→Close →打开GCP Tool Referense Setup 对话框(2-5)

ERDAS几何校正

图错误!文档中没有指定样式的文字。-1 Erdas几何校正 几何校正的具体过程: 1 显示待校正影像和参考影像 在ERDAS图标面板中点击Viewer图表两次,打开两个视窗(Viewer1/Viewer2),并将两个视窗平铺放置,操作过程如下: ERDAS图表面板菜单条:Session→Title Viewers 然后,在Viewer1中打开需要校正的影像 Viewer2中打开参考的影像 2启动几何校正模块 Viewer1菜单条:Raster→Geometric Correction →打开Set Geometric Model对话框 根据影像类型和要求选择不同的校正模型。 图一Set Geometric Model对话框 →选择多项式几何校正模型:Polynomial→OK →同时打开Geo Correction Tools对话框(图二)和Polynomial Model Properties对话框(图三)。 在Polynomial Model Properties对话框中,定义多项式模型参数以及投影参数:

定义多项式次方(Polynomial Order):2 定义投影参数:不需定义 →Apply→Close →点击打开GCP Tool Referense Setup 对话框(图四) 图二图三 图四图五 3 启动控制点工具 首先,在GCP Tool Referense Setup对话框(图四)中选择采点模式:→选择视窗采点模式:Existing Viewer→OK

→打开Viewer Selection Instructions指示器(图五) →在显示作为参考图像的Viewer2中点击左键 →打开reference Map Information 提示框(图六);→OK →此时,整个屏幕将自动变化为如图七所示的状态,表明控制点工具被启动,进入控制点采点状态。 图六 图七 4 采集地面控制点 GCP的具体采集过程: 在图像几何校正过程中,采集控制点是一项非常重要和繁重的工作,具体过程如下: 1、在GCP工具对话框中,点击Select GCP图表,进入GCP选择状态; 2、在GCP数据表中,将输入GCP的颜色设置为比较明显的黄色。 3、在Viewer1中移动关联方框位置,寻找明显的地物特征点,作为输入GCP。 4、在GCP工具对话框中,点击Create GCP图标,并在Viewer3中点击左键定点,

遥感图像预处理实验报告

实验前准备:遥感图像处理软件认识 1、实验目的与任务: ①熟悉ENVI软件,主要是对主菜单包含内容的熟悉; ②练习影像的打开、显示、保存;数据的显示,矢量的叠加等。 2、实验设备与数据 设备:遥感图像处理系统ENVI4.4软件; 数据:软件自带数据和河南焦作市影响数据。 3、实验内容与步骤: ⑴ENVA软件的认识 如上图所示,该软件共有12个菜单,每个菜单都附有下拉功能,里面分别包含了一些操作功能。 ⑵打开一幅遥感数据 选择File菜单下的第一个命令,通过该软件自带的数据打开遥感图像,可知,打开一幅遥感影像有两种显示方式。一种是灰度显示,另一种是RGB显示。 Gray(灰度显示)RGB显示 ⑶保存数据 ①选择图像显示上的File菜单进行保存; ②通过主菜单上的Save file as进行保存

⑷光谱库数据显示 选择Spectral > Spectral Libraries > Spectral Library Viewer。将出现Spectral Library Input File 对话框,允许选择一个波谱库进行浏览。点 击“Open Spectral Library”,选择某一所需的 波谱库。该波谱库将被导入到Spectral Library Input File 对话框中。点击一个波谱库的名称, 然后点击“OK”。将出现Spectral Library Viewer 对话框,供选择并绘制波谱库中的波谱曲线。 ⑸矢量化数据 点选显示菜单下的Tools工具栏,接着选择下面的第四个命令,之后选择第一个命令,对遥感图像进行矢量化。点击鼠标左键进行区域选择,选好之后双击鼠标右键,选中矢量化区域。 ⑹矢量数据与遥感影像的叠加与切割 选择显示菜单下的Tools工具,之后点选第一个 Link命令,再选择其下面的第一个命令,之后 OK,结束程序。 选择主菜单下的Basic Tools 菜单,之后选择 其中的第二个命令,在文件选择对话框中,选择 输入的文件(可以根据需要构建任意子集),将 出现Spatial Subset via ROI Parameters 对 话框通过点击矢量数据名,选择输入的矢量数 据。使用箭头切换按钮来选择是否遮蔽不包含在 矢量数据中的像元。 遥感图像的辐射定标 1、实验目的与任务: ①了解辐射定标的原理; ②使用ENVI软件自带的定标工具定标; ③学习使用波段运算进行辐射定标。 2、实验内容与步骤: ⑴辐射定标的原理 辐射定标就是将图像的数字量化值(DN)转化为辐射亮度值或者反射率或者表面温度等

遥感图像几何精校正实验报告

遥感图像几何精校正 实验名称:遥感图像的几何精校正。 实验目的:1.了解和熟悉envi软件的几何校正的原理 2.熟悉和掌握envi软件的几何校正的功能和使用方法; 3.对自己的图像先找到投影,再另存一幅图像,去掉投影,在其它软件中旋转一 角度,用原先的图像作为参考对旋转后的图像进行几何校正,使得其比较精确。实验原理:几何校正,主要方法是采用多项式法,机理是通过若干控制点,建立不同图像间的多项式控件变换和像元插值运算,实现遥感图像与实际地理图件间的配准,达 到消减以及消除遥感图像的几何畸变。 多项式几何校正激励实现的两大步: 1. 图像坐标的空间变换: 有几何畸变的遥感图像与没有几何畸变的遥感图像,其对应的像元的坐标是不一 样的,如下图1右边为无几何畸变的图像像元分布图,像元是均匀且不等距的分 布。为了在有几何畸变的图像上获取无几何畸变的像元坐标,需要进行两图像坐 标系统的空间装换。 图1:图像几何校正示意图 在数学方法上,对于不同二维笛卡儿坐标系统间的空间转换,通常采用的是二元 n次多项式,表达式如下: 其中x, y为变换前图像坐标, u, v为变换后图像坐标, aij , bij为多项式系数, n = 1, 2, 3, ?。 二元n次多项式将不同坐标系统下的对应点坐标联系起来, ( x, y )和( u, v )分别应 不同坐标系统中的像元坐标。这是一种多项式数字模拟坐标变换的方法,一旦有 了该多项式,就可以从一个坐标系统推算出另一个坐标系统中的对应点坐标。 如何获取和建立二元n次多项式,即二元n次多项式系数中a和b的求解,是几何 校正成败的关键。数学上有一套完善的计算方法,核心是通过已知若干存在于不 同图像上的同名点坐标,建立求解n次多项式系数的方程组,采用最小二乘法,得出 二元n次多项式系数。 不同的二元n次多项式,反映了几何畸变的遥感图像与无几何畸变的遥感图像间的 像元坐标的对应关系, 其中哪种多项式是最佳的空间变换模拟式,能达到图像间 坐标的完全配准,是需要考虑和分析的。 在二元n次多项式数字模拟中,从提高几何校正精度的角度考虑,需要兼顾的因素

实验四 几何校正

实验四 几何校正 图像校正主要指辐射校正和几何校正。 1、辐射校正包括传感器的辐射校正、大气校正、照度校正以及条纹和斑点的判定和消除。 2、几何校正就是校正成像过程中造成的各种几何畸变,包括几何粗校正和几何精校正。 图像几何校正的一般步骤: 掌握遥感图像几何校正的方法,利用Data Preparation 模块通过采集控制点实现图像的几何校正 第一步 数据和校正模型的准备 1.请直接点击桌面IMAGINE 图标,等待Viewer1出现。 2.由Viewer1打开开启欲校正的图像 C: \linchuan\linchuang-TM.img 。 3.点击 Viewer 图标,等待Viewer2出现。读取已校正过的参考图像C: \linchuan\linchuang-geo.img 。 数据准备 输入显示数字影像 确立校正变换模型 确定输出影像范围 像元空间坐标变换 像元的灰度重采样 输出纠正数字影像

4.由主菜单中,点击"DataPrep"图标,选择其中之Image Geometric Correction。 5.在出现的Set Geo Correction Input File对话框中,点击"Select Viewer"此选项,然后将鼠标光标对Viewer1点一下。 6.接着屏幕会出现Set Geometric Model之对话框,请选取第二项Polynomial,按 7.在出现的Polynomial Model Properties对话框中,如下图设置

因为参考图像panAtlanta.img 已经含有投影参数,故不再需要定义投影参数 请按下apply,然后按下close关闭对话框。 8.接着会出现GCP Tool Reference Setup之对话框,选择校正参考坐标的来源,由于先前已开启Viewer2图像,故于此选取第一项Existing Viewer,按下OK,屏幕会出现Viewer Selection Instructions对话框,

基于ERDAS航空影像单片正射纠正指南

基于ERDAS航空影像的单片正射纠正 王芳霞 摘要:介绍利用遥感图像处理系统ERDAS IMAGINE对单张航片进行正射纠正的原理、方法,纠正过程中应注意的问题及用此方法生产正射影像的优点和应用前景。 关键词:正射影像纠正遥感影像处理 1引言 随着数字化技术和城市建设的飞速发展,传统的测绘产品已不能满足各个行业高速发展的需要,高分辩率实时性好的数字测绘产品已逐渐替代了原先的传统测绘产品。数字正射影像地图就是其中一种重要的产品,所谓正射影像,指改正了因地形起伏和传感器误差而引起的像点位移的影像。数字正射影像不仅精度高,信息丰富,直观真实,而且数据结构简单,生产周期短,能很好的满足社会各行业的需要。利用摄影测量的方法生产正射影像,要求有准确的外业控制资料,且耗时耗力,遇到某地区没有现成的DEM,又没有带高程信息的地形图可供利用时,其不失为一种很好的方法。但若有现成DEM可供利用,则可采用单片数字正射纠正方案,该方案不仅可以省一道很费人力物力的工序,而且还可根据相片本身的重叠度,进行隔片纠正,从生产成本和速度上大大提高了生产效率,本文就是基于这种思路,介绍基于遥感图像处理系统ERDAS来进行航片正射纠正生产正射影像的方法。 2正射纠正的原理和方法

2.1原理: 数字正射纠正的实质就是将中心投影的影像通过数字元纠正形成正射投影的过程,其原理是将影像化为很多微小的区域,根据有关的参数利用相应的构像方程式或按一定的数学模型用控制点解算,求得解算模型,然后利用数字元高程模型对原始非正射影像进行纠正,使其转换为正射影像。注意纠正时尽量利用影像中心区域的影像,而避免利用影像边缘的影像。 2.2正射纠正的工艺流程如下:

相关文档
最新文档