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
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
|
"""New tests of writing feature collections."""
import pytest
from .conftest import requires_gdal33
import fiona
from fiona.crs import CRS
from fiona.errors import DriverError
from fiona.model import Feature
def test_issue771(tmpdir, caplog):
"""Overwrite a GeoJSON file without logging errors."""
schema = {"geometry": "Point", "properties": {"zero": "int"}}
feature = Feature.from_dict(
**{
"geometry": {"type": "Point", "coordinates": (0, 0)},
"properties": {"zero": "0"},
}
)
outputfile = tmpdir.join("test.geojson")
for i in range(2):
with fiona.open(
str(outputfile),
"w",
driver="GeoJSON",
schema=schema,
crs=CRS.from_epsg(4326),
) as collection:
collection.write(feature)
assert outputfile.exists()
for record in caplog.records:
assert record.levelname != "ERROR"
@requires_gdal33
def test_write__esri_only_wkt(tmpdir):
"""https://github.com/Toblerity/Fiona/issues/977"""
schema = {"geometry": "Point", "properties": {"zero": "int"}}
feature = Feature.from_dict(
**{
"geometry": {"type": "Point", "coordinates": (0, 0)},
"properties": {"zero": "0"},
}
)
target_crs = (
'PROJCS["IaRCS_04_Sioux_City-Iowa_Falls_NAD_1983_2011_LCC_US_Feet",'
'GEOGCS["GCS_NAD_1983_2011",DATUM["D_NAD_1983_2011",'
'SPHEROID["GRS_1980",6378137.0,298.257222101]],'
'PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],'
'PROJECTION["Lambert_Conformal_Conic"],'
'PARAMETER["False_Easting",14500000.0],'
'PARAMETER["False_Northing",8600000.0],'
'PARAMETER["Central_Meridian",-94.83333333333333],'
'PARAMETER["Standard_Parallel_1",42.53333333333333],'
'PARAMETER["Standard_Parallel_2",42.53333333333333],'
'PARAMETER["Scale_Factor",1.000045],'
'PARAMETER["Latitude_Of_Origin",42.53333333333333],'
'UNIT["Foot_US",0.3048006096012192]]'
)
outputfile = tmpdir.join("test.shp")
with fiona.open(
str(outputfile),
"w",
driver="ESRI Shapefile",
schema=schema,
crs=target_crs,
) as collection:
collection.write(feature)
assert collection.crs_wkt.startswith(
(
'PROJCS["IaRCS_04_Sioux_City-Iowa_Falls_NAD_1983_2011_LCC_US_Feet"',
'PROJCRS["IaRCS_04_Sioux_City-Iowa_Falls_NAD_1983_2011_LCC_US_Feet"', # GDAL 3.3+
)
)
def test_write__wkt_version(tmpdir):
"""https://github.com/Toblerity/Fiona/issues/977"""
schema = {"geometry": "Point", "properties": {"zero": "int"}}
feature = Feature.from_dict(
**{
"geometry": {"type": "Point", "coordinates": (0, 0)},
"properties": {"zero": "0"},
}
)
target_crs = (
'PROJCS["IaRCS_04_Sioux_City-Iowa_Falls_NAD_1983_2011_LCC_US_Feet",'
'GEOGCS["GCS_NAD_1983_2011",DATUM["D_NAD_1983_2011",'
'SPHEROID["GRS_1980",6378137.0,298.257222101]],'
'PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],'
'PROJECTION["Lambert_Conformal_Conic"],'
'PARAMETER["False_Easting",14500000.0],'
'PARAMETER["False_Northing",8600000.0],'
'PARAMETER["Central_Meridian",-94.83333333333333],'
'PARAMETER["Standard_Parallel_1",42.53333333333333],'
'PARAMETER["Standard_Parallel_2",42.53333333333333],'
'PARAMETER["Scale_Factor",1.000045],'
'PARAMETER["Latitude_Of_Origin",42.53333333333333],'
'UNIT["Foot_US",0.3048006096012192]]'
)
outputfile = tmpdir.join("test.shp")
with fiona.open(
str(outputfile),
"w",
driver="ESRI Shapefile",
schema=schema,
crs=target_crs,
wkt_version="WKT2_2018",
) as collection:
collection.write(feature)
assert collection.crs_wkt.startswith(
'PROJCRS["IaRCS_04_Sioux_City-Iowa_Falls_NAD_1983_2011_LCC_US_Feet"'
)
def test_issue1169():
"""Don't swallow errors when a collection can't be written."""
with pytest.raises(DriverError):
with fiona.open(
"s3://non-existing-bucket/test.geojson",
mode="w",
driver="GeoJSON",
schema={"geometry": "Point"},
) as collection:
collection.writerecords(
[
Feature.from_dict(
**{
"id": "0",
"type": "Feature",
"geometry": {"type": "Point", "coordinates": (1.0, 2.0)},
}
)
]
)
|
