Nos interesa el "polo de inaccesibilidad" ( https://en.wikipedia.org/wiki/Pole_of_inaccessibility ) de nuestro país.
Nuestro enfoque es utilizar QGIS, y descargar la línea de costa ([ngdc.noaa.gov]), y crear repetidamente "curvas de nivel de igual distancia" hasta que haya un punto que sea el más alejado de la línea de costa.
¿Sabes cómo generar un polígono que sea un polígono interior y tenga la misma distancia que el exterior?
Un ejemplo de curvas de nivel de igual distancia en la página de Wikipedia:
Tras un recordatorio de @underdark, por fin tengo el resultado. Aquí están mis pasos:
-
Obtenga el Shapefile de la costa continental
$ rizar el rizo https://www.ngdc.noaa.gov/mgg/shorelines/data/gshhg/latest/gshhg-shp-2.3.5-1.zip -o gshhg-shp.zip
[acción] utilizar QGIS para añadir un campo
value
con valor 1 al polígono de la costa continental -
Rasterizar el polígono de la costa continental
$ gdal_rasterize -a valor -tr 0.003 0.003 coastline.shp mainland.tif
-
Proximidad a la tierra firme rasterizada
$ gdal_proximity.py -values 0 -distunits GEO -of GTiff mainland.tif proximity.tif
-
Contorno basado en el resultado de proximidad
$ gdal_contour -a distancia -i 0,005 -f "ESRI Shapefile" proximity.tif proximity_contour.shp
¡Ahora tienes las líneas de contorno de la inaccesibilidad!