Statistical computing language and software environment.

learn more… | top users | synonyms

0
votes
0answers
12 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
0answers
21 views

Raster data in R to merge to a time series

I have gridded data files about 2000 ( see sample -each file consist of monthly rainfall and is label pYYYYMM) https://www.dropbox.com/sh/63z166tjxyu12s5/AAAs3Ccn1zdVoBYMj8Y1o303a?dl=0 ...
1
vote
0answers
22 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 ...
-1
votes
1answer
22 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 ...
0
votes
1answer
29 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 ...
0
votes
3answers
30 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
1answer
61 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 ...
0
votes
1answer
29 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
136 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
79 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) ...
0
votes
0answers
32 views

Comparing point pattern distributions using Geostatistics with R?

I have a data set consisting of about 1,500 points of an unknown distribution. About 60 of the points represent special conditions like infekted people. How can I determine if the 1,500 points and ...
2
votes
1answer
50 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 ...
0
votes
0answers
35 views

Convert Eps to geographic distance using DBSCAN

I am using the dbscan cluster (package fpc) in R to find clusters on a set of latitudes/longitudes coordinates. I want to find an eps distance that corresponds to a meaningful geographic distance ...
0
votes
0answers
24 views

weight a raster by another for a model ( R/ arcgis/maxent)

I have an idea in mind for my project but I can't figure out how to do it in R and Arcgis. Basically I would like to do a weight a raster by another for my maximum entropy model. -I got a raster ...
0
votes
1answer
17 views

TopologyException error during fortify SpatialPolygonsDataFrame

I tried to reproduce a recent question on spatial data visualization using ggplot2 and received the following error message during fortify(world, region = "ISO3"): Error: TopologyException: Input ...
2
votes
1answer
29 views

Perform affine transformation on shapefile with OGR, GDAL, or other command line tools?

Transformations with WKT or EPSG code are made easy with ogr2ogr, but are there tools for me to define my own affine transformation to apply to the coordinates of a shapefile? Or, in general, is ...
0
votes
1answer
30 views

Add Polygons to SpatialPolygons via loop / iteration using R?

Just out of interest: I'm creating a SpatialPolygon (R with package "sp") with a number of Polygons within. Is there a way to loop this task? I'm doing it manually as follows. The Polygons consist of ...
2
votes
1answer
41 views

Extracting point information from raster in R

I'd like to work with climate data available from here and extract values of certain locations (for which I have lat and lon). This is what I've done so far: library(SDMTools) library(raster) test ...
1
vote
2answers
54 views

How can I loop through a list to create rasters and write them to .tif files?

I received 3 .mat files that each contain several daily maps. I need to take them and get them into ArcMap (10.2 if that matters) as a .tif with spatial information. I have already used the R.matlab ...
0
votes
0answers
36 views

Rasterize and union problem

i've done the union first and then rasterized it. But now i want instead of that, to rasterize each one and then merge it as one. How do i do that? I mean without the union first off. ...
2
votes
1answer
78 views

How can I parse modis MOD13Q1 quality layers in R?

I'd like to parse MOD13Q1 VI quality layers in R. My original code must have been incorrect, because it gave lots of "deep ocean" flags in the middle of Brazil. After the changes suggested by ...
2
votes
2answers
152 views

Classifying DEM according to elevation in R?

r<-raster("DEM.tif") How to do the classification based on every ten meters for example and have different colors too?
0
votes
0answers
61 views

Error in downloading MODIS image using R script

I am trying to download "MOD13A3" MODIS image using R script (http://r-gis.net/?q=ModisDownload) but each time it's ending in error. How to make it work ? ...
3
votes
1answer
44 views

Plotting data without geographical coordinates besides the main map in ggplot2

I am trying to plot some data as a bubble chart on the world map using ggplot2 package in R. I am able to proceed as follows. # Load necessay packages library(ggplot2) library(rworldmap) ...
0
votes
0answers
26 views

Location allocation optimization using R?

I need to allocated zip codes to offices. Each zip code has certain demand from 0 to 120 clients and each office has certain capacity from 500 to 2500. How do I do this so that: All the zip codes ...
2
votes
1answer
75 views

Difference between raster() and readGDAL()?

I am new in working with rasters and after having trouble working with them i'd like to know the difference between those two import functions. for example i did this library(raster) ...
3
votes
1answer
91 views

Reading and plotting shapefiles with a custom function in R?

I have this piece of code that want to turn into a function r<-readShapeSpatial("roads") b<-readShapeSpatial("buildings") lc<-readShapeSpatial("landcover.shp") here is the preparation ...
1
vote
1answer
32 views

error of randomforest classification with Landsat data

I want perform randomforest classification on Landsat data. Here is some code I found from a paper. But when I run into the "classification" part, there is an error said "Error in eval(expr, ...
1
vote
0answers
41 views

Matching (joining) two slightly incompatible polygon shapefiles

I have two shapefiles that represent two slightly different zonings of the same area (municipality zoning and transport zones), about 3000 polygons each. Several smaller municipalities may form a ...
1
vote
2answers
29 views

Include non-raster object in overlay function in R

I have two raster objects and I would like to calculate values of the first conditional on values in the second. I am currently using the overlay function, but unfortunately I can't seem to pass non ...
2
votes
1answer
52 views

Defining projection using readOGR if the shapefile doesn't have it?

I think that by using readOGR to import a shapefile renders worthless to put proj4string to the script since if any coords are available in the shapefile, it will be transported through the readOGR. ...
1
vote
1answer
66 views

encoding error between maptools and arcgis

I am having a trouble of encoding at writing shape file by maptools (0.8-34) of R by MacOSX and reading that by ArcMAP 10.2.1 in windows 8.1 through parallels desktop. Using an original shape file, I ...
1
vote
0answers
56 views

What is the most user friendly OSM route planning API?

I ran a workshop last week on GIS for Transport Applications. In Nick Bearman's tutorial (available from the previous link) we used Google's routing API (via ggmap) and CycleStreets.net routing ...
0
votes
0answers
53 views

How to set raster projection in R

I am new to R and grib files. library(raster) I have a "grd" file in and I am trying to set the coordinates to polar stereographic projection r <- raster("p190001.grd") # class : ...
2
votes
0answers
37 views

How to create geometric shapes in R? [closed]

I want to plot some geometrical shapes but I dont know what values to put in a vector so they can be plotted. How do you know what numbers should be inserted in a vector so they can create certain ...
0
votes
1answer
35 views

Rasterize gives weird image

ro<-readShapeSpatial("roads") ra<-readShapeSpatial("railways") co<-readShapeSpatial("LIE_adm0") u<-union(union(extent(ro),extent(ra)),extent(co)) rast<-raster(u,res=1) ...
1
vote
1answer
25 views

Shapefile One-to-Many Relationship in R

Is there any way to append a shapefile in a one-to-many relationship with its attributes that R can read? I understand the fix that resides in ArcInfo by creating through MakeQueryTable, but as this ...
1
vote
0answers
52 views

Put output of SAGA algorithm of Processing model in R script

I'm trying to combine the Variogram Cloud model of SAGA with an R script inside a Processing model. First, the R script is: ##Vector processing=group ##showplots ##Layer=table ##X=Field Layer ...
0
votes
0answers
46 views

How to get time-zones from latitude and longitude in raster?

I have a raster with longitudes and latitudes for the US with meteorological data (for about 100k grid points). I would like to get the time zone for each point in order to transform the time to local ...
1
vote
0answers
54 views

How to avoid “Unable to compute a transformation between pixel/line”?

I am trying to re-project from Equi-rectangular to the NASA global EASE grid (Equal Area Scalable Earth) grid 1383 586 with a ~25 km cylindrical projection (ESPG=3410). The file is provided as HDF5. ...
3
votes
2answers
54 views

How to compute layer by layer after stack?

I have some rasters that I want to compute the moving average and then rewrite the rasters again. r1 <- r2 <- r3 <- r4 <- r5 <- r6 <- raster(nrows=10, ncols=10); r1 <- ...
2
votes
0answers
35 views

How to re-grid a raster using area-overlap weighted averaging?

I have two rasters r and r2 with different resolutions. I want to re-sample (re-grid) r2 to be the same as r using area-overlap weighted averaging:If an r grid cell is covered by more than 1 ...
0
votes
0answers
25 views

How to add a plot with a different x axis to an existing plot?

I am plotting a raster but I want combine to things in the same plot but in two different axis. library (raster) r <- raster(nrows=10, ncols=10);r <- setValues(r, 1:ncell(r)) extent(r) = ...
1
vote
1answer
60 views

How to compute the number of none NA-values for each row?

I have a matrix(raster) that I am computing the the mean of each row in this raster as: library (raster) r <- raster(nrows=10, ncols=10);r <- setValues(r, 1:ncell(r)) # The x-values will ...
0
votes
1answer
59 views

Using R to combine long/lat data into UK region polygons

I've a CSV file with geoinfo: JobSeekerId LocationId UpdatedDt EndDt Latitude Longitude Country Name JobTitle CompanyName ...
0
votes
0answers
35 views

Transforming GIS based XY data (polygons) so that they overlap in R

So I'm attempting to come up with the average shape of these polygon features as a part of my masters research. To do this, I want to run a shape analysis of GIS-based X-Y data using the shapes ...
3
votes
3answers
188 views

R: How to get latitudes and longitudes from a RasterLayer?

I am an absolute beginner of geographic data, so please, forgive me if the question is not appropriate. I downloaded data from NCDC NARR and managed to load into R using the raster package. I would ...
2
votes
1answer
57 views

How to extract values based on names of rasters?

I have a raster r that contains values for each grid cell (which vary from one cell to another). I have other 24 rasters where the names of the rasters are: the name of r1 is data-00.img, the name ...
2
votes
1answer
39 views

Georeferencing in R Based on the Four Corners

I have a raster image (photograph) that was taken from a helicopter. Along with that raster image, the system also tells me the positions of the four corners of the image in a different file. I have ...
0
votes
2answers
43 views

Count objects in polygons

I have a .txt file containing the coordinate of the two polygons in the picture below as well as coordinates for each of the lines inside the polygons. Could anyone help to count each line inside ...