similar to: map axis on projected shapefiles

Hi everyone, I wonder if anyone knows the reason why the outputs of the same reprojection in r and arcgis are different?. The magnitude of the change can be up to 40 km in the poles. Basically, I have a database of points equally separated by one degree over the globe. In ARCGIS, I am projecting the data in GCS-WGS-1984 and then reprojected it to Berhmann to ensure equal area

I imported the attached tiff file and converted the coordinate system to long lat and graphed it: californiatiff<- readGDAL("california1.tif") proj4string(californiatiff) > rasterprojection <- spTransform(californiatiff, CRS("+proj=longlat") however, when using the plot command for rasterprojection, I get a blob. I can see the outline of the state of california and

I would like to create a geographical grid to have a sort of a reference grid for my georeferenced survey data. The grid should be in a xy format, wgs1984 with a 0.025 degree, alternatively 10km, resolution covering -14 to -24 S and 24 to -34 E (Zimbabwe). Additionally I need to be able to export it as a .dbf Hope someone can help an R- novice ! [[alternative HTML version deleted]]

Ok, now we have the mapping of color names to color codes (see below) and conversion to rgb (something like Ben Bolker's function), but how to convert rgb to hsv? Thanks to Brian Ripley, Peter Dalgaard and Ben Bolker Details below Regards Jens etc/colors.big maps 455 names to rgb in S syntax etc/rgb.txt maps 657 names to rgb in C syntax, but unlike colors() it has mixed upper and

Dear R users, I'm currently trying to re-project a geotiff in another coordinate system. For instance, I have a tif image in UTM 19 zone which I would like to reproject into UTM 18. I was wondering if it was possible in R. Furthermore, I looked into 'rgdal' package, but I can't really find out if I'm doing the right thing. So far, here is what I'm doing: library(rgdal)

Hi all, We followed some books and sample codes and did some EMM estimation, only to find it won't be able to generate forecast. This is because in the stochastic volatility models we are estimating, the volatilities are latent variables, and we want to forecast 1-step ahead or h-step ahead volatilities. So it is nice to have the system estimated, but we couldn't get it to forecast at

Saludos, estoy intentando usar la función reprojectHDF() ( http://r-gis.net/?q=ModisDownload). #### source('ModisDownload.R') library(raster) library(RCurl) path<-"~/R/MODIS/Data/Test" input<-list.files(path,"*.hdf",all.files=T,recursive=T,full.names=T) input outname<-paste(substr(input[1],40,45),".tif",sep='') outname

In the manual page for prcomp(), it says that sdev is "the standard deviations of the principal components (i.e., the square roots of the eigenvalues of the covariance/correlation matrix, though the calculation is actually done with the singular values of the data matrix)." ?However, this is not what I'm finding. ?The values appear to be the standard deviations of a reprojection of

> > Hi, > > I used > GDALinfo("MOD13Q1.A2001049.h13v11.005.2007002215512.250m_16_days_EVI.tif") and > got the results: > > rows 10 > columns 11 > bands 1 > origin.x 150701.4 > origin.y 7744897 > res.x 250 > res.y 250 > ysign -1 > oblique.x 0 > oblique.y 0 > driver GTiff >

Dear all, I am using a shape file containing the information regarding the borders of the PARISHES of Austria. I created a plot with different colours for different percentages of child care institutions. Now I would like to add the information of the COUNTY boundaries to this plot which I have got in another shape file. I would like these boundaries to be shown with bold lines so that they can

Hi, I am new to R and I want to display longitude and latitude for a projected map. map.axes won't do it since it only works for unprojected map. Can anyone help me with this? Thanks, Sybil -- View this message in context: http://r.789695.n4.nabble.com/use-map-axes-for-a-projected-map-tp3302918p3302918.html Sent from the R help mailing list archive at Nabble.com.

Dear all, I'm performing a detrended correspondence analysis on vascular plant community data (296 species), and I have a question on the species scores projected in the ordination diagram. When I run a ordiplot all species are projected in the output graph, but I'd like to restrict the number of species plotted in the final graph. Some species are so rare in the data, that no relevant

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

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

Very helpful. One (possible) correction inline below... On Thu, 30 Mar 2006, Ray Brownrigg wrote: > > From: nhy303 at abdn.ac.uk > > > > 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) > > > >

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

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

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