geospatial

0热度

1回答

空的SpatialPolygonsDataFrame对象已通过，并且将在Shiny中跳过

对于繁琐的代码，我很抱歉，但需要一些很好的帮助。我将从printscreens 所以开始，5出的6个过滤器在我闪亮的应用程序的工作完美：的问题出现了，当我尝试创建第六过滤器 - 所谓指示灯。并呈现与三个最后一个过滤器相同的数据和格式。基本上我想结合指标过滤前三个指标。并通过三个指标之一来过滤所有7个项目的圈子（基于一个指标）。是否有意义？但最后的代码不起作用。它让我错误警告在polygo

0热度

1回答

与PostGis ST_Area相比，为什么我会得到一个稍微不同的区域？

我试图自动计算OpenStreetMaps上显示的建筑物的面积。为了做到这一点，我通过立交桥API获取了建筑物多边形的坐标。例如 https://overpass-turbo.eu/ [out:json]; way(29858799); out ids geom; ，输出我 { "version": 0.6, "generator": "Overpass API",

0热度

1回答

地理空间加入雪花

我正在寻找某种地理空间加入雪花。目前，运行此查询需要一天时间。我有一个位置数据点（〜10亿条记录）和另一个POI位置表（〜6k）（均以纬度/经度坐标）表。我想绘制在POI的X米内看到多少个数据点。只能使用雪花或任何开源程序。任何帮助表示赞赏。我希望看到一个表与下面的列： POI_id半径数据点

0热度

1回答

Mapbox.js：在日期线两侧的区域上设置地图边界

Mapbox,fitBounds方法似乎不能很好地与大西洋两岸的点（该坐标从跳跃点359至0）。它怎么可能适合阿拉斯加的边界或从香港到SF的航班。使用中心点和自定义缩放级别来编写解决方法也许是可能的，但它似乎是一个过度工作的解决方法，围绕着对边界的不良处理。如果firBounds没有帮助，是否有另一种类似的方法可以实现类似的结果？周围阿拉斯加fitBounds的实施例中，由于阿拉斯加穿过子午

0热度

1回答

AttributeError的：“地图”对象有没有属性“circle_marker”

我应该如何解决这个问题？基本上，我想根据他们的地理信息创建推文地图。所有我收集到的鸣叫是在美国境内，也有地理信息。 Traceback (most recent call last): File "tweet_map.py", line 19, in <module> tweet_map.circle_marker(location=geo, radius=250) Attr

6热度

2回答

如何计算R的土地覆盖面积

基本上，我用ASCII的形式计算了一个全局分布概率模型，例如： gdpm。 gdpm的值全部为0和1之间然后我输入从形状文件的本地地图： shape <- file.choose() map <- readOGR(shape, basename(file_path_sans_ext(shape))) 下一步，我光栅化gdpm，并裁剪使用本地地图： ldpm <- mask(gdpm, ma

0热度

1回答

查找给定的经度，纬度对位于MongoDB中的任何一个多边形中

根据我目前对MongoDB的理解，我们通过传递一组[long，lat]对数组来搜索我们的集合，以检查该文档的[long， lat]对在坐标[[ [], [], [] ]]内退出。 I，E。， db.SomeCollection.find({ "location":{ $geoWithin: { $geometry: { type: "Polygo

2热度

1回答

ggmap + ggplot不会绘制某些值

我想要采用空间数据集，旋转它，并用ggplot/ggmap绘制它们两者。我已经包括了数据集，我的函数用于旋转关于原点的数据集，以及我希望使用的绘图方法。 library(ggplot2) library(ggmap) library(scales) source("DEMOFunctions.R") PlantLAT <- 39.28682 PlantLON <- -96.1172

1热度

1回答

聚焦统计信息（聚焦，“光栅”包，R）的自身函数给出错误的输出

我是R和编程的初学者。所以我觉得这个问题很简单，但我不能找到答案或解决它。我有栅格（100 * 100个单元格）。我需要在移动窗口处通过2D DFT获得谐波幅值的中值（例如，窗口大小= 21）。我发现了光栅包中的焦点功能。对于这个函数，我可以编写自己的函数，它取值为一系列值（窗口中的栅格值）并返回整个窗口的单个值。 r <- raster(matrix(rnorm(10000), nrow =

2热度

1回答

为什么作物有时会在分类栅格上引入NAs？

library(raster) r <- raster('glc2000_v1_1') # http://forobs.jrc.ec.europa.eu/products/glc2000/products/glc2000_v1_1_Grid.zip extent(r) # class : Extent # xmin : -180.0045 # xmax : 179.9955 #