V
A
JÊL
Fig. 3.
Versneden
vlakken met hun
centroïden
kaarten, ook van gebouwen die elkaar
slechts weinig overlappen, bijvoor
beeld als gevolg van de intrinsieke on
nauwkeurigheid van de twee kaarten;
zie fig. 2: gebouw '1193 (TOP 10-
vector) met gebouwen G272 en G437
(GBKN). De mate van overlap is dus
belangrijk, uit te drukken in twee
overlappingspercentages. Het ene per
centage is het overlappende gebied ge
deeld door de GBKN-oppervlakte, het
andere percentage is het overlappende
gebied gedeeld door de TOPlOvector-
oppervlakte. Als er een overlap wordt
gevonden tussen twee gebouwen,
maar beide overlappingspercentages
zijn lager dan een bepaalde drempel
waarde (bijvoorbeeld 5%), wordt aan
genomen dat deze gebouwen niet met
elkaar corresponderen. De overlap
pingspercentages worden berekend
door gebruik te maken van MAP-
OVER en TOPOL. Na het versnijden
van de twee kaarten en het recon
strueren van de topologie worden van
alle vlakken de bijbehorende polygo
nen berekend, alsmede alle centroïden
(fig. 3).
create table overlay jpolygons
gbkn_id char(10),
toplO_id char(10),
centroid point,
poly polygon(70))
De centroïden worden gebruikt om te
bepalen tot welk GBKN-gebouw of
TOPlOvector-gebouw het nieuwe
vlak oorspronkelijk behoorde, aange
nomen dat de berekende polygoon na
MAPOVER in de tabel „overlay_po-
lygons" staat:
update overlayjpolygons 0
from gbkn
set gbkn_id gbkn.bid
where inside(o.centroid, gbkn.geometry)
Eenzelfde update-qqery geldt voor de TOPlOvector. Met
behulp van de inhoud van tabel „overlay_polygons" is het
nu mogelijk alle correspondenties tussen gebouwen op te
sporen:
create table overlap_buildings
gbkn_id
topl0_id
gbkn_area
topl0_area
overlap_area
perc_gbkn
perc_topl0
char(10)
char(10),
float8,
float8,
float8
float8,
float8)
1
insert into overlap_buildings
(gbkn_id, topl0_id, overlap_area)
select gbkn_id, !topl0_id, sumareapoly
from overlayjpolygons
group by gbkn_icj, topl0_id
having gbkn_id ^s not null and topl0_id is not null
,,GBKN_area" en ,,TOP10_area" worden bepaald door
de oppervlakten te berekenen van de originele gebouwen
(deze query is hier niet getoond). De overlappingspercenta
ges kunnen nu worden berekend, bijvoorbeeld voor de
GBKN:
update overlap_buildings
set perc_gbkn o\/erlap_area gbkn_area 100
De combinaties geVpnden in de tabel „overlap_buildings"
moeten worden toegekend aan één van de correspondentie
categorieën (tabel 1). Dit gaat als volgt. Eerst worden com
binaties met lage pverlappingspercentages (minder dan
3%) uitgefilterd:
create view ob_5 as
select gbkn_id, topl0_id
from overlap_buildings
where perc_gbkrt> 5 or perc_topl0 5
De n-op-m-correspondenties worden gevonden door de
combinaties te seletteren uit ,,ob_5" waarvan zowel het
GBKN-gebouw als het TOPlOvector-gebouw aanwezig
zijn in andere combinaties:
create table n_tojn as
select distinct o.gbkn_id, o.topl0_id
from ob_5 o
where 1 (select count(*)
fróm ob_5 o2
where o.gbkn_j.d o2.gbkn_id)
and
1 (select count(*)
from ob_5 o3
where o.toplOTd o3.topl0_id)
113
GEODESIA
1998-3
-b-