Search Results

I have recently started to use SQL and geopackages queries in R. I know how to open it, access the individual layers (geometry and tables), and select and process individual layers. …

Here is a solution converting raster patches into polygons, and subsequently into lines, using this post library(terra) set.seed(123) # Create an empty raster and set all values to NA initially raste …

I have a set of points placed within one patch. I would like to measure the nearest distance between the point and the closest cell with NA values (eg. outside of my patch). Patch is defined as a cell …

I have a list of rasters of the same location for multiple years. The change in pixel value over time represents the time series of pixel. To further analyses, I need to extract the values over time …

I found another great answer here ! https://stackoverflow.com/questions/19833784/how-to-extract-values-from-rasterstack-with-xy-coordinates?rq=1 at he same time, so here it is ! No need to put raste …

The simplest answer... I had wrong position of the [i,] indexing. Instead of indexing the shapefile polys[i,]@data$open_edge <- TRUE the correct indexing is in shapefile attribute table slot polys@da …

FALSE, sep = "") my.tab.spdf<- SpatialPointsDataFrame(cbind(my.tab$Longitude, my.tab$Latitude), my.tab) and use Voronoi tesselation published here: http://carsonfarmer.com/2009/09/voronoi-polygons-with-r/ …

Here I post my answer: I missed out how to convert SpatialPolygon to SpatialPolygonDataFrame and difference between them. SpatialPolygon doesn't contain attribute table, SpatialPolygonDataFrame does. …

This will implement moving window in raster function "modal" fun=modal to keep the majority values of neighbour values movingFun() from the {Raster}is mostly for vectors R code for focal() # matrix( …

However I couldn't find any approach to do use this moving window for confusion matrix evaluation in R, normally it is used to values interpolation. …

I have a set of rasters. I crop them all using FOR loop. Then, I want to save new-created rasters in new directory, not in my actual directory. For one raster it is easy: writeRaster(raster1, filena …

I have created a list of rasters - "norm" and run a several loops. Now I would like to write every created raster from "norm" (containing 1 layer) into one raster, having its specific name. # I can g …

I am trying to check if my shapefile polygons have open edge or not. To have an open edge, the polygon either has no neighbors, or is >5 m taller than its neighbours. So, I want to add this new field …

OK, I found this after correcting my googling: To access indivudual attributes of an object, one can use attributes(x)$attributeName or attr(x, attributeName). So, in my case, to access the 'mean' col …

However, points attribute table is not transferred automatically in R to Voronoi tesselation polygons during calculation (as is the case in ArcGIS). …