Запитання з тегом «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 …

10 r  ggplot2  rgeos 

У мене є часовий ряд супутникових зображень (5 діапазонів) і хочу класифікувати їх по kmeans у Р. Мій сценарій працює чудово (перебирайте зображення, конвертуйте зображення у 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 

Зачинено. Це питання поза темою . Наразі відповіді не приймаються. Хочете вдосконалити це питання? Оновіть питання, щоб воно було актуальним для обміну стеками географічних інформаційних систем. Закрито 9 місяців тому . У мене є набір даних з 1000 різних широт-довгот. Я хочу отримати середню річну температуру та річну кількість опадів …

9 r  climate  extract 

Я перетворюю вектор у растр у Р. Однак процес був занадто довгим. Чи є можливість поставити скрипт в обробку багатопотокових чи графічних процесорів, щоб зробити це швидше? Мій сценарій для растрованих вектор. r.raster = raster() extent(r.raster) = extent(setor) #definindo o extent do raster res(r.raster) = 10 #definindo o tamanho do …

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 …

9 shapefile  r  ggmap 

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

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()функція цього (нижче повністю відтворюваного прикладу): # 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 шарів) із зазначеної папки в мій Р. Потім отримайте доступ до них як одиночні багатошарові растри. Так я і зробив: # 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) arguments: NA …

9 shapefile  r  coordinate-system