array_sum(
array_foreach(
generate_series(0, layer_property('MY_LAYER', 'feature_count'), 1),
if(attribute(get_feature_by_id('MY_LAYER', @element), 'MY_FIELD_TO_FILTER') = 'MY_FILTER_VALUE',
attribute(get_feature_by_id('MY_LAYER', @element), 'MY_FIELD_TO_AGGREGATE'),
0
)
)
)aggregate(
layer:='MY_LAYER',
aggregate:='sum',
expression:="MY_FIELD_TO_AGGREGATE",
filter:="MY_FIELD_TO_FILTER"='MY_FILTER_VALUE'
)
https://heatmap-external-a.strava.com/tiles-auth/ride/hot/{z}/{x}/{y}.png?Key-Pair-Id=LONGSTRING&Signature=LONGERSTRING
array_sum(
array_filter(
array_foreach(generate_series(1,num_geometries(nodes_to_points($geometry))),
z(geometry_n( nodes_to_points($geometry), @element+1))-
z(geometry_n( nodes_to_points($geometry), @element))
), @element>0
)
)
color_mix_rgb('blue', 'white', 0)'0,0,255,255'color_mix_rgb("color_name_code", 'white', 0.2)'DGM_merged'
array_min(
array_foreach(
generate_series(1, num_points($geometry), 1),
raster_value('DGM_merged', 1,
point_n(nodes_to_points($geometry), @element)
)
)
)
SELECT
t.Name AS Trail_Name,
g.GNR AS Parcel_NR,
SUM(LENGTH(INTERSECTION(t.geometry, g.geometry))) AS SUM_L_inP,
COLLECT(INTERSECTION(t.geometry, g.geometry)) as geom
FROM Trails AS t JOIN Parcels AS g ON INTERSECTS(t.geometry, g.geometry)
GROUP BY g.GNR
if(
regexp_substr(description, 'img src="!?.*img src="(.*)"') = '',
regexp_substr(description, 'img src="(.*)"'),
regexp_substr(description, 'img src="!?.*img src="(.*)"')
)
#!/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 would throw an error and discard attribute fields for all fieldnames that do not match the GPX XML definition (name, cmt, etc…).
## if you dont put the GPX_USE_EXTENSION=YES paramter..
## with <-t_srs epsg:4326=""> the target CRS, to which the input coordinates will be transformed, is given..
OIFS="$IFS"
IFS=$'\n'
for f in *.gpkg
do ogr2ogr -f GPX -dsco GPX_USE_EXTENSIONS=YES -t_srs EPSG:4326 -overwrite ${f%.*}.gpx $f
done
IFS="$OIFS"
read -p "Hit [Enter] to exit..."
array_foreach()array_filter()collect_geometries()
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)
generate_seriesaggregatearray_foreacharray_filtercollect_geometrieswith_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'
)
)
)
regexp_replace(aggregate(layer:='Prcl_Lyr', aggregate:='concatenate', expression:="Prp_Ownr",filter:=intersects( $geometry, geometry(@parent)), concatenator:='/'), '[/]{2,}', '/')aggregate(layer:='Corridors_86ef2ea9_d1a1_4d9f_9735_7b7b2fb54cb2',
aggregate:='concatenate',
expression:="Name", filter:=within($geometry, map_get( item_variables('Main Map'),
'map_extent')), concatenator:=', ')
aggregate(layer:='BIKFFH_Karwendel BIKFFH_PL', aggregate:='sum', expression:="LNUMMER",
filter:=intersects($geometry,geometry(@parent)))
from PyQt4.QtCore import *
from qgis.core import *
from qgis.utils import *
#===============================================
##[User scripts]=group
##Gst_Nr=string
##Name_Layer=string 81127GST_V2
#===============================================
# Gst_Nr is the field in which we are searching
# 81127GST_V2 is the default layer for searching
#===============================================
canvas = iface.mapCanvas()
# First zoom to desired scale
canvas.zoomScale( 2500 )
allLayers = canvas.layers()
n = len(allLayers)
for i in range(0, n):
if allLayers[i].name() == Name_Layer:
break
tarL = allLayers[i]
# Get a featureIterator from an expression
expr = QgsExpression( "\"GNR\"='" + Gst_Nr+ "'" )
it = tarL.getFeatures( QgsFeatureRequest( expr ) )
# Build a list of feature Ids from the result obtained above
ids = [j.id() for j in it]
# Select features with the ids
tarL.setSelectedFeatures( ids )
# Zoom to selected features
canvas.zoomToSelected( tarL )
SELECT t.gID AS gID, t.sp AS Sp, count(*) AS NrIndSp FROM (SELECT g.id AS gID, s.sp AS sp FROM grid AS g JOIN Sp_distr AS s ON within(s.Geometry, g.Geometry) ) as t GROUP BY t.gId, t.sp
Select t.gID, count(*) as NrInd From (SELECT g.id AS gID, s.sp AS sp FROM grid AS g JOIN Sp_distr AS s ON within(s.Geometry, g.Geometry) ) as t GROUP BY t.gID ORDER BY t.gID
SELECT v.gID, count(*) AS SpNr FROM (SELECT t.gID, t.sp FROM (SELECT g.id AS gID, s.sp AS sp FROM grid AS g JOIN Sp_distr AS s ON within(s.Geometry, g.Geometry) ) as t GROUP BY t.gId, t.sp ) as v GROUP BY v.gId
select
name as name,
makeLine(t.ptgeom, 0) as geom
from (
select
name as name,
DissolvePoints(Collect(ST_Reverse(geometry))) as ptgeom
from segments group by name )
as t
select
replace(replace(replace(replace(replace(replace(astext(Collect(t.geometry)), 'MULTILINESTRING((','§'), '))', '%'), '(', ''), ')', ''), '§', 'LINESTRING('), '%', ')'
) as geom
from (
select MultiLinestringFromText('MULTILINESTRING((-1 -1, 0 0), (1 1, 4 4))') as geometry
) as t
LINESTRING(-1 -1, 0 0, 1 1, 4 4)
LINESTRING(-1 -1, 0 0, 0 0, 4 4)
SELECT AsText(Linemerge(MultiLinestringFromText('MULTILINESTRING((-1 -1, 0 0), (0 0, 4 4))')))
LINESTRING(-1 -1, 0 0, 4 4)
select
g1.PKUID as gID,
count (distinct s.species) as sp_num_inbu,
g1.Geometry AS geom
from (
select g.*
from(select Gunion(geometry) as geom
from Sp_distr) as u, grid as g
where area(intersection(buffer(centroid(g.geometry), 500), u.geom)) > pow(500, 2)*pi()*0.8
) as g1 join Sp_distr as s on intersects(buffer(centroid( g1.Geometry), 500), s.Geometry)
group by gID
select
g.*
from(select Gunion(geometry) as geom
from MYPLGN) as u, grid as g
where area(intersection(g.geometry, u.geom)) = area(g.geometry)
SELECT
ID,
makeline(Geometry) AS geom
FROM Birds
GROUP BY ID
SELECT
t.region AS region,
t.species AS sp_dom,
count(*) AS sp_number,
max(t.sp_area) / 10000 AS sp_dom_area
FROM (
SELECT
g.region AS region, s.species AS species,
area(intersection(g.Geometry, s.Geometry)) AS sp_area
FROM Regions AS g JOIN Sp_Distribution AS s
ON INTERSECTS(g.Geometry, s.Geometry)
) AS t
GROUP BY t.region
ORDER BY t.region
SELECT
g.region AS regID,
s.species AS sp,
AREA(INTERSECTION(g.geometry, s.geometry)) AS sp_area
FROM Regions AS g JOIN Sp_Distribution AS s
ON INTERSECTS(g.Geometry,s.Geometry)
ORDER BY regID, sp_area ASC
SELECT g.plgnID AS "plgn_ID",
AVG(s.pollution) AS "Average Pollution",
MAX(s.pollution) AS "Maximum Pollution",
COUNT(*) AS "Number of Sensors"
FROM sensors AS s JOIN SHAPE1 AS g
ON contains(g.geometry, s.geometry)
GROUP BY g.plgnID
SELECT t.geom AS geom,
t.plgnID AS plgnID,
t.soiltype AS soiltype,
max(t.soil_area) AS MaxArea, count () AS n_soiltypes
FROM (SELECT
g.Geometry AS geom, g.plgnID AS plgnID, s.Soiltype AS soiltype,
AREA(INTERSECTION(g.geometryO, s.geometry)) AS soil_area
FROM SHAPE1 AS g JOIN Soils AS s
ON INTERSECTS(g.Geometry,s.Geometry)
) AS t
GROUP BY t.plgnID
ORDER BY t.plgnID
Set WshShell = CreateObject("WScript.Shell")
WshShell.Run chr(34) & "C:\OSGeo4W\bin\qgis.bat" & Chr(34), 0
Set WshShell = Nothing
gdal_contour -fl 900 "path/to/dem_raster.asc" "path/to/output.shp"
