Barry suggested converting lat long points e.g, by proj4stri

Some notes: Find UTM zone for conversion <p dir

Consider coordinate reference system about spocc HOT 1 CLOSED

ropensci commented on June 3, 2024

Consider coordinate reference system

from spocc.

Comments (1)

sckott commented on June 3, 2024

Some notes:

Find UTM zone for conversion

UTM reference map

This works for West longitudes only though (from SO)

long2UTM <- function(long) {
  (floor((long + 180)/6) %% 60) + 1
}
long2UTM(-149.9)

[1] 6

For all longitudes

long2utm <- function(lon, lat) {
  if(56 <= lat & lat < 64){
    if(0 <= lon & lon < 3){ 31 } else 
      if(3 <= lon & lon < 12) { 32 } else { NULL }
  } else 
  if(72 <= lat) {
    if(0 <= lon & lon < 9){ 31 } else 
      if(9 <= lon & lon < 21) { 33 } else 
        if(21 <= lon & lon < 33) { 35 } else 
          if(33 <= lon & lon < 42) { 37 } else { NULL }
  }
  (floor((lon + 180)/6) %% 60) + 1
}
long2utm(-60, 65)

21

Convert occ output lat/long coordinates to UTM

library(spocc)
dat <- occ(query='Accipiter striatus', from='gbif')
dat <- na.omit(dat@gbif@data)
coordinates(dat)=~longitude+latitude
proj4string(dat) = CRS("+init=epsg:4326")
spTransform(dat, CRS("+proj=utm +zone=11 +datum=WGS84 +units=m +no_defs +towgs84=0,0,0"))

Essentially would have to do this point by point

Since a data.frame we get back from a source we don't know what zones they are in ahead of time. Though I suppose we could do basic stuff, and just leave complicated stuff up to the user.

from spocc.

Recommend Projects