null栅格数据分析栅格数据分析邱岭
Email: gisql@sina.com
栅格数据分析Page *栅格数据分析 五、综合分析应用 四、栅格和矢量转换 三、栅格计算器使用本章
内容
财务内部控制制度的内容财务内部控制制度的内容人员招聘与配置的内容项目成本控制的内容消防安全演练内容
二、栅格的重分类一、栅格数据概念1、栅格数据概念Page *1、栅格数据概念栅格数据结构又称为网格数据,或则栅格数据
特点: 结构简单:只有行列,对应像元的值,属性明显,表达直观。
不连续性(离散):数据之间有间隔,图形无法实现无级放大。栅格数据分类Page *栅格数据分类 栅格数据有两种:一种是由单一的层构成的称为栅格,另一种是由多层复合而成的称为影像(比如,多波段的卫星照片)
而卫星影像等数据包含了多个波段的信息,是分层存储的,每个波段的数据存储在一层中,每层中记录了某一波段范围内电磁光谱的反射系数。通过给每一层赋不同的颜色(R,G,B),可以生成合成影像,增强某方面的信息,以利于评价土地覆盖、植被密度等因素。
简单说是栅格是单波段,
影像是多个波段。栅格表现Page *栅格表现 栅格在一种特殊的影像格式,一般影像格式有JPG,BMP,TIF,PNG等,在ArcGIS中一般表现为Grid,.img,.TIF,具体的应用为DEM,DOM等,最直观表现栅格为单个波段RasterBand,影像分析也是基于单个波段RasterBand,所以在ArcGIS影像分析也成栅格分析。
任何栅格分析(处理)都是信息有损的处理,原始数据一定要保存(保留)。栅格数据介绍-CELLSIZEPage *栅格数据介绍-CELLSIZE定义格格单元的大小时,我们需要平衡信息的精确性和数据量之间的矛盾。
栅格单元代表的尺度越小CellSIZE,表达的信息就越精确。
栅格单元代表的尺度越大,存储数据所需要的空间就更少,同时,表达的信息也就不精确。
在ArcGIS中改变CellSIZE,使用Resample 重采样
Resample:可以处理单波段,也可以处理多波段数据:12\resample改变分辨率Page *改变分辨率多段取一个波段Page *多段取一个波段按掩膜提取可以导出每个波段,不要输出扩展名数据:11_3d\dem3d\g14.jpg几个单波段合成波段Page *几个单波段合成波段Composite Bands波段合成使用上个幻灯片的数据栅格数据集像元的位深度容量Page *栅格数据集像元的位深度容量 像元的位深度(像素深度)决定着特定栅格文件可以存储的值的范围,该范围可根据公式 2n 计算得出(其中,n 表示位深度)。例如,一个 8 位的栅格可以具有 256 个不同的值(范围从 0 至 255)。2、栅格的重分类Page *2、栅格的重分类重分类(Reclassify)
将栅格图层的数值进行重新分类组织或者重新解释。
重分类的关键是确定原数据到新数据之间的对应关系。
重分类只能从(详)细到粗(略),不能相反、可逆操作。
Reclassify对多波段影像,按第一个波段处理2、栅格的重分类Page *2、栅格的重分类选择Raster,和对应字段(Values)
值的范围是含下不含上
右图是1900 – 2000
含义为value>=1900
<2000,“-”前后都有一个空格数据为:12\resample\tingrid2、栅格的重分类NoDataPage *2、栅格的重分类NoData NoData是空,什么也没有(NULL),空不等于0,设置为nodata,输出为空白(默认是透明的),从表面上颜色和背景色一样,使用信息查询Info查询对应位置数据返回的是Nodata两个不同的概念Page *两个不同的概念重分类Reclassify:是对单个波段,改变值的分布
重采样 Resample,是改变影像分辨率(每个像素点代表的矢量大小=cellsize),可以用于多波段3、栅格计算器使用(中级)Page *3、栅格计算器使用(中级) 栅格计算器Raster Calculator是栅格数据空间分析中数据处理和分析中最为常用的方法,应用非常广泛,能够解决各种类型的问
题
快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题
,尤其重要的是:建立复杂的应用数学模型的基本模块。
ArcGIS 提供了非常友好的图形化栅格计算器,利用栅格计算器,不仅可以方便的完成基于数学运算符的栅格运算,以及基于数学函数的栅格运算,而且它还支持直接调用ArcGIS 自带的栅格数据空间分析函数,并且可以方便的实现多条语句的同时输入和运行。3、栅格计算器Raster Calculator使用Page *3、栅格计算器Raster Calculator使用3、栅格计算器Raster Calculator使用Page *3、栅格计算器Raster Calculator使用注意事项:
图层名称前后加 双引号,9.3是中括号
运算符前后加空格 9.3时必须
大小写含义是不一样的,arcgis10新
通过name=表达式 可以修改输出栅格图层名称 9.3时使用
一次可以输出多个结果,以回车作为标志,上行的输出结果,可以直接后面的使用 9.3时使用
一行内容很长,使用“~” 放在上一行行尾作为续行标志 9.3时使用函数说明-算术函数Page *函数说明-算术函数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),反之返回0函数说明Page *函数说明三角函数包括:Sin(正弦函数)、Cos(余弦函数)、Tan(正切函数)、Asin(反正弦函数)、Acos(反余弦函数)、Atan(反正切函数)。
对数函数(Logarithms):对数函数可对输入的格网数字做对数或指数的运算。指数部份包括:Exp (底数e)、Exp10 (底数10)、Exp2 (底数2)三种;对数部份包括:Log (自然对数)、Log10 (底数10)、log2 (底数2)等三种。
幂函数(Powers):幂函数可对输入的格网数字进行幂函数运算。幂函数包括三种:Sqrt (平方根)、Sqr (平方)、Pow (幂)。函数说明-其他重要函数(高级)Page *函数说明-其他重要函数(高级)Con 语法如下
Con(
, , _ {, }, _ {, }, _ {false_expression})
如Con("tingrid" >= 2000,3,"tingrid">= 1900,2,1)
含义
if "tingrid" >= 2000 then 返回值3
else if "tingrid" >= 1900 then 返回值2
else 返回值1
使用Con可以实现栅格重分类,Con比较灵活,栅格重分类比较直观,简单。函数说明-其他重要函数Page *函数说明-其他重要函数2. SetNull 满足条件设置为空
SetNull("tingrid">= 2000, "tingrid") 含义
If tingrid" >= 2000 then 返回空Nodata
Else 返回自己"tingrid",返回小于2000的数据
SetNull ("tingrid" - 2000, "tingrid"),不等于2000.设置为空,只取值为2000的数据。
比较两幅影像的不同Page *比较两幅影像的不同SetNull("tingrid" == "tingrid2","tingrid") 比较("tingrid", "tingrid2"相同取空,不同的取"tingrid"
2. 直接相减"tingrid" - "tingrid2"空和0转换Page *空和0转换把空变成0
Con(IsNull("tingrid"),0,"tingrid")
把0变成空
Con("tingrid2" <> 0,"tingrid2" )
把某个值(2000)转换为空是类似的
Con("tingrid2" <> 2000,"tingrid2" )
4.矢量和栅格相互转换Page *4.矢量和栅格相互转换矢量转栅格:几种方法
是等值线,通过建立tin生成栅格,见第8讲
离散点,生成栅格见第8讲
分析工具中要素转栅格FeatureToRaster
转成图片格式等
①-③都是特定数据格式,不通用
栅格转矢量
4.1.矢量转图片栅格Page *4.1.矢量转图片栅格在Arcmap文件,加载地图,使用(文件File)中输出地图,格式有PNG,emf,eps,pdf,jpg,bmp,svg,gif,tif等,自己设置大小和格式,写坐标文件(write world file)是自动配准,可以调整分辨率,pdf,svg,emf默认是300,其他为96,数字越高,输出的文件越大,清晰度越高,图片大小地图窗口的尺寸*分辨率,该工具只是针对当前地图窗口的内容
另一种土方法:抓屏,printscrn,抓屏幕
alt+ printscrn抓当前窗口
数据:打开\12\pic\pdf.mxd,看导出地图。exe界面上有汉字输出pdf,出现乱码问题Page *界面上有汉字输出pdf,出现乱码问题 界面上有汉字输出pdf,出现乱码问题,
下面两个都选中数据:打开\12\pic\pdf.mxd5.1输出指定比例尺图片Page *5.1输出指定比例尺图片 在Arcmap中查看全图,进入布局窗口,在数据框中设置固定比例,如1:10000后确定,调整布局的大小,最后使用导出地图当前窗口在布局数据:打开\12\pic\pdf.mxd,看指定比例输出.exe4.2 已有栅格转矢量Page *4.2 已有栅格转矢量有两种
使用Arctoolbox中要素转栅格,
Arcscan
说明1、主要特定的栅格数据,2是真正意义上的栅格转矢量,1在数据分析中有很重要的意义,操作方法类似。
4.2 已有栅格转矢量栅格转面
raster to FeaturePage *4.2 已有栅格转矢量栅格转面
raster to Feature1、tingrid转面
直接也无法转换,先重分类reclassify,分级大小,决定了生成矢量和原始栅格影像准确程度,分级越多,越准确,生成的矢量记录越多,每一个矢量的格子就越小。
使用数据:12\fx\tingrid,栅格转面.exe5.影像(栅格)分析(高级)Page *5.影像(栅格)分析(高级)使用影像分析以下几个应用
分类统计区段面积
统计25度以上耕地的面积
计算图斑图的坡度,类似求一个地区平均海拔
DEM中某个点修改
影像分析,是基于一个波段的,对于多段,只能提取其中一个波段。5.1影像面积统计Page *5.1影像面积统计统计高程1900以下,1900-2000,和2000以上的面积
使用数据:12\fx\tingrid,操作看影像面积统计1.exe5.1影像面积统计Page *5.1影像面积统计用重分类,以表格显示分区统计。使用数据:12\fx\tingrid,操作看影像面积统计2.exe5.2统计25度以上耕地的面积Page *5.2统计25度以上耕地的面积数据:\12\fx\, 5.3计算图斑图的坡度Page *5.3计算图斑图的坡度 已有坡度图slope和图斑(nc_parcl)图计算,计算图斑图的坡度级别。
增加唯一值字段
以表格显示分区统计,
输出属性表,获得每个图斑的坡度。
属性关联
更新坡度
使用数据:12\fx\slope,图斑为12\fx\slope\Data下nc_parcel
计算图斑图的坡度.exe属性关联Page *属性关联设置对应的字段求一个地方的平均海拔Page *求一个地方的平均海拔已有dem和行政区划图,
求各行政区划的平均海拔
Dem使用:\12\fx\tingrid
行政区划:\12\fx\data\NCV_CANTON
操作看:\12\fx\平均海拔.exe5.4 DEM中某个点修改Page *5.4 DEM中某个点修改数据:\12\resample\dem.mxd,
要点:Con(IsNull("pp"),"tingrid","PP")
具体操作:\12\dem坏点.exe数据模型Page *数据模型本章掌握内容Page *本章掌握内容重采样, 重分类区别
栅格计算器使用方法
要素转栅格
栅格转要素
分区表格统计