Warning: file_get_contents(/data/phpspider/zhask/data//catemap/4/r/75.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181
R 如何从空间多边形获取坐标_R_Gis_Sp - Fatal编程技术网

R 如何从空间多边形获取坐标

R 如何从空间多边形获取坐标,r,gis,sp,R,Gis,Sp,我使用以下方法从drawpoly获取多边形: x<- ogrListLayers(file_map) shape=readOGR(file_map, layer=x) #will load the shapefile to your dataset. plot(shape) polygon<- drawPoly() 提前感谢嗯,我发现功能点击poly并使用: polygon <- clickpoly(add = TRUE, col = 2, lwd = 2) 为了得到多边形

我使用以下方法从drawpoly获取多边形:

x<- ogrListLayers(file_map)
shape=readOGR(file_map, layer=x) #will load the shapefile to your dataset.
plot(shape)
polygon<- drawPoly()

提前感谢

嗯,我发现功能
点击poly
并使用:

polygon <- clickpoly(add = TRUE, col = 2, lwd = 2)

为了得到多边形中的x,y点,我发现函数
单击poly
,并使用:

polygon <- clickpoly(add = TRUE, col = 2, lwd = 2)

要获取多边形中的x、y点

,如果添加MWE,则会有所帮助。这里有一个:

polygon = "POLYGON((186.547856450189 -46.4485209972227,182.523886136587 -87.3045772259105,170.606614118262 -126.59055812607,151.254014475596 -162.79672548354,125.209796759794 -194.531696591426,93.4748256519081 -220.575914307228,65.3272038227373 -235.621131625452,36.2687622716732 -244.435913500765,-9.05803921416017 -248.900213788796,-54.3848406999935 -244.435913500765,-97.9697586242894 -231.214573193132,-138.137849033921 -209.744281559306,-173.34547474097 -180.850129866096,-202.239626434179 -145.642504159047,-221.510885247941 -109.588514788964,-228.270643805406 -87.3045772259111,-232.294614119009 -46.4485209972234,-228.270643805406 -5.59246476853568,-216.353371787081 33.6935161316239,-197.000772144416 69.8996834890942,-170.956554428614 101.63465459698,-139.221583320728 127.678872312782,-103.015415963258 147.031471955448,-63.7294350630982 158.948743973773,-22.8733788344104 162.972714287376,17.9826773942774 158.948743973774,57.2686582944371 147.031471955449,93.4748256519074 127.678872312783,125.209796759793 101.634654596981,151.254014475596 69.8996834890953,170.606614118262 33.6935161316251,182.523886136587 -5.59246476853451,186.547856450189 -46.4485209972227))"
我不知道为什么这么难,但你来吧:

polygon <- lapply(polygon, FUN = function(x) readWKT(x))

poly1 <- sapply(polygon, FUN=function(x) x@polygons)
poly2 <- sapply(poly1, FUN=function(x) x@Polygons)
poly3 <- lapply(poly2, FUN=function(x) x@coords)
plot(poly3[[1]])

polygon如果添加一个MWE,它会有所帮助。这里有一个:

polygon = "POLYGON((186.547856450189 -46.4485209972227,182.523886136587 -87.3045772259105,170.606614118262 -126.59055812607,151.254014475596 -162.79672548354,125.209796759794 -194.531696591426,93.4748256519081 -220.575914307228,65.3272038227373 -235.621131625452,36.2687622716732 -244.435913500765,-9.05803921416017 -248.900213788796,-54.3848406999935 -244.435913500765,-97.9697586242894 -231.214573193132,-138.137849033921 -209.744281559306,-173.34547474097 -180.850129866096,-202.239626434179 -145.642504159047,-221.510885247941 -109.588514788964,-228.270643805406 -87.3045772259111,-232.294614119009 -46.4485209972234,-228.270643805406 -5.59246476853568,-216.353371787081 33.6935161316239,-197.000772144416 69.8996834890942,-170.956554428614 101.63465459698,-139.221583320728 127.678872312782,-103.015415963258 147.031471955448,-63.7294350630982 158.948743973773,-22.8733788344104 162.972714287376,17.9826773942774 158.948743973774,57.2686582944371 147.031471955449,93.4748256519074 127.678872312783,125.209796759793 101.634654596981,151.254014475596 69.8996834890953,170.606614118262 33.6935161316251,182.523886136587 -5.59246476853451,186.547856450189 -46.4485209972227))"
我不知道为什么这么难,但你来吧:

polygon <- lapply(polygon, FUN = function(x) readWKT(x))

poly1 <- sapply(polygon, FUN=function(x) x@polygons)
poly2 <- sapply(poly1, FUN=function(x) x@Polygons)
poly3 <- lapply(poly2, FUN=function(x) x@coords)
plot(poly3[[1]])

polygon没有太多的内容,但是
coordinates(shape)
能给你想要的吗?@GavinSimpson我想从drawPoly()函数中得到多边形的坐标,我可以使用polygon@polygon但无法从坐标(多边形)中获取从itsee这个问题来看:没有太多的东西要做,但是
坐标(形状)
能给你你想要的吗?@GavinSimpson我想从drawPoly()函数中得到多边形的坐标,我可以使用polygon@polygon但无法从坐标(多边形)中获取从itsee这个问题来看:没有太多的东西要做,但是
坐标(形状)
能给你你想要的吗?@GavinSimpson我想从drawPoly()函数中得到多边形的坐标,我可以使用polygon@polygon但无法从坐标(多边形)中获取,请参见以下问题: