遥感应用模型遥感反演土地变化检测

1、遥感觉用模型遥感反演土地变化检测遥感觉用模型遥感反演土地变化检测12/12遥感觉用模型遥感反演土地变化检测遥感觉用模型实习报告学院：班级：学号：姓名：日期：指导老师：Part1大冶研究区土地使用情况解析督查分类最大似然法2009年影像分类结果：分了4类：建筑，水体，裸地，植被2010年影像分类结果：栅格转矢量变化检测-叠加解析（1）先提取2009年和2010年分类图的建筑用地：2009年2010年（2）变化检测（叠加解析）Chang=2010-2009空间盘问变化图层Chang与建设用地红线进行盘问提取结果收效图：批而未用:用而未尽：正规使用：越界开发：未批先用将变化矢量图与2010年影像叠加

2、显示：Part2：遥感反演与建模数据预办理安装环境卫星数据办理补丁将ENVI_HJ1A1B_Tools.sav补丁放在ITTIDLIDL80productsenvi48save_add目录下。数据读取和定标选择主菜单-File-OpenExternalFile-HJ-1A/1BCCD，Inputpath选择环境卫星数据文件夹，点Tools，打开环境卫星数据办理补丁后，Search，设置输出路径，勾选“Calibration”“LayerStacking”，点Apply，以以下图：工程区裁剪由于整景影像范围太大了，进行几何校正从前，裁剪出我们需要的太湖及其周边地域。打开上一步办理好的数据：HJ1

3、B-CCD1-Calbrated_LayerStacking.img完成太湖及其周边地域的裁剪主菜单-File-SaveFileAs-ENVIStandard，弹出NewFileBuilder面板，单击ImportFile，弹出CreateNewFileInputFile面板，选中SelectInputFile面板中的数据，单击SpatialSubset，弹出SelectSpatialSubset面板，单击image弹出SubsetbyImage对话框，在其中裁剪出太湖及其周边地域，在几个对话框中单击ok，在NewFileBuilder面板中，单击choose，设置保存路径，输出文件名HJ1B

4、-CCD1-Cal-sub.img4.几何校正(1)打开基准影像TM_baseimage.img，选择Map-Registration-AutomaticRegistration：ImagetoImage，选择基准影像TM_baseimage.img的波段4作为般配波段选择被配准影像，选择band4作为般配波段，在提示可否手动选择同名点时，选择否后，弹出AutomaticRegistrationParameters面板(3)在AutomaticRegistrationParameters面板中，设置以下图参数单击ok执行基于像元值自动搜寻同名点。对于RMS高的点直接删除，水域中的点由于没有固定

5、参照物都删掉，湖岸线的点RMS高的点可以直接删掉，也许在两个影像的ZOOM窗口大将十字光标定位到正确位置，再点击Update进行微调，在缺少控制点的地方手动增加一些控制点。调整控制点，直到总的RMS小于1时，完成控制点的选择，点击GroundControlPointsSelection上的File-SaveGCPstoASCII，保存控制点。(4)在GroundControlPointsSelection面板上，选择Options-WarpFile（AsImageMap）选择校正文件。在校正参数面板中，投影参数默认。在X和Y的像元大小输入30，按回车，图像输出大小自动更正。选择输出路径和文件名

6、，点击ok进行ImageRegistration。(5)打开基准影像和校正影像，在显示校正后影像的窗口中，右键选择Link命令，选择显示图像的两个窗口，打开十字光标查察校正结果。Geographic5.大气校正环境小卫星供应了波谱响应函数，以文本形式供应第一列表示波长（nm），后边4列分别表示4个波段对应波长的波谱响应函值。为了做大气校正，需要制作波谱曲线来描述波谱响应函数，用于大气校正（波谱响应函数与CCD型号相对应）（1）制作波谱响应函数主菜单Window-StartNewPlotWindow，打开ENVIPlotWindow面板，选择File-InputData-ASCII，导入681_

7、HJ1BCCD1文本文件，波长单位选择Nanometers，单击ok。在绘制窗口出现了4条曲线，选择b1，b2，b3，b4，便于差异Edit-DataParameters，编写每条线的名称为选择LibraryFile-SavePlotAs-SpectralLibrary窗口中，单击SelectAllItems，单击，在ok。OutputPlotstoSpectral在OutputSpectralLibrary面板中，选择输出路径和文件名“HJ1BCCD1光谱响应”，单击ok，将波谱曲线保存为波谱库文件。2）FLASSH大气校正在选择主菜单BasicTools-ConvertData（BSQ、B

8、IL、BIPConvertFileParameters中选择BIL,选择ConvertInPlace）选择校正后的影像，：yes，单击ok。主菜单Spectral-FLAASH，打开FLAASH大气校正模块，点击InputRadianceImage，选择BIL格式的，在弹出的RadianceScaleFactors面板中，选择Usesinglescalefactorforallbands，在此Singlescalefactor选择默认：设置输出路径及文件名，传感器基本信息设置数据可以从数据头文件10，单击ok。.XML中读取大气模型选MLS，气溶胶模型选Rural，气溶胶反演方法选None，能

9、见度给40km单击MultispectralSetting按钮，在FilterFunctionFile谱响应曲线“HJ1BCCD1光谱响应”，单击ok。单击AdvancedSettingsa按钮，TileSize设置为100，点击导入从前做好的光ok。在大气校正模块中，单击apply，进行大气校正。校正完成后检查校正结果，分别加载校正前后的影像，将两幅影像进行地理链接，搬动到植被地域，在影像上右键选择ZProfile打开光谱曲线窗口，显示两幅图像同一地址的光谱曲线。左边为大气校正前，右边为大气校正后：地域裁剪打开大气校正后的图像dqjz.img，在image窗口中选择Overlay-Vecto

10、rs，打开VectorParameters窗口，选择File-OpenVectorFile，打开taihu_based.shp。选择正确的投影种类，弹出AvailableVectorsList面板，选择该矢量文件，点击LoadSelected，选择显示图像的Display，点击ok，将矢量叠加在影像上。在AvailableVectorsList面板，选择File-ExportLayerstoROI，在SelectDataFiletoAssociatewithnewROI面板中选择影像，点击ok，在弹出的对话框中选择ConvertallrecordsofanEVFlayertooneROI，点击

11、ok，将矢量转变成ROI；在图像窗口，选择Overlay-RegionofInterest，打开ROI面板，选择File-SubsetDataviaROIs，在SelectInputFiletoSubsetviaROI中，选择影像，单击ok；在SpatialSubsetviaROIParameters中选择太湖ROI，MaskpixelsoutsideofROI选择Yes；设置输出路径及文件名dqjz-taihu.img，单击ok。反演建模（1）获取采样星上数据使用BasicTool-Bandmath，在Enteranexpression下面输入表达式：float(b4)/b3，单击Addto

12、List，单击ok，在VariablestoBandsPairings面板中选择第3、4波段对应b3、b4，设置输出路径和文件名“b4除b3”点击ok，计算获取比值图像；在display中显示比值图像，选择overlay-RegionOfInterest，打开ROITool在ROITool中，选择ROI_Type-InputPointsfromASCII，选择文本格式的反演点.txt，x选择经度3，y选纬度2，Thesepointcomprise：IndividualPoints设置投影信息后，点击ok，将实地检查的点位信息加载到图像中。在ROITOOL中，选择File-OutputROIst

13、oASCII。选择b4/b3的图像，在OutputROIstoASCIIParameters面板中选择ROI点，单击EditOutputASCIIForm，在输出内容设置面板中，选择ID、经纬度（GeoLocation）、波段像元值（BandValues）模型参数反演：叶绿素反演:线性模型2.指数模型3.对数模型Part3QiuckBird数据交融后分类提取影像交融（1）多光谱影像重采样在MapGis中查察影像信息后，可知多光谱影像分辨率为2.4m，全色影像分辨率为0.6m，在进行影像交融前先对多光谱影像重采样，使其分辨率与全色影像的分辨率相同。重采样前：重采样后：（2）裁剪由于整幅影像太大，

14、办理起来计算机要赞同很长时间，于是先对影像裁剪出光谷周边的部分地域。记录下裁剪后的多光谱影像的行列值，对高分辨率的全色影像做相同地域的裁剪：（3）影像交融对影像进行加权交融法交融交融后影像：（4）影像交融精度议论督查分类(最大似然法)原始影像:采用AOI区：分类结果：精度议论:专题图:2.非督查分类(IsoData)分类结果：精度议论：面积统计：面向对象分类分类结果：精度议论：Part4总结本次实习主要分为三大部分，土地利用变化检测，叶绿素反演，分类。在这三块中主要用的的软件有MapGIS和ENVI4.5。在土地变化检测中，先对时相一和时相二的影像进行督查分类，提取出建筑用地，再进行叠加解析，得出建筑用地变化部分。与所获取的建筑用地红线进行比较，获取批而未用，用而未尽，越界开发，未批先用等违规现象。在这里面从影像里面提取建设用地的精度决定了最后结果的精度。主要锻炼了我们做督查分类，地图编写器的使用。实践了使用遥感的方式做土地利用变化的方法。在叶绿素的反演建模中，依照实习指导书上的操作，就可以进行叶绿素的反演，其实在这个过程中，操作不难，主若是经过这个过程，即可以熟悉反演这个过程。在第三部分的分类，主若是面向对象的分类。面向对象的分类，就主要分两大部分，发现对