如何使用rgdal包中的readOGR从shapefile中读取结果？
readOGR函数是rgdal包中的一个函数，用于从shapefile中读取空间数据。它可以读取矢量数据，包括点、线、面等几何对象，并将其转换为R中的Spatial类对象。
使用readOGR函数从shapefile中读取结果的步骤如下：
- 安装和加载rgdal包：
install.packages("rgdal")
library(rgdal)- 使用readOGR函数读取shapefile文件：
shapefile_path <- "path/to/shapefile.shp"
result <- readOGR(dsn = shapefile_path, layer = "layer_name")其中，dsn参数指定shapefile文件的路径，layer参数指定要读取的图层名称。
- 查看读取结果：
print(result)这将显示读取的空间数据的基本信息，包括几何类型、坐标系等。
readOGR函数的优势：
- 支持读取各种常见的矢量数据格式，如shapefile、GeoJSON、KML等。
- 能够处理大规模的空间数据，具有较高的性能和稳定性。
- 提供了丰富的空间数据处理和分析功能，如空间查询、空间拓扑关系分析等。
readOGR函数的应用场景：
- 地理信息系统（GIS）数据处理和分析。
- 空间数据可视化和地图制作。
- 地理空间分析和模型建立。
推荐的腾讯云相关产品和产品介绍链接地址：
- 腾讯云地理信息系统（GIS）服务：https://cloud.tencent.com/product/gis
- 腾讯云地图服务：https://cloud.tencent.com/product/tianditu
- 腾讯云大数据分析服务：https://cloud.tencent.com/product/dla
相关·内容
3回答
我想从shapefile中获取纬度和经度。到目前为止，我只知道如何读取shapefile。
library(rgdal)
centroids.mp <- readOGR(".","35DSE250GC_SIR")
但是如何从centroids.mp中提取经度和纬度呢？
浏览 0提问于2013-05-09得票数 1
回答已采纳
2回答
我有一些非常大的shapefile。我可以使用rgdal函数readOGR将它们读取到SpatialPolygonsDataFrame中，但每个文件都需要很长时间。实际上，我只对@data插槽中显示的data.frame感兴趣。有没有办法只读取数据，跳过资源密集型多边形？
示例代码：
## State of Alabama census blocks (152 MB compressed, 266 MB uncompressed)
shpurl <- "http://www2.census.gov/geo/tiger/TIGER2011/TABBLOCK/tl_2011_01_t
浏览 2提问于2012-11-15得票数 8
回答已采纳
1回答
我想要建立一个世界各国首都的距离矩阵。我已经建立了这样的国家中心线的距离矩阵。Shapefile是从派生的。
library(rgeos)
library(rgdal)
shapefile <- readOGR("./Map/World Map", layer = "TM_WORLD_BORDERS-0.3")
centroids <- gCentroid(shapefile, byid = TRUE, as.character(shapefile@data$UN)
dist_matrix <- as.data.frame(geosphere
浏览 1提问于2016-12-20得票数 0
回答已采纳
我试图使用ggplot2和瑞士的一个公共图形文件(它包含在：PLZO_PLZ.shp中)创建一个choropleth地图。
library(rgdal)
library(ggplot2)
library(ggmap)
library(broom)
# Load Shapefile
shapefile <- readOGR(file.choose())
# Next the shapefile has to be converted to a dataframe for use in ggplot2
shapefile_df <- fortify(shapef
浏览 70提问于2021-05-19得票数 0
2回答
我正在尝试合并shapefile和.csv文件来制作选举结果的地图。我可以在加载shapefile时绘制图形，但是一旦我将其与.csv合并，它就会显示“plot.window(...)中的错误:需要有限的'xlim‘值”。
我一直在网上阅读，我想也许我必须将csv文件合并到shapefile (我一直在将shapefile合并到csv文件)。但是，csv文件(包含选举结果)比shapefile (提供地区坐标)具有更多的值。如何创建更多的地区来匹配选举结果？这会解决我的问题吗，或者我还遗漏了什么？此外，数据是西班牙语的，但相关的值Distritos=districts，partidos
浏览 7提问于2019-01-16得票数 1
我在readShapePoly包中使用maptools读取了一个shapefile，但是不能用readOGR读取相同的文件。我希望有人能帮助我阅读readOGR的shapefile。
我从这里下载了文件orcounty.shp：
我还下载了相关文件：orcounty.shx、orcounty.sbx、orcounty.sbn和orcounty.dbf，并将所有五个文件放在文件夹中：c:/users/mark w miller/gis_in_R/shapefile_example/
下面的代码读取shapefile并显示一些属性：
library(maptools)
setwd('c:
浏览 4提问于2014-03-09得票数 14
回答已采纳
2回答
我必须在R中阅读一个巨大的ESRI shapefile中的一小部分，我要分两个步骤：
步骤1:我使用ogr2ogr将shapefile剪辑到我的边界框中：
ogr2ogr -clipsrc xMin yMin xMax yMax outfile.shp infile.shp
步骤2:我用rgdal将其读入R中
df = readOGR(dsn="/path", layer="outfile")
问题是，我必须对多个文件这样做，而且很难跟踪生成每个单独文件的OGR操作。有没有办法在R中输送ogr2ogr，所以第1步是在飞行中完成的？
浏览 1提问于2013-09-10得票数 2
回答已采纳
我正在从Arcmap导入一个多边形shapefile，该文件已经有一个投影集，并附带了它的所有文件(sbn、sbx、prj等)。
但是，在使用readShapePoly函数之后，当我进行总结时，投影信息看起来是空白的。该预测是否已被识别？
Object of class SpatialPolygonsDataFrame
Coordinates:
min max
x 35551.4 1585917
y 6318047.3 9408727
Is projected: NA
proj4string : [NA]
Data attributes:
我知道有一个proj4s
浏览 2提问于2013-10-13得票数 0
1回答
我有一个闪亮的应用程序，它包括到Server中的空间表(包含从ESRI shapefile上传的数据)的以下连接：
dsn <- paste0("MSSQL:server=host\\instance;",
"database=database;",
"UID=my_username;",
"PWD=my_password;",
"trusted_connection=no")
spdf <
浏览 0提问于2018-04-06得票数 0
回答已采纳
2回答
我想在特定的位置提取一个shapefile值。我使用的shapefile可以找到，点击Download IHO Sea Areas下载。形状文件包含所有可能的海洋。
我可以阅读它并用以下方法来绘制它：
require("maptools")
require(rgdal)
require(sp)
ogrListLayers("World_Seas.shp")
shape <- readOGR("World_Seas.shp", layer="World_Seas")
但是，我想提取特定位置的海价值，比如说
p <-
浏览 2提问于2015-12-14得票数 5
回答已采纳
1回答
我在R 4.0.4 - sp v 1.4-5，rgdal 1.5-23。我使用rgdal包加载了一个shapefile，然后尝试简单地应用spTransform()，但我得到了一组无法摆脱的荒谬的警告。有人对此有什么见解吗？
options("rgdal_show_exportToProj4_warnings"="none")
library(rgdal)
library(sp)
## Shapefile has a projection of WGS84
X <- readOGR(dsn=".", layer="myshapef
浏览 154提问于2021-03-01得票数 2
我已经使用readOGR导入了一个shapefile (从包'rgdal')，并获得了一个SpatialPolygonsDataFrame。当我使用'rasterize‘函数(包’raster‘的)时，我得到了
但是我只想栅格化边缘，这样我就可以获得如下所示的GeoTiff
浏览 0提问于2013-02-26得票数 0
回答已采纳
我有一个shapefile (多边形)，我只想绘制它的一部分，所以使用ggplot的geom_path，我用xlim和ylim限制了可见路径。Shapefile是。
library(rgdal)
library(ggplot2)
coast <- readOGR(dsn=".", layer="coast-rgdal")
coast.df <- fortify(coast, region="STATE_NAME")
p <- ggplot(coast.df, aes(x=long, y=lat, group=group))
浏览 0提问于2013-09-24得票数 4
回答已采纳
1回答
我需要在R中的一个项目中处理几个shapefiles。
因此，我想迭代readOGR函数，将这些文件加载到列表中。稍后，我可以使用plyr或构造一个循环来对列表中的每个shapefile执行相同的操作。
下面是一个最低可重现性的例子：
library("rgdal")
setwd("your.path.here")
download.file("http://www.naturalearthdata.com/http//www.naturalearthdata.com/download/10m/cultural/ne_10m_admin_0_count
浏览 4提问于2014-03-27得票数 1
回答已采纳
1回答
我有一个.shp文件，我想更改它的crs，我尝试过使用spTransform，但在我的情况下它不起作用。.shp文件可以在找到。
library(rgdal)
shpfile <- readOGR(dsn="D:/Map",layer = "gpr_000b11a_e")
crs(shpfile)
CRS arguments:
+proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0
spTransform(shpfile, CRS("+proj=longlat +datu
浏览 0提问于2018-05-04得票数 0
回答已采纳
我试图读入我的形状文件，然后绘制它，但似乎RStudio卡在某些东西上，无法完成打印函数的运行。现在我有： library(rgdal)
new_county_path <- paste(county_path, "tl_2014_us_county.shp",sep='/')
county1 <- readOGR(new_county_path)
plot(county1) 但它不会产生情节，而且它似乎一直被卡住了，因为它只说 >plot(county1) 在R控制台中。我是不是做错了什么导致了这种情况，有没有更好的方法来读取和绘制shap
浏览 17提问于2021-01-09得票数 1
1回答
我试着用rgdal库读取一个shapefile，但是没有运气。
当我尝试使用以下语法导入时：
geo <- readOGR("/path/to/layer","layer")
我遇到了错误
Error in stopifnot(is.list(srl)) : infinite label point
我如何去诊断和解决这个问题？非常感谢。
浏览 2提问于2015-03-10得票数 0
1回答
我是R的地理信息系统地图的新手，我正在尝试绘制越南地图，并用犯罪数字或其他东西“填充”它。我的代码如下
st <- c('HCMC', 'Hanoi','Danang', 'Thanh Hoa', 'An Giang')
lat <- c(10.78, 21.02, 16.06, 19.82, 10.38)
long <- c(106.69, 105.83, 108.24, 105.79, 105.42)
ba <- c(96,99,88,56,7)
bls <- c(0.85, 1.
浏览 0提问于2018-10-10得票数 0
1回答
我正在使用多边形shapefile从栅格中提取平均值。栅格和矢量的CRS是相同的，但提取函数也会返回NA值。
library(rgdal)
library(raster)
options(stringsAsFactors = FALSE)
Shapefile <- readOGR("Fishnet_geolocated.shp",layer="Fishnet_geolocated")
plot(Shapefile)
CHM_Napi<-raster("Napi_CRS.tif")
crs(Shapefile)
crs(CHM_Napi
浏览 16提问于2020-05-13得票数 0
回答已采纳
1回答
所以我有一些已经可以工作的R代码。这在方法中针对shapefile从数据和空间子集中提取了一堆点。
#data is a .csv file with lon lat points
data_points <- SpatialPoints(data)
proj4string(data_points) <- CRS("+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0")
data_ll <- spTransform(data_points, CRS("+proj=longlat
浏览 5提问于2015-10-10得票数 1
有没有办法在python中做类似的事情，就像我在R中做的那样？
library(maps)
library(rgdal)
mp = 'gadm41_USA_shp' #from gadm.org
mp = rgdal::readOGR(mp)
mymap = maps::SpatialPolygons2map(mp, namefield="ISO_1")
map.where(mymap, y = 40.725959, x = -73.978544)
#[1] "US-NY:88"
链接到shapefile：
在python中，我可以在shape文
浏览 9提问于2022-10-28得票数 0
1回答
我在栅格化shapefile以在0.5*0.5网格上生成点时遇到了问题。shapefile表示全球珊瑚礁的风险级别(低-0，中-100，高-1000，极高-1500)到综合威胁的分类。
我从另一个运行良好的示例中提取了代码，但是当我对我的数据尝试它时，我从plot函数中得不到任何东西。有关shapefile和我的代码的链接，请参见以下内容：
# Read shapefile into R
library(rgdal)
library(raster)
int.threat.2030 <- readOGR(dsn = "Global_Threats/Integrated_
浏览 7提问于2016-01-30得票数 3
回答已采纳
1回答
是否可以使用ggplot2 / ggmap在由shapefile定义的给定空间区域内随机绘制geom_points？
我考虑了geom_jitter，但是我需要这些地块是随机分布的，而不是跨越空间边界。
示例数据无耻地从@matthiash 借来。
library(rgdal)
library(ggmap)
# Get shapefile with Drammen municipality borders
tmpzip<-tempfile()
tmpdir<-tempfile()
dir.create(tmpdir)
download.file("http://www.k
浏览 3提问于2017-11-15得票数 0
回答已采纳
我正在使用{ raster }根据不规则的shapefile ( Amazon biome)来裁剪(或裁剪)栅格，但输出总是有一个矩形范围。但是，我需要在与shapefile完全相同的几何图形中输出。有什么建议吗？干杯。
library(raster)
library(rgdal)
myshp <- readOGR("Amazon.shp", layer="Amazon")
e <- extent(myshp)
myraster <- raster("Temperature.tif")
myraster.crop <-
浏览 0提问于2014-04-15得票数 24
回答已采纳
1回答
Spatstat中的样方
、、、
我现在正在做探索性分析，目标是绘制一个样本地图，确定是否存在完全的空间随机性(视觉上和使用卡方检验)等。然而，我在绘制样本地图时遇到了问题。我之前问过，这导致我重新投影了我的数据。
下面是我的代码：
library(rgdal) #Brings Spatial Data in R
library(spatstat) # Spatial Statistics
library(maptools)
library(raster)
# Load nyc zip code boundary polygon shapefile
s <- readOGR("/Users/my_name/D
浏览 3提问于2020-12-15得票数 0
如何显示在zip文件中的shapefile格式的地图。文件如下：
我试着做以下几件事：
library(rgdal)
temp <- tempfile()
download.file("https://geoftp.ibge.gov.br/organizacao_do_territorio/malhas_territoriais/malhas_municipais/municipio_2015/UFs/PR/pr_municipios.zip",temp)
data <- readOGR(unz(temp, "41MUE250GC_SIR.shp"
浏览 13提问于2022-06-05得票数 0
回答已采纳
我已经加载了多个包/尝试过许多读取.shp的方法。所有的方法都不起作用，我无法破解谷歌相关的建议。运行R中的代码，系统冻结，RStudio.cloud中相同的代码只返回各种错误消息。包和代码的列表如下。最终，我试图创建映射，加入另一个文件并生成一个结果表。如有任何建议，将不胜感激。
library(rgeos)
library(sp)
library(rgdal)
library(maptools)
library(sf)
require(sf)
area <- read_sf(dsn = ".", layer = "2020ATLclip")
req
浏览 16提问于2022-01-30得票数 0
1回答
我一直试图用hcmap从highcharter包中绘制一个choropleth地图；我从我自己的shapefile中获得了多边形，因为它是一个不在highmap收藏列表中的地图。
为此，我首先设法将shapefile转换为GeoJson文件，如下所述：
后来，我成功地使用包geosonio绘制了地图，如下所述：
但是，我不知道如何将数据和值合并到我的地图中绘制的多边形中。所有可用的示例都合并为data.frame格式的映射数据，在转换为GeoJson时我会丢失该格式。
到目前为止，这是我的代码：
library(rgdal)
library(geojsonio)
library(highcha
浏览 0提问于2019-06-17得票数 1
回答已采纳
我试图使用R包tmap绘制一个SpatialLinesDataFrame。然而，我一直收到这样的错误消息：
CPL_geos_is_empty中的错误(st_geometry(X))：计算错误: IllegalArgumentException:点阵必须包含0或>1个元素。
我第一次尝试用tmap绘制"sp“对象。然后，我将sp对象转换为"sf“对象，仍然得到相同的错误消息。我做了一点谷歌搜索，我认为这可能与SpatialLinesDataFrame中的线路没有完全连接的事实有关。
所以我给出了一组连在一起的线，它用tmap绘制得很好。
我把这里使用的形状文件放在
浏览 1提问于2019-08-28得票数 3
1回答
我试图将一个shapefile应用到ggmaps地图上，但它给了我非常奇怪的结果。讨论中的shapefile是澳大利亚统计局提供的中的“统计局部区域”(与邮政编码相似的组) shapefile。
通常，我可能会认为这是一个切断边缘点的问题，但是即使在缩放级别1(实际上它看起来更糟)，我也在点击它：
下面是我用来生成上面图表的一些代码：
library(tidyverse)
library(ggmap)
library(rgdal)
slas <- readOGR(dsn="SLA",layer="SLA11aAust")
aus4 &l
浏览 3提问于2017-10-05得票数 3
回答已采纳
2回答
在执行任务时循环遍历文件夹
、、、、
在读取R中的变量并将文件分配给变量时，我试图遍历目录中的文件夹。
我要分配给变量的文件是shapefiles，所以我使用来自rgdal包的函数readOGR。目的是以后合并属于特定物种的所有形状文件。目录的层次结构/结构为type1 > species > ids。shapefiles看起来像id.shp等等。
示例shapefiles可以下载和
#code
setwd("~/type1/")
#Extract ids belonging to $species. Later use in readOCR function
sp_id <- functio
浏览 2提问于2013-12-11得票数 0
回答已采纳
我正在尝试让R读取简单景观分类的shapefile (发布的)，但收到错误
shape1 <-readShapePoly("JustForest_BahiaSul_1990_Test.shp", force_ring=TRUE, delete_null_obj=TRUE)
Error in Polygon(coords = crds) : ring not closed
当我尝试使用rgdal作为替代方案时，也会出现同样的错误：
shape1 <-readOGR("JustForest_BahiaSul_1990_Test.shp", layer=&
浏览 5提问于2013-05-07得票数 0
1回答
我一直在使用ArcMap访问空间数据服务器上的地理信息系统数据。我想知道如何在R.
我知道如何将shapefiles读入R。我已经成功地使用了maptools和rgdal来打开和映射本地存储的shapefile(例如)。
我的问题是数据不是本地存储的，而是存储在Application上的。我相信这是一个Oracle数据库。我已经得到了关于1.server 2.实例(一个数字)3.数据库4.用户和5.密码的信息。通常，我会包括一个例子，但是外部用户是否能够访问服务器是值得怀疑的。
例如，下面是如何读取和绘制R中的本地文件
library(rgdal)
ogrInfo(".
浏览 2提问于2014-04-22得票数 2
回答已采纳
我正在尝试使用ggplot来绘制shapefile ()的多边形，以便可视化一些包含在其中的区域的数据。下面的代码会产生一个错误：
library(rgdal)
library(tidyverse)
library(maptools)
hsa <- readOGR(dsn = "Data/hsa_bdry", layer = "HSA_Bdry") %>% fortify(region = "HSA93")
ggplot() + geom_polygon(data = hsa, aes(x = long, y = lat, gro
浏览 2提问于2017-09-01得票数 0
回答已采纳
1回答
类似于state/territory，在R中是否有生成state/territory合唱团的本地方式？
我可以看到一些有文档记录的方法，但没有一个是完整的(例如，看起来有点有希望，但缺乏原始数据)
library(rgdal)
library(spdplyr)
library(geojsonio)
library(rmapshaper)
# Load Australian State and Territories shapefile data
aus_ste <- readOGR(dsn = "/Users/kannishida/Downloads/STE11aAust"
浏览 6提问于2020-04-15得票数 1
我有一个NetCDF文件，它有38年的降水数据。我使用我的区域的shapefile来提取我感兴趣的区域的数据。当我使用下面的代码时，我会得到与多边形内的网格点对应的数据。但是，我感兴趣的是根据shapefile的每个多边形提取数据(在本例中为12)，并使用多边形的名称保存数据帧。下面是我的示例代码
library(ncdf4)
library(rgdal)
library(raster)
library(maptools)
library(GISTools)
NC = brick("Daily_Pcp.nc")
Subbasin_SHP=readOGR("my_Sha
浏览 2提问于2019-06-25得票数 3
2回答
绘制UTM /坐标
、、、
我是R编程新手，我有一个csv/excel文件，其中包含了一个国家的20个城镇，其中包含以下格式，
城镇UTM协调UTM坐标
xxxxxxxx 1377777 249514
yyyyyyyy 142145 228942
我无法将它们绘制成地图，是否有人知道如何绘制这些UTM坐标。是否有可能用UTM在R编程中规划城镇？如果是的话有人能帮我一下吗。我也有这个国家的形状文件。但我不知道该怎么处理。
myfilepath <- file.choose()
Cordinates<-read.csv(myfilepath,header=TRUE)
Cordinates
str(Cordinat
浏览 6提问于2021-03-27得票数 2
回答已采纳
我正在尝试获得一块北极有冰存在的地块。我使用的是从获得的shapefile
library(maptools)
library(rgeos)
library(ggplot2)
library(rgdal)
library(plyr)
library(ggmap)
library(dismo)
library(raster)
library(plyr)
ICE <- readOGR(".", "masie_ice_r00_v01_2018094_4km")
plot(ICE)
ICE_WGS84 <- spTransform(ICE,
浏览 2提问于2018-04-06得票数 3
我试图从一个简单的shapefile加载并绘制一个图层。这是美国的国家边界。我可以把它装好：
> library("sp","rgdal")
> shape = readOGR("/home/username/data/share/mapnik/world_boundaries/states_EPSG4326.shp", layer="states_EPSG4326")
OGR data source with driver: ESRI Shapefile
Source: "/home/u
浏览 1提问于2012-12-02得票数 4
回答已采纳
1回答
在R中读取和绘制esri形状文件
、、、、
我在读.shp (esri形状文件)到R中有困难，我已经尝试了R中的几个选项，并试图将ArcMap中的形状文件转换为可以正确读取形状文件的内容，但都还没有起作用。(在ArcMap中，我更正了几何体，从单一多边形转换为多多边形，等等，这可能是不必要的或不相关的)
这可能与我的形状文件包含“区域”(多多边形)而不是“多边形”有关。
如何在R中正确读取图形文件类型以进行绘图？(它看起来像ArcMap中的普通形状)
在ArcMap中，形状文件如下所示：
(显示包含其他多边形中的多边形的形状文件)
在R中，它看起来像这样：
(显示形状文件，其中某些多边形填充错误)
我在R中使用了以下代码：
req
浏览 3提问于2010-01-15得票数 15
回答已采纳
我正在尝试使用sf和RGDAL包为一个项目创建一个掩码。我有一个shapefile，它是一个SpatiallinesDataFrame，它需要是一个SpatialpolygonsDataFrame才能继续。 我试过这些代码 fence1 <- as(fence1 ,"SpatialPolygonsDataFrame")
fence1 <- SpatialPolygonsDataFrame(fence1)
fence1 <- SpatialPolygonsDataFrame(fence1, match.ID = TRUE) fence1
浏览 16提问于2019-09-24得票数 0
有什么简单的方法可以在Windows XP上访问R中的PostGis数据吗？
我想我可以尝试rgdal，但没有简单的方法来添加PG驱动程序(参见file.show(system.file("README.windows", package="rgdal")) )。因此，这段简单的代码无法工作：
library(rgdal)
mylayer <-readOGR(dsn="PG:host=localhost user=MyUser dbname=MyDb password=Secret port=5432", layer = "MyLa
浏览 0提问于2013-03-07得票数 3
回答已采纳
我试图用tmap库绘制一张地图：
library("rgdal")
library("tmap")
area_dnipro <- readOGR("dnipro.shp", encoding = "UTF-8", use_iconv=TRUE)
tm_shape(area_dnipro) + tm_polygons()
但在策划过程中我得到了这样的信息：
Warning message:
In sp::proj4string(obj) : CRS object has comment, which is lost in ou
浏览 3提问于2021-04-21得票数 1
回答已采纳
2回答
我(在Win7上)试图获取我的shapefile数据(是数据文件)。使用leaflet包显示。但没有任何成功。我只得到背景地图瓦片，而不是我的数据。不过，我可以在QGIS中看到我的数据。我从得到了下面的代码和描述。
library(rgdal)
shapeData <- readOGR(".",'myGIS')
ogrInfo(".",'myGIS')
leaflet() %>% addTiles() %>% setView(lng = -106.363590, lat=31.968483,zoom=11) %
浏览 0提问于2015-10-10得票数 7
回答已采纳
3回答
我试图获得一组点的坐标，定义多边形内的网格(我有一个shapefile )。看起来最简单的方法就是创建一个点的网格，然后将这些点过滤到多边形内的点。我看了一下和，根据答案，我尝试了这样的方法：
library(rgdal)
city_bdry <- readOGR("Boundaries - City",
"geo_export_32ded882-2eab-4eaa-b9da-a18889600a40")
res <- 0.01
bb <- bbox(city_bdry)
gt <- GridTo
浏览 5提问于2017-09-05得票数 3
回答已采纳
1回答
The
密码
#
# This is code for mapping of CGE_Morocco results
#
# rm(list = ls(all = TRUE)) # don't use this in code that others will copy/paste
## Loading packages
library(rgdal)
library(plyr)
library(maps)
library(maptools)
library(mapdata)
library(g
浏览 0提问于2013-11-23得票数 2
回答已采纳
我有一个shapefile，http://census.cso.ie/censusasp/saps/boundaries/Census2011_Small_Areas_generalised20m.zip 并希望提取经度/纬度，但我不确定如何将正确的坐标映射到正确的小区域。 mycode是： require(ggplot2)
require(proj4)
require(rgdal)
a=readOGR(....shp)
dublin = a[a$NUS3NAME=='Dublin',]
dublin=spTransform(dublin,CRS('=proj=lo
浏览 18提问于2019-05-25得票数 0
1回答
我试图通过使用sf和rgdal软件包指导同事在R中使用Shapefile。我提供的示例代码：
dz_boundaries <- readOGR(dsn="~/Downloads/SG_DataZoneBdry_2011", layer="SG_DataZone_Bdry_2011")
example示例是从：下载的
我能够读取这个格式文件，并在家用Mac和Windows机器上操作它，而不存在任何问题。但是，我的同事在试图使用Mac遵循代码时遇到了错误：
Error in readOGR: no features found
ogrInfo:all feat
浏览 5提问于2020-07-07得票数 1
回答已采纳
1回答
我目前正在处理地理数据，以便绘制它们。我有两个不同的shapefile。我希望做的是在加拿大和美国之间增加一个大胆的边界，但我不知道如何做到这一点。
# Load packages ----------------------------------------------------------------
library(dplyr)
library(readr)
library(ggplot2)
library(rgdal)
library(ggmap)
# Load data ------------------------------------------------------
浏览 0提问于2017-04-12得票数 3
1回答
我正试图用传单来画许多多边形，但我不明白出了什么问题。
我使用的shapefile可以在这里找到：
library(leaflet)
library(rgdal)
df <- readOGR("C:/Users/me/codes_postaux","codes_postaux_region")
plot(df)
在我看来，shapefile似乎没问题，我使用的代码也相当简单。然而，我只得到了一个输出的地图，没有多边形。我已经为这个问题挣扎了很长一段时间，如果有人能在这里帮助我，我会非常感激的。
map <- leaflet(df) %>
浏览 0提问于2016-11-02得票数 3
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
云直播
对象存储活动推荐
腾讯云开发者
