search for: map2poly

Using the maptools function "map2poly" I have created a map/polylist object - a map of the Great Lakes. My goal is to passively plot additional data on this map. Unfortunately I am not able to change (distort) the relative scale of the axes for the map, and so the sample points do not line up correctly on the map. In the code bel...

Hi Guys Thanks for this! As I am a beginner, I bet I'm running into some really basic problems. Using the example from the Map2poly function in maptools: try <- read.shape("euadmll.shp") mappolys <- Map2poly(try) # this produces 14 warnings() of the type: # 1: level 2 circularity at 62 in: .mtlbuild(.mtafters(r1), rD) # etc # 14: From next release, default hole handling will change in: Map2poly(try) plot(mapp...

...web, and I am at the end of my own resources. I need to do analysis of a spatial point process occuring with four counties, using spatstat. I've had no trouble importing the shapefiles, but I can't seem to figure out what I need to do in order to use the boundary as an owin. Why I use the Map2poly command, it appears to work, but not when I attempt define the window using the shapefile; it also gives me errors when I try to define a ppp.object. I am using the latest release of Raqua on G4. My process looks like this: #begin code > library(maptools) > Shapefile<-read.shape("co...

Dear R experts, I am new to the list and R software. I need to convert arcview file to Winbugs having R has middle package. Got from friends that it is possible following the steps: 1. Converting arcview shapefile to "cgm clear text file" 2. Downloading "convert.r" to into R 2.1 use source("convert.r") 2.2 convert("filename_cgm") without the

...#________________________________________________________ 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 <- subset(mappolys, nc$att.data$NAME!="Honolulu, HI") nc$att.data=subset(nc$att.data, nc$att.data$NAME!="Honolulu, HI") nc$att.data[,1]=as.numeric(paste(nc$att.data$MSACMSA)) #attributes(nc$att.data) nc$att.data=merge(nc$att.data, dat, by.x="AREA", b...

...SparseM (0.71) loaded. To cite, see citation("SparseM") Loading required package: boot > library(maptools) > district.shp <- read.shape("~/Documents/SAHSU/MD/data/Carthography/ districts_ok.shp") Shapefile type: Polygon, (5), # of Shapes: 354 > district.pl <- Map2poly(district.shp, as.character(district.shp$att $DISTRICT_2)) > plot(district.pl) *** caught segfault *** address 0xc0000165, cause 'memory not mapped' Traceback: 1: polygon(coords[pFrom[i]:pTo[i], ], border = border, xpd = xpd, density = density, angle = angle) 2: polygonholes(x[[...

Hi! I have no idea if this is maybe an easy task utilizing R since I read there is geographical map data in some package: I have a huge number of geographical points with their coordinates in Germany. Now I want to determine for each point in which "Bundesland" = state it is located. Can anybody tell me if this is doable relatively easy in R and if so give me some hints or links

I am a very new user of R (ver 2.2.1 & Windows XP) I would like to see if I can plot some simple data (Number of people in certain regions of Canada and within Ontario and Quebec. Can anyone suggest a) a souce of Canadian maps and b) some simple worked examples of using maps. At the moment I have managed to load the map and mapdata libraries but that is about the extend of my knowledge.

...s as the shapefile and one column populated with numbers from 1 to 8 which represented an attribute for each of those states. The code then plotted those numbers as colours on the map so each state had the appropriate colour. This was the code: x<-read.shape("NGADMIN2.SHP") xpolys<-Map2poly(x) coverage<-read.table("coveragemar08.txt",header=T) colour1<-coverage[,2] palette(c("red4","red3","red2","orangered","orange","gold1","yellow","white")) plot(xpolys, axes=F, colour1) I've now come...

Hi, Is there any way to plot a Map with different shading patterns, instead of color intensity, in R2.0.1 over Linux ? Thanks. D ----------------------------------------- Stay ahead of the information curve. Receive GIS news and jobs on your desktop daily. Subscribe today to the GIS CafeNews newsletter. [ http://www10.giscafe.com/nl/newsletter_subscribe.php ] It's informative and

Hi: Is there any function in R which can convert the long/lat to UTM(Universal Transverse Mercator)? There are quite a few converters on Internet. However, the interface is designed as input->output which I can not convert lots of locations at the same time. Another question is whether there is a function in R which can tell the time zone from the location's lat/long? Thank you!

Dear All, In the moment i'm using the map and maptools package to read shapefiles and display the maps. I'm looking for the possibility to draw points (randomly positioned or positioned according to a grid) into the polygons instead of filling the polygons with colors. For example: a map (shapefile) with 10 countries, 15 points in the polygon of country A, 20 points in the

Hi, In another thread ("PBSmapping and shapefiles") I asked for an easy way to read "shapefiles" and transform them in data that PBSmapping could use. One person is exploring some ways of doing this, but it is possible I'll have to do this "manually". With package "maptools" I am able to extract the information I need from a shapefile but it is