Re: [R-sig-Geo] Question on histograms from Raster and RasterVIS packages

2014-11-19 Thread Nuno
. If that happens, then you can create the data for a histogram yourself by using the function 'freq', perhaps after 'reclassify' (and then probably 'barplot'). Robert On Tue, Nov 18, 2014 at 9:56 AM, Nuno Sá nunocesard...@gmail.com wrote: Hello! My aim is to add a Median or a Mean line

[R-sig-Geo] Question on histograms from Raster and RasterVIS packages

2014-11-18 Thread Nuno
Hello! My aim is to add a Median or a Mean line to an histogram plot in R. The problem is the following: The hist function from the raster package uses a maximum of 100 000 values for generating the histogram but allows me to edit the plot, so I can easily add a line within the plot using

Re: [R-sig-Geo] Question on histograms from Raster and RasterVIS packages

2014-11-18 Thread Nuno
: require(stats) set.seed(14) x - rchisq(100, df = 4) hist(x) meanx-mean(x) meanx [1] 3.734162 abline(v=meanx, col=red) Cheers! Tom On Tue, Nov 18, 2014 at 10:56 AM, Nuno Sá nunocesard...@gmail.com wrote: Hello! My aim is to add a Median or a Mean line to an histogram

Re: [R-sig-Geo] Question on histograms from Raster and RasterVIS packages

2014-11-18 Thread Nuno
in more detail, probably using the grid.lines function. There is also the gridBase package which helps grid and base graphics play nicely together, though it could be a bit of overkill for just adding a vertical line. On Tue, Nov 18, 2014 at 11:21 AM, Nuno Sá nunocesard...@gmail.com wrote: Hey

[R-sig-Geo] Landsat ETM+ Gap filling WELD data

2014-10-15 Thread Nuno
Hey So, the first question is - is there an implementation of the newer algorithms for Gap filling in any R package? If not, I would be available to work with an expert to try to develop it here - some of the actual IDL code is available through Dr. Xiaolin Zhu's page:

Re: [R-sig-Geo] Landsat ETM+ Gap filling WELD data

2014-10-15 Thread Nuno
Version 1.5 Product Documentation, subsection Product Map Projections, they use an Albers Equal Area projection. Thanks, Eric On Wed, Oct 15, 2014 at 5:16 AM, Nuno Sá nunocesard...@gmail.com wrote: Hey So, the first question is - is there an implementation of the newer algorithms for Gap

Re: [R-sig-Geo] Landsat ETM+ Gap filling WELD data

2014-10-15 Thread Nuno
and address them individually. It's too much to ask a general question about doing that, and there is a lot you could do first. Cheers, Mike. On Thu, Oct 16, 2014 at 2:29 AM, Nuno Sá nunocesard...@gmail.com wrote: Hello Eric Thank you, I had seen it though. Are you aware if the rgdal able

Re: [R-sig-Geo] Mann-kendall

2014-10-09 Thread Nuno
Change Institute University of Florida www.clas.ufl.edu/users/forrest On Wed, Oct 8, 2014 at 7:58 PM, Nuno Sá nunocesard...@gmail.com wrote: Hello! ah! that return(unlist(MannKendall(x))) might have been why i didn't manage to do it with the calc function, will try it out, thank you

[R-sig-Geo] Mann-kendall

2014-10-08 Thread Nuno
? Am I forced to do series expansions? would this not be to slow with raster? Or is there a quicker way to iterate pixel by pixel? Sorry for the long email and thank you in advance for any comments/corrections/recommendations! Best regards to all, Nuno Sá -- Nuno César de Sá +351 91 961 90 37

Re: [R-sig-Geo] Mann-kendall

2014-10-08 Thread Nuno
(2*n+1) } phi.aprox - phi-aprox*phi.integ p.value.rst.neeti - 2*(1-phi.aprox) --- The result falls within 0 to 1 which is a good indication. But all depends on everything I've asked before. Thank you! On 8 October 2014 22:23, Nuno Sá nunocesard...@gmail.com wrote: Hello

Re: [R-sig-Geo] Mann-kendall

2014-10-08 Thread Nuno
and Environmental Change Institute University of Florida www.clas.ufl.edu/users/forrest On Wed, Oct 8, 2014 at 6:02 PM, Nuno Sá nunocesard...@gmail.com wrote: Hello Using the integration by part summation was quite fast (I owe a beer to the ones who've developed the raster package