栅格数据分析

栅格数据分析1/41栅格数据分析五、综合分析应用四、栅格和矢量转换三、栅格计算器使用本章内容二、栅格重分类一、栅格数据概念22/411、栅格数据概念栅格数据构造又称为网格数据，或则栅格数据特点：

构造简单：只有行列，对应像元值，属性显著，体现直观。

不连续性（离散）：数据之间有间隔，图形无法实现无级放大。33/41栅格数据分类栅格数据有两种：一种是由单一层组成称为栅格，另一种是由多层复合而成称为影像（例如，多波段卫星照片）而卫星影像等数据包括了多种波段信息，是分层存放，每个波段数据存放在一层中，每层中统计了某一波段范围内电磁光谱反射系数。通过给每一层赋不一样颜色（R,G,B），能够生成合成影像，增强某方面信息，以利于评价土地覆盖、植被密度等原因。

简单说是栅格是单波段，影像是多种波段。44/41栅格体现栅格在一种特殊影像格式，一般影像格式有JPG，BMP，TIF，PNG等，在ArcGIS中一般体现为Grid,.img,.TIF，详细应用为DEM，DOM等，最直观体现栅格为单个波段RasterBand，影像分析也是基于单个波段RasterBand，因此在ArcGIS影像分析也成栅格分析。任何栅格分析（处理）都是信息有损处理，原始数据一定要保存（保存）。55/41栅格数据介绍-CELLSIZE定义格格单元大小时，我们需要平衡信息精确性和数据量之间矛盾。栅格单元代表尺度越小CellSIZE，体现信息就越精确。栅格单元代表尺度越大，存放数据所需要空间就更少，同步，体现信息也就不精确。在ArcGIS中变化CellSIZE，使用Resample重采样Resample：能够处理单波段，也能够处理多波段数据：12\resample66/41变化辨别率77/41多段取一种波段按掩膜提取能够导出每个波段,不要输出扩展名数据:11_3d\dem3d\g14.jpg88/41几个单波段合成波段CompositeBands波段合成使用上个幻灯片数据99/41栅格数据集像元位深度容量

像元位深度（像素深度）决定着特定栅格文献能够存放值范围，该范围可根据公式2n计算得出（其中，n表达位深度）。例如，一种8位栅格能够具有256个不一样值（范围从0至255）。1010/41栅格数据集格式技术规范-help1111/412、栅格重分类重分类（Reclassify）将栅格图层数值进行重新分类组织或者重新解释。重分类关键是确定原数据到新数据之间对应关系。重分类只能从（详）细到粗（略），不能相反、可逆操作。Reclassify对多波段影像，按第一种波段处理1212/412、栅格重分类选择Raster，和对应字段（Values）值范围是含下不含上右图是1900–2023含义为value>=1900<2023,“-”前后都有一种空格数据为:12\resample\tingrid1313/412、栅格重分类NoDataNoData是空，什么也没有（NULL），空不等于0，设置为nodata，输出为空白（默认是透明），从表面上颜色和背景色同样，使用信息查询Info查询对应位置数据返回是Nodata1414/41两个不一样概念重分类Reclassify：是对单个波段，变化值分布重采样Resample，是变化影像辨别率（每个像素点代表矢量大小=cellsize），能够用于多波段1515/413、栅格计算器使用(中级）栅格计算器RasterCalculator是栅格数据空间分析中数据处理和分析中最为常用办法，应用非常广泛，能够处理多种类型问题，尤其主要是：建立复杂应用数学模型基本模块。ArcGIS提供了非常友好图形化栅格计算器，利用栅格计算器，不但能够方便完成基于数学运算符栅格运算，以及基于数学函数栅格运算，并且它还支持直接调用ArcGIS自带栅格数据空间分析函数，并且能够方便实现多条语句同步输入和运行。1616/413、栅格计算器RasterCalculator使用1717/413、栅格计算器RasterCalculator使用注意事项：图层名称前后加双引号,9.3是中括号运算符前后加空格9.3时必须大小写含义是不同样，arcgis10新通过name=体现式能够修改输出栅格图层名称9.3时使用一次能够输出多种成果，以回车作为标志，上行输出成果，能够直接背面使用9.3时使用一行内容很长，使用“~”放在上一行行尾作为续行标志9.3时使用1818/41函数说明-算术函数1、ceil函数返回参数大于num最小整数（返回值大于等于原始数据）y=6.04;x=ceil(y);x为7.0.y=5;x=ceil(y);x为5.0.2、floor()小于或等于传入参数整数（返回值大于等于原始数据）。y=6.04;x=floor(y);x为6.0.y=5;x=floor(y);x为5.0.3、float（转换为浮点数）4、Int取整数5、Abs取绝对值6、IsNull（raster）如果为空，条件成立：返回值1（表示True），反之返回01919/41函数说明三角函数包括：Sin（正弦函数）、Cos（余弦函数）、Tan（正切函数）、Asin（反正弦函数）、Acos（反余弦函数）、Atan（反正切函数）。对数函数（Logarithms）：对数函数可对输入格网数字做对数或指数运算。指数部份包括：Exp(底数e)、Exp10(底数10)、Exp2(底数2)三种；对数部份包括：Log(自然对数)、Log10(底数10)、log2(底数2)等三种。

幂函数（Powers）：幂函数可对输入格网数字进行幂函数运算。幂函数包括三种：Sqrt(平方根)、Sqr(平方)、Pow(幂)。2020/41函数说明-其他主要函数（高级）Con语法如下Con(<Condition>,<true_expression>,_{<Condition>,<true_expression>},_{<Condition>,<true_expression>},_{false_expression})如Con("tingrid">=2023,3,"tingrid">=1900,2,1)含义if"tingrid">=2023then返回值3elseif"tingrid">=1900then返回值2else返回值1使用Con能够实现栅格重分类，Con比较灵活，栅格重分类比较直观，简单。2121/41函数说明-其他主要函数2.SetNull满足条件设置为空SetNull("tingrid">=2023,"tingrid")含义Iftingrid">=2023then返回空NodataElse返回自己"tingrid",返回不大于2023数据SetNull

("tingrid"-2023,"tingrid")，不等于2023.设置为空，只取值为2023数据。2222/41比较两幅影像不一样SetNull("tingrid"=="tingrid2"，"tingrid")比较("tingrid"，"tingrid2"相同取空，不一样取"tingrid"2.直接相减"tingrid"-"tingrid2"2323/41空和0转换把空变成0Con(IsNull("tingrid"),0,"tingrid")把0变成空Con("tingrid2"<>0,"tingrid2")把某个值（2023）转换为空是类似Con("tingrid2"<>2023,"tingrid2")2424/414.矢量和栅格互相转换矢量转栅格：几个办法是等值线，通过建立tin生成栅格，见第8讲离散点，生成栅格见第8讲分析工具中要素转栅格FeatureToRaster转成图片格式等①-③都是特定数据格式，不通用栅格转矢量2525/414.1.矢量转图片栅格在Arcmap文献，加载地图，使用（文献File）中输出地图，格式有PNG，emf，eps，pdf，jpg，bmp，svg，gif，tif等，自己设置大小和格式，写坐标文献（writeworldfile）是自动配准，能够调整辨别率，pdf，svg，emf默认是300，其他为96，数字越高，输出文献越大，清楚度越高，图片大小地图窗口尺寸*辨别率，该工具只是针对目前地图窗口内容另一种土办法：抓屏，printscrn，抓屏幕alt+printscrn抓目前窗口

数据：打开\12\pic\pdf.mxd，看导出地图。exe2626/41界面上有汉字输出pdf，出现乱码问题界面上有汉字输出pdf，出现乱码问题，下面两个都选中数据：打开\12\pic\pdf.mxd2727/415.1输出指定百分比尺图片

在Arcmap中查看全图，进入布局窗口，在数据框中设置固定百分比，如1：10000后确定，调整布局大小，最后使用导出地图目前窗口在布局数据：打开\12\pic\pdf.mxd，看指定百分比输出.exe2828/414.2已有栅格转矢量有两种使用Arctoolbox中要素转栅格，Arcscan说明1、主要特定栅格数据，2是真正意义上栅格转矢量，1在数据分析中有很主要意义，操作办法类似。2929/414.2已有栅格转矢量栅格转面

rastertoFeature1、tingrid转面直接也无法转换，先重分类reclassify，分级大小，决定了生成矢量和原始栅格影像精确程度，分级越多，越精确，生成矢量统计越多，每一种矢量格子就越小。

使用数据:12\fx\tingrid，栅格转面.exe3030/415.影像（栅格）分析（高级）使用影像分析下列几个应用分类统计区段面积统计25度以上耕地面积计算图斑图坡度,类似求一种地域平均海拔DEM中某个点修改影像分析，是基于一种波段，对于多段，只能提取其中一种波段。3131/415.1影像面积统计统计高程1900下列，1900-2023，和2023以上面积使用数据:12\fx\tingrid,操作看影像面积统计1.exe3232/415.1影像面积统计用重分类,以表格显示分区统计。使用数据:12\fx\tingrid,操作看影像面积统计2.exe3333/415.2统计25度以上耕地面积数据:\12\fx\,3434/415.3计算图斑图坡度

已有坡度图slope和图斑(nc_parcl)图计算，计算图斑图坡度级别。增加唯一值字段以表格显示分区统计，输出属性表，取得每个图斑坡度。属性关联更新坡度使用数据:12\fx\slope,图斑为12\fx\slope\Data下nc_parcel计算图斑图坡度.exe3535/41属性关联设置对应字段3636/41求一种地方平均海拔已有dem和行政区划图，求各行政区划平均海拔Dem使用:\12\fx\tingrid行政区划:\12\fx\data\NCV_