Statistical computing language and software environment.

learn more… | top users | synonyms

2
votes
1answer
680 views

Distance to nearest point for every point same SpatialPointsDataFrame in R

I have a SpatialPointsDataFrame that I'm working with in R. I want to add a new vector to the dataframe that contains, for every point, the distance to the nearest other point in the ...
0
votes
1answer
13 views

GDAL error in R path.expand(path): invalid 'path' argument

I would like to stick to using R since that's where a lot of my pre-processing and analysis will take place for a bigger project. I have downloaded a number of MODIS files for the MOD13A3 product, ...
-1
votes
1answer
48 views

Clustering of Spatial Data in R or Python

I am new to GIS. I have a large set of (latitude and longitude)spatial data. I am confused on which clustering method to adopt. I have came across two density based approach: DBSCAN and OPTICS. I am ...
0
votes
0answers
7 views

Using errorsarlm with two different listw objects

i am trying to estimate some spatial error models with "errorsarlm" from the spdep package. In my first regression everything is fine: sreg <- errorsarlm(map@data$avginc ~ ...
0
votes
2answers
15 views

R - Easy way to reproject an extent object

I have an extent object that I would like to re-project to lat and lon. >outer_extent outer_extent class : Extent xmin : -5559753 xmax : -4447753 ymin ...
0
votes
0answers
20 views

R-Studio0.99, How to compute correlation between two rasters and serial correlation

I am using R.Studio0.99 version installed on my HP Laptop 4th generation (8GB RAM) windows 8 operating system and I know how to use R studio and some of the library to incorporate spatial data ...
1
vote
1answer
25 views

“Hole” argument ignored in call to Polygon function

I am working with a dataset consisting of a list of coordinates for a set of polygons. Some of these polygons are holes. However, the dataset does not follow the convention of listing the coordinates ...
3
votes
1answer
25 views

Plot a world heat map with mollweide projection in R

I have a data set with global observations with lat/lon information. There are 288 longitudes from -179.375 to 179.375 with spacing 1.25, 180 latitudes from -89.5 to 89.5 with spacing 1. There is one ...
0
votes
0answers
11 views

Identify adjacent and no-adjacent polygons between multiple shapefiles using R?

My spatial data are recorded multiple years, providing one polygon for a year. (Polygon is derived from raster data, so I can use rasters instead if it will be easier to run my calculations). I ...
1
vote
0answers
19 views

raster::intersect not producing expected results

I am having trouble with R's raster::intersect method. I give intersect a SpatialPolygonsDataFrame and SpatialPoints as documented in the raster reference manual (p. 111). The result acts like it ...
2
votes
1answer
394 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 ...
0
votes
1answer
207 views

plot 3D graph based on raster pixel values

I am new to raster data. I have a raster tif file with population density value in each grid cell. I am interested in plotting a 3D graph based on raster values in ArcGIS or R. The similar graph is ...
5
votes
4answers
155 views

How to speed up raster to polygon conversion in R?

I've decided to process my Landsat data in R instead or ArcGIS - due to my missing knowledge of python and and because of (assumed) high computation capacities of R. I want to : import r1 raster to ...
2
votes
1answer
37 views

Getting long-lat coordinates from NARR LCC-gridded data using raster and rgdal in r

I am studying migration ecology of sea ducks on the Atlantic coast of North America. In my analysis, I require substantial environmental data over a large geographic range for the past 15 years and ...
0
votes
0answers
11 views

MODISTool R Error: Can't call method “prepare”

I have already used Modistool in R a few times, but lately I am facing the following error while using the function GetProducts(), GetBands(Product = "MOD13Q1") MODISTools::GetProducts() Error in ...
0
votes
0answers
26 views

R Levelplot overlaying two SpatialPolygonsDataFrame on a RasterLayer

I'm having some trouble plotting a RasterLayer and a SpatialPolygonsDataFrame in a level plot in R. Something is wrong with the projection, but I don't manage to understand what. Here you can find my ...
0
votes
1answer
28 views

Writeogr to PostgreSQL/PostGIS database with R

I'm producing some Spatialpolygondataframes with R that I would like to upload to my PostGIS Database. I've found out that the way to do this is via the writeOGR() function using the PostgreSQL ...
1
vote
2answers
56 views

Clustering geographical data based on point location and associated point values

Given data points with longitude, latitude, and a third property value of this point. How can I cluster points into groups (geographical sub-regions) based on the property value? I searched by google ...
1
vote
1answer
51 views

Creating a polyline with a specified slope

I would like to create a road (polyline) starting from a point, with a specified slope, until reaching the top of a hill (or going out of the DTM Raster). The dowstream slope of the DTM must always ...
1
vote
2answers
89 views

Leaflet output is grey

I would like to implement leaflet in shiny with R. I try to go with the simple examples I found online: library(leaflet) library(ggmap) somePlace <- ggmap::geocode("Vienna") # Information from ...
1
vote
1answer
27 views

Spatial correlograms for binary raster?

I have binary raster data (deforestation hotspot or not) that I would like to plot spatial correlograms for in R, to look at spatial autocorrelation over different spatial lags. I came across this ...
0
votes
1answer
61 views

Is there a R package that does spatial regularization?

By spatial regularization I mean the method described in Aitkinson and Curran (1995) "Defining an Optimal Size of Support for Remote Sensing Investigations", or Section 4.8 of Webster and Oliver ...
3
votes
2answers
30 views

Error to obtain coordinate reference system argument from a raster image in R

I got an error when i wanted to obtain the CRS from stacklayer. I don´t know what could it be. I need some advice or help to obtain it. This is my code. #LOAD RASTER LIBRARY >library(rgdal) ...
0
votes
1answer
246 views

How to convert list of stack images to time series files in R

I am using bfast package for time series analysis. But it requires the files to be read as time series. How to convert list of stacked images from 2000 to 2013 into a time series and given as input ...
2
votes
1answer
17 views

Read multiple layers of KML file using R

I have a KML that 3970 layers. I am trying to read it into R using the following script: library(rgdal) lyr <- ogrListLayers("file.KML") mykml <- readOGR("file.KML","Layer1") This works but ...
1
vote
1answer
28 views

rgdal / ggmap : rotate map data and clip map

I'm trying to plot a map of Toronto with ward boundaries and geocoded data on events over years. But the layout is such that I end up with mostly wasted space in the map. To increase resolution for ...
3
votes
2answers
103 views

Creating OUTSIDE_ONLY buffer around polygon using R?

I would like to create buffers in polygons, surrounding my original polygon. However, I want to keep only the "outside polygon", not include also original polygon. In Buffer in ArcGIS I can just set ...
0
votes
1answer
23 views

Goodness of fit of variogram model using geoR

How can I find a measure of goodness of fit (let's say an R2) of a model fitted to a variogram by the function variofit using the geoR package? The only value I have is the "minimised weighted sum of ...
0
votes
1answer
25 views

Rescale raster in R

What's the best way to rescale a raster with range 0, 200 to range -1 to 1 in R? I'm currently using gdal_translate, but this requires that I write out a new file with the updated range. If possible ...
4
votes
0answers
52 views

Is it possible to use a raster::click event within an interactive leaflet map?

Data Download: https://www.dropbox.com/s/f849ra8l6umo9d8/crops.tif?dl=0 I have various raster images of the world that I would ultimately like to plot within an interactive leaflet map; an example is ...
0
votes
1answer
19 views

Get value from list class (vector) in R [closed]

I have a simple question about corce(?) list class data into vector. I extract a tif file from polygon to achieve a list. Maybe the class of list is leaded by extract function. I check the list length ...
0
votes
1answer
51 views

How to read MODIS hdf file into R [closed]

I am having some problems with the processing MODIS HDF in R. I want to import MOD13A3 NDVI hdf file. Purpose is to read image cells value as a matrix. I have install MRT in my PC, I had planned to ...
3
votes
1answer
27 views

Extract rasterlayer from extent class

I used extract() function to get remote sensing image cell value.The same code to process tif file produce error message. Code as follows: r1 <- raster("2000_2605_1.1_km_monthly_NDVI.tif") ...
0
votes
1answer
67 views

R equivalent of GRASS r.grow?

What is a good R equivalent of the GRASS function r.grow? The closest I've been able to find is the script suggested by whuber in response to Creating randomly shaped clumps of cells in a raster from ...
11
votes
1answer
613 views

Multispectral image segmentation for natural-resources applications using R

R has the capacity for image segmentation, although all of the examples I have come across utilize a single band for the segmentation (example). I am interested in combining R's capacity for random ...
0
votes
1answer
206 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 ...
4
votes
2answers
145 views

Changing raster calculation values - shorten this script?

i'm doing something quite straightforward regarding raster algebra but am struggling either to find the correct function or make a function work correctly; Basically I have 2 rasters, representing ...
3
votes
1answer
83 views

Put points on a projected world map

With the help from this post Lines on reprojected sp objects with Mollweide projection, I draw a world map using Mollweide projection. I also have some points with latitudes and longitudes, and I ...
2
votes
2answers
33 views

Reducing number of attributes (classes of marks) in marked point pattern using R?

I have a marked spatial point pattern loaded in R (from Shapefile). This dataset has 4 different "classes" of marks: class A, class B, class C and class D. Each class of marks is categorical. One can ...
3
votes
1answer
3k views

Interpreting Moran's I results

I have a sample of data taken from regions of a country! i want to test if there are any spatial autocorrelation on my data using Moran Indice test. the null hypothesis: is that there is no spatial ...
1
vote
0answers
17 views

R: how to export my Kriging result to raster? [duplicate]

I want to export results of my Kriging calculation to raster or some other format. My aim is to take this converted file and intersect it with another SpatialPointDataFrame (location of the point ...
4
votes
1answer
223 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...
0
votes
1answer
218 views

How to download the MODIS data from R using MODISTools library

I want to download MODIS product MOD13Q1 to use the Scientific Data Sets (SDS) layers of "250m_16_days_EVI", "250m_16_days_pixel_reliability". MODISTools library in R software provides the capability ...
1
vote
1answer
28 views

Splitting polygon into multiple fixed-width polygons along edge in R?

I have a polygon and would like to spit into grid with fixed width along with one edge. The result will look like the figure below. Are there algorithm to achieve my target? My language is R, but ...
2
votes
1answer
27 views

Convert/display grayscale image as RGB

i want to display raster images using a color palette in a leaflet map. The images are tiled (gdal2tiles) from a single .tif to multiple .png files with values from 0 to 100. I add the image using ...
2
votes
1answer
16 views

Dealing with spatio-temporal data: sp::over function in R

I have a dataset called carcass, which contains occurrences of livestock carcasses in Norway in the last 5 years (each raw corresponds to a single carcass, which has an ID number, along with ...
1
vote
1answer
38 views

How to apply functions iteratively to SpatialLines in R?

Looking for a way to apply the 'rgeos::gIntersection' and 'rgeos::gDifference' function iteratively to SpatialLines, I created 3 SpatialLines as an example: #from the sp vignette: l1 = ...
1
vote
0answers
21 views

How to download GPM rainfall data automatically?

How to download GPM rainfall data automatically? I tried to use DownloadGPM.R Function, but I think I had problem with Axel download accelerator. I guess I need to declare Axel as environment ...
4
votes
2answers
2k views

How to convert a spatialpolygon to a SpatialPolygonsDataFrame and add a column to the atrribute table

coast<-readShapeSpatial("coastline.shp") landc<-readShapeSpatial("landcover.shp) ro<-readShapeSpatial("roads.shp") bc<-gBuffer(ro,width=100) landc$ratings=1 ...
0
votes
1answer
42 views

smoothing raster layer stack using R

I'm trying to do smoothing for .tif stack layer of 331 bands. I used this code in R : setwd("D:\\") raster_dir <- "D:\\Modis VI quality\\" VI <- ...