Statistical computing language and software environment.

learn more… | top users | synonyms

2
votes
0answers
100 views

does the census bureau provide clipped / water-bounded files as nice as certain cities?

this question says they do but when i compare what i can find in the cartographic boundary files section of the census bureau's gis page, they don't look as good as certain shapefiles provided by the ...
1
vote
1answer
64 views

R: Change RasterBrick dimensions from TXY to XYT

I am having difficulties in R getting climate data from a netCDF file into the right format for further analysis. I am working with gridded global data on Standardized Precipitation Index (downloaded ...
1
vote
0answers
16 views

Is there a replacement for ggmap::geocode? [duplicate]

I've found ggmap's geocode() function really useful over the years as it's got a 'fuzzy' search implemented by Google. However, I'm wondering if there are alternatives, because geocode() seems ...
1
vote
1answer
208 views

Example of simple kriging with varying local means in R

I'm looking for an example of how to use simple kriging with varying local means in R. I've tried google, a few books (for example Aplied Spatial Data Analysis with R) and I've found almost nothing. ...
-2
votes
1answer
66 views

Points don't show up on an overlaying shapefile?

library(ggmap) library(rgdal) cities<-c("city A","city B") country<-readOGR(dsn=".",layer="country") gc<-geocode(cities,source="google") gc$names<-c("city A","city B") coordinates(gc)<-...
4
votes
2answers
379 views

Alternatives of spTransform?

Let's say that we have a shapefile with a certain projection. s<-readOGR(dsn=".",layer="Spain") We also have the airports of Spain as points in another projection. a<-readOGR(dsn=".",layer="...
8
votes
3answers
661 views

Fastest way to convert big raster to polyline using R or Python?

I have a big raster file (129600 by 64800 pixel) with global water bodies (1 bit values 0 and 1) and try to extract ocean and inland water shorelines. I've tried with ArcGIS and QGIS to convert from ...
9
votes
2answers
259 views

Cleaning blocks pattern from urban LIDAR (blocks uplifted from streets)

We have a 1-meter LIDAR DEM from a city. A small subset can be downloaded from this link: This screenshot shows raw DEM with gray palette (darker belts are streets, and greyish and whitish ...
3
votes
3answers
1k 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 ...
1
vote
1answer
136 views

Pan sharpening using R

I want to pan sharpen my Landsat8 scene (according to FIHS algorithm, see here). According to the paper the following equation should do the pansharpening: MS_high <- MS_low + (pan - I) "...
0
votes
1answer
57 views

Split by attributes with 2 fields

I am using SplitLayerByAttributes toolbox. It’s a very useful tool for me to process certain activities. As you know, I can split the feature class as per 1 field. What if I want to split the features ...
1
vote
1answer
2k views

How to save R raster plot for ArcGIS raster file

I want to save my raster plot that produced with R to ArcGIS raster (.img, .tif and ASCII) file.
4
votes
2answers
1k views

R function for Thiessen Polygons

For a point data set with temperature data I have, I want to create their Thiessen Polygons in R. So far, the only algorithm I have found is that of Carson Farmer (http://carsonfarmer.com/2009/09/...
2
votes
1answer
231 views

R-bridge minimum requirement

The new R-bridge, unveiled at uc2015, has a dependency on ArcGIS 10.3.1. Has anybody tested it with older version of ArcGIS Desktop? I need to know why it depends on 10.3.1. I need to run it on 10.2....
1
vote
1answer
47 views

Adjusting sampling points to closest land surface QGIS

I have many sampling points in mangroves of Southeast Asia and am trying to extract climatic data for each using the WorldClim dataset. However, many of the points fall on areas classified as "water" ...
1
vote
0answers
287 views

Converting UNIX time to datetime in R [closed]

I am trying to convert a 13 character UNIX timestamp to datetime and I am getting a year that is in the 44000's. I'm not sure what I am doing wrong, hopefully someone can help. library(maptools) GMT ...
0
votes
1answer
99 views

variogramST in R

I've been following the tutorial in this site to compute the spatio-temporal variogram in R. The result from variogramST is the object called var, and it contains several NAs for np, dist, and gamma, ...
2
votes
1answer
67 views

How to separate out bands from hyperspectral data using R

I'm having one hyperspctral data with 426 bands, with dim 250*250*426. I want to seperate each bands. I tried like this img <- stack("name.tif") rr <- array(, dim= c(250,250,426)) for(i in 1:...
0
votes
2answers
377 views

how to remove cloud values of Landsat image from R [closed]

How can I remove cloud value of the Landsat8? How can I estimate cloud values for eg: modis cloud values are from 700 to 800. By using R how can I estimate>
6
votes
2answers
227 views

Create large amount of random points in binary raster

I want to create a point vector dataset of 10000 points (or larger) within a binary raster, where the points should be constrained to areas where the raster value is 1. I tried the following steps. ...
7
votes
1answer
1k views

Batch geocoding in R

I need to find a package in R or source code in R that performs batch geocoding using Bing, Yahoo, OpenStreetmap. Do you have any suggestions?
6
votes
1answer
111 views

How to make a wedge shapefile in R?

I'm quite okay at R programming, but fairly new in GIS (especially in R). I have a set of coordinates and I'd like to create a custom shape around them. What I want to do step-by-step: Create a ...
1
vote
0answers
26 views

Batch Geocoding in R over daily limit with pay as you go? [duplicate]

I have found the the geocode function in ggmap to be very helpful in the past. Now I am trying to geocode about 60,000 addresses. The google maps API daily limit is 2500. Google offers a "pay as you ...
-2
votes
1answer
500 views

How can open a shapefile using R maptools package? [closed]

library(maptools) mn.map.shp = readShapeSpatial("minnesota.shp") Results this error at second line: Error in getinfo.shape(fn) : Error opening SHP file
18
votes
3answers
4k views

Using R with ArcGIS Desktop

In response to comments below, a re-phrasing of the question. I've used the integration of R with QGIS, and found it a very uesful tool. Having taken a look at the R toolbox extension for 9 and the ...
5
votes
0answers
144 views

Interpreting semivariogram high nugget

I made a semivariogram in R using package gstat, variogram() function. I want to see if there is spatial autocorrelation in the residuals of my model (species abundance as a function of habitat, ...
6
votes
1answer
392 views

Performing cluster analysis in R with Point Datasets

I am a R-Beginner, but I want to use it for some spatial statistics. I want to check the relation between the occurrence of ill persons and the existence of water in a area of interest. To give you a ...
0
votes
1answer
44 views

Exporting “pennLC” object stored in R as a shapefile

I need to export the pennLC object from package:SpatialEpi as a shapefile. pennLC contains different components of lung cancer dataset for Pennsylvania state such as: $geo, a data frame $data, a ...
-1
votes
1answer
149 views

How to read metadata file from R [closed]

How to read metadata file parameters( like sun elevation,REFLECTANCE_MULT_BAND etc ) from the user interface or from R How can end user validate the parameter of metadat file(.mtl) How can I pass ...
3
votes
2answers
109 views

How to remove bounding box extents based on reflectance values of the image of Landsat8 by using R programming

I have taken the reflectance image from Landsat8 bands. Now the extents of bounding box showing negative values , to proceed to my next level of the process i need to remove this negative values and ...
2
votes
0answers
51 views

About shp for ggmap in the R

I am a new in R and GIS, so I kindly ask you on advice. I had a shp which I want to plot on google map. The shpfile is here: http://www.filedropper.com/exercisestanddata First I do this: ogrInfo(...
6
votes
1answer
4k views

How to intersect lines and polygons in R?

Following on from an answer about intersecting polygons with lines to chop up a polygon into smaller polygon units (in QGIS), I wanted to try the same thing in R. However, I cannot seem to find a ...
2
votes
1answer
1k views

Using Regression Trees in ArcGIS

I have a follow up question to Statistical analyses for survey against raster? I have a very similar problem as the initial poster David. I am creating a multivariate model from raster data in order ...
-2
votes
1answer
487 views

How to calculate statistical NDVI in R (for Landsat8) [duplicate]

Is there any function to calculate statistical NDVI in landsat8 package
-1
votes
2answers
592 views

unable to find an inherited method for function ‘writeRaster’ for signature ‘“SpatialGridDataFrame”, “character”’

Error in (function (classes, fdef, mtable) library(rgdal) x<-readGDAL("LC81360442014322LGN00_B5.TIF") z<-reflconv(x,2.0000E-05,-0.100000) writeRaster(z,filename="test.tif",format="GTiff",...
0
votes
0answers
93 views

Variogram to assess spatial auto-correlation in grouped surveys in R

I have a large dataset with coordinates of survey locations, and results for presence or absence of a species. These surveys are spatially clustered, with four in a line separated by 1km each, and the ...
1
vote
0answers
19 views

How can I quantify whether locations of one 'type' are closer to one another compared to a second 'type'?

This might have an obvious answer but I'm new to GIS and can't seem to find the right technique... I've got GPS points for a bunch of households and have calculated the distance (in km) between each ...
6
votes
1answer
153 views

Segmentation of a big raster

I'm actually trying to convert one raster into a SpatialPolygons. I've got a big raster (12 millions pixels) with a lot of zeros, and some negative values (representing some depression in a DEM). I ...
0
votes
1answer
45 views

Debugging spplot() and reflconv() on R Cannot allocate vector of size …?

I am trying to perform LANSAT data of Digital Numbers (DN) to reflectance in R(coding) using reflconv(). When I apply that function to example it shows below error: >library(rgdal) x <- ...
3
votes
1answer
134 views

Define cellsize from spsample in meters

I use a shape file to create a regular grid using spsample, from sp package. Two problems here: (1) I would like to define the size of my grid cells in meters, but I've checked the documentation of ...
1
vote
1answer
938 views

Saving output 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 (...
0
votes
0answers
91 views

Union overlapping shapefiles (poylgons) and merge the attributes where overlap occurs in R

I have 30 shapefiles of overlapping polygons and would like to union them into one shapefile. The attribute table of the 30 layers is exactly the same (two columns, Type and Origin). A new polygon ...
1
vote
1answer
318 views

Using Loop for calculating NDVI

I am using MOD09Q1 now and I would like to calculating NDVI whole year data. One year, I have 46 band 1 images and 46 band 2 images. list.files I creating list of band 1 and list of band 2. ...
1
vote
0answers
59 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 ...
2
votes
2answers
2k views

Opening .lyr file via rgdal/OGR

According to this answer I can read in a .lyr file via OGR as long as it was made with ArcGIS version 10. I have such a file (from the ArcGIS Maps & Data page) but I can't quite make R's readOGR ...
0
votes
1answer
100 views

How to extracting values of 40 plots from different bands of raster image in R?

I have 40 plots(samples) and I want to extract corresponding values from different bands of RapidEye image. So at the end I want to have totally 200 values of 5 bands in the format of txt. Till now I ...
3
votes
1answer
116 views

Find which latitude and longitude values are falling in which raster cells

I am an absolute beginner in analyzing geographic data with R. I have a database with longitude and latitude of around a million buildings in a country. I would like to find the average density of ...
1
vote
0answers
74 views

R raster stack 4 dimensions

I am trying to create a raster stack / brick from a list of ncdf files with the dimensions 160x160x24x1). r = stack(files) This gives me a raster stack but without the 4th dimensions which is ...
2
votes
2answers
84 views

How can I use a DBF layer in an R script in QGIS?

How can I use a DBF layer in an R script in QGIS? For example, I can input a vector by using ##Layer=vector, and a CSV table with ##Table=table, but it will fail if the table is based on a DBF file. ...
3
votes
1answer
3k views

merging 2 polygon adjacent to each other using R

I have several shp files which are the components (or tiles) of a map, I have joint all these tiles together to have a large shp file. Then I realized some of the polygon is now split into 2 halves, ...