[Gvsig_italian] Aspect - gvSIG e SEXTANTE

[Claudio]

Grazie Antonio,
era proprio ciò che desideravo sapere, ovvero come sextante  
restituisse l'esposizione rispetto al classico modello in gradi da 0 a  
360.
Ho provato sia manualmente in gvSIG 1.1.2 sia con il tuo modello in  
1.9 e funziona tutto perfettamente.
Spero che la versione 2.0 esca al più presto così svolgere tutte le  
vecchie esercitazioni in ambiente OS di gvSIG.
ciao
Claudio


Il giorno 09/dic/08, alle ore 18:39, Antonio Falciano ha scritto:

> Claudio ha scritto:
>> Buongiorno a tutti, sono Claudio e svolgo (attualmente) ricerca alo
>> IUAV sulle nuove tecnologie per il territorio e l'ambiente.
>> Mi sono avvicinato al mondo Open Source GIS perchè vorrei portarlo  
>> nei
>> laboratori GIS del mondo universitario.
>
> Benvenuto!
>
>> Dopo questa breve presentazione, vi pongo il primo quesito riguardo
>> gvSIG e SEXTANTE:
>> - ho provato ad utilizzare l'algoritmo Aspect ma non riesco a  
>> venire a
>> capo sulla attribuzione dei valori di cella rispetto alle diverse
>> esposizioni.
>> Qualcuno di voi si è già imbattuto in questo problema riuscendo a
>> codificare per le classi di valori la rispettiva esposizione dei
>> versanti?
>
> L'aspect, come è noto, esprime la direzione principale di pendenza  
> di ogni
> cella. E' generalmente espressa in gradi positivi da 0 a 360 misurati
> in senso orario a partire dal Nord. In Sextante è invece espressa in
> radianti ed è pertanto compresa tra 0 e 2 pi greco (PI).
> Per ottenere l'esposizione dei versanti, basterà quindi  
> riclassificare i
> valori ottenuti suddividendo opportunamente 2 PI in base al numero di
> classi di interesse. Ad es. se sei interessato alle canoniche 8 classi
> (N, NE, E, SE, S, SW, W, NW), occorre suddividere 2 PI in 8 angoli di
> uguale ampiezza (PI/4), ottenendo quindi:
>
> [0, PI/8) U [15/8 PI, 2 PI] --> 0 (N);
> [PI/8, 3/8 PI)                  --> 1 (NE);
> ...
> [11/8 PI, 13/8 PI)            --> 6 (W);
> [13/8 PI, 15/8 PI)            --> 7 (NW).
>
> Ho creato un piccolo modello in Sextante (reclassify_aspect.model in
> allegato) con la 1.9 alpha, in modo da non dover ogni volta ridefinire
> la tabella contenente i valori delle classi.
> Per visualizzare adeguatamente i risultati ottenuti, occorre poi
> definire ed associare una tabella dei colori al raster risultante.
>
> ciao
> Antonio
>
> <?xml version='1.0' encoding='ISO-8859-1' standalone='yes' ?>
> 	<model name="Reclassify aspect">
> 		<algorithm key="PROC0" class="class  
> es.unex.sextante.morphometry.aspect.AspectAlgorithm"  
> coords="85,98,57,58">
> 			<assignment key="METHOD" assigned_to="INNERPARAM0PROC0" />
> 			<assignment key="DEM" assigned_to="INPUT0" />
> 		</algorithm>
> 		<algorithm key="PROC1" class="class  
> es.unex.sextante.gridCategorical.reclassify.ReclassifyAlgorithm"  
> coords="162,176,72,58">
> 			<assignment key="INPUT" assigned_to="ASPECTPROC0" />
> 			<assignment key="LUT" assigned_to="INNERPARAM1PROC1" />
> 			<assignment key="METHOD" assigned_to="INNERPARAM0PROC1" />
> 		</algorithm>
> 		<input name="INPUT0" description="DTM" class="class  
> es.unex.sextante.parameters.ParameterRasterLayer"  
> coords="20,20,45,58">
> 			<attribute name="mandatory" value="true" />
> 		</input>
> 		<output key="RECLASSPROC1" class="class  
> es.unex.sextante.outputs.OutputRasterLayer"  
> description="Aspect_reclass" />
> 		<data_object key="INNERPARAM0PROC1" value="1"  
> description="Selection" />
> 		<data_object key="INNERPARAM0PROC0" value="0"  
> description="Selection" />
> 		<data_object key="INNERPARAM1PROC1" value="[[0.0|0.393|0.0],[0.393| 
> 1.178|1.0],[1.178|1.963|2.0],[1.963|2.749|3.0],[2.749|3.534|4.0], 
> [3.534|4.32|5.0],[4.32|5.105|6.0],[5.105|5.89|7.0],[5.89|6.283| 
> 0.0]]" description="Fixed Table" />
> 		<data_object key="ASPECTPROC0" value="class  
> es.unex.sextante.gui.modeler.elements.ModelElementRasterLayer,1"  
> description='"Aspect" / Process 0: Aspect' />
> 		<data_object key="RECLASSPROC1" value="class  
> es.unex.sextante.gui.modeler.elements.ModelElementRasterLayer,1"  
> description='"Reclassify" / Reclassify' />
> 	</model>
> <?xml version='1.0' encoding='ISO-8859-1' standalone='yes' ? 
> >_______________________________________________
> Gvsig_italian mailing list
> Gvsig_italian a lists.gvsig.org
> https://gvsig.org/lists/mailman/listinfo/gvsig_italian

Arch. Claudio Schifani
Phd in Pianificazione Urbana e Territoriale
Dipartimento Città e Territorio
via dei Cartari 19b,
90138 Palermo
tel.3398429091