similar to: maptools: Test if point is in polygon

Greetings! I am trying plot some data on a map in R. Here's the scenario. I have a variable called probworkinghealthy which contains a predicted probability of employment for every individual in my sample (about 100,000 observations). I have another variable, called a001ter, which contains the subject of residency in the Russian Federation (akin to a US state) for every individual in the

A few years back, I wrote some code to plot maps with the maptools package. Now I am trying to reproduce my results, only to find out that maptools has been updated and my code no longer works. I've been able to fix the first part of it by forcing spb <- maptools:::read.shape("/home/sasha/Documents/maps/spb.shp") (as read.shape has been deprecated) but now when I do plot(spb,

Hi, I meet a strange trouble writing shapefiles with write.pointShape() from maptools. > write.pointShape(mylines[,5:9], file="Sxtplinshp", mycoordspt) Error in write.pointShape(mylines[, 5:9], file = "Sxtplinshp", mycoordspt) : shapefile names must conform to the 8.3 format It seems that the function request a special shapefile name spelling (conformed to 8.3

cheers , I need help in shapefile manipulations. I have two shapefiles of my country Cameroon. The first contain 10 provinces and each province contains a certain number of administratives units. I dont have ESRI Arc view. I want to add a admistrative unit of the center province to the map with 10 province. Is it possible to to take a peace of shapefile and combining it to another ? Thanks !

Greetings! I attempted to read a STATA 10 file using the following syntax: library(foreign) data <- read.dta("~/Documents/agestandard/analysis.dta") I got the following error: *** stack smashing detected ***: /usr/lib/R/bin/exec/R terminated ======= Backtrace: ========= /lib/tls/i686/cmov/libc.so.6(__fortify_fail+0x48)[0xb7c13138]

The main limitation of the shapefiles package that I put together is that it does not create shapefiles from R objects - rather it only writes shapefiles that have been read into R and manipulated within the constraints of the existing file structure. By this I mean that for example you can change the coordinates of points and write them back out. Or you can add a bunch of blank columns in the

Hi, I want to use map("state") and have the ecoregion shape (please see link) file projected onto this. ftp://ftp.epa.gov/wed/ecoregions/cec_na/NA_CEC_Eco_Level3.zip ftp://ftp.epa.gov/wed/ecoregions/cec_na/NA_CEC_Eco_Level3.zip Could someone please show me how; I have never messed with this sort of thing. Thanks. -- View this message in context:

I am having difficulty loading a points shapefile (1600 records) I generated from ArcGIS as a PointsSpatialDataFrame using the maptools readShapePoints function. library(maptools) surveypts <- readShapePoints(system.file("C:/temp/survey.shp", package="maptools")[1] Error in getinfo.shape(filen) : Error opening SHP file I've tried numerous re-exports and renames of the

Greetings: I would like to kindly ask for a little help. The rough code is: #________________________________________________________ dat=data.frame(read.delim(file="all.txt", header = TRUE, sep = "\t", quote="\"", dec=".",na.strings = "NA")) nc=read.shape("astae.shp", dbf.data=TRUE, verbose=TRUE) mappolys=Map2poly(nc) submap

Dear all, I would like to produce a map with information about the patenting activity in German districts, by coloring districts with different degrees of patenting activity in different colors. I work with the packages maptools, maps and spdep. The map data is read from an external .shp file (+ the corresponding .shx and .dbf files). Plotting a map with the IDs or the patenting indicator itself

Hello. I'm new on that forum, so I would use some help... I'm programing with R program, have some exercise to do. So here it goes. 1. I have to upload two libraryes "sp" and "maptools" 2. This command is troubling me : towns <- readShapeSpatial("OB/OB.shp") 3. When I'm pressing f5, there is an error called: Error: could not find function

Hi, I have a problem that concerns combination of the package Maptools and Splancs I have 2 shapefiles that i want to manipulate (one of type point and one polygon).I import them in R using Maptools but then i can't estimate a quartic Kernel using Splancs. The package doesn't recognize the shapes (invalid points and poly argument).I don't know if this is an easy task but i have

Hi, I'd like to read a shapefile into a Map object. This is exactly what read.shape{maptools} is suppoed to do, according to the documentation I found in the link below. However, this doesn't seem to work library(maptools)x <- read.shape(system.file("shapes/sids.shp", package="maptools")[1])Error: could not find function "read.shape" I haven't

I'm not sure I'm barking up the right tree here, but would I need to make use of groups to plot two separate datasets within ONE panel in xyplot? The desired end result is a single xy plot of two separate (but similar in values and ranges). Full code follows, xyplot code at bottom #########Determine Frequencies ##########coastal_slope #needs the maptools package to read ESRI grid

Dear all, I would like to employ a European map in a presentation. My idea was to use: library(mapdata) map("worldHires", c("Austria", "Switzerland", "Germany")) where I include all countries from my analysis as a vector of character strings like in the example above. Unfortunately, I was unable to specify the Baltic States (Latvia, Lithuania, Estonia) or

Sorry, I did not make my question clear. Since I have a point theme with many points, some of them may clump together. the problems here are: 1. how to find clumps in a point theme? 2. the convex-hull extension I found only deal with all the points in a theme at each time? how to make each convex hull around each point clump automatically? Thanks. Xiaohua On 5/12/06, Bob Booth <bbooth

Hello R Community, I have imported a dataset which contain X Y coordinates and would like to recreate a shape file after some data analysis. What i have done is to import some taxlot data and join them based on some criteria. I want to check to see how well the joining went by reviewing the results in GIS. A couple things. I cant seem to import a shape file correctly using the maptools

Greetings. I'm trying to label some points using the identify() command. The labels and the coordinates are coming from a csv file. The file encoding is UTF-8. If I use labels in Latin characters, everything comes out OK. But when I switch to using labels in Cyrillic, I get garbage. To check that my file is not corrupt, I've attempted to create the labels using the text() command and

Hi, I need to compile Icecast on a kind of ancient AIX machine. I compiled all the required libraries, except libvorbis, and I can't make it compile at all. Is there a way to compile Icecast without Ogg Vorbis support? Thanks! -- Regards, Alek Andreev alek@zvuk.net

Hi, I had my stream running from ezstream 0.1.2 to icecast 2.2.0. It's been going for about a week, but two days ago I changed my icecast.xml. Since then, Icecast dropped my sources and stopped transmission twice. Everytime, it spew out xmlEncodeEntitiesReentrant : char out of range on STDERR and [2005-02-07 01:14:47] WARN source/get_next_buffer Disconnecting source due to socket