Statistical computing language and software environment.

learn more… | top users | synonyms

1
vote
2answers
116 views

Getting attribute table from WFS data source

I am attempting to get some WFS 2.0 INSPIRE data from the Finnish meteorological institute into a table format of some kind. Is there any software which can easily grab the attribute data for the WFS ...
0
votes
1answer
34 views

Calculate mean for a each object of SpatialPointsDataFrame

I have a list with several SpatialPointsDataFrame objects, and each object has certain slots. How can I calculate the meann value of a certain slot for every object in the list, and then save that ...
0
votes
1answer
63 views

Compare two rasters in R with the cell value from one replacing the cell value from another

I am trying to create a single raster from two rasters. The first raster is made up of all 0 values and the second is made up of the current date value (formatted in days from 001 to 365). The problem ...
0
votes
1answer
497 views

IDW in R with cross validation

I want to perform IDW interpolation using R, and more specifically, using the idw command from the gstat package. Also, I want to use the optimal number of neighbors (point locations) and distance ...
0
votes
2answers
80 views

GlobCover 2009 dataset - importing legend and colormap values into R

Does anyone know how to import a legend for a raster layer (loaded as a .tif) into R such that it can be attached to the raster file itself and used in plotting? In this case, I am using the GLC2009 ...
0
votes
1answer
315 views

Compare two rasters with raster overlay and replace cell value with variable

I want to compare 2 rasters and when there is a change from 2 to 1 create a new raster that replaces the cell where the change has been detected with the change variable. ## Finds Cells that have ...
4
votes
0answers
48 views

Why does @data@values slot for RasterLayer only contain logical(0) and not the actual values?

Trying to get to the bottom of why, when I read in a raster of NDVI, the @data@values slot does not contain the actual values until I set them manually. For example: NDVI <- ...
3
votes
1answer
119 views

Processing large amounts of Rasters in R

I have 6500 processed and clipped rasters that indicate snow levels on a daily basis in the northern hemisphere. I need to create a raster that represents a sum of the snow for the season. The winter ...
1
vote
2answers
62 views

Aggregate SpatialPoint to form SpatialPolygon

Is it possible to convert SpatialPointDataFrame to SpatialPolygonDataFrame by some kind of generalization or aggregation function in R. how can i form polygons from points of similar color/classified ...
1
vote
1answer
102 views

Trying to limit my shapemap to state boundary lines in R

I have a feeling I'm either trying to do this incorrectly, or I'm just not searching for the correct terms, but I cannot figure out how to do what I want to do. I am trying to plot the ranges of fish ...
1
vote
1answer
28 views

Workflow for updating raster using polygons of different attributes

I require a specific output but it is proving difficult for me to work out a workflow. I have: 1) polygons of new housing developments (356, non overlapping) 2) a 25m raster of land cover I want ...
2
votes
1answer
41 views

Last SEXTANTE (Processing) update results in missing RScripts examples

I have installed today the new version of the Processing plugin on QGIS 2.8 (on Windows 7 x64) and I cannot anymore find the RScripts example that were available (and functioning) before the udpate. ...
2
votes
1answer
27 views

Mosaic rasterstacks based on minimum of certain layer?

I am trying to mosaic 42 remote sensing rasterstacks (with 250 bands) based on the criterion that in overlapping areas, the pixel should be taken that has the most nadir viewing angle Beside my ...
1
vote
0answers
30 views

Using errorsalm() from spdep with multilevel data?

I have a dataset that has multiple observations per county, and I am interested in estimating the effect of a county-level variable on individual-level outcomes. I would like to account for adjacency ...
0
votes
1answer
150 views

Switch an R script from “MODIS Reprojection Tool” to “ModisSwathTool”

This R script works with MRT (MODIS Reprojection Tool), and its performance is unbelievable for land data. The only problem is when you use Atmospheric data from Modis, like "MOD05_L2". The problem ...
1
vote
3answers
633 views

R Spatial Kernel Density Estimation

I am new on spatial kernel density estimation with r and need some suggestions. Saying, I like to estimate the density for some event occurring at a location, for example, the probability of ...
2
votes
2answers
149 views

Use GRASS through R

I want to use GRASS GIS 7 through R 3.1.3 in kubuntu operational system. Could someone help me connect GRASS with R, and use GRASS commands through R (not the other way round!).
0
votes
2answers
170 views

How to get average values per coordinates in ArcMap/R?

I have a point shapefile. It has multiple points at different spatial locations. Within each point there are multiple values that I want to get the average of. In this case the new file would have one ...
1
vote
1answer
51 views

How to create a Planer Point Pattern in R using a Shapefile of point features

I have just started to learn spatial analysis of point pattern with R. On page 20 of this book, the "Swedishpines" dataset has been used for inspection purposes. However, I want to use my own dataset ...
0
votes
0answers
39 views

Total area of polygons over latitudes

I'd like to calculate the total area of several polygons along a north-south axis. A bit like this http://www.closr.it/canvas/1307/#/spot There is a similar question with R : Spatial graphing, ...
0
votes
0answers
30 views

how to assign the same centroid to polygons grouped by ID , R

I have a SpatialPolygonDataFrame object from a shapefile of different zip codes in New York state, many zip codes are represented by several polygons( for different reasons) so I want to assign the ...
1
vote
0answers
97 views

Simulate water runoff to get water height over each pixel of DEM

Here is my problem that I am trying to deal with. I have given a height map of around 512X512 resolution(It is Geo-referenced and I have the meta-data for this). And I have given that this area ...
1
vote
0answers
141 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
0answers
18 views

Creating Concentric circles using the ggmap function from GGMAP package R [duplicate]

I would like to create two concentric circles using the ggmap package and ggmap function, one being 1/4 mile and the other 1 mile radius around the 16 metro rail stations in Los Angeles California. It ...
1
vote
0answers
68 views

Effects of heterogeneous vegetation in mapping canopy cover using randomForest in R

I have a 5.8 m resolution satellite image (LISS IV by IRS-P6, in bands 4, 3 and 2) of a densely forested region in the Eastern Himalayas. My objective is to predict regions of canopy cover. Regions of ...
1
vote
1answer
116 views

The stability of randomForest models after increasing predictor variables

With reference to the posts: Incorporating terrain data to predict canopy cover using randomForest in R and Random-Forest Classification of 10cm Imagery for species-distribution in R (no point-shapes) ...
2
votes
1answer
118 views

Incorporating terrain data to predict canopy cover using randomForest in R

This is a subjective question. Thank you bearing with me. I have a 5.8 m resolution satellite image (LISS IV by IRS-P6, in bands 4, 3 and 2) of a densely forested region in the Eastern Himalayas. My ...
1
vote
0answers
41 views

Assign Z value to new coordinates based on interpolated response surface in R

I've created a response surface of sorts based on interpolated z values (bottom of post). I'm looking for the most efficient way to assign the respective z value to a new point (actually, a dataset ...
2
votes
2answers
234 views

gDistance to determine closest point to polygon: does changing projection change results qualitatively?

I'm using gDistance from the rgeos package for R to find the nearest point to a polygon. The projection of both points and the polygon is +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 ...
0
votes
0answers
45 views

Creating a for.loop to randomly extract data from a raster stack within shapefile extents

I'm currently comparing a MODIS time series of about 15 years to determine the effects of pasture management practices on the "greeness" of pastures in a water catchment. So far I've been able to ...
2
votes
2answers
129 views

R apply function in turn to each subset of a spatialPolygonsDataFrame

I'm trying to get the area of all the polygons in a SpatialPolygonsDataFrame: parcProj. library("rgeos") gArea(parcProj[2,]) works fine but: apply(parcProj,1,gArea) or ...
1
vote
0answers
19 views

weighting scheme of gstat for fitting models to variograms

I use R and the gstat package to construct variograms of temperature, and I would like to ask If anyone knows what weighting scheme does the gstat package use to fit a model to the sample variogram?
0
votes
0answers
15 views

fit polynomial function on spectral data

I want to fit a polynomial function (max. 3rd order) on each raster cell over all my spectral bands (Landsat 1-7) creating a new raster(stack) representing the coefficients. I got my data (including ...
3
votes
0answers
201 views

Python vs R for analzing spatial data [closed]

I am relatively new to both the python and R programming languages, having a small amount of exposure to both. Which language has the edge when working with spatial data and when would you use R over ...
2
votes
0answers
62 views

How to clip or erase a csv file with a non-rectangular shapefile?

I can clip a csv file with a ogr2ogr having 4 extents of a rectangular polygon. Is it possible to clip a csv file with a non rectangular shapefile?
2
votes
2answers
259 views

How to make a movie from a ascii grid files?

Lets say I have several ascii grid files for each time step for a particular location. I would like to spit out the image for each of these ascii grid files and legend on the side bar. I also would ...
0
votes
0answers
50 views

Merging two regions of several into one, to get output SpatialPolygonsDataFrame file

I've been using the large SpatialPolygonsDataFrame file with outdated border of regions. I want to correct boundaries, but also hope to get in the end Polygons with Data Frame again. I used the ...
12
votes
2answers
506 views

How does spatial polygon %over% polygon work to when aggregating values in r?

I'm working on an environmental epidemiology project where I have point exposures (~ 2,000 industrial hog operations - IHOs). These IHOs spray on nearby fields, but the feces water droplets and smell ...
0
votes
0answers
21 views

R to ArcMap, erroneous change of spatial reference

I have SPOT5 imagery from Northern Sweden / Finland. Data is downloaded from SACCESS service and is provided in Sweref99 TM projection. When opened in ArcMap the data is shown in Sweref99 TM as ...
0
votes
1answer
202 views

Extracting point data for each lat and long from raster file into a dataframe

I have gridded data files ( see sample -each file consist of monthly precip and is label pYYYYMM) https://www.dropbox.com/sh/63z166tjxyu12s5/AAAs3Ccn1zdVoBYMj8Y1o303a?dl=0. I am trying to extract ...
1
vote
0answers
117 views

moving window regression on raster stack in R

I want to perform a moving window regression on every pixel of two raster stacks representing Band3 and Band4 of Landsat data. The result should be two additional stacks, one representing the ...
0
votes
0answers
95 views

How can I convert a geojson file to topojson in R?

I have tried the geojsonio package but it does not allow me to convert my geojson to topojson format. I know there are web applications available online to convert from geojson to topojson, but my ...
1
vote
1answer
141 views

Function to merge geojson and csv data in R

I have been using merge() to merge my geojson file with csv data in R. My codes are as follows: map@data <- merge(map@data, csvdata, by.x="ID", by.y="ID") However, after exporting the geojson ...
3
votes
3answers
670 views

Unable to read raster into R: 'TIFFReadDirectory:Failed to read directory at offset'

I'm getting a weird error when trying to read a raster into R. >r <- raster('myraster.tif') Error in .local(.Object, ...) : TIFFReadDirectory:Failed to read directory at offset 2606110 ...
3
votes
2answers
331 views

Strange spatial interpolation results from ordinary kriging

I'm using Ordinary Kriging to interpolate a set of spatial data. I'm working on different sets which represent different days of measurement. I get good results, apart from one specific set. Will post ...
11
votes
2answers
316 views

Point overplotting on maps in R (e.g. with ggplot2) - how to push points to the side, mark original positions, combine close-by points …?

I've got a database containing languages, their longitudes and latitudes and a feature value (either category 1, category 2 or both - in the plot these are marked red, blue and green respectively). ...
0
votes
1answer
320 views

How can I export a GeoJSON file from R?

The write. function allows one to export a JSON file, but it does not seem to work for a GeoJSON file. This is the code i ran: write(map,"Map Layer2.GeoJSON") and this was the error produced: ...
1
vote
3answers
635 views

How to calculate mean aspect in circular statistics?

Measurements of angles (e.g. compass bearings, aspect) are circular in nature (e.g. 0 and 360 degrees coincide) and therefore cannot be meaningfully described by ordinary arithmetic means. I have a ...
1
vote
3answers
158 views

Spatial distribution of an invasive specie across a landscape : what analysis?

I'm studying the spatial distribution of an invasive specie introduced in a natural area. This plant, that thrives in sunlight, grows in thickets (patches). The population of my study site (10x10 km) ...
3
votes
1answer
333 views

How do I set up the R Script folder in QGIS 2.8.1 Processing Toolbox (Advanced Interface)?

I am currently lacking the R Script folder in my processing toolbox under the advanced interface tab. I suspect this could be a result of the firewall in my company's wifi. I already have R ...