Grass

Adquisicin de datos y la integracin
Recursos: resumen GRASS GIS y manual
La comprensin de los comandos de GRASS GIS

Puede ejecutar GRASS varias formas, por medio de su interfaz grfica de usuario (wxGUI,
GUI) y la interfaz de lnea de comandos (CLI):
RECOMENDADO: (. D *) para la visualizacin de comandos utilizan wxGUI,

para otros comandos escriba el nombre del comando en la consola de comandos
(ver pestaa de la parte inferior del Administrador de capas) y Enter para abrir (R *,
g *, v *...) su interfaz grfica de usuario;
ejecutar todo a travs de wxGUI: iniciar el comando usando un botn
correspondiente o seleccionarlo de un men desplegable (usar "Mdulo de
bsqueda" para encontrar el comando);
(. D *) para la visualizacin de comandos utilizan wxGUI, para otros comandos
escribir o pegar (pasta o Shift Insertar) de comandos con todas las opciones en la
consola de comandos y escriba (til cuando se est acabando el tiempo, pero usted
aprender menos); puede utilizar flechas arriba / abajo en el teclado para recuperar
el comando ejecutado anteriormente, si desea volver a ejecutarlo con diferentes
opciones;
cuando se ejecuta en Linux o Mac, pegar el comando con las opciones en la cscara
(se puede copiar / pegar varios comandos a la vez)
Tutorial basado en pantalla muestra estas opciones: Introduccin a SIG GRASS

utilizando wxGUI
Instrucciones tanto para wxGUI y el enfoque de lnea de comandos se proporcionan en las
asignaciones iniciales, pero las instrucciones sern en su mayora de lnea de comandos ms
adelante. No dude en utilizar el enfoque que ms le convenga mejor. Para una referencia
rpida, aqu estn los equivalentes wxGUI para el d * comandos.:
En la barra de herramientas del Administrador de capas SIG GRASS
d.rast : Aadir trama capa de mapa

d.vect : Aadir capa de mapa vectorial
d.shade : Aadir varias capas de mapas raster> Aadir un mapa de relieve
sombreado
En la barra de herramientas de GRASS GIS Map
D.ZOOM : opciones de zoom

d.out.file : Save visualizacin para presentar
d.barscale : Aadir elementos del mapa> Mostrar / ocultar barra de escala
d.legend : Aadir elementos del mapa> Mostrar / ocultar leyenda

d.northarrow : Aadir elementos del mapa> Mostrar / ocultar la flecha norte
El d. * Comandos que aparecen en las asignaciones indican qu capas de mapas deben ser
incluidos en el mapa de salida, d.out.file indica el mapa debe ser guardado para el
informe. Es posible que ya tienen las capas de mapas correctos en su Administrador de
capas, por lo que no tiene que cargar de nuevo. Las instrucciones de la CLI tambin se
pueden utilizar para ejecutar las tareas en Linux y Mac como scripts.
ubicacin SIG GRASS y el directorio de mapas

Antes de comenzar con GRASS GIS es importante entender que GRASS GIS utiliza el
concepto de ubicacin, y la regin mapset computacional para apoyar el anlisis y
modelado eficiente sin la necesidad de un subconjunto de datos y volver a muestrear. La
ubicacin es un espacio de proyectos que tiene el sistema de coordenadas comn. Cada
ubicacin puede tener varios directorios de mapas, uno se llama PERMANENTE para
almacenar datos de referencia, otros son fijados por los usuarios. Puede modificar los datos
slo en su mapset actual, pero se puede ver o copiar datos desde otros directorios de mapas.
Ver tambin la estructura de base de datos de GRASS .
NO MUEVA los archivos en el directorio ubicacin utilizando herramientas que no
GRASS.
flujo de trabajo sugerido es crear un nuevo directorio de mapas para cada asignacin y
hacer todos los clculos all.
Regin de clculo se establece por g.region mdulo (comando) y especifica el alcance y la
resolucin con la que se realiza la operacin de trama.
Iniciar GRASS GIS

Iniciar GRASS - hacer clic en el icono de hierba o de tipo Linux:
grass70
En Pannel de arranque establecido SIG Directorio de datos de ruta de acceso a bases de

datos, por ejemplo, en MS Windows, C: \ Users \ mi_nombre \ grassdata . Por
Ubicacin del proyecto seleccione nc_spm_08_grass7 (Carolina del Norte, del plano del
estado, metros) y para el directorio de mapas Accesible crear un nuevo directorio de
mapas (por ejemplo, llamada HW_data_acquisition ) y haga clic en Inicio GRASS .
datos de visualizacin en 2D
En primer set regin que trabajan en GUI: Ajustes> Regin> Regin Set> Establecer
regin actual de la regin llamada> selecciona swwake_30m> Ejecutar> Cerrar o,
alternativamente, la pasta en la consola de comandos:
g.region swwake_30m -p
Visualizar mapa de elevacin:

En el Administrador de capa> botn Aadir capa de trama > selecciona trama "elevacin">
Aceptar.
En Visualizacin del mapa> Zoom a la regin de medida computacional .
Asegrese de que el botn Render (esquina inferior derecha) para la representacin
automtica de mapas se comprueba sucesivamente. Tambin puede utilizar el primer botn
de mapa de visualizacin en pantalla del mapa para hacer mapas.
Alternativamente, se puede aadir capa de trama pegando este comando en la consola de
comandos:
elevacin d.rast
Pantalla de elevacin como un mapa con sombra de color:

En el Administrador de Capa> Aadir varias capas de los mapas raster > Aadir capa de
mapa de relieve sombreado > seleccione "elevation_shade" como mapa de relieve
sombreado y "elevacin", como mapa de colores.
Bajo Opcional conjunto de pestaas por ciento para iluminar al 30> OK.
lnea de comandos equivalente: pegar en la consola de comandos y escriba:
d.shade sombra de color elevation_shade = = = elevacin iluminar 30
Nota: dependiendo de la versin del csped instalado, los nombres de los d.shade mdulo
pueden ser diferentes. Abra la d.shade de dilogo para encontrar las opciones correctas.
Ahora compruebe la gama de valores de elevacin de todo el mapa. Ya sea pegando el
siguiente comando en la consola de comandos o hacer clic derecho sobre la capa raster de
elevacin y la seleccin de metadatos en el men contextual.
elevacin r.info -r
Visualizacin de vectores lineales y puntuales mapas:

En Administrador de Capas: Aadir botn de capa de mapa vectorial > seleccione
"streets_wake"> Aceptar. Ajustar los colores, anchos de lnea, smbolos haciendo clic en
las pestaas correspondientes. Aadir mapas adicionales ( "roadsmajor", "schools_wake").
Estos son los equivalentes de lnea de comandos:
d.vect streets_wake
d.vect de color rojo = ancho roadsmajor = 2
d.vect schools_wake tamao de los iconos bsica = / = 10 fill_color
crculo azul =
Guarde el mapa que se muestra:

En Visualizacin de mapa, haga clic en Guardar pantalla de archivo grfico .
myMap d.out.file
datos de visualizacin en 3D
Esto es slo una vista previa rpida: Vamos a explorar la visualizacin 3D en tercera
asignacin.
En Administrador de capas marque la capa elevation_shade, pero mantener el raster de
elevacin sucesivamente. En Mapa interruptor de visualizacin de la vista 2D a 3D. Una
vez cargado el modo 3D, aumentar la altura de visualizacin con control deslizante, rotar
disco hacia el sur. Haga clic en la ficha de datos y establecer la resolucin del modo fino a
1, para obtener la superficie mxima resolucin. Guardar la imagen haciendo clic en
Guardar pantalla de archivo grfico (se requiere el apoyo TIFF) o tomar captura de
pantalla.
vista de datos y consulta

Para familiarizarse con las herramientas de consulta y zoom, establecer la vista de nuevo a
2D. La consulta de datos significa la bsqueda de los valores en los mapas raster o atributos
en mapas vectoriales en un lugar seleccionado por el ratn. En Administrador de capas,
seleccione la capa (s) para consultar el mapa haciendo clic en l. En Visualizacin de mapa
clic el botn de consulta mapa (s) de trama / vector y haga clic en un mapa en lugares en
los que desea conocer los valores / atributos.
Visualizacin del mapa en explorar las opciones para hacer un zoom: zoom-in, zoom de
alejamiento, zoom al mapa seleccionado, hacer zoom a la regin computacional . Tambin
puede acercar a la capa de mapa especfico haciendo clic derecho sobre la capa y
seleccionando zoom al mapa (s) seleccionado.
Ahora mostrar DEM en todo el estado, los lmites del condado y las estaciones climticas:
Retirar todas las capas del mapa del Administrador de capas.
En el men de administrador de capas: Ajustes> Regin> Selecciona la regin .
En g.region de dilogo: Conjunto regin actual de la regin llamada > seleccionar
"nc_500m"> Ejecutar > Cerrar .
En Visualizacin del mapa: Zoom a la regin computacional .
En Administrador de capas: Aadir capa de trama > seleccione "elev_state_500m".
Aadir capa de mapa vectorial "precip_30ynormals", el juego de smbolos a bsico /
marcador.
Aadir capa de mapa vectorial boundary_county, para mostrar slo lmite del condado, en
la pestaa Seleccin de apagar reas. Reorganizar las capas mediante el ratn, si es
necesario.
La lnea de comandos equivalente est por debajo, a excepcin de la eliminacin de todas

las capas del mapa del Administrador de capas . Al usar interfaz grfica de usuario y su
consola de comandos , puede desactivar las capas de mapa o completamente eliminarlos en
el Administrador de capas . Cuando se utiliza la lnea de comandos del sistema y d.mon ,
puede utilizar d.erase .
g.region nc_500m -p
d.rast elev_state_500m
precip_30ynormals d.vect icono = bsica / marcador
d.vect tipo boundary_county = lmite
nc_precip_stations d.out.file
La gestin de los datos

Vamos a aprender cmo encontrar ms informacin acerca de los datos, copiar y cambiar el
nombre de los mapas.
Lista de datos raster, vector y de la regin (extensiones guardados) en la base de datos
proporcionada:
En la GUI: Archivo> Administrar mapas> Lista.
O pegue el siguiente comando en la consola de comandos:
g.list f = tipo raster, vector, regin
Lista nicos mapas vectoriales que comienzan con "s":

En la GUI: Archivo > Administrar mapas > Lista .
En la pestaa Requerido -> seleccione el tipo de mapa (s) de vectores y en la pestaa
Patrn Tipo * s de patrn de Mapa Bsqueda Nombre .
patrn de tipo vectorial g.list = = "s *"
Encontrar el trabajo del sistema y el alcance espacial de coordenadas:

En la GUI: Ajustes > Las proyecciones cartogrficas > Pantalla de proyeccin de mapa
En la GUI: Ajustes > Regin > Pantalla Regin
g.proj -p
g.region -pl
Para saber cmo copiar, renombrar y eliminar mapas:

interfaz grfica de usuario: Archivo > Administrar mapas > Copiar > selecciona / tipo en
los nombres de, al
GUI: Archivo > Administrar mapas > Renombrar > selecciona / tipo en los nombres de, al
GUI: Archivo > Administrar mapas > eliminar > seleccione el tipo y el nombre y
comprobar Fuerza de eliminacin .
Tenga en cuenta que slo se puede eliminar mapas del directorio de mapas de trabajo
actual.
g.copy trama = geology_30m, mygeology
g.rename trama = mygeology, mygeology_30m
raster g.list
Tipo g.remove = Nombre de trama = -f mygeology_30m
Encontrar informacin sobre mapas individuales capas y sus atributos:

En la GUI: haga clic derecho capa de trama "elevacin", seleccione Metadatos , haga clic
derecho en la capa de vector "streets_wake", seleccione Metadatos , Mostrar datos de
atributos .
elevacin r.info
streets_wake v.info
v.info -c streets_wake
Todos los datos (mapas) se almacenan en el directorio de la base de datos de GRASS.

Puede guardar su proyecto en un espacio de trabajo HIERBA archivo * .gxw bajo Archivo
> Espacio de trabajo > Guardar como .
Importacin y exportacin de datos raster y vectoriales

En esta seccin vamos a cubrir la importacin y exportacin de raster y vector de datos que
estn en la misma proyeccin como la ubicacin de GRASS que trabajamos. Los casos en
que los difiere de proyeccin estn cubiertos en una gua separada .
Importacin y exportacin de archivos de forma
Download del archivo con los puntos geodsicos: geod_pts_spm.zip . Guardar y
descomprimirlo en el escritorio o en un directorio en el que tenga permiso de lectura /
escritura.
Para importar el archivo en uso GUI:
File > datos vectoriales Importar > formatos de importacin comunes .

En el cuadro de dilogo, seleccione Tipo de fuente > Directorio .
Utilice Examinar botn para encontrar la ruta de acceso al directorio
"geod_pts_spm" para llenar el nombre de la fuente de campo.
A continuacin, haga clic en Importar botn.
En caso de que los datos del punto importado no se muestra automticamente, utilice
Agregar vectorial para mostrar los datos y hacer zoom a la capa vectorial.
Comparar el vector con el vector importado "streets_wake" para ver cmo se distribuyen.
Exportar el mapa como un archivo roadsmajor forma:
File > Exportar > Los formatos ms comunes

Seleccionar roadsmajor como Nombre de la entrada de mapa vectorial
Tipo C: \ mydirectory \ roadsmajor ruta como OGR nombre de origen de datos de
salida (ruta debe ser a un directorio en el que tiene permiso de escritura)
haga clic en Ejecutar
A continuacin, puede ver la roadsmajor en ArcGIS.

Importacin y exportacin de vectores cuando las secuencias de comandos
Para flujos de trabajo de secuencias de comandos de una lnea de comandos de importacin
/ exportacin puede ser utilizado.
Se requiere establecer adecuadamente ruta de los archivos externos. En este caso, "gdc.shp"
sera sustituido, por ejemplo, "C: \ path_to_my_file \ gdc.shp" en MS Windows y por
"/path_to_my_file/gdc.shp" en los sistemas de tipo Unix. Y la ruta completa, por ejemplo,
"C: \ path_to_my_file \ roadsmajor.shp", se utilizaran en lugar de "roadsmajor.shp".
v.in.ogr de entrada = salida gdc.shp = geodetic_points
v.out.ogr de entrada = salida roadsmajor = roadsmajor.shp
En MS Windows a menudo es ms fcil hacer la importacin a travs de interfaz grfica de

usuario para obtener la ruta correcta a sus datos.
En Mac, Linux y otros sistemas se puede usar el directorio de trabajo actual y / o la funcin
de autocompletar en la lnea de comandos del sistema para evitar tener que escribir el ruta
completa; por supuesto, puedes usar interfaz grfica de usuario, si lo prefiere.
archivo de importacin de trama proporcionado en formato GeoTIFF
Descargar una trama de uso del suelo en formato GeoTIFF (imagen TIFF referencia):
lc96ras_cut.tif . Guardarlo en el escritorio o en un directorio de su eleccin.
Para importar el archivo en uso GUI: Archivo > Importar datos raster > formatos de
importacin comunes . Mantenga el valor predeterminado GeoTIFF como el formato y
utilizar Navegar botn para encontrar el camino hacia descargado "lc96ras_cut.tif" archivo.
Asegrese de que el archivo est desprotegido en la Lista de capas GDAL . Haga clic en
Importar .
En MS Windows, la importacin es a menudo ms fcil de hacer, porque a travs del GUI
de conseguir el camino correcto a los datos. En general, se puede utilizar la alternativa de
lnea de comandos a continuacin.
Tenga en cuenta que se requiere para proporcionar una ruta completa de los datos de
entrada, por lo que tendra que sustituir "lc96ras_cut.tif" por algo como "C: \
path_to_my_file \ lc96ras_cut.tif "o" /path_to_my_file/lc96ras_cut.tif "dependiendo de la

ubicacin de los datos y la plataforma que est utilizando.
r.in.gdal de entrada = salida lc96ras_cut.tif = landuse96_subset
Exportacin de archivos de mapa de bits como rejilla ASCII y GeoTIFF

Ahora, los archivos de mapa de bits de exportacin como rejilla ASCII y GeoTIFF.
Utilice Archivo > Mapa exportacin trama > formatos de exportacin comunes .
No se olvide de especificar la ruta completa de los archivos de salida de lo contrario es
posible que tenga dificultades para encontrar ellos, especialmente en MS Windows.
Como alternativa, puede utilizar la lnea de comandos, pero de nuevo, no se olvide para
especificar la ruta completa al archivo de salida o configurar el directorio de trabajo a la
posicin deseada.
r.out.gdal de entrada = salida elev_ned_30m = formato elev_ned_30m.ascii
= AAIGrid
r.out.gdal de entrada = salida basin_50K = basin_50K.tif
Opcional: Aadir capa WMS

La visin y descarga de ortofotos a travs del servicio de cartografa web. Requiere buena
conexin a Internet. Si esto no funciona, informe del problema incluyendo cualquier
mensaje de error en el informe.
Regin primer cambio en la pequea zona rural y visualizar la ortofoto proporcionado.
g.region rural_1m -p
d.rast ortho_2001_t792_1m
En la barra de herramientas Administrador de capas Aadir capa de servicios web .

Pegue el enlace con el servicio en el campo de servidor y pulse conectar:
http://raster.nationalmap.gov/arcgis/services/Orthoimagery/USGS_EROS_Orth
o_1Foot/ImageServer/WMSServer?
Espere hasta cambios en la interfaz. Sea paciente, puede tardar hasta un minuto.
En la lista de capas, no debera haber "Orthoimagery_All", haga clic en eso.
Pulse el botn Aadir capa y esperar.
El cuadro de dilogo se puede cerrar despus.
Para guardar la capa:
Haga clic derecho sobre la capa de servicio Web en el Administrador de Capa - elige
Guardar capa de servicios web . En el dilogo, establecer el nombre de la capa ortho_new
y pulse Guardar capa .
modelos de datos geoespaciales

recursos:
Informacin general del SIG GRASS y manual

Recomendaciones y tutorial de cmo utilizar wxGUI de la primera asignacin
Iniciar GRASS GIS

Iniciar GRASS - hacer clic en el icono de hierba o de tipo
grass70

datos, por ejemplo, en MS Windows, C: \ Users \ mi_nombre \ grassdata .
Por Ubicacin del proyecto seleccionar nc_spm_08_grass7 (Carolina del Norte, del plano
del estado, metros) y para el directorio de mapas accesible a crear un nuevo directorio de
mapas ( por ejemplo, se llama HW_data_models).
Haga clic en Inicio GRASS .
Si usted prefiere trabajar en la interfaz grfica de usuario, que debe ser capaz de encontrar
t mismo los equivalentes GUI para las tareas siguientes.
Algunos consejos para la interfaz grfica de usuario se incluyen, pero a partir de ahora, la
mayor parte de las instrucciones que se proporcionan como lnea de comandos.
Pista para correr la mayor parte de los comandos en la interfaz grfica de usuario -. escriba
o pegue el nombre del mdulo en la consola de comandos y luego pulsa Enter para abrir el
dilogo de interfaz grfica de usuario Lea la pgina de manual para cada comando que est
utilizando por primera vez para aprender lo que est haciendo y lo que significan los
parmetros.
Volver a muestrear con una resolucin ms alta

Volver a muestrear el mapa raster dado a alta y baja resolucin (30m-> 10m, 30m-> 100 m)
y compara el muestreo si el vecino ms cercano con bilineal y el mtodo bicbico.
En primer lugar, establecer la regin de clculo medida a nuestra rea de estudio y ajustar
la resolucin a 30 metros. La regin computacional (regin para abreviar) se establece
mediante g.region mdulo. Aqu por conveniencia que utilizan named regin que define
tanto la extensin como la resolucin. Esta regin llamada est incluida en los datos
(ubicacin) que estamos utilizando, pero es posible crear nuevas regiones nombradas y
utilizarlos para marcar diferentes reas de estudio.
regin g.region = -p swwake_30m
El -p bandera para g.region se utiliza para imprimir la regin apenas fijamos.

A continuacin, presentamos una resolucin de 30 metros NED raster de elevacin.
d.rast elev_ned_30m
Para volver a muestrear a una resolucin de 10 metros, establecido por primera vez la
regin de clculo para la resolucin de 10 metros, a continuacin, volver a muestrear la
trama utilizando el mtodo del vecino ms cercano.

Sugerencia: Para abrir la r.resamp.interp en interfaz grfica de usuario, escriba o pegue el
nombre de mdulo en la consola de comandos , a continuacin, Intro para abrir el panel de
interfaz grfica de usuario, no se olvide de configurar el mtodo en virtud cercano Opcional
pestaa.
g.region res = 10 -p
r.resamp.interp elev_ned_30m cabo = = elev_ned10m_nn mtodo ms cercana
Mostrar el mapa resampled aadiendo "elev_ned10m_nn" al Administrador de capas en

caso de que no lo tiene en el Administrador de capas ya.
Como alternativa, el uso en lnea de comandos lo siguiente:
d.rast elev_ned10m_nn
El mapa de elevacin "elev_ned10m_nn" tiene el mismo aspecto que el original, por lo que
ahora, mira la superficie de elevacin resampled cmo usar el mapa de aspecto:
r.slope.aspect aspecto elev_ned10m_nn = aspect_ned10m_nn
Mostrar el mapa resampled aadiendo "aspect_ned10m_nn" al Administrador de capas o

en la lnea de comandos usando:
d.rast aspect_ned10m_nn
Guarde el mapa que se muestra y explicar lo que est pasando en el informe y en qu se

diferencia del aspecto calculada a partir del mapa de elevacin original?
Para guardar el mapa, haga clic en Visualizacin del mapa para en el botn Guardar
pantalla en el archivo grfico " o, como alternativa, el uso el siguiente comando:
d.out.file aspect_nn
Ahora, usando DEM reinterpolate bilineal y la interpolacin bicbica. Compruebe las

superficies de elevacin interpolados usando mapas de aspecto.
r.resamp.interp elev_ned_30m cabo met = = elev_ned10m_bil bilineal
r.resamp.interp elev_ned_30m cabo met = = elev_ned10m_bic bicbica
r.slope.aspect aspecto elev_ned10m_bil = aspect_ned10m_bil
aspecto elev_ned10m_bic r.slope.aspect = aspect_ned10m_bic
aspect_ned10m_bil d.rast
aspect_ned10m_bic d.rast
Guarde el mapa que se muestra y, en su informe, comparar el resultado con el resultado del
vecino ms cercano calculado previamente.
En Visualizacin del mapa de clic de botn Save visualizacin de archivos de grficos , o
utilizar el siguiente en la lnea de comandos:
aspect_bic d.out.file
Volver a muestrear con una resolucin ms baja

Volver a muestrear a una resolucin ms baja (30m -> 100 m).
En primer lugar, mostrar los mapas originales de elevacin y el uso de la tierra:
d.rast elev_ned_30m
landuse96_28m d.rast
A continuacin, cambiar la resolucin de la regin y la elevacin volver a muestrear (que es

un campo continuo) y el uso de la tierra (que tiene categoras discretas). Explicar la
seleccin del mtodo de agregacin. Se puede utilizar tambin para promedio de uso de la
tierra? Qu significa el modo?
r.resamp.stats elev_ned_30m cabo = mtodo elev_new100m_avg = Media
d.rast elev_new100m_avg
Antes del siguiente clculo, eliminar todas las capas del mapa desde el Administrador de
capas debido a que no es necesario ni verlos.
r.resamp.stats landuse96_28m cabo = mtodo landuse96_100m = Modo
Eliminar o desactivar los mapas de uso del suelo, elevacin y de aspecto.
La conversin entre tipos de datos vectoriales

Convertir los bloques censales polgonos a puntos con sus centroides (til para interpolar
una superficie de tendencia densidad de poblacin):
Tipo v.to.points census_wake2000 = centroide a cabo el uso = =
census_centr vrtice
Lmites del censo de visualizacin utilizando la GUI: Aadir vectorial "census_wake2000"

Seleccin > Tipo de accidente > lmite (desconectar los otros tipos). Guarda el mapa que se
muestra en la pantalla del mapa de clic de botn Guardar pantalla de archivo grfico .
Como alternativa, utilice los siguientes comandos para controlar la visualizacin.
Tenga en cuenta que, tanto en la lnea de comandos y la interfaz grfica de usuario debe
introducir la ruta completa al archivo que est guardando la imagen en, o que debe conocer
el directorio de trabajo actual.
d.vect icono census_centr = / = crculo fcol bsica verde tamao = 10
d.vect census_wake2000 col = Tipo rojo = lmite
vect_to_points d.out.file
Convertir a curvas de nivel (puntos tiles para calcular los contornos de DEM):
v.to.points de entrada = salida elev_ned10m_cont10m = elev_ned_contpts
type = = uso de la lnea de vrtice
Mostrar los "puntos" elev_ned_contpts del vector y zoom en la zona de muy pequeo para
ver los puntos reales.
d.vect elev_ned_contpts tamao co = = icono marrn bsico / punto = 3
Convertir a partir del vector a raster

Convertir datos del vector a raster para su uso en el anlisis basado en la trama. En primer
lugar, ajuste la regin de clculo de la resolucin 200:
g.region swwake_30m res = 200 -p
A continuacin, retire todas las capas del Administrador de capas .

Convertir puntos vectoriales "escuelas" de la trama. Como valor para la columna uso trama
atributo "CORECAPACI" para la capacidad bsica.
Para aadir leyenda en el uso GUI Aadir elementos del mapa > Mostrar leyenda / ocultar
y seleccione "schools_cap_200m".
schools_wake d.vect
v.info -c schools_wake
v.to.rast schools_wake cabo = uso schools_cap_200m = attr = attrcol tipo
CORECAPACI = punto
schools_cap_200m d.rast
d.vect streets_wake co = gris
d.legend schools_cap_200m a = 70,30,2,6
Ahora convertir lneas en el vector "calles" de la trama. Ajuste la resolucin y el uso de

atributos 30m lmite de velocidad.
v.to.rast streets_wake cabo = uso streets_speed_30m = attr = attrcol tipo
VELOCIDAD = lnea
Si no lo ha hecho, agregue eliminar todas las otras capas de mapas de administrador de

capas y aadir "streets_speed_30m" capa de trama.
Aadir leyenda para la trama "streets_speed_30m" usar interfaz grfica de usuario en
pantalla del mapa : Aadir leyenda > Establecer opciones > avanzado > Lista de discreta
nmeros de gato y el tipo de los lmites de velocidad 25,35,45,55,65; mover la leyenda con
el ratn, segn sea necesario.
Como alternativa, utilice los siguientes comandos:
d.rast streets_speed_30m
d.legend streets_speed_30m a = = 5,30,2,5 uso 25,35,45,55,65
Guarde el mapa visualizado. En Visualizacin del mapa de clic de botn Save visualizacin
de archivos de grficos , o utilizar el siguiente.
vect_to_rast d.out.file
Convertir de la trama al vector
Convertir lneas de trama de lneas vectoriales.

En primer lugar, establecer la regin y quitar las capas de mapa del Administrador de capas
. A continuacin, realice la conversin.
Explicar por qu estamos usando r.thin mdulo . Es posible que desee eliminar todas las
capas utilizadas previamente desde el Administrador de capas antes de iniciar estos nuevos
clculos.
g.region trama = -p streams_derived
d.rast streams_derived
r.thin streams_derived salida = streams_derived_t
r.to.vect streams_derived_t salida = Tipo = lnea streams_derived_t
comparar visualmente el resultado con corrientes de fotografas areas digitalizadas.

d.vect de color azul = streams_derived_t
d.vect arroyos color = red
Guarde el mapa que se muestra (en la visualizacin del mapa de clic de botn Save
visualizacin de archivos de grficos ).
streams_compare d.out.file
Convertir las reas de mapa de bits que representan las cuencas a polgonos vectoriales.
Utilice el valor de trama como el nmero categora (-v) y polgonos vector Ver llenas de
colores al azar.
En GUI: Aadir vectorial> Colores> Encender colores al azar. Es posible que desee
eliminar todas las capas utilizadas previamente desde el Administrador de capas antes de
iniciar estos nuevos clculos.
g.region trama = -p basin_50K
d.rast basin_50K
r.to.vect -sv basin_50K salida = Tipo = rea basin_50Kval
d.vect -c basin_50Kval
d.vect arroyos color = blue
Guarde el mapa que se muestra ya sea mediante GUI o usando el siguiente en caso de que
se est trabajando en la lnea de comandos.
cuencas d.out.file
visualizacin de datos y visualizacin

recursos:

Capturas de pantalla para la visualizacin 3D con wxnviz en GRASS6.4.3RC2 (similar
a GRASS7)
o La visualizacin en la hierba I: superficie (71MB)
o La visualizacin en el csped II: drapeado puntos y lneas (45 MB)
o Visualizacin en la hierba III: planos de corte (62 MB)
Iniciar GRASS GIS

estado, metros) y para el directorio de mapas Accesible crear un nuevo directorio de mapas
(por ejemplo, llamada HW_data_visualization) y haga clic en Inicio GRASS .
grass70
Visualizacin en perspectiva 3D
Establecer la regin a la trama "elevacin".
g.region trama = -p elevacin
Interactiva ver los datos de "elevacin" raster y vectoriales ( "corrientes", "usgsgages") en

vista 3D, ajuste la posicin, propiedades de la superficie de visin e iluminacin para resaltar
caractersticas.
Eliminar o desactivar cualquiera de las capas de mapa en el Administrador de capas. A
continuacin, aadir la trama "elevacin", vector Aadir "corrientes", "vector Aadir
usgsgages".
Zoom a la regin computacional y Cambie a la vista 3D.
Siga la visualizacin de video I y II para tareas similares (el enlace est arriba) y guardar las
imgenes 2-3 para su informe (guardar como TIFF o tomar captura de pantalla, si no tiene el
apoyo tiff).
La visualizacin de mltiples superficies
Visualizar mltiples superficies (tierra desnuda y la superficie con vegetacin y estructuras),
analizar su relacin mediante crossections generados por planos de corte interactivo
Se recomienda dejar de fumar GRASS antes de iniciar la tarea a continuacin, ya que vamos a
trabajar con una regin ms pequea de alta resolucin.
GRASS7 comenzar con su mapset anterior. En primer set regin para "rural_1m", entonces
interpolar superficie con vegetacin desde mltiples puntos de retorno lidar utilizando el
mdulo de v.surf.rst (explicaremos ms adelante interpolacin).
v.surf.rst elev_lidrural_mrpts elevacin = elev_lidruralmr_1m npmin = 120 =

25 segmax tensin = 40 = 1,0 lisa
En Administrador de capas raster Aadir "elev_lidruralmr_1m" y "elev_lid792_1m" (tierra

desnuda). En Visualizacin de mapa zoom a la regin computacional.
Asegrese de que tiene slo los dos mapas de elevacin activada en el Administrador de
capas. Cambie a la vista 3D y seguir la visualizacin de vdeo III (el enlace est arriba).
Guarde 2-3 imgenes interesantes para su informe, incluir al menos una crossection.
operaciones de visualizacin 2D bsica

Interfaz grfica de usuario se recomienda para las tareas abajo, ver los equivalentes GUI para
los comandos seleccionados d * arriba, la lnea de comandos es para indicar el flujo de trabajo
y la produccin.
Visualizacin subconjuntos de datos
Explorar visualmente relacin entre las reas desarrolladas y la topografa.
Uso de la tierra visualizacin categoras 1, 2 (tierra desarrollada) ms de topografa
sombreada. Ajustes > Regin > regin de Ajuste > Selecciona la regin para que coincida con
esta trama > SELECT "landuse96_28m". Aadir trama > seleccione "elevation_shade".
Aadir trama > seleccione "landuse96_28m"> seleccin > Lista de gatos para que aparezca >
1,2. zoom a la regin computacional . Haga clic derecho sobre "landuse96_28m" capa de
mapa y cambiar la opacidad para mostrar la topografa mezclado con el uso del suelo.
Guardar pantalla de archivo grfico. Es posible que desee eliminar todas las capas utilizadas
previamente desde el Administrador de capas antes de empezar a visualizar los nuevos.
g.region trama = -p landuse96_28m

d.rast elevation_shade
Los valores landuse96_28m d.rast = 1,2
landuse_elev d.out.file
Cambiar los colores de los mapas raster

Hay muchas maneras de cmo hacerlo ver r.colors manual).
En primer lugar crear su propia copia del mapa de elevacin (consulte el comando a
continuacin).
Exhbala por Agregar trama> "myelev".
g.copy trama = elevacin, myelev
myelev d.rast
Establecer la tabla de colores de uno predefinido, por ejemplo "bgyr" (azul, verde, amarillo,
rojo).
Opciones de la GUI: haga clic derecho en la capa de trama "myelev" en el Administrador de
capa> Establecer tabla de colores > Definir > Nombre de la tabla de colores > "bgyr">
Ejecutar .
Si no ve los nuevos colores - haga clic en el segundo botn de visualizacin del mapa para
volver a dibujar.
r.colors myelev color = bgyr
myelev d.rast
Cambiar la tabla de colores de forma interactiva utilizando reglas:

En el Administrador de capas manu clic Raster > Administrar colores > Administrar reglas de
color de forma interactiva .
Seleccionar mapa raster "myelev".
Haga clic en las casillas de color para cambiar el color, haga clic en Vista previa .
Si se ve bien - haga clic en Aceptar.
Para crear la tabla de colores para futuros usos mltiples escribir las reglas en un archivo de
texto plano mycolor.txt usando un editor de texto (edicin de textos, el Bloc de notas), por
ejemplo:
50 azul
70 aqua
90 verde
110 amarillo
130 de naranja
160 marrn
A continuacin, asigne la siguiente manera:

Haga clic derecho sobre "myelev" capa de mapa> conjunto de tabla de colores > Definir >
Ruta al archivo rules > Ejecutar .
Si es necesario, volver a dibujar con el segundo botn de visualizacin del mapa. Tambin
puede crear y asignar la tabla de colores por encima escribiendo o pegndolo en los valores de
entrar de forma interactiva ventana o ejecutarlo desde la lnea de comand (asegrese de que
tiene la ruta correcta, en el siguiente ejemplo mycolor.txt debe sustituirse por ejemplo, /
Usuarios /john/mycolor.txt en Mac OS X):
r.colors myelev reglas = mycolor.txt
Tenga en cuenta que se pueden combinar RGB y las definiciones de nombres, as como% y el
gato / val para crear tablas de colores personalizados complejos y almacenarlos en un archivo
de texto para uso futuro, vanse los ejemplos de r.colors pgina del manual.
Comparar el uso de intervalos iguales e histograma ecualizado tabla de colores de la
pendiente
Para agregar el uso leyenda Aadir elementos del mapa botn de visualizacin del mapa.
g.copy trama = pendiente, myslope

r.colors myslope color = bgyr
myslope d.rast
myslope d.legend
myslopecolor d.out.file
Para explicar la diferencia entre los dos mapas, se puede generar un histograma. En el mapa de
visualizacin, haga clic en Analizar mapa > Crear histograma para abrir la herramienta de
histograma y guardar los resultados en archivos de grficos. O utilizar este comando que se
sumar un histograma como una capa de Administrador de capas.
myslope d.histogram
Ahora ajuste el histograma igual tabla de colores, y guardar los nuevos mapas de pendiente y
el histograma imgenes. Describe el efecto del histograma igualaron tabla de colores.
r.colors -e color de myslope = bgyr
d.out.file myslopecolorequalized
myslope d.histogram
Modificar la leyenda, escala y rejilla

Para cambiar el tamao de la leyenda de myslope, haga clic derecho en la leyenda, seleccione
la leyenda cambiar el tamao y cambiar el tamao con el ratn. Como alternativa, puede
cambiar el tamao mediante la seleccin del indicador del modo de visualizacin del mapa,
hacer doble clic en la leyenda para iniciar el dilogo leyenda y definir las opciones>
opcionales > Colocacin > 50,90,4,7.
Los nmeros son de abajo, arriba, izquierda, derecha como porcentaje de la pantalla
coordenadas.
Aadir unidades a la leyenda: Aadir elementos del mapa > Aadir capa de texto > escriba
deg> OK.
Aadir barra de nivel de Visualizacin de mapa: Aadir elementos del mapa > Mostrar /
ocultar la flecha norte (doble clic sobre ella para cambiarla).
Nota: puede utilizar leyendas horizontal mediante el uso de nivel al = 6,10,2,30 o simplemente
estirar horizontalmente con el ratn.
d.barscale
d.legend myslope a = 50,90,4,7
Apagar "myslope", la pantalla de la trama "myelev" y cambiar la leyenda de la "myelev"

trama. Aadir rejilla para el estado de las coordenadas planas a 5000m con las garrapatas, a
1.000 m. Tambin agregar una cuadrcula de latitud / longitud en el intervalo de 2 minutos de
arco.
En el Administrador de Capa> Aadir varios datos > Aadir capa de rejilla > Tamao 5000>
Ejecutar.
Para dibujar nica frontera con las garrapatas: Cambia el tamao a 1000 y haga clic en
Desactivar > dibujo de rejilla Desactivar .> Ejecutar
Apagar el dibujo de rejilla Desactivar .
Para dibujar rejilla de latitud / longitud : conjunto Requerido > Tamao de 0:02> haga clic en
Dibujar > Dibujar cuadrcula geogrfica > Ejecutar.
alternativa de lnea de comandos:
myelev d.rast
d.legend myelev
tamao d.grid = 5000 = de color marrn
d.grid tamao -n = 1000
d.grid -g size = 0: 02 color = negro
myelevmap d.out.file
Pantalla color compuesto

En Administrador de Capas: Aadir varias capas de mapas raster > Aadir capa RGBmap .
d.rgb rojo = verde = lsat7_2002_30 lsat7_2002_20 azul = lsat7_2002_10
d.vect de color amarillo = roadsmajor
mylandsat d.out.file
Exportacin NC estaciones de precipitacin de datos a kml y post en Google

Earth
Si est ejecutando el comando en la interfaz grfica de usuario no se olvide de hacer clic en el
punto para el tipo de datos. Para ver los datos en Google Earth slo tienes que arrastrar y soltar
el archivo KML en la ventana de Google Earth.
nc_500m g.region
v.out.ogr precip_30ynormals salida = c: \ temp formato \ precipitation.kml
= tipo KML = punto
parte opcional
Aadir capa WMS
La visin y descarga de ortofotos a travs del servicio de cartografa web. Requiere buena
conexin a Internet. Si esto no funciona, informe del problema incluyendo cualquier mensaje
de error en el informe.
Regin primer cambio en la pequea zona rural y visualizar la ortofoto proporcionado.
En la barra de herramientas Administrador de capas Aadir capa de servicios web .

Pegue el enlace con el servicio en el campo de servidor y pulse conectar:
http://raster.nationalmap.gov/arcgis/services/Orthoimagery/USGS_EROS_Ortho_
1Foot/ImageServer/WMSServer?
Espere hasta cambios en la interfaz.

En la lista de capas, no debera haber "Orthoimagery_All", haga clic en eso.
Pulse Aadir botn de la capa y esperar.
El cuadro de dilogo se puede cerrar despus.
Para guardar la capa:
Haga clic derecho sobre la capa de servicio Web en el Administrador de Capa - elige Guardar
capa de servicios web . En el dilogo, establecer el nombre de la capa ortho_new y pulse

Guardar capa .
mapa Flagelo
Utilice mapswipe para comparar la nueva ortofoto con el orto proporcionada en el conjunto de
datos: Archivo > Mapa golpe fuerte . Aadir trama> seleccione "ortho_2001_t792_1m" y
"ortho_new".
Mover deslizador para comparar imgenes - Qu cambios observas? Guardar una captura de
pantalla de la ventana mapswipe que muestra la comparacin de ortos. Ver ms opciones aqu
.
Anlisis Geoespacial I:, operaciones zonales y focales,

el lgebra de mapas mundial
recursos:

Para ejecutar MapCalc expresiones, puede ejecutar el comando completo dentro de la

consola de interfaz grfica de usuario del sistema, o en caso de cualquier problema, escriba
o copie la expresin en el SIG de la trama de clculo del mapa de gramneas que se puede
iniciar desde la barra de herramientas del Administrador de capas.
Los archivos de texto con las reglas de color:
srtmneddiff_color.txt
Iniciar GRASS GIS

mapas (por ejemplo, llamada HW_map_algebra) y haga clic en Inicio GRASS .
grass70
Cambiar el directorio de trabajo: Ajustes > entorno de trabajo HIERBA > directorio de
trabajo Cambiar > Seleccionar / crear cualquier directorio o escribe cd (siglas de cambio
de directorio) en la consola de comandos de interfaz grfica de usuario y pulse Enter:
discos compactos
Descargar todos los archivos de texto con las reglas de color (vase ms arriba) en el
directorio seleccionado. Ahora puede utilizar los comandos de la tarea que requiere el
archivo de texto sin necesidad de especificar la ruta completa al archivo.
resmenes de cmputo
reas de cmputo para cada categora en dos resoluciones diferentes.
Son resultados iguales? Explica (vase tambin la Conferencia 1) . Copia y pega el informe
de la ventana de salida o guardar el informe en un archivo de texto: ventana de salida >
Guardar . Utilice fuente de ancho fijo (por ejemplo, Courier, Andale Mono en su informe
para conservar el formato).
g.region trama = landuse96_28m res
unidad landuse96_28m r.report = c,
g.region trama = landuse96_28m res
unidad landuse96_28m r.report = c,
= 12 -ap
h, p
= 30 -ap
h, p
reas de cmputo para cada categora de uso de la tierra para cada cdigo postal. Comparar
27601 27511 Raleigh con Cary. Incluir slo la parte correspondiente de la tabla en su
informe.
cdigos postales r.report, unidad landuse96_28m = h, p
Calcular estadsticas zonales mapas que representan inclinacin para cada cuenca.
Aadir leyendas utilizando Aadir elementos del mapa en la barra de herramientas
Visualizacin del mapa. Recordatorio: d.out.file significa guardar en el archivo de grficos
para su informe.
g.region trama = -p pendiente
base de la cubierta r.stats.zonal = = basin_50K mtodo de la pendiente
media de salida = = slope_avgbasin
r.colors slope_avgbasin color = gyr
slope_avgbasin d.rast
d.legend slope_avgbasin a = 90,50,5,8
d.vect arroyos color = 15: 25: 110
slope_avgbasin d.out.file
Calcular estadsticas zonales mapas que representan el uso del suelo ms comn para cada
cuenca.
Al aadir la leyenda, asegrese de definir la lista de nmeros de categora discretos .
base de la cubierta R.Mode = basin_50K = salida landuse96_28m =
landuse96_modebasin
landuse96_modebasin d.rast
corrientes d.vect
d.legend landuse96_modebasin en 40,20,2,5 = -n
landuse96_modezip d.out.file
Aplicar los operadores de vecinos

Utilice el operador barrio para calcular el uso del suelo mapa diversidad. Cun diversa es
la utilizacin del suelo en la zona de NCSU? Antes de comenzar a mostrar los nuevos
resultados, retire todas las capas del mapa aadido anteriormente del Administrador de
capas .
r.info -g landuse96_28m
r.neighbors landuse96_28m salida = mtodo lu_divers = tamao de la
diversidad = 7
lu_divers d.rast
d.legend lu_divers en 70,15,5,10 = -v
d.vect streets_wake
r.report unidad lu_divers = p
Utilice el operador barrio para suavizar el mapa de elevacin SRTM y comparar las
medidas estadsticas globales para el DEM original y alisada. Cmo sera el tamao de la
zona de influir en el resultado?
g.region trama = -p elev_srtm_30m
r.neighbors elev_srtm_30m salida = elev_srtm30m_sm5 method = tamao medio
= 5
d.rast elev_srtm_30m
d.rast elev_srtm30m_sm5
elev_srtm_30m r.univar
r.univar elev_srtm30m_sm5
Parche mltiples capas de mapa de bits en una sola trama

Azulejos parche de trama para lidar base, 6m res. DEM para Centennial Campus. Antes de
mostrar nuevos datos, eliminar las capas ms antiguas del mapa Administrador de capas .
g.region trama = el_D793_6m, el_D783_6m, el_D782_6m, el_D792_6m -p
r.patch de entrada = el_D793_6m, el_D783_6m, el_D782_6m, la produccin
el_D792_6m = elevlidD_6m
r.colors elevlidD_6m trama = elevacin
d.rast elevlidD_6m
lgebra de mapas
Ver r.mapcalc pgina del manual de sintaxis y funciones.  Si usted est recibiendo en error al ejecutar r.mapcalc en la
consola de comandos, ejecutar la versin GUI desde la barra de herramientas del
Administrador de capas.
Diferencia Normalizada calcular ndices de vegetacin (NDVI). Explicar la diferencia entre
el punto flotante y de enteros en el manejo ndvi1, ndvi2 y el resultado ndvi3 .
g.region trama = -p lsat7_2002_40

r.mapcalc "ndvi1 = (lsat7_2002_40 - lsat7_2002_30) / (+ lsat7_2002_40
lsat7_2002_30)"
r.mapcalc "ndvi2 = 1.0 * (lsat7_2002_40 - lsat7_2002_30) / (+
lsat7_2002_40 lsat7_2002_30)"
r.mapcalc "ndvi3 = flotacin (lsat7_2002_40 - lsat7_2002_30) / flotador
(lsat7_2002_40 + lsat7_2002_30)"
r.info -r ndvi1
r.info -r ndvi2
r.info -r ndvi3
r.colors de color ndvi3 = NDVI
d.rast ndvi3
d.out.file ndvi3
Explorar la diferencia entre el SRTM DSM y basada en lidar NED DEM. Calcular el mapa
de diferencias de elevacin.
g.region trama = -p elev_ned_30m
r.mapcalc "srtm_ned30_dif = elev_srtm_30m - elev_ned_30m"
Crear una tabla de colores personalizada para distinguir los valores negativos y positivos:
r.info -r srtm_ned30_dif
FP: Alcance de los datos es -142,24 a 86,19 ... ...

Asignar encargo color de la tabla srtmneddiff_color.txt .
GUI: Haga clic derecho sobre la capa> Propiedades > Establecer color de la tabla >
Colores > Ruta del archivo de reglas .
r.colors srtm_ned30_dif reglas = srtmneddiff_color.txt
Zoom a la regin computacional y apagar capas de mapas anteriores. Mostrar la capa de

mapa de diferencias:
srtm_ned30_dif d.rast
d.legend srtm_ned30_dif a = 70,15,5,10
d.vect streets_wake
srtm_ned30_dif d.out.file
elev_ned_30m r.univar
Es la SRTM sobre todo mayor o menor que elev_ned? Cuyo resultado va a utilizar para
responder a la pregunta anterior - el mapa srtm_ned30_dif o nmeros proporcionada por
r.univar ?
Trabajar con las instrucciones if
Crear un mapa de las zonas urbanas (de alta densidad y de clase baja densidad) con 0 clase
en otro lugar.
r.mapcalc "urban1_30m = si (== landuse96_28m 1,1,0) + si (==

landuse96_28m 2,2,0)"
urban1_30m d.rast
d.legend urban1_30m a = 10,30,5,8
Como alternativa con "o" operador lgico y nula en otra parte:

r.mapcalc "urban2_30m = si (landuse96_28m == == 1 || landuse96_28m 2,
landuse96_28m, nula ())"
urban2_30m d.rast
El manejo de valores nulos

Reemplazar los lagos en el mapa landclass96 con representacin ms precisa.
A partir de mapas raster lagos, lagos que dividir por 1000 porque el mapa lago tiene
grandes nmeros que habran requerido tabla de colores personalizados.
r.mapcalc "luselakes = si (IsNull (lagos), landclass96, lagos / 1000)"
luselakes d.rast
luselakes d.out.file
Crear mscara para las reas desarrolladas bajas (por ejemplo, vulnerable a las
inundaciones) con elevacin entre 60 y 100 metros y el uso de la tierra 1 o 2.
Establecer la regin, visualizar los mapas de entrada y crear una mscara.
Antes de iniciar nuevos clculos, eliminar o desactivar anterior las capas de mapa en el
pesebre de capa . Tambin puede ampliar a regin computacional en Visualizacin del
mapa una vez que establezca una nueva.
elevacin d.rast
r.mapcalc "low_elev_developed = if ((elevacin <100 && elevacin> 60) &&
(landuse96_28m == == 1 || landuse96_28m 2), 1, nula ())"
r.mask trama = low_elev_developed
Comando r.mask crea un mapa raster "mscara" en su directorio de mapas.

Eliminar la capa de "low_elev_developed" si se aadi.
Mostrar las cuencas hidrogrficas para ver el efecto de mscara:
d.rast basin_50K
d.out.file basin_masked
Desactivar la mscara, y mostrar basin_50K nuevo para mostrar que la mscara se elimin.
r.mask -r
d.rast basin_50K
El uso de coordenadas de la ventana que se mueve en el lgebra de mapas
Vuelva a colocar la seccin del DSM SRTM con elevacin DEM NED. Trata de explicar
cmo este r.mapcalc funciona expresin.
r.mapcalc "elev_combined = si (y () <224274. && x ()> 637455., elevacin,
elev_srtm_30m)"
plano inclinado
Crear un mapa que representa la trama plano inclinado (por ejemplo, falla geolgica):
r.mapcalc "tiltplane = 0,2 * (0,1 * fila () + columna ()) + 50"
r.mapcalc "tiltpl_under = si (tiltplane <elev_lid792_1m + 2, tiltplane,
nula ())"
Ver la superficie de elevacin y el plano del subsuelo en 3D. Apagar todas las capas en el
administrador de capas excepto por elev_lid792_1m y tiltpl_under . Cambio de
visualizacin en Vista 3D, ajuste la posicin de visualizacin para una vista desde el sur.
Guardar una imagen para su informe.
diversas tareas tiles adicionales opcionales Utilice el lgebra de mapas para crear subconjuntos de mapas.
Cambiar regin para el azulejo foto area 792 y la resolucin de 10 metros. Puesto que
vamos a trabajar en diferentes son, es una buena idea para eliminar todas las capas de los
mapas utilizados previamente desde el Administrador de capas .
g.region trama = ortho_2001_t792_1m res = 10 -ap
Crear un subconjunto del mapa de elevacin para esta subregin.

r.mapcalc "elevation_792_10m = elevacin"
d.rast elevation_792_10m
Alejar el zoom para ver que es un subconjunto.

Trabajar con NULL y la mscara. Ajustar la mscara y comprobar su efecto.
elevacin d.rast
d.vect streets_wake
r.mask maskcats urbanas = 55
elevacin d.rast
Retire la mscara:
r.mask -r
Calcular el% de rea para cada categora en cada cdigo postal. Qu cdigo postal tiene el
desarrollo ms alta densidad?
r.stats cdigos postales -pl, landuse96_28m
Trabajar con coordenadas relativas.

Escribe la expresin en una sola lnea y sin \
De nuevo, es una buena idea para eliminar las capas de los mapas utilizados previamente
antes de empezar a trabajar en una nueva tarea.
g.region trama = -p elev_srtm_30m
r.mapcalc "elev_srtm30m_smooth = (elev_srtm_30m [-1, -1] \
+ Elev_srtm_30m [-1,0] + elev_srtm_30m [-1,1] \
+ Elev_srtm_30m [0, -1] + elev_srtm_30m [0,0] \
+ Elev_srtm_30m [0,1] + elev_srtm_30m [1, -1] \
+ Elev_srtm_30m [1,0] + elev_srtm_30m [1,1]) / 9 ".
Asignar el mapa resultante de la misma tabla de colores que el original. Comparacin de las
estadsticas globales
r.colors elev_srtm30m_smooth trama = elev_srtm_30m
elev_srtm30m_smooth r.univar
d.rast elev_srtm_30m
d.rast elev_srtm30m_smooth
Anlisis Geoespacial II: tampones, superficies de

costos, menos ruta de coste
recursos:

Los archivos de texto con reglas de recodificacin, reglas de color y ubicacin de los sitios:
noise_cats.txt
friction_rules.txt
friction_color.txt
fire_pt.txt
lostperson.txt
Iniciar GRASS GIS

mapas (por ejemplo, llamada HW_buffers_cost) y haga clic en Inicio GRASS .
grass70
Eliminar mscara en caso de que se queda de asignacin anterior (comando de uso r.mask
-r ).
discos compactos
Descargar todos los archivos de texto con reglas de recodificacin, reglas de color y
ubicacin de los sitios (vase ms arriba) en el directorio seleccionado. Ahora puede
utilizar los comandos de la tarea que requiere el archivo de texto sin necesidad de
especificar la ruta completa al archivo.
Los tampones
Encuentra las reas desarrolladas cerca de los lagos
Selecciona la regin, crear tampones:
r.buffer lagos salida = distancias lakes_buff = 60120240500
d.rast lakes_buff
Categoras de la Lista en la tierra mapa de uso para identificar los nmeros de gato para las
reas desarrolladas. A continuacin, ejecute r.mapcalc :
landuse96_28m r.category
r.mapcalc "developed_lake = si (landuse96_28m == == 1 || landuse96_28m 2,
lakes_buff, nula ())"
r.support developed_lake trama = lakes_buff
Mostrar resultados:
En la GUI: se apagan mapa anterior, utilice Agregar trama, Aadir iconos vectoriales.
developed_lake d.rast
d.vect de color gris = streets_wake
lagos d.rast
d.legend uso developed_lake -c = 2,3,4,5 a = 5,25,2,5
mylakesbuffmap d.out.file
Halla el rea total dentro de los tampones y el rea desarrollado en ha:

unidades lakes_buff r.report -n = h
r.report unidades developed_lake -N = h
Encuentra las reas desarrolladas afectadas por el ruido de las carreteras

Regin de establecer y crear memorias intermedias a lo largo de las carreteras principales.
Intersect reas de mapa de uso del suelo con tampones carretera desarrollados.
Transferir las etiquetas de categora y calcular el rea afectada.
Ejecutar r.category de interfaz grfica de usuario y busque el archivo txt, o ejecutar

r.category comando (asumiendo noise_cats.txt est en el directorio de trabajo).
r.buffer roadsmajor salida = roads_buffers distancias = 250,500,2500
r.mapcalc "ruido = si (landuse96_28m == == 1 || landuse96_28m 2,
roads_buffers, nula ())"
r.colors color de ruido = ryg
normas de ruido r.category = = separador noise_cats.txt:
r.report -n unidades de ruido = p, h
Mostrar los resultados. En GUI apagar mapas anteriores.

ruido d.rast
d.vect streets_wake
ruido d.legend -c a = 5,25,2,5
Encontrar escuelas afectadas por altos niveles de ruido

Convertir las escuelas a raster usando atributo CORECAPACI (capacidad de la escuela).
Utilice el lgebra de mapas para superponer con amortiguadores de impacto de ruido y
calcular el nmero de estudiantes expuestos a ruido (ver resultado de r.univar ).
v.to.rast salida schools_wake = uso schoolscap_10m = attr = attrcolumn
tipo CORECAPACI = punto
r.mapcalc "schools_noise = si (int (schoolscap_10m) && roads_buffers ==
2, int (schoolscap_10m), null ())"
schools_noise r.univar
r.to.vect schools_noise salida = tipo schools_noise = punto
d.vect schools_wake tamao de los iconos bsica = / = 10 fill_color
crculo azul =
d.vect schools_noise size = 14 = cian fill_color icono = bsica / crculo
de color negro =
mynoisemap d.out.file
superficies de costes
Calcular la superficie costo acumulado a un lugar del accidente determinado,
basndose en los lmites de velocidad
Utilice 5mi lmite de velocidad / hr para las reas fuera de la carretera (valores nulos).
Para la leyenda establece los valores lmite de velocidad discretos en la leyenda de dilogo
GUI:
pestaa Requerido > streets_speed, subconjunto > Lista de nmeros de categora discretos
v.info -c streets_wake
v.to.rast salida streets_wake = uso streets_speedtmp = attrcolumn attr =
Tipo = lnea VELOCIDAD
r.mapcalc "streets_speed = si (isnull (streets_speedtmp), 5,
streets_speedtmp)"
streets_speed r.info
r.colors de color streets_speed = gyr
streets_speed d.rast
d.legend streets_speed a = = 5,40,2,5 uso 5,25,35,45,65
myspeedmap d.out.file
Asignar el tiempo de viaje para cruzar una celda de cuadrcula de 30 metros en horas.
Tenga en cuenta que el costo de GRASS GIS se define como el tiempo de viaje por cada
clula .
Calcular la superficie costo acumulado hasta el punto dado. Dnde est la 0.018641
constante viene? Puede modificar la expresin para obtener el tiempo en minutos.
r.mapcalc "streets_travtime = 0.018641 / streets_speed"
r.cost salida streets_travtime -k = streets_cost start_coordinates =
634886,224328
Importar la ubicacin del punto de accidentes contenidas en fire_pt.txt para que pueda
visualizarlo en el mapa.
Ejecutar v.in.ascii de interfaz grfica de usuario, utilice el explorador de archivos para
proporcionar la ruta completa al archivo o introducir valores interactiva pegando las
coordenadas en el proporcionado espacio, o ejecutar el comando suponiendo que el archivo
est en el directorio de trabajo.
v.in.ascii de entrada = salida fire_pt.txt = separador de espacio =
fire_pt
Calcular iscronas y mostrar el mapa de superficie costo acumulado.

r.contour streets_cost paso de salida = 0.04 = streets_cost_04
r.colors streets_cost color = elevacin
d.rast streets_cost
d.vect de color rojo tamao fire_pt = = icono bsica / marcador = 20
d.vect streets_cost_04
d.legend streets_cost a = 5,50,2,5
mycostmap d.out.file
Encuentra costo (tiempo de viaje) de firestations seleccionados

En primer lugar hacer su propia copia del mapa firestations y la lista de atributos.
A continuacin, consultar la superficie costo acumulado en el lugar firestations.
El tiempo de viaje en horas se almacenar en el CVLAG columna de atributos.
g.copy vector = firestations, myfirestations
v.info myfirestations -c
v.what.rast myfirestations trama = = columna streets_cost CVLAG
Aadir myfirestations al Administrador de capas y haga clic derecho sobre el botn Mostrar
datos de atributos.
Haga clic en CVLAG que ordenar y encontrar el costo ms bajo (menor tiempo)> 0
(firestations con coste 0 estn fuera de la regin).
Usted debe ver una tabla como la siguiente .
Exportacin firestations con tiempo de trnsito de menos de 0,1 hr:
v.out.ascii myfirestations separador de espacio = precisin = 3 columnas

= ID, ubicacin, donde CVLAG = "CVLAG <0,1 y CVLAG> 0"
Usted debe obtener algo como lo siguiente:

635775.565
635940.262
637386.831
633178.155
228121.693
225912.796
222569.152
221353.037
19
20
21
52
19 4021 Distrito Dr. 0,076

0 5001 Western Blvd 0,037
0 1,721 0,071 Trailwood Dr.
27 6000 Holly Springs Rd 0.060
Encuentra la ruta de menor costo para las dos estaciones ms cercanas:

r.drain -n entrada = salida streets_cost = route_20Westernb
start_coordinates = 635940.3,225912.8
r.drain -n entrada = salida streets_cost = route_52Hollyb
start_coordinates = 633178.2,221353.0
r.colors route_20Westernb de color gris =
r.colors route_52Hollyb de color gris =
Mostrar los resultados en interfaz grfica de usuario.

Puede que ya tenga algunos de estos mapas en rbol de capas interfaz grfica de usuario,
por lo que slo encenderlos y fuera y re-orden, segn sea necesario:
d.rast streets_cost
d.vect de color rojo tamao fire_pt = = icono bsica / marcador = 20
d.rast route_20Westernb
d.rast route_52Hollyb
mylcpmap d.out.file
Imprimir la longitud de la ruta en clulas (se multiplican por 30 metros para llegar aprox.
M).
Debe tener el tiempo en horas ya partir del mapa costo. A qu velocidad media [km / h]
necesita el recorrido del carro para llegar all en prevista hora? Es el tiempo, la velocidad y
la distancia realista?
r.describe route_20Westernb
r.describe route_52Hollyb
Si desea visualizar la superficie streets_cost en 3D como se muestra en la conferencia, es

necesario multiplicar la superficie en 1000 utilizando r.mapcalc porque los nmeros en
horas son muy bajos.
A continuacin, apague todas las capas excepto por el streets_cost_1000 y cambiar la
pantalla a 3D.
Calcular mapa de accesibilidad para ayudar en la bsqueda de la persona perdida
Crear un mapa de friccin basado en la ocupacin del suelo.
Clases de cobertura de la primera pantalla de la tierra:
landclass96 r.category
Recodificar el mapa de uso del suelo al mapa de friccin usando las reglas en
friction_rules.txt .
r.recode landclass96 cabo = reglas de friccin = friction_rules.txt
Aadir las calles al mapa de friccin - que se pierden por el mapa de uso del suelo.
r.mapcalc "friction2 = si (streets_speed> 6, 0,1, friccin)"
r.colors friction2 reglas = friction_color.txt
d.rast friction2
d.legend friction2 a = 5,50,2,5
Recordatorio: Tiempo al tiempo, se debe eliminar o al menos desmarque utilizado

previamente capas de mapa en el Administrador de capas , por lo que las capas no se rindi
cuando no los necesita ms.
Calcule el mapa de costos y generar iscronas:
r.walk -k = elevacin de la friccin elev_ned_30m = salida friction2 =
walkcost start_coordinates = 635576,216485 lambda = 0,5 max_cost = 10000
r.contour salida walkcost = paso walkcost = 1000
Para mostrar la entrada, las coordenadas de importacin del punto en que la persona
perdida fue visto por ltima vez dan en el archivo lostperson.txt .
v.in.ascii de entrada = salida lostperson.txt = separador lostperson =
coma
Mostrar los resultados:

lagos d.rast
d.rast walkcost
d.vect streets_wake
d.vect de color marrn walkcost =
d.vect de color rojo donde walkcost = = = nivel de 6000 width = 2
tamao d.vect lostperson = 20 color = negro = fill_color icono negro =
bsica / marcador
d.legend walkcost a = 5,50,2,5
mywalkcostmap d.out.file
superficies de costes a las entidades de lnea

Calcular la distancia ms corta a una carretera principal.
Establecer regin, convertir vector de hoja de ruta para la trama:
v.to.rast salida roadsmajor = uso roadsmajor = tipo val = lnea
Calcule la distancia en el mapa y la superficie de coste a las carreteras:

r.mapcalc
r.cost de
r.mapcalc
r.mapcalc
"area_one = 1"
entrada = salida area_one = dist_toroad start_rast = roadsmajor
"dist_meters = dist_toroad * (ewres () + () nsres) / 2."
"dist_class = int (dist_meters / 500)"
dist_class d.rast
d.vect roadsmajor
dist_roadsmaj d.out.file
r.cost -k de entrada = salida streets_travtime = cdist_toroadk start_rast
= roadsmajor
r.colors cdist_toroadk color = bgyr
cdist_toroadk d.rast
d.vect de color negro = streets_wake
d.vect de color rojo = roadsmajor
cdist_toroadk d.legend a = 5,50,2,5
cost_roadsmaj d.out.file
interpolacin espacial y la aproximacin I: mtodos

recursos:

Iniciar GRASS GIS

grass70

mapas ( por ejemplo, se llama HW_interpolation_1).
Calcular polgonos de Voronoi

Mostrar los polgonos con centroides.
Encontrar el nombre de la columna donde se almacena z y convertir los polgonos de la
trama.
Calcule aspecto a evaluar la geometra de la superficie.
v.voronoi elev_lid792_randpts salida = elev_vor
d.vect elev_vor size = 1 Tipo -c = rea, centroide
v.info -c elev_vor
v.to.rast elev_vor salida = elev_vor_1m attrcolumn = valor de uso = attr
r.colors elev_vor_1m color = elevacin
r.slope.aspect aspecto elev_vor_1m = asp_vor_1m
Mostrar el mapa de elevacin resultante y mapa de aspecto.

d.rast elev_vor_1m
d.rast asp_vor_1m
d.vect elev_lid792_randpts size = 1 color = red
elev_voronoi d.out.file
Opcionalmente, se puede ver en perspectiva 3D (apagar todas las capas excepto por
elev_vor_1m y cambiar a la vista 3D).
Interpolacin utilizando IDW

Establecer regin y resolucin, encontrar un nombre de columna en la que se almacena z.
Interpolar DEM usando IDW, comprobar el resultado usando aspecto.
v.info -c elev_lid792_randpts
v.surf.idw elev_lid792_randpts salida de la columna = elev_idw_1m = valor
r.colors elev_idw_1m color = elevacin
r.slope.aspect aspecto elev_idw_1m = asp_idw_1m
d.rast elev_idw_1m
d.rast asp_idw_1m
d.out.file elev_idw
experimento de diseo que aclara el impacto de los parmetros IDW en la superficie, se

centran en el impacto de:
exponente por ejemplo, potencia = 0,5, 1, 5 (2 es el valor predeterminado)

nmero de puntos vecinos por ejemplo, npoint = 1, 5, 20, 60 (12 es el valor
predeterminado)
Incluir imgenes seleccionadas (por ejemplo, sombreado o aspecto) y correspondientes

estadsticas (por ejemplo, media, mnimo, mximo desde r.univar , histograma) que ponen
de relieve las diferencias en las superficies resultantes en su informe.
Compruebe la superficie interpolada con parmetros por defecto utilizando la vista 3D.
No se olvide de apagar todo a excepcin de las elevaciones interpolados y establecer la
resolucin fina a 1.
Puede utilizar el color constante de la superficie para poner de relieve su estructura.
Guardar una imagen para su informe.
Calcular DEM de contornos

Calcular DEM de contornos utilizando la interpolacin lineal entre isolneas:
v.to.rast salida elev_lid792_cont1m = = el_lid792_cont1m attrcolumn uso

level = attr
r.surf.contour el_lid792_cont1m salida = el_rcont
r.colors de color el_rcont = elevacin
Compruebe el resultado usando un mapa 2D o aspecto el_rcont en 3D.

En 3D fijada vista desde SE y la luz del NW para revelar sutil geometra.
aspecto el_rcont r.slope.aspect = asp_rcont
d.rast el_rcont
d.rast asp_rcont
d.vect elev_lid792_cont1m col = blanco
d.out.file asp_rcont
Opcional: crear el modelo TIN

Convertir z-valor almacenado como "valor" atributo a la coordenada z.
Compute estao y visualizarla como datos vectoriales 3D:
Mantener solamente, "elev_idw_1m" y "elev_rand_tin" "level90" encendido.
Para "elev_rand_tin"> Propiedades en Seleccin Tab apagar todo a excepcin de lmite.
cambiar la vista de 2D a 3D y quitar elev_idw_1m.
v.to.3d elev_lid792_randpts salida de la columna = elev_lid792_randpts3d
= valor
v.delaunay salida elev_lid792_randpts3d = elev_rand_tin
r.mapcalc "level90 = 90"
interpolacin espacial y la aproximacin II: estras

recursos:

Los archivos de texto con las reglas de color:
deviations_color.txt
precip_color.txt
Iniciar GRASS GIS

grass70

discos compactos
Descargar todos los archivos de texto con las reglas de color (vase ms arriba) en el
archivo de texto sin necesidad de especificar la ruta completa al archivo.
Interpolar raster de elevacin de puntos con estras y con diferentes

tensiones
Calcular aspecto simultneamente con la interpolacin y evaluar el impacto de la tensin
mediante el uso de la tensin = 40 (por defecto), la tensin y la tensin = 10 = 160.
g.region rural_1m res = 1 -p
v.surf.rst elev_lid792_randpts elevacin = elev_rstdef_1m zcolumn = valor
de relacin = asp_rstdef_1m segmax = 30 = 140 npmin
v.surf.rst elev_lid792_randpts elevacin = = aspecto elev_rstt10_1m
asp_rstt10_1m zcolumn = tensin de valor = 10 = 30 segmax npmin = 140
v.surf.rst elev_lid792_randpts elevacin = = aspecto elev_rstt160_1m
asp_rstt160_1m zcolumn = tensin de valor = 160 segmax = 30 = 140 npmin
Comparacin de las superficies de elevacin interpolados usando mapas de aspecto.

Cambie la tabla de colores de aspecto gris de aspecto.
Guardar imgenes para su informe.
r.colors asp_rstdef_1m color = aspecto
r.colors asp_rstt10_1m color = aspecto
r.colors asp_rstt160_1m color = aspecto
d.rast elev_rstdef_1m
d.rast asp_rstdef_1m
d.out.file asp_rst_t40
d.rast asp_rstt10_1m
d.rast asp_rstt160_1m
O usar las vistas 3D de elev_rstdef_1m, elev_rstt10_1m, elev_rstt160_1m, asegrese de

apagar los raster de aspecto y guardar las 3 imgenes para su informe.
raster de elevacin de clculo y desviaciones mapa vectorial de puntos

Para diferentes valores de suavizar las estadsticas de desviacin comparar para suavizar
0,1 y 10.
Bsqueda medio de la raz cuadrada rmse desviacin.
v.surf.rst elev_lid792_randpts elevacin = elev_rstdef_1mb zcolumn =
valor lisa = 0,1 desviaciones = elev_rstdef_devi segmax = 30 = 140 npmin
v.build elev_rstdef_devi
v.surf.rst elev_lid792_randpts elevacin = elev_rstsm10_1mb zcolumn =

valor lisa = 10 = desviaciones elev_rstsm10_devi segmax = 30 = 140 npmin
v.build elev_rstsm10_devi
v.info -c elev_rstdef_devi
v.univar columna elev_rstdef_devi = tipo Flt1 = punto
elev_rstdef_1mb r.info
v.info -c elev_rstsm10_devi
v.univar columna elev_rstsm10_devi = tipo Flt1 = punto
elev_rstsm10_1mb r.info
Calcular y mostrar las desviaciones mapas utilizando misma tabla de colores. Es necesario
utilizar tabla de colores personalizados para ver bien los resultados.
Tenga en cuenta que estamos aqu interpolar las desviaciones, no las elevaciones dadas.
v.surf.rst elev_rstdef_devi elevacin = elev_rstdef_devi zcolumn = Flt1
segmax = 30 = 140 npmin
v.surf.rst elev_rstsm10_devi elevacin = elev_rstsm10_devi zcolumn = Flt1
segmax = 30 = 140 npmin
Aplicar la tabla de colores descargado deviations_color.txt a la trama de desviacin.

De forma opcional, para ver los resultados en 3D utilizar "elev_rstdef_1mb" para la
elevacin (desconectar todo lo dems) y la cada de los mapas como desviaciones de color.
r.colors elev_rstsm10_devi reglas = deviations_color.txt
r.colors elev_rstdef_devi trama = elev_rstsm10_devi
d.rast elev_rstdef_devi
d.rast elev_rstsm10_devi
d.legend elev_rstsm10_devi a = 2,50,2,6
d.out.file elev_rstsm10_devi
Calcular error predictiva de interpolacin

Calcular el error de prediccin de interpolacin para cada punto utilizando la validacin
cruzada (sin salida de trama, slo los puntos con errores pred.).
v.surf.rst -c elev_lid792_randpts zcolumn = valor cvdev = elev_rstdef_cv
npmin = 120 = 35 segmax
v.build elev_rstdef_cv
columna elev_rstdef_cv v.univar = tipo Flt1 = punto
Calcular un mapa de la trama de los errores de prediccin y de identificar los lugares donde
el muestreo es insuficiente.
De forma opcional, para ver el resultado en 3D utilizar "elev_rstdef_1mb" para la elevacin
(desconectar todo lo dems) y cubra el mapa de validacin cruzada "elev_rstdef_cv" como
el color.
elevacin elev_rstdef_cv v.surf.rst = elev_rstdef_cv zcolumn = Flt1
r.colors elev_rstdef_cv trama = elev_rstsm10_devi
d.rast elev_rstdef_cv
d.vect elev_rstdef_cv size = 2
d.legend elev_rstdef_cv a = 2,50,2,6
elev_rstdef_cv d.out.file
Interpolar precipitacin con influencia de la topografa

Establecer la regin 3D (lea la pgina del manual de g.region ).
Hemos establecido tbres a alto valor - tenemos un solo nivel, porque no estamos calculando
la trama 3D (vase la presentacin para ms detalles).
g.region trama = -p elev_state_500m
g.region t = 2000 b = 0 tbres = 2000 = 500 res3 -P3
Calcular la precipitacin mapa raster y sin influencia de la elevacin (con segmax = 700
segmentacin no se lleva a cabo lo que la funcin de interpolacin se calcula utilizando
todos los puntos a la vez).
Vamos a usar la mscara durante la interpolacin.
r.mask trama = ncmask_500m
v.info -c precip_30ynormals
precip_30ynormals v.surf.rst elevacin = precip_annual_500m zcolumn =
segmax anual = 700
Utilice la tabla de colores descargado el precip_color.txt .

Zoom a la regin de clculo cuando se muestra el resultado.
r.colors precip_annual_500m reglas = precip_color.txt
precip_annual_500m d.rast
d.legend precip_annual_500m a = 2,30,2,5 rango = 970,2400
Calcular la precipitacin mapa raster con la elevacin.

No es tanto la produccin voxel 3D y salida de trama 2D - queremos que la salida de trama
2D (cross_output).
De forma opcional para ver los resultados en 3D, apague todo a excepcin de
elev_state_500m y precip_30ynormals_3d,
cambiar a 3D, ajuste (escribir) de altura espectador en 300000, z-exag a las 6, cosa fina = 1,
el uso precip_anntopo_500m para el color, ajustar el tamao del icono de los puntos esfera, 5000.
visualizar el resultado y guardar la imagen de la memoria.
v.info -c precip_30ynormals_3d
v.vol.rst precip_30ynormals_3d cross_input = elev_state_500m cross_output
= precip_anntopo_500m maskmap = elev_state_500m wcolumn = anual zScale =
90 = 700 segmax
r.colors precip_anntopo_500m trama = precip_annual_500m
precip_anntopo_500m d.rast
Tratar de explicar cmo fue la elevacin utilizado para la interpolacin de trama

precipitacin.
Despus de que haya terminado, retire la mscara.
r.mask -r
Geomorphometry I: modelado del terreno

Recursos: resumen GRASS GIS y manual , herramientas libLAS para las conversiones de
datos lidar
Descargar los datos lidar Raleigh 2013 como archivos LAS y ortofoto:
azulejo LAS 0793_016 en Metros NC State Plane

geotif ortofoto (mosaico en 1 pie resolucin)
Iniciar GRASS GIS

grass70

mapas ( por ejemplo, se llama HW_terrain_modeling).
discos compactos
Descargar todos los archivos (vase ms arriba) en el directorio seleccionado. Ahora puede
utilizar los comandos de la tarea que requiere el archivo sin la necesidad de especificar la
ruta completa al archivo.
Analizar tierra desnuda y mltiples propiedades de los datos de retorno

lidar por hurgar en la basura
Importar los puntos utilizando v.in.lidar . Podemos especificar qu clase y que devuelven
(primero, segundo, apellido) que desee importar.
Podemos ver la clasificacin, ya sea en los metadatos distribuido con los datos LIDAR o
puede ser representada con lasinfo herramienta (en caso lasinfo comando no est
disponible, Saltarlo.):
tile_0793_016_spm.las lasinfo
Clase 2 representa los puntos de tierra desnuda.

Ahora nos importan los puntos de tierra desnuda y puntos de retorno primeros separado:
v.in.lidar -t entrada = tile_0793_016_spm.las salida = elev_lid793016_be
class_filter = 2
v.in.lidar -t entrada = tile_0793_016_spm.las salida = =

elev_lid793016_1r return_filter primera
Nota: si tiene algn problema con los archivos, consultar las instrucciones anteriores o usar
interfaz grfica de usuario de dilogo v.in.lidar para navegar para obtener la ruta de los
archivos.
Establecer la regin desde el archivo de puntos importada con una resolucin de 1 metro.
Calcular mapas raster (con r.in.lidar ) que representan el nmero de puntos por celda de la
cuadrcula 1 m.
Comparar las densidades de punto de tierra desnuda, primera vuelta.
g.region vect = elev_lid793016_1r res = 1 -p
r.in.lidar de entrada = tile_0793_016_spm.las salida = lid_be_binn1m
class_filter = 2 method = n
r.in.lidar de entrada = tile_0793_016_spm.las salida = = lid_1r_binn1m
return_filter primer mtodo = n
r.colors lid_be_binn1m, color lid_1r_binn1m = -e bcyr
d.rast lid_be_binn1m
d.legend lid_be_binn1m a = 2,50,2,9
unidad lid_be_binn1m r.report = p, c
lid_be_binn1m r.univar
d.rast lid_1r_binn1m
d.legend -d -s lid_1r_binn1m a = 2,50,8,12
unidad lid_1r_binn1m r.report = p
lid_1r_binn1m r.univar
Calcular un mapa de trama que representa la elevacin media de cada celda de 1 m.

Tendr agujeros.
r.in.lidar de entrada = tile_0793_016_spm.las salida = lid_be_binmean1m
class_filter = 2 = media mtodo
r.in.lidar de entrada = tile_0793_016_spm.las salida = = lid_1r_binmean1m
return_filter primer mtodo = media
r.colors lid_be_binmean1m color = elevacin
r.colors lid_1r_binmean1m color = elevacin
mapa r.info = lid_1r_binmean1m
d.rast lid_be_binmean1m
d.legend lid_be_binmean1m a = 2,40,2,5
d.rast lid_1r_binmean1m
d.legend lid_1r_binmean1m a = 2,40,2,5
mylidmrmean1m d.out.file
Calcular un mapa de trama que representa la elevacin media de cada celda de 2 m.

El resultado es casi lo suficientemente bueno para el 1 de vuelta, pero todava hay muchos
agujeros de tierra desnuda.
g.region res = 2 -ap
r.in.lidar de entrada = tile_0793_016_spm.las salida = lid_be_binmean2m
class_filter = 2 = media mtodo
r.in.lidar de entrada = tile_0793_016_spm.las salida = = lid_1r_binmean2m
return_filter primer mtodo = media
r.colors lid_be_binmean2m color = elevacin
r.colors lid_1r_binmean2m color = elevacin

d.rast lid_be_binmean2m
d.legend lid_be_binmean2m a = 2,40,2,5
d.rast lid_1r_binmean2m
d.legend lid_1r_binmean2m a = 2,40,2,5
mylidmrmean2m d.out.file
Calcular gama de valores de elevacin en cada celda de la cuadrcula:

r.in.lidar de entrada = tile_0793_016_spm.las salida = lid_be_binrange2m
class_filter = 2 = mtodo de rango
r.in.lidar de entrada = tile_0793_016_spm.las salida = lid_1r_binrange2m
return_filter = = primer mtodo gama
d.rast lid_be_binrange2m
d.legend lid_be_binrange2m a = 2,40,2,5
d.rast lid_1r_binrange2m
d.legend lid_1r_binrange2m a = 2,40,2,5
Importacin y ortofoto pantalla, se apagan todas las capas excepto por ortofoto.
r.in.gdal ortho_0793_016_ncspm.tif salida = ortho_2013_0793
d.rast ortho_2013_0793
Identificar las caractersticas que se asocian con valores grandes de rango.

Buscar slo los altos valores de rango. Lo que la cobertura del suelo se asocia con una gran
gama de mltiples datos de retorno?
Los valores lid_1r_binrange2m d.rast = 10.-200.
mylid_1rrange2m d.out.file
Ahora sabemos cmo densa son los datos y cul es el rango dentro de la clula.
Si necesitamos un DEM resolucin 1m o DSM para nuestra aplicacin de este anlisis nos
dice que necesitamos para interpolar desde la nube de puntos. Qu pasos empezar cuando
procesamiento de los datos del punto de enturbiamiento no est familiarizado con?
Interpolar DEM y DSM

Para interpolar DEM y DSM utilizamos parmetros por defecto, excepto para el nmero de
puntos utilizados para la segmentacin y la interpolacin
(segmax y npmin y una mayor tensin para mltiples datos de retorno).
Puede configurar dmin = 1 para hacer la interpolacin funcione ms rpido (consulte el
manual).
Sea paciente, puede tardar unos minutos en ejecutarse dependiendo de la potencia de los
ordenadores.
g.region res = 1 -ap
v.surf.rst elev_lid793016_be elevacin = elev_lid793016_be_1m npmin = 120
= 25 segmax dmin = 1
v.surf.rst elev_lid793016_1r elevacin = elev_lid793016_1r_1m npmin = 120
= 25 segmax tensin = 100 = 0,5 lisa dmin = 1
r.colors elev_lid793016_be_1m color = elevacin
r.colors elev_lid793016_1r_1m color = elevacin
d.rast elev_lid793016_be_1m
d.rast elev_lid793016_1r_1m
Ocultar leyenda y apagar todas las capas del mapa a excepcin de los ltimos 2
interpolados queridos.
Utilice la vista 3D con planos de corte para comparar la tierra desnuda y la superficie del
terreno.
Asegrese de que la resolucin fina se establece en 1 para todas las superficies.
Asignar a cada color constante de la superficie, aadir avin constante de referencia.
Shade la crossection usando el color de la opcin superficie inferior.
Si no recuerda esto, ver el vdeo de captura de pantalla para la vista 3D.
Guardar la imagen para su informe.
Geomorphometry II: anlisis espacial y temporal del

terreno
Recursos: resumen GRASS GIS y manual , GRASSbook .
Descargar tablas directorio de mapas y de color:
Descargar series de tiempo Nagshead y copiarlo en su ubicacin nc_spm_08_grass7

(debe ser el mismo directorio nivel que PERMANENTE). No deje que su programa
de descompresin crear directorio de nivel adicional con el mismo nombre! Si no
est seguro acerca de la estructura de base de datos SIG GRASS ledo en el manual
.
stddev_color.txt
regrslope_color.txt
Iniciar GRASS GIS

grass70

mapas ( por ejemplo, se llama HW_terrain_analysis).
discos compactos
Descargar todos los archivos de texto (vase ms arriba) en el directorio seleccionado.

Ahora puede utilizar los comandos de la tarea que requiere el archivo de texto sin necesidad
de especificar la ruta completa al archivo.
Calcular los parmetros topogrficos bsicos: pendiente y orientacin

elevacin r.slope.aspect = pendiente de elevacin = = aspecto myslope
myaspect
Pantalla resultante mapas con la leyenda usando el GUI.

myslope d.rast
d.legend myslope a = 2,40,2,6
myslope d.out.file
d.rast myaspect
d.legend myaspect a = 2,40,2,6
myaspect d.out.file
. Mostrar impacto de valores enteros en metros sobre rasante y el patrn aspecto

Calcular DEM nmero entero y deriva su pendiente y la orientacin.
Utilice interfaz grfica de usuario para visualizar el histograma: en Visualizacin del mapa
> Analizar > Crear histograma :
r.mapcalc "elev_int = int (elevacin)"
r.slope.aspect aspecto elev_int = pendiente aspect_int_10m =
slope_int_10m
slope_int_10m d.histogram
myslope d.histogram
aspect_int_10m d.histogram
d.histogram myaspect
myslope d.rast
Zoom en la zona suroeste de la regin actual (rea relativamente plana cerca grande de
intercambio). Puede explicar la diferencia en los mapas de pendientes derivado de nmero
entero (m resolucin vertical) y de punto flotante (resolucin mm vertical) DEM?
slope_int_10m d.rast
slope_int d.out.file
Calcular la pendiente a lo largo de la carretera

Calcular la pendiente de una ruta de autobs. valores de primer extracto de capa de mapa de
pendientes:
res g.region trama = elevacin = 20
v.to.rast entrada = tipo busroute11 = salida de lnea = = dir uso
busroute11
r.mapcalc "route_slope = si (busroute11, myslope)"
A continuacin, calcular la pendiente en la direccin de la ruta:
v.to.rast entrada = tipo busroute11 = salida de lnea = = dir uso

busroute11_dir
r.mapcalc "route_slope_dir = abs (atan (tan (myslope) * cos (myaspect busroute11_dir)))"
r.colors mapa = route_slope, color route_slope_dir = pendiente
Pantalla, zoom si es necesario, y calculan las estadsticas univariantes y comentar sobre la

diferencia de dos resultados.
route_slope r.univar
route_slope_dir r.univar
d.shade sombra de color elevation_shade = = = 20 myslope iluminar
r.contour de entrada = salida elevacin = contornos paso = 2
contornos d.vect
route_slope d.rast
route_slope d.out.file
d.rast route_slope_dir
route_slope_dir d.out.file
curvaturas
Calcular la pendiente, aspecto y curvaturas simultneamente con la interpolacin.
Usted puede hacer los ejemplos a continuacin los datos de tierra desnuda solamente
(primer ejemplo),
ejemplo de retorno mltiple es opcional (si usted es curioso cmo se diferencia de BE).
v.surf.rst elev_lid792_bepts elevacin = = aspecto elev_lid_1m asp_lid_1m
pcurvature = pc_lid_1m tcurvature = tc_lid_1m npmin = 120 = 25 segmax
elevacin = elev_lidmr_1m aspecto v.surf.rst elev_lidrural_mrpts =
asp_lidmr_1m pcurvature = pc_lidmr_1m tcurvature = tc_lidmr_1m npmin =
120 = 25 segmax tensin = 300 = 1 lisa.
Mostrar los resultados como imgenes en 2D o en 3D.

Para la vista 3D, apague todo a excepcin de superficie de elevacin que desea ver.
Coloque los mapas raster parmetros topogrficos sobre DEM como el color.
d.rast elev_lid_1m
pc_lid_1m d.rast
profcurvature1 d.out.file
d.rast elev_lidmr_1m
pc_lidmr_1m d.rast
Los mapas reflejan la curvatura patrn de encuesta en lugar de las caractersticas

topogrficas.
As que reducen la tensin y aumentar el suavizado.
Puede utilizar varias pantallas para comparar los resultados de lado a lado.
elevacin = elev_lidt15_1m aspecto v.surf.rst elev_lid792_bepts =
asp_lidt15_1m pcurvature = pc_lidt15_1m tcurvature = tc_lidt15_1m npmin =
120 = 25 segmax tensin = 15 = 1 lisa.
elevacin = elev_lidmrt15_1m aspecto v.surf.rst elev_lidrural_mrpts =

asp_lidmrt15_1m pcurvature = pc_lidmrt15_1m tcurvature = tc_lidmrt15_1m
npmin = 120 = 25 segmax tensin = 15 = 1 lisa.
d.rast elev_lidt15_1m
pc_lidt15_1m d.rast
profcurvature2 d.out.file
d.rast tc_lidt15_1m
d.rast elev_lidmrt15_1m
pc_lidmrt15_1m d.rast
Accidentes geogrficos
Extraer formas de relieve en diferentes niveles de detalle ajustando el tamao de la ventana
en movimiento.
Establecer subregin rural a 1 m de resolucin, calcular formas de relieve utilizando 9m y
45m barrio: leer el manual para obtener ms informacin. Explicar los tipos de accidentes
geogrficos y el papel del tamao de la vecindad.
salida elev_lid792_1m r.param.scale = tamao feature9c_1m = 9 = mtodo de
funcin
salida elev_lid792_1m r.param.scale = tamao feature45c_1m = 45 = mtodo
de funcin
Mostrar con la leyenda, guardar las imgenes de la memoria.

Opcionalmente mostrar los mapas de caractersticas cubri sobre elev_lid792_1m como el
color.
d.rast feature9c_1m
d.legend feature9c_1m a = 2,20,2,6
d.rast feature45c_1m
d.legend feature45c_1m a = 2,20,2,6
d.vect de color marrn elev_lid792_cont1m =
anlisis de series temporales de la trama

Para este ejercicio vamos a utilizar NagsHead_series Mapset que ha descargado.
Tienes que hacer primero el directorio de mapas accesibles.
En interfaz grfica de usuario: men Ajustes -> entorno de trabajo HIERBA -> Mapset
acceso o usando un comando:
operacin g.mapsets = aadir directorio de mapas = NagsHead_series
Si no ve el directorio de mapas, asegrese de que ha descargado y descomprimido

correctamente.
Ejecutar el anlisis de series y explicar los resultados: el mapa que est ncleo y
envolvente? Qu formas de relieve tienen una desviacin estndar alto y qu significa?
g.region trama = -p NH_2008_1m
d.rast NH_2008_1m
r.series NH_1999_1m, NH_2001_1m, NH_2004_1m, NH_2005_1m,

NH_2008_1m cabo = mtodo NH_9908_min = mnimo
NH_2008_1m cabo = mtodo NH_9908_max = mximo
NH_2008_1m cabo = mtodo NH_9908_mintime = min_raster
NH_2008_1m cabo = mtodo NH_9908_maxtime = max_raster
NH_2008_1m cabo = mtodo NH_9908_range = rango
NH_2008_1m cabo = mtodo NH_9908_avg = Media
NH_2008_1m cabo = mtodo NH_9908_stddev = STDDEV
NH_2007_1m,
NH_2007_1m,
NH_2007_1m,
NH_2007_1m,
NH_2007_1m,
NH_2007_1m,
NH_2007_1m,
r.colors NH_9908_stddev gobierna = stddev_color.txt

d.rast NH_9908_stddev
d.rast NH_9908_range
Usar plano de corte en la vista 3D para mostrar el ncleo y envolvente.

Aadir plano de elevacin constante en -1m de referencia, fijado en algn lugar zexag 3-5
(el valor predeterminado es demasiado alto).
Asignar superficies de color constante, utilice la superficie superior o inferior para el color
crossection .
Cuando se utiliza la parte superior de color, bajar la fuente de luz para hacer la parte
superior oscura superficie y resaltar la crossection.
Anlisis del potencial de la energa solar y la cuenca

visual
recursos:

GRASSbook
Los archivos de texto con ubicaciones de los sitios:
viewshed_points.txt
Iniciar GRASS GIS

mapas (por ejemplo, llamada HW_viewshed_solar) y haga clic en Inicio GRASS .
grass70
discos compactos
Descargar todos los archivos de texto (vase ms arriba) en el directorio seleccionado.

Ahora puede utilizar los comandos de la tarea que requiere el archivo de texto sin necesidad
de especificar la ruta completa al archivo.
anlisis de cuenca visual

Calcular cuenca visual de una nueva torre de 32 pisos ubicado en el centro.
g.region trama = elev_ned_30m -ap
v.in.ascii -z de entrada = salida viewshed_points.txt = puntos de vista
separador =, z = 3
salida elev_ned_30m r.viewshed = tower_165_los coordina = 642212,224767
observer_elevation = 165 MAX_DISTANCE = 10000
Mostrar el resultado en relieve sombreado:

d.his hue = tower_165_los intensidad = elevation_shade
d.vect streets_wake
d.vect tamao de puntos de vista = 10 color = icono naranja = bsica /
marcador
d.barscale a = 3,6
d.legend tower_165_los a = 50,90,2,6
towerview_angle d.out.file
Averiguar cul es el uso de la tierra dentro de la vista usando el lgebra de mapas:

r.mapcalc "tower_los_lu30m = si (tower_165_los, landclass96, nula ())"
r.colors tower_los_lu30m trama = landclass96
r.category trama tower_los_lu30m = landclass96
r.report -n tower_los_lu30m unidad = p, h
Mostrar slo las siguientes capas y guardar el resultado de:

tower_los_lu30m d.rast
d.vect streets_wake
d.legend tower_los_lu30m a = 25,50,1,3
towerviewlu d.out.file
Tambin podemos hacer la visibilidad desde antigua sede RedHat:

salida elev_ned_30m r.viewshed = redhat_25_los coordina = 638898,224528
observer_elevation = 25 = 10000 MAX_DISTANCE
Mostrar slo las siguientes capas y guardar el resultado de:

redhat_25_los d.rast
d.vect streets_wake
d.vect tamao de puntos de vista = 10 = col icono rojo = bsica /
marcador
d.out.file RHview
Utilice mapalgebra para calcular el uso del suelo en la vista, asignar los colores de uso de
suelo de uso de la tierra mapa visibles y etiquetas de las categoras utilizando r.colors y
r.category (vase el ejemplo anterior) y utilizar r.report -n para comparar la composicin
de tamaos y uso de la tierra dentro de cuenca visual de la torre de glbulos rojos y la sede
de humedad relativa .
anlisis de la radiacin solar

Establecer la regin y poner el edificio proyectado para el DEM, vamos a utilizar este
nuevo DEM para los anlisis.
Eliminar todas las capas y hacer zoom a la regin.
r.mapcalc "elevfacility_1m = si (isnull (instalacin), elev_lid792_1m,
138)."
r.colors elevfacility_1m color = elevacin
d.rast elevfacility_1m
Preparar mapas de entrada (pendiente y orientacin):

r.slope.aspect aspecto elevfacility_1m = pendiente aspect_elevfac_1m =
slope_elevfac_1m
ngulo de incidencia y sombras

Calcular la posicin del sol el 22 de diciembre a las 2:25 pm, EST (sin salida de mapa de
espera):
r.sunmask -s elevfacility_1m ao = 2,001 mes = 12 das = 22 hora = 14
minutos = 25 s = 0 = -5 zona horaria
. Calcular los ngulos de incidencia incluyendo sombras

Asignamos histograma igual tabla de colores - Puede explicar por qu? (Pista: tratar la
misma tabla de colores sin -e). Cul es el valor en el techo? Cmo se relaciona al da /
hora?
r.sun elevfacility_1m aspecto = pendiente aspect_elevfac_1m =
slope_elevfac_1m incidout = incid_elevfac_1m da = 356 = 14,416667 tiempo
incid_elevfac_1m r.info
r.colors -e incid_elevfac_1m co = bcyr
d.rast incid_elevfac_1m
d.legend incid_elevfac_1m a = 25,50,1,3
incid_angle d.out.file
Se extrae la zona de sombra proyectada de 14,4 horas y calcular y extraer rea de sombra
para el 7,5 hr:
r.mapcalc "shadow_1m = si (isnull (incid_elevfac_1m), 1, nula ())"
r.colors shadow_1m de color gris =
d.rast elevfacility_1m
shadow_1m d.rast
r.sun elevfacility_1m aspecto = pendiente aspect_elevfac_1m =
slope_elevfac_1m incidout = incid_elevfac7_1m da = 356 = 7,50 momento
r.mapcalc "shadow7_1m = si (isnull (incid_elevfac7_1m), 1, nula ())"
shadow7_1m d.rast
sombras d.out.file
Radiacin solar
Calcular mundial (haz + difusa + refl) de radiacin para todo el da durante el verano y el
solsticio de invierno.
Visualizar los mapas de radiacin y tambin el tiempo de insolacin mapas insol_time
utilizando mismos comandos.
Opcionalmente mostrar los mapas de radiacin se cubre ms elevfacility_1m elevacin en
la vista 3D para ver relacin entre el terreno la geometra y la radiacin solar.
r.sun elevfacility_1m aspecto = pendiente aspect_elevfac_1m = da
slope_elevfac_1m = 356 = g356 glob_rad insol_time = its356
r.colors g356 Co = -e gyr
r.sun elevfacility_1m aspecto = pendiente aspect_elevfac_1m = da
slope_elevfac_1m = 172 = glob_rad G172 insol_time = its172
r.colors G172 Co = -e gyr
d.rast g356
g356 d.legend a = 25,50,1,3
solar_winter d.out.file
d.rast G172
d.legend G172 = 25,50,1,3 al rango = 8800,8867
solar_summer d.out.file
Calcular la radiacin e insolacin tiempo solar directa para una regin ms grande.
Trate de encontrar buenas tablas de colores (a medida, hist. Igualado, para ver el patrn).
r.slope.aspect aspecto elev_ned_30m = pendiente asp_ned_30m = slp_ned_30m
r.sun aspecto elev_ned_30m = pendiente asp_ned_30m = slp_ned_30m
linke_value = 2,5 albedo_value = 0,2 beam_rad = B356 diff_rad = D356
refl_rad = R356 insol_time = it356 da = 356
Zoom a la nueva regin de clculo y mostrar las siguientes capas:

B356 d.rast
B356 d.legend a = 2,30,2,6
d.out.file beamrad356
d.rast it356
d.legend it356 a = 2,30,2,6
insol_time d.out.file
modelado de la erosin
recursos:

Breves tericos fondo, ecuaciones, unidades
Los valores C factoriales
Factor K
Los archivos de texto con reglas de recodificacin y reglas de color:
Descripcin del uso del suelo categora ( lu_labels.txt )

Tabla C-factor de recodificacin ( cfac_rules.txt )
Tabla C-factor de colores ( cfac_color.txt )
Tabla de colores prdida de suelo ( soilloss_color.txt )
La erosin / deposicin tabla de colores ( erdep_color.txt )
Clases de erosin / deposicin ( erdep_class.txt )
Etiquetas de la erosin / deposicin de clase ( erdep_label.txt )
Iniciar GRASS GIS

mapas (por ejemplo, llamada HW_erosion) y haga clic en Inicio GRASS .
grass70
discos compactos
Descargar todos los archivos de texto con reglas de recodificacin, reglas de color y
ubicacin de los sitios (vase ms arriba) en el directorio seleccionado. Ahora puede
utilizar los comandos de la tarea que requiere el archivo de texto sin necesidad de
especificar la ruta completa al archivo.
Calcular desprendimiento de suelo utilizando USLE3D

Compute potencial topogrfico (factor LS)
Comparar el impacto de los exponentes de la funcin de encendido en el patrn de erosin.
Qu ecuacin representa las condiciones cuando el flujo de agua tiene un mayor impacto y
qu modelos ms fuerte influencia de la pendiente?
g.region trama = -p elev_lid792_1m
r.slope.aspect elev_lid792_1m pendiente = = aspecto slope_1m aspect_1m
r.flow elev_lid792_1m flowaccumulation = flowacc_1m
r.mapcalc "lsfac3d_1m = 1.2 * pow (flowacc_1m * 1./22.1,0.2) * pow (sin
(slope_1m) /0.09,1.3)"
r.mapcalc "lsfac3d_s1_1m = 1,5 * pow (flowacc_1m * 1./22.1,0.5) * pow
(sin (slope_1m) /0.09,1.0)"
r.colors lsfac3d_s1_1m color = -e gyr
r.colors lsfac3d_1m trama = lsfac3d_s1_1m
Mostrar capas y guardar los resultados:

lsfac3d_1m d.rast
d.vect elev_lid792_cont1m
d.legend lsfac3d_1m a = 2,50,2,6
lsfac_s13 d.out.file
d.rast lsfac3d_s1_1m
d.vect elev_lid792_cont1m
d.legend lsfac3d_s1_1m a = 2,50,2,6
lsfac_s10 d.out.file
Calcular el desprendimiento del suelo para la cobertura del suelo variables en el

espacio
Establecer la regin de la zona rural y landcover_1m recodificar a cFactor utilizando el
r.recode mdulo. Asignar tabla de colores especiales y visualizacin de resultado.
r.recode salida landcover_1m = reglas cfactor_1m = cfac_rules.txt
r.colors mapa = reglas cfactor_1m = cfac_color.txt
d.rast cfactor_1m
Calcular la ecuacin USLE3D usando lgebra de mapas, cfactorbare_1m es el mismo que

cfactor_1m, cfactorgrow_1m tiene el uso del suelo recodificada para condiciones de
verano.
Especificar las unidades de los mapas raster utilizando el r.support mdulo. Comparacin
de las tasas de erosin y de distribucin para el invierno (desnudo) y las condiciones de
verano.
r.mapcalc "soillossbare_1m = 270. * * soils_Kfactor lsfac3d_1m *
cfactorbare_1m"
r.mapcalc "soillossgrow_1m = 270. * * soils_Kfactor lsfac3d_1m *
cfactorgrow_1m"
r.colors soillossbare_1m reglas = soilloss_color.txt
r.colors soillossgrow_1m trama = soillossbare_1m
mapa r.support = unidades soillossbare_1m = "ton / (acre.year)"
mapa r.support = unidades soillossgrow_1m = "ton / (acre.year)"
Eliminar las capas anteriores y mostrar los nuevos:

d.rast soillossbare_1m
d.legend soillossbare_1m a = 2,50,2,6
soillossbare d.out.file
d.rast soillossgrow_1m
d.legend soillossgrow_1m a = 2,50,2,6
soillosgrow d.out.file
soillossbare_1m r.univar
soillossgrow_1m r.univar
Calcular nuevo DEM con erosin labrada en

r.mapcalc "elev_erodedb_1m = elev_lid792_1m- (soillossbare_1m / 100.)"
Elev_erodedb_1m pantalla en 3D y cuelgue de soillossbare_1m como el color.

Para verla en 3D apagar todo excepto elev_erodedb_1m.
Coloque soillossbare_1m como el color y no se olvide de establecer la resolucin fina a 1.
El uso de iluminacin de SW, z-exag al menos 2,0
mapas de la red de erosin / deposicin de computacin (utilizando USPED)
Calcular mapas netas de erosin / deposicin como la divergencia del flujo de sedimentos
(USPED).
En primer flujo de sedimentos de cmputo y sus componentes en direcciones x, y:
r.mapcalc "sedflow_1m = 270. * * soils_Kfactor cfactorgrow_1m *
flowacc_1m * sen (slope_1m)"
r.colors sedflow_1m trama = soillossbare_1m
d.rast sedflow_1m
r.mapcalc "QSX = sedflow_1m * cos (aspect_1m)"
r.mapcalc "QSY = sedflow_1m * sen (aspect_1m)"
Calcular el cambio del flujo de sedimentos en las direcciones X e Y y luego cambiar en la

direccin del flujo, utilizando divergencia.
r.slope.aspect QSX dx = qsx_dx
r.slope.aspect dy = QSY qsy_dy
r.mapcalc "erdep = qsx_dx + qsy_dy"
r.info -r erdep
r.colors erdep reglas = erdep_color.txt
erdep d.rast
d.legend erdep a = 2,50,2,6 rango = -100100
Elev_lid792_1m pantalla en 3D y la cada sobre erdep como el color (apagar todas las
capas excepto por elev_lid792_1m).
resumen estadstico de cmputo
Utilice r.recode para clasificar la erosin / deposicin y r.category aadir etiquetas (estable,
alta erosin, etc.) en las categoras individuales:
r.recode erdep salida = reglas erdep_class = erdep_class.txt
normas erdep_class r.category = erdep_label.txt sep =:
r.report erdep_class unidad = p, h, una
Ejemplo de salida:
[...]
| # | Descripcin | % | hectreas | acres |
| -4 | erosin severa. . . | 0,19 | 0.101300 | 0.25031 |
| -3 | alta erosin. . . . | 1.34 | 0.701600 | 1.73365 |
| -2 | la erosin moderada. . | 3.89 | 2.042600 | 5.04726 |
| -1 | bajo la erosin. . . . | 19.74 | 10.366000 | 25.61438 |
| 0 | estable. . . . . . . | 61.32 | 32.192000 | 79.54643 |
| 1 | baja deposicin. . . | 8,40 | 4.407600 | 10.89118 |
| 2 | deposicin moderada | 2.49 | 1.307500 | 3.23083 |
| 3 | alta deposicin. . | 1.29 | 0.676900 | 1.67262 |
| 4 | deposicin grave. | 0,24 | 0.126100 | 0.31159 |
| * | Ningn dato. . . . . . . | 1.10 | 0.578400 | 1.42922 |
| ------------------------------------------------- - |
| Total | 100,00 | 52.500000 | 129.7275 |
Opina sobre ventajas, desventajas y riesgos de la clasificacin de los datos de erosin /

deposicin en categoras.
Calcular las estadsticas univariantes:

erdep r.univar
Busque la lnea con suma:

suma: -2374.473814
Las unidades son toneladas / (acre.year), pero nuestras clulas son de 1 m 2 . Por lo tanto
tenemos que dividir por 4046 (1 acre = 4046 m 2 ) para obtener tonelada anual total
transportado desde la cuenca. Puede utilizar la terminal de Python do la divisin.
flujo de agua superficial y la erosin de simulacin

A continuacin, se calcula la escorrenta superficial y el transporte de sedimentos utilizando
el modelo de muestreo ruta existente.
recursos:

libro HIERBA
Para la animacin en GRASS GIS 7, consulte las instrucciones de vdeo o instrucciones

para el manejo de datos espacio-temporal y la visualizacin en el taller GRASS GIS
Iniciar GRASS GIS

grass70

la profundidad de flujo superficial y la descarga

g.region rural_1m res = 2 -p
Calcular derivadas parciales para definir el gradiente de la superficie de elevacin:

v.surf.rst entrada -d = elev_lid792_bepts elevacin = pendiente
elev_lid792_2m = = aspecto dx_2m tensin dy_2m = 15 = 1,5 lisa npmin =
150
Tenga en cuenta que las derivadas parciales pueden ser tambin calcularon utilizando
r.slope.aspect .
Calcular patrn espacial de la profundidad de flujo superficial y la descarga mediante la
ejecucin del modelo de SIMWE implementado en el r.sim.water mdulo. Vamos a utilizar
el exceso de lluvias uniforme, la infiltracin y la ocupacin del suelo.
elevacin r.sim.water -t = elev_lid792_2m dx = dy = dx_2m dy_2m
rain_value = 50 = 0 infil_value man_value = 0.05 = profundidad de
descarga wdp_2m = nwalkers disch_2m = 100000 niterations = 120 = 20
output_step
Mostrar los resultados, aadir la leyenda de trama individual y lo que cree que es
apropiado.
d.rast wdp_2m.120
d.rast disch_2m.120
Opcionalmente, la animacin de la serie de tiempo: Archivo> herramienta de animacin

escorrenta pico con direccin de flujo predefinido
Calcular la escorrenta pico con la direccin del flujo predefinida a lo largo de la corriente
inluding alcantarilla en virtud de la carretera.
En primer lugar, la direccin de clculo (aspecto) de los flujos dados:
v.to.rast flujos de salida = = dir uso streams_dir_2m
Calcule la corriente dx y dy utilizando la direccin y pendiente igual a 2 grados:

r.mapcalc "dxstr_2m = tan (2.) * cos (streams_dir_2m)"
r.mapcalc "dystr_2m = tan (2.) * sin (streams_dir_2m)"
Calcular DEM combinado y dx y dy corriente:

r.mapcalc "dxdemstr_2m = si (isnull (dxstr_2m), dx_2m, dxstr_2m)"
r.mapcalc "dydemstr_2m = si (isnull (dystr_2m), dy_2m, dystr_2m)"
Ejecutar el modelo:
elevacin r.sim.water -t = elev_lid792_2m dx = dy = dxdemstr_2m
dydemstr_2m rain_value = 50 = 0 infil_value man_value = 0.05 =
profundidad de descarga wdpstr_2m = nwalkers dischstr_2m = 100000
niterations = 120 = 20 output_step
d.rast dischstr_2m.120
Utilice la herramienta de animacin para animar la descarga.

Escorrenta de la cubierta vegetal variables en el espacio y el exceso de precipitaciones
Calcular el escurrimiento de la cobertura del suelo espacialmente variable representada por

Manning espacialmente variables y la infiltracin.
Vamos a utilizar la variable coeficiente de Manning se define mediante la reclasificacin de
clase de cobertura. Este es el contenido de un archivo llamado land_to_mannings.txt con
reglas de reclasificacin para el r.recode mdulo.
1: 1: 0,9: 0,9
2: 2: 0,5: 0,5
3: 3: 0,01: 0,01
4: 4: 0,03: 0,03
5: 5: 0,01: 0,01
6: 6: 0,05: 0,05
7: 7: 0,1: 0,1
8: 8: 0,1: 0,1
9: 9: 0,9: 0,9
10: 10: 0,03: 0,03
11: 11: 0,5: 0,5
r.recode de entrada = salida landcover_1m = reglas mancover_2m =
land_to_mannings.txt
Del mismo modo, vamos a crear la trama con tasas de exceso de lluvia espacialmente
variables en funcin de las clases de cobertura del suelo. Utilizamos archivo llamado
land_to_rain.txt para especificar las tarifas de las clases individuales.
1: 1: 50: 50
2: 2: 5: 5
3: 3: 40: 40
4: 4: 35: 35
5: 5: 50: 50
6: 6: 40: 40
7: 7: 25: 25
8: 8: 15: 15
9: 9: 50.: 50.
10: 10: 40: 40
11: 11: 10: 10
Una vez ms, se utiliza el archivo de normas relativas al r.recode mdulo.

r.recode de entrada = salida landcover_1m = reglas raincover_2m =
land_to_rain.txt
Ejecutar el modelo:
r.sim.water
dydemstr_2m
de descarga
niterations
elevacin -t = elev_lid792_2m dx = dy = dxdemstr_2m

lluvia = raincover_2m infil_value = 0 = hombre = profundidad
mancover_2m wdpstrcov_2m = nwalkers distrcov_2m = 100000
= 120 = 20 output_step
Mostrar los resultados estticos.

d.rast wdpstrcov_2m.120
d.legend wdpstrcov_2m.120
d.rast distrcov_2m.120
d.legend distrcov_2m.120
Animacin de la serie de tiempo con la funcin de animacin.
velocidad de flujo de sedimentos, erosin y deposicin

Clculo de la tasa de flujo de sedimentos y la erosin / deposicin neta utilizando parte del
transporte de sedimentos del modelo implementado en el SIMWE la r.sim.sediment
mdulo. Para hacer los clculos ms rpido, fijado regin justo a la parte alta de la cuenca:
g.region s = s + 290
Calcular la capacidad de transporte de entrada y mapas de coeficiente de desprendimiento:

r.mapcalc "tranin = 0,001"
r.mapcalc "DETIN = 0,001"
Calcule entrada esfuerzo crtico de corte:

r.mapcalc "tauin = 0,01"
Ejecutar el modelo utilizando la ltima profundidad de la ejecucin anterior:

g.copy rast = wdp_2m.120, wdp_2m
elevacin r.sim.sediment = elev_lid792_2m dx = dy = dx_2m dy_2m
water_depth = wdp_2m detachment_coeff = DETIN transport_coeff = tranin
shear_stress = tauin man_value = 0.05 = 1000000 nwalkers niterations =
120 = transport_capacity tcapacity tlimit_erosion_deposition = erdepmax
sediment_flux = sedflow erosion_deposition = erdepsimwe
Mostrar estos resultados despus de unos segundos:

tcapacity d.rast
d.rast erdepmax
Mostrar los resultados finales:

d.rast sedflow
d.legend sedflow
erdepsimwe d.rast
d.legend erdepsimwe
El flujo de enrutamiento y el anlisis de las cuencas

hidrogrficas
recursos:

archivo de texto con el punto de salida de las coordenadas:
coordenadas del punto ( outlet_point.txt )
Iniciar GRASS GIS

mapas (por ejemplo, llamada HW_hydrology) y haga clic en Inicio GRASS .
grass70
discos compactos
Descargar todos los archivos de texto con la ubicacin de los sitios (vase ms arriba) en el
archivo de texto sin necesidad de especificar la ruta completa al archivo. Iniciar GRASS
Calcular la direccin del flujo, la acumulacin de flujo y subcuencas

Calcular la direccin del flujo, la acumulacin de flujo y subcuencas con aprox. tamao de
10000 clulas de 30m NED.
r.watershed -s elev_ned_30m umbral de acumulacin = 10000 = = accum_10K
drenaje draindir_10K cuenca = basin_10K
Extraer las corrientes ms detallados de la trama acumulacin de flujo:

r.mapcalc "streams_der_30m = if (abs (accum_10K)> 100, 1, nula ())"
Convertir al formato y mostrar los resultados a lo largo del vector con las corrientes
oficiales del Condado de Wake (rojo):
r.to.vect -s salida basin_10K = tipo basin_10K = rea
r.thin salida streams_der_30m = streams_der_30m_t
r.to.vect -s streams_der_30m_t cabo = Tipo = lnea streams_der_30m
Generar un mapa con sombra y la pantalla: De qu manera las corrientes procedentes

comparan con el mapa oficial de corriente?
Cmo se puede modificar la expresin MapCalc para hacer orgenes de transmisin
encajan mejor con el mapa oficial de corriente?
d.his hue = basin_10K intensidad = elevation_shade iluminar = 40
d.vect basin_10K type = lmite
lagos d.rast
d.vect de color azul streams_der_30m =
d.vect arroyos color = red
mystreams d.out.file
Crear un mapa de depresiones DEM

Llenado depresin es a menudo necesario para ciertos algoritmos de enrutamiento flujo
pero puede alterar los datos de elevacin de manera significativa. Averiguar qu tan extensa
las depresiones estn en nuestra DEM.
Tenga en cuenta que r.watershed no necesita ninguna depresin llenado gracias a su
algoritmo subyacente que utiliza ruta de coste mnimo para superar depresiones.
r.fill.dir de entrada = salida elevacin = = elev_fill1 direccin reas
dir1 = unres1
entrada r.fill.dir = salida elev_fill1 = = elev_fill2 direccin reas
Dir2 = unres2
entrada r.fill.dir = elev_fill2 las zonas de salida = elev_fill3
direccin Dir3 = = unres3
r.mapcalc "depr_bin = if ((elevacin-elev_fill3) <0, 1, nula ())"
Retirar todas las capas utilizadas previamente desde el Administrador de capas y mostrar
los nuevos resultados, comparar las depresiones derivadas con lagos reales:
d.rast elevacin
d.vect roadsmajor
d.rast depr_bin
d.vect lagos type = rea fill_color = aqua
d.out.file depresiones
Derivar rea que contribuye para una salida dada

Establecer la regin de la zona de estudio alta resolucin y el zoom a la misma:
r.watershed -como elev_lid792_1m umbral = 5000 = acumulacin accum_5K
drenaje = draindir_5K cuenca = basin_5K
Eliminar las capas anteriores. Mostrar ms corrientes extrados imagen area:

d.rast accum_5K valores = 1-1000000
Identificar toma de corriente de la corriente extrada.

Crear un mapa del vector con el punto al este del norte = 638.845,52 = 220.085,26
(descargar el archivo de texto) que tiene accum_5K = 224.510.
v.in.ascii de entrada = salida outlet_point.txt = = separador outletA30
espacio
d.vect de color amarillo = outletA30
Delinear el rea que contribuye asociado con esta salida y convertirlo a formato vectorial:
r.water.outlet de entrada = salida draindir_5K = basin_A30 coordina =

638845.52,220085.26
r.to.vect -s basin_A30 salida = Tipo = rea basin_A30
Pantalla lmite de la cuenca, junto con contornos.

Calcular el rea de la cuenca (categora 1)
Tipo d.vect basin_A30 = color de lmite de ancho = verde = 2
r.contour salida elev_lid792_1m = paso elev_lid792_cont_1m = 1 = 104
MinLevel
d.vect de color blanco elev_lid792_cont_1m =
d.out.file watershedA30
unidad r.report basin_A30 = h, una
Evaluar y mitigar el impacto de la carretera en flowrouting

Comparacin de las corrientes extradas (Acum> 800) con los datos oficiales de la
secuencia:
d.rast
d.rast
d.vect
d.vect
ortho_2001_t792_1m
accum_5K valores = 800-1000000
de color rojo = streets_wake
arroyos color = verde
Tallar un canal predefinido determinado por los datos de la secuencia en

DEM
raster r.carve = vector elev_lid792_1m = width = 2 corrientes de
profundidad = 0,8 salida = elev_lidcarved_1m
r.colors elev_lidcarved_1m trama = elev_lid792_1m
d.rast elev_lidcarved_1m
Extraer corrientes del DEM tallada y comparar con el mapa oficial corrientes. Cul es la
diferencia entre accum_5K, accumc_5K1m y arroyos? Explicar las ventajas y desventajas
de la talla.
umbral -como r.watershed elev_lidcarved_1m = 50000 = acumulacin
accumc_5K1m
d.rast accumc_5K1m valores = 1500-10000000
corrientes d.vect
d.out.file streamcarved
Calcular patrn de flujo refinado usando D-inf

Comparacin de las lneas de flujo de ladera y ladera abajo: En qu tipo de accidente
geogrfico (cresta, valle) que convergen?
Eliminar o desactivar todas las capas anteriores.
g.region trama = -p elev_lid792_1m
r.flow elev_lid792_1m lnea de flujo = lneas de flujo flowlength =
flowlg_1m flowaccumulation = flowacc_1m
r.flow -u elev_lid792_1m flowlength = flowlgup_1m flowaccumulation =

flowaccup_1m
Mostrar mapas junto con contornos para ver relacin al terreno:

d.rast flowacc_1m
d.vect elev_lid792_cont_1m color = red
d.rast flowaccup_1m
d.vect elev_lid792_cont_1m color = red
de lneas de flujo d.vect
d.out.file lneas de flujo
Comparar el resultado mltiple direccin de flujo con D8 y Dinf

Comparar el resultado de accum_5K y accum_mfd_5K1m. Por qu son diferentes, lo que
es ms preciso y por qu?
r.watershed -a elev_lid792_1m umbral = 50000 = acumulacin accum_mfd_5K1m
d.rast accum_mfd_5K1m
d.vect corrientes
d.out.file stream_mfd
ndice de humedad de cmputo

r.topidx elev_lid792_1m salida = wetness_1m
d.rast wetness_1m
Elev_lid792_1m pantalla en 3D y cuelgue de wetness_1m como el color.

Nota: apague todas las capas excepto por elev_lid792_1m.
Crear un mapa de rea inundada

Crear un mapa de la zona inundada para un determinado nivel de agua y punto de siembra:
r.lake elev_lid792_1m water_level = 113,5 lago = flood1 coordenadas =
638728,220278
d.rast elev_lid792_1m
d.rast flood1
d.out.file floodedarea
Aumentar el nivel de agua para 113.7m y 114.0m y crear mapas de reas inundadas en
estos dos niveles.

Grass

Cargado por

Copyright:

Formatos disponibles

Grass

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Grass

Cargado por

Copyright:

Formatos disponibles

Adquisicin de datos y la integracin

Recursos: resumen GRASS GIS y manual

La comprensin de los comandos de GRASS GIS

RECOMENDADO: (. D *) para la visualizacin de comandos utilizan wxGUI,

Tutorial basado en pantalla muestra estas opciones: Introduccin a SIG GRASS

d.rast : Aadir trama capa de mapa

En la barra de herramientas de GRASS GIS Map

D.ZOOM : opciones de zoom

d.legend : Aadir elementos del mapa> Mostrar / ocultar leyenda

ubicacin SIG GRASS y el directorio de mapas

Iniciar GRASS GIS

En Pannel de arranque establecido SIG Directorio de datos de ruta de acceso a bases de

Visualizar mapa de elevacin:

Pantalla de elevacin como un mapa con sombra de color:

Visualizacin de vectores lineales y puntuales mapas:

Guarde el mapa que se muestra:

vista de datos y consulta

La lnea de comandos equivalente est por debajo, a excepcin de la eliminacin de todas

La gestin de los datos

Lista nicos mapas vectoriales que comienzan con "s":

Encontrar el trabajo del sistema y el alcance espacial de coordenadas:

Para saber cmo copiar, renombrar y eliminar mapas:

Encontrar informacin sobre mapas individuales capas y sus atributos:

Todos los datos (mapas) se almacenan en el directorio de la base de datos de GRASS.

Importacin y exportacin de datos raster y vectoriales

File > datos vectoriales Importar > formatos de importacin comunes .

File > Exportar > Los formatos ms comunes

A continuacin, puede ver la roadsmajor en ArcGIS.

En MS Windows a menudo es ms fcil hacer la importacin a travs de interfaz grfica de

path_to_my_file \ lc96ras_cut.tif "o" /path_to_my_file/lc96ras_cut.tif "dependiendo de la

Exportacin de archivos de mapa de bits como rejilla ASCII y GeoTIFF

Opcional: Aadir capa WMS

En la barra de herramientas Administrador de capas Aadir capa de servicios web .

modelos de datos geoespaciales

Informacin general del SIG GRASS y manual

Iniciar GRASS GIS

En Pannel de arranque establecido SIG Directorio de datos de ruta de acceso a bases de

Volver a muestrear con una resolucin ms alta

El -p bandera para g.region se utiliza para imprimir la regin apenas fijamos.

trama utilizando el mtodo del vecino ms cercano.

Mostrar el mapa resampled aadiendo "elev_ned10m_nn" al Administrador de capas en

Mostrar el mapa resampled aadiendo "aspect_ned10m_nn" al Administrador de capas o

Guarde el mapa que se muestra y explicar lo que est pasando en el informe y en qu se

Ahora, usando DEM reinterpolate bilineal y la interpolacin bicbica. Compruebe las

Volver a muestrear con una resolucin ms baja

A continuacin, cambiar la resolucin de la regin y la elevacin volver a muestrear (que es

Eliminar o desactivar los mapas de uso del suelo, elevacin y de aspecto.

La conversin entre tipos de datos vectoriales

Lmites del censo de visualizacin utilizando la GUI: Aadir vectorial "census_wake2000"

d.vect elev_ned_contpts tamao co = = icono marrn bsico / punto = 3

Convertir a partir del vector a raster

A continuacin, retire todas las capas del Administrador de capas .

Ahora convertir lneas en el vector "calles" de la trama. Ajuste la resolucin y el uso de

Si no lo ha hecho, agregue eliminar todas las otras capas de mapas de administrador de

Convertir de la trama al vector

Convertir lneas de trama de lneas vectoriales.

comparar visualmente el resultado con corrientes de fotografas areas digitalizadas.

visualizacin de datos y visualizacin

Informacin general del SIG GRASS y manual