<div dir="ltr">Hi Shiva,<br><br><div class="gmail_quote"><div dir="ltr">El jue., 12 jul. 2018 a las 2:45, Shiva Khanal (<<a href="mailto:shiva_khanal@hotmail.com">shiva_khanal@hotmail.com</a>>) escribió:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">




<div dir="ltr">
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="font-family:Calibri,sans-serif;font-size:11pt">Hi,</span></div>
<div>
<div class="gmail-m_1725168727447190363WordSection1">
<p class="MsoNormal" style="margin:0in 0in 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">
 </p>
<p class="MsoNormal" style="margin:0in 0in 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">
I am interested to run the spatial-mode PCA on raster time series. In R it is implemented as an ‘EOF’ function in spacetime Package. The function transposes the raster stack and applies ‘prcomp’ function.  I would like to know if there is any possible way to
 implement i.pca similar to this approach.</p>
<p class="MsoNormal" style="margin:0in 0in 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">
 </p>
<p class="MsoNormal" style="margin:0in 0in 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">
Interestingly, I found a very good tutorial on GRASS raster time series at <a href="https://grasswiki.osgeo.org/wiki/Temporal_data_processing/GRASS_R_raster_time_series_processing" target="_blank">
https://grasswiki.osgeo.org/wiki/Temporal_data_processing/GRASS_R_raster_time_series_processing</a> This includes the example using spacetime package to read GRASS raster time series into R.</p>
<p class="MsoNormal" style="margin:0in 0in 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">
</p></div></div></div></blockquote><div><br></div><div>there's also an alternative way to export raster time series from GRASS that was not there at the moment of creating the wiki: the add-on <a href="http://t.rast.out.xyz">t.rast.out.xyz</a> [0] that produces a table with xy coordinates and all dates in columns. You can then import that into R, transpose and run prcomp, EOF or others. Do you want it for gap-filling?<br></div><div><br></div><div>If you dare to implement EOF in GRASS, please go ahead!! I will be very happy to test and provide feedback, I've been wishing for such a thing (and also DINEOF) for a while :) <br></div><div><br></div><div>Maybe the devs can provide hints and advice on how to address it.  <br></div><div><br></div><div>Best, <br></div><div>Vero</div><div><br></div><div>[0] <a href="https://grass.osgeo.org/grass74/manuals/addons/t.rast.out.xyz.html">https://grass.osgeo.org/grass74/manuals/addons/t.rast.out.xyz.html</a></div><br></div></div>