Statistical computing language and software environment.

learn more… | top users | synonyms

1
vote
1answer
80 views

Connecting Grass with R

I need help to connect R with Grass. In RStudio when I send the line: loc <- initGRASS("C:/Program Files/R/R-3.1.0/bin/i386", home=tempdir(), override=TRUE) It says: Error in ...
0
votes
0answers
70 views

Unserialize socklist error when using rasterEngine in r

Basically, I am running the "rasterEngine" function in R. I want to use parallelization to change values in a raster image (I have around 400 images that I need to work on). However, the rasterEngine ...
1
vote
1answer
21 views

RasterStack extract() for duplicated dates

I have a a Large RasterStack of temperatures and a 10570 points in a Large SpatialPointsDataFrame and I'm using extract() to assign the temperatures to each point in R. ...
0
votes
0answers
31 views

How to get specific coordinate system to get arcGIS point shapefile in r?

I have a Point shapefile from arcGIS 10. I want to get it in r and also want to get it a specific predefined coordinate systen for my data. I can read the shapefile in R but summary in r Shows ...
3
votes
1answer
48 views

writeOGR column limit error: “Creating Name field failed”

I have a SpatialPointsDataFrame with 1057 elements and 4760 columns and want to export it to a shapefile. When I try to do this with writeOGR() I get an error: > ...
0
votes
0answers
16 views

Creating a loop to predict over small areas

Based on a SDM, I am I am trying to use the predict function (raster package) to predict the model to South Africa. It runs off memory and gives me an error when I predict for the entirely country, ...
0
votes
0answers
29 views

How do I peform a true GIS Clip of polygons layer using a polygon layer in R

I'd like to do a true GIS Clip in R of soils polygons using a series of single boundary polygons, but I cannot find an R function to properly do it. It should work just like the CLIP function in ...
0
votes
0answers
48 views

How to open a raw binary raster in QGIS or R

I have a raw binary raster without a file ending, and no header file. I would like to open this in either QGIS or R and then convert to geotiff, but opening it seems to be a challenge. Here is a ...
0
votes
1answer
111 views

Saving kriging information from automap package in the ESRI ASCII raster format

I would like to save the information generated from an ordinary kriging (processed with the function autoKrige from the automap R package), into a file extension readable in ArcGis for desktop 9.3.1 ...
2
votes
1answer
55 views

Averaging XY points locations if in close proximity automatically over a large dataset in R

I have a point dataset (XYZ values) containing a few million points covering a 5km x 5km extent representing maxima locations. However there a great deal of these points which exist in 0.5m of others. ...
0
votes
0answers
38 views

calculate distance between points crossing the antimeridian

How can I calculate the distance between two points (lon/lat) if any of the points exceeds the antimeridian (+-180 longitude)? Im using the distVincentyEllipsoid {geosphere} to calculate the shortest ...
6
votes
2answers
213 views

Calculating spatial lags per year in R

At the moment I am having some difficulties with calculating a spatial lag in R. I know how to calculate the lag in space-wide format but am unable to do it in long form, i.e. have repeated ...
2
votes
1answer
149 views

Splitting strings into new rows, retaining geometry

UPDATED BELOW I have 9 shapefiles, each shapefile has 700,000+ rows and I have roughly 12 attribute fields in each shapefile. There are three attribute fields, common to every shapefile, that I am ...
4
votes
2answers
113 views

Converting shapefile to raster: Character attributes

I've got a shapefile that I would like to convert into a raster, but the trick is that the attribute I want a raster version of is character data. I'm familiar with geoprocessing using GDAL (in the ...
0
votes
1answer
81 views

How do I generate a correlogram from a GeoTiff using R?

I'm trying to evaluate the scale of spatial autocorrelation for a bunch of rasters I have as GeoTIFFs. I'd like to be able to do this entirely using R, in order to maintain smooth workflow. What's ...
0
votes
1answer
39 views

Point-sampling rasters using R? [duplicate]

I need to be able to sample many raster files (imported from GeoTIFFs) using a previously defined set of spatial points, using R. What's the best way to do this?
3
votes
2answers
336 views

Extract Raster from Raster using Polygon shapefile in R

I am new to R and using the raster package. I have a problem extracting polygons from an existing raster file. If I use extract(raster, poly_shape) function on the raster it always creates a list ...
4
votes
0answers
27 views

Projecting polygon when calculating its area? [duplicate]

I have a broad range of polygons (in vector format and in SpatialPolygonsDataFrame class), some going east-west, others going south-north, while even others span polar regions or very large ranges as ...
1
vote
2answers
94 views

How to create a time series with a NetCDF raster from Ensembles E-OBS using GDAL, R, GRASS GIS or SAGA GIS

I downloaded a NetCDF file that contains daily precipitation values from 1-1-1950 to 30-6-2013. The raster has 23191 bands. http://eca.knmi.nl/download/ensembles/download.php My goal is to analyze ...
2
votes
1answer
76 views

Compute the area of a given polygon (in vector format) lying between specified lines of latitude

I have species distributions (in vector format) and I would like to calculate how much of the area for a given species distribution is within certain lines of latitude, such as temperate and tropic ...
1
vote
1answer
51 views

How to deal with empty corners in raster when calculating fragmentation index

I have clipped buffer area around coordinate from one big raster (forest cover). Such created raster was loaded in R where PatchStat function (library SDMTools) gives fragmentation index. After ...
2
votes
1answer
52 views

Plotting spatial data when two spatial objects have different CRS

I have a spatial polygon object and a spatial points object. The latter was created from xy latlong data vectors (called latitude and longitude, respectively) in a dataframe, while the former was ...
2
votes
3answers
102 views

R - color SpatialLines by attribute

I have a SpatialLinesDataFrame dataset of rivers in R. The dataset consists of lines (river segments) with a "Turbidity" attribute. I want to create a map and legend, where: Turbidity between 0 and ...
1
vote
1answer
55 views

How to infer a relationship between two spatial points using R

I am quite new to the area of spatial statistics, but I'm very interested. For learning and demo purposes, I've created three datsets. Dataset - Persons: This describes individuals at a certain ...
7
votes
1answer
121 views

How can I conduct Geographically Weighted Analysis using an exponential decay function?

I have some points (representing 30 study sites) and I want to calculate the weighted mean of several variables (landscape layers) using a negative exponential decay (weights function) to give more ...
0
votes
0answers
22 views

formula for GWPCA cross validation on GWmodel package on R software

i want to know the formula of for GWPCA cross validation to get the best bandwuth on GWmodel package on R software. thanks
6
votes
1answer
293 views

Route analysis in R based on the OSM transport network: possible?

I'm analysing the likely paths taken by commuters on their way to work, based on Census flow data, job density distributions and the transport network. At the moment I'm using pgRouting and the ...
5
votes
1answer
758 views

Shortest distance from point to line (big dataset)

What is the best method to calculate the distance from the nearest road for a series of address points in QGIS, having a very big dataset? (QGIS 2.0.1, Win XP) I read the answers to previous similar ...
4
votes
1answer
107 views

Is it possible to use RPy/RPy2 in a Geoprocessing Script that is exposed as a GP Task on ArcGIS Server?

Is it possible to publish a Script tool that uses the RPy/RPy2 module (or more generally any Python-to-R link) to ArcGIS Server 10 SP2? If so, what specifically do I need to do, i.e. Install R on the ...
1
vote
3answers
120 views

How to access shapefile table in QGIS Modeler R scripts?

I need to develop .rsx scripts for QGIS 2.0.1 Modeler to do statistical elaborations on data tables of loaded vector layers. In this case I am not interested in accessing shapefile geometries, I'd ...
0
votes
0answers
46 views

Failure to display R script result plots in QGIS

I'm trying to run R scripts (3.0.2) from QGIS (2.0.1) on Windows 7 x64. There seems to be a known bug whereby the output results do not display in the QGIS GUI, but rather just an html path ending in ...
2
votes
1answer
146 views

MODIS Level 1B batch process

I am a new MODIS Level 1B user. I have been using the MODIS Swath tool to process my images. However, I have over 1000 images to process. I was wondering if there a script available in IDL or R (very ...
2
votes
0answers
84 views

R raster Package Moran's I interpretation

I am using raster package to calculate the local Moran's I . The example gives the range of Moran's I between - 1 to 2.47. On my own data I have see the value range -3.070423 - 7.228558 How can ...
3
votes
3answers
318 views

Convert 800 .hdf4 files to Tiff using GRASS/QGIS/SAGA/GDAL/R

I've have MODIS ndvi data (MODIS_Grid_16 days NDVI) from the iberian peninsula, from February of 2000 to August of 2012 (really big data). How do I batch convert about 800 .hdf MODIS files to .tif ? ...
2
votes
1answer
77 views

Plot a Raster in R with white space for NA values

I want to produce a .png of a raster with background values plotted as white space. Below you can see the current state of the plot. As you can tell, the pixels which should contain NA values, ...
2
votes
2answers
69 views

R - What is the purpose of .GAL lattice file?

Could someone please explain to me what the purpose of a .GAL lattice file is? I want to see if my data are spatially autocorrelated and I've been following Bivand (2013). I have loaded my polygon ...
3
votes
0answers
106 views

R Raster package loop over pixels

I want to perform a trend analysis for each pixel over it's time dimension. Therefore I wrote a loop over each pixel doing the calculation but this is extremely slow. I am sure there must be a faster ...
6
votes
2answers
846 views

Resources on using R in QGIS for R users?

This is a very basic question, but I'm having a heck of a time figuring it out because so much of the QGIS documentation is oriented toward use of the GUI. I want to be able to run R scripts in QGIS, ...
2
votes
1answer
56 views

Matching two point clouds without loss (nearest neighbor, but more complicated)

I have two point clouds, one with demographic and one with electoral data. Points in both clouds represent the same empirical entity (say, a village, or an urban neighbourhood), but are not in the ...
0
votes
0answers
46 views

Resources for mapping and displaying spatial data in R

I'd like this to be a repository for all resources available specifically for mapping and displaying spatial data in R. I've seen (and contributed) to this post, but would like this post to focus ...
23
votes
6answers
6k views

Tutorials to handle spatial data in R

Does anyone know some good tutorials to handle spatial data in R? I know that there are some packages like rgdal and maptools but I haven't actually found usefull tutorials to get in touch with the ...
1
vote
1answer
140 views

R - Dealing with missing data in SpatialPolygonDataFrames for Moran.test

I'm trying to use Moran.test on a SpatialPolygonDataFrame consisting of 7194 elements in R. I know that there is around 150 polygons with NA values. First I generate a spatial weights matrix: ...
1
vote
1answer
448 views

How to process MODIS ndvi data using QGIS/GRASS/R?

I have MODIS ndvi data (MODIS_Grid_16 days NDVI) from the iberian peninsula (3 tiles), from February of 2000 to August of 2012 (really big data). I want to know how to deal with this kind of files. My ...
2
votes
0answers
191 views

Mapping Geocoded addresses onto Census shapefiles in R

I have a list of addresses that I have geocoded to get the corresponding longitude and latitude. > AddressDetail <- geocode(match_address) > coordinates(AddressDetail) <- c("lon", "lat") ...
0
votes
0answers
87 views

Defining the boundary of a catchment area

I digitized contour lines from topographic maps and created a DEM. With the DEM, I estimated the boundary of a catchment area. Estimated positional accuracy of topographic maps is RMSE = 65 m, with a ...
5
votes
1answer
2k views

R : convert a spatial polygon objet to spatial polygon data frame

My goal is to modify an existing shapefile by merging certain polygons. After importing the shapefile and using the UnionSpatialPolygons command, I get the polygon outline that I want. However, this ...
6
votes
3answers
464 views

Overlay a spatial polygon with a grid and check in which grid element specific coordinates are located

How can one use R to split a shape file in 200 meter squares/sub-polygons, plot this grid (incl. ID numbers for each square) over the original map below, and evaluate in which square specific ...
8
votes
3answers
209 views

R, rasterVis and RGL: how to prevent reflective surfaces?

I have been experimenting with digital terrain model files in ESRI ASCII grid raster format. Despite zero experience with this type of data I found it very easy to load into R, as per code below. The ...
2
votes
2answers
496 views

R Raster package: write NetCDF with time dimension

I am writing a rasterbrick with multiple layers to a netcdf file using the raster package in R: writeRaster(gridfile, "Data", "CDF", overwrite=TRUE, ...
2
votes
0answers
138 views

Per-pixel calculations on a raster stack using GDAL

In R, it's relatively trivial to perform per pixel calculations based on a raster stack (e.g., get std.dev for each pixel on a 12 layer GeoTIFF). Unfortunately, the speed is less than desirable when ...