Statistical computing language and software environment.

learn more… | top users | synonyms

2
votes
0answers
19 views

How to analyse point data as covariate for areal data

I have two data sets: a) a set of polygons with areal data variables. b) point measures of various co-variates of interest (note I do not necessarily have a point for each polygon) How do I ...
2
votes
2answers
256 views

Iterating polygons on a SpatialPolygonsDataFrame in R with readWKT

I have fetched some data (39 objects) from my PostGIS database and converted this data into a matrix: rs <- dbSendQuery(con, "SELECT ST_AsText(polygon_accumulation) AS ShapeWKT from ...
0
votes
1answer
74 views

Processing HDF4 data in R

I am trying to process MODIS files in R, which are in HDF4 format (link to data). I have experience reading in HDF5 and netCDF files in R, but HDF4 seems to be a common problem in R. I am using a MAC ...
2
votes
0answers
11 views

Error in unionSpatialPolygons using R

I am migrating from ArcMAP to opensource GIS. I've succesfully installed R and many packages in a Ubuntu 12.04 machine. I am used to linux and programming in C, but not to R and its syntax. I am ...
0
votes
1answer
36 views

Conducting geographically weighted regression on a raster

I have a relatively big raster file on which I want to conduct geographically weighted regression. Size of raster ~ 2 GB and number of rows is around 100,000 One way might be to convert from raster ...
0
votes
0answers
11 views

error in using akima package in r for irregular point data interpolation

I am tring to interpolate my irregular Point data with an Elevation value. I found akima package for Interpolation and tring to use it but it Shows error message like this: Error: could not find ...
2
votes
2answers
47 views

Select polygons that do not intersect any other polygon using QGIS, GRASS, SAGA or R

I would like to select the polygons of the layer below that do not intersect with any other polygon (within the same layer). How can I do this using QGIS, GRASS, SAGA or R?
0
votes
1answer
22 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
6 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
17 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
11 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
26 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
6 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
16 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
22 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 ...
1
vote
1answer
88 views

Interpolation with netcdf4 file

I have a netcdf4 data set that contains array and information about latitude and longitude which can be accessed with the following R command: ncvar_get(dataname, 'lon'); ncvar_get(dataname, 'lat'); ...
4
votes
1answer
150 views

How to compare spatial patterns (crime cases vs population)?

I am looking for a way to statistically compare two spatial patterns in shape polygon data (not point pattern analysis). Let's say I have a city map divided in neighborhoods (barrios), in shape ...
0
votes
1answer
35 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
28 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
23 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 ...
0
votes
0answers
11 views

display followers count on google maps using R

How can I plot the count of followers, locations wise on google maps- Im using R language my data is as below for(i in 1:25){print(friends_obj[[i]]$location)} [1] "London & Glasgow/Edinburgh, UK" ...
6
votes
2answers
179 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
109 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 ...
1
vote
1answer
399 views

GDAL loads in Python, but not IPython (Linux)

I am trying my hardest to utilize the work of the good folks at the GeoDa center, but I cannot seem to get all of the components to load. GDAL, in particular, has been giving me some difficulty. The ...
4
votes
2answers
87 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
30 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
28 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?
2
votes
2answers
81 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
24 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
55 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
41 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
44 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
32 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
56 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 ...
0
votes
1answer
163 views

How to access each raster in a stack with a loop

I have a stack of rasters and I want to access each layer, combine each later with another file and finally export the resulting file as a csv file. Is there any way to do this with a loop? So far I ...
1
vote
1answer
41 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
92 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
13 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
4
votes
1answer
242 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 ...
1
vote
1answer
71 views

How to analyze and combine overlaping elevation maps into one?

I have 30 overlapping maps that contain tree density heights (measures between 1-100%) of a large area, and I need to get the highest point by pixel and create a new map with all these highest points ...
4
votes
1answer
391 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 ...
0
votes
1answer
138 views

Calculate area covered by a particular raster layer within polygons in QGIS

I am trying to extract information on vegetation cover within a set radius of point locations (my sites). I have succeeded in creating a new vector layer of circles with a 5 km radius around my sites. ...
4
votes
1answer
106 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
85 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
26 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
130 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
46 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
148 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 ? ...
5
votes
2answers
63 views

join grid or polygon layers in R

I have two grids. Cells in each grid have an attribute. Size of grid cells differs between the two grids and the two grids do not overlay each other perfectly. I would like to create a new layer, ...
2
votes
1answer
34 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, ...