Re: [R-sig-Geo] Intersect & Plot a Polygon and wkbLineString

2009-02-16 Thread Jim Burke
Hi Everyone, Sorry folks, my email program, Netscape Communicator, cannot send text email correctly. So this attempt goes out in both html and text. I have a nice solid SpatialPologonsDataFrame and a new line file (wkbLineString) called "major roads.shp". MY QUESTIONS. I would appreciate exam

[R-sig-Geo] Intersect & Plot a Polygon and wkbLineString

2009-02-16 Thread Jim Burke
I have a nice solid SpatialPologonsDataFrame and a new line file (wkbLineString) called "major roads.shp". MY QUESTIONS. *I would appreciate examples of* *a)How to merge the intersection of a **SpatialPologonsDataFrame and a wkbLineString ** b)How to plot the result c)How to display the s

Re: [R-sig-Geo] Point pattern analysis

2009-02-16 Thread Michel Barbosa
Dear Tom, Actually, R complains about duplicate points. I'm pretty sure some restaurants are really close to each other, or perhaps in the same building. Could be some geocoding errors as well. If I were to 'wiggle' the data, would this be dependent on the resolution of the dataset? I'm still figu

Re: [R-sig-Geo] Point pattern analysis

2009-02-16 Thread Michel Barbosa
Dear Nicholas, But from your description, it seems like there might be other questions > that should guide your analysis. > Context should drive your exploration of the data. > The most important context is the user's current location in my case since I'm researching the usefulness of SDA for Loc

Re: [R-sig-Geo] Point pattern analysis

2009-02-16 Thread Michel Barbosa
Dear Virgilio, Interesting problem. Let me know if you need help collecting data. ;) Thanks :) Actually, I'm busy with developing a Location-Based Service (a restaurant finder to be precise) utilising SDA. The goal of my research is to integrate SDA in an LBS. For this purpose, I've gathered abo

Re: [R-sig-Geo] Point pattern analysis

2009-02-16 Thread Tom Petersen
Michael, Although your points might seem cluttered (are there reasons to believe that they shouldn't be?), you should still check what your spatial resolution is, and if there are several restaurants occupying the same coordinates - in that case, only the first one will be considered on import. I

Re: [R-sig-Geo] SpatialFiltering in spdep: How to create an 'nb' object from an imported spatial weights matrix

2009-02-16 Thread Timo Mitze
Dear Roger, it works now, many Thanks! Timo -Ursprüngliche Nachricht- Von: Roger Bivand [mailto:roger.biv...@nhh.no] Gesendet: Sonntag, 15. Februar 2009 23:00 An: Timo Mitze Cc: r-sig-geo@stat.math.ethz.ch Betreff: Re: AW: [R-sig-Geo] SpatialFiltering in spdep: How to create an 'nb' o

[R-sig-Geo] GEOSTAT 2009 Summer School, 3-10 May 2009 MEDILS, Split, Croatia

2009-02-16 Thread Hengl, T.
1st call: GEOSTAT 2009 Summer School "Spatio-temporal data analysis with R + SAGA + Google Earth" Period: 3 – 10 May 2009 Location: Mediterranean Institute For Life Sciences (MEDILS), Split, Croatia The objective of the GEOSTAT Summer School is to promote use of open source tools for the analy

Re: [R-sig-Geo] Point pattern analysis

2009-02-16 Thread Nicholas Lewin-Koh
Hi, Yes, if you need help rating restaurants, put me in your grant too :-) Seriously, there are many ways to skin a cat. I don't think cartograms will help you much in this particular case. If you have data besides your point pattern, eg postal codes, census data, zoning, ... You could look for the

Re: [R-sig-Geo] Masking interpolations

2009-02-16 Thread Paul Hiemstra
Hi, Could you provide us with some reproducible code? In addition, look at the output of traceback() to get more info on what exactly went wrong. cheers Paul, Wesley Roberts wrote: Hi Paul, I have managed to get the code to run and have figured out what the function does. Many thanks, it s

Re: [R-sig-Geo] Masking interpolations

2009-02-16 Thread Wesley Roberts
Hi Paul, I have managed to get the code to run and have figured out what the function does. Many thanks, it should do exactly what I need. Unfortunately I have now run into another problem. The code appears to run properly until the error below appears. . . . . 60963 (-76170.4, -3300950) 20.06

Re: [R-sig-Geo] Masking interpolations

2009-02-16 Thread Paul Hiemstra
Hi, If you want to use the standard convex hull of automap, just omit the new_data argument. That will stuff 5000 points into the convexhull on a regular grid. If you want to specify the pixelsize you need to adapt the create_new_data function like this: create_new_data = function(obj) { # F

[R-sig-Geo] Personal invitation from srinivasa raghavan

2009-02-16 Thread srinivasa raghavan
Hi, On 2/2/2009 3:10:07 AM, you were invited to join srinivasa raghavan's UNYK address book so he/she would always have access to your contact info and you to his/hers. To accept his/her request, Click here. http://www.

Re: [R-sig-Geo] Masking interpolations

2009-02-16 Thread Wesley Roberts
Excellent, that will help alot. I have the newest version of automap and am running an experiment now. Just to confirm I need to omit the statement 'grd(new_data argument)' which in my case is 'gridded(grd) <- TRUE'. That way the resolution is defined in x.range <- . and y.range <- . s

Re: [R-sig-Geo] Masking interpolations

2009-02-16 Thread Paul Hiemstra
Hi, A convex hull should do the trick, if you have the latest automap version (>0.5-2) you don't have do to anything to get a convex hull. Just omit grd (new_data argument) and automap will automatically use a convex hull. cheers, Paul Wesley Roberts wrote: Hi Paul, Many thanks for the re

Re: [R-sig-Geo] Masking interpolations

2009-02-16 Thread Wesley Roberts
Hi Paul, Many thanks for the reply. I am not sure about the convex hull approach as I am not sure how to implement it as part of my program. Is the code you wrote below replacing the following statements? Would I then pass grd to the autoKrige statement as shown in the original code below? ***

Re: [R-sig-Geo] Masking interpolations

2009-02-16 Thread Paul Hiemstra
...in addition: You could also use AcrGIS, QGis or another GIS package to make a polygon that delineates the areas in which you want to interpolate. You can then use spsample() like I showed you in the previous mail. Making this polygon can also be done in R using spplot and grid.locator, if y

Re: [R-sig-Geo] Masking interpolations

2009-02-16 Thread Paul Hiemstra
Hi Wesley, You could take a look at using a convex hull. I'm not sure if this will fix your problem as we cannot see how exactly your points are irregular. The latest version on my website (0.5-2) uses a convex hull off the data if you don't pass a new_data object. You could try this. The func