# finite-element-reprojection.R

#### mdsumner

#### Tue Jun 6 12:23:37 2017

    prj <- "+proj=omerc +lonc=165 +lat_0=-22 +alpha=23 +k=0.99984 +x_0=0 +y_0=0 +no_uoff +gamma=23 +ellps=WGS84 +units=m +no_defs"

    #install.packages("sfdct")
    library(sfdct)
    library(sf)

    ## Linking to GEOS 3.5.1, GDAL 2.1.2, proj.4 4.9.3

    data("wrld_simpl", package = "maptools")
    amap <- st_buffer(st_as_sf(wrld_simpl[, "NAME", drop = FALSE]), dist = 0)

    ## Warning in st_buffer.sfc(st_geometry(x), dist, nQuadSegs): st_buffer does
    ## not correctly buffer longitude/latitude data, dist needs to be in decimal
    ## degrees.

    ## kludges to get properties of the data as they are
    ## without geographic or metadata assumptions about what
    ## I'm trying to do
    areas_native <- as.numeric(st_area(st_set_crs(amap, NA)))
    w <- ct_triangulate(amap, a = 1)

    wp <- st_cast(st_transform(w, prj))

    ## Warning in st_cast.sf(st_transform(w, prj)): repeating attributes for all
    ## sub-geometries for which they may not be constant

    library(dplyr)

    ##
    ## Attaching package: 'dplyr'

    ## The following objects are masked from 'package:stats':
    ##
    ##     filter, lag

    ## The following objects are masked from 'package:base':
    ##
    ##     intersect, setdiff, setequal, union

    ## note that st_length is native length because it's not longlat
    ## and so not geoid-aware so we good ...
    wpf <- dplyr::filter(wp[as.numeric(st_length(wp)) < 3e7, ])
    plot(wpf, border = NA)