¿Cómo obtener un archivo de forma oceánica a partir de un archivo de forma de área terrestre muy detallado?

Estoy buscando la posibilidad de "invertir" un archivo de forma. Lo que quiero hacer es convertir un archivo de forma con los continentes e islas en un archivo de forma que tenga polígonos de agua. El archivo de forma oceánica existente de la Tierra Natural es demasiado inexacto.

Como colaborador activo de OSM, pensé que usar los shapefiles existentes de OSM estaría bien para mí. Estos son mucho más precisos pero muy defectuosos desde una perspectiva SIG. Sin embargo, el servidor de mapas puede manejarlo, así que pensé que esto debería estar bien. Estos archivos de forma, algunos de los cuales puede encontrar aquí http://metro.teczno.com/#coastline, son muy grandes.

He creado un archivo de forma con un polígono con límites mundiales (-180 a 180, -90 a 90). Luego he tratado de recortar las áreas de mi archivo de forma de límites mundiales donde las áreas de tierra están en el otro archivo de forma. Lo probé con QGIS, herramientas vectoriales -> función de diferencia. Después de 12 horas con 100% de uso de CPU, lo detuve. Siempre mostraba un progreso del 100%. La RAM no fue un problema.

Estoy buscando otra posibilidad para recortar las áreas terrestres de un "polígono mundial", preferiblemente la línea de comandos de Linux.

¡Rasters al rescate! Puede realizar esta operación convirtiendo su archivo de forma en un ráster, realizando la operación 'invertir' y convirtiéndolo de nuevo en vector

• Ráster -> Convertir -> Rasterizar. Capa de salida ráster1
• Calculadora de campo abierto, ingrese la expresión raster1 @ 1 = 0
• Ráster -> Convertir -> Poligonalizar, Capa de salida ocean_polygons
• Abra la tabla de atributos, busque DN = 0, cambie la edición, elimine las funciones seleccionadas y guarde. Lo que te quedará serán polígonos para las áreas oceánicas.

SAGA GIS (y muchos otros probablemente también) tiene un comando llamado "intersección poligonal". Cree un polígono que contenga todo el mundo y reste sus polígonos terrestres. Entonces, el método 2 o 3 en esa biblioteca.

Usage: saga_cmd libshapes_polygons "Polygon Intersection" -SHAPES_A <str> [-FIELD_A <str>] -SHAPES_B <str> [-FIELD_B <str>] -SHAPES_AB <str> [-METHOD <str>] [-SPLITPARTS]
  -SHAPES_A:<str>   Layer A
    Shapes (input)
  -FIELD_A:<str>    Attribute A
    Table field
  -SHAPES_B:<str>   Layer B
    Shapes (input)
  -FIELD_B:<str>    Attribute B
    Table field
  -SHAPES_AB:<str>  Intersection
    Shapes (output)
  -METHOD:<str>     Method
    Choice
    Available Choices:
    [0] Complete Intersection
    [1] Intersection
    [2] Difference (A - B)
    [3] Difference (B - A)
  -SPLITPARTS       Split Parts
    Boolean

Un ráster puede tener demasiados datos si necesita detalles en las líneas costeras. ¿Podría ser una opción el mismo principio que el sugerido por los pensamientos espaciales pero usar un TIN de Delaunay?

Estos archivos de forma pueden ser beneficiosos. Son del Centro de Mapeo de ESRI y contienen datos hidrológicos y de características marinas bastante precisos.

No he tenido mucho éxito con la inversión de archivos de formas terrestres OSM debido a su complejidad, sin embargo, es posible generar archivos de formas terrestres y oceánicas a partir de los datos sin procesar de OSM. http://openstreetmapdata.com/data ofrece descargas de ambos que encuentro extremadamente útiles.