Analisis Espacial Con Datos Raster en Spatial Analyst-Libre
Analisis Espacial Con Datos Raster en Spatial Analyst-Libre
Analisis Espacial Con Datos Raster en Spatial Analyst-Libre
NDIC E
1. E l Mo d e lo d e Da to s Ra ste r. 2. Fue nte s d e o rig e n d e d a to s ra ste r. 3. Bre ve intro d uc c i n al tra ta mie nto y func io ne s d e a nlisis ra ste r: Anlisis L o c al, Fo c al, Zo nal y Glo b al. 4. Alma c e namie nto d e Ra ste r e n Ge o d a ta b a se . 5. Intro d uc c i n a la Exte nsi n Spa tia l Analyst pa ra Arc GIS De skto p. 6. Esp e c ific a c i n d e las o pc io ne s de a nlisis d e Spa tia l Ana lyst. 7. Ob te nc i n d e l MDE me d ia nte la he rramie nta 3D Ana lyst. 8. Ob te nc i n d e lo s pa rme tro s mo rfo m tric o s d e rivad o s d e l MDE c o n Sp a tia l Analyst. 9. Func io ne s d e a nlisis Ra ste r c o n Spa tia l Analyst. 10. E l uso d e Arc T o o lb o x y Mo d e l Build e r p a ra e l a n lisis e spa c ial d e d a to s ra ste r. 11. Bre ve re fe re nc ia d e a plic a c io ne s re ale s.
Pgina 2
Suelen englobar distintos modelos (raster/vectorial/TIN) y tienen mltiples usos. Tabla 1. Modelos de dat os en GI S.
Pero por qu es t an difcil saber si el fenm eno que vam os a analizar es m ej or analizarlo m ediant e un m odelo de dat os rast er o un m odelo de dat os vect orial? Es aqu donde ent ran en j uego los lm it es del fenm eno analizado, as desde un punt o de vist a concept ual, en general los elem ent os paisaj st icos nat urales, norm alm ent e no present an bordes m arcados sino zonas de t ransicin, por lo
Pgina 3
D e fin icin de Rast e r : Se t rat a de un Modelo de dat os en el cual las ent idades geogrficas son represent adas usando celdas, generalm ent e cuadrados ( aunque algunos sist em as utilizan ot ras figuras geom t ricas com o t ringulos o hexgonos) , ordenados conform ando una grilla regular. Un rast er es esencialm ent e com o una m at riz bidim ensional ( filas y colum nas) cuyo origen se sit a en el ext rem o superior izquierdo, al cont rario de lo que ocurre con el Sist em a de Coordenadas Cart esiano cuyo origen se encuent r a en el punt o inferior izquier do ( v ase Fig. 2) . Un r ast er puede t ener una o m s bandas de inform acin ( ent indase Fot ografa Area con t res bandas una por cada color visible RGB; im genes de Sat lit e m ult iespect rales o hiperespect rales, et c.) y la m nim a ent idad en la que se descom pone dicho rast er o grilla es el pxel que perm it e alm acenar t res t ipos de variables ( X o Longit ud, Y o Lat it ud, y Z o variable que represent am os m ediant e dicho rast er) . As el m odelo de dat os rast er se cent ra m s en las propiedades y variaciones de las variables en el espacio que en la repr esent acin discret a de los elem ent os que lo conform an.
Fig. 1. Form at o Rast er a. Sim ple, b. Multiespect ral. Todo ello conlleva un problem a derivado de la falt a de exact it ud a la hora de delim it ar los elem ent os represent ados, y la nica aproxim acin que se puede llevar a cabo es m ediant e el aum ent o de la resolucin o una clasificacin del fenm eno a nivel de subpixel. Por lo que est am os ant e un pr oblem a no t ant o de exactit ud posicional, sino de resolucin del r ast er.
Pgina 4
Fig. 2. Relacin ent re el sist em a de coordenadas rast er y cart esiano. Existen diversas clasificaciones de los dat os rast er que t ienden a difer enciar los ut ilizando cr it erios com o los que det allo a cont inuacin: Segn o o Segn Segn o o o o el valor del pix el. Cont inuos. Discret os o t em ticos. el nm ero de bandas o variables alm acenadas en el rast er. el tipo de dat os alm acenado y la fuent e de procedencia I m genes de Sat lit e o Rast er espect rales ( incluye fot ografa area) . Modelos Digit ales del Ter reno. Rast er t em ticos. Fot ografa convencional ( podra incluirse dent ro del prim er grupo)
A) Segn el valor del pxel: Cont inuos: se refiere a dat os que se encuent ran dist ribuidos de form a cont ina com o su nom bre indica en cualquier lugar de la superficie t errest re ( por ej em plo: t em perat ura, precipit acin, alt it ud, pr esin at m osfrica, et c.) y por t ant o cualquier pix el t endr un v alor de dicha v ar iable. Discret os o t em ticos: por ot ro lado t enem os los rast er discret os que proceden de una esquem at izacin o discret izacin del fenm eno est udiado de acuerdo a un crit erio, o se t rat a de represent aciones de fenm enos cuyos lm it es quedan perfect am ent e est ablecidos (com o es el caso de un rast er de usos del suelo, t ipologa de suelos, et c.) . El proceso de discretizacin de una variable cont inua pasa por la agrupacin de los pxeles que com ponen el rast er en det erm inadas unidades y la post erior asignacin de un valor com n a t odos los pxeles que com ponen cada clase o int ervalo, lo que se ha venido a denom inar RECLASI FI CACI N de un rast er cont inuo para la obt encin del un rast er discret o.
Pgina 5
Fuente: The Remote S ensing Tutorial: http:/ / www.fas.org/ irp/ imint/ docs/ rst Fig. 3. El espect ro elect rom agnt ico y las im genes m ult iespect rales ( im genes de sat lit e y fot ografas areas) e hiperespect rales.
C) Segn el t ipo de dat os alm acenado y la fuent e de procedencia: Siguiendo est e crit erio podem os diferenciar t em t icam ent e los rast er de acuerdo a la fuent e de infor m acin ut ilizada as com o la finalidad del anlisis ( conv ersin v ect or a r ast er, ex t raccin de inform acin o reconocim ient o de pat rones, et c.). As t enem os:
Pgina 6
Fig. 4. I m agen de Satlit e Landsat 7 ETM+ . Zona del Torcal de Ant equera ( Mlaga) o Modelos Digit ales del Terr eno. Rast er cuya car act erstica pr incipal es la represent acin de un fenm eno cont inuo sobre la superficie terrest re m ediant e una m alla de pxeles que se com ponen de una com ponent e espacial ( t am ao del pixel y localizacin del m ism o en coordenadas XY) , y com ponent e t em t ica o de at ribut os ( el valor de la variable alm acenada) . A part ir de est e t ipo de dat os rast er podr em os deriv ar ot r o t ipo de infor m acin m ediant e la aplicacin de las funciones de anlisis espacial.
Fig. 5. Modelo Digit al del Terreno ( Represent acin de un MDE) o Rast er t em t icos. Proceden de la agrupacin t em tica de los dat os cont inuos de acuerdo a un crit erio preest ablecido.
Pgina 7
Fig. 6. Rast er t em tico con t res valores de alt it ud ( m enor de 200 m et ros, de 200 a 700 y m s de 700 m et ros) pr ocedent e de la r eclasificacin del m odelo digit al de elev aciones de la Fig. 5.
Fot ografa convencional ( podra incluirse dent ro del prim er grupo) . Verem os un ej em plo de est e t ipo de dat os que puede alm acenarse no com o un rast er georreferenciado, sino com o un at ribut o rast er en el m odelo de dat os de geodat abase, o com o un hiper enlace en cualquier clase de ent idad (Shapefile) .
Centrndonos m om ent neam ent e en los rast er m ult iespectrales debem os at ender a la organizacin de dichos rast er de acuerdo a la siguient e ilust racin, donde se int roduce ot ro ej e cart esiano derivado de la seal det ect ada en los pxeles de cada banda de inform acin.
Pgina 8
Fig. 8. Ej es en el espacio n- dim ensional de im genes m ult iespect rales. As de izquierda a derecha y de arriba abaj o t endram os respect ivam ent e las colum nas y filas que com ponen cada una de las bandas, y cada una de st as variarn en funcin del at ribut o represent ado en el ej e Z que nos m uest ra la variabilidad t em t ica de un conjunt o rast er Mult iespect ral, que ira desde 1 hast a n bandas ( 3 en el caso de una fot ografa area o fot ografa conv encional, hast a cent enares de bandas par a im genes hiperespect r ales.) .
ESTRUCTURA Y CARACTER STI CAS MS SI GN I FI CATI VAS DE UN RASTER Com o se ha com ent ado un fichero rast er es sim ilar a una m alla de pxeles en cuya posicin se alm acena un v alor ( usualm ent e un Nivel Digit al [ ND] dada una codificacin especfica) . Est e v alor suele asociarse por convencin a un color conform ando la apariencia visual de dicho fichero en un ordenador, pero en ocasiones est o no es del t odo ciert o ya que en det erm inadas ocasiones un pixel no alm acena un valor o caract erst ica visual, sino ot ro elem ent o, a veces abst ract o com o puede ser la alt it ud o la t em perat ura, no percept ible visualm ent e ( pero s su m anifest acin) . Generalm ent e el alm acenam ient o de dicha inform acin en las com put adoras se realiza recorriendo dicha m alla de izquier da a derecha y de arriba abaj o desde la esquina superior izquierda del rast er, para post eriorm ent e reconst ruir dicha im agen el soft ware, asignando un color de acuerdo a una palet a de colores especificada, aunque exist en ot ros t ipos de alm acenam ient o rast er com o los que se det allan a cont inuacin: Est r uct ur a s de Rast er Sim ple s : - En u m er a cin Exh a ust iva . Se t rat a de alm acenar uno por uno el valor de cada celda, de acuer do con una secuencia, generalm ent e fila a fila a par tir de la celda superior izquier da. Fig. 9. - Codifica cin por gru pos de longit u d va r iable ( Run Le n gt h Encodin g) . La codificacin no se hace valor a valor, sino por grupos de valores iguales seguidos. Fig. 10. Est r uct ur a s de Rast er Jer rqu ica s: - Ar boles Cu at e rn a r ios . Una derivacin de las est ruct uras convencionales para ordenar y codificar dat os rast er es la de rboles cuat ernarios ( quadt rees) . Se t rat a de operar en una m ism a capa rast er con dist int os t am aos celdas, lo que significa t rabaj ar con una resolucin variable, algo poco frecuent e. Dicha est ruct ura es el result ado de una divisin del espacio basada en est ruct uras j errquicas de 4 cuadrant es. El espacio inicial es una m at riz regular de
Pgina 9
Pgina 10
Respect o a st e lt im o punt o de la referenciacin m ediant e ficheros world vam os a punt ualizar com o se crea un fichero w orld de acuerdo a una im agen area.
Pgina 11
Fig. 12. Ej em plo de fichero world. La denom inacin y alm acenam ient o de dicho ficher o se har siguiendo las siguient es conv enciones de nom enclat ura del fichero: Usar el m ism o nom bre de la Fot ografa Area u Ort ofot ografas, I m agen de Sat lit e o fichero raster. Com o ext ensin usarem os la m ism a pero agregarem os la let ra w al final de la ext ensin y opcionalm ent e elim inarem os la segunda let ra de la ext ensin original. Quedando por t ant o del siguient e m odo par a los ficheros siguient es: o JPG = JGW. o TI F = TFW. o SI D = SDW. o
Llegados a est e punt o en el prxim o epgrafe verem os cules son los m t odos m s ut ilizados para la obt encin de dat os raster.
Pgina 12
Fig. 13. Escner Fot ogram trico y de sobrem esa. I m gen es de Sa t lit e o fot ogr afa a rea ( c m a ra s m t rica s o digit a le s) . La t eledet eccin se considera com o la ciencia que perm it e obt ener inform acin acerca de la superficie t errest re sin est ar en cont act o con ella. Para ello se det ect a y graba la energa em it ida o reflej ada en sensores aerot ransport ados para post eriorm ent e enviarlas a cent ros de t rat am ient o especficos para su uso final dent ro de las aplicaciones GI S. Est e proceso conlleva la int egracin de varios elem ent os ent re los que dest acan: o o o o Fuent e de Energa. At m sfera. El objet o que int er acta con dicha r adiacin. El sensor ( que puede ser elem ent o act ivo o pasivo) .
Pgina 13
Obt e ncin de m ode los digit a le s de eleva cion es m edia nt e t r iangu la cin y post e r ior r ast e r iz acin . Est e proceso se basa en la obt encin de un Modelo Digit al de Elevaciones siguiendo uno de los procesos m s generalizados ent re los tcnicos GI S dada su facilidad de obt encin, as com o la facilidad de acceso a los dat os de origen ( Curvas de Nivel) . El proceso que se det allar m s adelant e (vase epgrafe Obt encin del MDE m ediant e la herr am ient a 3D Analy st ) , est basado en las siguient es fases principalm ent e: o o o Obt encin de las curvas de nivel a una equidist ancia dada de una zona concret a. Triangulacin de las Curvas de Nivel obt eniendo as el TI N o Red I rregular Triangulada. Conversin del TI N a un m odelo de dat os Rast er o Teselar.
D a t os Lidar ( Ligh t D e t ect ion a n d Ran gin g ). El sist em a Lidar em it e pulsaciones de luz que reflej an en el t erreno o cualquier obj et o en alt ura y la pulsacin de regreso se conviert e de fot ones a im pulsos elctricos. La resolucin de los dat os posicionales ( x, y) puede llegar a precisiones de 1 m et ro en la horizont al y 15 cm en alt ura, y las posibilidades que br inda perm it en r ealizar m ediciones en condiciones clim at olgicas y de ilum inacin que no seran fact ibles m ediant e la r ecoleccin de fot ogr afas areas, posibilit ando t am bin la obt encin de un m odelo digit al de Superficie incluy endo elem ent os am bient ales y ant rpicos ( m asas forest ales, lechos fluviales, edificaciones, et c.)
Com o vem os son m lt iples las fuent es de dat os para obt ener inform acin en form at o rast er. En los prxim os epgrafes se har uso principalm ent e de curvas de nivel, es decir, dat os en form at o vect orial que se t riangularn y a partir de los cuales obt endrem os un Modelo Digit al de Elevacin.
Pgina 14
Fig. 14. Tipologa de Anlisis y Funciones Locales. Fu n cion e s u Ope r a dor e s Foca le s : se denom inan as o de vecindad, ya que con ellas obt endrem os para cada pixel un valor en funcin de los valores alm acenados en dichos pxeles y en los vecinos o aledaos al m ism o. Para est e t ipo de anlisis usam os lo que se denom ina una m at riz de recorrido que puede est ar com puest a por los pxeles m s prxim os ( operadores de vecindad inm ediat a) , o m s alej ados al pixel analizado ( operadores de vecindad ext endida) y el rango de la m at riz ser algo que debam os determ inar nosot ros, aunque por defect o se suele ut ilizar una m at riz de 3x3 pxeles. Com o ejem plos de este t ipo de funciones t enem os la obt encin de est adst icas respect o a los pxeles vecinos, aplicacin de filt ros sobre im genes ( m ediana, m oda, det eccin de bordes, et c.) as com o el clculo de parm et ros com o la pendient e, orient acin, som breado, a part ir de un MDE. Para el caso de las operaciones de vecindad ext endida, t endrem os en cuent a un m ayor nm ero de pxeles cont iguos com o es el caso del anlisis de visibilidad o dist ancia en lnea rect a o de cost o.
Fig. 15. Tipologa de Anlisis y Funciones Focales. Fu n cion e s u Oper a dore s Zona le s: en est e t ipo de operaciones el valor del pixel en una capa result ant e del anlisis depender del valor del pixel en la capa de origen, as com o el valor de los
Pgina 15
Fig. 16. Tipologa de Anlisis y Funciones Zonales. Fu n cion e s u Ope r a dor e s Globa les : por lt im o t enem os las operaciones de m bit o global que t endrn en cuent a la t ot alidad de los pxeles de la im agen o rast er, y ent re st as se encuadran los m t odos de int erpolacin.
Fig. 17. Tipologa de Anlisis y Funciones Globales. Adem s de est os operadores exist en ot ro t ipo de anlisis que se realizan en funcin de la et apa dent ro del proceso de int egracin de los dat os en un GI S, com o son las operaciones de relocalizacin de una im agen rast er en coordenadas ( georreferenciacin) o las correcciones de sus niveles digit ales, o la realizacin de m osaicos, que aunque bien podran incluirse dent ro de las operaciones ant eriores, han venido a ident ificarse dent ro de las operaciones de pre- procesado rast er. Com o aclaracin debem os decir que no t odos los pxeles t endrn asignado un valor, bien porque han sido m odificados, bien por falt a de recepcin por el equipo recept or ( im genes de sat lit e, Lidar, et c.) , indicando que no hay m edicin alguna para dicho rast er.
Pgina 16
Pgina 17
Est e t ipo de dat os se puede alm acenar en varias m odalidades dent ro de una geodat abase lo cual redunda en m lt iples beneficios asociados al m t odo de alm acenam ient o. Ent re los beneficios que podem os encont rar del uso y alm acenam ient o de dat os rast er en geodat abase t enem os: Uso de rast er t ant o a nivel de geodat abase personal com o m ult iusuario, cada una con sus part icularidades en cuant o a gest in. Mant enim ient o de grandes volm enes de dat os. Posibilidad de uso de dat asets individuales o en m osaico, o cat logos de dat os raster. Rpido despliegue de dat os rast er dada la posibilidad de crear pirm ides ( m t odo de com presin visual de dat os rast er) en nuest ros rast er. Posibilidad de ext raccin de form a sencilla de ot ros dat os a part ir de dichos rast er. Posibilidad de com presin con/ sin prdida. Et c.
A veces no result a t an obvio saber el m odo de alm acenar, organizar y analizar los dat os rast er en una base de dat os geogrfica, es por ello que darem os un breve repaso a los t ipos de organizacin y alm acenam ient o de dat os rast er en geodat abase y verem os algunos ej em plos y su organizacin.
Pgina 18
Los t r es tipos de alm acenam ient o r ast er en geodat abase inclu yen: Ra st e r D a t a se t s ( Conj unt o de Dat os Rast er): se denom ina rast er Dat aset o conj unt o de dat os rast er a cualquier form at o rast er vlido organizado en una o m s bandas, siendo cada banda un Array o m at riz de pxeles ( celdas) donde cada uno de st os dispone de un valor asociado a la t em t ica represent ada por dicho rast er. I ndependient em ent e del nm ero de bandas, un rast er Dat aset puede anexarse horizont alm ent e a m odo de m osaico o perm anecer individualm ent e. Ra st er Ca t a log ( Cat logos Rast er) : se denom ina com o cat logo rast er a la coleccin de rast er dat aset s en una t abla dent ro de la geodat abase en la que cada regist ro r epresent a un Dat aset rast er individual, y aunque en la m ay ora de los casos se usa para alm acenar rast er individuales ( uno por regist ro) dent ro de un m osaico ( im genes cont iguas con una zona de solape) , aunque es posible alm acenar t am bin dist int as im genes de una m ism a zona pero de periodos t em porales dist int os, as com o ot ro t ipo de im genes. Dent ro de dicho cat logo rast er cada im agen t endr su propio sist em a de referencia espacial, y dichos rast er dat aset s podrn ser alm acenados de dos form as dist int as dent ro de una Geodat abase: Rast er Gest ionado por la Geodat abase: convirt iendo las im genes en form at o Erdas I m agine ( I MG) y alm acenndolas direct am ent e en una car pet a al m ism o nivel que la geodat abase. Com o beneficio t enem os la gest in cent ralizada por la geodat abase, per o com o cont rapart ida t enem os el aum ent o del peso de nuest ra geodat abase sobre t odo cuando de reas ext ensas se t rat a. Rast er No Gest ionado por la Geodat abase: no se conviert en las im genes a un form at o nativo com o es ERDAS I m agine ( I MG) sino que se m ant ienen en la ubicacin y con el form at o original, y slo se alm acena en una t abla ( cam po NAME) de la geodat abase las referencias a la ubicacin de cada rast er Dat aset . Adem s de la incorporacin de t odos los rast er en el cat logo nos perm it e ver las im genes o un ndice vect orial cuy a finalidad es la posibilidad de que podam os cargar slo las im genes que necesit em os a una escala dada, o que podam os realizar selecciones t em t icas y espaciales dent ro de ese ndice vect orial com o si de ent idades vect oriales se t rat ase. At r ibu t os Ra st er : los at ribut os rast er se alm acenan en la geodat abase a m odo de enlace ent re un cam po de t ipo Rast er en una Feat ure Class y la im agen alm acenada int ernam ent e en dicha geodat abase. Es sim ilar al hiperenlace en Shapefile pero en el que la referencia es int erna y no se alm acena el enlace a una im agen ext erna a la geodat abase. A cont inuacin verem os brevem ent e los pasos para crear cada uno de est os alm acenes de dat os rast er. A la hora de saber qu t ipo de m odalidad de alm acenam ient o seleccionar debem os t ener en cuent a los siguient es aspect os: Si slo se va a hacer uso de un ra st er D a t ase t deberem os alm acenarlo com o t al, no dentro de un cat logo rast er. Si querem os t rabaj ar con un m osaico rast er y no necesit am os disponer de las zonas de solape podem os opt ar por un ra st e r Da t a set a m odo de m osa ico . Si nuest ra necesidad es disponer de un m osaico rast er pero necesit am os conservar las zonas de solape individual opt arem os por un cat logo r ast er . Si querem os disponer de un conjunt o de rast er dat aset s organizados t em poralm ent e ( por fechas de capt ura) pero referidos a una m ism a zona usarem os un ca t logo r a ste r . Si necesit am os disponer de un conjunt o de im genes cont iguas pero no es necesario car gar t odas la im genes al m ism o t iem po podem os usar un ca t logo r a ste r , ya que
Pgina 19
Para ver com o se alm acenan y crean los rast er dat aset s, cat logos y at ribut os rast er en geodat abase vam os a crear un ficher o de est e t ipo ( Geodat abase Personal en nuest ro caso) . Los pasos a seguir son: Pulsar bot n derecho del rat n sobre la carpet a de dest ino New Personal Geodat abase
Fig. 19. Creacin de una geodat abase personal. Una vez disponem os de nuest r a geodat abase procedem os a renom brarla com o cream os conv enient e, en nuest r o caso hem os usado Alm a ce n a m ie n t o Ra st er .m db .
Cr e acin de u n Rast e r Da t a se t .
Para crear un rast er Dat aset t enem os dist int as opciones: Desde Ar cCa ta log . Bot n derecho sobre la geodat abase New Rast er Dat aset .
Fig. 20. Cr eacin de un Rast er Dat aset desde ArcCat alog. Desde Ar cToolbox . Dat a Managem ent Tools Rast er Creat e Rast er Dat aset .
Pgina 20
Desde la Ln e a de Com a n dos . Bast a con t eclear creat erast erdat aset y proseguir rellenando los dat os solicit ados. Bast ar con especificar en el cuadro de dilogo los siguient es parm et ros: Ou t put Loca t ion ( Localizacin del rast er Dat aset ) . Por defect o la geodat abase sobre la que se ha pulsado bot n derecho. Ra st er D a t a se t N a m e w it h Ex t en sion ( Nom bre del Rast er Dat aset con extension) . Aqu especificarem os el nom bre del rast er Dat aset y las ext ensiones disponibles son ( im g, t if, bm p, j pg, et c.) . Cellsize ( Opcional) . Es el t am ao de celda que t endr nuestr o r ast er Dat aset . Pix el Type ( t ipo de dat os soport ado por el pixel) . Aqu especificarem os dent ro de los est ndares de alm acenam ient o de dat os rast er el t ipo de valor que podrem os alm acenar. Spa t ial Refe re n ce for Rast e r ( Opt ional) . Sist em a de Referencia Espacial de nuest ro rast er dat aset . N u m be r of Ban ds ( Nm er o de Bandas) . Nm er o de bandas que com pondr n el r ast er ( 1 a n) . Adem s de est os parm et ros exist en ot ros parm et ros de configuracin del rast er dat aset s que no se suelen especificar, salvo casos excepcionales, y que pueden consult arse en la docum ent acin incluida con ArcGI S Deskt op.
Pgina 21
Para im port ar rast er dat aset s en Geodat abase t enem os las siguient es opciones: Desde Ar cCat a log . Bot n derecho del rat n sobre la geodat abase I m port Rast er Dat aset s / Rast er Dat aset s ( Mosaico)
Fig. 22. I m port acin de Rast er Dat aset s individuales o a m odo de m osaico en ArcCat alog Desde Ar cToolbox : Conversion To Geodat abase Rast er t o Geodat abase ( Mult iple)
Fig. 23. I m port acin de Rast er Dat aset m lt iples a Geodat abase o Rast er Cat alog. Desde la Ln e a de Com an dos : int roduciendo la sent encia rast ert ogeodat abase e int roduciendo el rest o de par m et ros.
Pgina 22
Para crear un Cat logo Rast er t enem os dist int as opciones: Desde Ar cCa ta log . Bot n derecho sobre la geodat abase New Rast er Cat alog.
Fig. 24. Generacin de un Cat logo Rast er desde ArcCat alog. Desde Ar cToolbox . Dat a Managem ent Tools Rast er Creat e Rast er Cat alog.
Pgina 23
Fig. 26. Ej em plo de cat logo rast er a m odo de m osaico con ArcCat alog.
Pgina 24
Para crear un Atribut o Raster debem os proceder del siguient e m odo: En una Feat ure Class creada de ant em ano vam os a agregar un cam po de t ipo Rast er, pulsando bot n derecho del rat n sobre la Fea t ure Cla ss Propieda des Fie lds . Ah agregarem os el cam po y especificarem os com o t ipo de dat os Rast er. ( Vase Fig. 26) .
Fig. 27. Creacin de un cam po de t ipo Rast er en una Feat ure Class m ediant e ArcCat alog. Seguidam ent e vam os a ArcMap y en m odo edicin pulsam os el bot n At ribut os y en el cam po que hem os agregado pulsarem os el bot n izquierdo del rat n y seleccionarem os m ediant e las opciones que se m uest ran en la Fig. 27 la im agen a cargar com o at ribut o.
Fig. 28. Adicin de una fot o com o at ribut o r ast er. Adem s de est e conj unt o de operaciones disponem os de otro t ipo de operaciones de m ant enim ient o de dat os rast er com o son los cam bios de resolucin ( D a ta M a n a ge m e n t Tools Rast e r Re sa m ple ) , elabor acin de m osaicos pr opiam ent e dichos ( D a t a M a na ge m e n t Tools Raste r M osa ic t o N e w Ra st er / M osa ic ) que debido a la ext ensin que conllevara com ent ar t odas las operaciones sobrepasara la ext ensin de est a gua. En fut uras versiones se incluirn dichas operaciones.
Pgina 25
Adem s de las herram ient as incorporadas en la propia barra de herram ient as de Spat ial Analyst , se incorporan a las herram ient as de ArcToolbox, num erosas funciones de anlisis ent re las que dest acan, oper ador es arit m t icos, oper ador es booleanos, anlisis hidrolgicos, y un lar go et ct era. Operables desde cualquiera de los com ponent es del ent orno de geoprocesam ient o de ArcGI S Deskt op t ales com o las pr opias herram ient as de ArcToolbox , Model Builder o la Lnea de Com andos.
Operadores de Distancia (Operadores Focales) Operadores de obtencin de Densidades (Operadores Focales) Operadores de Interpolacin de Superficies (Operadores Globales) Operadores de Obtencin de Superficies (Operadores Focales) Estadsticas de C eldas (Operadores Locales) Resmenes Estadsticos de Pxeles Vecinos (Operadores Focales) Estadsticas Zonales (Operadores Zonales) Reclasificacin (Operadores Locales) Calculadora Raster Conversiones (Raster a Vector Vector a Raster) Opciones (Configuracin de las herramientas de S patial Analyst)
Fig. 29. Barra de Herram ient as de Spat ial Analyst para ArcGI S Deskt op y agrupacin de funciones de Anlisis Espacial de acuerdo a la clasificacin del capt ulo 3.
Pgina 26
1 2
Fig. 30. Calculador a Rast er de ArcGI S Spat ial Analy st . La m ayor part e de las operaciones que se realizan son arit m t icas com o la sum a de un rast er con ot ro, la ext raccin de pxeles de un rast er basado en un crit erio especfico, la com paracin de un rast er con ot ro, et c. En las siguient es lneas se puede ver algunos ej em plos de la funcionalidad de dicha herr am ient a: [ Ra ste r ] + 2 5 [ Ra ste r 1 ] / [ Rast e r2 ]
Resultado: Adiciona 25 unidades a los pxeles del raster de ent rada. Resultado: Divide los pxeles de un raster respect o al ot ro cada pxel con su respect ivo en la ot ra capa. Resultado: Conviert e a valores enteros los pxeles de un raster cuyos dat os son alm acenados m ediant e dat os de t ipo Float . Resultado: obt iene un rast er donde los valores son la raz cuadrada del raster de ent rada.
I n t ( [ Ra st er ] )
Sqr t ( [ Ra st er ] )
Pgina 27
Com o h em os com ent ado adem s de est as funciones arit m t ico-lgicas y ot ras m s avanzadas ( logar t m icas, exponenciales, t r igonom t ricas o booleanas) t enem os ot ras funciones predet erm inadas para clculos m s especficos, ent re las que se encuent ran algunos de los ejem plos que det allo a cont inuacin ( para m s inform acin acerca de las herram ient as disponibles Consult ar la Docum ent acin de Ay uda de Ar cGI S Deskt op en Ext ensions Spat ial Analy st Spat ial Analy st Funt ional Refer ence) : M osa ic ( [ Rast er 1 ] , [ Ra ste r 2 ] , [ Rast e r3 ] [ Ra ste rN ] )
Result ado: obt enem os un m osaico con t odas los rast er de ent rada a part ir de una zona de solape.
Aspect ( [ Ra st er ] )
Resultado: obtenem os un raster de orient acin de la pendient e a part ir de un m odelo digital del terreno ( frecuentem ent e un MDE).
Com o vem os la m ayora no son com plicadas pero det erm inadas funciones pueden requerir de dist int os parm et ros que perm it en personalizar el result ado obt enido com o por ej em plo la siguient e operacin: Operador Mxim a Pr obabilidad o v er osim ilit ud (Max im um Likelihood Classifier) Funcin genrica: M LCla ssify( < st ack > , < sign a t u re > , { a _ pr iori_ f ile} , { o_ re j ect _ gr id} ) { r ej ect _ fr a ct ion } , { EQUAL | SAM PLES | FI LE} ,
Aplicacin: m lcla ssify ( redla nds, redl8 , # , f ile, w e ight s_ file , r ej ect 2 _ gr d ) Donde los par m et r os especificados nos perm it en per sonalizar la obt encin de una clasificacin supervisada por el m t odo de m xim a verosim ilit ud, algorit m o am pliam ent e ut ilizado en Teledet eccin, y donde com o vem os algunos argum ent os que no queram os especificar al ser opcionales nos apar ecer n con el sm bolo # ( en est e caso el parm et ro { rej ect _fract ion} o la porcin de la celda que perm anecer sin asignarse a una clase debido a un a probabilidad baja de asignacin correct a. Para int roducir est os dat os bast a con seleccionar el rast er u operando m ediant e doble clic en la zona 1, y seleccionar los operador es en la zona 2 pulsando un sim ple clic sobr e el operador a ut ilizar , pudiendo com binarse de form a correct a hast a com plet ar la funcin o algorit m o de anlisis deseado, con la posibilidad que nos brinda de realizar varios anlisis en una sim ple funcin com binada en la zona 3.
Pgina 28
Fig. 31. Opciones de Configuracin de Spat ial Analyst Pest aa General. Ficha Ext ent : en est a ficha slo t endrem os que especificar los lm it es de anlisis de nuest ro t rabaj o con Spat ial Analyst m ediant e la especificacin del parm et ro siguient e:
Pgina 29
Fig. 32. Opciones de Configuracin de Spat ial Analyst Pest aa Ext ent .
Ficha Cell Size: por lt im o, est a pest aa nos perm it ir especificar el t am ao de celda de salida de nuest ro anlisis rast er, pudiendo especificar algunas de las siguient es opciones dent ro del parm et ro siguient e: Ana lysis Cell Size (t am ao de celda o pixel de anlisis) donde podem os seleccionar: Max im um of I nput s: el valor de t am ao de pixel m s grande usado por los rast er de ent rada. Minim um of I nput s: el valor de t am ao de pixel m s pequeo usado por los rast er de ent rada. Sam e as Below: usarem os un t am ao de pix el que especificarem os en la casilla Cell Size. Por t ant o, con est os parm et ros t am bin se act ualizarn los valores del nm ero de filas y colum nas del rast er de salida.
Pgina 30
Fig. 33. Opciones de Configuracin de Spat ial Analyst Pest aa Cell Size. Una vez especificados est os parm et ros est arem os en condiciones de usar est a ext ensin para cualquier anlisis rast er, pero ant es de com enzar y ya que la capa que obt endrem os en el prxim o epgrafe ser usada para proyect os post eriores, debem os hacer uso de ot ra ext ensin disponible con ArcGI S Deskt op ( disponible m ediant e adquisicin de licencia apar t e) denom inada 3D Analy st que nos va a perm it ir, ent re ot ras cosas, derivar a part ir de una capa vect orial de curvas de nivel y punt os acot ados un TI N y post eriorm ent e una superficie de alt itudes, conocida norm alm ent e com o Modelo Digit al de Elevaciones ( MDE) un subt ipo de los Modelos Digit ales del Ter reno ( MDT) .
Pgina 31
Pasem os a la descripcin det allada de los pasos que debem os llevar a cabo para realizar nuest ro MDE en form at o raster. Ge n e ra cin de la s Cu r vas de N ive l . Aqu slo m encionar que podem os digit alizarlas m ediant e digit alizacin m anual o asist ida m ediant e Ar cScan. Tria ngu lacin ( TI N) : para realizar la t riangulacin usarem os el m dulo 3D Analyst a t ravs de las funciones de t riangulacin localizadas en el m en de la barra de herram ient as 3D Analyst Creat e/ Modify Tin Creat e Tin from Features.
Fig. 34. Her ram ient a par a t riangular a part ir de dat os vect oriales ( punt uales, lineales y poligonales) Una v ez seleccionada dicha herr am ient a nos solicit ar los siguient es dat os (v ase Fig. 35) : La ye r s ( Capas): son las capas de ent rada, bien disponibles en la Tabla de Cont enidos, o seleccionables dent ro de una unidad de disco o red. Una v ez seleccionadas podrem os especificar los parm et ros de configuracin para cada una de st as. Se t t in gs for Se le ct e d layer ( Parm et ros de configuracin para la capa seleccionada) : donde especificarem os par a cada capa, dependiendo de la funcin que t enga dent ro del
Pgina 32
Fig. 35. Dat os solicit ados par a la Triangulacin. El TI N obt enido se puede observar en la siguient e ilust racin com o quedara aproxim adam ent e:
Pgina 33
Fig. 36. Obt encin del TI N t ras int roducir los dat os solicit ados por la herram ient a Creat e Tin from Feat ures. De fondo podem os observar las curvas de nivel de las que procede. Ra st eriza cin del TI N ( MDE) : una vez obt enido el TI N, pasarem os a calcular la superficie cont inua referida a los valores de alt it ud del t errit orio analizado. Para ello usarem os las herram ient as de conversin t ridim ensional, disponibles en 3D Analyst Convert Tin t o Rast er. Donde especificarem os los siguient es parm et ros: I nput TI N ( TI N de ent rada) : aqu especificarem os el TI N obt enido en el paso ant erior. At t ribut e ( At ribut o) : aqu especificam os el valor de los at ribut os que cont endr el MDE en form at o rast er ( cada pixel t endr un valor dependiendo de est a seleccin) , que puede ir desde valores de elevacin, hast a orient aciones, pasando por valores de pendient e ( en grados o porcent aj e) del TI N. Z Fact or: es un fact or de exageracin del relieve usado para zonas donde no es posible discrim inar de form a adecuada las diferencias alt it udinales y m orfogrficas. Por defect o se dej a 1,0 ya que sera com o dej ar el valor de cada pixel com o est ( ya sabem os un valor m ult iplicado por 1,0 queda igual) . Cell Size ( Tam ao de Celda) : por defect o el program a int roduce un valor que en m uchas ocasiones hay que m odificar. En m i caso he seleccionado 10 m et ros pero ya depender de vuest ras necesidades, y de los dat os de ent rada. Exist en ciert a discrepancia en cuant o a la t cnica o m et odologa para la seleccin de un t am ao de pixel de acuerdo a num erosos elem ent os com o son la equidist ancia, la escala de nuest ro m apa, et c. pero com o prim era aproxim acin podram os usar o bien la equidist ancia, o bien usar m t odos para ext raer la m nim a dist ancia ent re los vrt ices que conform an las curvas de nivel, un proceso algo t edioso pero que suele dar buenos result ados. Nunca debem os usar una resolucin dem asiado grosera ya que est aram os obviando inform acin, ni usar una resolucin dem asiado det allada que no aum ent e el det alle de nuest ros dat os ya que t endram os la m ism a inform acin repet ida en num erosos pxeles. Al asignar est e parm et ro aut om t icam ent e se act ualizarn los valores de las filas y colum nas del rast er de salida,
Pgina 34
Fig. 38. MDE obt enido m ediant e Triangulacin y post erior conversin a rast er.
Pgina 35
Pgina 36
Fig. 39. Herram ient a para la obt encin de cont ornos o I solneas y cuadro de dilogo donde especificarem os los parm etros. Par m et r os par a la obt encin de los cont ornos: I nput Sur face ( Rast er o superficie de ent r ada) : nuestr o MDE inicial.
Pgina 37
En base a est os parm et ros nuest r a capa de isolneas se configur ara segn los parm et ros est ablecidos en la part e inferior de dicho cuadro de dilogo. El r esult ado sera sim ilar al de las cur vas de nivel. Pe n die n te (Slope): la pendient e se considera com o la rat io m xim a de cam bio del valor de Z a t ravs de una celda y sus cont iguas. Para ello se hace uso de un filt ro o Kernel que va recorriendo el rast er para obt ener el valor de pendient e del pixel respect o a sus vecinos ( usando un filt ro de 3x3 pxeles por defect o) . Dicha herram ient a se encuent ra en Spa t ia l Ana lyst Surfa ce Ana lysis Slope .
Fig. 40. Localizacin de la funcin Slope ( Pendient e) en la bar ra de herr am ient as Spat ial Analyst y cu adr o de dilogo con los par m et ros de configuracin. El algorit m o ut ilizado se describe br ev em ent e a cont inuacin: Dado que el recorrido de los pxeles m ediant e un filt ro o Kernel se realiza de izquierda a derecha y de arriba abaj o correlat ivam ent e y t eniendo en cuent a que dicho filt ro const a de una m at riz de 3x3 pxeles, t endrem os en cuent a para cada valor del rast er de ent rada ( alt it udes en nuest ro caso) , los 8 pxeles vecinos al pixel dado, y siguiendo el t eorem a de Pit goras podrem os obt ener los siguient es elem ent os que ut iliza el algorit m o de form a int erna: Parm et ros de dist ancia ( base del t ringulo) : obt enido bien por la cabecera del fichero, bien por la pest aa Source de nuest ro Layer (Vanse las Propiedades) . Para el caso de las diagonales har falt a aplicar el m ism o t eorem a y as obt ener la dist ancia ent re el pix el cent ral y los cuat ro diagonales. Parm et ros de diferencia de alt ura: bast a con rest ar el valor del pixel cent ral y los 8 rest ant es, ese valor est disponible en el rast er de ent rada m ediant e el valor del pixel de alt it udes.
Pgina 38
dEO= [ (Zi+ 1,j+ 1 + 2Zi+ 1,j + Zi+ 1,j- 1) - ( Zi- 1,j+ 1 + 2Zi- 1,j + Zi- 1,j- 1) ] / 8dX dNS = [ (Zi+ 1,j+ 1 + 2Zi,j + 1 + Zi- 1,j+ 1) - ( Zi+ 1,j- 1 + 2Zi,j - 1 + Zi- 1,j- 1) ] / 8dy
Donde dx = la dist ancia est e- oest e a t ravs de la celda ( anchura de la celda) dy = la dist ancia nort e- sur a t ravs de la celda (altura de la celda)
Fig. 41. Ncleo de anlisis para la obtencin del rast er de pendient es.
Los parm et ros ut ilizados son los siguient es: I nput Surface ( Superficie o Rast er de ent rada): aqu especificam os el MDE inicial t al y com o hicim os para los cont ornos. Out put Measurem ent : valores de m edicin de la pendient e ya sea en grados o porcent aj e. Z Fact or ( Fact or de exageracin) : se t rat a del fact or de exageracin del relieve en caso de t rat arse de zonas con escasa variacin. Por defect o t enem os 1.0. Out put Cell Size ( Tam ao de pixel por defect o): dej arem os el m ism o que el MDE incial, ser el valor de r esolucin que t endr nuest r o rast er de pendient es final. Out put Rast er ( Rast er de Salida) : ser el nom bre del fichero rast er en form at o ESRI Grid. Ya sabem os las convencin de nom enclat ura exist ent e para este t ipo de ficheros.
Pgina 39
Fig. 42. Rast er de pendient es obt enido m ediant e Spat ial Analy st . Or ie n t acin ( Aspect ) : se t rat a del valor de la orient acin respect o al Nort e que van de 0 a 360 gr ados de or ient acin. El valor de cada celda indica la or ient acin de las lader as. El valor -1 indicar que no hay or ient acin y se t r at a de una super ficie plana. Par a obt ener dicho valor se or ient a el MDE de acuerdo a la figura siguient e. En cuant o a la localizacin de la herram ient a y sus parm et ros se m uest ran en la siguient e figura (Vase Fig. 43) . La herram ient a la podem os encont rar en Spa t ia l Ana lyst Surfa ce Ana lysis Aspect .
Grados de Orientacin
Fig. 43. Localizacin de la herr am ient a Aspect (Orient acin) y su cuadr o de dilogo.
Pgina 40
Som br e a do Topogr fico (Hillshade) : superficie que m uest ra el grado de ilum inacin de los pxeles de acuerdo a dos posiciones relat ivas, una respect o a la orient acin o ngulo acim ut al, y ot ra respect o al ngulo de elevacin solar o cenit al, que por defect o suelen usarse 45 de Ele vacin Sola r y 315 de Orie nt a cin . Dicha herram ient a se puede usar desde Spa t ia l Ana lyst Surfa ce Ana lysis H illsha de . Est e rast er result ant e que se m uest ra en la figura 46, podr ser usado bien com o com ponent e cart ogrfico que facilit ar la represent acin del relieve, o bien para clculos relacionados con la insolacin.
Orientacin
Elevacin S olar
Fig. 45. Localizacin de la herr am ient a Hillshade en Spat ial Analyst par a Ar cGI S.
Pgina 41
Fig. 46. a. Som breado t opogrfico obt enido a part ir del MDE y b. con el MDE superpuest o. Para superponer el MDE sobre el Som breado t opogrfico creando el efect o visualizado en la im agen ant erior sit uam os la capa MDE de la t abla de cont enidos sobre la del som breado, y post eriorm ent e nos dirigim os hacia las propiedades del MDE, pulsando bot n derecho del rat n sobre la capa MDE, y seleccionando Propiedades Pest aa Display Transparency , donde especificarem os un porcent aj e de t r anspar encia ( 0 100% ) .
re as de Visibilida d o Cue ncas V isu a les ( Viewshed) : se t rat a de un rast er que represent a los pxeles del rea visible desde un punt o de observacin. Para ello debem os disponer de un rast er de elev aciones y un punt o o lnea de observacin que debe est ar localizado tridim ensionalm ent e. Para poder disponer de un punt o de observacin en 3D debem os crear en prim er lugar el punt o de observacin m ediant e una Feat ure Class o un Shapefile, y convert irlo de Feat ure Class a 3D Feat ure Class m ediant e 3 D An a lyst Con ve r t Fe a tu r e s To 3 D , funcin que requiere com o parm et ros de ent rada el Shapefile o la Feat ure Class, y un MDE o superficie de referencia a part ir de la cual t om ar el v alor de alt it ud.
Fig. 47. Localizacin de la herr am ient a Viewshed y los par m et r os de configur acin.
Pgina 42
Va cia do / Re llena do ( Cut / Fill) : se t rat a de una funcin que calcula el valor de vaciado o relleno de un MDE o cualquier ot ro m odelo digit al del t erreno ent re dos fechas o sit uaciones. Para ello se basa en una operacin local ( pixel a pixel) ent re am bas sit uaciones, apareciendo com o result ado los valores de acum ulacin o vaciado con dist int os colores ( ganancias en roj o y con valores negat ivos, y prdidas con colores azules, y valores posit ivos; para el caso de los valores est ables sin prdidas o ganancias se m ost rarn en un color neut ro o gris, y su valor ser 0) respect ivam ent e. Los parm etros necesarios se lim it an al rast er de ent rada y salida ( Rast e r I n pu t y Ra st er Ou t pu t )
Fig. 49. I lust r acin de la funcin de Cut / Fill de 3D Analy st y Spat ial Analyst .
Fuente: Sistem a de Ayuda ( ESRI ArcGI S Deskt op)
Todas las capas obt enidas m ediant e est os procesos de anlisis superficial podrn ser usados post eriorm ent e para cualquier anlisis am bient al, bien para la obt encin de par m et ros derivados, com o par a su uso en un anlisis int egrado del t errit orio.
Pgina 43
Pgina 44
Fig. 51. Localizacin de la funcin para obtener la distancia en lnea rect a y sus parm et ros de configuracin. Este pr oceso nos perm it e ext raer la dist ancia exist ente desde cada punt o del t erritorio hast a el punt o de dest ino Pun t oD est ino.shp . Para ello requiere de los siguient es dat os: D ist a nce To ( Dist ancia a) : aqu especificarem os la capa que cont iene el punt o o localizacin de dest ino o referencia. M a xim u m Dist a n ce ( Dist ancia Mxim a) : es la distancia que lim it ar nuest ros clculos, a part ir de la que no obtendrem os ningn valor de dist ancia. Out put Ce ll Siz e ( Tam ao de Celda de Salida) : t am ao de salida de los pxeles del raster de distancia ( opcionalm ent e tam bin para los rast er de direccin y asignacin) . Cr e a t e D ir e ct ion ( Crear Rast er de Direcciones) : aqu especificarem os opcionalm ent e un fichero de t ipo rast er ESRI Grid que cont endr los valores de dir eccin a la que t ienen que dirigirse cualquier punto o pixel del t errit orio para llegar al punt o m s cercano, en este caso el punt o de referencia. Los valores van de 1 a 8 segn la direccin que t om en. Crea t e Alloca t ion (Crear Rast er de Asignaciones) : raster opcional que contendr un det erm inado valor relacionado con el punt o m s cercano.
El result ado para el ej em plo se m uest ra a cont inuacin: obviando el rast er de allocat ion o asignaciones dado que slo tendram os un valor al exist ir un nico punt o de referencia o dest ino.
Fig. 52. a. Raster relat ivo a la dist ancia m nim a en lnea recta desde cada pixel al punt o de dest ino y b. Rast er de direcciones.
Pgina 45
C lcu lo de la D ist an cia de Cost e : Para el clculo de la dist ancia ponderada por una superficie de coste es debem os ir a Spa t ia l Ana lyst D ist a nce Cost W eight ed . Para el clculo de la distancia ponderada por una o varias superficies de coste debem os disponer, al m enos, de los siguientes datos: Lugar de origen y dest ino. En nuest ro caso dos punt os sit uados encim a de la superficie a analizar, y que posteriorm ente nos servirn para obtener la ruta m s corta basada en dicha dist ancia. Superficie de coste o friccin: usarem os la pendient e com o principal elem ento de im pedancia, debidam ente ponderado. Superficie sobre la que t rabaj arem os. Nuest ro MDE.
Fig. 53. Localizacin de la herr am ient a par a la obt encin de la dist ancia ponder ada por el coste de t rnsit o, y parm et ros necesarios para su configuracin.
Lo prim er o que debem os obtener antes de realizar nuest ras operaciones de clculo de distancias ponderadas, es obtener el fact or o factores de ponderacin. En nuest ro caso la pendiente. Para ello procedem os a calcular la pendient e en porcent aj e o en grados, y a part ir de st a reclasificarla en clases segn su grado de afeccin sobre nuest ro m ovim ient o a t ravs de los pxeles, en definit iva, a t ravs del t errit orio. Para ello harem os uso de la funcin de Reclasificacin que se ver m s adelante en este m ism o epgrafe. Una vez reclasificada y ponderada convenientem ente se asignar al Rast er de Coste o I m pedancia. Los parm et ros solicit ados por dicha funcin son los siguient es: D ist a nce To ( Distancia a) : aqu nos pide el elem ent o de referencia a part ir del cual vam os a calcular la distancia de cost e ponderada por la pendient e (o cualquier ot ra superficie) para cualquier pixel o localizacin del t errit orio analizado. Cost Rast er ( Rast er de I m pedancia o cost e) : ser nuest ro raster de pendient e, el que posibilit e u ofrezca cualquier resistencia a nuest ro m ovim ient o. M a x im um D ist a nce ( Dist ancia Mxim a de Anlisis) : aqu podem os especificar opcionalm ent e la distancia m xim a a part ir de la cual no se t om arn m s pxeles para el clculo. Cr e a t e D ir e ct ion (Crear raster de direcciones) : podem os especificar opcionalm ent e un rast er de direcciones que t om arn los pxeles para dirigirse al punt o m s prxim o de referencia. Crea t e Alloca t ion (Crear rast er de Asignaciones) : raster de asignacin del lugar o lugares m s cercanos a un determ inado pxel. Out p ut Ra st er ( Raster de Salida) : raster cuyos valores tendrn el valor de distancia m nim o ponderado por el coste de im pedancia de la pendiente.
A cont inuacin vem os los parm et ros usados ( MDE) , Pendientes, y punt os de origen y dest ino.
Pgina 46
Fig. 54. Dat os necesarios para la evaluacin de la distancia de cost e a un punto dado (a. Localizaciones puntuales, y b. Superficie de friccin Pendient e) . El result ado de dicha operacin se puede observar a cont inuacin ( capas de direccin y distancia ponderadas por la pendient e).
Fig. 55. Capas result ant es de la operacin de distancia de cost e (ponderada por la pendiente) . a. Distancia de coste y b. Raster de Direcciones. Not a. Se ha t om ado com o referencia el punt o de origen o par t ida.
C lculo de la Ru t a m s cor t a : Una vez obt enidas las superficies de coste de desplazam ient o y los ot ros parm et ros (direccin y asignacin) necesarios para el clculo de la ruta m s cort a, podrem os calcular la ruta m s corta ent re los puntos de origen y dest ino basndonos en los parm et ros ant eriorm ent e calculados. La funcin la podem os localizar en Spa t ia l Ana lyst D ist a nce Shor t est Pa t h , y los parm etros a especificar son los siguient es: Pat h To ( Rut a hacia) : aqu especificam os el punt o de dest ino donde querem os llegar. Cost D ist a n ce Rast e r ( Rast er de Dist ancia de Cost o) : raster que obt uvim os en el clculo anterior. Cost Dir ect ion Ra st er ( Rast er de Direccin) : rast er de direcciones que se usar para el clculo de la rut a. Pat h Type : aqu podrem os seleccionar si querem os que el clculo se haga para cada zona de asignacin, para cada celda, o la m ej or. Out p ut Fea t ures : Aqu obtendrem os la ruta o rut as finales en form at o Shapefile o Feat ure Class.
Pgina 47
Fig. 56. Localizacin de la herram ient a para la obtencin de la rut a m s cort a ent re dos punt os dados y los parm et ros de configuracin.
Fig. 57. Rut a obtenida ent re el punt o de inicio y el punt o final ponderada por la pendient e o superficie de coste.
Opcionalm ent e podem os obt ener las zonas de asignacin t al y com o hicim os ant eriorm ent e respect o a una capa de punt os o ent idades de referencia, pero de form a individual y no integrada en el clculo de dist ancias para ello vam os a Spa t ial An a lyst D ist an ce - Allocat ion .
Pgina 48
Los parm et ros a especificar son los siguientes: I nput D a t a ( Datos de Ent rada) : Popu la t ion Fie ld (Cam po de dat os) : dat os a usar para evaluar su grado de concent racin o dispersin. D e nsity Typ e ( Tipo de anlisis) : m ediant e un ncleo o filt r o, o sim ple. Est ando la diferencia principal ent re am bos m t odos en que el Kernel perm it e valorar m s los punt os m s prxim os al ncleo del Kernel que los m s alej ados, m ient ras que para el sim ple, se evalan t odos los elem ent os del m ism o m odo dent ro del rea de anlisis, dndonos el m t odo Kernel superficies m s suavizadas. Se ar ch Radius ( Radio de Bsqueda) : Radio de anlisis en las unidades especificadas para el anlisis. Ar e a Unit s (Unidades de rea) : unidades de m edicin del rea de m uest reo. Out put Ce ll Siz e ( Tam ao de celda de salida) : t am ao que tendr el raster de salida, por defect o el m ism o del MDE que podam os usar posteriorm ente. Out p ut Rast er ( Rast er de Salida) : fichero rast er de form at o ESRI Grid que contendr los dat os finales de densidad.
En la siguient e figura podem os observar un ej em plo donde hem os obt enido la densidad de puntos de cota por rea de anlisis. Com o se puede observar exist e una m ayor concent racin de punt os en la zona noreste del t rm ino m unicipal, dndose casos de escasez de punt os en la zona suroest e, por lo que sera recom endable obt ener alguna fuent e adicional para las zonas que disponen de pocos punt os acotados que puedan det ectar las form as del relieve en dicha zona.
Pgina 49
Fig. 59. Mapa de densidad de punt os de cot a usados para gener ar el MDE y verificar las zonas con escasa densidad de punt os.
Pgina 50
Fig. 60. Mt odos de I nterpolacin Det erm inst icos. Mt odo I DW y parm et ros de configuracin. El resultado podem os observarlo a cont inuacin:
Fig. 61. I nterpolacin result ant e m ediant e el m t odo I DW. I n t e rpola cin Splin e ( Regularizado o de I sotensin) : Este m t odo est im a los valores usando una funcin m atem t ica que reduce al m nim o la curvat ura de la superficie tot al, dando com o result ado una superficie que pasa a t ravs de los puntos m uest reados. Hay dos m t odos Splines: regularizado o t ensin. - Regu lar iza do. Este m todo crea una super ficie suave, de m aner a gr adual com binando la super ficie con los valores que pueden est ar fuera de la gam a de dat os de la m uest ra. - Te nsin. El m t odo de tensin suaviza la rigidez de la superficie segn el carcter del fenm eno m odelado. Esto crea una superficie lisa con valores est recham ent e ligados a los valores de la m uest ra.
Pgina 51
Fig. 62. I nterpolacin m ediante Spline y sus parm et ros. Dej ando los parm et ros por defect o que nos facilit a el soft ware podem os obt ener una superficie com o la que se m uest ra a cont inuacin:
Fig. 63. Result ado de la I nt er polacin m ediant e Spline. I n t e rpola cin Kriging ( Universal u Ordinario) : est e m t odo de interpolacin global t iene en cuenta un variogram a para la obt encin de la superficie final, Para ello se resuelve un conj unt o de ecuaciones ut ilizando la inform acin procedente del v a riogra m a y las dist ancias relat ivas ent re los dat os y la posicin del punt o donde el valor int erpolado es solicit ado. Est as ecuaciones cont ienen la cov a rian z a ent re el punt o a ser est im ado y los datos y las covarianzas ent re los dat os en s m ism os. La herram ienta se encuent ra en Spa t ial An alyst I n t er polat e To Rast er Kr iging , y donde podem os seleccionar los siguientes parm et ros:
Pgina 52
Fig. 64. Localizacion de la herram ient a Kriging y los parm et ros a especificar.
Fig. 65. I m agen del raster procedent e de la funcin de int erpolacin m ediante Kriging.
Pgina 53
Para realizar est e t ipo de anlisis debem os disponer com o m nim o dos bandas o capas raster referidas a un m ism o territ orio y georreferenciadas en un m ism o sistem a de referencia espacial, para poder ext raer los valores est adst icos de un pxel respect o a todas las bandas im plicadas. La her ram ient a se puede localizar en Spa t ial An a lyst Ce ll St a t ist ics , y los parm et ros requeridos son los siguientes:
Fig. 66. Localizacin de la herram ient a de Est adst icas basadas en Celdas o Pxeles y sus parm et r os de anlisis.
Los parm et ros son los siguient es: La ye r s ( Capas) : Capas disponibles para seleccionarlas e im plicarlas en el anlisis. I nput Ra st ers ( Rasters de Ent rada) : capas usadas para el anlisis. Ove r lay St a t ist ics ( Estadst icas de Solape) : funcin estadst ica que calcularem os para cada pixel basndonos en su localizacin en t odas las bandas o capas rast er de ent rada. Out p ut Rast er ( Raster de Salida) : raster que alm acenar los dat os est adst icos resultantes. Ent re los valores est adst icos que se pueden obtener t enem os: Maj orit y, m axim um , m ean, m edian, m inim um , m inorit y, range, st andard deviat ion, sum , variet y. As por ejem plo si disponem os de una gam a de rast er de t em perat ura o precipit aciones para t odo un ao dividido por m eses, podram os obtener valores de variabilidad, concent racin y dispersin con respecto al valor m edio, et c. Post eriorm ent e se podrn realizar correlaciones ent re raster, redundando en un anlisis clim t ico m s det allado. En el ej em plo aplicado hem os calculado el valor m edio de t odos los pxeles de una im agen Landsat usando las bandas con com binacin RGB ( 321 o Color Natural) .
Fig. 67. Result ado de la funcin de resum en estadst ico local ( Media de los valores iniciales) .
Pgina 54
Los parm et ros de configuracin son los siguient es: I nput D a t a ( Datos de Ent rada) : Raster a analizar m ediante la funcin. Fie ld ( Cam po) : Cam po de Referencia (al ser Raster Value). St a t ist ic Type ( Tipo de Estadst ico) : t ipo de estadst icas a obt ener. N e ighbor hood ( Crit erio de Proxim idad) : figura geom t rica de referencia. N e igh bor hood Se t t ings ( Parm et ros de Proxim idad) : son los parm et ros que se especifican post eriorm ente para obtener los valores de los pxeles o las reas (en unidades m t ricas) m s prxim as. Ou t pu t Cell Size ( Tam ao de celda de salida) : ser el t am ao del raster por defect o. Out p ut Rast er ( Rast er de Salida) : aqu escribim os el ESRI Grid de salida que alm acenar los valores de salida de las est adst icas focales.
Fig. 68. Localizacin de la herram ient a Neight borhood St at ist ics y sus parm et r os de uso.
Las form as geom t ricas que podem os usar para el clculo son:
Fig. 69. Form as geom t ricas de anlisis de vecindad. As para las dist int as figuras t enem os dist int os parm et ros de configuracin que nos perm it irn personalizar la salida raster:
Pgina 55
Rect a ngle ( Rectngulo) : Anchura y Alt ura. Cir cle ( Crculo) : Radio. Annulus ( Anillo) : Radio int erior y exterior del anillo. W edge (Cua) : ngulo inicial, ngulo final y radio. Las unidades se podrn especificar en nm ero de pxeles o en unidades de m apa. Com o result ado de la aplicacin podem os obt ener parm et ros que nos perm it an ver el grado de var iabilidad a pequea escala ent re un punt o y su rea m s prxim a, pudiendo det erm inar zonas con un elevado gradient e o var iabilidad de una det erm inada var iable ( niveles de cont am inacin, focos t rm icos, et c.)
Fig. 70. Localizacin de la herram ient a Zonal Stat ist ics y parm et ros de uso.
Pgina 56
I nput Ra st er ( Rast er de Ent rada) : rast er a r eclasificar ( por ej em plo de pendient es) Re class Field (Cam po) : cam po de base para la reclasificacin. Se t Value s t o Re classify : aqu podem os especificar los valores que se usarn para la reclasificacin y que post eriorm ent e poblarn t ant o el rast er reclasificado, com o el vect orial de salida ( m ediante el cam po GRI DCODE) en caso de convert ir a vect orial. Ent re las opciones que podem os usar dent ro de est e grupo de elem ent os t enem os Classify, si no hem os usado est e m t odo previam ent e para clasificar los datos cont inuos en intervalos; Unique, para establecer valores nicos com o clases del raster, es decir, un valor por cada valor del pixel; Add Ent ry, para agregar nuevas clases; Delete Ent ries, para borrar clases asignadas ant eriorm ente; Precision, parm et ro de configuracin usado para la divisin de las clases y que int erviene en el nm ero de elem ent os que ent ra en una u ot ra clase. Si disponem os de una tabla LUT ( Look Up Table) con la configuracin podr em os cargarla Load o Salvar Save nuest ra tabla existente para ot ro anlisis sim ilar. Ch ang e M issing v alues N oD a t a : est a opcin nos perm it e cam biar los valores sin asignacin de nuest ro raster a valores NoData. Out p ut Rast er ( Rast er de Salida) : el nom bre del raster en form at o ESRI Grid que contendr las clases tem t icas.
Fig. 71. Parm et ros y localizacin de la Her ram ient a de Reclasificacin Rast er.
Pgina 57
Fig. 72. Result ado de la reclasificacin del m apa de pendient es en 5 int ervalos.
Pgina 58
Fig. 73. Localizacin del m en de conversin y t ipos de conversin. Parm et ros de conversin Rast er Vect or y Vect or Rast er.
Aqu podem os observar com o ejem plo la conversin de rast er a vect or de las cuencas visuales:
Fig. 74. Conversin Rast er a Vect or de las Cuencas Visuales. Ot ro grupo de herram ient as nos perm it en realizar anlisis m s especficos com o son las herram ientas de Anlisis Hidrolgicos que nos provee dicha extensin y que puede usarse m ediant e ArcToolbox, Model Builder, Lnea de com andos o la Calculadora Raster, as com o Anlisis de Aguas Subterrneas, ent re ot ros. Este t ipo de procesos sern am pliados en una segunda edicin de esta gua.
Pgina 59
Fig. 75. Her ram ient as de Geoprocesam ient o de ArcToolbox Spat ial Analyst Tools. Ent r e est as herr am ient as t enem os los siguient es grupos: lgebr a de m a pa s : com puest a por los grupos Condit ional, Map Algebra, Mat h, Overlay, et c. H e rr a m ient a s de Gene ra liza cin : Gener alizat ion. H e r r a m ien t a s de e x t r accin de in for m acin : Ext r act ion, Densit y , Dist ance, Local, Zonal, Neighbor hood, et c. H e rr a m ient a s a plica da s : Hydrology, Groundwat er, Surface, Solar Radiation, et c.
De t odas est as herram ient as las que m ayorm ent e usarem os desde ArcToolbox sern las de Ext r accin de ent idades de una superficie en form at o rast er ( enm ascarado o ext r accin de un rast er por m scara) o vect orial ( ext raccin por punt os por ej em plo) . En el caso de la ext raccin del Rast er a part ir de una m scara podem os encont rar dicha herram ient a en Spat ial Analyst Tools Ext ract ionExt ract By Mask, donde nos lim it arem os a especificar el rast er de ent rada, la m scara ( bien sea rast er o vect orial) y un fichero rast er de salida. Debido a la ex t ensin que conllevara com ent ar t odas las herr am ient as de Anlisis Espacial Rast er desde ArcToolbox, dej arem os est e apart ado para la am pliacin de est a gua que se act ualizar prxim am ent e. M odeliz a cin de An lisis Espa cia l m edia nt e M odel Bu ilder . Com o hem os com ent ado las herram ient as de anlisis rast er nos perm it en m odelar la superficie t errest re en base, principalm ent e, a dat os raster.
Pgina 60
Fig. 76. Funcin prim ar ia de geopr ocesam ient o. Vam os a realizar un m odelo de geoprocesam ient o que nos va a perm it ir m ost rar gran part e de la funcionalidad de dicha herram ient a. En est e m odelo usarem os un MDE para generar a part ir de st e, las superficies de pendient e. A part ir de st os usando una operacin de lgebr a de m apas obt endrem os las superficies que cum plan con la condicin siguient e: Crit erios: Elevaciones por encim a de los 200 Met ros y cuya pendient e sea superior al 30% . Lo prim ero que vam os a hacer es crearnos una Toolbox en ArcToolbox, para ello bast a con pulsar bot n derecho del r at n sobr e la cabecer a de ArcToolbox y seleccion ar New Toolbox. La renom bram os com o pr efiram os. En su int erior podrem os agregar un Modelo de Geoprocesam ient o, pulsarem os para ello bot n derecho del rat n New Model. A cont inuacin se puede observar el m odelo pr opuest o.
Fig. 77. Creacin de la herram ient a de geoprocesam ient o con Model Builder.
Pgina 61
Est os elem ent os podr em os r enom br arlos pulsando bot n derecho del rat n sobre el obj et o Renam e ( Renom br ar) . Una vez que disponem os de st os en el m odelo vam os a enlazarlos. Para ello pulsarem os el bot n Add Conn e ct ion y unirem os cada elem ent o con su funcin de anlisis de la barra de herram ient as pr incipal del m odelo.
Herramienta de Conexin
Fig. 78. Herr am ient a de conexin de objet os y Modelo definit ivo ant es de su ej ecucin. Lo nico que t endrem os que especificar es la frm ula en la funcin de lgebra de m apas:
( mde > 200) AND ( raster pendiente > 30)
Fig. 79. Adicin de la frm ula de lgebra de m apas para el clculo especifico.
Pgina 62
Fig.80. Result ado del proceso del m odelo de geopr ocesam ient o.
Pgina 63
Una vez realizado el anlisis procederem os a int erpret ar los result ados y a present ar debidam ent e la infor m acin obt enida, m ediant e m apas, est adst icas, inform es, etc.
Con est e epgr afe llegam os al final de est a gua m et odolgica de Anlisis Espacial m ediant e herr am ient as de anlisis Rast er con Spat ial Analy st para ArcGI S. Recordad que est a gua est en const ant e evolucin, y act ualizacin por lo que en breve habr ot ra versin m s act ualizada con
Pgina 64
Gr acias nuevam ent e Gabr iel Or tiz. ht t p: / / ht t p: / / www.gabrielort iz.com ht t p: / / foro.gabrielort iz.com
Referencias recom endadas: Gut irrez Puebla, Javier; Gould Michael ( 2000) . SI G: Sistem a s de I nform a cin Geogr fica . Ed. Snt esis. Prim er a r eim presin. ESRI (2003) . Using Spa t ia l Ana lyst . Envir onm ent al Syst em Resear ch I nstit ut e. Moreno Jim nez, A. ( 2005, coord.) : Sist e m as y a n lisis de la in for m a cin geogr fica . M a nua l de a u t oa pre n diza j e con ArcGI S. Edit orial Ra-Ma. Pea Llopis, Juan ( 2006) : Sist e m as de I n for m a cin Geogr fica a plica dos a la gest in de l t e r r it or io . Depart am ent o de Ecologa. Univ ersidad de Alicant e.
Not a: Est a gua ha sido realizada y sus im genes se han obt enido con Licencia Original de ArcGI S Deskt op ( Ar cView Nodo Localizado) y ver sin dem o de las Ext ension es por un periodo lim it ado.
Pgina 65