# Data accessed using command:

wget -O {output_dir}/settlements.geojson "https://kort.nunagis.gl/refserver/rest/services/Kortportal/Byer_og_bygder/MapServer/0/query/?f=json&where=OBJECTID+is+not+null&outFields=*&orderByFields=OBJECTID+ASC" && wget -O {output_dir}/towns.geojson "https://kort.nunagis.gl/refserver/rest/services/Kortportal/Byer_og_bygder/MapServer/1/query/?f=json&where=OBJECTID+is+not+null&outFields=*&orderByFields=OBJECTID+ASC"

# Data processed using the following steps:

ogrmerge.py -single -o {output_dir}/populated_places.gpkg {input_dir}/*geojson

ogr2ogr -lco ENCODING=UTF-8 -t_srs EPSG:3413 -clipdst {assets_dir}/greenland_rectangle.geojson -makevalid -nln populated_places -dialect sqlite -sql "SELECT geom, Ny_grønlandsk as 'New Greenlandic', Ny_grønlandsk as 'label', Gammel_grønlandsk as 'Old Greenlandic', Dansk as Danish, Alternativt_stednavn  as 'Alternative placename', Indbyggertal_2016 as 'Population 2016' FROM merged" {output_dir}/final.gpkg {input_dir}/populated_places.gpkg