Extraer el punto máximo de un ráster para una rejilla de 4x4 con los métodos de GDAL Python en QGIS 3

Esto puede realizarse con el método 'ReadAsArray' del módulo Python de GDAL empleando una ventana de 4x4 y almacenando cada punto donde ocurren los máximos en una memory layer. Para probar esta aproximación se creó un ráster de resolución de 0.5x0.5, con una amplitud (filas, columnas) de 20x20 y valores entre 1 y 1000 para evitar valores repetidos. El código desarrollado se expone a continuación y contiene varios print con propósito de corroboración.

Creando un ráster de distancia a la línea costera con PyQGIS 3 y GDAL/OGR

Con PyQGIS, cada centro de celda de un ráster puede ser calculado en un doble loop como un punto y, su distancia a la línea costera, evaluada mediante el método 'closestSegmentWithContext' de la clase QgsGeometry. Posteriormente, la librería GDAL puede ser usada para crear dicho ráster con estos valores de distancia en un array de filas x columnas.

Reproyectando y guardando un shapefile tipo polígono con GDAL/OGR

El codigo siguiente corrige los problemas de ejecución que presenta este otro para reproyectar un shapefile con base en el sistema de referencia espacial que es tomado de un ráster.