我们付出一些成本,时间的或者其他,最终总能收获一些什么。
【FME-HOW-TO系列】13 通过重新采样修改栅格像元大小
除了使用RasterResampler转换器进行重采样的操作外,还需要了解不同的插值方法,各方法大概的不同。
可以参考ArcGIS栅格重采样方法介绍
Nearest Neighbor(最邻近法):最邻近法是把原始图像中距离最近的像元值填充到新图像中;该方法的优点是不会改变原始栅格值,而且处理速度快,但是该方法会有半个像元大小的位移。适用于表示分类或某种专题的离散数据,如土地利用、植被类型等
Bilinear(双线性内插法):双线性插值法是基于三次线性插值的方法,将输入栅格数据集中的4个最邻近像元(4邻域)的像元值进行加权平均计算出新的像元值,并将其赋予输出栅格数据集的相应像元。其中,权值是由4邻域中每个像元的中心与内插点之间的距离决定的。该方法的重采样结果会比最邻近法更平滑,但会改变原来的栅格值。适用于表示某种现象分布、地形表面的连续数据,如 DEM、气温或降雨量分布、坡度等,这些数据本来就是通过采样点内插得到的连续表面。
Bucubic(三次卷积插值法):三次卷积内插法是基于五次多项式插值的方法,将输入栅格数据集中的16个最邻近像元(16邻域)的像元值进行加权平均计算出新的像元值,并将其赋予输出栅格数据集的相应像元。其中,权值是由16邻域中每个像元的中心点与内插点之间的距离决定的。
三次卷积内插法通过增加邻近点来获取最佳插值函数,可以进一步提高内插精度,算法较为复杂,计算量大,处理时间较长。由于该方法使用16邻域进行加权计算,处理结果会更加清晰,栅格数据的边界会有锐化的效果。该方法同样会改变原来的栅格值,且有可能会超出输入栅格的值域范围。适用于航片和遥感影像的重采样。
Average 4,Average 16:平均4和平均16是对于像DEM这样的数字栅格有用的插值类型,并且产生比双线性和最近邻的更高质量的图像。(大致原理没有搜)
还可以对重采样相近的名词进行知识拓展,重采样包括上采样和下采样。
回到FME中RasterResampler转换器的使用,主要任务是Average4的采样方式,将原始的DEM像元数减少一半。
RasterResampler转换器说明:使用所需的维数重采样一个输入的栅格,期望的单元格大小是地面单位或是单元格大小的百分比。
添加写模块Esri ASCII Grid,将“TAB文件生成”设置为“是”
TAB文件用于存储控制点,坐标系和用户属性
【FME-HOW-TO系列】14 按地理要素更改栅格值
读模块添加一个矢量面数据和一个GeoTiff栅格数据,传入到MapnikRasterizer转换器,修改相关参数,前面对矢量数据进行栅格化时,已使用过该转换器。
平滑Gamma值为1使用模糊逻辑,会使得边界模糊,此值设为0
存在的问题是,FME中的显示是有四个波段红、绿、蓝、透明度,ArcMap中只有红、绿、蓝三个波段,原来的geoTiff(0,0,0)与周围的nodata一致的情况
【FME-HOW-TO系列】15 减少不必要的坐标点
主要使用Generalizer转换器,通过设置参数的方式,对图形进行简化
Generalizer转换器说明:用特定的算法变换和测量要素。
这里有四类算法:
【FME-HOW-TO系列】16 来自点数据的TIN(不规则三角网)
主要使用TINGenerator转换器进行三角网的创建,使用Scaler转换器设置Z比例因子为0.1
TINGenerator转换器说明:基于输入的点和断裂线构建一个Delaunay三角网。表面模型可以用多个表示形式输出:不规则三角网(TIN),TIN顶点,TIN边和三角网。
Scaler:放大或者缩小要素。
读模块读取:ASPRS激光雷达数据交换格式(LAS)格式
写模块为shapefile
【FME-HOW-TO系列】17 通过克里金法将点插值转换为栅格
使用到的转换器
Reprojector转换器说明:把要素从一个坐标系统重投影到另一个坐标系统。
CoordinateExtractor转换器说明:按指定的索引检索坐标点XYZ值存储到属性字段中。
后续使用的RCaller需要安装geoR和spR模块,该模块没有安装,后续需要了解一下
TempPathnameCreator转换器:保留并返回在转换完成后将被删除的临时文件或目录路径。