NDVI (MOD13Q1) 5 de la serie MODIS: flujo de procesamiento NDVI

Prefacio: (Consejo personal)

1. Realice el trabajo de procesamiento de datos, porque la cantidad de datos suele ser relativamente grande. Por lo tanto, en las computadoras personales, es necesario clasificar datos básicos, datos de procesos crecientes y datos de resultados, lo que facilita su búsqueda.

2. Haga una copia de seguridad de los datos importantes, porque en el proceso de procesamiento de datos, los datos importantes se perderán debido a errores, no permita que los datos importantes se agoten. Por lo tanto, generar o guardar archivos comprimidos es una buena opción.

3. Cuando la cantidad de datos es grande y se requiere procesamiento por lotes. El documento o código de plantilla debe evitar la aparición de errores de detalle. Es fácil de completar, lento e ineficiente.

Herramienta de procesamiento: ArcMap

Procesamiento de datos: archivo de producto terminado NDVI (.tif) de MOD13Q1

Archivo vectorial del área de estudio (tomando la provincia de Henan como ejemplo)

Flujo de procesamiento:

1. Importe el archivo vectorial del área de estudio (realice el procesamiento hueco)

2. Importe el archivo .tif a ArcMap (puede ser de varias hojas, principalmente en función de la capacidad de la computadora para resistir).

Ya sea para construir una pirámide: sí

3. Proyección de definición de lote (MRT lo ha hecho en procesamiento por lotes, pero doble seguro y este paso es simple) (ArcMap tiene muchos métodos de proyección de lote único o múltiple, este artículo solo presenta uno)

Haga doble clic en la capa.

Haga clic en el sistema de coordenadas.

Haga clic en el sistema de coordenadas proyectadas

Haga clic en UTM

Haga clic en WGS1984

Hemisferio norte

Este artículo utiliza WGS 1984 UTM ZONE 49N.

Haga clic en Aplicar y luego en Aceptar. (Si hay coordenadas de uso común, haga clic derecho y haga clic en las coordenadas para seleccionar el favorito de la carpeta)

4. Extraer por máscara (el cultivo también funciona, pero se recomienda extraer por máscara) (este paso solo se puede procesar uno por uno, el procesamiento por lotes de Python aún no se ha actualizado)

Haga clic en ArcToolbox para abrir ArcToolbox, como se muestra en la figura a continuación para extraer por máscara, haga doble clic para extraer por máscara

Haga clic en ok

Para mayor comodidad, se puede eliminar el archivo .tif original Los resultados son los siguientes:

Si se requiere código para este paso, el método de código es el siguiente:

Corte por lotes el mismo archivo de máscara:

Importación ArcPy
 de ArcPy importación env
 de arcpy.sa importación * # Conjunto entorno 
arcpy.env.workspace r = " E: FAN InterpolationlKrigO " 
de Tramas . = ArcPy ListRasters ( " * " , " GRID " ) # get lista de archivos # ruta de salida 
outRasterPath = " E :. EANInterpolationKrigO_ extractlKrigO_ " para ras en de rásters:
     Imprimir outRasterPathtras # la hora de salida las licencias necesarias cualesquiera







    
    arcpy CheckOutExtension ( " Spatial " )   
    
    # Variables locales 
    inRaster ras 
    inMaskData = " E: FUJIANModelBoundary.shp "  # 掩膜 文件 
    outRaster-outRasterPath-ras + " .tif " 
    
    # Proceso: Extraer por máscara 
    arcpy.gp. ExtractByMask_ sa (inRaster, inMaskData, outRaster) 

print ( " OK " )