Dates incorrect for ESRI feature service

[bengannon-fc]

Date fields do not appear to read in correctly for ESRI feature services. The same data is read in correctly by the sf package.

# Read in current fire perimeters from the National Interagency Fire Center
poly_url <- 'https://services3.arcgis.com/T4QMspbfLg3qTGWY/arcgis/rest/services/WFIGS_Interagency_Perimeters_Current/FeatureServer/0/query?where=1=1&returnGeometry=true&outFields=*&f=geojson'

# With sf as reference
library('sf')
service_sf <- st_read(poly_url)

# With terra to illustrate date issue
library('terra')
service_terra <- vect(poly_url)

# sf and terra dates differ - appear to use different origins
summary(service_sf$poly_CreateDate)
summary(service_terra$poly_CreateDate)
plot(service_sf$poly_CreateDate,service_terra$poly_CreateDate,
     xlab='sf date number',ylab='Terra date number')

# sf dates are correct
dates_sf <- as.POSIXct(service_sf$poly_CreateDate/1000,origin='1970-01-01')
dates_sf[1:5]
# [1] "2024-09-16 09:31:07 PDT" "2024-08-01 15:12:03 PDT"
# [3] "2024-07-07 16:52:58 PDT" "2024-07-23 09:22:27 PDT"
# [5] "2024-09-18 05:30:00 PDT"

# terra dates are incorrect
dates_terra <- as.POSIXct(service_terra$poly_CreateDate/1000,origin='1970-01-01')
dates_terra[1:5]
[1] "1969-12-30 19:50:14 PST" "1970-01-03 18:33:57 PST"
[3] "1969-12-09 20:14:52 PST" "1969-12-25 12:44:21 PST"
[5] "1970-01-01 15:49:07 PST"

Terra version information:

packageVersion("terra")
#[1] ‘1.7.46’
terra::gdal(lib="all")
#    gdal     proj     geos 
# "3.6.2"  "9.2.0" "3.11.2"

[rhijmans]

Thank you for reporting that. Fixed:

poly_url <- 'https://services3.arcgis.com/T4QMspbfLg3qTGWY/arcgis/rest/services/WFIGS_Interagency_Perimeters_Current/FeatureServer/0/query?where=1=1&returnGeometry=true&outFields=*&f=geojson'

library('sf')
service_sf <- st_read(poly_url)
library('terra')
service_terra <- vect(poly_url)
dates_sf <- as.POSIXct(service_sf$poly_CreateDate/1000,origin='1970-01-01')
dates_sf[1:5]
#[1] "2024-09-16 09:31:07 PDT" "2024-08-01 15:12:03 PDT" "2024-07-07 16:52:58 PDT" "2024-07-23 09:22:27 PDT" "2024-09-18 05:30:00 PDT"

dates_terra <- as.POSIXct(service_terra$poly_CreateDate/1000,origin='1970-01-01')
dates_terra[1:5]
#[1] "2024-09-16 09:31:07 PDT" "2024-08-01 15:12:03 PDT" "2024-07-07 16:52:58 PDT" "2024-07-23 09:22:27 PDT" "2024-09-18 05:30:00 PDT"