请选择 进入手机版 | 继续访问电脑版
2123 8

[问答] 请教空间插值的问题 [推广有奖]

  • 0关注
  • 0粉丝

大专生

6%

还不是VIP/贵宾

-

威望
0
论坛币
0 个
通用积分
7.4120
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
353 点
帖子
24
精华
0
在线时间
28 小时
注册时间
2022-4-8
最后登录
2023-12-14

相似文件 换一批

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

感谢您参与论坛问题回答

经管之家送您两个论坛币!

+2 论坛币
已知一组点数据如下:
longi        lati        AMM
124.2500         27.0833         0.717763911
124.4167         27.0833         0.565288714
124.2500         27.2500         0.715888277
124.4167         27.2500         0.723441372
125.2500         27.4167         0.902380367
125.4167         27.4167         0.789868463
125.5833         27.4167         0.80516665
125.7500         27.5833         0.867977949
125.2500         27.7500         0.931652271
127.4167         31.9167         0.164210845
122.4167         26.5833         0.834836401
122.5833         26.5833         0.865860023
122.9167         26.7500         0.647638529
123.9167         26.7500         0.790268065
122.5833         27.0833         0.767969355
124.4167         27.0833         0.887081911
123.5833         27.2500         0.920635416
123.9167         27.5833         0.861132427
124.0833         27.7500         0.964586229
124.4167         27.7500         0.892144971
122.5833         26.0833         0.432547089
122.9167         26.0833         0.457594433
122.4167         26.2500         0.631898365
122.4167         26.4167         0.570757814
122.5833         26.4167         0.566880461
122.4167         26.5833         0.626675953
122.4167         26.7500         0.617377349
122.9167         26.9167         0.68981381
122.0833         27.7500         0.406830771


请问如何对其进行空间插值形成面域绘制成图
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

关键词:long LAT ATI

本帖被以下文库推荐

llb_321 在职认证  发表于 2022-5-5 17:28:12 |显示全部楼层 |坛友微信交流群
library(RandomFields)
这个包,可以用空间点数据,实现克里金法插值,然后可以利用ggplot2绘图
已有 2 人评分论坛币 收起 理由
cheetahfly + 10 热心帮助其他会员
jiangbeilu + 5 精彩帖子

总评分: 论坛币 + 15   查看全部评分

使用道具

llb_321 发表于 2022-5-5 17:28
library(RandomFields)
这个包,可以用空间点数据,实现克里金法插值,然后可以利用ggplot2绘图
老师您好,library这个包,显示了错误如下:
WARNING: Rtools is required to build R packages but is not currently installed. Please download and install the appropriate version of Rtools before proceeding:

https://cran.rstudio.com/bin/windows/Rtools/
Warning in install.packages :
  package ‘RandomFields’ is not available for this version of R

A version of this package for your version of R might be available elsewhere,
see the ideas at
https://cran.r-project.org/doc/manuals/r-patched/R-admin.html#Installing-packages

使用道具

llb_321 在职认证  发表于 2022-5-5 21:51:06 |显示全部楼层 |坛友微信交流群
学R语言的星星 发表于 2022-5-5 18:14
老师您好,library这个包,显示了错误如下:
WARNING: Rtools is required to build R packages but is  ...
抱歉,我刚查了一下,这个包已经不在CRAN中了。
貌似可以用fields包,不过这个包我没用过。

使用道具

llb_321 发表于 2022-5-5 21:51
抱歉,我刚查了一下,这个包已经不在CRAN中了。
貌似可以用fields包,不过这个包我没用过。
好嘞,谢谢老师!

使用道具

llb_321 在职认证  发表于 2022-5-6 16:33:13 |显示全部楼层 |坛友微信交流群
学R语言的星星 发表于 2022-5-6 16:10
好嘞,谢谢老师!
交流而已,不是什么老师。
那个fields包,我试了一下,可以插值,但是surface效果不理想。
您在CRAN的packages列表里搜一下Kriging,有一些包可能有用。

使用道具

llb_321 在职认证  发表于 2022-5-6 19:59:54 |显示全部楼层 |坛友微信交流群
可以看看LatticeKrig包。
假设你的数据为dd
  1. library(LatticeKrig)
  2. res <- LatticeKrig(dd[,1:2], dd[,3], nlevel = 1) #插值
  3. surface(res) #绘制面域
  4. points(dd[,1:2]) #绘制原始数据点
复制代码

这是最简单的用法,那个nlevel参数是插值的阶数。
这个包挺复杂的,但效果我个人感觉没有randomFields好,没办法,那个包被清掉了,github上也没有了。

使用道具

llb_321 在职认证  发表于 2022-5-7 13:09:15 |显示全部楼层 |坛友微信交流群
这是我用soil数据,fields包Tps()插值的结果,途中黑点是原始数据点。
2332232.PNG

使用道具

llb_321 发表于 2022-5-7 13:09
这是我用soil数据,fields包Tps()插值的结果,途中黑点是原始数据点。
特别感谢!!

使用道具

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

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-3-29 19:48