Displaying 1 result from an estimated 1 matches for "spgeom1".
Did you mean:
spgeom2
2020 Oct 20
1
sp:gIntersection warning message about projection
...: +proj=utm +zone=16 +datum=WGS84 +units=m +no_defs
variables : 2
names : x, y
value : 560615.564407584, 4698904.71208346
Please observe that both files have identical projections.
p1 <- gIntersection(buff, pln)
Warning message:
In RGEOSBinTopoFunc(spgeom1, spgeom2, byid, id, drop_lower_td,
unaryUnion_if_byid_false, :
spgeom1 and spgeom2 have different proj4 strings
This warning is baffling since both files have same projection:
crs(buff)
CRS arguments:
+proj=utm +zone=16 +datum=WGS84 +units=m +no_defs
crs(pln)
CRS arguments:
+proj=utm +zone=1...