Edzer Pebesma

Results 708 comments of Edzer Pebesma

Why is ``` gdal(warn = 1) ``` not the default?

I was looking into `b.nc` created above, and noted it has ``` // global attributes: :Conventions = "CF-1.4" ; :created_by = "R packages ncdf4 and terra (version 1.7-37)" ; ```...

Yes, [this](https://github.com/OSGeo/gdal/blob/master/frmts/netcdf/netcdfdataset.h#L284-L288) is where GDAL looks for "valid" time unit signatures, and `years since` is not in it - no idea if this is an omission or on purpose.

> Years and months are valid units, but tricky, I read that in the udunits docs, and agree that "year" is a difficult unit, however specified dates defined as e.g....

Do `terra` objects btw have a way to register whether, e.g. in case of yearly time steps, the values are associated with a time instance (Jan 1st 00:00) or with...

Please give the entire output you see after giving the `install.packages("mapview")` command.

On a fresh colab notebook I can install `sf`, then `mapview`, but when loading mapview I get the following error: ![Screenshot from 2023-01-01 20-48-33](https://user-images.githubusercontent.com/520851/210183029-d229566e-07c0-4d42-a7a9-b6b479113c26.png) `htmltools` is required by `IRkernel`, which...

I need some kind of reprex to be able to look at this; if you want you can share data or a link through email.

This is because `template` has values over 12000; if you use `tm

Good idea ```r ... # Create the raster with the value in grid observations raster bind_cols(grid_points[,c("X", "Y"),]) |> st_as_sf(coords = c("X", "Y")) |> st_rasterize(template, align = TRUE) ``` now does...