similar to: from long/lat to UTM

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

I'm hoping someone can help me. I have downloaded the proj4R.zip and under my version of R (2.3.1) I install the package from local zip file. This worked great. I then type library(proj4R) to load the library and I get the error: Error in library(proj4R) : 'proj4R' is not a valid package -- installed < 2.0.0? I have read through the install documentation and have downloaded and

Hi, I'm trying to use transform my mercator locations into utm but I'm doing something wrong because only x is transformed ... see > xy.sp SpatialPoints: loni lati [1,] -8.85 38.16 [2,] -9.19 37.99 [3,] -9.11 37.97 [4,] -9.06 38.15 [5,] -9.03 37.87 [6,] -9.14 37.81 [7,] -9.09 37.70 [8,] -8.95 37.45 [9,] -9.17 37.37 [10,] -9.00 37.34 Coordinate Reference

Hi: I have two vectors of data, x and y and I want to get the "polynomial" expansion of (x+y)^p with any integer power p in R. Suppose p=2, then I want a matrix of five vectors, namely, x y x^2 y^2 x*y. The coefficient of the polynomial is not needed. I can write it manully if p is small. But I want it in the case of p=10 or even bigger, is there any function in R can do that

Hi: I have some problems when I use the function "solve" function in a loop. In the following code, I have a diagonal martix "ttt" whose elements change in every iteration in a loop. I defined a "dpoMatrix"class before the loop so I do not need to define this class every time in the loop. The reason is to save some computing time. The code is below. The inverse

Dear R users, Could somebody tell me why the grey color ramp is repeated in this graph, eventhough the ramp values go from 0 to 1? I must be missing something obvious, but I can not see it! z <- c(0.064329041,0.117243316,0.161565116,0.19923015,0.231642175,0.259835539,0.284571226, 0.038507288,0.094184749,0.140959431,0.180803984,0.215159105,0.245096084,0.271412845,

Hi R helpers, I tried to convert a list of LatLong coordinates (in DD format) into UTM zone 11 NAD 27. I first tried this from PBSmapping: library(PBSmapping) LatLong<-cbind(c(56.85359, 56.85478),c(-118.4109, -118.4035)) colnames(LatLong)<-c("X","Y") attr(LatLong, "projection") <- "UTM" attr(LatLong, "zone") <- 11

Hi R-Users, I have plotted a region whose polygon coordinates are given in shp format ED50 UTM (zone=30) ) using "readShapePoly" in library(maptools). Now I need to plot a set of points in that region (my.dataframe, with X and Y geographic coordinates), which have been read using GPS in Longitud-Latitud form (using WGS84 system), so I first need to convert these Longitud-Latitud data

Dear all, I have been trying to convert coordinates from longitude/latitude to utm but I got an error. As soon as the longitude coordinate is greater than 90, I get the folloowing error message: "error in pj_transform: latitude or longitude exceeded limits" Here is what I did: SP<-SpatialPoints(cbind(126.59,-14.30),proj4string=CRS("+proj=longlat"))

Dear R users, The Sweave code below runs fine, as it is. However, an error occurs when the line 'library(xtable)' is uncommented: Error: chunk 1 Error in "label<-"(`*tmp*`, value = "month") : no applicable method for "label<-" Is anybody aware of this and knows a workaround? Thanks, Sander. *******************

Dear all, I am trying to project my LongLat-maps to a plane. The ultimate purpose is to do a search of points in vicinity of other points using overlay-commands (sp) with radius in km. I am applying spTransform (package rgdal) and it gives my some curious results. An example. Let's take a point lying somewhere in Germany, zone=32U x <- 8.968735 y <- 49.454735 After conversion I

Hi -- I have searched the documentation and archives on graphing capabilities in R for the past couple of hours, but I haven't been able to find anything directly related to my problem. I want to create a plot with several lines displayed on it. I want each line to be displayed in a different color, so it seems the only way to do this is to create a plot and then use the lines() function.

Hi, is there a function in R to convert data read with read.shape and which is originally in UTM coordinates into longitude / latitude coordinates? I found the convUL() function from the PBSmapping package but I have no idea how I could apply that to the read.shape object. Many thanks, Werner __________________________________________________ Do sragenden Schutz gegen Massenmails.

A non-text attachment was scrubbed... Name: not available Type: text/enriched Size: 905 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/20010420/61ada248/attachment.bin

Dear all, We have a large data set with temperature data for weather stations across the globe (15000 stations). For each station, we need to calculate the number of days a certain temperature is exceeded. So far we used the following S code, where mat88 is a matrix containing rows of 365 daily temperatures for each of 15000 weather stations: m <- 37 n <- 2 outmat88 <- matrix(0,

hello, to use the function plot3d, i should use the package R.basic! plot3d {R.basic} If people know exactly a site to load this package, please give me the URL! Thanks Sabine --------------------------------- [[alternative HTML version deleted]]

Does MGCP transverse NAT? Seeing as the only decent yet cheap IP phone is the Swissvoice, it would be rather helpful.

Dear all, I have a data frame containing the results of an experiment. Like this: a<-seq(1,4,by=1) b<-seq(1,2,by=1) test<-expand.grid(b,a,a) colnames(test)<-c("replicates","bins", "groups") test$abc <- rnorm(32) test$def <- rnorm(32) test$ghi <- rnorm(32) test The following code snippet aggregates the data for one variable and then draws a

Hello all, I noticed that both Windows and OS X version of R have a GUI (Rconsole). Is there a GUI for Linux? I'm running Debian on which the CLI for R works just fine. Regards, - Robert

Dear expe-R-ts, I try some test data for a factorAnalysis (resp. pca) in the sense of Prof. Ripley's MASS ? 11.1, p. 330 ff., just to prepare myself for an analysis of my own empirical data using R (instead of SPSS). 1. the data. ## The test data is (from the book of Backhaus et al.: Multivariate ## Analysemethoden. Springer 2000 [9th ed.], p. 300 ff):