0

I'm running a fresh R installation with RStudio. Installed packages are:

library(randomForest)
library(randomForestExplainer)
library(raster)
library(rgdal)
library(sf)

When I try loading any raster with South America Albers Equal Area Conic projection (ESPG: 102033) I receive the following warning:

Warning message: In showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj = prefer_proj) : Discarded datum Unknown based on Australian Natl & S. Amer. 1969 ellipsoid in Proj4 definition

Keeping the projection is important to me, cause I need to be able to collect point samples from that raster and I need to make sure the points and rasters are projected in the same coordinate system. Albers seems to be the most accurate for my study area.

I have packages updated to the last version availble @ cran and I cannot seem to find a solution for this.

Projection details are (quoting QGIS 3.4.4 defs):

Extent: Extent not known Proj4: +proj=aea +lat_1=-5 +lat_2=-42 +lat_0=-32 +lon_0=-60 +x_0=0 +y_0=0 +ellps=aust_SA +units=m +no_defs

I need a way to import rasters with albers projections without having them discarded.

4
  • No error there, only a warning. See gis.stackexchange.com/questions/364667/… and gis.stackexchange.com/questions/372692/…
    – JRR
    Commented Feb 25, 2021 at 11:13
  • 102033 is an ESRI projection code, not an EPSG one. The warning looks like the CRS has an ellipsoid but no datum definition, and PROJ is loathe to assume a datum. The unofficial epsg.io site has an entry ESRI:102033 which gives a datum that text-matches to two possible datums in the official epsg.org registry. Even so its possible your PROJ doesn't have datum info for those datum codes. Can you supply a raster with the coordinate info? (one with all zeroes as content will do if the data is confidential)
    – Spacedman
    Commented Feb 25, 2021 at 11:26
  • @JRR you are right, I'm not receiving an error, it's only a warning. Theoretically, I could still run my code, but there will be a point in my workflow where i'll need to point-sample my rasters with a point shapefile and I need to ensure that everything has a coordinate so the dots fall in the right place. It seems like a warning serious enough to me. Is it possible I'm overreacting?
    – Eric Lino
    Commented Feb 25, 2021 at 14:34
  • 1
    @Spacedman Weird. I did see 102033 being called both ESRI and ESPG projections and to be honest, I was a little unsure which one was right. QGIS says it's an ESPG code (imgur.com/a/dW0YmyY). As for the raster, you can find it here. Link will last 7 days: we.tl/t-bj1j6Lb8k6
    – Eric Lino
    Commented Feb 25, 2021 at 14:38

0

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Browse other questions tagged or ask your own question.