问题标签 [maptools]

我一直在使用 ArcMap 来访问空间数据服务器上的 GIS 数据。我想弄清楚如何在 R 中做同样的事情。

我知道如何将 shapefile 读入 R。我已成功使用 maptools 和 rgdal 打开和映射本地存储的 shapefile（例如 http://www.nceas.ucsb.edu/scicomp/usecases/ReadWriteESRIShapeFiles）

我的问题是数据不是存储在本地，而是存储在应用程序服务器上。我相信这是一个Oracle数据库。我已经获得了有关 1. 服务器 2. 实例（一个数字） 3. 数据库 4. 用户和 5. 密码的信息。通常，我会举一个例子，但外部用户能否访问服务器是值得怀疑的。

例如这里是如何在 R 中读取和绘制本地文件

这 ”。” 指向本地目录。我将如何更改它以访问服务器上的数据？代码的实际语法会很有帮助。

我有一个包含生物群及其地理信息的数据集：

结果2

我使用 maptools 在地图上绘制这些点，这里的代码：

但我不知道如何为每个组划定多边形，并更改所研究组的颜色选项。

我有一些坐标 lat、lon 和值 Z 的数据点（数百万个数据点）。每个坐标点代表一个大小为 0.0005 的网格的中心。我想用颜色代表Z值来绘制它们。但是，我的数据不是网格数据，它们是分散的数据点。我知道另一个函数“geoshow”也可以制作类似的地图，但用于网格数据。使用 matlab 将散点数据转换为网格数据时内存不足。所以我使用 pcolor 实现了这一点。下面的代码：

对于小数据集，它工作得很好。但是，它将花费太多时间来处理更大的数据集。是否可以在没有网格的情况下更快地绘制散点数据？任何功能 pcolor ， geoshow 或其他东西都可以。谢谢你的帮助！

我有一堆包含在类的对象中的 shapefile：

我想以km2 获得这些shapefile 的面积。由于这些 shapefile 分布在全球各地，因此我使用“等面积”投影对 shapefile 进行了投影，即：兰伯特圆柱等面积“+proj=cea +lon_0=-80”，以确保这些区域不会受纬度影响而变形。

我这样做是这样的：

我以这种方式计算面积：

但是，我获得的区域很奇怪……值如下：7.861879e-02 1.819334e+00 4.510745e+00。他们只是没有意义...

我究竟做错了什么？我怎样才能以平方公里为单位获得这些区域？（考虑到这些 shapefile 分布在全球...）

任何帮助将不胜感激。非常感谢您提前。

蒂娜。

我猜有一个简单的解决方案可以解决我一直遇到的问题，但是我遇到了一些麻烦。

我正在尝试转换以下map对象：

使用函数进入SpatialPolygon对象：map2SpatialPolygons

我不断收到以下错误：

经过一番研究，在应用该函数后，ID 的长度似乎为 63，而map对象的长度为 169 .NAmat2xyList(cbind(map$x, map$y))（我无法找到其来源）。

有人有想法么？这是usa地图对象的结构：

这与我之前提出的一个问题有关：

从 Shiny (R) 下载 png

我现在已经创建了多个闪亮的图并下载了它们，但是对于使用 maptools 包的图，我只得到一个空的 png。

这是 Shiny 中的错误还是我的代码有问题？

以下是我的服务器文件的相关摘录：

我确信这个问题已经在其他地方得到了回答，但我无法通过搜索提出这个问题。

我有代表一个国家内的城市以及每个城市的人口的点。我还有一个县的多边形文件。我想找到每个县内最大城市的位置。

如何才能做到这一点？

这是一些数据

结构（列表（国家= c（“我们”，“我们”，“我们”，“我们”，“我们”，“我们”，“我们”，“我们”，“我们”，“我们”，“我们“,
“我们”, “我们”, “我们”, “我们”, “我们”, “我们”, “我们”, “我们”, “我们”, “我们”, “我们”, “我们”, "us", "us"), City = c("cabarrus", "cox store", "cal-vel", "briarwood townhouses", "barker heights", "davie
十字路口”、“蟹点村”、“杜鹃花”、“切斯特菲尔德”、“查尔斯蒙特”、“康纳”、“三叶草花园”、“科里赫高地”、“卡利森”、“克雷斯特维尤英亩”、“克莱格”、“迦南”公园”，“尚蒂伊”，“贝尔格莱德”，“布里斯十字路口”，“虚张声势”，“巴特纳”，“底部”，“班迪”，“博斯蒂安高地”），AccentCity = c（“Cabarrus”，“Cox Store” , “Cal-Vel”, “Briarwood Townhouses”, “Barker Heights”, “Davie Crossroads”, “Crab Point Village”, “Azalea”, “Chesterfield”, “Charlesmont”, “Connor”, “Clover Garden”, “ Corriher Heights”、“Callisons”、“Crestview Acres”、“Clegg”、“Canaan Park”、“Chantilly”、“Belgrade”、“Brices Crossroads”、“Bluff”、“Butner”、“Bottom”、“Bandy”、“Bostian Heights”）、地区= c（“NC”，“NC”，“NC”，“NC”，“NC”，“NC”，“NC”，“NC”，“NC”，“NC”，“NC”，“NC” 、“NC”、“NC”、“NC”、“NC”、“NC”、“NC”、“NC”、“NC”、“NC”、“NC”、“NC”、“NC”、“ NC"), 人口 = c(NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, A_integer_,NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_), Latitude = (35.2369444, 35.275, 36.4291667, 35.295, 35.3111111, 35.8319444, 34.7602778, 35.58, 35.81, 5.9341667, 35.7419444, 36.1883333, 35.5605556, 35.0841667, 35.0213889, 35.8655556, 36.2761111, 36.3016667, 34.88, 34.8186111, 35.8377778, 36.1319444, 36.4747222, 35.6419444, 35.7544444), Longitude = c(-80.5419444, -82.0352778, -78.9694444, -81.5238889, -82.4441667, -80.535 , -76.7305556, -82.4713889, -81.6611111, -81.5127778, -78.1486111, -79.4630556, -80.635, -76.7255556, -80.5427778, -78.8497222, -79.7852778, -76.1711111, -77.2352778, -78.1016667, -82.8580556, -78.7569444, - 80.7741667, -81.09, -80.9294444)), .Names = c("Country", "City", "AccentCity", "地区”, “人口”, “纬度”, “经度”), row.names = c(544L, 889L, 551L, 434L, 190L, 975L, 894L, 147L, 717L, 700L, 831L, 773L, 862L, 559L, 915L, 753L, 584L, 695L, 262L, 437L, 372L, 537L, 406L, 178L, 02L), 类别 = "data.frame")

以及在北卡罗来纳州阅读的一些代码

我想找到每个县内人口最多的城市。对不起，我没有可重复的例子。如果我这样做了，我会得到答案！

我正在开发一个原始的闪亮应用程序，它将映射来自苏格兰开放数据项目的一些数据。我制定了 SPARQL 查询，这些查询制造了一个类似于下面提供的提取的数据框

我有强化的 shapefile，我正在尝试使用以下代码映射它们：

代码返回错误

如果我理解类似的讨论（1、2），则问题与值类有关。但是，我的理解是as.factor()应该解决这个问题，但显然情况并非如此。该应用程序的完整代码包含在下面。当然，我会很感激任何帮助。

服务器

用户界面

我的代码中有这一行来为点密度图创建点。我滥用的函数是 maptools 包中的 dotsInPolys 。当我运行它时出现此错误，我不确定它是什么意思。任何人都可以帮忙吗？

NSWdots <- dotsInPolys(NSW, as.integer(datajoin$pop.10))

NSW 是一个 shapefile，datajoin$pop.10 是一个数字向量。

给定的错误是：

Error in if (x[i] > 0) { : missing value where TRUE/FALSE needed

我有一个情节的一部分wrld_simpl

我想给一些国家上色，比如说德国为红色，丹麦为绿色。我该怎么做？我尝试了这样的事情，但它没有按照我想要的方式工作：

谢谢！