to 实验十二 Image 遥感图像几何精校对 Map形式

本文目录导航：

• 实验十二 遥感图像几何精校对：Image to Map形式
• 实验三 遥感图像空间测量
• matlab图像解决教程

实验十二 遥感图像几何精校对：Image to Map形式

一、实验目的

经过ENVI的Image to Map Registration配置命令经常使用，学会用带天文坐标的地形图（Map）去校对几何畸变的遥感影像，使后者也成为一幅无畸变的带天文坐标的遥感图像，从而加深对遥感图像几何校对的了解。

二、实验内容

①桂林市TM 影像的Image和Map管理点GCPs选用；②Image和Map的管理点坐标提取和输入；③基于Image to Map的图像空间配准；④插值算法选用；⑤坐标投影确定；⑥Wrap操作；⑦遥感图像坐标网笼罩。

三、实验要求

①课前预备好地形图上的管理点坐标数据；②管理点选用必定是其所在位置具备明晰可辨的点位识别特色，如点、两线交点和复线锐角拐点等；③管理点在影像上要平均散布；④全图区内的管理点数量不能少于20个；⑤对校对前后的图像启动比拟剖析。

四、技术条件

①微型计算机；②桂林市TM 卫星影像数据；③ENVI软件；④Photoshop软件（ver.6.0以上）和ACDSee软件（ver.4.0以上）。

五，实验步骤

（1）图像输入。

在ENVI主菜单栏当选用“File>Open Image File”，产生文件目录窗口，将桂林市TM 1～7波段的数据调入“Available Bands List”窗口，从波段列表当选用TM1显示之。

（2）几何校对打算选用。

在ENVI主菜单栏中，选用“Map”，在其下拉菜单中，选用“Map>Registration>Select GCPs: ImagetoMap”，如图12-1所示。

（3）坐标投影确定。

在选用“Select GCPs:Image to Map”之后，产生“Image to Map Registration”窗口。

此框中列出了各种坐标投影，从当选用本实验驳回的投影及地图分带“GK Zone18”（高斯克吕格投影18度带）；并点击关上“Datum..”，产生“Select Geogprhic Datum”（选用天文坐标系统窗口），从当选用“WGS-84”；并在“Image to Map Registration”窗口中，对影像空间分辨率尺寸“X Pixel Size”和“Y Pixel Size”都输入所解决遥感影像的空间分辨率。

这里是TM为30m。

而后，对此两个窗口区分点击【OK】按钮分开，如图12-2所示。

图12-1“Image to Map”解决打算选的操作界面

图12-2 几何校对坐标投影确定的操作界面

（4）输入管理点坐标数据。

在口头完步骤（3）后，产生“Ground Control Points Selection”（选管理点）窗口，如图12-3所示。

在此窗口中，左边的高低两个空白框是要输入的管理点的W GS-84坐标值（E对应于X坐标，N对应于Y坐标），左边的“Image X”和“Image Y”区分为与该管理点对应的遥感影像图像元样（Sample）、行（Line）坐标值（S,L）。

只能是手工输入（包含可以将WORD表格中的数据拷贝到此粘贴）。

参与新点输入或许修负数据可以经过操作【Add Point】按钮成功。

图12-3 几何校对坐标投影确定的操作界面

（5） Warp（纠正）解决图像选用。

在成功步骤（4）的管理点坐标数据输入之后，不用分开“Ground Control Points Selection”界面，间接在该窗口上方的菜单中，关上“Options”，产生下拉菜单，菜单最下面一项是对正在显示的图像作纠正，第二项是对磁盘文件中的图像作纠正，本实验选用第一项，如图12-4所示。

图12-4“Warp”（纠正）解决图像选用的操作界面

（6）口头Warp（纠正）解决。

在口头步骤（5）后，产生“Reigstration Parameters”（校对参数）框，如图12-5所示。

在该框中，左侧为前面各项操作构成的数据，普通不用再改变。

右侧的四行中，从上到下依次是：校对方法（Method）、重采样（Resampling）、背景值（Background）和输入校对图像的存储位置（Output Result to），对“Method”采取多项式（2次），对“Resampling”可选用“Nearest Neighbor”（最临近）、“Bilniear”（双线性插值）或“Cubic Convolution”（三次卷积）。

对存储位置，假设是实验解决，可以存在Memory；假设是正式结果，就可以作为文件存到磁盘上，如图12-5所示。

设置完这些校对参数后，点击【OK】按钮，成功Image to Map形式几何精校对。

图12-5 口头“Warp”（纠正）解决的操作界面

（7）用Photoshop软件将经过几何校对的图像与原始图像并列关上，观察比拟几何校对前后影像有何差异，用WORD文档记载，取名为《几何校对前后差异剖析》，存入自己的上班文件夹。

六、实验报告

（1）简述实验环节。

（2）回答疑问：①“Image to Map”是何含意？②运用“Image to Map”配置须要做哪些相应的预备上班？③参考实验十一的②题和③题。

实验报告格局见附录一。

实验三 遥感图像空间测量

一、实验目的

经过用像元计算地物之间的距离，了解遥感数字图像像元与空间分辨率的相关，以及遥感数字图像的几何构成及空间结构，把握基于像元的图像空间测量计算方法。

二、实验内容

（1）运用像元坐标及多边形面积计算公式测算恣意形态影像区域的面积；

（2）运用ENVI影像测量工具（Measurement Tool）测算恣意形态影像区域的面积。

三、实验要求

预习本实验，仔细观摩教员展示。

学会ENVI的Cursor Location/Value和Measurement Tool配置经常使用，实际了解像元含意，能灵敏运用手工测算法和ENVI影像测量工具法计算遥感影像上的距离和面积。

测量结果存档。

编写实验报告。

四、技术条件

①微型计算机；②桂林市TM 1～7波段数据；③ENVI软件；④Photoshop软件（ver.6.0以上）和ACDSee软件（ver.4.0以上）。

五、实验步骤

（1）建设用自己名字命名的实习公用文件夹。

（2）数据输入。

选用“File> Open Image File”，产生文件目录窗口，找到寄存桂林市TM 1 ～7波段数据的子目录并关上。

（3）影像遨游。

经过拉宽主窗口使影像区域获取最大限制的显示，用鼠标移动Scroll窗口中的白色矩形框，可以使主窗口中影像跟着移动，从而可以观察到整个影像区域的遥感影像。

（4）读取图像坐标。

用鼠标从上方横条配置菜单当选“Tools> Cursor Location/Value…”，屏幕上产生“Cursor Location/Value”对话框，其中显示四行数据，如图3-1所示。

Disp#1

——鼠标箭头所在处遥感数字图像的图像坐标，即Sa mple——样（有的也称列），Line— 行—。

图3-1 光标位置/数值窗口

Scrn:R：

——该坐标处三个色彩图像的灰度值（假设是单波段图像，此三值相反）。

Projection: Geographic Lat/Lon— 遥—感影像的地图投影：经纬度天文坐标。

只要经过几何校对的遥感影像，才正确显示出该点的经纬度天文坐标。

否则，只显示为：LL:0？′0.00″N,0？′0.00″E。

——该点的影像灰度值，数值范畴0～255。

在“Cursor Location/Value”形态下，鼠标箭头移动到图像区任何位置，这四组数据都会追随显示，因此，就可以随便读出感兴味指标位置的准确坐标。

（5）运用像元坐标及多边形面积计算公式手工计算恣意影像区域的面积：在成功上述步骤后，用“Tools/Cursor Location/Value”配置，对桂林市TM 假黑白分解影像图（图3-2）中红框的一切角点坐标启动量度，量度结果记入表3-1中。

在所有角点量度成功之后，将这些角点坐标值依照顺序代入公式（3-1），就可以计算出桂林市城区的面积。

表3-1 桂林市TM影像城市多边形区域角点坐标测算表

续表

图3-2 桂林市TM5、TM4和TM3波段假黑白分解增强影像图

据此，计算其南北向和物品向长度（结果用千米示意）及区域的面积（结果用平方千米示意）。多边形面积（S）计算公式为

遥感地质学实验教程

（6）运用ENVl的影像测量工具（Measurement Tool）启动面积测量。

1）在主窗口上方菜单下选用“Tools>MeasurementTool”，产生“Display Measurement Tool”对话框。

在“Display”文本框内，输入想测量的显示号。

选用适当的图像窗口切换按钮，包含主窗口（Image）、滚动窗口（Scroll）、缩放窗口（Zoom），要在任何期间禁用测量配置，选用【Off】切换按钮，如图3-3所示。

图3-3 显示测量工具对话框

2）选用“Type＞所需测量的区域形态”。

其中，Polygon为多边形；Polyline为折线；Rectangle为矩形；Ellipse为椭圆。

在图像显示内，经过用鼠标左键点击，绘制所须要的形态。

经过点击鼠标右键，闭合多边形或成功线段。

要删除形态，再次点击鼠标右键。

◎关于多边形形式，顶点间的距离被列出，当多边形闭合时，周长和总面积被记载。

◎关于折线形式，列出了顶点间的距离，当折线成功时，给出总距离。

◎关于矩形或椭圆形式，用鼠标左键点击并拖曳到所须要的形态大小。

如须要绘制一个正方形或圆，在矩形或椭圆形式下，按住鼠标中键的同时，点击并拖曳。

3）选用“Units＞所需测量单位”。

若图像的像元大小没有保留在文件头中，当产生“Input Display Pixel Size”对话框时，选用除“Pixel”之外的任何单位，在“X Pixel Size”和“Y Pixel Size”文本框中，输入图像的像元大小（注：TM 像元尺寸为30m ×30m）。

4）保留测量消息。

在“Display Measurement Tool”对话框内经常使用“File”菜单（图3-3），选用“File>Save Points to ASCII”，在“OutputMeasurement Directory”对话框内，输入一个输入文件名。

点击【OK】按钮，把测量消息保留到一个ASCII文件，格局为txt。

六、实验报告

（1）简述实验环节。

（2）回答疑问：①经过本次实验观察到的理想，说明遥感数字图像模型的构成。

②口头“Coursor Location/Value”操作可以获取遥感数字图像的哪些技术参数？Location操作框中显示的数值单位是什么？③比拟用手工量度加公式（3-1）测算的城市面积和用ENVI影像测量工具（Measurement Tool）测算的城市面积，对两者测量的分歧性和差异启动便捷剖析。

④可以经过对TM 影像像元坐标的量度，就能粗略计算出地物间的距离和面积，而不用依赖于天文坐标，为什么？

实验报告格局见附录一。

matlab图像解决教程

楼主你太狠了，5分要他人做这么多！1.图像的读入、显示及消息查问：(1)I=imread () %图像读入imshow(I)%图像显示(2)inf=imfinfo() % 图像消息查问2.图像的罕用解决语句：(1) X=rgb2gray(I) ; imshow(X) %黑白图像转灰度图像 (2)X2=grayslice(I,64) ; imshow(X2,hot(64)) %将灰度图像转为索引色图像(3) X3=im2bw(X1) ; imshow(X3) %将图像转二值图像3.图像滤波：clear allg0=imread()g0 = g0(:,:,2);%三维转二维figure(1);imshow(g0) ;title(原图) %如图 5g1=imnoise(g0,salt & pepper,0.2)g1=im2double(g1);figure(2);imshow(g1);title(参与椒盐噪声)%如图 6h1=fspecial(gaussian,4,0.3)g2=filter2(h1,g1,same)figure(3);imshow(g2);title(启动高斯滤波)%如图 7h2=fspecial(sobel)g3=filter2(h2,g1,same)figure(4);imshow(g3);title(启动sobel滤波)%如图 8h3=fspecial(prewitt)g4=filter2(h3,g1,same)figure(5);imshow(g4);title(启动prewitt滤波)%如图 9h4=fspecial(laplacian,0.5);g5=filter2(h4,g1,same);figure(6);imshow(g5);title(启动拉普拉斯滤波);%如图 10h5=fspecial(log,4,0.3);g6=filter2(h5,g1,same);figure(7);imshow(g6);title(启动高斯拉普拉斯滤波);%如图 11h6=fspecial(average);g7=filter2(h6,g1,same);figure(8);imshow(g7);title(启动均值滤波);%如图 12h7=fspecial(unsharp,0.3);g8=filter2(h7,g1,same);figure(9);imshow(g8);title(启动含糊滤波);%如图 13h8=[0 -1 0;-1 5 -1;0 -1 0];g9=filter2(h8,g1,same);figure(10);imshow(g9);title(启动高通高斯滤波);%如图 14h9=g1;g10=medfilt2(h9);figure(11);imshow(g10);title(启动中值滤波);%如图 15