Krigeage avec barrières dans ArcGIS Desktop?

10

J'essaie d'interpoler dans une zone côtière avec des points de données dans plusieurs rivières se vidant dans une baie et sur les côtés opposés des péninsules.

J'ai créé une couche côtière pour l'utiliser comme barrière, mais je ne sais pas comment procéder dans ArcGIS 10.

Je ne vois aucun objet barrières dans aucun des outils de krigeage.

J'ai essayé de faire le krigeage puis de le couper au sol, mais j'obtiens des résultats irréalistes où les points de données rapprochés mais séparés par le sol créent des problèmes.

Je vois la documentation en ligne pour ce faire dans les scripts de géotraitement en 9.3 mais rien de similaire pour 10.

— elayne927
source

Pour ajouter des liens pour afficher les modifications, comparez la documentation d' ArcGIS 9.3 à ArcGIS 10

— Mike T

3

Les rivières coulent-elles généralement dans une direction spécifique? Si tel est le cas, vous pouvez utiliser le regroupement directionnel le long des lignes centrales des rivières pour limiter les semi-variogrammes empiriques uniquement aux paires de points s'étendant sur la longueur des rivières. (Voir ici pour plus de détails: http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/Exploring_directional_autocorrelation/0031000000m2000000/ )

Une autre option serait de partitionner les échantillons ponctuels en ensembles de données séparés et d'exécuter votre krigeage sur chaque ensemble séparément. Cela pourrait être plus approprié de toute façon, car chaque rivière pourrait avoir des modèles statistiques différents; mais ne fonctionnerait pas si la taille de votre échantillon est petite pour chaque rivière.

Avez-vous besoin d'utiliser une méthode géostatistique ou une méthode déterministe fonctionnerait-elle? Je pense que vous pourriez être en mesure d'atteindre ce que vous voulez en utilisant l'interpolation de diffusion avec des barrières. Vous pouvez utiliser votre terrain comme une barrière extrêmement coûteuse. http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//00310000002s000000.htm

En outre, voici l'interpolation de diffusion avec des barrières absolues comme outil de script de géotraitement. Votre littoral ferait une barrière absolue parfaite. http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/Diffusion_Interpolation_With_Barriers/003000000005000000/

— blord-castillo
source

0

Avez-vous essayé l'outil Topo to Raster? Vous pouvez spécifier le type des entités en entrée (altitude du point, ruisseau, puits, contour, limite, lac.) Vous pouvez essayer de définir les entités de votre couche de ligne de côte sur "TopoLake"

En outre, le résultat des outils Topo vers raster est un raster "hydologiquement correct"

Citations de l'aide: TopoLake Classe d'entités surfaciques qui spécifie l'emplacement des lacs. Toutes les cellules raster en sortie d'un lac seront affectées à la valeur d'élévation minimale de toutes les cellules le long du rivage.

— Jakub Sisak GeoGraphics
source

Je pense que plutôt que de construire une représentation de la topographie, le PO souhaite créer une surface de mesures, telles que la qualité de l'eau, effectuées à l'intérieur des plans d'eau. Topo2Raster ne fera pas ça. (Ce n'est pas Krige non plus.)

— whuber

Oui, j'essaie de faire exactement cela.

— elayne927

Vous avez absolument raison Bill, Topo to Raster ne permet pas le krigeage. J'avais juste l'impression qu'Elayne faisait référence à l'interpolation en général. (J'entends souvent des gens regrouper toute interpolation sous krigeage.)

— Jakub Sisak GeoGraphics

Elayne, je pense que je sais ce que vous entendez par résultats irréalistes. Vous pouvez essayer de convertir vos lacs en une grille de points régulièrement espacés avec la valeur d'élévation du lac affectée à chaque point et de l'ajouter en entrée à l'outil Topo vers raster en tant qu'élévation de point. Cela ne donnera pas non plus un résultat parfait, surtout le long du rivage mais apparaîtra plus réaliste. Avez-vous essayé le TIN et les lignes de rupture, puis converti le TIN résultant en raster? Pouvez-vous publier le lien vers les scripts de géotraitement auxquels vous faites référence?

— Jakub Sisak GeoGraphics

0

Je suis un adopteur tardif, je ne peux donc que vous dire ce qui fonctionne en 9.3 avec l'extension Spatial Analyst - publier juste au cas où la méthode fonctionne toujours en 10. Tout d'abord, créez une chaîne fermée représentant votre domaine d'intérêt (sonne comme l'inverse de votre littoral ici). Spécifiez votre "conteneur" dans la boîte de dialogue des options de Spatial Analyst sous l'onglet d'étendue. Effectuez votre analyse en sélectionnant Krigeage dans Interpoler vers raster dans le menu Spatial Analyst.

Bonne chance!

— rec.thegeom
source