-
0 votesanswersviews
（R）“收缩空间多边形” - 基于核密度估计阈值提取栅格/空间网格
我正在尝试使用R将我的shapefile从LHS“转换”为RHS，如下所示 . 原始数据（如LHS所示）取自2016年加拿大人口普查 . 它包括加拿大不列颠哥伦比亚省萨里市的所有传播区域（最小的基地之一） . 形状文件可以下载here . 将其转换为RHS的基本原理是因为典型的区域插值函数，例如 st_interpolate_aw 来自 sf 假设值在单位面积上是常数 . 我希望为不同的地理单位... -
0 votesanswersviews
子集化空间权重矩阵.listw时出错
我试图运行空间自回归模型（这里的一般形式）： SAR1=lagsarlm(log(dependentvariable)~X1+X2+X3, data=DATANAME, listw=SW, na.omit, type="lag", method="eigen",quiet=NULL,interval=NULL,tol.solve=1.0e-10,trs=NU... -
2 votesanswersviews
使用gIntersect查找多个SpatialPolygons的成对交集？
我想生成一个空间对象，其中包含空间多边形列表的所有成对交叉点 . 手动使用gIntersect只需两层的交叉点很容易，但我想同时找到所有成对交叉点 . 我有13个多边形，因此有156种不同的对组合来检查重叠 . 无论是lapply还是for循环似乎都可以工作，但我想我需要一个包含所有可能的空间多边形组合的矩阵 . 以下是数据子样本的要点： https://gist.github.com/dwwol... -
0 votesanswersviews
R-计算点和线之间的最小距离
我有一大堆包裹的多边形形状文件（50个县），我想循环遍历每个包裹并计算从每个包裹的中心到最近的铁路的距离，全国存储在线形文件中 . 我提出了两种不同的措施：来自geospheres包的point2Line和来自rgeos包的gDistance . 每个人都有不同的问题 . 对于point2line，我的代码看起来像这样： rail_projection <-spTransform(rail_... -
1 votesanswersviews
带有splm的空间面板数据：单变量回归时出错
我正在尝试使用R的“splm”包来估计具有面积数据的空间自回归模型 . 当尝试在拟合多变量之前执行固定效应单变量回归时，我偶然发现了以下错误： colnames <-（* tmp *，value =“log（pcap）”）出错：尝试在小于二维的对象上设置'colnames' 使用随机效果或多个自变量时不会发生这种情况 . 根据splm的包文档，这也可以使用示例数据重现： library(... -
1 votesanswersviews
R空间面板固定效果：lag.listw中的错误，对象长度不同
我正在使用R包“splm”估计空间面板数据的固定效应模型 . 强有力 balancer 的小组由77个县和31年的观察组成 . 因变量是一种作物的县级产量 . 当我对随时间变化的趋势和单个固定效果进行回归时，我得到了对象长度不同的错误，即使我有完美 balancer 的面板并且使用相同的一组县来创建空间权重矩阵 . 下面我生成了一个类似于我的数据的“面板数据集”，并且还生成了具有随机数的县的经度和... -
0 votesanswersviews
R中的空间面板回归：spgm出错
我对下面显示的空间回归代码有一个问题（ Code Example ） . 运行回归后，我得到以下错误 Error in listw %*% as.matrix(ywithin) : Cholmod error 'X and/or Y have wrong dimensions' at file ../MatrixOps/cholmod_sdmult.c, line 90 当我从回归中删除空... -
1 votesanswersviews
R中面板数据的空间回归
我有一个面板数据集，有几百个区域，大约10年，区域的空间数据 . 我用 spdep 包创建了一个权重矩阵（通过标准方式，然后， nb2listw ） . 因此，我有一个矩阵，每个区域都有权重（相对于其他区域） - 但每个区域只表示一次 . 我想从 spdep 包中运行一些空间回归（ lagsarlm ， errorsarlm ），但是我收到一个错误： Error in subset.listw(l... -
2 votesanswersviews
Geoserver无法识别空间列
我正在使用Geoserver和SQL Server 2008 . 我有一张表有[地理]类型的列 . 我能够看到Geoserver中显示的表格（我已经安装了SQL服务器的扩展名），但它将列识别为byte []类型而不是地理位置 . 有人可以建议我解决这个问题 . 从上一篇文章开始，我尝试多次删除并重新附加Geoserver中的数据库 . 此外，SQL服务器显示数据是地理类型并正确映射它 . 只是Ge... -
3 votesanswersviews
空间权重：不对称邻接矩阵？
我正在创建一个邻接矩阵，用于在R中进行空间分析 . 数据是美国大陆的所有县 . 我有来自美国人口普查Tiger文件的县空间多边形 . 我能够创建邻居列表，它是对称的 . 但是当我将其转换为邻接矩阵时，它不是对称的 . 这是一个问题，因为我的目标是使用 ngspatial::autologistic 运行空间自动模型，并且我得到一个错误，我必须提供对称的二元邻接矩阵 . 这是我创建邻接矩阵的R代码：... -
3 votesanswersviews
R中的空间面板回归：不一致的空间权重？
我试图用splm包在R中运行空间面板回归 . 所以我有多边形，随着时间推移汇总数据，我想看看因变量如何受其他变量的影响，这些变量也随着时间的推移而变化 . 我有546个区域有许多变量，但为了测试它是如何工作的，我把3个多边形的数据子集，包括用于计算权重的shapefile和数据 . https://drive.google.com/file/d/0B4SK0f2zZUKxZ0dDU2lnclB2M... -
15 votesanswersviews
刻面ggplot的不同图例和填充颜色？
很抱歉没有包含我的问题的任何示例数据 . 我找不到轻松生成示例形状文件的方法 . 希望有经验的 ggplot 用户可以从下面的描述中看到我想做的事情 . 我有： A data.frame X，包含有关样本图的信息（ plotid ， var1 ， var2 ， var3 ， var4 ，...） 多边形shapefile Y ，包含样本图的空间信息 将shapefile Y （... -
4 votesanswersviews
如何使用R中的ggplot2在投影 Map 上绘制插值数据
我想使用ggplot2在投影 Map 上绘制一些插值数据，我已经在这个问题上工作了几个星期 . 希望有人可以帮助我，非常感谢 . shapefile和数据可以在https://www.dropbox.com/s/8wfgf8207dbh79r/gpr_000b11a_e.zip?dl=0和https://www.dropbox.com/s/9czvb35vsyf3t28/Mydata.rdata... -
1 votesanswersviews
提取SpatialPolygonsDataFrame（索引）的特定多边形的栅格值
我有一个包含120个多边形和一些相关数据的SpatialPolygonsDataFrame . 现在，我想分别提取每个多边形内栅格上的值的平均值 . 我成功地绘制了单个多边形： plot(SpatialPolygons(SPdataframe@polygons)[i]) 但它不能以相同的方式提取值： extract(raster, SpatialPolygons(SPdataframe@poly... -
1 votesanswersviews
如何在R中使用gstat拟合半变异函数模型？
我有一个csv文件文件包含3月1日下午12点的大气PM10浓度数据 . please, download.我想在R中使用gstat包绘制一个半变异函数 . 我试着在R中编写这些代码 . 但是对于这些数据，我不能适应模型 . library(sp) library(gstat) seoul3112<-read.csv("seoul3112.csv") ... -
1 votesanswersviews
如何选择krige插值的时空变异函数参数
最近我用R进行空间 - 时间插值 . 我在包 gstat 中选择函数krigeST，它在时空数据中使用克里金法 . 这是我的步骤： 首先，我创建一个名为stfdf的STFDF对象： stfdf<-stConstruct(observation,space = list(values = 1:ncol(observation)),time = Node$Time,SpatialObj = ... -
7 votesanswersviews
使用spplot过度绘制两个SpatialPolygonsDataFrame
我有一堆数据，我在县一级绘制，without borders . 我想加入州界 . 我有一个状态shapefile（多边形），但 spplot 似乎没有任何方法可以添加到上一个 Map 之上 . 有没有办法做到这一点，而不重写面板功能，以采取两个SPDF（这似乎非常专业的可能是其他人的问题）？ 这是一个可重复的例子： library(sp) Srs1 = Polygons(list(Polygon... -
59 votesanswersviews
通过R中的属性对SpatialPolygonsDataFrame（即删除多边形）进行子集化的简单方法
我想根据@data数据框中的相应属性值从SpatialPolygonsDataFrame对象中删除一些多边形，以便我可以绘制简化/子集化的shapefile . 到目前为止，我还没有找到办法做到这一点 . 例如，假设我要删除面积小于30000的world shapefile中的所有多边形 . 我将如何进行此操作？ 或者，同样，我如何删除Antartica？ require(maptools) g... -
25 votesanswersviews
裁剪SpatialPolygonsDataFrame
我有两个 SpatialPolygonsDataFrame 文件：dat1，dat2 extent(dat1) class : Extent xmin : -180 xmax : 180 ymin : -90 ymax : 90 extent(dat2) class : Extent xmin ... -
0 votesanswersviews
当行数不相等时覆盖SpatialPointsDataFrame和SpatialPolygonsDataFrame [关闭]
我想知道这是否可能在R. 我有一个SpatialPointsDataFrame，其中包含分配了属性的点 . 我想用SpatialPolygonDataFrame覆盖此SpatialPointsDataFrame . 以下是我正在使用的代码 . 只是为了澄清多边形和点数据框都在同一个投影协调系统中 . buptpol <- SpatialPolygonsDataFrame(bu_proj_po... -
6 votesanswersviews
在SpatialPolygonsDataFrame上绘制SpatialPointsDataFrame
我有一个SpatialPolygonsDataFrame，我使用以下命令绘制： spplot(milanoNILNew, "percPolVita",at = c(0,0.0696,0.08979,0.0974,0.1116,0.181),col.regions =pal,main="Quota Clienti Unipol con Polizze Vita &quo... -
3 votesanswersviews
为什么可以使用`$`来子集SpatialPolygonsDataFrame？
当我有一个SpatialPolygonsDataFrame对象时，我知道我可以通过两种方式获取数据： spatial_df@data$column spatial_df$column 但是我不明白为什么第二种方式是可能的 . 我以为我必须使用 @ 访问 data 插槽？这是关于 SpatialPolygonsDataFrame 类的独特之处，还是关于S4对象的一般情况？ 一个可能的答案是在sp ... -
1 votesanswersviews
R：根据多边形合并两个SpatialPolygonsDataFrame
我有2个SpatialPolygonsDataFrames（data1和data2） . 多边形代表建筑物 . data1是所有建筑物，data2是一些特定的建筑物 . 我想将data2 @ data中的所有变量添加到data1 @ data（对于属于SpatialPolygonsDataFrame的建筑物） . 我已经找到了一种通过多边形聚合或合并两个SpatialPolygonsDataFr... -
1 votesanswersviews
将SpatialPointsDataFrame与SpatialPolygonsDataFrame组合错误：超出最大返回密集矩阵大小
我试图在SpatialPolygonsDataFrame（info）上组合一个1000x1000m正方形的SpatialPointsDataFrame（网格）来聚合每个网格方块内的点的所有信息 . 我试过代码： combined <- intersect(info, grid) 但我记得这个错误： Error in RGEOSBinPredFunc(spgeom1, spgeom2, byi... -
2 votesanswersviews
R中的邻近空间滤波
我有一个物种的出现点，我想消除潜在的采样偏差（某些地区的点密度可能比其他地区大得多） . 实现此目的的一种方式是最大化不小于彼此的特定距离X的点的子集 . 从本质上讲，我会防止点太靠近彼此 . 是否有任何现有的R功能可以做到这一点？我搜索了各种空间包，但没有找到任何东西，也无法弄清楚如何自己实现 . 可以下载示例事件点数据集here . 谢谢！ -
-1 votesanswersviews
用于空间点模式分析的Kcross - 在R中
我'm trying to use Kcross for spatial point pattern analysis. I wish to compare between species A and species B, from my dataset called '鸟' . 感兴趣的变量是 species_name . 它只有2个级别 - species A 和 species B ... -
0 votesanswersviews
用于geohash的Oracle Spatial 12C MD.HHENCODE无法正常工作
我试图在Oracle Spatial 12C中使用MD.HHENCODE（或MDSYS.HHENCODE）空间函数 . 基于我能找到的唯一文档，我无法获得编译功能 . 这是我的功能代码： CREATE OR REPLACE FUNCTION linear_key ( p_shape in mdsys.sdo_geometry, ... -
0 votesanswersviews
如何过滤/剪辑WMS提供的栅格 - Geoserver
我正在尝试过滤/剪辑Geoserver提供的Web Map 服务提供的栅格 . 我不想基于多边形进行剪辑;我想根据一个值进行过滤，使得低于此值的所有光栅像素都是黑色或透明的 . 据我所知，我不能使用cql_filter，因为它仅用于WFS功能请求，而不能用于WMS Raster图像 . 也许你知道如何解决我的问题 . -
0 votesanswersviews
R：计算两个顶点之间的单个最短路径
目前，我正在开展一个涉及纽约出租车数据的项目，在这个项目中，我可以获得一个人在网络中被接送的地方 . 我正在使用 ESRI shapefile ，我可以使用 shp2graph 包将其作为 igraph 对象加载到R中;我需要利用Dijkstra算法（或类似的最短路径算法）来找到两个给定顶点之间的单个最短路径 . 我认为 igraph 包的 get.shortest.paths() 方法将是我... -
2 votesanswersviews
使用randomForest，Caret和factor变量预测栅格时出错
我试图用randomForest和插入符号包预测栅格图层，但在引入因子变量时失败 . 没有因素，一切正常，但一旦我带来一个因素，我得到错误： Error in predict.randomForest(modelFit, newdata) : Type of predictors in new data do not match that of the training data. 我在下面创建了...