«r» 태그된 질문

ggplot2를 사용하여 공간 데이터를 플로팅하는 데 문제가 있습니다. spplot을 사용하여 플롯하면 맵이 잘 보이므로 요새화 단계에서 찢어짐이 발생한다고 가정합니다. 코드는 다음과 같습니다. #install the packages library(rgdal) library(mapproj) library(raster) library(rgeos) library(ggplot2) library(plyr) if (!require(gpclib)) install.packages("gpclib", type="source") gpclibPermit() setwd("C:/Users/My Documents") #read in laa to regional mapping #must aggregate to higher level regions …

10 r  ggplot2  rgeos 

나는 일련의 위성 이미지 (5 밴드)를 가지고 있으며 R에서 kmeans로 분류하고 싶습니다. 내 스크립트가 정상적으로 작동합니다 (이미지를 반복하고 이미지를 data.frame으로 변환하고 클러스터링 한 다음 다시 래스터) : for (n in files) { image <- stack(n) image <- clip(image,subset) ###classify raster image.df <- as.data.frame(image) cluster.image <- kmeans(na.omit(image.df), 10, iter.max = 10, …

10 r  classification 

나는 levelplot멋진 rasterVis패키지 의 기능을 사용하여 값이 0으로 갈라지는 맵을 만듭니다. 빨간색에서 파란색 팔레트를 사용하여 래스터 파일을 플로팅하고 있지만 맵에서 회색을 0으로 값을 할당하려고합니다. 특히이 그림의 색상을 재현하고 싶습니다. 빨간색에서 파란색으로 스케일이 사용되지만 0 값은 회색으로 표시됩니다. 현재 내지도는 다음과 같습니다. 그리고 이것은 그것을 재현하는 코드입니다 ( https://www.dropbox.com/s/cypfdu1eaz2fuok/r.annual.tif?dl=0에있는 파일 ) …

9 raster  r  legend  visualisation  color-ramp 

닫은. 이 질문은 주제에 맞지 않습니다 . 현재 답변을받지 않습니다. 이 질문을 개선하고 싶습니까? 지리 정보 시스템 스택 교환에 대한 주제가 되도록 질문을 업데이트하십시오 . 휴일 7 개월 전 . 위도 경도가 1000 개인 데이터 세트가 있습니다. 이 좌표 각각에 대한 평균 연간 온도와 연간 강수량을 추출하고 싶습니다. 이러한 데이터는 …

9 r  climate  extract 

R에서 벡터를 래스터로 변환하고 있지만 프로세스가 너무 길었습니다. 더 빠르게하기 위해 스크립트를 멀티 스레드 또는 GPU 처리에 넣을 수 있습니까? 래스터 화 된 벡터에 대한 스크립트입니다. r.raster = raster() extent(r.raster) = extent(setor) #definindo o extent do raster res(r.raster) = 10 #definindo o tamanho do pixel setor.r = rasterize(setor, r.raster, 'dens_imov') …

9 raster  r  rasterization  rgdal  parallel 

ggmap을 사용하여 일부 위치 지점이 포함 된 맵의 셰이프 파일에서 지방 자치 경계 (폴리곤)를 포함하고 싶습니다. 이 스크립트는 다각형 플로팅을 제외한 모든 작업을 수행합니다. library(rgdal) library(ggmap) # Get shapefile with Drammen municipality borders tmpzip<-tempfile() tmpdir<-tempfile() dir.create(tmpdir) download.file("http://www.kartverket.no/Documents/Kart/N50-N5000%20Kartdata/33_N5000_shape.zip",tmpzip) unzip(tmpzip, exdir=tmpdir) kommune <- readOGR(dsn=tmpdir, layer="NO_AdminOmrader_pol") kommune<-kommune[kommune$NAVN=="Drammen",] kommune<-spTransform(kommune, CRS("+init=epsg:4326")) # Get location point …

9 shapefile  r  ggmap 

패키지를 choropleth()사용하여 R에서 올바른 맵 을 만들려고 GISTools합니다. choro.legend()범례를 표시 하는 데 사용 합니다. 그러나 올바른 범례 배치를 만들 수 없으며 함수 fmt에서 매개 변수 의 역할을 이해하지 못합니다 choro.legend(). fmt내 전설의 색과 설명 사이의 공간을 줄일 수있는 것 같습니다 . 내가 발견 choro.legend 도움말 페이지 국가는 fmt이다 "choropleth 클래스 …

9 r  legend  choropleth 

이진 래스터 내에 10000 포인트 이상의 포인트 벡터 데이터 세트를 만들고 싶습니다. 여기서 래스터 값이 1 인 영역으로 포인트를 제한해야합니다. 나는 다음 단계를 시도했다. 다각형 래스터 QGIS : 벡터-> 연구 도구-> 임의의 점 이것은 최대 2000 포인트까지 잘 작동하지만 위의 모든 것은 QGIS 충돌을 유발합니다. 이진 래스터 (또는 다각형 버전)에 …

9 qgis  arcgis-desktop  python  r  point-in-polygon 

R에 data.frame을 만들었습니다. 결과로 SpatialLinesDataFrame을 가져 오려면 data.frame을 선으로, 선을 SpatialLines로, SpatialLines을 SpatialLinesDataFrame으로 변환해야한다는 것을 이해하고 싶습니다. 나는 이걸했다 filedata.frame=data.frame(matrix(file),ncol=14, byrow=T)) file=Lines(filedata.frame) 원하는 data.frame을 얻지 만 선은 얻지 못합니다. 나는 기본이지만 확실하지 않은 이유는 무엇입니까?

9 r  line 

래스터를 투영하려고합니다. R에는 이것에 대한 projectRaster()기능이 있습니다 (완전히 재현 가능한 예 아래). # example Raster require(raster) r <- raster(xmn=-110, xmx=-90, ymn=40, ymx=60, ncols=40, nrows=40) r <- setValues(r, 1:ncell(r)) projection(r) # project to newproj <- "+init=epsg:4714" # using raster package to reproject pr1 <- projectRaster(r, crs = CRS(newproj), method = 'bilinear') …

9 raster  gdal  r  coordinate-system  gdalwarp 

.adf 파일을 R로로드하고 싶습니다. 데이터는 다음 페이지에 있습니다. http://www.fao.org/geonetwork/srv/en/metadata.show?id=14057 인터넷에서 조사 한 결과 다음 코드를 시도했습니다. 문제는 RasterLayer 클래스에서 거기에 있으면 안되는 부정적인 가치를 얻는다는 것입니다. 왜 이런 일이 발생하는지 모르겠으므로 누군가가 나를 도울 수 있기를 바랍니다!? 암호: library(rgdal) library(RColorBrewer) dpath<- path... x <- new("GDALReadOnlyDataset", dpath) getDriver(x) getDriverLongName(getDriver(x)) xx<-asSGDF_GROD(x) r …

9 raster  data  r  esri-grid-format 

표시된 폴더에서 여러 개의 모자이크 래스터 (각 모자이크는 7 개의 레이어로 구성됨)를 내 R 로 가져오고 싶습니다. 그런 다음 단일 멀티 레이스 터 래스터로 액세스하십시오. 그래서 나는했다 : # read all mosaics named "mos....img" in R raster_data<-list.files(path=getwd(), pattern="mos.*.img$") # read files as rasters s <- stack(raster_data) # check my imported …

9 raster  r  remote-sensing  landsat 

점 집합이 속하는 다각형을 식별하기 위해 공간 오버레이를 수행하고 싶습니다. 그러나 먼저 다각형을 버퍼링하고 분해하여 병합 된 다각형 내에 있지만 구멍 내에 있지 않은 지점에 오버레이 절차에 의해 유사한 레이블이 지정되도록합니다. 불행히도, 내가 사용하는 버퍼링 및 / 또는 용해 프로세스는 SpatialPolygons객체를 다중 다각형으로 줄이는 것 입니다. 를 사용 gBuffer하면의 경우 …

9 r  sp 

SpatialLinesR에서 두 개의 객체를 만들었습니다 . 이러한 객체는 다음과 같이 생성되었습니다. library(sp) xy <- cbind(x,y) xy.sp = sp::SpatialPoints(xy) spl1 <- sp::SpatialLines(list(Lines(Line(xy.sp), ID="a"))) 이제 어떻게 든 이것이 회전하고 뒤집힌 동일한 선이며 그 차이가 0과 같다는 결론을 내립니다. 이를 위해 maptools패키지 를 사용 하고 라인 # 1을 회전 시킬 수 있습니다 . …

9 r  geometry  similarity 

덴마크 지방 자치 단체의 셰이프 파일 중심을 찾으려고 노력한 후 그 사이의 운전 시간을 찾습니다. 내가 사용 R의 readShapePoly에서 기능 maptools와 결합 gCentroid에서 기능 rgeos및 모든 작품. 그러나 다음과 같은 공간 점을 얻습니다. SpatialPoints: x y 1 571860.7 6225016 Coordinate Reference System (CRS) arguments: NA 여행 시간을 잡기 위해 Google에서 …

9 shapefile  r  coordinate-system