Agus, That is indeed because of an unintended difference between the
released version (Barry's) and the development version on R-Forge (yours). I
have fixed the development version. Thanks for pointing this out. Robert
On Fri, Jul 15, 2011 at 12:12 PM, Agustin Lobo wrote:
> Barry,
>
> I tried tha
Barry,
I tried that first, but by some reason (my raster version?), the
result here is different than yours:
> CODEQP <- raster(extent(PNMlimit))
> show(CODEQP)
class : RasterLayer
dimensions : 10, 10, 100 (nrow, ncol, ncell)
resolution : 2586.382, 1886.275 (x, y)
extent : 438048.4,
# You could do this:
dem100 <- resample(dem, CODEQP, method='bilinear')
# But as the resolutions are rather different, I think it is better to first
aggregate such that they are about the same:
temp <- aggregate(dem, 10)
# And then use resample:
dem100 <- resample(temp, CODEQP, method='biline
>> I've a raster stack with 3 layers.
>>
>> I need to classify it assigning one class to each unique combination of
>> values.
> Perhaps the below works. It may be very slow if unq has many rows.
That was not quite right, I think it works when you do:
unq <- unique(as.matrix(data_stack))
Robert,
While you're re-opening the code, I'd like to ask about the col parameter. I
personally had trouble with setting the right range for my palette as zlim
is overwritten. For instance I have a raster which range isn't symmetric
around 0, but I have a palette which I'd like to be symmetric aro
On Fri, Jul 15, 2011 at 11:10 AM, Agustin Lobo wrote:
> I have 2 raster objects
>
>> show(CODEQP)
> class : RasterLayer
> dimensions : 10, 10, 100 (nrow, ncol, ncell)
> resolution : 100, 100 (x, y)
> extent : 438048.4, 439048.4, 4633415, 4634415 (xmin, xmax, ymin, ymax)
> coord. re
On Fri, Jul 15, 2011 at 12:00 PM, Agustin Lobo wrote:
> Hi!
>
> Is it possible to create a raster object using extent and resolution,
> thus having
> raster() to calculate the required nb. of rows and columns? Something like,
>
> a = raster (ext=miextent, res=100)
>
> (Obviously the nb. of pixels
Hi!
Is it possible to create a raster object using extent and resolution,
thus having
raster() to calculate the required nb. of rows and columns? Something like,
a = raster (ext=miextent, res=100)
(Obviously the nb. of pixels would have to be rounded)
Thanks!
Agus
I have 2 raster objects
> show(CODEQP)
class : RasterLayer
dimensions : 10, 10, 100 (nrow, ncol, ncell)
resolution : 100, 100 (x, y)
extent : 438048.4, 439048.4, 4633415, 4634415 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=utm +zone=31 +ellps=intl +units=m +no_defs
values : i
to remove duplicates in DF you need to replace it:
DF = remove.duplicates(DF, zero=0, remove.second = TRUE)
zerodist(DF)
On 07/15/2011 10:23 AM, Matevž Pavlič wrote:
> Hi, thanks for your time
>
> I am sure the problem is in same locations...but i dont know why
> remove.duplicates() does no
Dear Robert
i have tested using our example and it seems to work fine on mac
Best
On 15 Jul 2011, at 09:51, lancelot wrote:
> Dear Robert,
>
> Working with:
>
> > sessionInfo()
> R version 2.13.1 (2011-07-08)
> Platform: i386-pc-mingw32/i386 (32-bit)
>
> locale:
> [1] LC_COLLATE=French_France.
Hi, thanks for your time
I am sure the problem is in same locations...but i dont know why
remove.duplicates() does not remove duplicates...
When i try this :
>coordinates(DF)<-~X+Y
> zerodist(DF)
[,1] [,2]
[1,] 83 84
[2,] 242 243
[3,] 266 267
[4,] 431 432
[5,] 439 440
[6,]
Thank you Ashton for your response. I mistyped in the email; it wasn't a
problem with the capitalization.
Robert, I closed everything out and began anew with your first set of
instruction. That worked perfectly.
Thank you so much for your help!
Best,
Rebecca
--
View this message in context:
Dear Robert,
Working with:
> sessionInfo()
R version 2.13.1 (2011-07-08)
Platform: i386-pc-mingw32/i386 (32-bit)
locale:
[1] LC_COLLATE=French_France.1252 LC_CTYPE=French_France.1252
[3] LC_MONETARY=French_France.1252 LC_NUMERIC=C
[5] LC_TIME=French_France.1252
attached base packages:
[1] sta
14 matches
Mail list logo