Michael Lawton wrote:

> I am trying to create a pixel image from the marks of a marked
> spatial point pattern for use with the Kinhom function in the
> spatstat library and am having some difficulty.  Is there a way to
           ^^^^^^^
        ***package*** !!!
        (NOT ``library'')

> create such an object using x-y coordinates and values at those
> coordinates?

        It is not clear to me what you actually want to do, but there
        are tools in spatstat which may accomplish your purpose.  In
        particular, the function as.im() (see the on-line help) will
        take a function of the form ``function(x,y,...)'' (and a
        window W) and convert that function into an image defined on W.

        So if you have or can create a function which will calculate
        the ``values at those coordinates'' then there is no problem.

        Otherwise you will need to explain more clearly what it is
        that you want.

        If your are thinking of an ``image'' where the value is, say,
        10, at points of the pattern whose marks are of type 1, and
        say -10, at points of the pattern whose marks are of type 2,
        and say 0 at ``background'' points, then this is readily
        accomplished using nearest.raster.point().  But I don't think
        you gain anything by representing the point pattern as an image
        in this way, and this is not the sort of thing that one would
        feed to Kinhom().

        Quite often an input to Kinhom() would be calculated by
        ``smoothing'' the field of ``delta functions'' corresponding
        to a collection of points.  The function density.ppp() is
        a good way to accomplish this.  Choosing the ``bandwidth''
        parameter sigma for the smooth is however a bit of a mystery,
        apparently even to the cognoscenti.

                                cheers,

                                        Rolf Turner
                                        [EMAIL PROTECTED]

P. S.  A library, as someone --- Martin Maechler??? --- has repeatedly
tried to explain to the list, is a *collection of packages*.  One
loads a package *from a library* by using the library() function.

                                        R. T.

______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to