我有一个任务是为一些空间数据产生一个反距离加权估计。我发现,当只有一个已知点和指定的最大距离时,使用idw会导致R(版本3.3.2和3.4.2)崩溃。
library(sp)
library(gstat)
# THIS IS OK MORE THAN ONE KNOWN POINT
X<-c(15,4,20,0)
Y<-c(10,10,10,10)
P<-c(1,2,3,4)
known.df<-data.frame(X,Y,P)
coordinates(known.df)<- ~ X + Y
X<-c(10)
Y<-c(10)
我的问题是:我如何能够掩蔽或裁剪结果,使用R,IDW插值的结果,只包含原始的一组数据点?
在下面的示例中,使用20个随机点在gstat中使用IDW函数插值曲面。从点集得到凸包,并在插值映射上绘制。但是,我想要裁剪地图,所以只有点云内的区域显示插值结果。
谢谢你的建议!
###
# IDW in R
#
library(gstat)
#
set.seed(1234)
x <- rnorm(20)+10
y <- rnorm(20)+10
z <- rnorm(20)+10
#
xyz <-
我有一个包含R2和pvalue的lm模型结果,我在一个条形图中绘制了它们。然后,我用两个离散变量对它们进行了剖分。我想将*放在条形图的顶部,以标记统计意义(pvlue <= 0.05),如下图左下角所示。
我还没有找到关于如何做到这一点的有洞察力的教程。
有办法吗,拜托?
这是我用过的一些code
> head(res_all_s2)
WI aggre_per Season yield_level slope Intercept r.squared
1 R IDW2 Dec Season2 Region II -7.06
我在检查输入是否正确的if语句中遇到了问题。当我输入M或F时,它会打印出“错误的输入”,我不是很清楚。
def check(aw,iw):
if abs(aw-iw)<5 :
print "Your weight is normal."
elif abs(aw-iw)>5 and abs(aw-iw)<15 :
print "Your weight is about normal."
else:
print "Your weight is not normal.
我正在尝试将通常在oracle服务器上运行的查询转换为使用SQL server (使用Teradata SQL Assistant客户端),但我在理解语法更改是什么方面遇到了问题。
目前针对Oracle,它生成具有相关计划工作的“父”设备列表,同时还生成与“子”设备it相关的任何计划工作。
SELECT
EquipmentSurfaceAssetBK AS "DEPARTMENT",
ParentEquipmentCodeBK AS "PARENT",
EquipmentCodeBK AS "POSITION",
EquipmentDescri
我正在学习空间数据格式。我已经创建了一个简单的例子,它使用逆距离加权(IDW)根据初始的点集在网格上插值值。现在我想得到每个初始点的插值值。我认为这应该比较容易,但很难找到一种优雅的方法。换句话说,我想在"geog2.idw“中为40个点(x,y)中的每个点找到最近的网格值。
代码
#adjusted from example created by Brian J. Knaus
#packages
library(gstat)
require(gstat)
require(lattice)
require(raster)
require(sp)
##### ##### #####
我正在尝试编写一个循环,将一个函数应用到一个列表列表中。我想用gstat包的idw命令插入不同国家的温度(最高平均温度/MAXMEAN)。 idw函数需要坐标(经度、纬度)和变量MAXMEAN的信息。这些被组合在数据帧列表"temperatures.coordinates“(看起来像这样:https://imgur.com/cr3PquB)中,这是一个包含37个国家的列表,其中每个列表又被分成12个月(子列表长度为12)。该函数还需要来自grd.list数据帧的信息,该数据帧是一个长度为37的列表,没有子列表。 应用于单个国家/地区和单个月份,代码如下所示: # Create a g
我试图将IDW (逆距离加权)应用于数据库中的不同组。我试图使用dplyr对每个组应用这个函数,但是我在拆分-应用-组合中犯了一个错误。当前函数为每组10个观察返回10个值,但目前dplyr试图在每个变异单元中插入10个返回值,而不是一个新的变异单元格值。
问题可能与函数无关,但不幸的是,我找不到一个显示同样错误的更简单的函数。
我得到一个错误消息,即dataframe已损坏,新列中填充了值。
group N Lat Long Obs idw_val
1 A 1 49.43952 20.42646 11 <dbl[10]>
2 B
我试图修改一个查询,以生成对未来工作的预测,以总结分配给某些行业的工作量,但是当我添加SUM()函数时,我仍然会得到一个3504错误。我不知道我错过了什么?
编辑:包括组和选择,不确定为什么他们没有复制第一.
注: WS.ActivityEstimatedHours的列类型为十进制
SELECT
WO.WorkOrderCodeBK AS "WORKORDER",
WO.WorkOrderPMCodeBK AS "PM",
WO.WorkOrderDescription AS "DESCRIPTION
我在试着提取一个圆球,我一直在犯这个错误
tar: mediatek/custom/common/lk/lcm/bp101wx1_n/bp101wx1_n.c: Cannot open: Not a directory
tar: mediatek/custom/common/lk/lcm/hsd070idw1/hsd070idw1.c: Cannot open: Not a directory
tar: mediatek/custom/common/lk/lcm/nt35582/nt35582.c: Cannot open: Not a directory
tar: mediatek/cust