¿Cómo interpolar mientras se proporciona la geometría de recorte? (ArcGIS 10)

Question

¿Hay alguna forma de interpolar (IDW) mientras se proporciona la geometría de recorte? Mi objetivo es que quiero interpolar los datos de temperatura pero no quiero que interpuele toda la extensión del archivo de puntos sino que ajuste la salida a la extensión de un polígono concreto. Gracias.

Answer 1

Ir a arctoolbox => herramientas de análisis 3d => IDW => entornos => Configuración de análisis de trama => Máscara

configuración general => extensión

o

puede utilizar extracto por máscara La herramienta .extract by mask no es una buena opción para usar la interpolación.

Answer 2

No creo que se pueda especificar una geometría de recorte como parte del proceso IDW, sin embargo se puede realizar una operación de recorte por separado ya sea en los datos de entrada o de salida y automatizar este proceso utilizando comandos de python.

Me gustaría advertir que si recortas un subconjunto de puntos y luego realizas un IDW en el subconjunto obtendrás una superficie diferente que si realizas un IDW y luego recortas el resultado a tu polígono.

A continuación se presentan ejemplos de comandos python adaptados de los archivos de ayuda de IDW , Clip y Extracto de la máscara

Clip antes de IDW

import arcpy
from arcpy import env
from arcpy.sa import *

env.workspace = "C:/data"
arcpy.Clip_analysis([POINTFIELD], [POLYGON], "studypts.shp")

outIDW = Idw("studypts.shp", "temperature", 2000, 2, RadiusVariable(10, 150000))
outIDW.save("C:/output/idw.tif")

Clip después de IDW

import arcpy
from arcpy import env
from arcpy.sa import *

env.workspace = "C:/data"

outIDW = Idw([POINTFIELD], "temperature", 2000, 2, RadiusVariable(10, 150000))
outIDW.save("C:/output/idw")

outExtractByMask = ExtractByMask("C:/output/idw", [POLYGON])
outExtractByMask.save("C:/output/maskextract")