Has anyone figured out how to compensate for the changes in the PROJ projection engine?
I use crs and spTransform throughout my R code and reproject to reference. I am trying to Wraptoreference but I get a huge error code (listed at the end of post).
I was told to read these two articles about the changes to PROJ:
But I don’t know how to fix the problem.
Sorry if this is very vague, I can go into more details if necessary. I am still a novice at R.
Thank you, any help is appreciated.
Error in warpToReference(source = inPaths$temp, destination = outPaths$temp, :
An error might have occured. The function returned: ERROR 1: PROJ: proj_create_from_database: SQLite error on SELECT name, type, coordinate_system_auth_name, coordinate_system_code, datum_auth_name, datum_code, area_of_use_auth_name, area_of_use_code, text_definition, deprecated FROM geodetic_crs WHERE auth_name = ? AND code = ?: no such column: area_of_use_auth_nameERROR 1: PROJ: proj_create_operations: SQLite error on SELECT name, type, coordinate_system_auth_name, coordinate_system_code, datum_auth_name, datum_code, area_of_use_auth_name, area_of_use_code, text_definition, deprecated FROM geodetic_crs WHERE auth_name = ? AND code = ?: no such column: area_of_use_auth_nameERROR 6: Cannot find coordinate operations from `GEOGCRS["GRS 1980(IUGG, 1980)",DATUM["unknown",ELLIPSOID["GRS80",6378137,298.257222101,LENGTHUNIT["metre",1,ID["EPSG",9001]]]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9122]]],CS[ellipsoidal,2],AXIS["latitude",north,ORDER[1],ANGL
In addition: Warning messages:
1: In getProjectionRef(x, OVERRIDE_PROJ_DATUM_WITH_TOWGS84 = OVERRIDE_PROJ_DATUM_WITH_TOWGS84, :
Discarded datum unknown in Proj4 definition: +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs,
but +towgs84= values preserved
2: In getProjectionRef(x, OVERRIDE_PROJ_DATUM_WITH_TOWGS84 = OVERRIDE_PROJ_DATUM_WITH_TOWGS84, :
Discarded datum unknown in Proj4 definition: +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs,
but +towgs84= values preserved
3: In system(command = command, intern = TRUE, wait = TRUE) :
running command '"gdalwarp" -co TFW=YES -t_srs "+proj=aea +lat_0=23 +lon_0=-96 +lat_1=29.5 +lat_2=45.5 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs" -overwrite -te 1155355.08746798 1609595.57403268 2258335.08746798 3013055.57403268 -tr 30 30 -r bilinear