Importacion automatica por lotes de archivos (batch scripts) de imagenes raster a GRIDs y su procesamiento usando ArcInfo
Descripcion de herramientas disenadas para conversion rapida y procesamiento de datos raster en ArcInfo
Обсудить в форуме Комментариев — 0
Introduccion
Es posible utilizar ArcInfo Workstation para importar rapida y confiablemente datos raster que estan originalmente en formatos TIF, BSQ/BIL/BIP a un formato GRID, usado para el analisis y calculos en ArcInfo, ArcView y ArcGIS. La generacion de archivos por lotes (Batch scripts o simplemente scripts) para realizar esta tarea puede ser util tambien para ejecutar series de operaciones en GRIDs de manera automatica. Los archivos por lotes son la creacion de una lista de instrucciones para conversion o procesamiento que son interpretadas y ejecutadas automaticamente por un programa especifico, en este caso el ArcInfo.
A continuacion se ofrecen una serie de archivos por lotes que ejemplifican la capacidad de ArcInfo para realizar tareas automaticas (en Arc info a este tipo de archivos se les conoce como aml). Para ejecutar estos archivos todos los archivos raster (cualquier numero de ellos) en formato BIL o TIF deben colocarse en el mismo directorio.
Para utilizar los scripts mencionados aqui, debe cambiarse en el archivo de entrada la direccion del folder en donde todos los archivos deben ser procesados asi como el folder donde los archivos de salida se almacenaran. Si el nombre del folder de salida no se agrega, el aml usara el nombre del folder de entrada para almacenar los archivos de salida. Es tambien posible utilizar rutas relativas (Path) para el folder de salida. Por ejemplo: Si el folder de entrada esta configurado para ser c:\work\, entonces escribir result\ en el archivo por lotes significa que la ruta completa del archivo de salida seria c:\work\result\. Ambos folders, de entrada y salida, deben tener diagonal invertida al final ("\")
Ejecutando arc y configurando el ambiente de trabajo
Para comenzar ArcInfo es necesario escribir arc en una ventana de comandos o bien seleccionar Start\Programs\ArcGIS\ArcInfo Workstation\Arc o algo similar. Despues de esto se abre la consola de arc, en esta consola es necesario escribir todos los comandos manualmente, todos los comandos se escriben despues del prompt arc: . En todos los ejemplos que a continuacion mostramos se omite este prompt. Para comenzar a trabajar es mejor especificar la direccion completa del folder, donde todos los aml se localizan, al hacer esto tu debes especificar la direccion completa a estas herramientas cada vez que desees ejecutarlas. Esto puede hacerse con el comando &workspace.
Para asegurarse de que la instruccion que le dimos a arc se ejecuto correctamente se puede usar la letra w como comando en el prompt, el resultado seria similar a este: Current location: c:\work\myscripts.
Es recomendable almacenar todos los scripts en un folder, esto hace mas facil correrlos, uno des pues del otro o por separado.
Scripts
Para comenzar a trabajar con datos en formato GRID, es frecuentemente necesario convertir datos originalmente almacenados en otros formatos raster a formatos GRID, Es sencillo hacerlo con el script image2grid.aml. Para hacerlo solo hay que escribir este comando en la consola de ArcInfo Workstation:
o
&run image2grid.aml c:\work\ d:\result\
Si las imagenes raster de entrada tienen fomato 16-bit signed, entonces tendra que hacerse mas de un paso para corregirlas despues de importarlas usando con2.aml:
o
&run con2.aml c:\work\ d:\result\
El resultado de este script sera un nuevo grupo de GRIDs con la letra "c" como prefijo en el nombre, aquellos archivos temporales producidos por el script antes de este paso de importacion (usando image2grid) seran borrados del disco.
La operaciomn contraria podria ser exportar de GRID a BIL nuestros mapas. Esto tambien puede ser automatizado usando el script grid2bil.aml. Los archivos en formato BIL pueden utilizarse, por ejemplo, para cargar datos en TIMESAT.
o
&run grid2bil.aml c:\work\ ..\result\
Algunas veces es necesario enmascarar algunos valores de nuestras imagenes raster haciendolos nulos (NODATA), esta operacion es muy frecuente al manejar imagenes provenientes de sensores remotos, donde ciertos valores negativos son utilizados para los bordes de las imagenes y no tienen por tanto ninguna iformacion util. Los formatos raster mas comunes (TIF, BIL entre otros) no tienen la habilidad de almacenar datos nulos, mientras que GRID si puede. Para convertir cierto valor o intervalo a NODATA, se puede utilizar el script setnull.aml.
o
&run setnull.aml c:\work\ d:\result\
Por default este script hace nulos todos los valores menores a -3,000. Se puede cambiar esta condicion en el script usando un editor de texto y modificando la linea que contiene la condicion %current% <= - 3000, por ejemplo, si modificamos esta linea y ponemos en su lugar %current% == -3000, el resultado tendra valores nulos para los valores que originalmente eran exactamente -3000: %current% == -3000
El resultado de este script sera un nuevo grupo de GRIDs con la letra "z" como prefijo en el nombre, aquellos archivos temporales producidos por el script antes de este paso de importacion (usando image2grid.aml) seran borrados del disco.
o
%output% = setnull(%current% == -3000, %current%)
o
%output% = setnull(%current% <= -3000 and %current% >= -5000, %current%)
Para recortar todos los GRIDs usando una cobertura vectorial se puede utilizar el script clip-all.aml. Para hacerlo funcionar, la cobertura vectorial que cortara los datos necesita estar disponible y contener al menos un poligono que pueda ser usado para cortar los datos. Para mayor facilidad, esta cobertura puede colocarse en el mismo folder que los GRIDs a recortar, en este caso seria suficiente escribir en la consola directamente el nombre de la cobertura vectorial sin la ruta al archivo.
o
&run clip-all.aml clipper c:\work\data\ ..\result\
Si la cobertura vectorial esta localizada en un directorio diferente al de los GRIDs, entonces es necesario agregar la ruta completa del archivo, por ejemplo:
o
&run clip-all.aml c:\temp\clipcover c:\work\data\ ..\result\
Despues de realizar varias operaciones automaticas sobre nuestros archivos, sus nombres pueden lucir complejos y/o carecer de significado. Para cambiar todos los nombres de un grupo de GRIDs de una vez suando cierta substitucion para los nombres de los archivos se puede usar el script rename.aml.
o
&run rename.aml c:\work\data\ znc
Para ejecutar este script es necesario especificar tres cosas: 1) el folder de origen, donde estan localizados los GRIDs a renombrar 2) el patron de busqueda de archivos, fragmento de texto que utilizara el sistema para encontrar GRIDs por su nombre 3) el patron de substitucion, texto que se utilizara para substituir el patron de busqueda de archivos, este ultimo parametro puede omitirse, el resultado sera borrar el patron de busqueda del nombre del (los) GRID(s).
Estos scripts son relativamente faciles de modificar, cambiandolos es posible hacer cualquier otra operacion sobre cualquier numero de GRIDs o archivos.
Scripts disponibles: image2grid.aml, con2.aml, grid2bil.aml, setnull.aml, grid2bil.aml, clip-all.aml.
Обсудить в форуме Комментариев — 0
Ссылки по теме
Последнее обновление: September 09 2021
Дата создания: 08.08.2006
Автор(ы): Максим Дубинин
© GIS-Lab и авторы, 2002-2021. При использовании материалов сайта, ссылка на GIS-Lab и авторов обязательна. Содержание материалов - ответственность авторов. (подробнее).