Thanks for the hint Michael
there must have been a change in spelling the function name. The maptool
function is spelt "ContourLines2SLDF" rather than "contourLines2SLDF" as
it was spelt before. Anyway I am happy it is still there and it works.
Cheers
Karl
Michael Sumner <[EMAIL PROT
Hello list,
I used to be able to apply the function contourLines2SLDF for the
following operation on a SpatialGridDataFrame object.
> CL <- contourLines2SLDF(contourLines(as.image.SpatialGridDataFrame
+ (as(ht.idw["var1.pred"],"SpatialGridDataFrame"
The last time I ran this function
Hi Edzer,
your suggestion worked well. I saved the plot as an object named z1. In
the legend variable I found the settings for symbol size. In my case there
were 4 entries which I changed from 0.1 to 1
using the following statement.
z1$legend$inside$args$key$points$cex <- c(1,1,1,1)
plot(z1)
R
Hello list,
I was wondering if it were possible in spplot() to influence the size of
the plotting symbol appearing in the legend independently of the size of
the plotting symbol in the graph.
I have a situation where the symbol in the graph needs to be very small,
i.e. cex=0.1. Consequently, the
Hello list
I am stuck with a problem regarding the manipulation of two data sets.
The first data set contains electromagnetic induction (EMI) readings of a
field obtained with a high sampling density,14965 records in a field 1000
by 500 m. Each reading is geo-referenced.
The second set contains s
When performing ordinary kriging under gstat in R I received the following
error after the program had stopped.
---R output
vert.ok <- krige(vert~1, em38a, SPDF, model=vert.fit)
[using ordinary kriging]
"chfactor.c", line 130: singular matrix in function LDLfactor()
gstat caught an error that oc
As suggested by Roger Bivand and Edzer Pebesma I have included the
output of str(t1) and str(SPDF) below.
The error persited when I tried to do the kriging in two steps, using gstat
() and project().
Kriging a different variable, eg. "band1", results in the same error.
Converting t1 to a SpatialP
I am encountering a error message when attempting to do "inverse distance
kriging" or variogram calculations. The error is always the same
ERROR: dimensions do not match: locations 1513600 and data 611452.
What puzzles me is the fact that I am using the same code as I have
successfully
used with
Hello
I was wondering if there is a way in "R" or "sp" to calculate shortest
distance from a given object such as a river, similar to the data example
"meuse.grid" under the gstat package. I was wondering how the "dist"
variable in the "meuse.grid" file was orignially calculated? I would like
|-+>
| | [EMAIL PROTECTED]|
| | nl |
| ||
| | 29/06/2006 17:28 |
| ||
|-+>
>
|-+>
| | [EMAIL PROTECTED]|
| | nl |
| ||
| | 27/06/2006 19:12 |
| ||
|-+>
>
ur
> lines from such information and overlay it onto the existing interpolated
> graph?
>
> Cheers
>
> Karl Sommer
>
Simpler than I thought it would be: add contour = TRUE to the argument list
of spplot.
If you want the actual contours, there is
density as the sampling locations
indicated with a grey x in fig08. Is there a method to generate contour
lines from such information and overlay it onto the existing interpolated
graph?
Cheers
Karl Sommer
_
Karl J Sommer,
Department of Primary Industries,
Catchment
Roger
of course, the problem was, as you suggested, that I did not adjust the
cell dimensions.
Thanks again
Karl
|-+>
| | [EMAIL PROTECTED]|
| | no |
| ||
| |
Yes, the grid coordinates reflect the c(0,5, 0,5)
[1,] 617400.0 6190061.0
[2,] 617400.5 6190061.0
[3,] 617401.0 6190061.0
[4,] 617401.5 6190061.0
[5,] 617400.0 6190060.5
[6,] 617400.5 6190060.5
But for some reason, when using the smaller cellsize my spplot of the
interpolated field all of
after some fiddling with the GridTopology() method I eventually succeeded
in making it work and now can interpolate onto the grid.
Tremendous, thanks.
One problem I encountered when specifying a number smaller than 1 for the
cellsize argument eg. c(0.5, 0.5) instead of c(1,1), I get erroneous
r
thanks Roger Bivand, this is really powerful and much less cumbersome than
I feared.
I ran an inverse distance interpolation on some of the screened data. In
order to do so I created a destination file from scratch with a 2 x 2 m
grid. The destination file is of rectangular shape and is larger
Hello list
I have been trying to figure out if it is possible to select a subset of
data in a SpatialGridDataFrame without prior conversion to a data.frame
using the method as.data.frame().
I have imported an ESRI " *. bil " file containing different spectral bands
using rgdal. The file imported
Hello list
thanks to the helpful advice given by list members I have been able to read
and re-project vector based shape files and attributes data using readOGR()
under rgdal for further analysis (kriging, k-means clustering) and
plotting.
Using the sp.layout under spplot allowed me to precisely
thanks for the prompt reply
the method spTransform() worked
I now have an attribute file of class SpatiaPointsDataFrame with EM38
readings. I also have a separate shape file of class
SpatialPolygonsDataFrame which I imported form an ESRI shape file using the
readOGR() method of rgdal. I would
Hello list
I have been trying to do a projection conversion from lat lon to a UTM grid
I followed an example in the April discussion list provided by Roger Bivand
and adapted it to my own data.
Everythings seemd to work well up until the transform() method which failed
to execute. The failure se
the opposite
ie. demote the same data set from SpatialPointsDataFrame back to a
data.frame?
Regards
Karl Sommer
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
22 matches
Mail list logo