GDAL

Conversion

Batch convert coordinates

cat swiss.csv | gdaltransform -s_srs EPSG:2056 -t_srs EPSG:4326 -output_xy &> swiss.dat

Change CRS

gdalwarp swissalti3d_2021_2666-1244_0.5_2056_5728.tif erdmannli_terrain.tif -s_srs EPSG:2056 -t_srs EPSG:4326

Build pyramid

gdal_retile.py -v -r bilinear -levels 4 -ps 2048 2048 -co "TILED=YES" -co "COMPRESS=JPEG" -targetDir bmpyramid Piri-Reis-1513A_cut_modifiziert.tif

Merge

gdal_merge.py -o erdmannli_ortho.tif swissimage-dop10_2021_266*.tif

Clip

gdalwarp -te 8.314 47.347 8.315 47.348 swissimage_4326.tif swissimage_cut.tif

XYZ to tiff

# Remove leading spaces
sed 's/^ *//' 003881_XYZ_2018_M151_EM122_Azores_V1.xyz > 003881_XYZ_2018_M151_EM122_Azores_V1_converted.xyz
# multiple spaces to comma
tr -s '\t ' ',,' < 003881_XYZ_2018_M151_EM122_Azores_V1_converted.xyz > 003881_XYZ_2018_M151_EM122_Azores_V1_converted2.xyz