Вопросы с тегом «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") …

10 r  ggplot2  rgeos 

У меня есть временной ряд спутниковых изображений (5 полос), и я хочу классифицировать их по kmeans в R. Мой скрипт работает нормально (просматривайте мои изображения, конвертируйте изображения в data.frame, группируйте их и конвертируйте обратно в растр): for (n in files) { image <- stack(n) image <- clip(image,subset) ###classify raster image.df …

10 r  classification 

Я использую levelplotфункцию удивительного rasterVisпакета, чтобы создать карту, значения которой расходятся около нуля. Я строю растровый файл, используя красную палитру синего цвета, но я пытаюсь присвоить серым нулевые значения на карте. В частности, я хочу воспроизвести цвета этой фигуры: Обратите внимание, что там, где используется шкала от красного до синего, …

9 raster  r  legend  visualisation  color-ramp 

Закрыто. Этот вопрос не по теме . В настоящее время он не принимает ответы. Хотите улучшить этот вопрос? Обновите вопрос так это на тему для географической информационной системы Stack Exchange. Закрыто 9 месяцев назад . У меня есть набор данных с 1000 разных широт-долгот. Я хочу извлечь среднегодовую температуру и …

9 r  climate  extract 

Я конвертирую вектор в растр в R. Однако процесс был слишком долгим. Есть ли возможность поместить скрипт в многопоточную или графическую обработку, чтобы сделать это быстрее? Мой скрипт в растеризованный вектор. r.raster = raster() extent(r.raster) = extent(setor) #definindo o extent do raster res(r.raster) = 10 #definindo o tamanho do pixel …

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 data …

9 shapefile  r  ggmap 

Я пытаюсь создать правильную choropleth()карту в R, используя GISToolsпакет. Я использую choro.legend()для отображения моей легенды. Однако я не могу создать ни правильное размещение легенды, И / ИЛИ Я не понимаю роль fmtпараметра в моей choro.legend()функции. Мне кажется, это fmtможет уменьшить пространство между цветами моей легенды и описанием. Я нашел страницу …

9 r  legend  choropleth 

Я хочу создать точечный векторный набор данных из 10000 точек (или больше) внутри двоичного растра, где точки должны быть ограничены областями, где значение растра равно 1. Я попробовал следующие шаги. Многоугольный растр QGIS: Вектор -> Инструменты исследования -> Случайные точки Это прекрасно работает до 2000 очков, но все, что выше, …

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

Я создал data.frame в R. В результате я хочу получить 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()функция to this (ниже полностью воспроизводимого примера): # 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 …

9 raster  gdal  r  coordinate-system  gdalwarp 

Я хочу загрузить файлы .adf в R. Данные с этой страницы: http://www.fao.org/geonetwork/srv/en/metadata.show?id=14057. Я попробовал следующий код, который я нашел после некоторых исследований в Интернете. Проблема в том, что в классе RasterLayer я получаю отрицательные значения, которых там быть не должно. Я не знаю, почему это происходит, так что, надеюсь, кто-нибудь …

9 raster  data  r  esri-grid-format 

Я хотел бы выполнить пространственное наложение, чтобы идентифицировать полигон, в который попадает набор точек. Однако сначала я хочу буферизовать и растворить многоугольники, чтобы точки, падающие где-либо внутри объединенных многоугольников (но не внутри каких-либо отверстий), были помечены аналогично с помощью процедуры наложения. К сожалению, процесс буферизации и / или растворения, который …

9 r  sp 

Я хотел бы импортировать все множественные растры мозаики (каждая мозаика состоит из 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) …

9 raster  r  remote-sensing  landsat 

Я создал два SpatialLinesобъекты в R: . Эти объекты были созданы таким образом: library(sp) xy <- cbind(x,y) xy.sp = sp::SpatialPoints(xy) spl1 <- sp::SpatialLines(list(Lines(Line(xy.sp), ID="a"))) Теперь я хочу как-то сделать вывод, что это одна и та же линия, повернутая и перевернутая, и что разница между ними равна 0 (т.е. форма равна). …

9 r  geometry  similarity 

Я пытаюсь найти несколько центроидов шейп-файла для датских муниципалитетов, а затем найти время в пути между ними. Я использую R«s readShapePolyфункцию из в maptoolsсочетании с gCentroidфункцией от rgeos, и все работает. Тем не менее, я получаю пространственные точки, такие как SpatialPoints: x y 1 571860.7 6225016 Coordinate Reference System (CRS) …

9 shapefile  r  coordinate-system