Comment obtenir un fichier de formes océaniques à partir d'un fichier de formes de zone terrestre très détaillé?

10

Je recherche une possibilité "d'inverser" un fichier de formes. Ce que je veux faire, c'est convertir un fichier de formes avec les continents et les îles en un fichier de formes contenant des polygones d'eau. Le fichier de formes océanique existant de la Terre naturelle est trop inexact.

En tant que contributeur OSM actif, je pensais que l'utilisation des fichiers de formes existants d'OSM me conviendrait. Ce sont beaucoup plus précis mais très bogués du point de vue SIG. Cependant, mapserver peut y faire face, donc j'ai pensé que cela devrait aller. Ces fichiers de formes, dont certains que vous pouvez trouver ici http://metro.teczno.com/#coastline sont très volumineux.

J'ai créé un fichier de formes avec un polygone avec des limites du monde (-180 à 180, -90 à 90). Ensuite, j'ai essayé de découper les zones de mon fichier de formes des limites du monde où se trouvent les zones terrestres dans l'autre fichier de formes. Je l'ai essayé avec QGIS, outils vectoriels -> fonction de différence. Après 12 heures avec 100% d'utilisation du CPU, je l'ai arrêté. Il affiche toujours une progression de 100%. La RAM n'était pas un problème.

Je cherche une autre possibilité de découper les zones terrestres d'un "polygone mondial", de préférence en ligne de commande Linux.

shapefile polygon clip Thomas
la source
2
Dans ArcGIS, l'opération que vous décrivez est connue sous le nom d'effacement. Il semble que vous recherchiez une solution open source, donc peut-être que la recherche sur ce terme fera monter quelque chose. MapWindow semble avoir un outil d'effacement mais la discussion ici suggère qu'il peut aussi avoir des limitations: mapwindow.org/phorum/read.php?4,20317
PolyGeo
2
Vous souhaiterez diviser votre fichier de formes "polygone avec limites du monde (-180 à 180, -90 à 90)" en plus petits morceaux. Ensuite, vous pouvez découper les zones. Vous rencontrez deux problèmes probables, l'un est que vous traitez simultanément des données sur l'hémisphère nord et l'hémisphère sud. Deuxièmement, des problèmes potentiels avec les limites -180/180. Essayez de diviser vos limites mondiales en 4 petits morceaux tels que: en haut à gauche, en haut à droite, en bas à gauche, en bas à droite du monde.
Michael Markieta
Je suppose, lorsque vous avez créé votre polygone mondial que vous l'avez projeté dans le même système que les données terrestres?
djq

Réponses:

6

Des rasters à la rescousse! Vous pouvez effectuer cette opération en convertissant votre fichier de formes en raster, en effectuant l'opération «inverser» et en la reconvertissant en vecteur

  • Raster -> Convertir -> Rasteriser. Raster de couche de sortie 1
  • Calculatrice de champ ouvert, entrez l'expression raster1 @ 1 = 0
  • Raster -> Convertir -> Polygoniser, Couche de sortie ocean_polygons
  • Ouvrez la table des attributs, recherchez DN = 0, basculez l'édition, supprimez les fonctionnalités sélectionnées et enregistrez. Ce qui vous restera sera des polygones pour les zones océaniques.

entrez la description de l'image ici

pensées spatiales
la source
1

SAGA GIS (et bien d'autres aussi probablement) a une commande appelée "intersection de polygones". Créez un polygone contenant le monde entier et soustrayez-y vos polygones terrestres. Donc méthode 2 ou 3 dans cette bibliothèque.

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
johanvdw
la source
0

Un raster peut être trop lourd en données si vous avez besoin de détails sur les côtes. Le même principe que celui suggéré par les réflexions spatiales mais en utilisant un TIN Delaunay pourrait être une option?

AnserGIS
la source
0

ces fichiers de formes pourraient être bénéfiques. Ils proviennent du Centre de cartographie ESRI et contiennent des données hydrologiques et marines assez précises.

Craig
la source
0

Je n'ai pas eu beaucoup de succès en inversant les fichiers de formes terrestres OSM en raison de leur complexité, mais il est possible de générer des fichiers de formes terrestres et océaniques à partir des données OSM brutes. http://openstreetmapdata.com/data propose des téléchargements des deux que je trouve extrêmement utiles.

ajashton
la source