R에서 지리 참조 정보가 포함 된 래스터 생성

이전 질문을 참조하여 생성 된 래스터 이미지를 R을 사용하여지도에 레이어로 표시 할 수있는 형식으로 "지오 코딩"/ 변환하는 방법은 무엇입니까? 요컨대, 상용 SpatialKey ( 스크린 샷 ) 소프트웨어 의 결과를 재현 하려고 합니다.

GeoServer 인스턴스를 사용하여 맵에 대한 데이터를 호스팅하고 있지만 필요한 투영 및 좌표 정보를 이미지에 통합하여 실제 지리적 위치에 올바르게 정렬되도록하려면 어떻게해야합니까? 밀도 열지도는 공간 데이터를 사용하여 생성되었습니다.

raster r georeferencing

— 미모
소스

이것은 프로그래밍 문제 일 뿐이므로 실제로는 주제가 아닙니다. spatstat 패키지 정보와 짧은 코스 pdf를 살펴 보았는지 확인하십시오. 나는 이것이 전에 올 것이라고 확신하고 빠른 Google 검색은 이 관련 R-help answer을 나타 냅니다. 더 직접적인 솔루션과 검색 기능이 더 있으면 놀라지 않을 것입니다.

— Andy W

생성 된 래스터 이미지를 어떻게 만들었습니까? 다른 응용 프로그램에서 나왔거나 R 내에서 만들어 졌습니까?

— Simbamangu

image () 함수를 사용하여 R 내에서 생성됩니다.

— Mimo

x의 수평 및 수직 해상도가 다릅니다. 이러한 데이터는 Arc-ascii 형식으로 저장할 수 없습니다. My x의 치수는 4684, 1950이며 해상도에 대해 불평합니다. 나는 그것을 동일하게 만들었다. Buggy

답변:

래스터 패키지를 사용하면 임의의 래스터 (크기, 투영)를 만들거나 행렬과 같은 기존 객체를 사용할 수 있습니다. 그런 다음 투영 및 범위를 지정할 수 있습니다. 귀하의 image()기능을 직접 사용할 수 행렬 인수를 :

library(raster)    
## Create a matrix with random data & use image()
xy <- matrix(rnorm(400),20,20)
image(xy)

# Turn the matrix into a raster
rast <- raster(xy)
# Give it lat/lon coords for 36-37°E, 3-2°S
extent(rast) <- c(36,37,-3,-2)
# ... and assign a projection
projection(rast) <- CRS("+proj=longlat +datum=WGS84")
plot(rast)

그런 다음 writeRaster()다양한 형식으로 사용할 수 있습니다 ( 예 :

writeRaster(rast, "~/myraster.asc", format = "ascii")

— 심바 만구
소스

여러 가지 방법이 있지만 Windows에서 오픈 소스 도구를 사용하려는 찔러 세 가지 옵션을 제안 할 수 있습니다.

GDAL (오픈 소스 지오 래스터 조작 라이브러리) 에서 gdal_translate 도구를 사용하십시오 . -gcp pixel line easting northing elevation래스터 ( http://www.gdal.org/gdal_datamodel.html ) 에서 픽셀과 선이 x / y 이고 지리적 좌표를 동쪽 / 북쪽으로하는 매개 변수가 필요할 수 있습니다. 또한 투영을 지정해야합니다.-a_srs srs_def
GDAL에 대한 r 바인딩 을 사용하고 동일한 작업을 수행하십시오.
더 쉬운 "GUI"방법은이 튜토리얼 을 따르고 GDAL과 함께 QGIS를 사용하는 것입니다. 이것의 단점은 QGIS와 GDAL을 설정하지 않았다면 설정해야한다는 것입니다.

어떤 경로를 선택하든 내려 가기 전에 필요한 각 소프트웨어 비트의 종속성을주의 깊게 살펴 보겠습니다.

— Stev_k
소스