maptools

    0热度

    1回答

    我想使用dplyr来计算一组lon/lat/timestamp坐标的日出时间,使用maptools的sunriset函数。这是一个可重现的例子。 library(maptools) library(dplyr) pts <- tbl_df(data.frame( lon=c(12.08752,12.08748,12.08754,12.08760,12.08746,12.08748),

    3热度

    1回答

    我是R新手,它是包系统,刚写了我的第一个包,意图在OpenCPU中使用它。 在执行此功能: #' Create a PNG which shows interpolated senseMap Data #' #' @export #' @import jsonlite #' @import sp #' @import gstat #' @import rgeos #' @import

    0热度

    1回答

    我已经创建了一个SpatialPolygonsDataFrame,其中的国家与区域相关联。这里是什么样子https://cloudstor.aarnet.edu.au/plus/index.php/s/RpYr3xyMrmhaGKA(有数据EDU链接): 由于对象太大用,并作为我只是在信息的区域,而不是国家的相关工作,我要合并各地区内的国家。换句话说,我希望为每个区域创建一个多边形。 我试图运行下

    1热度

    1回答

    我有加载为SpatialPolygonsDataFrame在R个空间数据: library(sp) library(tmap) d <- readRDS("data.rds") qtm(d) 对于该示例,从我http://www.gadm.org/download使用的数据为德国。 下面,你会看到由qtm(d)绘制的地图。但是,我想显示自己的数据。我有经度和纬度的位置,我想要计算映射多边

    0热度

    1回答

    我有以下两个shape文件: > summary(precincts1) Object of class SpatialPolygonsDataFrame Coordinates: min max x -74.25545 -73.70002 y 40.49613 40.91540 Precinct Shape_Leng Shape_Area Min.

    0热度

    1回答

    我正在为区域内的各个学校创建地图,并且想要在学校边界周围创建阴影区域。我可以使用QGIS手动完成,但是希望能够在R中生成类似的东西。我现在可以在地图上绘制一些人口统计变量,添加学校位置和边界。有一个循环会自动为每所学校创建一张新地图,并将其插入到knitr报告中。 我无法弄清楚的是如何创建一个图层来遮蔽边界以外的区域。这是我的地图使用R: 但是想使它看起来更像是这样的:但是它是基于内核 我试着从G

    0热度

    1回答

    我试图在R中绘制我的研究区域。到目前为止,我已经下载了人口普查OA边界数据,并且已经将属性数据合并到了我希望查看的OA中。 因此可以从总共1000个OA中说出。我只想看看这些OA中的500个,为此我有外部数据集。 土地是包含1,000个OA的我的空间多边形数据框。它有一列是OA名称“地理编码”。然而,它有它的坐标,因为你运行: coordinates(land) 它产生所有的坐标。 House

    10热度

    2回答

    不妨让猫从袋子里走出来,并说我一直在尝试约一天的时间来找出发生了什么,当我呼叫plot为类SpatialPolygons*的对象。我的直觉是寻找plot.SpatialPolygonsDataFrame,这种方法适用于我见过的任何其他示例&遇到,但这失败了。 简单的例子: US States .shp file这里。 library(maptools) us.states<-readShapeP

    3热度

    2回答

    我正在使用R与maps/mapproject/maptools包来绘制一些地图,并注意到一种行为似乎违反直觉我实际上限制了我能做的事情。 绘制地图欧洲的(与来自ETRS89/ETRS-LCC取出,所以没有冰岛并且还限幅在东方的限制)不指定任何投影: library(maps) map("world",xlim=c(-10.67,34.5),ylim=c(31.55,71.05), interio

    1热度

    1回答

    我有我的情节和传说正确的颜色范围的问题。 这是代码使用: data.ch4 <- read.csv2("v42_CH4_1970_TOT.txt",skip = 3,stringsAsFactors = FALSE, header = F) num_data <- data.frame(data.matrix(data.ch4)) library(maptools) library(lat