Statistical computing language and software environment.

learn more… | top users | synonyms

3
votes
1answer
153 views

problem with polygon smoothing in R

With R I am looking to smooth my Polygon: inter2.csv code: inter1 <- read.table("c:/inter2.csv", header=TRUE) # add a category (required for later rasterizing/polygonizing) inter1 <- ...
3
votes
1answer
99 views

make a raster from dataset (csv)

I was trying to make a raster from an dataset (https://www.dropbox.com/s/jsm1x06lnzfbqty/data.csv?dl=0) code: inter1= read.table("c:/data.csv", header=TRUE) #convert to spatial points ...
3
votes
1answer
269 views

Polygon fill colour specified by column R spplot

I am trying to plot a geological map using the spplot function in R. For every unit I already assigned a specific color in the col-column of the SpatialPolygonsDataFrame. Now I would like to plot the ...
2
votes
1answer
147 views

How to change position of legend in levelplot of rasterVis in R? [closed]

By default legend is plot right side and I want to change it to bottom of plot. levelplot(CHL1, xlim=c(30,130), ylim=c(-20,30), xlab="", ylab="", margin=FALSE, col.regions=pal, ...
3
votes
2answers
237 views

Reprojection of shapefile in R and final extraction of values

I need to extract values from some raster files (stack) using a shape file. to bring both of them in the same CRS, I first transformed the shapefile using the EPSG codes. Opening it in an GIS software ...
2
votes
1answer
215 views

Count overlapping polygons in single shape file with R

I have a shapefile of species ranges for ~ 346 species. Some species have multiple polygon ranges so there are 1,078 polygons and many polygons are multi-part. It looks something like this: I ...
2
votes
1answer
82 views

Extracting data, not colors, from geotiffs in R

I'm an experienced R user new to GIS, but finding my way. I've recently run into an issue that I haven't been able to find any satisfactory documentation on (but could just be searching poorly). I'm ...
0
votes
1answer
75 views

How to create an interrupted Mollweide map in R?

I am wondering if it is possible to create an interrupted Mollweide map in R. This question is similar to Discontinuous or Interrupted map projections? which does not really have a satisfactory ...
3
votes
2answers
138 views

Extract values from raster based on points and pixels adjacent to the point

I have an issue. I have a set of points (n=~2000) and I would like to extract the pixel values from multiple rasters (a single-raster solution is also ok) for location of that point as well as that of ...
1
vote
0answers
61 views

R package spatstat: How to use point process model covariate as factor when pixel image values are numeric

I am trying to model a point process with an image covariate using the ppm() function in the spatstat package in R. I convert my raster to an im object for use with spatstat, and I run into a problem ...
4
votes
2answers
1k views

Converting osm file to shapefile (or data frame) in R

I was looking for a road shapefile for Australia today. I ended up going to Geofabrik. Open Street Map has data in .osm format. But they do not have one in .shp for Australia. If possible, I want to ...
3
votes
1answer
76 views

How to plot a full hemispheric orthographic projection in R

I want to make a plot like that using R: All I want is the country borders as seen from outer space, no fancy colors needed. It seems I need to use the spTransform function, but it gives me a "non ...
3
votes
2answers
241 views

Ripley's K function in QGIS Chugiak 2.4

I'm a noob user and I need to apply the Ripley's K-function on a point pattern in QGIS 2.4 Chugiak, but I don't know how to do. I think I need to use some plugin or extension (maybe GRASS, R or ...
2
votes
2answers
258 views

Convert shapefile with lists to data frame in R

I would like to ask some help for converting a shapefile to data frame. I downloaded data related to roads from Natural Earth. I imported the shapefile using shapefiles package like the following. ...
1
vote
0answers
118 views

Using r, readShapeSpatial fails with simplified mapshaper.org shape

So I'm mapping out the postal regions of my country and it's an enormous, 5.4million row dataframe. I don't need every single little twist and turn of the boundries, so I was super excited to ...
1
vote
1answer
88 views

How to convert raster to im object for point process model covariate in R

I am trying to convert a raster object to an .im object for use with a point process model in the spatstat package in R. I begin by creating the raster from a tiff file using the raster() package. No ...
1
vote
1answer
173 views

How to find highest kernel density point

I've been using spatstat to analyse point patterns in R. Here's an example to put this in context: library(spatstat) library(maptools) # to convert sp class to ppp Some point data: set.seed(1985) ...
0
votes
1answer
119 views

Subset netCDF based on last valid variable by level

I have a netCDF file (in this LINK) which contains monthly average values for 4 parameters at 11 different depth levels (0-200m). I would like to extract the monthly value for temperature "votemper" ...
1
vote
0answers
28 views

How to convert shapefiles to rasters? [duplicate]

I have my environmental variables such as elevation, soil, vegetation and others as shape-files. I want to use a raster command to upload them as tiff files for running a boosted regression tree in ...
0
votes
0answers
64 views

R PostGIS convert SpatialPolygon to PostGIS table

I've been having some success in scraping some webdata and converting it from HTML to spatial data as originally asked here: Parse PHP/HTML to shapefile or similar I've been using R to both fetch the ...
1
vote
2answers
255 views

Transforming and Projecting Shapefiles in R & ArcGIS for Desktop?

Perhaps this is an issue with one of my original data sources, but I've searched extensively to address this problem both through ArcGIS and R. I've got 2 separate sources of a county layer SourceA ...
1
vote
1answer
92 views

Clump and eliminate classes in a raster in ArcMap or R

after interpolation of sediment classes, I obtained a raster which looks like the image below (sample of the river). I would like to clump and eliminate areas within the raster that are below a ...
3
votes
0answers
59 views

Resampling air temperature gridded data to sub-grid resolution

I'd like to obtain a 30-m resolution gridded dataset of air temperature from the ~1km gridded data of WorldClim, of Hijmans et al., 2005) for a mountainous region. The easiest procedure (I thought ...
1
vote
1answer
137 views

Problem in creating neighbors, area and perimeter in spatial econometrics using R

I’m working on a spatial cross section model with R-programming-language. I’m running these codes which refers to anselin data (spatial_columbus) in spdep package. I found them from these links: ...
0
votes
0answers
49 views

Doing joint hypothesis test in spatial regression using R?

I'm doing estimation of spatial lag/error model using R package "spreg"/"spdep" but I cannot find any method to do anhypothesis test after regression. For example, I want to test whether two ...
1
vote
1answer
156 views

How to open and plot Globcolour chlorophyll image in R?

Level 3 globcolour chlorophyll images (download from hermes) are in netCDF format. Those files can be imported via ncdf package to R 3.1. ch.nc = ...
5
votes
1answer
288 views

Fixing orphaned holes in R

I'm trying to perform a union on a common field after merging two adjacent shapefiles. The shapefiles end up with at least one thin sliver of space between them. When I attempt a union I get the ...
5
votes
4answers
235 views

Create a rectangular buffer from point shapefile using QGIS/GRASS/SAGA GIS/R

I have a points shapefile and I want to create a rectangular buffer (22 x 16 km) around each point. My goal is to create a 16x22 km grid where each of the points is the centroid of the grid. Seed the ...
1
vote
2answers
337 views

Avoiding R for windows gui front-end has stopped working error?

I am trying to load huge number of points (shape file) in R. But when I load it, the program displays following error message and just stops. "R for Windows GUI front-end has stopped working" My OS ...
2
votes
1answer
89 views

How to create header files for multiple rasters in R

I have nearly 1000 generic binary images, i.e., headerless, band sequential and binary (16-bit signed integer). I have been succesfully able to load these with package 'raster' only when I create an ...
1
vote
1answer
100 views

Adding a grid and selecting random cells from spatial polygon in r

I have a fairly similar question to this: Similar Question I am working with a map of Denmark or more a subset of a Denmark Map. This file comes from the Global administrative areas site GADM site ...
2
votes
0answers
89 views

Plotting Bare Soil Line on landsat data in RStudio?

I am using RStudio with landsat package for computing and plotting the Bare Soil Line (BSL) on landsat tif image running on my HP Pavilion i7 laptop. The codes for running the BSL() function is as ...
4
votes
1answer
146 views

Why covariance values are different in R and ArcGIS?

I've calculated covariance between 7 rasters using R and ArcGIS. It gives me different results. Each times covariance values from R are approximately 1.874 times higher than ArcGIS values. R - ...
0
votes
1answer
107 views

Extract single datapoint from netCDF file in R using ndcf

I'm trying to extract a single variable from a 3-dimensional netCDF file. [1] "file blah.nc has 3 dimensions:" [1] "longitude Size: 200" [1] "latitude Size: 100" [1] "time Size: 2000" [1] ...
3
votes
2answers
134 views

How to iterate through features in a shapefile using R?

I am interested in iterating through a point shapefile to perform tasks on a point by point basis using R. In Python, this can be accomplished using the following workflow: import arcpy shp = ...
3
votes
2answers
90 views

Counting number of points in particular polygon for list of individuals using R?

I have distribution points for a set of individuals and a file of polygons-with matching IDs. What I would like to do is to go through and for each individual find the total number of points and the ...
4
votes
2answers
99 views

Map location of unknown point with distances to two known points

I'm trying to develop a function in R that will locate a point (A) on a map given distances between that point A and two other points (B and C), for which exact coordinates are known. The problem I'm ...
0
votes
1answer
47 views

How to batch export using r.out.gdal command in GRASS

I created a list of files (NPPFILES) using g.mlist command in GRASS 7. Now I would like to export these files as .tiff files. I tried this command in GRASS prompt for files in $NPPFILES; do ...
1
vote
0answers
65 views

Calculating surface and lines in 3D projection in R

For my internship research I am trying to compare the surfacearea and length of geometric shapes and lines in one coordinate system to another. I have a SpatialPolygonDataFrame in the European equal ...
5
votes
1answer
233 views

Convert ESRI Shapefile to Raster using R - biggest polygon area to raster cell value?

I want to create a raster from an ESRI shapefile. Multiple polygons may exist in one raster cell. If so, I want to assign the attribute variable which occupies the biggest area of the raster cell. ...
1
vote
0answers
67 views

Rounding off coordinates in an SpatialPolygonDataFrame in R

EDIT AFTER COMMENTS After a comment of Whuber (Thank you!) I realized it will probably be most easy to incorporate a round(...) statement within my code to be able to work with less decimals. ...
0
votes
0answers
43 views

How to reference a number variable in r.mapcalc within QGIS graphical modeler

I need to run the following formula in r.mapcalc: (xor(abs(A-B),abs(A-C))>120)*A I would like to however use a user defined number input variable instead of the 120. I already have this variable ...
2
votes
1answer
72 views

Is this a known issue in gaps between masked raster and spatial polygon in R with mask command?

This is my codes raster file: RAD polygon(shape): ws_sfc RAD.sub <- crop(RAD, extent(ws_sfc)) plot(RAD.sub) RAD.sub <- mask(RAD.sub, ws_sfc) plot(RAD.sub) plot(ws_sfc, add = TRUE) ...
0
votes
0answers
21 views

Kernel Density Analysis using R and ArcGis [duplicate]

I am doing some Kernel Density Estimations, and like everyone elese I face difficulties with setting the bandwidth to create georeferenced rasters. ArcGis (10.1) does it just fine, but it has no ...
2
votes
1answer
147 views

Plotting species distribution on map in R

I'd like to plot the distribution of a species on a map of Switzerland, following this tutorial. library(maps) library(mapdata) map('worldHires', 'Switzerland', xlim=c(5, 12), ylim=c(45,48)) From ...
1
vote
0answers
132 views

Unable to access R from QGIS 2.4

I have been unable to run R from within QGIS on my computer. I have R version 3.03 installed and QGIS 2.4 on Windows 7. R is located at C:\Program Files\R\R-3.0.3\bin\and this is the path e specified ...
3
votes
0answers
72 views

How to calculate 90% confidence interval in a raster stack

I have a stack of grid files (23) in ASC format from a MODIS sensor and I need to calculate a 90% confidence interval, mean, min, etc, for each one in R, how can I do that? specially the interval. ...
1
vote
0answers
40 views

QGIS using different R lib

I am working with a colleague who is learning QGIS and has a question about the R plugin. His library path (.libPaths()) has two folders: "C:/Users/RDCERWAW/Documents/R/win-library/3.0" "C:/Program ...
2
votes
0answers
89 views

Using R to subtract the polygon area of one shp file from another

I'm looking to produce a map from two separate maps in R. My Maps files are as follows: ftp://ftp.cpc.ncep.noaa.gov/GIS/us_tempprcpfcst/610prcp_20140822.zip #Map1 ...
1
vote
2answers
180 views

How to produce “SpatialPixelsDataFrame” from stacked RasterLayer?

How to produce "SpatialPixelsDataFrame" from staccked RasterLayer ? I have a stacked rasterlayer and I want to convert them to SpatialPixelsDataFrame. RasterLayer information is > RAD.all[[1]] ...