Elizabeth,
Creating a raster object from a nc file is probably easiest done like this:
r - raster(filename, band=1)
Or to get the all:
b - brick(filename)
then:
ba - aggregate(b, )
Normally, you should be able to do this to save it as ncdf:
ba - writeRaster(ba, filename='out.nc')
Or
Recently, Robert kindly added a function to raster package (raster_1.5-14
from RForge) called lineValues. After loading the new package, I am unable
to access (or use) the function. Any ideas what I might be doing wrong?
rst.cln - raster::lineValues(lns = rst.lines, rst.cellnum)
Error:
Hello again,
okay, i found the ncf package and the correlog function.
# correlog(x, y, z, w = NULL, increment, resamp = 100, latlon =
FALSE, na.rm = FALSE, quiet = FALSE)
my short question is:
if latlon = T and the positions are "degree,
Jens,
in your code snippet below, the number 6370 rings a bell. According to
http://en.wikipedia.org/wiki/Earth_radius it is close to the average
earth radius. In km. Given that there's just one such number, it seems
the code assumes a sphere is close enough to approximate the earth.
On
I want to read Iknos image and then want to see display in R. Just read and
display nothing else.
I try using raster package
rasterfile-system.file(C:/iknos/edwards_airbase_usa_1m_tc.tif,package=raster)
Then plot with the help of this command
plot(rasterfile)
Error in plot.window(...)
Hi Malik,
To a look at the rgdal package to read your tiff file and spplot from
the sp package to plot it.
cheers,
Paul
On 10/06/2010 01:36 PM, Malik Shahzad wrote:
I want to read Iknos image and then want to see display in R. Just read and
display nothing else.
I try using raster package
I believe there is a mistake in your command. You skipped the raster
command itself. From the raster help files here is the correct way:
r - raster(system.file(external/test.grd, package=raster))
On Wed, Oct 6, 2010 at 7:36 AM, Malik Shahzad geoma...@live.com wrote:
I want to read Iknos
Indeed, and as Malik is not using a system file (files packaged with
R to serve as examples), you can should leave that bit out and do:
r -raster(C:/iknos/edwards_airbase_usa_1m_tc.tif')
plot(r)
Robert
On Wed, Oct 6, 2010 at 4:44 PM, kapo coulibaly kmcou...@gmail.com wrote:
I believe there is
Dear thanks for help,
In first line there is no problem
rasterfile-system.file(C:/iknos/edwards_airbase_usa_1m_tc.tif,package=raster)
This works fine without any error, It mean image has been loaded into variable
raster file, but now how i can display raster file in R. I mean I want to
Hello
I am trying to use the correlog function to estimate a spatial correlogram
for the residuals of a logistic regression and I have run accross the
following error.
summary(binom1 - glm(Use~X20mslop+X20mdem+soilsst, family=binomial,
+ data=M60m2000NE_1.df))
correlog1.1 -
Hello Jenn
Unfortunately I can not offer any advice but I am facing very similar
challenges and actively searching for answers. I have presence/absence data
for a colonial species and I am trying to compare the performance of a
logistic model with spatial models: 1) a glm including an
Greetings,
I have two questions/problems related to spplot:
1. I am trying to plot one SpatialPolygonsDataFrame over another, and
fill the background polygon with grey. I can't seem to fill the
background polygon without having it cover over my supposed foreground
polygon, even if I specify
12 matches
Mail list logo