similar to: Plotting shapefiles on existing maps

Dear All, This is probably a very basic question but: I have plotted a map of the Barents Sea and surrounding coastline using: map('worldHires',ylim=c(50,85),xlim=c(5,65),fill=T,resolution=0) map.axes() map.scale(x=30,metric=T) Next, I imported a shapefile with depth contours for the sea: contours<-read.shape("D://My Documents/BarentsSea.shp",dbf.data=T) (This is in

Dear R-users, I have imported a shapefile with depth contours for a sea: depths<-read.shape("D://My Documents/BarentsSea.shp",dbf.data=T) (This is in mercator projection) **Is there a way to convert this shapefile into a format that it may be plotted on a contour plot?** I wish to add these contours onto a map (already coded using 'maps' package) to map the sea contours

On Thu, 24 May 2018 at 18:41 Mark R Payne <markpayneatwork at gmail.com> wrote: Hi, I seem to be having a problem adding the following two raster objects together - one is a rasterLayer, the other is a rasterBrick. The extent, resolution, and origin are the same, so according to my understand it should work. The objects look like so: > obs.clim class : RasterLayer dimensions : 60, 200,

This is now fixed in development on RForge, you can try it out by installing from there, or from the Github mirror with devtools::install_github("rforge/raster/pkg/raster"). (To get fixes into raster email the maintainer directly - you might not get a response but it'll be addressed). Cheers, Mike. On Thu, 24 May 2018 at 20:08 Michael Sumner <mdsumner at gmail.com> wrote:

Hi I have a netCDF file of volumetric soil water content at four different soil layers and want to convert each soil layer in the netCDF file to a GeoTIFF layer. This code converts the netCDF file to one GeoTIFF layer, i.e. unclear which soil depth. file.nc <- "C:/Soil_Weather_data/Agro/VMC21/VMC21.nc" file.tiff <- "C:/Soil_Weather_data/Agro /VMC21/VMC21" importnetcdf

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 ????????

Hi Ahmed, When reading from a ncdf file you can use the 'varname', 'lvar' and 'level' arguments - see the 'Details' section in the docs https://www.rdocumentation.org/packages/raster/versions/2.6-7/topics/raster <https://www.rdocumentation.org/packages/raster/versions/2.6-7/topics/raster> We can't tell what is in the ncdf file from what you report other

I just uploaded two packages to CRAN. shapefiles_0.1.tar.gz - functions to read and write ESRI shapefiles (including dbfs) emme2_0.1.tar.gz - functions to read binary data from an EMME/2 databank data (EMME/2 is a transportation modeling program) Please let me know if you find any bugs or have some suggestions. Thanks. Regards, Benjamin Stabler Transportation Planning Analysis Unit Oregon

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 !

Dear List, I have examined the sp, maptools, and shapefiles packages (Windows) but cannot seem to figure out a solution for writing shapefiles for 'curvy shapes'. The scenario is that I generate a series of polygons, circles, and ellipses on a plot and would like to convert them to shapefiles. Since the circles and ellipses are do not return a coordinate list during the plotting

Hi, I have overlayed 2 projected shapefiles using the "plot" function. When I plot them the numbers next to the x-axis ticks are printed in scientific format (e.g. -4e+05, -3e+05,...,), and the problem is these dont match up to the projected point locations I also wish to overlay from my dataset. For example, according to the x-axis ticks the prime meridian (0?) is in totally the

Now available on CRAN is a revised version of the shapefiles package for reading and writing shapefiles in R. New additions, courtesy of others, include the ability to convert a simple R data frame of points, polylines or polygons to a shp format list, which can then be written out to a shapefile with write.shp. There is also a function to convert the read.shp shp format list to a simple data

I am attempting to write what MATLAB does as : [S,A]= shaperead() which returns an N by 1 geographic array structure S containing geometric information, and a paralell N by 1 attribute structure array, A containing feature attribute information. How can this be done with R code I tried using shp<-readShapePoints() but this doesnt return the same structure that I am looking for. Please

Hi, I'm a SAS veteran and R newbie. Can anyone point me to where I can best find examples on reading in shapefiles with R, and then merging in measurement info for the different polygons? I'm struggling a bit finding a good toehold in the documentation, as I am not yet familiar with navigating it yet. Thanks, Matt [[alternative HTML version deleted]]

I'm hoping that someone could guide me in how to extract data from shapefiles. I want to extract data from a shapefile (classed as "SpatialPolygonsDataFrame") and more specifically, the data is contained within the slot called "coords" from the class "polygons" within this file and despite my best efforts (and much reading!), I'm at a lost. Many thanks, Kate

Hi, Is there a way, preferably with R, to read shapefiles and transform them in a format that I could then use with package PBSmapping? I have been able to read such files into R with maptools' read.shape and plot it with plot.Map, but I'd like to bring the data to PBSmapping and plot from there. I also looked at the package shapefile, but it does not seem to do what I want

Dear All, I'm new with mapping and using shapefiles in R. I use the version: R 3.2.1 GUI 1.66 with Snow Leopard Mac OS (6956) I want to build a map of the countries of Africa, in order to add points of collecting species of beetles, for later publishing. library(maptools) library(mapdata) When I put the extreme coordinates of Africa, the map began only at longitude 0? and thus lacks the

Hi all, I'm looking for zip code shapefiles for Italy and France.Does anyone know where I might be able to find these? thank you! [[alternative HTML version deleted]]

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

Dear List, I'm trying to import a multiple polygon shapefile in splancs. I've found the following command but it seems to work only for single polygons: Name_of_splancs_polygon<-getPolygonCoordsSlot (getPolygonsPolygonsSlot(getSpPpolygonsSlot(name_of_sp_polygon) [[1]]) [[1]]) I've tried to change the numbers but it seams to choose polygons instead of importing them