-
1 votesanswersviews
使用gstat进行局部变异函数的局部块克里格法
我一直无法使用R中的gstat包找到本地变异函数的局部块克里金特有关的信息 . 澳大利亚精确农业中心有一个名为VESPER的免费软件可以做到这一点，并且从我读过的内容应该可以在R中，我可以使用一些帮助将一个for循环放在一起，使gstat函数在本地工作 . 使用meuse数据集作为示例，我已经能够计算并将全局变异函数拟合到数据集： library(gstat) data(meuse) ... -
0 votesanswersviews
使用R将空间点数据集与空间栅格数据集合并（主数据集采用SP点格式）
我正在研究使用R的空间数据集 . 数据描述 我的主数据集采用 SpatialPointsDataFrame 格式，并且每个月都有表面温度数据（列名称 - "ruralLSTday"，01231898） . 数据摘要如下所示： 主数据 - （以SpatialPointsDataFrame格式） TOWN_ID ruralLSTday ruralLSTnight year mon... -
4 votesanswersviews
如何使用gstat预测进行普通克里金法
我试图在R中编写一个使用gstat库的代码来创建插值 . 我已经阅读了gstat手册，基于互联网上的一些例子，我设法编写了这段代码（这只是一部分）： g <- gstat(id="tec", formula=TEC ~ 1, data=data) ##I create an object v <- variogram(g) # plot the empirica... -
1 votesanswersviews
来自克里金数据的ggmap上的热图
我想在ggmap上绘制热图 . library(ggmap) turku<-get_map('turku', zoom=13) turkumap<-ggmap(turku, extent="device", legend="topleft") turkumap turkumap+geom_density2d(mapping=aes... -
1 votesanswersviews
如何在R中使用gstat拟合半变异函数模型？
我有一个csv文件文件包含3月1日下午12点的大气PM10浓度数据 . please, download.我想在R中使用gstat包绘制一个半变异函数 . 我试着在R中编写这些代码 . 但是对于这些数据，我不能适应模型 . library(sp) library(gstat) seoul3112<-read.csv("seoul3112.csv") ... -
0 votesanswersviews
R：差异变异函数拟合，克里金结果差
我试着在雅加达湾做一个克里金 . 我有一套测量点，有适当的坐标和属性（pH，盐度，......） 为了做克里金法，我首先需要为变差函数找到一个模型 . 当我使用"variogram"函数输出不完美但是应该没问题，但是当我尝试拟合变异函数时，我得到一个警告信息：在fit.variogram中（ph.vgm，model = vgm（0.12，"Sph"，0.1）... -
0 votesanswersviews
Automap包：负值和变异函数拟合
两天前我开始使用automap软件包进行自动变异函数建模并玩弄它（手动操作很繁琐而且没有太多优点......虽然我觉得很奇怪自动化似乎真的很喜欢Matérn函数用于对变差函数建模） . 我想到的是数据集中存在负值，而这些值未在此包生成的 Map 上显示 . 这是我刚刚粘在一起的数字：Nr . 1是使用autoFitVariogram和autoKrige后的输出 . 它似乎没有负面 Value ... -
6 votesanswersviews
如何反向转换正态分数转换后的数据
我每天有61个测量站的降雨量，为期12年（8000平方公里） . 目标是创造5公里和25公里分辨率网格日降雨量产品 . 由于车站的数量很少，并且即使在雨季也不是所有车站都有雨，我选择使用气候变异图 . 典型日（irain）的雨量计测量如下，具有由NA表示的少量缺失值 . 7.8 4.4 15.4 19.1 5.8 0 42 6.4 21 21 0 0 ... -
0 votesanswersviews
提取变差函数模型的拟合值
如何提取变异函数模型的拟合值（不是参数！）？我打算在Excel软件中绘制理论变差函数 . 我在gstat包中使用了fit.variogram，如下所示： library(sp) data(meuse) coordinates(meuse) = ~x+y vgm1 <- variogram(log(zinc)~1, meuse) fit.variogram(vgm1, vgm(1, "...