Statistical computing language and software environment.

learn more… | top users | synonyms

0
votes
1answer
21 views

Raster Mosaic fault using R

I'm trying to mosaic two raster images using R. I used this script: setwd("D:\\Modis EVI\\") raster_dir <- "D:\\Modis EVI\\" r1 <- brick(paste(raster_dir,"evi8_stack.tif", sep="")) r2 <- ...
1
vote
1answer
23 views

Remove raster pixel in R

This is the 32 year trend NDVI (raster data) I would like to remove the green grids data by using R (the green pixel value is larger than 0.03) and this one is R code that i use but i cannot ...
0
votes
0answers
20 views

How to define the Northern Lambert Conformal Conic CRS?

The dataset I got my hands on uses the the Northern Lambert Conformal Conic CRS? I was not able to find a ESPG number using http://prj2epsg.org/search so I decided to define a custom CRS in R but I ...
4
votes
2answers
612 views

LSA-SAF data reprojection with GDAL and R

The Land Surface Analysis Satellite Applications Facility (LSA SAF) distributes several products derived from the EUMETSAT satellites data. I am working with R to read and process the data (the file ...
-1
votes
0answers
23 views

How to convert .HDF to .Tiff and keep the number of the output files exactly as the input?

I have 332 .HDF files that I have processed in R to convert them to .Tiff. It has worked fine, however, the number of output files are not the same as the number of the input files! Input is 332, ...
1
vote
1answer
22 views

Error in rasterToPolygon trying to get several coverage layers

I have a very strange problem: I am trying to get a raster with coverage from 0 to 1 in each cell for types of vegetation, I started with a raster of vegetation type, but I needed a raster with 12 ...
1
vote
0answers
26 views

R script for find the minimum distance [closed]

I need an R code to select which polygon in a shapefile has got the minimum distance from a given point.
0
votes
1answer
32 views

how to create raster stack using R

I have 300 .tiff files, that i need to create raster stack of them. Can you please help me with R code for that? I can do in Arcgis, Envi, R. I have tried: setwd("F:\\MODIS\\Modis EVI\\HDF8 EVI") ...
1
vote
1answer
218 views

way in R or Python to add image/video in map plot marker click popup/infowindow

I have many different (lat,long data) points of world associated with unique place names and corresponding each point specific image or video data. Now I want to create an html file where if user ...
2
votes
2answers
730 views

Is there an R equivalent of NEAR in ESRI ArcGIS?

I have 54 points in one data set A I have 300+ data points in a dataset B. I wish to find the distance for each point in dataset-set B, to points in dataset A. Nearest Neighbor packages appear to ...
9
votes
2answers
560 views

Help choosing most suitable bivariate choropleth mapping coloring scheme

I am trying to overlay a nominal variable (2 categories - rainfed vs irrigated) in a choropleth map over a numeric variable. In order to highlight the difference in pattern among the two ...
0
votes
0answers
19 views

How can I create a NetCDF from precipitation rasters data in R? [duplicate]

I have a lot of precipitation raster data (more 5000 files GeoTIF), created by kriging interpolation. All grid cells have the same size and the maps the same spacial extent I want merge all data in ...
2
votes
1answer
434 views

How to extract percentage of land cover from a raster according to multiple shapefiles?

I have a raster map that represents land cover of a region, and different shapefiles (total 120) in the same region. I want to know the percentage of every land cover for each shapefile. For example, ...
1
vote
1answer
20 views

How to work with Baccini et al.'s biomass map from ArcGIS?

Baccini et al. (2012) has produced one of the best biomass maps available for the world. One can find the map from one of ArcGIS repositories here. I downloaded the .pkinfo file for that map but I ...
1
vote
0answers
14 views

R Freezes when Plotting Spatial Data Frames

I have run into a problem where anytime I try and plot a Spatial DataFrame in R, the console freezes, regardless of the size of the file. I imagined this must be due to an issue with the path to one ...
0
votes
0answers
50 views

Call gdal_proximity() from R

I have a raster including eight land cover types (pixel values ranged from 1 to 8): > r class : RasterLayer dimensions : 9887, 30532, 301869884 (nrow, ncol, ncell) resolution : 10, 10 ...
1
vote
0answers
18 views

Group location points that have a cumulative value 2000

I have a set of 400 locations (points) joined to a table having population attributes. I want to form groups of nearby locations such that the cumulative population is 2000, group can have maximum 15 ...
1
vote
1answer
66 views

How to digitize 3D geospatial data using FOSS applications

I need to digitize soil horizon data from photos to be later used in a 3-dimensional GIS workflow. For each photo I have DGPS coordinates of the top left and right corner and the depth of the ...
2
votes
0answers
17 views

TryCatch for if two extents?

I have some spatial processing in R that I want to do. After selecting out the appropriate polygon (clip1) I want to crop the raster on that polygon with the command clip1 <- crop(imported_raster, ...
0
votes
0answers
13 views

locate sample geo trace data

I'm working on a project where I need to form clusters of GPS coordinates that are separate from noisy geo-data points. I am using DBSCAN/OPTICS XI algorithms in R on randomly generated test data, but ...
2
votes
1answer
54 views

Recommendation for getting TIF raster from modis HDF, for processing in R

I've downloaded a set of modis mod13q1 HDF files, e.g. MOD13Q1.A2016097.h11v05.006.2016114040325.hdf. I'd like to load the 250m_16_days_EVI and 250m_16_days_pixel_reliability bands into R as rasters. ...
1
vote
0answers
29 views

Polygons created with R are corrupted when imported into ArcMAP

I created some polygons using R and to get a better visualization and some other stuff I need to import those into ArcMap. However, when I imported the polygons I realized there are some holes in some ...
1
vote
1answer
22 views

R: Krovak projection in R - identical in ArcGIS, different in R?

I have two files: shapefile and raster, both in the same coordinate system in ArcGIS: EPSG 5514. Check information on projection in ArcGIS for raster: S-JTSK_Krovak_East_North WKID: 5514 Authority: ...
2
votes
1answer
64 views

Error in In predict.gstat in R

I am trying to do kriging by gstat pcakge through R. Here is my Code: krg<-krige(var ~ 1,data,grid,model=vgm,nmax=30) But each time I recive this warning: 1: In predict.gstat(g, newdata = ...
0
votes
3answers
540 views

Trilateration algorithm for n amount of points in R?

I read the page Trilateration algorithm for n amount of points and I have to do similar kind of trilateration in R. I am not that expert in R, can anyone help, if there is similar function/package ...
1
vote
0answers
30 views

Spatial AND temporal join between two point data sets

I was hoping I could get some pointers in how to approach this issue, which is tricky version of a spatial join. I have been searching and I can't find a particular solution to what I need. I have ...
0
votes
1answer
177 views

Plotting an eia shapefile in google maps with R?

I'm trying to plot a eia shapefile onto google maps. I have a code working with the map funciton but can't get it to work with the ggmap function. Where am I going wrong? Working code: ...
1
vote
1answer
25 views

Proximity filter

The aim of what I'm currently messing around: Creating a distance buffer around points in a SpatialPointsDataFrame If the value a this point is > 0, always keep it If the value a this point is = 0 ...
5
votes
1answer
94 views

Krovak projection EPSG 5514: Correct data transfer between R and ArcGIS?

I have a problem which many ArcGIS, QGIS and R users from Czech Republic and Slovakia may encounter. In ArcGIS, I am working EPSG: 5514, Projected Coordinate System, defined as: ...
0
votes
1answer
215 views

using raster.predict for Cubist model, factor levels list

I would like to thank in advance any help or advice. Let me know if I need to provide more reproducible data/code. I am working on a Cubist prediction model for soil texture, using continuous and ...
1
vote
0answers
37 views

Creating an ArcGIS toolbox with an R script

I posted this over on another website but since I think it might be an issue with ArcGIS, I decided to post it here also. So I have the following script in R. in_params[[1]] represents the input CSV ...
0
votes
1answer
22 views

How to extract single band from multi-band radiometry raster using R?

I am having trouble extracting single band from multiple multi-band raster. I have no problem doing it on single raster file. Is there a way to loop through all the raster images to extract a single ...
2
votes
0answers
34 views

Failed to create a raster-brick from a netCDF file (.nc)

I am trying to build a raster-brick from a netCDF file (network Common Data Form) leave this link for more information () but to make a brick (which is part of the package raster) throws me the ...
2
votes
1answer
29 views

What is a UK boundary shapefile available for download via persistent URL (for cURL) with parish-level detail?

I'm in the process of generating a series of choropleth maps for a journal article and am writing this up in R as an .Rmd. My aim is to set this up as reproducible research, that is in a form which ...
0
votes
2answers
53 views

R: creating an error matrix from two rasters

I have two land use rasters (6 classes each) and I'd like to create an error matrix (errors of omission and commission) from them both. The premise is that I have a very fine 'base' raster and have ...
0
votes
2answers
22 views

mapping with spplot

This question should be easy! I am using spplot for visualizing a map with suitable results. However, I would like to know how to set the limits of z, because the values are not being discriminated. ...
0
votes
1answer
39 views

Use of Leaflet plugins in 'Leaflet for R' to connect with ArcGIS online

I was wondering if it is possible to use the plugins for the normal (non-R) version of Leaflet, and if, how I can do that. reason I need to do that is, that some of the data I need to show are hosted ...
2
votes
0answers
19 views

Creating a reference raster in R

I need to create global reference raster (100km grid, Eckert IV). Later, I will project other various rasters (with different resolutions and projections) into the resolution/projection of the ...
1
vote
2answers
190 views

projection differences in R

I am trying to reproject some climate data that I downloaded from http://gdo-dcp.ucllnl.org/downscaled_cmip_projections/dcpInterface.html#Welcomein The BCSD-CMIP5 data is found under the tab ...
8
votes
2answers
4k views

Clip Spatial object to bounding box in R

Given a Spatial object in R, how do I clip all of its elements to lie within a bounding box? There are two things I'd like to do (ideally I'd know how to do both, but either is an acceptable solution ...
2
votes
2answers
460 views

Clump and eliminate classes in a raster in ArcMap or R

after interpolation of sediment classes, I obtained a raster which looks like the image below (sample of the river). I would like to clump and eliminate areas within the raster that are below a ...
0
votes
0answers
28 views

Why ArcGIS crashes when I load raster file created in R? (HFA format)

I was processing my raster operation in R, as I think it's easier then in ArcGIS: - read shp + raster1 - crop them - classify raster1 by location - create classified raster2 - export the ...
0
votes
1answer
29 views

rasterVis plot3D in R: obtain consistent x-y axes with extent() / projection() shown in RasterLayer object

I have a RasterLayer object z: class : RasterLayer dimensions : 100, 100, 10000 (nrow, ncol, ncell) resolution : 0.06534037, 0.03165736 (x, y) extent : -96.63908, -90.10505, ...
0
votes
1answer
22 views

Finding neighbors for non adjacent polygons in R [closed]

I'm looking for a method or function in r-project that create a neighborhood matrix from a group of polygons that don't touch each other, like polygons of set of blocks separated by streets in a town. ...
0
votes
1answer
47 views

Creating SpatialPolygons in R

I am trying to do kriging with my data. But there is a problem here. I want to use krige command in gstat package. But I do not have "Spatial object with prediction/simulation locations". How I can ...
2
votes
1answer
8 views

Retain “region” elements for when using fortify() on a SpatialLinesDataFrame

I'm layering multiple plots using ggplot2 to create a map of Alaska that combines school location and the major road network. I've learned that despite being labeled "major roads" there is wide ...
11
votes
2answers
4k views

Aggregating points to grid using R

I have a question with regard to spatial aggregation in R. What I am trying to do is aggregate a point dataset to a grid. I am unsure however how to do this as I have little experience with this sort ...
0
votes
3answers
67 views

Converting position data created in R to GPX format?

Does anyone know if there is a way to convert position data (latitude and longitude) created in R to GPX format to use with navigation software? I have been searching and havent found an answer. ...
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 ...
3
votes
0answers
197 views

Using R to map matching GPS data? [closed]

I have many vehicles GPS data points, which consist of "timestamp, longitude, latitude, fuel consumption and distance" and the data is OSM(open street map) based. I want to match the data points to ...