Pythonを用いてShapeFileをSVGに変換できる方法

## ライブラリをインポート
import geopandas as gpd
from bs4 import BeautifulSoup
import random
import textwrap

## ヨーロッパ評議会から取得したデータをGeoPandasデータフレームとして読み込む
fp = "geodata/TR_2011_LAU1.shp"
data = gpd.read_file(fp)

## shapefileの行をpathに変換する関数を定義
def process_to_svg_group(row):
	rd = row.to_dict()
	del rd["geometry"]
	
	# パスの属性としてデータフレームの列を保持する
	to_add = []
	for key, val in rd.items():
					to_add.append('data-{}="{}"'.format(key, val))

	# geometry列をSVGパスに変換
	ps = BeautifulSoup(row.geometry._repr_svg_(), "xml").find_all("path")
	paths = []
	for p in ps:
		new_path = f"""<g {' '.join(to_add)}>{str(p)}</g>"""
		paths.append(new_path)

	return "\n\n".join(paths)

## データフレームの各行をSVGパスに変換
processed_rows = []
for i, row in data.sample(frac=1).iterrows():
	p = process_to_svg_group(row)
	processed_rows.append(p)

## SVGの属性を設定
props = {
	"viewBox": f"{data.total_bounds[0]} {data.total_bounds[1]} {data.total_bounds[2] - data.total_bounds[0]} {data.total_bounds[3] -data.total_bounds[1] }",
	"xmlns": "http://www.w3.org/2000/svg",
	"xmlns:xlink": "http://www.w3.org/1999/xlink",
	"transform": "scale(1, -1)",
	"preserveAspectRatio":"XMinYMin meet"
}

## 属性を追加
template = '{key:s}="{val:s}"'
attrs = " ".join([template.format(key=key, val=props[key]) for key in props])

## SVGを作成
raw_svg_str = textwrap.dedent(r"""
	<?xml version="1.0" encoding="utf-8" ?>
	<svg {attrs:s}>
	{data:s}
	</svg>""").format(attrs=attrs, data="".join(processed_rows)).strip()

## SVGを保存
with open("map.svg", "w") as f:
	f.write(raw_svg_str)