R-경계 상자를 만들고 다각형 클래스 및 플롯으로 변환

나는이 NE lat lng와 SW lat lng. 내 목표는 간단합니다 .RStudio를 사용하고 위의 두 위도 및 경도에서 경계 상자를 만들고 다각형을 생성하고 사각형 다각형에 임의의 점을 추가하고 싶습니다. data.frame으로 가져올 수있는 CSV 파일에 NELat / Lng 및 SW-Lat / Lng이 있습니다.

나는 시작하기 위해 약간의 지침이 필요합니다 .R 언어에 매우 익숙합니다. 내가 익숙한 것과는 다릅니다.

지금까지 나는 이것을 가지고있다 :

coords = cbind(78.46801, 19.53407)
coordsmax = cbind(78.83157, 19.74557 )
sp = SpatialPoints(coordsmax)
sp2 = SpatialPoints(coords)

r1 = rbind(coords, coordsmax[1, ])  # join
P1 = Polygon(r1)
Ps1 = Polygons(list(P1), ID = "a")
plot(Ps1)

방금 다른 예제를 보았지만 4 모서리의 다각형을 그릴 수는 없습니다.

코드가 약간 변경되었습니다.

먼저 포인트 생성을 중단했습니다. 을 사용하지 않고 다각형을 형성 할 수 있습니다 SpatialPoints. 많은 점이 관련된 경우 점에서 다각형을 만드는 것이 좋습니다.

둘째, 아래 매트릭스에 5 개의 좌표를 썼습니다. 각 좌표 부부는 경계 상자의 한 구석을 가리키고 다섯 번째는 첫 번째 점을 반복합니다. 즉, 매트릭스는 [ (x_min, y_min) , (x_max, y_min), (x_max, y_max), (x_max, y_min), (x_min, y_min) ]을 포함한다

마지막으로, 나는 지리적 맥락에서 플롯 가능한 객체를 형성하는 데 사용 SpatialPolygons했습니다 espg:4326.

library(sp)

coords = matrix(c(78.46801, 19.53407,
               78.46801, 19.74557,
               78.83157, 19.74557,
               78.83157, 19.53407,
               78.46801, 19.53407), 
             ncol = 2, byrow = TRUE)


P1 = Polygon(coords)
Ps1 = SpatialPolygons(list(Polygons(list(P1), ID = "a")), proj4string=CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs"))
plot(Ps1, axes = TRUE)

코드를 플롯하면 이런 일이 발생합니다.

그리고 이것은 여기에 제시된 코드 수정 후입니다 :

extent래스터 라이브러리에서 함수 를 중첩하여 asSpatialPolygons 객체를 만들 수 있습니다.

library(sp)
e <- as(raster::extent(78.46801, 78.83157, 19.53407, 19.74557), "SpatialPolygons")
proj4string(e) <- "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs"
plot(e)