14得票1回答
使用sf::st_read读取shape文件无法捕获UTF8编码。

我想读取一个以UTF8编码的shape文件。使用rgdal::readOGR读取时,它可以正常工作,但是sf::st_read无法获取正确的编码。有什么建议可以解决这个问题吗? 为了提供可重现的示例,您可以从这里下载我正在尝试读取的shape文件。 使用rgdal::readOGR进行读取...

14得票3回答
R如何在一定容差范围内叠加点和多边形?

我希望使用R叠加空间点和多边形,以便将地理区域的属性分配给这些点。 通常我使用sp包中的over命令来实现。但问题是我正在处理全球范围内大量的地理位置事件,在某些情况下(特别是在沿海地区),经纬度组合略微超出了国家/地区的边界。 以下是一个可重现的示例,基于这个非常好的问题。## exampl...

14得票1回答
将泰森多边形和地图结合起来

我想将 Voronoi 多边形与地图结合起来,以便稍后用于空间分析。我有一些点和形状文件要组合,然后保存为形状文件/空间多边形。为了获得 Voronoi 多边形,我使用此主题中的函数。 以下是我的代码: coords<-data.frame(LONG=c(16.9252,16.936...

14得票1回答
从SpatialPolygons和其他sp类中提取特征坐标

包sp提供了许多用于不同空间概念(点、线、多边形)的类。对于某些类,例如SpatialLines,访问要素坐标很简单。所有示例都来自于各自类别的帮助页面。l1 = cbind(c(1,2,3),c(3,2,2)) l1a = cbind(l1[,1]+.05,l1[,2]+.05) l2 = ...

14得票3回答
在R中将起始坐标和结束坐标转换为空间线

我有一组看起来像这样的开始和结束坐标: begin.coord <- data.frame(lon=c(-85.76,-85.46,-85.89), lat=c(38.34,38.76,38.31)) end.coord <- data.frame(lon=c(-85.72,-8...

14得票1回答
如何从GEOMETRYCOLLECTION返回MULTIPOLYGON?

我有一个全球国家的数据集,希望能够在本初子午线上进行分割,并重新调整数据以关注太平洋。 我尝试使用简单要素(sf)进行此操作,但遇到了无法解决的对象类型问题。 为了拆分数据,我尝试了以下方法: st_wg84 <- "+proj=longlat +datum=WGS84 ...

14得票2回答
R中的GIS地图自动标签放置

我正在使用sf包(及相关包)在R中制作GIS地图,以读取shapefile文件,并使用ggplot2(及相关工具)进行绘图。这很好用,但是我找不到自动/编程创建河流和道路等要素的标签位置的方法。这些要素通常是线串,形状不规则。请参见维基媒体上的示例图像。 ggrepel包可以很好地自动标...

14得票3回答
在R中创建网格以进行gstat中的克里金插值

lat long 7.16 124.21 8.6 123.35 8.43 124.28 8.15 125.08 考虑这些坐标,这些坐标对应于测量降雨数据的气象站。 R中gstat软件包的介绍使用了meuse数据集。在本教程https://rpubs.com/nabil...

13得票3回答
MySQL对POINT列的INSERT/UPDATE操作

我正在尝试用我的国家的地理位置来填充我的数据库。我的其中一个表有4个字段:ID[主键]、纬度、经度和geoPoint。EDIT `SCDBs`.`Punto_Geografico`; SET @lat = 18.469692; SET @lon = -63.93212; SET @g =...

13得票2回答
NHibernate.Spatial和Sql 2008地理类型

我目前正在处理使用SQL Server 2008地理类型的项目。 由于该项目的大部分部分使用NHibernate作为ORM,我想知道使用/扩展NHibernate.Spatial与此类型会有多少工作量(就我所看到的而言,Nh Spatial仅支持几何类型)。 顺便问一下,还有其他支持此SQL类...