rgdal

    3热度

    2回答

    我想在R中的光栅文件上绘制shapefile,但我无法使它们完美重叠:光栅似乎逆时针旋转几度。这是投影问题吗? 请考虑以下MWE library(raster) library(rgdal) # Download from http://biogeo.ucdavis.edu/data/gadm2/shp/ITA_adm.zip shape_file = "ITA_adm1.shp" #

    0热度

    1回答

    我有一个包含900K个随机采样坐标点的数据集,其中每个坐标点的值都与每个点相关联。我想在地图上制作一个网格,并为每个单元格指定位于单元格内的所有点的平均值。最后,我想将网格绘制为热图。 电网的解决方案可能是100平方公里。小区达10000平方公里。

    1热度

    2回答

    我想从shapefile中获取纬度和经度。直到现在,我只知道如何阅读shapefile。 library(rgdal) centroids.mp <- readOGR(".","35DSE250GC_SIR") 但我怎么能从centroids.mp中提取经纬度?

    0热度

    2回答

    我从1983年开始在国家平面坐标系统中具有以下数据集。所有坐标位于长岛区(3104)。 dput(example) structure(c(1008031L, 1000852L, 1001869L, 1005306L, 986887L, 998031L, 1018703L, 1014319L, 1016186L, 1006977L, 1006891L, 1000883L, 1001403L,

    4热度

    1回答

    我试图在Rstudio(0.97.332)中使用Ubuntu 12.10(quantal)和R 2.15.1(烤棉花糖)安装rgdal库。 我在Ubuntu的安装以下关于@中: $ apt-cache search gdal dans-gdal-scripts - GDAL contributed tools by Geographic Information Network of Alaska

    7热度

    1回答

    我试图做同样的事情在这个问题Cartogram + choropleth map in R问,但从一个SpatialPolygonsDataFrame开始,并希望结束了相同类型的对象。 我可以将对象保存为shape文件,使用scapetoad,重新打开它并将其转换回来,但我宁愿将它全部放在R中,以便过程完全可重现,并且可以自动编写几十个变体。 我已经在Github上分叉了Rcartogram代码,

    2热度

    1回答

    我成功导入了一个带有readOGR的多边形shapefile文件,但未导入.prj文件中的信息。 测试 类:SpatialPolygonsDataFrame n功能:19407 程度:35551.4,1585917,6318047,9408727(XMIN,XMAX,YMIN,YMAX) 坐标。 REF。 :NA 摘要(测试)类SpatialPolygonsDataFrame 的 物体坐标: 最小

    1热度

    1回答

    我想知道如何根据特定颜色的矩阵编写一个GTiff?这可能是颜色不代表数量变量,而是定性空间现象,如人口统计集群。如果我有例如一个1000x1000的颜色值矩阵(像'#FF0000'这样的字符串),我可以很容易地用grid.raster(mtrx, interpolate=FALSE)来绘制这些矩阵,但这很难与其他形状文件叠加。我也可以将它们转换为一个SpatialPixelsDataFrame,字

    11热度

    4回答

    我发现了一个相当简单的例子,说明如何做到这一点,但我不能让它为我工作。我很新的R library(rgdal) xy <- cbind(c(118, 119), c(10, 50)) project(xy, "+proj=utm +zone=51 ellps=WGS84") [,1] [,2] [1,] -48636.65 1109577 [2,] 213372.05 554

    0热度

    2回答

    我有一个来自UTM区域33北部的大纬度经度信息。 我尝试下面的命令来这​​个地理信息转换为UTM坐标(我的数据集对象最初称为S3km): library(rgdal) UTM33N<-"+proj=utm+zone=33+north" UTM33N<-paste(UTM33N,"+ellps=WGS84",sep="") UTM33N<-paste(UTM33N,"+datum=WGS84"