登記所備付地図データを使ってみる

import zipfile
from bs4 import BeautifulSoup
from shapely.geometry import Point
import geopandas as gpd

# 解凍するファイルとCRSの設定（先に調べてあるものとして）
zip_file1 = '08224-0527.zip'
zip_file2 = '08224-0527-75.zip'
crs = 'EPSG:6677'

# zipファイルの解凍
with zipfile.ZipFile(zip_file1) as zf:
    with zf.open(zip_file2) as ff:
        with zipfile.ZipFile(ff) as zff:
            with zff.open(zff.namelist()[0]) as fff:
                data = fff.read()
                soup = BeautifulSoup(data, "lxml-xml")

# 空間属性＜GM_Point＞の解析（idと座標の辞書を作成）
point = {}
for member in soup.find_all('GM_Point'):
    point[member['id']] = [float(member.find('X').text), float(member.find('Y').text)]

# 主題属性＜基準点＞の解析
pieces = []
for member in soup.find_all('基準点'):
    r = {}
    # 属性部分
    r['名称'] = member.find('名称').text
    r['基準点種別'] = member.find('基準点種別').text
    r['埋標区分'] = member.find('埋標区分').text
    
    # ジオメトリの作成
    key = member.find('形状')['idref']
    r['geometry'] = Point((point[key][1], point[key][0])) # GISではxyが逆になる

    pieces.append(r)
    
# GeoDataFrame作成
gdf = gpd.GeoDataFrame(pieces, crs=crs)

# 保存
gdf.to_file('sample1.shp', encoding='cp932')

import zipfile
from lxml import etree as ET
import io
from shapely.geometry import Point
import geopandas as gpd

# 解凍するファイルとCRSの設定（先に調べてあるものとして）
zip_file1 = '08224-0527.zip'
zip_file2 = '08224-0527-75.zip'
crs = 'EPSG:6677'

# zipファイルの解凍
with zipfile.ZipFile(zip_file1) as zf:
    with zf.open(zip_file2) as ff:
        with zipfile.ZipFile(ff) as zff:
            with zff.open(zff.namelist()[0]) as fff:
                data = fff.read()
                tree = ET.parse(io.BytesIO(data))
                root = tree.getroot()

# 空間属性＜GM_Point＞の解析（idと座標の辞書を作成）
point = {}

for child in root.findall('.//{http://www.moj.go.jp/MINJI/tizuzumen}GM_Point'):
    id = child.attrib['id']
    x = float(child.find('.//{http://www.moj.go.jp/MINJI/tizuzumen}X').text)
    y = float(child.find('.//{http://www.moj.go.jp/MINJI/tizuzumen}Y').text)
    
    point[id] = [x, y]

# 主題属性＜基準点＞の解析
pieces = []
for child in root.findall('.//{http://www.moj.go.jp/MINJI/tizuxml}基準点'):
    r = {}
    # データ取得部分
    for c in child:
        if '形状' in c.tag: # 形状の要素のみidrefを取得する必要がある
            v = c.attrib['idref']

            # ジオメトリ作成
            r['geometry'] = Point((point[v][1], point[v][0])) # GISではxyが逆になる      
        else:
            v = c.text

        k = c.tag.replace('{http://www.moj.go.jp/MINJI/tizuxml}','')
        r.setdefault(k, v)

    pieces.append(r)
    
# GeoDataFrame作成
gdf = gpd.GeoDataFrame(pieces, crs=crs)

# 保存
gdf.to_file('sample2.shp', encoding='cp932')