Edzer Pebesma

> And that is why it would be good to do this test? The test doesn't tell you whether the coordinates are to be treated as unknown, Cartesian.

> The test would tell you that you cannot transform the data. That is equivalent to having a crs that is NA in the context of making a map. That...

of course - ``` > st_crs("ESRI:104905") Coordinate Reference System: User input: ESRI:104905 wkt: GEOGCRS["GCS_Mars_2000", DATUM["D_Mars_2000", ELLIPSOID["Mars_2000_IAU_IAG",3396190,169.894447223612, LENGTHUNIT["metre",1]]], PRIMEM["Reference_Meridian",0, ANGLEUNIT["degree",0.0174532925199433]], CS[ellipsoidal,2], AXIS["geodetic latitude (Lat)",north, ORDER[1], ANGLEUNIT["degree",0.0174532925199433]], AXIS["geodetic longitude (Lon)",east, ORDER[2], ANGLEUNIT["degree",0.0174532925199433]],...

> I've noticed that shapefiles are used quite a lot in docs of sf The original reason for that was that at the time of first submission, certain CRAN platforms...

I really should migrate this blog to quarto, I no longer have good tools to check things locally before committing / merging. Let's see how this comes out!

the warning is now gone. Now the images...

OK like this, or do you want me to update the date?

I'm volunteering - let me know when is a good time.

Is this idea still current?

Cc: @Enchufa2 - About the benefits: see [here](https://journal.r-project.org/archive/2016/RJ-2016-061/index.html): Using units may eliminate a whole class of potential scientific programming mistakes. If NetCDF files come with units, why would an R...