Statistical computing language and software environment.

learn more… | top users | synonyms

2
votes
0answers
44 views

Failed to create a raster-brick from a netCDF file (.nc)

I am trying to build a raster-brick from a netCDF file (network Common Data Form) leave this link for more information () but to make a brick (which is part of the package raster) throws me the ...
1
vote
0answers
28 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 ...
5
votes
1answer
98 views

Krovak projection EPSG 5514: Correct data transfer between R and ArcGIS?

I have a problem which many ArcGIS, QGIS and R users from Czech Republic and Slovakia may encounter. In ArcGIS, I am working EPSG: 5514, Projected Coordinate System, defined as: S-...
0
votes
0answers
29 views

Why ArcGIS crashes when I load raster file created in R? (HFA format)

I was processing my raster operation in R, as I think it's easier then in ArcGIS: - read shp + raster1 - crop them - classify raster1 by location - create classified raster2 - export the ...
0
votes
1answer
32 views

rasterVis plot3D in R: obtain consistent x-y axes with extent() / projection() shown in RasterLayer object

I have a RasterLayer object z: class : RasterLayer dimensions : 100, 100, 10000 (nrow, ncol, ncell) resolution : 0.06534037, 0.03165736 (x, y) extent : -96.63908, -90.10505, 40....
0
votes
2answers
30 views

mapping with spplot

This question should be easy! I am using spplot for visualizing a map with suitable results. However, I would like to know how to set the limits of z, because the values are not being discriminated. ...
0
votes
0answers
23 views

Problems to export raster files in R software

I'm working with R software with raster and netcdf4 packages adn netcdf files to get information from SST. I got the information and now I'd like to convert in ascii or geotiff files to work in ArcMap ...
1
vote
0answers
27 views

Polygon to Raster: Cells with multiple values

I'm trying to use R to rasterize a number of shapefiles and I'm having trouble dealing with cells on the edge of shapes, in other words the if a cell is partially covered by an NA value and partially ...
1
vote
1answer
26 views

Proximity filter

The aim of what I'm currently messing around: Creating a distance buffer around points in a SpatialPointsDataFrame If the value a this point is > 0, always keep it If the value a this point is = 0 ...
0
votes
1answer
52 views

Creating SpatialPolygons in R

I am trying to do kriging with my data. But there is a problem here. I want to use krige command in gstat package. But I do not have "Spatial object with prediction/simulation locations". How I can ...
2
votes
1answer
15 views

Retain “region” elements for when using fortify() on a SpatialLinesDataFrame

I'm layering multiple plots using ggplot2 to create a map of Alaska that combines school location and the major road network. I've learned that despite being labeled "major roads" there is wide ...
0
votes
1answer
35 views

Finding neighbors for non adjacent polygons in R [closed]

I'm looking for a method or function in r-project that create a neighborhood matrix from a group of polygons that don't touch each other, like polygons of set of blocks separated by streets in a town. ...
0
votes
0answers
15 views

Why the krige.cv gstat function does not ask for a semivariogram model for every training set?

I'm performing kriging leave-one-out-Cross-validation with the krige.cv gstat function: Is it possible to specify the semivariogram model for every training set? or is there a reason why the model has ...
1
vote
1answer
41 views

csv points into raster using R

I have got the forest cover data of a region in africa (raster layer) and some data collected inside this region by a survey (.csv). In the .csv file I have got x,y coordinates and 3 other columns ...
0
votes
0answers
44 views

north arrow outside map

I am plotting 5 maps using spplot having a 2*3 layout. The space not used for maps should contain cartographic features such as a north arrow. Could someone guid me in the right direction? I know how ...
0
votes
0answers
26 views

R: how to correctly set number n for sample grid for IDW?

I'm trying to make IDW in R. In my tutorial, I found this approach which is easily understandable, expect of one point - why they used n = 6000 ??? how did they established the value to create a grid ...
0
votes
1answer
61 views

R: Voronoi tesselation from long/lat data?

I want to create a SpatialPointDataFrame and Voronoi Tesselation from my .csv data downloaded from EPA, available here: http://uloz.to/xboQc7cu/no2-calif-csv I convert my data in SPDF format in 2 ...
-3
votes
1answer
57 views

How to colour a SpatialLinesDataFrame with Rainbow

how to color a SpatialLinesDataFrame according to some given value (Rainbow) FRA=readShapePoly("FRA_adm0.shp") grid <- raster(extent(FRA)) # Choose its resolution. I will use 2 degrees of latitude ...
2
votes
0answers
35 views

Filtering Small Areas of Shapefiles Using R

I'm trying to use the method proposed by the @ialm in the Question "How to speed up the plotting of polygons in R?" There, he loads the maps using getData from the raster package, I tested and it ...
0
votes
0answers
17 views

Use SVM to make prediction based on two Raster input in R

Right Now I have two raster images, one is elevation another one is slope. And I already generate a SVM model in R. Right now I want to use this SVM model to make a classification prediction using the ...
2
votes
1answer
75 views

Create GPX track in R. writeOGR function

I have a gpx file where someone exported two different activities in two separate locations to a single track. I need to split them and avoid getting the erroneous line in the track layer that ...
2
votes
0answers
22 views

Creating a reference raster in R

I need to create global reference raster (100km grid, Eckert IV). Later, I will project other various rasters (with different resolutions and projections) into the resolution/projection of the ...
1
vote
0answers
29 views

Specifying the “idp” and the “nmax” in the krige.cv function

I'm performing IDW interpolation using the idw R function from the gstat package. For the purpose of determining what combination of "idp" and "nmax" produces the lowest RMSE in leave-one-out-cross ...
1
vote
1answer
36 views

How to deal with two dimensional lat and lon variables in a NetCDF file in R?

I'm trying to convert a NetCDF file to a dataframe in R using the RNetCDF package. The NetCDF file is very different from what I have experience with, in that the lat and lon values each have two ...
1
vote
1answer
36 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 ...
0
votes
2answers
54 views

R: creating an error matrix from two rasters

I have two land use rasters (6 classes each) and I'd like to create an error matrix (errors of omission and commission) from them both. The premise is that I have a very fine 'base' raster and have ...
0
votes
0answers
14 views

Raster Values converted to 3.4e+38 in R

I have read a .grd file into R as RasterLayer which should only contain 1, 0 and NA as cell values. Turns out all the NA cells contain 3.4e+38 as value now. How can I change that and why does this ...
2
votes
0answers
41 views

Error mapping data frame in R: is.data.frame(map) is not TRUE [closed]

I am trying to generate a map using the following data and code: fileURL <- c("https://drive.google.com/open?id=0B-jAX5hT2D3hNnVtLVhROENKRGs") require(data.table) require(tidyr) require(dplyr) ge....
1
vote
1answer
78 views

dynamic visualization of spatiotemporal data on a grid

I am handling spatiotemporal data and it is a quite new area for me. Basically I have data related to the position of cars at a given time in Montreal. I would like to dynamically visualize the ...
1
vote
0answers
38 views

Raster package limitations

I am very new to R and to raster package and I was hoping that someone could provide some guidance. I am trying to run a function that averages a series of rasters, 300 to be exact, they are in geotif ...
2
votes
1answer
34 views

What is a UK boundary shapefile available for download via persistent URL (for cURL) with parish-level detail?

I'm in the process of generating a series of choropleth maps for a journal article and am writing this up in R as an .Rmd. My aim is to set this up as reproducible research, that is in a form which ...
7
votes
1answer
104 views

How to decrease the distance between legend's items in R?

I'm trying to create correct choropleth() map in R using GISTools package. I use choro.legend() to display my legend. However I can't create either correct legend placement, AND/OR I don't understand ...
0
votes
0answers
36 views

R: why are missing choro.legend() items when displayed in zoom()? RStudio

I wish to create map of Oregon, using choropleth(). Using locator() trick I've revealed coordinates for my choro.legend() position. It seems correctly placed in my plotting window, however my legend ...
0
votes
2answers
46 views

How to check my EPSG and units in R when it is not defined in proj4string?

After reading my .shp in R using readOGR(), I tried to figure out what is my EPSG and what are my units in my projections? when I tape my.shp = readOGR(dsn=setwd(dsn = dsn, my....
0
votes
0answers
27 views

Extract raster values from RASTERSBRICK using XY coordinates falling within each polygon R

I have a rasterbrick s (10000 raster layers) that actually is a time series raster. Accompanying s is a SpatialPolygonsDataFrame with many sub-polygons which can be identified by names and/or IDs etc. ...
1
vote
1answer
82 views

R kriging in QGIS

I'm trying to work with the R Kriging script in QGIS (processing toolbox). I have installed R before directly to C:. I also have installed the packages raster, rgdal, sp and automap. When I'm starting ...
0
votes
1answer
29 views

Calculating the number of days above freezing in R

Im wondering if anyone has come across any data source that can be used to calculate the number of days where temperatures are greater than 0 degrees C, 10 C and 20 C at a global scale? I have a site ...
3
votes
2answers
131 views

Plot tilted map in R

I would like to plot a tilted map, similar to the image below. I would prefer to do this using ggplot2 but I'm glad if there is a solution using base R. Any ideas on how to do this? Here is a dataset ...
3
votes
1answer
35 views

Is it possible to start the GRASS GUI from R?

I am using the rgrass7 package in R in order to work with GRASS GIS. After having set the needed enviornment variables, I can now launch a GRASS session from R and and work on my GRASS data using ...
2
votes
0answers
24 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 ...
0
votes
1answer
29 views

Get feature type of a layer

I have to assign the feature type of a layer. Does anyone know an function, or way to access this information. I have noticed that the message of the function readOGR (package rgdal) contains an ...
0
votes
1answer
25 views

histogram of raster coordinates in R

I'm trying to generate an histogram of the y coordinates where raster values are 1. Mysteriously the following, different rasters, give me the same histogram. What could I be doing wrong? here are ...
1
vote
1answer
44 views

Using 'sp' method to subset SpatialPointsDataFrame with SpatialPolygon returns error in parLapply()

I have numerous very large data sets of ship tracking data (Automated Information System; AIS). These are originally in a file geodatabase and can be downloaded here (this is a smaller file for ...
0
votes
1answer
58 views

Use of Leaflet plugins in 'Leaflet for R' to connect with ArcGIS online

I was wondering if it is possible to use the plugins for the normal (non-R) version of Leaflet, and if, how I can do that. reason I need to do that is, that some of the data I need to show are hosted ...
1
vote
0answers
52 views

Clipping raster brick object in R

I have MODIS NDVI TIFF files for a year I want to first clip it with a irregular polygon (not by its extent only) of my study area and save it as a multi band raster in TIFF format The TIFF image are ...
0
votes
1answer
31 views

Add information from a raster nc-file to a polygon shapefile

DATA Download from this website the following ZIP file: NUTS_2006_60M_SH.zip and save all the NUTS_RG_60M_2006 files. In what follows I refer to this as "NUTS3 regions". Download from this website ...
1
vote
1answer
32 views

Deciphering coord. ref. in R – very basic question

I'm trying to make sense of the coord. ref. parameters for rasters in R. For example, for a given raster, crs() gives the following: +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 ...
0
votes
0answers
70 views

collocated cokriging with gstat in R

I am trying to perform cokriging with temperature as the target variable and elevation as a secondary variable. Elevation (the secondary variable) is collocated with temperature (the target variable) ...
0
votes
0answers
29 views

Using the extract function in raster package for R with polygons

I have a shape file that contains some polygons and an unprojected raster. I'd like to extract the cells within the polygons using R, so I use code like this: # Load my polygons MyPoly <- ...
3
votes
2answers
80 views

How to generate a buffer outside polygons having a given pixel value?

For example, I would like to build a 20-m buffer (represented in yellow in the image below) outside polygons having a value of 2 (represented in green) and to assign the value of 3 to buffer cells. ...