span8
span4
简介:RCaller入门|上一篇:RCALLER:树高和树宽是否相关?
本文也是教程的一部分在FME里我该怎么做|下一步:从点生成等高线数据
的运输机transformer允许进行更多的统计分析统计计算器变压器。这个例子展示了如何使用RCaller从点创建栅格。克里格法是一种插值方法,用于从细胞信号点数据创建一个估计的表面。
CSV格式的手机信号点数据。
使用数据检查器查看的源csv数据
下载克里金FME2016.1.FMWT是完整的工作区。如果您想自己创建工作空间,请下载数据点.csv,然后按照下面的步骤创建它。
请注意,工作区要求用户安装“geor”和“sp”r模块。
1.读取源数据
使用CSV阅读器读取数据点。CSV。
2.责备
使用投影仪将数据从lat/long(度)重新投影到UTM83-10(米)
三。提取坐标值
CoordinateExtractor提取点的x和y坐标值,并将它们添加为名为easting (x)和northing (y)的属性。东向和北向将被用于在RCaller中创建地图和kriging。
四。为光栅输出创建临时文件
其中一个特性用于使用TempPathnameCreator生成临时文件。来自RCaller的光栅图像将被临时保存到这个位置。
5个。添加呼叫器
R脚本将点插入到光栅中,并在临时文件位置创建光栅png。对于R脚本,请指定以下内容:
库(sp)库(geoR) map = SpatialPointsDataFrame(data=CellSignals, coords=cbind(CellSignals$easting,CellSignals$northing)网格< -扩大。grid( seq(min(CellSignals$easting),max(CellSignals$easting),50),seq(min(CellSignals$northing),max(CellSignals$northing),50) ) ml=likfit(coords=coordinates(map), data=map$Quality, ini=c(10000,10000), nug=10) kc <- krige.conv(coords=coordinates(map), data=map$Quality, loc=pred.grid, krige=krige.control(obj.m=ml)) png(CellSignals$pathname[1]) image(kc, loc = pred.grid, col=gray(seq(1,0,l=30)), xlab="Coord X", ylab="Coord Y") dev.off() fmeOutput <- data.frame(pathname=CellSignals$pathname[1])
调用方参数
6.读取光栅PNG
FeatureReader读取RCaller中创建的光栅png,从而允许我们查看输出数据。
第七章。查看结果
在数据检查器中检查输出光栅。结果表明,RCaller可以用来从点数据生成栅格。请注意,此示例在统计上可能不准确。
数据检查器中显示的克里格输出
可选:直接从rcaller读取r光栅对象
在FME 2019.0 beta(从构建19074开始)中,您可以直接从RCaller输出一个R光栅对象,以便在FME中进一步使用。对于那些希望使用R包进行栅格处理,然后在FME工作流中继续工作的用户来说,这个场景非常理想。亚搏在线
关于这个过程的详细文档将在春季发布,但是现在,您可以找到一个显示该方法的附加模板。基本流程是:
所附模板kriging_fme2019_光栅_对象.fmwt包含此示例。
在早期版本的fme中,技术上也可以得到类似的结果,但您将无法从rcaller直接输出光栅对象。相反,您必须编写一个临时光栅(使用例如光栅包的writeRaster()函数)然后使用featurereader和文件路径属性将其读回工作区。
你好,欧文,
我添加了一个简短的部分并上传了一个模板(克里格-fme2019-光栅-对象.fmwt)显示如何输出地理参考光栅而不是打印。请注意,您需要FME 2019.0测试版才能轻松使用此功能。这在旧版本中是可能的,但它有点复杂。请看一下模板,看看它是否适合您的需要。我很高兴回答你对它使用的任何问题。
我们希望在spring中发布一个更全面、更出色的例子。
©2019安全亚搏在线软件公司法律