Questions tagged [rgdal]

R Bindings for the Geospatial Data Abstraction Library

Filter by
Sorted by
Tagged with
1
vote
2answers
55 views

How to convert MTM coordinates to WGS84 with R?

I have MTM coordinates in the following format df <- data.frame(xcentro = c(250605,250770,250871),ycentro = c(5186087,5185610,5186140)) and need them to be converted in the WGS84 format. I expect ...
0
votes
0answers
13 views

Are these sp files equal in terms of proj4string? Can I do some transformation of the proj4string?

summary(SCat_sp_E6) Object of class SpatialPolygonsDataFrame Coordinates: min max x -74.095072 -74.018988 y 4.647323 4.832934 Is projected: FALSE proj4string : [+proj=longlat +...
6
votes
3answers
130 views

sp::proj4string equivalent in the PROJ6 framework?

I am confused about the shift from PROJ4 to PROJ6 in the (relatively) recent versions of the spatial packages of R (and gdal for that matter). I used to use sp::proj4string() to get the projection ...
0
votes
0answers
36 views

How do I specify the CRS in light of the rgdal update?

My data (df) are in Africa Albers equal area projection: x_ y_ t_ id 704109. -3029239. 2015-05-12 09:28:00 A 703680. -3029412. 2015-05-12 10:00:00 A 703043. -3030028. ...
0
votes
1answer
27 views

Reprojecting GOES-16 data in R

I'm trying to plot GOES-16 data with world boundaries. Data comes from GOES-16 files such as this. Data comes with "GOES-R ABI fixed grid projection" reference system, which to my better ...
0
votes
0answers
24 views

Can not load spatial data from MSSQL on Linux

I have a table on the MS SQL server. I am trying to import the data in R using the rgdal package. I can import the table on Windows without any problems using the following code: readOGR(dsn = "...
1
vote
0answers
21 views

Installing rgdal - cannot solve dependencies

I am trying to build a docker image based on osgeo/gdal:ubuntu-XXXXX. However it keeps failing when trying to install rgdal. I've tried this on gdal-3.1.1/proj-7.1.0 and gdal-3.0.3/proj-6.3.0, as well ...
1
vote
0answers
54 views

Using rdgal on R- error- No GDAL installation found

I trying to map burnt area using hdf data. For this i have all the neccessary packages, including the rdgal pacakge. However when i type in my burnt area code, i get the response I have downloaded ...
1
vote
1answer
150 views

Error “projects.h file not found” while installing rgdal in R on Mac

I cannot install rgdal on MacBook with Catalina. I installed r using Brew and r is up to date. I have the gdal framework installed. I've tried to install rgdal from source, but I keep getting an ...
2
votes
1answer
87 views

r - Setting custom projection in PROJ >=5 paradigm

Before PROJ >=5.0, one could set a custom projection centered anywhere in the world by adjusting the +lon_0=and +lat_0= arguments in proj4string specification. For example, a projection centered on ...
1
vote
0answers
117 views

R can't work with ecw while gdal/qgis does. Ubuntu 20.04 config

I'm an Ubuntu 20.04 user with all configured to use .ecw raster format. I can read these files in QGIS (3.10) and work with them without a problem. This is what I get in console: $ gdalinfo --...
1
vote
0answers
76 views

Updating PROJ GEOS and GDAL in R

I have been having trouble linking to GEOS 3.8.1, GDAL 3.0.4, PROJ 6.3.0. I have tried reinstalling sf, rgdal, raster, and rgeos packages. I still couldn't link to GEOS 3.8.1, GDAL 3.0.4, PROJ 6.3.0. ...
0
votes
0answers
81 views

GDAL conflicts between GDAL and PostGIS versions in R

I'm trying to learn postgis and I'm using R as an interface to make queries under Linux. I'm using GEOS 3.8.0, GDAL 3.0.4 and PROJ 7.0.0. When everything is up-to-date my R ecosytem (sp, sf, rgdal) is ...
2
votes
1answer
321 views

Install latest GDAL on docker image of CentOS

I am trying to install gdal latest version on a docker image of CentOS. Below is the minimal code. # Start and connect . docker run --rm -it centos:latest bash yum -y upgrade yum -y install wget yum ...
0
votes
0answers
28 views

Get Nearest Grid Data Using Leaflet and R

I am trying to write a piece of code in R using the Leaflet, sf and rgdal libraries which will identify the closest grid cell to the user's click and retrieve the data associated with that grid cell. ...
0
votes
1answer
67 views

Erreur gdal_polygonize.py package ARfun with R

I'm a french student, in advance I'm sorry for my english. I'm installed the package APfun on R to transform Raster to Polygons with Appolygonize function. But I have a error : crownsPoly <- ...
0
votes
2answers
100 views

Using writeOGR on a list of SpatialPolygons

This question has been asked before, but I find the answer confusing since there are no reproducible examples that link the question and the answer, and the answer doesn't actually use the writeOGR ...
2
votes
1answer
170 views

R spatial: st_overlaps to subset also fully covered polygons?

I want to subset all polygons in fc1 completely or partially overlapping with fc2. But, the st_overlaps subsets only polygons that overlaps with the border of fc2, but skips the fully covered polygons....
0
votes
1answer
108 views

Mass conversion of .JP2 images to .TIF format in R

I need to convert a list of .jp2 files to .tif. I cannot do it manually since it is a large number of files. So I can't use: gdal_translate(path_vhr_imag_1, "splt1.tif") gdal_translate(...
7
votes
1answer
221 views

R spatial: Erase one polygon from another: correct use of `st_difference`?

I am trying to erase one shapefile from the other, similar to ArcGIS Erase tool: This worked pretty well with rgeos::gDifference(spdf1, spdf2), which needs as input geometry of ...
1
vote
1answer
45 views

Use geometry binary predicate list to subset simple feature object?

I am trying to create a new simple feature from two overlapping features. I want to subset all polygons that are fully or partially located within my buffer (red border). However, the st_overlaps ...
1
vote
1answer
26 views

Keep getting NAs when I run sp's over function on LatLon points on shapefile polygons

I am not sure why i am keep getting NA whenever I run the Over function with Latitude and Longitude point on the polygon from shapefile. Please note that this is first time for me doing the spatial ...
2
votes
1answer
38 views

Check if points overlap with polygon

I have data that maps lat lon of a storm for every day as follows (I have removed the day ID so it has only year and month) lat_lonDat <- structure(list(Latitude = c(9.1, 9.15755, 9.2, 9.21496, ...
7
votes
3answers
234 views

Find maximum extent out of list of shapefiles (in projected coordinate system) in R

I have a script where I crop a list of rasters by the extent of one of many shapefiles. So far all of by shapefiles have been such that I could manually select the one with the largest extent (i.e. ...
2
votes
2answers
125 views

Loading shapefiles for reprojection using sf

I'm trying to automate my workflow, but having trouble with one step. My data come as processed shapefiles saved in .RDA files; I load them to a PostGIS system, cleaning up little issues (spTransform()...
1
vote
0answers
118 views

Attributes encoding problem with rgdal::writeOGR

I performed an intersection of two SpatialPolygonsDataFrame on R (3.5.1) using raster::intersect. Here's the structure of the resulting SPDF attributes table > str(shp@data) 'data.frame': 944 ...
1
vote
1answer
345 views

Keep spatial lines intersecting spatial polygons using R?

I am trying to keep spatial Lines objects intersecting a spatial polygon in R. The intersect function from raster package is doing actually a crop. But what I need is to keep the lines intersecting a ...
0
votes
1answer
2k views

GDAL 2.4.2 from UbuntuGIS unstable broken because of missing libogdi.so.3.2? [closed]

After an unthoughtful apt-get update today, GDAL no longer works on my Ubuntu Machine (bionic) because libogdi.so.3.2 is missing: daniel@gin-nuest:~$ gdalinfo gdalinfo: error while loading shared ...
1
vote
1answer
46 views

Read Postgis geom available as string variable with R

I have a Postgis Database from which I query data. I know that I can use the rgdal::readOGR() function with the PG driver to get geoms from the database into R. However, I am also querying lots of ...
0
votes
1answer
432 views

How to get WMTS tile coordinates?

What can I do with the given information: I don't know where I am in the World having some particular tile. It's only an ocean (or sea). "https://services.arcgisonline.com/ArcGIS/rest/services/...
0
votes
1answer
550 views

Identify polygon grid in which the point belongs to?

I have a dataframe with 3 columns , which has my custom lat / long with a failure rate associated with each Lat / long head(df) Latitude Longitude Failure.rate 51.3 ...
1
vote
1answer
194 views

Trying to make calculation with gdal_calc.py?

I have a tif image that I'm trying to run gdal_calc.py. In the source file the values are between 1,118-9,704 and I want to divide all the values by 10,000 to get results between 0-1. This is the ...
0
votes
0answers
302 views

Convert DXF file to RasterLayer

I am an R user and have a .dwg file of city buildings that I wish to transform into a RasterLayer using the heights as elevation data. The metadata for the DXF file states: DXF File Metadata: ...
0
votes
0answers
49 views

GeoTools and PostGIS coordinate transformation gets different result [duplicate]

I've got a problem transforming coordinates. The result differs in PostGIS and GeoTools. I want to transform from EPSG:4326 to EPSG:3398. PostGIS: indicia=> select ST_AsText(ST_Transform(...
3
votes
1answer
256 views

gIntersection error: linearRing not created

I am trying to do something relatively simple but with large shapefiles. I have a shapefile of protected areas in the USA: and watershed boundaries: I want to calculate the percentage and ...
0
votes
1answer
40 views

Unable to plot in R some locations over a shapefile

I have a shapefile of polygons that I got from here, and a CSV file with longitude, latitude data (i.e. lon= -71.4; lat=42.5). The shapefile has the following CRS information: "+proj=lcc +lat_1=41....
1
vote
1answer
203 views

How to access GDAL utilities from command line using Rgdal installation? [closed]

I have just installed rgdal from within R install.packages("rgdal"). I can now successfully do library(rgdal) with the following message: Loading required package: sp rgdal: version: 1.3-9, (SVN ...
3
votes
1answer
2k views

Reading shapefile in R with shapefile function?

I want to read a shapefile(polyline) into R, using rgdal and raster packages. The code looks like: d<-file("/Volumes/Environment_layers/original/WaterRoad/Road.shp") o<-shapefile(d) but I got ...
2
votes
0answers
124 views

Warning message in raster::rotate(): “this does not look like an appropriate object for this function”

I'm attempting to rotate a stacked raster using the raster package in R. The raw stacked raster looks like this: When I attempt to rotate the raster to change the longitudinal scale from 0-360 to -...
0
votes
0answers
2k views

Get Google Maps driving directions between points in a KML file?

I'm trying to improve a workflow that currently involves manually composing a multi-destination driving route in Google Maps. I have the ability to export a KML file containing the points we want via ...
1
vote
0answers
39 views

Resolving topology errors after merging polygons

I have two SpatialPolygonDataFrames, and I have two merge them. SPDF1 contains watershed are falling within the country, SPDF2 has watershed are falling outside the country, so the both should match ...
0
votes
0answers
528 views

Extract raster cell values by polygon field

I have a raster as a RasterLayer (raster) and a polygon layer (bound) with 38 individual polygons as a SpatialPolygonsDataFrame. I used the polygon to mask from a larger raster. Each polygon contains ...
1
vote
0answers
528 views

Possible to drop NA geomType in readOGR?

I have a gdb file which is giving me trouble. f = '/path/to/folder.gdb' layers = ogrListLayers(f) geom_list = lapply(layers, readOGR, dsn = f) Gives an error on one of the layers: Source: "f&...
0
votes
1answer
510 views

batch merge error: Error in match.arg(method) : 'arg' must be of length 1

I am attempting to read in a series of polygon shapefiles from a directory and batch merge them into a single shapefile. They all have the exact same field names. I am using this code: all_data <-...
2
votes
1answer
208 views

Binding multiple polygons in R

I have a directory that contains around 700 .shp files. I want to create a map in R that displays all of them. The easiest way i could think to do this was to bind them all together and then use ...
0
votes
2answers
780 views

Extracting non-spatial data from imported KML file using R?

I have a KML file which was created using Google's My Maps. The original file can be downloaded here: Google My Maps Using R, I can import this using the "readOGR" function of the rgdal library ...
1
vote
1answer
221 views

Parallel rasterizing with GetCover option? [closed]

I'm trying to rasterize a shapefile in parallel using the solution from this post: Processing vector to raster faster with R However, this fails to give the correct answer when the GetCover = T in ...
0
votes
1answer
41 views

manipulate ASCII grid

I have an ASCII grid (Data information in this box link https://app.box.com/s/jpdalfjnyfhqbvrghd8csu3pz5b2fl2j) with the following projection information Projection POLAR Units METERS ...
2
votes
0answers
99 views

Points and bioclim layers plot correctly in R, but not in Maxent.jar

I am trying to use Maxent.jar. This program uses a CSV file of points and ascii bioclim layers (basically, a map of the Americas with climate data) and overlays them on top of one another. I find ...
0
votes
1answer
270 views

Pattern matching with grepl on SpatialPolygonsDataFrame object in R

I am trying to use grepl as part of an if control statement in R. I want to subset a SpatialPolygonsDataFrame object in my list (i.e. googleProjectedGroundTruth or googleProjectedFallow) only if its ...