Statistical computing language and software environment.

learn more… | top users | synonyms

1
vote
2answers
203 views

Leaflet output is grey

I would like to implement leaflet in shiny with R. I try to go with the simple examples I found online: library(leaflet) library(ggmap) somePlace <- ggmap::geocode("Vienna") # Information from ...
1
vote
1answer
39 views

Spatial correlograms for binary raster?

I have binary raster data (deforestation hotspot or not) that I would like to plot spatial correlograms for in R, to look at spatial autocorrelation over different spatial lags. I came across this ...
0
votes
1answer
68 views

Is there a R package that does spatial regularization?

By spatial regularization I mean the method described in Aitkinson and Curran (1995) "Defining an Optimal Size of Support for Remote Sensing Investigations", or Section 4.8 of Webster and Oliver (...
3
votes
2answers
44 views

Error to obtain coordinate reference system argument from a raster image in R

I got an error when i wanted to obtain the CRS from stacklayer. I don´t know what could it be. I need some advice or help to obtain it. This is my code. #LOAD RASTER LIBRARY >library(rgdal) >...
2
votes
1answer
25 views

Read multiple layers of KML file using R

I have a KML that 3970 layers. I am trying to read it into R using the following script: library(rgdal) lyr <- ogrListLayers("file.KML") mykml <- readOGR("file.KML","Layer1") This works but ...
1
vote
1answer
51 views

rgdal / ggmap : rotate map data and clip map

I'm trying to plot a map of Toronto with ward boundaries and geocoded data on events over years. But the layout is such that I end up with mostly wasted space in the map. To increase resolution for ...
3
votes
2answers
173 views

Creating OUTSIDE_ONLY buffer around polygon using R?

I would like to create buffers in polygons, surrounding my original polygon. However, I want to keep only the "outside polygon", not include also original polygon. In Buffer in ArcGIS I can just set ...
0
votes
1answer
36 views

Goodness of fit of variogram model using geoR

How can I find a measure of goodness of fit (let's say an R2) of a model fitted to a variogram by the function variofit using the geoR package? The only value I have is the "minimised weighted sum of ...
0
votes
1answer
52 views

Rescale raster in R

What's the best way to rescale a raster with range 0, 200 to range -1 to 1 in R? I'm currently using gdal_translate, but this requires that I write out a new file with the updated range. If possible ...
0
votes
1answer
26 views

Get value from list class (vector) in R [closed]

I have a simple question about corce(?) list class data into vector. I extract a tif file from polygon to achieve a list. Maybe the class of list is leaded by extract function. I check the list length ...
3
votes
1answer
41 views

Extract rasterlayer from extent class

I used extract() function to get remote sensing image cell value.The same code to process tif file produce error message. Code as follows: r1 <- raster("2000_2605_1.1_km_monthly_NDVI.tif") ...
0
votes
1answer
74 views

R equivalent of GRASS r.grow?

What is a good R equivalent of the GRASS function r.grow? The closest I've been able to find is the script suggested by whuber in response to Creating randomly shaped clumps of cells in a raster from ...
11
votes
1answer
686 views

Multispectral image segmentation for natural-resources applications using R

R has the capacity for image segmentation, although all of the examples I have come across utilize a single band for the segmentation (example). I am interested in combining R's capacity for random ...
4
votes
2answers
161 views

Changing raster calculation values - shorten this script?

i'm doing something quite straightforward regarding raster algebra but am struggling either to find the correct function or make a function work correctly; Basically I have 2 rasters, representing ...
3
votes
1answer
96 views

Put points on a projected world map

With the help from this post Lines on reprojected sp objects with Mollweide projection, I draw a world map using Mollweide projection. I also have some points with latitudes and longitudes, and I ...
2
votes
2answers
39 views

Reducing number of attributes (classes of marks) in marked point pattern using R?

I have a marked spatial point pattern loaded in R (from Shapefile). This dataset has 4 different "classes" of marks: class A, class B, class C and class D. Each class of marks is categorical. One can ...
4
votes
1answer
3k views

Interpreting Moran's I results

I have a sample of data taken from regions of a country! i want to test if there are any spatial autocorrelation on my data using Moran Indice test. the null hypothesis: is that there is no spatial ...
1
vote
0answers
17 views

R: how to export my Kriging result to raster? [duplicate]

I want to export results of my Kriging calculation to raster or some other format. My aim is to take this converted file and intersect it with another SpatialPointDataFrame (location of the point ...
4
votes
1answer
287 views

A faster polygon intersection tool than rgeos gIntersection in R?

I'm looking for a R function able to retrieve the polygon intersecting between two SpatialPolygonsDataFrame. I need to use R and the rgeos::gIntersection() is too slow in comparison with ArcGis...
1
vote
1answer
46 views

Splitting polygon into multiple fixed-width polygons along edge in R?

I have a polygon and would like to spit into grid with fixed width along with one edge. The result will look like the figure below. Are there algorithm to achieve my target? My language is R, but ...
2
votes
1answer
23 views

Dealing with spatio-temporal data: sp::over function in R

I have a dataset called carcass, which contains occurrences of livestock carcasses in Norway in the last 5 years (each raw corresponds to a single carcass, which has an ID number, along with longitude,...
1
vote
1answer
43 views

How to apply functions iteratively to SpatialLines in R?

Looking for a way to apply the 'rgeos::gIntersection' and 'rgeos::gDifference' function iteratively to SpatialLines, I created 3 SpatialLines as an example: #from the sp vignette: l1 = cbind(c(1,2,3,...
1
vote
0answers
30 views

How to download GPM rainfall data automatically?

How to download GPM rainfall data automatically? I tried to use DownloadGPM.R Function, but I think I had problem with Axel download accelerator. I guess I need to declare Axel as environment ...
5
votes
2answers
3k views

How to convert a spatialpolygon to a SpatialPolygonsDataFrame and add a column to the atrribute table

coast<-readShapeSpatial("coastline.shp") landc<-readShapeSpatial("landcover.shp) ro<-readShapeSpatial("roads.shp") bc<-gBuffer(ro,width=100) landc$ratings=1 landc$ratings[landc$LANDUSE_ID==...
0
votes
1answer
60 views

smoothing raster layer stack using R

I'm trying to do smoothing for .tif stack layer of 331 bands. I used this code in R : setwd("D:\\") raster_dir <- "D:\\Modis VI quality\\" VI <- brick(paste(raster_dir,"...
0
votes
0answers
25 views

What is the approach for comparison two spatial point sets on the global scale?

I have to point sets, both of them dispersed over my study area, which counts several districts. Now I want to find a way how to compare them between each other using spatial statistics on a global ...
1
vote
1answer
398 views

Warning: implicit list embedding of S4 objects is deprecated

From this code: na_cells <- as.data.frame(extract(r_temp, poly_contour, cellnumbers = TRUE)) Why do I get this warning? Warning in `[<-`(`*tmp*`, cnt, value = p@polygons[[i]]@Polygons[[j]]) :...
1
vote
2answers
40 views

How to overlay a polyshape to faceted maps in r with ggplot?

So i am quite new to the use of r and ggplot to produce high quality output maps for publications. I have managed to create a faceted map from my data following the approach proposed here http://...
1
vote
0answers
29 views

writeOGR R and add an icon symbol [closed]

I am wondering if there is a way to add an icon symbol to the R code I am using to export data in GPX format. In the navigation program I import the GPX file into all of the of coordinates are ...
0
votes
0answers
27 views

Smooth raster stack

How can I smooth a raster stack of 300 bands using R/Arcgis/Envi? I used this code: setwd("D:\\") raster_dir <- "D:\\Modis VI quality\\" x <- brick(paste(raster_dir, "laikipia_VI_quality....
2
votes
0answers
27 views

Path Matching for Journeys Using XY Coordinates [closed]

I am trying to find some suggestions on ways or algorithms to determine path similarity. An illustration of 3 paths using made up data is given below. The title should be finding similar paths instead ...
44
votes
11answers
8k views

What is the most useful spatial R trick?

R is becoming pretty strong tool for handling and analyzing spatial data. I learned some useful things through questions like these at SO and thought it might be useful to have something simmilar, but ...
1
vote
1answer
26 views

Swapped axes in R raster brick plot from netcdf

I am using the R raster package to import a tiny sample netcdf file, available here: https://www.dropbox.com/s/fht4g52p948vxva/persiann.ncml.nc?dl=0. The main data website is: https://www.ncdc.noaa....
6
votes
2answers
233 views

Transforming geostationary satellite image to lon/lat

The ICARE Data and Services Center distributes the SEV_AERUS-AEROSOL-D3 product in HDF5 format (http://www.icare.univ-lille1.fr/archive/?dir=GEO/MSG+0000/SEV_AERUS-AEROSOL-D3.v1.03/). The main ...
2
votes
0answers
36 views

Kernel Density with vector output and/or possible automation?

I'm currently working on a project that will have me calculate the kernel density of a number of shapefiles with varying parameters. My current workflow goes like this: Import shapefiles in QGIS 2....
0
votes
1answer
28 views

Raster Mosaic fault using R

I'm trying to mosaic two raster images using R. I used this script: setwd("D:\\Modis EVI\\") raster_dir <- "D:\\Modis EVI\\" r1 <- brick(paste(raster_dir,"evi8_stack.tif", sep="")) r2 <- ...
1
vote
1answer
37 views

Remove raster pixel in R

This is the 32 year trend NDVI (raster data) I would like to remove the green grids data by using R (the green pixel value is larger than 0.03) and this one is R code that i use but i cannot ...
0
votes
0answers
29 views

How to define the Northern Lambert Conformal Conic CRS?

The dataset I got my hands on uses the the Northern Lambert Conformal Conic CRS? I was not able to find a ESPG number using http://prj2epsg.org/search so I decided to define a custom CRS in R but I ...
4
votes
2answers
670 views

LSA-SAF data reprojection with GDAL and R

The Land Surface Analysis Satellite Applications Facility (LSA SAF) distributes several products derived from the EUMETSAT satellites data. I am working with R to read and process the data (the file ...
0
votes
0answers
49 views

How to convert .HDF to .Tiff and keep the number of the output files exactly as the input?

I have 332 .HDF files that I have processed in R to convert them to .Tiff. It has worked fine, however, the number of output files are not the same as the number of the input files! Input is 332, ...
1
vote
1answer
23 views

Error in rasterToPolygon trying to get several coverage layers

I have a very strange problem: I am trying to get a raster with coverage from 0 to 1 in each cell for types of vegetation, I started with a raster of vegetation type, but I needed a raster with 12 ...
1
vote
0answers
34 views

R script for find the minimum distance [closed]

I need an R code to select which polygon in a shapefile has got the minimum distance from a given point.
0
votes
1answer
45 views

how to create raster stack using R

I have 300 .tiff files, that i need to create raster stack of them. Can you please help me with R code for that? I can do in Arcgis, Envi, R. I have tried: setwd("F:\\MODIS\\Modis EVI\\HDF8 EVI") ...
2
votes
2answers
912 views

Is there an R equivalent of NEAR in ESRI ArcGIS?

I have 54 points in one data set A I have 300+ data points in a dataset B. I wish to find the distance for each point in dataset-set B, to points in dataset A. Nearest Neighbor packages appear to ...
9
votes
2answers
599 views

Help choosing most suitable bivariate choropleth mapping coloring scheme

I am trying to overlay a nominal variable (2 categories - rainfed vs irrigated) in a choropleth map over a numeric variable. In order to highlight the difference in pattern among the two categories,...
0
votes
0answers
21 views

How can I create a NetCDF from precipitation rasters data in R? [duplicate]

I have a lot of precipitation raster data (more 5000 files GeoTIF), created by kriging interpolation. All grid cells have the same size and the maps the same spacial extent I want merge all data in ...
2
votes
1answer
466 views

How to extract percentage of land cover from a raster according to multiple shapefiles?

I have a raster map that represents land cover of a region, and different shapefiles (total 120) in the same region. I want to know the percentage of every land cover for each shapefile. For example, ...
1
vote
1answer
25 views

How to work with Baccini et al.'s biomass map from ArcGIS?

Baccini et al. (2012) has produced one of the best biomass maps available for the world. One can find the map from one of ArcGIS repositories here. I downloaded the .pkinfo file for that map but I ...
1
vote
0answers
17 views

R Freezes when Plotting Spatial Data Frames

I have run into a problem where anytime I try and plot a Spatial DataFrame in R, the console freezes, regardless of the size of the file. I imagined this must be due to an issue with the path to one ...
0
votes
0answers
56 views

Call gdal_proximity() from R

I have a raster including eight land cover types (pixel values ranged from 1 to 8): > r class : RasterLayer dimensions : 9887, 30532, 301869884 (nrow, ncol, ncell) resolution : 10, 10 (...