免费视频|新人指南|投诉删帖|广告合作|地信网APP下载

查看: 1789|回复: 8
收起左侧

【FME-HOW-TO系列】17 通过克里金法将点插值转换为栅格

[复制链接]

665

主题

2万

铜板

34

好友

传奇会员

Rank: 30Rank: 30Rank: 30Rank: 30Rank: 30Rank: 30Rank: 30Rank: 30

积分
20649

精华勋章爱心勋章地信专家组贡献勋章冰雪节勋章

发表于 2021-3-26 15:04 | 显示全部楼层 |阅读模式
介绍

RCaller转换器允许比StatisticsCalculator转换器多得多的统计分析。此示例显示如何使用RCaller从点创建栅格。Kriging是一种插值方法,用于根据单元信号点数据创建一个拟合曲面。

下载

datapoints.csv

kriging-fme20161.fmwt

来源数据

CSV格式的单元信号点数据。

使用Data Inspector查看的源CSV数据

说明

下载的kriging_FME2016.1.fmwt是完整的工作空间。如果您想自己创建工作空间,请下载DataPoints.csv,然后按照以下步骤创建它。

请注意,工作空间要求用户安装“geoR”和“sp”R模块。

1.读取源数据

使用CSV读模块读取DataPoints.csv。

2.重投影

使用Reprojector将数据从纬度/精度(度)重新投影到UTM83-10(米)

3.提取坐标值

CoordinateExtractor提取点的x和y坐标值,并将它们添加为名为easting(x)和northing(y)的属性。Easting和northing将用于在RCaller中创建地图和克里金。

4.为栅格输出创建临时文件

一个要素用于使用TempPathnameCreator生成临时文件。来自RCaller的栅格图像将临时保存到此位置。

5.添加RCaller

R脚本将点插入到栅格中,并在临时文件位置创建栅格png。对于R脚本,请指定以下内容:

·
library(sp)

·  library(geoR)

·   

·  map = SpatialPointsDataFrame(data=CellSignals, coords=cbind(CellSignals$easting,CellSignals$northing))

·  pred.grid <- expand.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])


RCaller参数

6.读取栅格png

FeatureReader读取在RCaller中创建的栅格png,它允许我们查看输出数据。

7.查看结果

在Data Inspector中检查输出栅格。结果表明RCaller可用于从点数据生成栅格。请注意,此示例可能在统计上不精确。

克里金输出显示在Data Inspector


本文转载来自CSDN作者:fmechina  版权归作者所有

原文链接:https://blog.csdn.net/fmechina/article/details/102458513?spm=1001.2014.3001.5502

10

主题

2万

铜板

34

好友

黄金会员

Rank: 23Rank: 23Rank: 23Rank: 23Rank: 23Rank: 23Rank: 23

积分
4395

精华勋章灌水勋章活跃勋章贡献勋章

发表于 2021-3-26 15:27 | 显示全部楼层
谢谢分享
回复

使用道具 举报

2

主题

337

铜板

1

好友

技术员

Rank: 3Rank: 3

积分
27
发表于 2021-3-27 19:53 | 显示全部楼层
:ai:ai:ai
回复

使用道具 举报

206

主题

1万

铜板

25

好友

传奇会员

Rank: 30Rank: 30Rank: 30Rank: 30Rank: 30Rank: 30Rank: 30Rank: 30

积分
38667

灌水勋章宣传勋章贡献勋章活跃勋章

QQ
发表于 2021-3-29 08:10 | 显示全部楼层
谢谢分享
回复

使用道具 举报

6

主题

1万

铜板

15

好友

传奇会员

Rank: 30Rank: 30Rank: 30Rank: 30Rank: 30Rank: 30Rank: 30Rank: 30

积分
20795
发表于 2021-3-29 08:31 | 显示全部楼层
谢谢分享谢谢分享
回复 支持 反对

使用道具 举报

0

主题

3万

铜板

9

好友

黄金会员

Rank: 23Rank: 23Rank: 23Rank: 23Rank: 23Rank: 23Rank: 23

积分
4457
发表于 2021-12-12 12:56 | 显示全部楼层
谢谢分享
回复

使用道具 举报

0

主题

1047

铜板

4

好友

助理工程师

Rank: 5Rank: 5

积分
399
发表于 2022-8-12 16:52 | 显示全部楼层
感谢分享
回复

使用道具 举报

2

主题

3万

铜板

6

好友

钻石会员

Rank: 26Rank: 26Rank: 26Rank: 26Rank: 26Rank: 26Rank: 26

积分
5665
发表于 2024-3-1 14:13 | 显示全部楼层
感谢楼主分享
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

在线客服
快速回复 返回顶部 返回列表