1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
|
# Copyright (c) 2020, Manfred Moitzi
# License: MIT License
import pytest
from ezdxf.entities import factory
from ezdxf.render.forms import square, translate
from ezdxf.lldxf import const
from ezdxf.addons import geo
shapely_geometry = pytest.importorskip("shapely.geometry")
def test_shapely_geo_interface():
point = shapely_geometry.shape(
{
"type": "Point",
"coordinates": (0, 0),
}
)
assert (point.x, point.y) == (0, 0)
def validate(geo_proxy: geo.GeoProxy) -> bool:
return shapely_geometry.shape(geo_proxy).is_valid
def test_resolved_hatch_with_intersecting_holes():
hatch = factory.new("HATCH")
paths = hatch.paths
paths.add_polyline_path(square(10), flags=const.BOUNDARY_PATH_EXTERNAL)
paths.add_polyline_path(
translate(square(3), (1, 1)), flags=const.BOUNDARY_PATH_DEFAULT
)
paths.add_polyline_path(
translate(square(3), (2, 2)), flags=const.BOUNDARY_PATH_DEFAULT
)
p = geo.proxy(hatch)
# Overlapping holes already resolved by fast_bbox_detection()
polygon = shapely_geometry.shape(p)
assert polygon.is_valid is True
p.filter(validate)
assert p.root["type"] == "Polygon"
assert len(p.root["coordinates"]) == 2
def test_valid_hatch():
hatch = factory.new("HATCH")
paths = hatch.paths
paths.add_polyline_path(square(10), flags=const.BOUNDARY_PATH_EXTERNAL)
paths.add_polyline_path(
translate(square(3), (1, 1)), flags=const.BOUNDARY_PATH_DEFAULT
)
paths.add_polyline_path(
translate(square(3), (5, 1)), flags=const.BOUNDARY_PATH_DEFAULT
)
p = geo.proxy(hatch)
polygon = shapely_geometry.shape(p)
assert polygon.is_valid is True
p.filter(validate)
assert p.root != {}
def test_proxy_given_hatch_when_force_line_string_generates_valid_multi_line_string() -> None:
hatch = factory.new("HATCH")
paths = hatch.paths
paths.add_polyline_path(square(10), flags=const.BOUNDARY_PATH_EXTERNAL)
paths.add_polyline_path(
translate(square(3), (1, 1)), flags=const.BOUNDARY_PATH_DEFAULT
)
p = geo.proxy(hatch, force_line_string=True)
feature = p.__geo_interface__
assert feature["type"] == "MultiLineString"
# According to the __geo_interface__ specification:
# https://gist.github.com/sgillies/2217756
assert isinstance(feature["coordinates"][0][0], tuple)
if __name__ == "__main__":
pytest.main([__file__])
|
