similar to: Bubble plot on shapefile - projection issues?

Dear all: I have a shapefile in UTM coordinate system and I would like to transform it to Lat/Log coordinates (WSG84). The package PBSmapping contains function convUL to transform between the two coordinate systems when data is in the form of a data frame with attributes specifying the coordinate system. However, when shapefiles are imported using function read.shape (package maptools), a list

Hi, ? I am trying to write a shapefile with projection. I have my data in a data.frame called try and consists in xy coordinates and a numerical attribute value z1. ? Libraries loaded are: sp, rgdal, raster, maptools ? head(try) ???????? x?????? ??? y??? ? ? ? ? ? ? ??? z1 1 610237.1 ???????? 3375751 ???????????? 8.221 2 610236.1 ???????? 3375750 ???????????? 8.153 3 610236.1 ????????

Hello I wonder if anybody can help, I am using the package adehabitatHR to estimate the potential distribution of a species using the command "domain" In the example given in the AdehabitatHS manual a map containing elevation information is loaded (class= spatial pixels data frame) as well as the GPS points of the animals being tracked, these are then plotted on each other and

I have several large shapefiles that I need to use to draw a map. I want to simplify the files and then save the output so that I can use them again. The methods that I have tried do not allow me to save the simplified output as a shapefile. I have tried: thinnedSpatialPoly in the package maptools, but I understand that you cannot save the output as a shapefile because the topographical features

I'm new to R and I'm working on point shapefiles. Is there a way that you could interpolate a shapefile via kriging in R using an attribute? All examples on the internet are using txt files and CSVs. Thanks a lot. [[alternative HTML version deleted]]

Hello I'm working with NetCDF files in an animation and am trying to superimpose a shapefile on the image as it is generates the html pages. If I take out the lines for the shapefile, it works correctly, however I'm having difficultly including the shapefile. If anyone has any ideas I would greatly appreciate the assistance. R Version 2.15.0 (2012-03-30) Platform i386-pc-mingw32/i386

My sincerest apologies, as this is a very elementary problem, but I have searched through archives and FAQs and on the 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

Hi, I have done the interpolation for my data and I was able to create the contours in multipanel with the help of Pascal. Now, I want to clip the contour with the shapefile. I want only the portion of contour to be displayed which falls inside the boundary of the shapefile. The data mydata.csv can be found on https://www.dropbox.com/s/khi7nv0160hi68p/mydata.csv The data for shapefile can be

Hi all, who knows how to get an ESRI Shapefile for the NUTS-2 Regions of the enlarged European Union? Particularly I want to draw maps of Germany, Poland, Czech Republik, Hungary and Austria. I've found Shapefiles for the US, Russia and other countries elsewhere in the web, but for Europe it seems really difficult. With many thanks for any hint Albrecht

I am trying to estimate home range size using the plug-in method with kernel density estimation in the kernel smoothing (ks) package. Unless there is another way I am not familiar with, in order to calculate spatial area under the space I need to convert my kde () object into a spatial object somehow in order to calculate its spatial area. Could someone demonstrate how this might be done? --

I would like to create a voronoi map that can show me point values like the ones in colum street_cle. I was able to do it in tileplot(latticeExtra) like the ones attached, and with the following commands: xy=list(x=x2$X, y=x2$Y) xy$street_cle <- x2$street_cle tileplot(street_cle~ x*y, xy) But this is a trellis. Is it possible to be converted to a shapefile or image so I can plot it using

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

Hi, I am working on a library of map functions for R and I have started with code to read shapefiles (the format is not great but it is ubiquitous). I am not a windows person, but where I am working I don't have access to any UNIX machines easily. So here's the problem I wrote C code to read the information out of the shapefile header and I compiled it using the script given in Ripley and

Dear R users, I would like to know how to indentify the polygons that are located on the border of a map (i.e.shapefile). Do you have any suggestion on how to do it? Thank you very much, Leo Monasterio. [[alternative HTML version deleted]]

Hi everybody: Could anybody help how I solve the next problem?. I'm doing interpolation maps of tropospheric ozone of my region, and after create it using IDW, and kriging methods, I want from shapefiles (*.shx, *.shp, *.dbf, *.sbx & *.sbn ) create contour over the interpolation maps. Could anybody tell me how do it?. Thanks in advance. ###################################### Antonio

Scotland is the second feature in the UK data, so get it and split this one MULTIPOLYGON feature into individual POLYGONS scot = st_cast(the_uk$geometry[2],"POLYGON") # which is the largest polygon? which.max(st_area(scot)) [1] 1 # the first one. ok... plot(scot[[1]]) # mainland # add the rest of the islands for context, in grey, maybe to show they're outside our study area:

I'm not sure if this is what you want. But the function coordinates() in sp package gives you the coordinates of SpatialObjects. Regards. Rodrigo. 2010/6/9 Nikhil Kaza <nikhil.list@gmail.com> > You need to execute gpclibPermit() to enable gpclib. > > library(maptools) should have issued a warning to that effect. > > > Nikhil Kaza > Asst. Professor, > City and

Hopefully the attachment will make it this time... Hi: I am practicing with the attached shapefile and was wondering if I can get some help. Haven't used 'rgdal' and 'maptools' much but it appears to be a great way bring map data into R. Please take a look at the comments and let me know if I need to explain better what I am trying to accomplish. library(rgdal)

Is it possible to convert an hexbin object like hexbin(rep(1:100,each=100),rep(1:100,100),xbins=10) and convert it to a vectorial (shp) format? Paulo ---------------------------------------------------------------------------- ------------------------------------------------- This electronic mail transmission including any attachment hereof, contains information that is private and confidential,

Hello, I have an ESRI ArGis geo database directory which I would like to convert to geojson or a shape file (or anything else that I can read into R). Unfortunately that does not work out of the box with rgdal, because it does not come with the fileGDB or openfileGDB driver. I could successfully install gdal and the fileGDB driver/extension, but it seems that i can use gdal only to convert the