search for: readshapespati

...ile 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 package and the readShapeSpatial. I have tried Building=file("data/input/BuildingShape/Building.shp") Bldg<-readShapeSpatial(fn=data/input/BuildingShape/Building,proj4string=NAD83) #---------------------- Bldg<-readShapeSpatial(data/input/BuildingShape/Building,proj4string=NAD83) #--------------- Building=file...

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 "readShapeSpatial" So, I would use some help. Thanx in advance. Gaja -- View this message in context: http://r.789695.n4.nabble.com/Problem-with-uploading-library-tp3170455p31...

....jpg (but without the top right area, just the different grayscale kriging area). library(xlsReadWrite) library(maps) library(mapdata) library(maptools) Pontos<-read.xls('Pontos.xls',sheet=1,rowNames=T) Prof<-read.xls('Pontos.xls',sheet=2,rowNames=T) Batimetria<-readShapeSpatial(‘Batimetria.shp') Municipios<-readShapeSpatial('Municipios.shp') Recifes<-readShapeSpatial('Recifes.shp') par(mar=c(4,0,0,0)) map('worldHires','brazil',ylim=c(,),xlim=c(,),type='n') plot(Batimetria,ylab='',xlab='',border='...

Hi ALL, Can you please help me on how to determine the prediction error for ordinary kriging?Below are all the commands i used to generate the OK plot: rsa2 <- readShapeSpatial("residentialsa", CRS("+proj=tmerc +lat_0=39.66666666666666 +lon_0=-8.131906111111112 +k=1 +x_0=0 +y_0=0 +ellps=intl +units=m +no_defs")) x2 <- readShapeSpatial("ptna2", CRS("+proj=tmerc +lat_0=39.66666666666666 +lon_0=-8.131906111111112 +k=1 +x_0=0 +y_0=0...

...s Dados<-odbcConnectExcel(''Campos.xls'',readOnly=T) Coord<-sqlFetch(Dados,''CoordMed'',colnames=F,rownames=''Ponto'') odbcClose(Dados) N<-pixmap(''Norte.png'',nrow=166,ncol=113) # Carregar pontos e shapes Batimetria<-readShapeSpatial(''C:/Users/Rodrigo/Documents/UFPR/Micropaleontol ogia/Campos/ShapeFiles/Batimetria_BC.shp'') Estados<-readShapeSpatial(''C:/Users/Rodrigo/Documents/UFPR/Micropaleontologi a/Campos/ShapeFiles/Estados_Sudeste.shp'') Faciologia<-readShapeSpatial(''C:/U...

Hi R users I need the map of France?s ? communes ? (towns) to build a map Is there a way to get it? More generally: How to do to get the map of a country and its different geographical levels? Best regards -- View this message in context: http://old.nabble.com/map-of-a-country-and-its-different-geographical-levels-tp26225645p26225645.html Sent from the R help mailing list archive at

...uild choropleth. Basically i am trying to colour regions on the map by population. I possess the shape file of the country, and also the population data, however, i am having trouble to create the plot, below is my code: population=read.csv("nz2.csv") population names(population) nz=readShapeSpatial((sprintf('nz-geodetic-marks',tmp_dir))) nz$land_distr mapnames = nz$land_distr mapnamesnew=as.character(mapnames) popnames =as.character(mapnamesnew) pop =population$People pop2=as.numeric(pop) popdich = ifelse(pop2 < 100000, "red", "blue") popnames.lower = to...

Dear all, I am currently working with the spatstat package, using windows and pixel images. First: My aim is to transform a shapefile (see attached) into a pixel image. My idea is to start transforming the shapefile into a Spatial Polygon file: x <- readShapeSpatial("200001441.shp") y <- as(x, "SpatialPolygons") z <- as.owin(y) Given z, I want to identify each polygon with a single constant value. This is like adding marks to the SpatialPolygons file. Then I want to convert these polygons into an image, such that the value of each...

Dear all, I am trying to create a map for central asian countries(Kazakhstan, Uzbekstan, Kyrgyzstan, Turkmenstan, and Tajikstan). I tried google for shapefiles and Rdata files of central asia, but I can't seem to find them. There is only a world map and individual maps for each countries. If I use individual maps to create this map, how do I put them on a same graphic with different colours?

...id=115&parent=category Problem: I want to extract centroid coordinates for each State, but there is some coding problem with the shapefiles that prevents this. #Code: #After extracting the shapefiles from the india_state.zip file, then: library(maptools) library(sp) library(rgeos) india <- readShapeSpatial("india_state.shp") #Some states are made up of more than one polygon. #State names are in the NAME column of the dataframe india at data plot(india[india$NAME == "Tamil Nadu",], col = 1:8) #Extract centroid for a State (this works): gCentroid(india[india$NAME == "Tamil...

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 currently working with the spatstat package, using windows and pixel images. > > First: > > My aim is to transform a shapefile (see attached) into a pixel image. > > My idea is to start transforming the shapefile into a Spatial Polygon file: > > x <- readShapeSpatial("200001441.shp") > y <- as(x, "SpatialPolygons") > z <- as.owin(y) > > Given z, I want to identify each polygon with a single constant value. This is like adding marks to the SpatialPolygons file. Then I want to convert these polygons into an image, such tha...

Dear list, I have a question regarding the included polygon in the khat function of the splancs library defining the area where points appear. I have not only one simple polygon included, my map includes several islands. I read my map which was a shape file by readShapeSpatial of the maptools library. Then I split up the SpatialPolygonsDataFrame into the included polygons by coor <- NULL for (i in 1:[number of polygons]) coor<- rbind(coor,map@polygons[[1]]@Polygons[[i]]@coords) I am using now this newly generated polygon "coor" in the khat command....

...2 2 US002 8 3 US003 3 4 US004 5 5 US005 6 6 US006 7 7 US007 9 8 US008 9 9 US009 10 10 US010 11 11 US011 13 12 US012 15"),head=TRUE,stringsAsFactors=FALSE)) library(maptools) library(ggplot2) library(gpclib) gpclibPermit() setwd("C:/Documents") us_counties.shp <- readShapeSpatial("uscounties.shp") us_counties.shp.p <- fortify.SpatialPolygonsDataFrame(us_counties.shp, region="CNTY_FIPS") us <- merge(us_counties.shp.p, us_counties.shp, by.x="id", by.y="CNTY_FIPS") p <- ggplot(data=us, aes(x=long, y=lat, group=group)) + ge...

...elow data will work but you will have to set your directory for the shape file. Never did this before so hopefully this works. ResProbs is just supposed to be a value between 0-1, sorry if that more complicated than it needed to be. TazFile <- "*directory*/TAZ.shp" TazShape <- readShapeSpatial(TazFile) TazShape<-TazShape[order(TazShape$TAZ_NUM),] ResTaz<-25:666 ResProbs<-rnorm(642,0:1) ResProbs[ResProbs>1]=.5 ResProbs[ResProbs<-1]=.2 ResProbs<-abs(ResProbs) ResProbs..<-data.frame(ResTaz,ResProbs) names(ResProbs..)<-c("Taz","SFsubM") TazShap...

Dear all, in GADM map there are three levels (nation, province and precinct) for each country of the world but for all of them you are never able to plot only one part of a chosen country. To be sure, I am trying to plot only one region of “Italy” and colour the different precincts in it. So far I am able to colour only the part of my interest but the programme still plot the whole country. Is

Dear list, I'd like to use the function smooth.map{maps} to smooth out aggregated spatial data. The first argument to this function must be a map object. I have a shape file which I read with readShapeSpatial{maptools} and I got an object of class "SpatialPolygonsDataFrame". Is it possible to create a map object from a SpatialPolygonsDataFrame object? Alternatively, is there any other function similar to smooth.map that will work with a SpatialPolygonsDataFrame? Thanks for any help, Axel....

Dear all, I am currently working with the spatstat package, using windows and pixel images. First: My aim is to transform a shapefile (see attached) into a pixel image. My idea is to start transforming the shapefile into a Spatial Polygon file: x <- readShapeSpatial("200001441.shp") y <- as(x, "SpatialPolygons") z <- as.owin(y) Given z, I want to identify each polygon with a single constant value. This is like adding marks to the SpatialPolygons file. Then I want to convert these polygons into an image, such that the value of each...

...----------------------------------- #configuramos donde setwd("D:/mapas/") #Carga de librerias library (sp) library (RColorBrewer) library (maptools) # carga de mapa politico y datos datos_manies<-read.csv("D:/mapas/corrientes/tabla.csv", sep=",") cartografia<-readShapeSpatial("D:/mapas/corrientes/corrientes.shp",proj4string=CRS("+proj=longlat +datum=NAD27")) # comprobamos si el mapa cargo, si aparece el grafico esta bien # si no es asi hay que revisar las librerias y dependencias plot(cartografia) # para ver la paleta de brewer # display.brewer.al...

By any chance has anyone put together the R code for reproducing the Firefox Downloads Analytics Map in R? Here is a link to an example: http://www.knitwareblog.com/wp-content/uploads/2008/06/firefox-3-download-map.jpg I'll try to track down the inputs, but I am interested in putting the code to produce this or something similar. Any hints about the best packages, or references to R Blog