lawn (version 0.6.0)

lawn_intersect: Intersection

Description

Finds the intersection of two data-Polygon's and returns just the intersection of the two

Usage

lawn_intersect(poly1, poly2, lint = FALSE)

Arguments

poly1

A data-Polygon.

poly2

A data-Polygon.

lint

(logical) Lint or not. Uses geojsonhint. Takes up increasing time as the object to get linted increases in size, so probably use by default for small objects, but not for large if you know they are good geojson objects. Default: FALSE

Value

data-Polygon, data-MultiLineString, or undefined

Details

Polygons with just a shared boundary will return the boundary. Polygons that do not intersect will return NULL.

See Also

Other transformations: lawn_bezier(), lawn_buffer(), lawn_concave(), lawn_convex(), lawn_difference(), lawn_merge(), lawn_simplify(), lawn_union()

Examples

# NOT RUN {
poly1 <- '{
 "type": "Feature",
 "properties": {
   "fill": "#0f0"
 },
 "geometry": {
   "type": "Polygon",
   "coordinates": [[
     [-122.801742, 45.48565],
     [-122.801742, 45.60491],
     [-122.584762, 45.60491],
     [-122.584762, 45.48565],
     [-122.801742, 45.48565]
    ]]
 }
}'

poly2 <- '{
 "type": "Feature",
 "properties": {
   "fill": "#00f"
 },
 "geometry": {
   "type": "Polygon",
   "coordinates": [[
     [-122.520217, 45.535693],
     [-122.64038, 45.553967],
     [-122.720031, 45.526554],
     [-122.669906, 45.507309],
     [-122.723464, 45.446643],
     [-122.532577, 45.408574],
     [-122.487258, 45.477466],
     [-122.520217, 45.535693]
     ]]
 }
}'
lawn_intersect(poly1, poly2)

view(poly1)
view(poly2)
lawn_intersect(poly1, poly2) %>% view()

x1 <- lawn_buffer(lawn_point(c(-122.6375, 45.53)), 1500, "meters")
x2 <- lawn_buffer(lawn_point(c(-122.6475, 45.53)), 1500, "meters")
lawn_intersect(x1, x2)
structure(x1, class = "featurecollection") %>% view()
structure(x2, class = "featurecollection") %>% view()
lawn_intersect(x1, x2) %>% view()

# not overlapping
x3 <- lawn_buffer(lawn_point(c(-122.6375, 45.53)), 1500, "meters")
x4 <- lawn_buffer(lawn_point(c(-122.6975, 45.53)), 1500, "meters")
structure(x3, class = "featurecollection") %>% view()
structure(x4, class = "featurecollection") %>% view()
lawn_intersect(x3, x4)
# }