<div dir="ltr">hi Veronica,<br><br>On Thu, May 1, 2014 at 3:17 PM, Veronica Andreo <<a href="mailto:veroandreo@gmail.com">veroandreo@gmail.com</a>> wrote:<br>> Dear all,<br>><br>> I'm using the following command to run r.hants on a series of 506 maps of<br>
> chlorophyll. It seems to start thinking and a few seconds after i get an<br>> error:<br>><br>> r.hants file=lista_hants nf=6 fet=0.1 dod=5 delta=0.1 base_period=46<br>><br>> ERROR: no temp files available<br>
<br>probably it runs out of memory? How many pixels does one map have?<br><br>> As stated before, the total number of images is 506 (11 years) but the base<br>> period is 1 year, hence base_period=46 images, and i'm passing a plain text<br>
> file with the map names (one per line) as stated here:<br>> <a href="http://grass.osgeo.org/grass70/manuals/addons/r.hants.html">http://grass.osgeo.org/grass70/manuals/addons/r.hants.html</a><br><br>You can simply use g.mlist for that, much easier (see below).<br>
<br>> The files are registered in a spatio-temporal raster data set, if that means<br>> something, dunno.<br><br>I guess that this will be simply ignored, so no problem.<br><br>-------------<br><br>I have tried myself, using the new dataset prepared by Anna:<br>
<br>Climate dataset for North Carolina from publicly available data, with <br>monthly temperature average and monthly precipitation from 2000 to 2012:<br><a href="http://courses.ncsu.edu/mea592/common/media/02/nc_climate_spm_2000_2012.zip">http://courses.ncsu.edu/mea592/common/media/02/nc_climate_spm_2000_2012.zip</a><br>
<br><font face="courier new, monospace">#########<br><br>GRASS 7.1.svn (nc_climate_spm_2000_2012):~ ></font><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace"># set computational region</font></div>
<div><font face="courier new, monospace">g.region rast=2004_03_tempmean -p<br><br></font></div><div><font face="courier new, monospace"># check for the maps, especially how many</font></div><div><font face="courier new, monospace">g.mlist rast pattern=20??_??_tempmean | wc -l<br>
raster map(s) available in mapset <climate_1970_2012>:<br>156<br><br># try out at lower resolution for quick assessment:<br>g.region -p res=5000</font></div><div><font face="courier new, monospace"><br></font></div>
<div><font face="courier new, monospace"># HANTS</font></div><div><font face="courier new, monospace"># just wildly guessing the parameters for a test run:<br>r.hants input=`g.mlist rast pattern=20??_??_tempmean sep=comma` \</font></div>
<div><font face="courier new, monospace"> nf=6 fet=0.1 dod=5 delta=0.1 base_period=12<br>Harmonic analysis of 156 input maps...<br> 100%<br><br># assign reasonable color tables<br>for map in `g.mlist rast pattern="*tempmean_hants"` ; do</font></div>
<div><font face="courier new, monospace"> r.colors $map color=celsius</font></div><div><font face="courier new, monospace">done<br><br># verify with one of the 156 results:<br># remember that we had reduced the resolution</font></div>
<div><font face="courier new, monospace"># by factor 10 to this quick assessment<br>r.mapcalc "diff = 2000_06_tempmean - 2000_06_tempmean_hants"</font></div><div><font face="courier new, monospace"><br>r.colors diff color=differences</font></div>
<div><font face="courier new, monospace">d.mon wx0</font></div><div><font face="courier new, monospace">d.rast diff</font></div><div><font face="courier new, monospace"><br>r.univar diff -g<br>n=5066<br>null_cells=5040<br>
cells=10106<br>min=-0.0899336115228095<br>max=0.359362050140941<br>range=0.449295661663751<br>mean=0.188579838052468<br>...<br><br># see original as animation (still at 5km resolution)<br># note: color table is different from standard "celsius"<br>
g.gui.animation rast=`g.mlist rast pattern="*tempmean" sep=comma`<br><br># see HANTS time series as animation<br>g.gui.animation rast=`g.mlist rast pattern="*tempmean_hants" sep=comma`</font><br></div>
<div><br></div><div>Looks all fine!</div><div><br></div><div>Best</div><div>markusN</div></div>