带有项目的R警告中的spTransform()
在R语言中，带有项目的R警告中的spTransform()是一个函数，用于在空间数据分析中进行坐标系转换。它是sp包中的一个函数，该包提供了处理空间数据的功能。
spTransform()函数的作用是将空间数据对象的坐标系从一个参考系统转换为另一个参考系统。它接受两个参数：sp对象和目标坐标系。sp对象可以是SpatialPoints、SpatialPointsDataFrame、SpatialPolygons、SpatialPolygonsDataFrame等类型的对象。
该函数的优势在于可以方便地将空间数据对象从一个坐标系转换为另一个坐标系，以满足不同分析需求。它可以用于各种空间数据分析任务，如地理信息系统（GIS）、地图制作、空间统计分析等。
在腾讯云的产品中，与空间数据分析相关的产品是地理大数据分析平台（GeoDT）。GeoDT提供了丰富的地理数据处理和分析功能，包括坐标系转换、空间数据查询、空间分析等。您可以通过以下链接了解更多关于GeoDT的信息：https://cloud.tencent.com/product/geodt
需要注意的是，本回答仅提供了腾讯云的相关产品作为参考，其他云计算品牌商也提供类似的产品和服务。
相关·内容
1回答
我有一张传单地图，我在上面画了几个多边形。我正在尝试使用rmapshaper::ms_simplify()来减少地图的加载时间。 polygons_ll <- (
spTransform(polygons, CRS("+proj=longlat +datum=WGS84"))
)
polygons_ll <- ms_simplify(polygons_ll) 这一切都在一个R闪亮的应用程序中，它在遇到警告后不会运行： Warning in sp::proj4string(sp) :
CRS object has comme
浏览 50提问于2020-12-08得票数 0
回答已采纳
1回答
是否有办法在保留R的同时改变LiDAR点云的投影(比如sf包中的LiDAR)？
使用lidR包，您可以设置但不能修改点云投影。我知道来自lastools的las2las可以完成这项工作(并且这个部分是开放的)，但是我希望保持它简单，在R内，并且不需要很难地复制这个过程的输入/输出。
有什么建议吗？
浏览 5提问于2020-09-26得票数 0
回答已采纳
1回答
我从level4下载了TDWG的形状。这个形状有一个WGS84投影。我需要将这一层重新投影到欧洲的LAEA或CRS("+init=epsg:3035")。
我使用了以下代码：
TDWG4 <- readOGR("D:/GIS/Administrative/world/TDWG/level4/level4.shp", layer="level4")
TDWG4.LAEA <- spTransform(TDWG4, CRS("+init=epsg:3035"))
并得到以下错误：
.spTransform_Polygon
浏览 6提问于2017-09-25得票数 1
回答已采纳
我从我的同事那里得到了一个地图图形文件，我还得到了一些观测数据，他们的坐标记录在WGS84中。我读取了R中的shapefile：
# read shapefile
> shpmap <- readOGR(dsn = "soil_type", layer = "soil_type_data")
OGR data source with driver: ESRI Shapefile
Source: "soil_type", layer: "soil_type_data"
with 18542 features
It ha
浏览 0提问于2015-10-13得票数 1
我在R中加载了SpatialPointsDataFrame对象，范围以米为单位。我要他们转换成度小数。请帮帮我怎么做。
下面是我的数据在控制台中的样子：
class : SpatialPointsDataFrame
features : 23
extent : 351912.5, 457807, 3236835, 3367232 (xmin, xmax, ymin, ymax)
crs : +proj=utm +zone=44 +datum=WGS84 +units=m +no_defs
variables : 1
names
浏览 8提问于2022-05-15得票数 0
我是R中的GIS新手，我试图在shapefile中添加lat/long或UTM坐标。我从这里下载了芝加哥市(City_Boundaries.shp)的边界：
我加载了maptools和rgeos库：
library(maptools)
library(rgeos)
library(rgdal)
我将数据输入R&试图为区域16T添加一个UTM代码：
city<- readShapeSpatial("C:/Users/Luke.Gerdes/Documents/GIS Files/Chicago/City_Boundary.shp")
proj4string(city
浏览 3提问于2015-02-25得票数 1
回答已采纳
我正试图创建一个选定的加拿大省/地区和选定的美国州的地图。到目前为止，最好的地图似乎是用GADM数据生成的地图：。
然而，我无法在同一张地图上绘制美国和加拿大的地图，也无法只绘制选定的省/地区和州。例如，我对阿拉斯加、育空、西北、不列颠哥伦比亚省、艾伯塔和蒙大拿等地感兴趣。
此外，美国地图似乎是沿着国际日期线划分的。
有谁能帮我：
datelineoverlay 将上述省份/地区和州绘制在一张地图上，避免美国沿着国际地理区域的一个纬度-经度网格选择一个特定的投影，可能是多边形。。
也许spplot不允许用户指定投影。我没有看到在see页面上选择投影的选项。我知道如何在地图包中选择带有映射函数的
浏览 4提问于2012-05-26得票数 12
回答已采纳
首先，我是新来的R，所以请容忍我。
我最终要完成的是展示一张带有传单的阿姆斯特丹交互式地图。为此，我使用RGDAL读取shapefiles。
这个包含阿姆斯特丹的shapefiles。
我使用下面的代码来读取shapefiles并显示地图。
amsterdam <- readOGR(".", layer = "sd2010zw_region", verbose = FALSE)
leaflet(amsterdam) %>%
addProviderTiles("CartoDB.Positron", options= provide
浏览 5提问于2016-06-01得票数 7
回答已采纳
我在R中遇到了一些我不明白的投影问题。
我下载了以下全局数据集：
然后，我用不同的投影来制作地图，教授空间投影的概念。
我成功地映射了lat/lon WGS84中的数据，并将其重新投影/映射到了Robinson
library(rgdal)
library(ggplot2)
setwd("~/Documents/data")
# read shapefile
worldBound <- readOGR(dsn="Global/Boundaries/ne_110m_land",
layer="ne_110m_la
浏览 6提问于2016-02-24得票数 2
回答已采纳
1回答
我对R中的映射函数很陌生。我正试图在地图上绘制点，然而，我需要他们在学位(他们目前在脚)。我试图找出坐标转换为度的代码，有人知道这个坐标的代码或公式吗？
浏览 0提问于2017-03-22得票数 0
2回答
从Python运行R脚本
、、、、
我想从Python脚本运行R脚本。在不同坐标系中投影经度坐标需要R脚本。我考虑了两个选项来做到这一点。在第一个选项中，我喜欢将经度和经度坐标解析为R脚本，如下所示。最后，我希望R脚本将x和y返回给python脚本，但我不知道如何做到这一点。
project<-function(lat,lon){
library(sp)
library(rgdal)
xy <- cbind(x = lon, y = lat)
S <- SpatialPoints(xy)
proj4string(S) <- CRS("+proj=longlat +ellps=WGS84 +da
浏览 0提问于2013-05-09得票数 4
回答已采纳
昨天很好，今天当我清理或构建项目时，它显示了这个警告：
警告:忽略Android伪影com.google.android:android:2.1_r1用于appDebug警告:忽略Android伪影com.google.android:android:2.1_r1用于appRelease
我认为这是因为升级了android的gradle插件。
有人知道这意味着什么吗？
浏览 1提问于2016-05-19得票数 32
我有使用纬度和经度指定的多边形数据，以及在向东和向北时的点数据。我的问题是如何将东向和北纬数据转换为经度和纬度点来绘制多边形上的图？我的数据如下：
feature_easting feature_northing
<dbl> <dbl>
1 209420 922820
2 260123 666308
3 337004 795573
4 259210 723869
浏览 2提问于2018-06-06得票数 0
回答已采纳
1回答
空间数据自动克里格
、、、
我试图用克里格函数来创建海洋剖面中化学参数的垂直图，我现在很难开始工作。
我的数据如下所示：
horiz=rep(1:5, 5)
depth=runif(25)
value = horiz+runif(25)/5
df <- data.frame(horiz, depth, value)
autoKrige包中的automap函数看起来应该为我完成这项工作，但它需要一个SpatialPointsDataFrame类的对象。据我所知，包rgdal中的函数rgdal创建了SpatialPointsDataFrame对象，但是有两个问题：
来自CRAN的这个的OSX二进制文件，我在OXS1
浏览 2提问于2012-06-13得票数 3
回答已采纳
我正在尝试使用库(Tmap)在R中创建一张带有全色RGB蒙版陆地卫星图像的地图。然而，NAs显示为黑色。这是我所做的。
使用库(Sf)，我计算了5个多边形的质心，并将它们缓冲了5000米。之后，我使用库(光栅)通过缓冲质心来遮罩陆地卫星图像。代码看起来像这样，并且工作得很好。
# Read in the data
polys <- st_read("nybb.shp")
rast <- brick("LC08_L1TP_013032_20171027_20171027_01_RT.tif")
# Transform polys, crop rast
浏览 1提问于2017-11-03得票数 2
我有一些在SWEREF 99TM中协调的，我不能使用R中的proj4包转换成WGS84。
temp <- data.frame (x = c(598223, 598812, 598824, 598232, 597614, 597629),
y = c(7095460, 7095426, 7094827, 7094227, 7094821, 7095433))
ptransform(temp,
src.proj = '+proj=utm +zone=33 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +n
浏览 6提问于2019-04-09得票数 0
1回答
我在Windows8.1中安装rgdal时遇到了问题，我安装了OSGeo4W (这是强制性的吗？)，我明白我必须安装proj4，但不知道在上安装什么，我读了一些文档，但是我觉得很无知，无法让它正常工作。
我想通过R-studio通过rgdal访问R-studio，但是如何告诉r-studio在哪里查找gdal和proj4？我真的很困惑，显然需要一些文档，但不知道到哪里去找，因为我所读到的东西似乎都没有用。
rgdal是根据R安装的，但是当我使用spTransform时，它说
Error in spTransform(locations, CRS("+proj=longlat +ellp
浏览 4提问于2014-06-29得票数 2
回答已采纳
我在Behrmann等面积圆柱投影中有一个200 in x 200 in的网格单元的shapefile，覆盖了世界上的陆地区域。我的目标是将shapefile转换为+proj=longlat WGS84格式，这样我就可以将其与常用投影中的地图进行匹配，例如地图工具中的wrld_simpl。然而，我并没有成功，希望在这方面能得到一些帮助。 rm(list = ls())
library(RCurl)
library(raster)
library(maptools)
library(rgdal)
data("wrld_simpl")
tmp <- tempfile()
浏览 23提问于2019-01-19得票数 0
1回答
我对GIS和R还比较陌生，我正在做一个项目，把犯罪数据映射到城市的特定社区。数据文件采用CSV格式，具有事件ID、经度和纬度。邻域多边形由shapefile定义。以下是我在R中的代码：
require(sp)
require(rgdal)
require(maps)
crime<-read.csv("crime sample.csv")
crime<-subset(crime,!is.na(crime$Latitude))
coordinates(crime)<-c("Longitude","Latitude")
Neig
浏览 1提问于2013-11-18得票数 2
回答已采纳
1回答
我想知道是否有人知道如何在R中实现地图和栅格的平等地球投影？
我找到了这个教程：
但我不能复制它。看起来"+proj=eqearth“是不可用的，尽管它现在看起来应该是可用的。
library(rgdal)
library(maps)
library(maptools)
RobinsonProj="+proj=robin +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs"
EqualEarthProj= "+proj=eqearth +lon_0=0 +x_0=0 +y_0
浏览 6提问于2020-04-21得票数 0
我用这种方式创建一个SpatialPolygon
#make spatial points
#assign Original CRS WGS84 EPSG:4326 (DATUM --> on sphere)
cornersEPSG4326 <- SpatialPoints(coords=cbind(x,y), proj4string = CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs"))
#transform to EPSG3857 (Web Mercator PROJECTION)
cornersE
浏览 1提问于2015-09-03得票数 1
回答已采纳
我刚刚开始从事一个旧项目，在这个项目中有许多Xcode警告。
其中一个警告是R.swift警告："R.swift project.pbxproj在内部不一致。“
在构建目标-日志中，它说：“警告: R.swift project.pbxproj在内部不一致。
不使用XCRemoteSwiftPackageReference (FF0493CB254AF9C600793461)
不使用XCSwiftPackageProductDependency (FF0493CC254AF9C600793461)
不使用XCRemoteSwiftPackageReference (FF04
浏览 11提问于2022-08-31得票数 0
回答已采纳
我想预测2012年秘鲁人口普查的坐标，以便与另一个具有下列坐标参考系统的光栅层一起工作：
"+proj=utm +zone=18 +south +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs“。
帖子的答案在我的情况下行不通。
crs(census)
CRS参数：+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0
newCRScensus=CRS("+proj=utm +zone=18 +south +ellps=WGS84 +
浏览 0提问于2018-03-26得票数 1
回答已采纳
1回答
我有一个空间多边形对象和一个空间点对象。后者是由xy数据向量(分别称为纬度和经度)在数据中创建的，而前者则直接使用rgdal读取到R中。我的代码如下：
boros <- readOGR(dsn = ".", "nybb")
rats <- read.csv("nycrats_missing_latlong_removed_4.2.14.csv", header = TRUE)
coordinates(rats) <- ~longitude + latitude
此时，两个空间物体都没有投影。如果我按以下方式投影这些对象：
pro
浏览 4提问于2014-04-03得票数 1
回答已采纳
我在Intellij Idea 13中创建了一个Griffon项目，Griffon版本是1.5.0。在整个项目中，我可以看到许多带有下划线的关键字。当我将光标移到下划线上时，下面的文字会显示为警告“无法解析符号”。当我运行应用程序时，它工作得很好。请帮助我解决这个错误。有关详细信息，请参阅下图。
浏览 1提问于2014-09-03得票数 0
回答已采纳
在Android Studio版本3.1中，依赖编译word被替换为实现。android studio 3.1中带有警告的依赖项：
dependencies {
compile fileTree(dir: 'libs', include: ['*.jar'])
compile 'com.android.support:appcompat-v7:27.1.0'
compile 'com.android.support.constraint:constraint-layout:1.0.2'
浏览 0提问于2018-04-07得票数 1
我正在使用Sentinel 2图像数据进行NDVI监测。
我在R中有一个导入为SpatialPolygonsDataFrame的区域的村级边界。当我试图裁剪出感兴趣的区域时，当我到达NDVI光栅层后，它显示范围不重叠。
我在我的NDVI layer(utm)中发现我的区域变量的CRS是different(longlat)。
> crs(ndvi)
CRS arguments:
+proj=utm +zone=42 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0
> crs(region)
CRS argume
浏览 1提问于2018-08-22得票数 0
1回答
我想在R中建立一个以太平洋为中心的罗宾逊投影.
使用spTransform将shapefile转换为罗宾逊CRS并将经度设置为150度是行不通的，fortify错误也不起作用。
就像这样：
world_robin <- spTransform(world2, CRS("+proj=robin +lon_0=150 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs"))
plot(world_robin, col = "khaki", bg = "azure2")
在稍后的过程
浏览 0提问于2015-09-15得票数 1
我试图在R中绘制一条路由，.kml文件来自Marble，一个虚拟的地球仪应用程序。看起来是这样的：
<?xml version="1.0" encoding="UTF-8"?>
<kml xmlns="http://earth.google.com/kml/2.2" xmlns:gx="http://www.google.com/kml/ext/2.2">
<Document>
<visibility>1</visibility>
<Folder&g
浏览 2提问于2014-01-29得票数 1
回答已采纳
1回答
编辑:我稍微改写了标题问题，并调整了文本以回应@DWin的评论。
将预测的地理层与未投影的地理层组合起来可能是一项挑战。通常情况下，由于地理层来自不同的产品和出版商，因此似乎有必要进行一些转换。
我知道R有几个工具来执行地理转换。例如：
对于sp包中的类Spatial*对象，可以使用rgdal包中的spTransform()函数；对于raster包中的Raster*类对象，可以使用projectRaster()函数。<code>H 211</code><code>G 212</code>
下面是我想在R: Transform to UTM网格区域
浏览 0提问于2012-02-10得票数 3
回答已采纳
我有一个带有复选框的可扩展列表视图。当我单击子对象时，会出现一个警告对话框，我会选择数量，然后子对象的textView就会改变。但是当我向下滚动列表时，这个孩子从视图中消失了，列表忘记了更改的文本视图，而是设置了旧的。理由是什么呢？
浏览 0提问于2015-11-12得票数 1
1回答
输入
我有两个形状文件，我输入到R，所以我结束了。一种包含公共汽车路线的空间数据。一个包含公交车站的空间点。
划出一个给定的路线和它的停站如下所示。
样本数据
此包含两个表单文件，作为一个zip下载，其中包含两个路径示例。
目标
我的目标是计算每一对车站之间的地理距离(以米为单位)：停1到停2，停2到停3，等等。
我发现的大多数计算欧几里德距离的方法，或者说是“乌鸦苍蝇”，在这里都是行不通的。
提到了PBSmapping，它有一个calcLength函数，它可以很好地计算路由的总距离，但是我无法找到一种方法来匹配它与停止对的情况，也无法找到一种方法来根据它的属性对shapefile进
浏览 0提问于2018-04-10得票数 2
我想在带有wgs84坐标的线状文件周围创建一个缓冲区。
我准备了一个带有单行段和Datum: D_WGS_1984的shapefile。之后，我使用‘.shp’命令将readOGR加载到R中。
之后，我从rgeos包中尝试了用于计算缓冲区的gBuffer方法：
gBuffer(l2, width=1.0, quadsegs=5, capStyle="ROUND", joinStyle="ROUND", mitreLimit=0.01))
Warning:
In gBuffer(l2, width = 1, quadsegs = 5, capStyle = &#
浏览 2提问于2012-03-16得票数 4
回答已采纳
1回答
我是R的新手，在为一些点模式分析设置我的数据时遇到了一些麻烦。 我想做的是:对纽约市逮捕数据进行点模式分析，看看逮捕人数和新冠肺炎案件之间是否存在空间相关性。 到目前为止，我所做的是:以shapefile的形式下载数据 https://data.cityofnewyork.us/City-Government/Borough-Boundaries/tqmj-j8zm (邮政编码边界) https://www1.nyc.gov/site/nypd/stats/crime-statistics/citywide-crime-stats.page (按邮政编码列出的纽约市逮捕人数年初至今数据) 代
浏览 18提问于2020-12-14得票数 0
我有一个拉丁美洲的.TIFF文件(栅格)和巴西市政当局的.shp文件(形状文件)，我想重叠这两个，但我想和所有的信息，一个栅格。
在我的例子中，我有一个关于光污染的TIFF文件。光栅中的每个点代表一个光信息，范围从0到63。shapefile中只有几何体。
最后，我希望市政当局的每个信息都返回所有栅格点的总和(或平均值)。
有人知道这在R软件中是可能的吗？
非常感谢
我不能举一个例子，因为TIFF文件太大了，但我发送了一个
我想覆盖像这样的地图，但尺寸较小
浏览 32提问于2019-11-01得票数 0
回答已采纳
2回答
自从R# 2017.1.1以来，我在自动生成的GetHashCode()函数中得到警告。让我解释一下这个函数是如何创建的：
如果重写类的Equals函数，R#建议您创建相等的成员。如果让R#生成那些相等的成员，它也会重写GetHashCode()函数。但是，在那里，它使用了我类的所有属性。由于这些属性并不都是只读的，R#告诉我应该让它们成为只读的，并显示一个警告。
因此，我的问题是，我是否应该将GetHashCode()函数保留为空(或者删除那些带有警告的部分)，或者是否应该使属性成为只读的，这样R#就不再警告我了。
这是我的密码：
public override bool Equals(ob
浏览 3提问于2017-04-12得票数 2
回答已采纳
我尝试使用将NAD83状态平面坐标中的一组x-y坐标转换为以度为单位的常规Lan/Lat坐标。我可以复制这篇文章中给出的例子，但它不能给我的集合正确的答案！下面是我尝试过的和我得到的错误答案。
library(rgdal)
nad83_coords <- data.frame(x=c(577430), y=c(2323270)) # My coordinates in NAD83
coordinates(nad83_coords) <- c('x', 'y')
proj4string(nad83_coords)=CRS("+init=esri
浏览 0提问于2016-09-15得票数 3
2回答
我很难理解这个查询的工作原理。它应该返回已订购所有商品的客户的姓名。
R指的是客户所做的商品订单表，其中包含客户id (cid)和项目id (iid)。
I引用可以排序的项表，其中包含项id。
C是带有customer id的客户表
SELECT cname
FROM Customer C
WHERE NOT EXISTS
( (SELECT I.iid
FROM Item I)
EXCEPT
(SELECT R.iid
FROM Order R
WHERE R.cid=C.cid))
根据我的假设，带有SELECT R.iid的底层嵌套查
浏览 1提问于2013-02-27得票数 3
回答已采纳
2回答
我有两个空间对象，一个是空间多边形对象，另一个是.csv文件，我将其转换为空间点对象。第一个是智利政府为其一个公社创建的官方形状文件，另一个是通过使用HERE API地理编码创建的，这里是同一公社的街道地址。
首先，我使用readOGR从加载了空间多边形对象：
quilpue <- readOGR( dsn= getwd() , layer="quilpue-rgdal",
encoding = "UTF-8")
然后，我将.csv文件加载到R中，并使用sp包中的coordinates()函数将其转换为空间点对象。
po
浏览 7提问于2017-05-09得票数 2
分配给我的任务是从.nc文件中从.tif文件中剪辑光栅。
编辑(从评论)：我想提取临时。来自.nc的信息，因为我需要检查特定地区的年平均气温。要进行比较，就必须在完全相同的区域上进行比较。.nc文件比以前检查过的区域要大，所以我需要将它“剪辑”到.tif的程度。.tif数据为0\1格式，其中它是0(或者.tif小于.nc)，.nc数据应该是“剪裁”的。最后，我希望保留.nc数据，但保留.tif的范围，同时仍然保留其分辨率和投影。(.tif和.nc有不同的投影和像素大小)
现在，通常情况下，这不会是一个问题，因为我可以使用raster::crop。但这并不能处理不同的投影和不同的像素大小/分辨率
浏览 1提问于2020-08-13得票数 0
1回答
我在R中使用以下代码将我的UTM坐标转换为十进制度
utmcoor<-SpatialPoints(cbind(starling2$Easting,starling2$Northing), proj4string=CRS("+proj=utm +zone=51"))
longlatcoor<-spTransform(utmcoor,CRS("+proj=longlat"))
如何将结果('latlongcoor‘输出)作为两个新列添加到原始dataframe 'starling2’中？
浏览 2提问于2017-07-21得票数 1
最近我研究了R，使之成为1km×1km的网格。所以，我已经做了下面的代码。但shapefile(korea, seoul and WGS84)无法准确地将(或覆盖)添加到网格中。所以我需要你的帮助来解决这个问题。
shapefile有类似于这个+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84+towgs84=0,0,0的CRS
为了使网格分辨率达到1kmx1km，我一直在使用spTransform函数，如下所示。seoullayer_t <- spTransform(seoullayer, CRS("+init=epsg:3
浏览 2提问于2019-07-13得票数 0
回答已采纳
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
回答已采纳
所以昨天我更新了我的android studio到最新的版本(4.2)，在构建我的项目后，我收到了这个警告 Please remove usages of `jcenter()` Maven repository from your build scripts and migrate your build to other Maven repositories.
This repository is deprecated and it will be shut down in the future.
See http://developer.android.com/r/tools/jcent
浏览 63提问于2021-05-10得票数 0
回答已采纳
我需要加载一个光栅地图和一个.shp文件。
两者都以下列方式加载到全局环境中
library(raster)
library(maptools)
library(maps)
topo <- raster("topo Europe.tif")
roads <- readShapeSpatial("roman_roads_v2008.shp")
现在我在全球环境中都有。我可以独立地画出来
plot(plot(topo, col=gray.colors(100), axes=F, box=F, legend=F))
plot(roads)
然而，我只能看到
浏览 0提问于2018-03-16得票数 0
回答已采纳
我在我的Android项目中添加了几个依赖项(我相信都是与mapbox相关的)，现在我在代码中得到了以下警告：
如何删除这方面的所有警告？我知道我可以添加@SuppressLint来删除这个函数的警告，但是我想从我的应用程序中删除所有这些警告。
浏览 1提问于2019-07-10得票数 8
回答已采纳
3回答
我想用sf包在都柏林机场周围画一个110海里的圆圈。(稍后，我将通过st_intersect与ADS的航班位置报告进行交叉。)
我为NM定义了一个新的单位如下：
library(units)
library(tidyverse)
library(sf)
NM <- make_unit("NM")
install_conversion_constant("NM", "km", 1.852)
然后定义都柏林机场坐标：
# DUB/EIDW location, see
# https://skyvector.com/airport/EIDW/D
浏览 3提问于2017-10-12得票数 9
回答已采纳
1回答
我试着用Rstudio上的传单制作一张地图。我尝试了许多例子和自我选择，它们都正常工作。
但我想使用我在工作中的形状文件。它与所有的软件(QGIS，R spplot等)一起工作。但是当我在传单中使用形状文件时，它就会在地图的顶部出现一条蓝线。
下面是我使用的代码：
dat<-readOGR(dsn="shape/Shape ER 2015"，layer=filename，编码=‘UTF-8’)传单(Dat)%>%addTiles()%>% addPolgons(data=dat，weight=2)
下面是代码的结果：
这是用spplot获得的
浏览 3提问于2015-10-26得票数 1
回答已采纳
1回答
我以tif文件的形式下载了新西兰的海拔数据，其中附有一个"world file“(tfw)，
80.000000000000000
0.000000000000000
0.000000000000000
-80.000000000000000
1062571.750000000000000
6282711.000000000000000
我从wikipedia的理解是，第一个和第三个数字表示比例尺(像素/公里)，最后两个数字位于UTM坐标中(如果我正确理解的话)。
现在，我想用一个不同的数据集对齐这些光栅数据，该数据集由GPS坐标跟踪的路径组成，这给了我一个lat/long格式。
浏览 1提问于2018-09-26得票数 2
我已经将我的google-services插件更新到了4.3.0
当我同步我的项目时，我得到以下警告。
WARNING: API 'variant.getMergeResources()' is obsolete and has been replaced with 'variant.getMergeResourcesProvider()'.
It will be removed at the end of 2019.
For more information, see https://d.android.com/r/tools/task-configur
浏览 0提问于2019-07-04得票数 3
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
即时通信 IM活动推荐
腾讯云开发者
