GDAL: Bash Script for Converting All GPKG-Files From Directory To GPX-Format

Sinve lately QGIS (I use QGIS 3.12) comes equipped with a nice tool ("Split Vector Layers"), that will split your layers based on an attribute an name your files accordingly. THe only drawback is, that you can not choose different output formats, but need to go with the one and only GPKG-format. Now, I needed GPX in my case and had to find a way of converting the gpkg-files in batch mode. This a achieved by writing a small bash script looping over the files with GDAL's ogr2ogr command, and convert each file to GPX.

#!/bin/bash
#first two lines and last line prevent whitspace problems in filenames (see: https://unix.stackexchange.com/questions/9496/looping-through-files-with-spaces-in-the-names
#then loop over all files with gpkg extension an convert to gpx format preserving the old filename!
#the script will throw an error and discard attribute fields for all fieldnames that do not match the GPX XML definition (name, cmt, etc…). 
#with <-t_srs epsg:4326=""> the target CRS, to which the input coordinates need to be transformed, is given..

OIFS="$IFS"
IFS=$'\n'
for f in *.gpkg
	do ogr2ogr -t_srs EPSG:4326 -overwrite ${f%.*}.gpx $f
done
IFS="$OIFS"
read -p "Hit [Enter] to exit..."

QGIS 3: Symbology with Geometry Generator - Draw One Convex Hull For All Features With Same Attribute

I have a layer name "Trails" with an attribute "Trail ID", which contains uniquie, consecutive Feature/Trail IDs, and an attribute "Schwierigkeit" (trail difficulty). For all features with the value "black" for the attribute "Schwierigkeit" I want to render one convex hull. The below code will select the last element in of all the features in the layer and apply the code for drawing the polygon only once. The first

array_foreach()

in the code will create an array of all features (series generated from 1 to feature count number), over this array, the second array_foreach() will apply the geometry finction on each element which meets the condition of the

array_filter()

function. The

collect_geometries()

function finally puts all those single geometries within the resulting array into one multiline geometry, for which I then draw the hull.

if($id = maximum($id),
convex_hull(
collect_geometries(
with_variable('my_arr', 
array_foreach(
generate_series(1,  layer_property( 'trails', 'feature_count'), 1),
get_feature('trails', 'Trail ID', @element)
),
array_foreach(
array_filter(@my_arr, attribute(@element, 'Schwierigkeit')='black') , geometry(@element)))))
, NULL)

QGIS 3: Geometry Generator Expression for Points at Line Intersections

A solution for symbolizing line intersections between the current layer's features and another layer's features ('Tracks Subset=Detail'). The

generate_series

is used with an

aggregate

which yields the number if features in the other layer used for the intersection. The

array_foreach

iterates over all features of the other layer and intersects each line with the input layer's current feature. The

array_filter

is needed to filter out the epmty geometries that could result after the intersection. If you wouldn't do this, the code would break when calling

collect_geometries

, which is needed to finally convert the array of geometries into a valid, digestible geometry collection, fed into the geometry generator with checking Point/Multi-POint as geometry type..

with_variable ('my_series', generate_series(1, aggregate(layer:='Tracks Subset=Detail', aggregate:='max', expression:="id")),
	collect_geometries(
		array_filter(
		array_foreach(@my_series, if(is_empty_or_null( 
			intersection(geometry(get_feature('Tracks Subset=Detail', 'id', @element)),
			$geometry)), 'x', intersection(geometry(get_feature('Tracks Subset=Detail', 'id', @element)),
			$geometry))
			), 
		@element != 'x'
	)
  )
)