Johannes Frölich pcace

Convert .xyz elevation models to GeoTiff

In Germany, more and more state agencies allow free access to high resolution elevation models. However, these are often released as xyz tables, which are not easily used in GIS environments. A standard method to convert this format to raster formats (eg. GeoTiff) is the GDAL function gdal_translate [1, 2]. However, converting 1M lines takes dozens of seconds and is not trivial to parallize.

This gist, using highly optimized Python library Pandas does the conversion on the same computer in less then 1 second.

	#!/bin/bash
	set -e

	# Install build dependencies
	sudo apt update
	sudo apt install -y build-essential cmake wget unzip python3-dev \
	libsqlite3-dev sqlite3 libtiff-dev libcurl4-openssl-dev \
	libxml2-dev libexpat-dev libxerces-c-dev libpq-dev \
	libfontconfig1-dev libcairo2-dev libpython3-dev swig \
	python3-pip python3-numpy

	{
	"captainVersion": "2",
	"documentation": "Based on https://www.digitalocean.com/community/tutorials/how-to-set-up-a-private-docker-registry-on-ubuntu-18-04",
	"displayName": "Docker Container Registry",
	"dockerCompose": {
	"version": "3.3",
	"services": {
	"$$cap_appname": {
	"image": "registry:2",
	"containerHttpPort": "5000",