[R-sig-Geo] UTM problem

2010-10-13 Thread Werner W.
Hi, this questions is not really about R but it would be great if somebody could point me to a forum or so where I can get help with this. I only know a few basics about GIS and have a problem with some UTM coordinates from a survey. Coordinates are given as, e.g. 3258843 32780 which are

Re: [R-sig-Geo] UTM problem

2010-10-13 Thread Sarah Goslee
Hi, I double-checked your results with an online conversion tool (http://www.rcn.montana.edu/resources/tools/coordinates.aspx) and got the same answer. That means that your R code is fine, and there's something wrong with the original data. If possible, you should check with the source. There

Re: [R-sig-Geo] UTM problem

2010-10-13 Thread Werner W.
Hi Sarah, many thanks for checking this. I looked at the documentation again and this is the info about the GPS receiver setting they used: POSITION FRMT: hddd.d MAP DATUM: WGS 84 CDI: +-0.25 UNITS: METRIC HEADINGS: GRID E000 DEGEES I'm not sure if this is any useful. Any idea what I

Re: [R-sig-Geo] UTM problem

2010-10-13 Thread Ageel Ibrahim Bushara
Hi, It seems the problem with your data. It seems there is inconsistency in your data. Here: 3258843 32780 You have 7 digits for x and 5 digits for y And here are the coordinates of Kampala: 459199 33160 on WGS84 You have 6 digits for x and 5 digits for y Ageel

Re: [R-sig-Geo] UTM problem

2010-10-13 Thread Sarah Goslee
It looks like Uganda is mostly but not entirely in 36N. It's possible the data switches zones, or that it keeps the same zone throughout and uses northing and/or easting values outside the usual range. Plotting the entire dataset as x, y coordinates would answer that question very quickly. Here's

[R-sig-Geo] Creating density heatmaps for geographical data

2010-10-13 Thread Karl Ove Hufthammer
Dear list members, Does anybody have suggestions for the best way of creating density heatmaps for geographical data? With ‘density heatmaps’ I’m not thinking of heatmaps as in R’s ‘heatmap’ function, but more along the lines of http://www.heatmapapi.com/ or

Re: [R-sig-Geo] UTM problem

2010-10-13 Thread Werner W.
Many thanks for your replies Ageel and Sarah! I still have no clue what to do with the data to make it usable. The northings are in the range of 26237 to 38149 which seem to make sense but the eastings are in the range of 3252446 to 3264111 and these don't seem to be in Uganda. If I use one of

Re: [R-sig-Geo] UTM problem

2010-10-13 Thread Barry Rowlingson
On Wed, Oct 13, 2010 at 4:14 PM, Werner W. pensterfuz...@yahoo.de wrote: Coordinates are given as, e.g. 3258843 32780 which are supposedly zone 36N on system WGS84 (according to the documentation). The point should be close to Kampala which has the coordinates 0.3°, 32.63° and UTM 36N

Re: [R-sig-Geo] reading .grd file in R?

2010-10-13 Thread Thiago Veloso
  My suggestion is to convert your .grd files to self-descriptive ones, for example netcdf. Considering that you have some experience with GrADS, a helpful tutorial is available on 

[R-sig-Geo] DBF file size

2010-10-13 Thread kapo coulibaly
I imported a DBF with more than 256 columns using read.dbf. But when i write it back out with write.dbf, only the first 256 columns are written out. Is there a way to export dbf files with more columns. (I actually need about 1200 columns) thanks [[alternative HTML version deleted]]

[R-sig-Geo] Convert asc grid to points

2010-10-13 Thread Arnaud Mosnier
Dear GeoR users, Is there a function to convert ascii grid (created with adehabitat kernel function) into a points format (ppp or others). Thanks Arnaud [[alternative HTML version deleted]] ___ R-sig-Geo mailing list

Re: [R-sig-Geo] DBF file size

2010-10-13 Thread Alex Mandel
On 10/13/2010 12:32 PM, kapo coulibaly wrote: I imported a DBF with more than 256 columns using read.dbf. But when i write it back out with write.dbf, only the first 256 columns are written out. Is there a way to export dbf files with more columns. (I actually need about 1200 columns)

Re: [R-sig-Geo] DBF file size

2010-10-13 Thread kapo coulibaly
I'm doing Groundwater Modeling and that is the appropriate input format for the GUI I'm using (Groundwater vistas) for what I'm trying to do. Funny enough I was able to create a dbf with 512 columns in ArcGIS but as I said I needed more. Thanks On Wed, Oct 13, 2010 at 4:11 PM, Alex Mandel

[R-sig-Geo] legend does not appear in all cases in a plot?

2010-10-13 Thread govindas
Hello all, Here I am with a plot issue - when I try to add legend it does not work at all times. Can someone tell me what's wrong with my code? the following legend does not work! plot(ts(tot_avg1[,1], freq=12, start=c(1981,12), end=c(2009,12)), ylim=c(min(tot_avg1[,1]), max(tot_avg1[,3])),