我需要R语言的帮助
从我的代码:
inter1= read.table("C:/inter.csv", header=TRUE)
inter1$xx<-inter1$long
inter1$yy<-inter1$lat
coordinates(inter1) = ~long + lat
#Plot the results:
plot(inter1)
Run Code Online (Sandbox Code Playgroud)
我有这个情节:http://i.stack.imgur.com/98aTf.png
我现在正在寻找情节上绘制多边形的每一组点,我不知道这个过程我必须更加努力才能到达那里,谢谢你的帮助
inter.csv:
long lat var1.pred
1 4.2 19 31.8216045615229
2 4.3 19 31.913824396486
3 4.4 19 32.0090783396173
4 4.5 19 32.1067681024233
5 4.6 19 32.2061094352961
6 4.7 19 32.3061148156713
7 4.8 19 32.4055837134796
8 4.9 19 32.503104196147
9 5 19 32.5970697606984
10 5.1 19 32.6857147918646
11 5.2 19 32.767170733855
12 5.3 19 32.8395428348418
13 5.4 19 32.9010042955024
14 5.5 19 …Run Code Online (Sandbox Code Playgroud) 我正在寻找一种方法来计算列表空间多边形中每个多边形的重心:
我以为用了一个循环,但他让我为第一个多边形,我不知道的方式,我是R的新手,有人可以帮助我代码:
for ( i in 1:length(polys1_T)) {
xx=mean(coordinates(polys1_T[[i]])[,1])
yy=mean(coordinates(polys1_T[[i]])[,2])
aa<-as.data.frame(cbind(xx,yy))
}
Run Code Online (Sandbox Code Playgroud)
编辑:
码:
inter1 <- read.table("c:/inter1.csv", header=TRUE)
# add a category (required for later rasterizing/polygonizing)
inter1 <- cbind(inter1,
cat
= rep(1L, nrow(inter1)), stringsAsFactors = FALSE)
# convert to spatial points
coordinates(inter1) <- ~long + lat
# gridify your set of points
gridded(inter1) <- TRUE
# convert to raster
r <- raster(inter1)
# convert raster to polygons
sp <- rasterToPolygons(r, dissolve = T)
plot(sp)
# addition transformation to distinguish well …Run Code Online (Sandbox Code Playgroud) 我找不到在R中的数字小数点后面只显示2个数字的方法,我不想要对数字进行舍入,我只想显示数字2.
例如:
1.29964我想显示:1.29