[GRASS-user] i.pca with scrips

Markus Metz markus.metz.giswork at googlemail.com
Mon Jul 4 14:49:15 EDT 2011


On Mon, Jul 4, 2011 at 8:43 PM, Gaspar <gos47 at hotmail.com> wrote:
> Dear Grass users
>
> I am trying to run the principal components analysis with "i.pca" for 12
> raster.
>
> I made a script to process three times the "i.pca" with three sets
> (twelve raster by each  sets) of different parameters (temperature,
> chlorophyll, Insolation), but when I run the scripts this only takes 11
> images by sets. I tried i.pc manually and this takes twelve raster
> correctly. Thank you in advance.
>
> EXAMPLES WHIT SCRIP:
>
> GRASS 6.4.1 > i.pca_.sh
> .
> .
> .
> Computing means for band 11...
>  100%
> Computing row 1 (of 11) of covariance matrix...
>
> EXAMPLES MANUALLY FOR ONE SET:
>
> GRASS 6.4.1 >
> input=tmedia_1,tmedia_2,tmedia_3,tmedia_4,tmedia_5,tmedia_6,tmedia_7,tmedia_8,tmedia_9,tmedia_10,tmedia_11,tmedia_12 output=PCA_tmedia_
>
> Computing means for band 12...
>  100%
> Computing row 1 (of 12) of covariance matrix...
>
>
> EXAMPLE SCRIPS
>
> #!/bin/bash
>
> # Analisis de componentes principales
>
> for map in tmedia_ aston_ clor_
> do
>
> B1=1
> B2=2
> B3=3
> B4=4
> B5=5
> B5=6
> B7=7
> B8=8
> B9=9
> B10=10
> B11=11
> B12=12
>
> i.pca rescale=0,0 input="$map$B1,$map$B2,$map$B3,$map$B4,$map$B5,$map
> $B7,$map$B8,$map$B9,$map$B10,$map$B11,$map$B12" output=PCA_$map
>

Is $map$B6 really missing in the script or is it just a typo above?


More information about the grass-user mailing list