Overlay med GeoPandas¶
Viser bruk av GeoPandas overlay-funksjon
In [1]:
import pandas as pd
import geopandas as gpd
import fiona as fi
FKB ledning¶
In [2]:
ledning_gdb = "Basisdata_1133_Hjelmeland_5972_FKB-Ledning_FGDB.gdb"
print(fi.listlayers(ledning_gdb))
In [3]:
ledning = gpd.read_file(ledning_gdb, layer = 'fkb_ledning_beliggenhet')
In [4]:
ledning.head()
Out[4]:
FKB AR5 - Arealtype 30 (skog)¶
In [5]:
ar5_gdb = "Basisdata_1133_Hjelmeland_25832_FKB-AR5_FGDB.gdb"
print(fi.listlayers(ar5_gdb))
In [6]:
ar5 = gpd.read_file(ar5_gdb, layer = 'fkb_ar5_omrade')
In [7]:
skog = ar5[(ar5["arealtype"] == "30")]
skog.head()
Out[7]:
Datasjekk¶
In [8]:
ledning.is_valid.value_counts()
Out[8]:
In [9]:
# Denne gir feilmelding:
#skog.is_valid
In [10]:
# Denne koden viser at det er tre polygoner som ikke er gyldige:
for index,value in enumerate(skog.geometry):
try:
if value.is_valid:
pass
except:
print('index: ', index)
Oppretting av feil i datagrunnlag¶
Triks for å få etablert geometrien for ugyldige polygoner på nytt:
Kjører en buffer-operasjon med bufferbredde 0
Dette er inspirert av dokumentasjonen for geopandas.overlay make_validSkriver tilbake til GeoDataFrame med pandas.DataFrame.iat-funksjon
In [11]:
# Finner kolonne-indeks for geometry
skog.columns.get_loc('geometry')
Out[11]:
In [12]:
for index,value in enumerate(skog.geometry):
try:
if value.is_valid:
pass
except:
skog.iat[index, 26] = value.buffer(0)
print('fixed geometry at index: ', index)
Problem-polygoner etter buffer(0)-operasjon¶
In [13]:
skog.iloc[513].geometry
Out[13]:
In [14]:
skog.iloc[518].geometry
Out[14]:
In [15]:
skog.iloc[5289].geometry
Out[15]:
Overlay ledning + skog¶
In [16]:
ledning_skog = gpd.overlay(ledning, skog, how='intersection', keep_geom_type=False, make_valid=True)
ledning_skog.shape
Out[16]:
In [17]:
ledning_skog.geom_type.value_counts()
Out[17]:
In [18]:
base = skog.plot(column='treslag', cmap='tab20c', figsize=(25,25))
ledning.plot(ax=base, color='cyan');
ledning_skog.plot(ax=base, color='red');
In [ ]: