diff options
| -rwxr-xr-x | gis-observation-map | 26 |
1 files changed, 22 insertions, 4 deletions
diff --git a/gis-observation-map b/gis-observation-map index 0d26ce0..5b75b89 100755 --- a/gis-observation-map +++ b/gis-observation-map @@ -229,6 +229,7 @@ def idx_intersects(idx_path, geometries): locFieldIdx = layer.FindFieldIndex('location', True) locations = [] + union = None feature = layer.GetNextFeature() while feature is not None: geometry = feature.GetGeometryRef() @@ -238,9 +239,20 @@ def idx_intersects(idx_path, geometries): loc = feature.GetField(locFieldIdx) locations.append(loc) break + if union is None: + union = geometry.Clone() + else: + union = union.Union(geometry) feature = layer.GetNextFeature() + + contains = True + for source_geometry in geometries: + if not union.Contains(source_geometry): + contains = False + break + ds = None - return locations + return locations, contains def find_qlr(basedir): if not basedir.is_dir(): @@ -743,6 +755,7 @@ getObservations(taxonLists, taxonRedlistCategories, searchFilter) topo_maps = ['Topografi 10', 'Topografi 50', 'Topografi 100'] if args.topo_basedir is not None: + idxContains = {} for topo in topo_maps: qlr_paths = [] d = args.topo_basedir.joinpath(topo) @@ -756,7 +769,7 @@ if args.topo_basedir is not None: idx = d.joinpath('index.geojson') if idx.is_file(): - tiles = idx_intersects(idx.as_posix(), geometries) + tiles, idxContains[topo] = idx_intersects(idx.as_posix(), geometries) for tile in tiles: t = d.joinpath(tile) qlr_path = find_qlr(t) @@ -804,7 +817,12 @@ if args.topo_basedir is not None: maxScale = 1 for topo in topo_maps: minScale = None - if topo == 'Topografi 10': + rootGroup = QgsProject.instance().layerTreeRoot().findGroup(topo) + if topo in idxContains.keys() and not idxContains[topo]: + # if the tileset doesn't entirely covers the area of interest, + # then uncheck the layer group and disable scale-based visibility + rootGroup.setItemVisibilityChecked(False) + elif topo == 'Topografi 10': minScale = 50000 elif topo == 'Topografi 50': minScale = 100000 @@ -812,7 +830,7 @@ if args.topo_basedir is not None: minScale = 10000000 if minScale is not None: - for lt in QgsProject.instance().layerTreeRoot().findGroup(topo).findLayers(): + for lt in rootGroup.findLayers(): lyr = lt.layer() lyr.setScaleBasedVisibility(True) lyr.setMinimumScale(minScale) |