Michael Sumner

mesh_plot from raw data

1

Nice example here https://gis.stackexchange.com/questions/321144/convert-xyz-to-stars-curvilinear-grid-and-plot#321144 ``` r f print.quosures rlang ```  ``` r im Warning in spex_to_pt(spx): loc CRS is not set, assuming longlat #> Warning in project_spex(x, "+proj=merc +a=6378137...

```R library(quadmesh) library(ggplot2) library(raster) r

```R library(raster) #' Resample a raster, fast #' #' Equivalent to `raster::aggregate` where `fun` returns a central value, but #' without taking any care for missing values or summarizing any...

Remove with a warning. Occurs with ```R rgdal::project(cbind(-179.9999, 0), "+proj=laea") ``` came up as a problem with mars.tif, when projected to "laea" (or geocent) [mars.zip](https://github.com/hypertidy/quadmesh/files/2976113/mars.zip)

https://gist.github.com/fogleman/60ace074221806986bf9290d997e3a24 See basis in examples/stl which has some geo and some saddle function examples. Python/cairo examples here: https://twitter.com/FogleBird/status/1091538515496701953

Quadmesh is obviously a good way to do this: https://github.com/cran/sp/blob/master/src/surfaceArea.c Does hinge a little on the quad vs. triangle interpretation, so just make that clear.

Consider applying a maximum to how much work is going to occur, mesh_plot also needs to be aware that it may be adding to a plot - so the maximum...

Here's etopo2 with recent sst textured on, and the north and southern ice edge. (uses lazyraster because impatience) ```R rasterRGB

pure index mesh

2

To be more generally useable `quadmesh` should provide the index, and allow the triangulated version as an alternative. So a ncol*nrow grid of any kind has two index options: quads,...

hurrah https://github.com/tylermorganwall/raybevel