A couple of hints:
asc files are very slow to process in R. I would write it back to the disk
as either a native grd, or even geoTiff, before reclassifying (in my
experience, perhaps a 100x increase in throughput).
If you have enough memory and are running 64-bit R, load the raster into
memory rat
I found some problems with this code, but no worries, I got subplot() working
and finished my figure. The resolution and sharpness of the figure is
definitely superior to a version I had been trying to make with ArcMap.
1.R-3.0.1 didn't find subplot() in TeachingDemos, but it did in Hm
Here is an example using the subplot function:
library(TeachingDemos)
library(maptools)
data(wrld_simpl)
plot(wrld_simpl,border=NA,col='blue',axes=TRUE, xlim=c(100,130),
ylim=c(-40,30))
tmp <- subplot(
plot(wrld_simpl,border=NA,col='green',bg='white',axes=FALSE),
'bottomleft', size=c(2,1),
Hi,
I would like to calculate kernel areas using kernel.area function from
adehabitat package.
My original coordintes are in latidude/longitude. I suppose that as an
input for the function the coordinates should be in other fomat - I assume
it is UTM coordinates (seeing the example in help, I ha
Thank you Oscar, I will try it!
By the way, I attached the dput(limit_NBver10) on the end of my previous
email (not as an attachment to he email).
Zuzana
On 10 September 2013 15:25, Oscar Perpiñán Lamigueiro <
oscar.perpi...@gmail.com> wrote:
> Hello,
>
> The next code combines three differen
Dear Roger Bivand,
thank you for your answer. Indeed yesterday I tried to post the question
however I received a warning telling me that I was not allowed to post
(I did not register yet to the mailing list), so I assumed the post did
not go through. I apologize for the double post.
I am act
You asked the same question yesterday (but simply added it to the original
thread from 2009). Do read the instructions for posting and the posting
guide. This is not a question concerned with the use of software, indeed,
had you used the software, you could have examined the questions you ask
e
Dear list,
I found a message asking same kind of things I am wondering.
Unfortunately I dont find proper answers and thus would like to update
the topic. Maybe Xingli could you share what your learn from the authors
with us to the questions below?
Regarding the weights, is it imperative for
Dear list,
What parameter set allows one to perform automatically Indicator Kriging in
automap?
Dr. Francesco Carotenuto (Ph.D.)
Department of Earth Sciences
Federico II University
Largo S. Marcellino 10,
80138 NAPLES (ITALY)
http://francesco-carotenuto.blogspot.it/
http://3dpaleontology.blogsp
On Fri, Sep 13, 2013 at 8:06 AM, Lionel Hertzog wrote:
> Dear Nevil Amos,
>
> Maybe the function 'reclassify' in the raster package is what you are
> looking for.
This should be trivially parallelizable, so if you have a cluster
handy there should be a way to split the raster and farm out the tas
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
In the script, you set the wrong CRS to the long/lat points. Set the
right CRS, then spTransform to the new one. Check that they match
after reprojection by plotting with axes, eg.
plot(muni.sp,axes=T)
points(muni.sp, col = 'red')
Second, when you wo
Dear Nevil Amos,
Maybe the function 'reclassify' in the raster package is what you are
looking for.
Yours,
Lionel
On 13/09/2013 08:36, nevil amos wrote:
I have a large raster (4gb as a .asc file) with 197000 values, that I need
to reclassify by values in an large lookup table, 197,000rowsx8
12 matches
Mail list logo