Questions tagged [rgdal]

R Bindings for the Geospatial Data Abstraction Library

Filter by
Sorted by
Tagged with
0
votes
1answer
39 views

Raster subset error in R: “Error: object of type 'S4' is not subsettable”

I am running this code in R: country_template <- raster() dimensions <- extent(124000, 463000, 2953600, 3127600) country_template <- setExtent(country_template, dimensions) res(...
1
vote
1answer
25 views

projectraster in R cannot complete transformation

I have a few .tif rasters (ELEV and POP) projected in long/lat, but I would like the CRS and extent to be the same as a shapefile I have (STUDYAREA), for example, my study area lies mostly in UTM16 so ...
1
vote
1answer
40 views

Creating adjacent matrix with polygons in shapefile with R

I´m very new at these geospatial data analysis and manipulation in R. Right now I have a shapefile that you can find here: http://sina.conagua.gob.mx/sina/tema.php?tema=cuencas&ver=mapa&o=0&...
0
votes
0answers
14 views

A way to add Terrain type coefficient (as a raster) in functions of movecost package R?

Im using movecost R package, it allows to add changes in Terrain Type factor as a coefficient (N). I need this to vary between the several terrain types in my study area. So, can N be a raster with a ...
-1
votes
1answer
21 views

“.shp” extension missing in R while other layers are present

I am trying to read a shapefile in R using the command: rivnet <- readOGR(dsn = 'sa_riv_30s',layer = 'sa_riv_30s') However, the '.shp' layer does not seem to be recognized by R. The layers present ...
0
votes
1answer
23 views

How do I get a raster layer containing cell centroid coordinates (x,y) of another raster using R?

I have a rasterStack object in R and a set of polygons that represent different landcover classes in that area. For all the cells that fall in those polygons, I have extracted pixel values for each ...
0
votes
1answer
20 views

read_sf is truncating/rounding coordinates

I am using this function read_sf from package sf to read a shapefile in CRS 'ETRS89 / UTM zone 29N'. But coordinates are rounded to the meter after reading it. Trying to make a reproducible example I ...
0
votes
0answers
24 views

Getting a rasterized object from irregular lat and long data

I have a dataframe (I've named it ozone2raster) with latitude and longitude points, and values for ozone concentration that I want to use for making a rasterized map. The dimensions of this dataframe ...
1
vote
1answer
43 views

Error in writeOGR when trying to export spatial polygons as shapefile in R

I am trying to export polygons that i have created for home range data but when i go to export them get this error message: Error in writeOGR(colkud95[1, ], ".", "colonyhomerange", ...
1
vote
1answer
26 views

Converting a SpatialPointsDataFrame into km units from m, in R

I read in a UTM shapefile to R using readOGR. I want to use the resulting SpatialPointsDataFrame to build a mesh in the INLA package, but I've been told that to help with numerical stability, I should ...
0
votes
1answer
34 views

How to merge (mosaic) images that have some similar parts in their file names in R

I have multiple GeoTIFF files each representing a single tile of MOD13A3 NDVI data. My study area is covered by 2 tiles and I have a time series spanning 60 months. I want to mosaic(merge) every two ...
1
vote
1answer
33 views

R: Should I transform my Raster to spatial objects here?

This question is a follow-up based on some of the advice provided in the answers and comments at Converting raster to spatial points without loss of information in R. First, to provide a bit of ...
2
votes
1answer
33 views

How do I change extent of shapefile into lat lon of NetCDF file?

I have loaded a shapefile and NetCDF into R, and I want to extract data based on the extent of the shapefile. I have done the following, but when I want to extract it I get NA's as the extent differs ...
1
vote
1answer
40 views

Create new polygon files from existing vector layer in R

I have an existing vector shapefile with all the counties in the United States (called counties.shp) I also have a dataset that contains many species-county presence data points (species.csv) which ...
2
votes
1answer
253 views

Error importing TIFF file using raster in R: “Failure during raster IO”

I've been trying to import a TIFF formatted raster image into R using the raster::raster() function with no luck, and I would like to know if the file is corrupted, or if there is something I am ...
0
votes
1answer
35 views

Eurostat shapefiles plotting with errors

I am trying to plot interactive maps at the NUTS3 level. Initially I was using shapefiles that I constructed myself for the data in QGIS, however, when I tried to use them in R they rendered poorly ...
3
votes
1answer
63 views

error running readLAS() in lidR: 'new_proj_and_gdal' is not an exported object from 'namespace:rgdal'

I've just updated lidR to 3.1.1 and am now unable to load a .las file using readLAS(): las <- lidR::readLAS(paste(sfm_dir,"/","...
1
vote
0answers
21 views

Proper way to assign CRS in rgdal 1.5-23 [duplicate]

Perhaps my question has been already answered but if it did, I haven´t be able to find it. I´ve updated the rgal package to version 1.5-23 (GDAL 3.2.1 & PROJ 7.2.1) and when I run the following ...
0
votes
1answer
169 views

ogrListLayers cannot open data source

I recently updated R and all of it's packages, and my code for reading in the layers of a geodatabase no longer works. It works fine on my windows desktop, but doesn't work on my macbook (I change the ...
0
votes
0answers
15 views

Changing hasNoDataValue from TRUE to FALSE by using rgdal

I was wondering if there is a way in R to change hasNoDataValue from TRUE to FALSE. Basically I have a raster that, when using rgdal::GDALinfo('myraster.tif'), I can see that hasNoDataValue" = ...
0
votes
0answers
69 views

Running GDAL function in R to read File Geodatabase

I was able to install GDAL through Conda and run readOGR with the following steps. Programs needed: R (download available here with RStudio: https://rstudio.com/products/rstudio/download/) Anaconda3 (...
0
votes
1answer
177 views

Subscript out of bounds error while mosaicking TIFFs together

I have just downloaded NDVI imagery from Google Earth Engine. The image was too big so it was divided into about 35 .tif files. Each tile has approximately 950 bands of NDVI values. I am now trying to ...
1
vote
0answers
29 views

Is it better to calculate area of a global raster using long/lat or using an equal area projection (e.g. Galls-Peters)?

I am trying to calculate the area of values in a global spanning raster with values of NA, 1, 2 or 3. I would like to know the total area of each category, however, the map is currently in lat long (...
1
vote
1answer
422 views

Unable to library “rgdal” and “gdalUtils” after upgrading to macOS big sur

I just upgraded my laptop to macOS big sur, but then I could not library "rgdal" and "gdalUtils". It threw me an error: Error: package or namespace load failed for ‘rgdal’ in dyn....
0
votes
1answer
462 views

Unable to find inherited method for function ‘crop’ for signature ‘“numeric” using R

I'm trying to iterate a crop and then a resampling process by using a for loop but I've faced some problems because I'd obtain an error message. Below is my code, then the error message #find all tifs ...
0
votes
0answers
373 views

Why does my rgdal discard datums?

I'm using rgdal to process rasters and shapefiles that contain CRS information [CH1903+ LV95 (EPSG2056) or CH1903 LV03 (EPSG:21781)]. For some reason part of the CRS information is always discarded (...
1
vote
1answer
75 views

Reading or creating empty shapefile in R

I have several polygon shapefiles that I need to merge into one using a loop in R. I want to start with an empty shapefile and merge it with the existing shapefiles. All shapefiles will be chacked ...
2
votes
2answers
392 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
18 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 +...
8
votes
4answers
698 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
191 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
191 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 ...
2
votes
1answer
64 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 = "...
2
votes
0answers
110 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
287 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
294 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
141 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
299 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
228 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. ...
2
votes
1answer
910 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
1answer
109 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
269 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
765 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
323 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(...
8
votes
1answer
699 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
118 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
29 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
81 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
403 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
189 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()...