similar to: Calculating distance between longitude,latitude of 2 points

Deemed colleagues I would appreciate your help with a sentence to calculate the linear distance between 2 geographical points (coordinates in UTM). In advance thnks for your attention, -- Ricardo Bandin Llanos rbandin@udec.cl Estudiante - Magíster Cs. m. Pesquerías Universidad de Concepción, Región del Bio-Bio, Chile Celular: (0056-41) 97949957 [[alternative HTML version deleted]]

Hi, I drew an ellipse with the package ellipse. Now I would like to know if a point is inside the ellipse. Is any R functions to do it without computing the equation of the ellipse manually? Thanks. For example, if I do "plot(ellipse(0.8), type = 'l')", I would like to know if (0,1) belongs to the drawn ellipse. Regards, Alain -- Alain Guillet Statistician and Computer

Hello all, I recently used the Vegan library quite extensively (in the context of text similarity assessment) on an Ubuntu 6.06 LTS system with R version 2.2.1 (2005-12-20 r36812). The Vegan lib is version 1.6-10. I hit on a problem yesterday, though, when trying to install R and Vegan on two further computers - one Windows XP and one further Ubuntu 6.06 machine, taking either R version 2.4.0

I have a question that I have been trying to figure out and I imagine there is a very simple answer to it. I am trying to use the distAB function in the clim.pact package I have a dataframe with 4 columns in which A is ref pt and B is site 1. longitude of A 2. Latitude of A 3. longitude of B 4. Latitude of B Problem is the columns are unequal in length column 1 and 2 are length=29 and 3 and 4

Hi All, I am unable to get the plot geomap in googleVis package. data is as follows > head(index.ret) country ytd 1 Argentina -10.18 2 Australia -3.42 3 Austria -2.70 4 Belgium 1.94 5 Brazil -7.16 6 Canada 0.56 > map1 = gvisGeoMap(index.ret,locationvar = 'country', numvar = 'ytd') > plot(map1) But it just displays a blank page, showing an

I'm having trouble understanding the output from as.windrose(). For one thing, data on a boundary between sectors seem to be left out of the counts. I assume that explains the missing point in the output below (angle 45). Shouldn't one side of each sector interval be open, to include values such as my 45 in the example? Also, why does the angle 180 in my input apparently not result in

Dear All, I would wish to make multiple plots and give title/ headings same time, is there a simpler/ tidier way compared to below(, especially the headings as they are missing)? See output attached. #####################################Making multiple Plots start######################## plot(mybrick9, y = 2, col = terrain.colors( length(seq(0,1, by = .2))-1),breaks= seq(0, 1, by = 0.2), axes =

Hello I will make an interpolation of data which represents azimuth direction ( angle from north in clockwise direction) values. But there is a problem. Say, for instance, while 1 and 359 indicate somewhat same direction, interpolation puts values in the range from 1 to 359. What can I do to solve the problem ? Anything you offer ? thanks in advance Ahmet Temiz General Directory of

Thanks Mike. My data has longitude and latitude coords and I used distAB {clim.pact} then mst {ape} to calculate my minimum spanning tree. The nodes are telecoms sites from all over Australia. My goal is to determine the minimum cost of linking them via cabling, and I'm starting by calculating the distance "as the crow flies", but will probably eventually need to calculate the

Is there packaged code to convert geographical coordinates (e.g., longitude, latitude, elevation) to Cartesian coordinates in 3-space? I can see how to do this using 1. a spherical-to-Cartesian conversion like pracma::sph2cart(tpr) http://cran.r-project.org/web/packages/pracma/ 2. a geographical-to-spherical conversion. This seems to involve (in roughly increasing order of difficulty or

Hi all, I have been searching all sorts of documentation, reference cards, cheat sheets but can't find why R's crossprod(A, B) which is identical to A%*%B does not produce the same as Matlabs cross(A, B) Supposedly both calculate the cross product, and say so, or where do I go wrong? R is only doing sums in the crossprod however, as indicated by (z <- crossprod(1:4)) # = sum(1 +

Dear all, I have an xts object , t.xts with 4 columns: "v1" "DD1" "v2" "DD2" and created a data frame : t <- as.data.frame(t.xts) I would like to extract data and create a new data frame for when the values in column DD1 falls between 0 and 30 and extract the corresponding v1 value. How can I do this? Thanks. -- View this message in context:

Hi, I am very new to R so these questions may seem simple! I have a huge 2 sets of data(matrix 5x20000++) in the following formats , for example "data.txt" and "data2.txt": Date Time X Y 03/03/1983 20:00 0.1 990 I would like to recreate a new matrix which filters through "data.txt" and

Hi all, My apologies for the confusion ... but I am re-naming this thread again because the proposal doesn't come from Ambisonia. As the founder of Ambisonia I have only been 'mediating' between some members of the ambisonic community ... so it would be wrong to associate the spec with Ambisonia. To clarify. This is what the spec proposes: - that Mapping = 1 means the contents are

On 2/28/07, Ivo Emanuel Gon?alves <justivo@gmail.com> wrote: > On 2/28/07, Ralph Giles <giles@xiph.org> wrote: > > Well, there are todo pages at wiki.xiph.org, but I meant more in the > > community folklore sense. My point is a roadmap doesn't help much unless > > there are people committed to making things happen. That's been the > > problem with a

On 2/28/07, Ivo Emanuel Gon?alves <justivo@gmail.com> wrote: > On 2/28/07, Ralph Giles <giles@xiph.org> wrote: > > Well, there are todo pages at wiki.xiph.org, but I meant more in the > > community folklore sense. My point is a roadmap doesn't help much unless > > there are people committed to making things happen. That's been the > > problem with a

Hello R list, I need to compute empirical variograms using data from a large geographic area (~10^6 km2). Although I could not find a specific reference, I assume that both geoR and gstat calculate distances among data points assuming points are on a flat surface (using the Pythagorean Theorem). Because the location of my data is large and located near the pole, assuming that latitude and

Hello R-helpers, What I have: I am working with vectors not matrice: #Basic equations A <-c(-20,-9,-2) x <-c(0.17,0.22,0.61) B <- crossprod(A,x) # R matrix multiplication works with vectors A%*%x==B # Is true... Question: If x is unknown and A and B are known, how do I solve for x in R? solve(A,B) won't work because A is not a square matrix solve(A,B) Error in solve.default(A,

I''m new to Ruby, Rails, and this list. I''ve tried searching for an answer but it all seems to be over my head. I have all the books but again they seem to leave me on my own to figure out the stuff that really matters. Where is the best place to find out the details? Exactly what are these much celebrated "naming conventions"; i.e. where do I find an exhaustive

I have spatial data on a sphere (the Earth) for which I would like to run an gls model assuming that the errors are autcorrelated, i.e. including a corSpatial correlation in the model specification. In this case the distance metric should be calculated on the sphere, therefore metric = "euclidean" in (for example) corSpher would be incorrect. I would be grateful for help on how to