similar to: Map projections - converting latitude/longitude to eastings and northings

Hi, I have done the interpolation for my data and I was able to create the contours in multipanel with the help of Pascal. Now, I want to clip the contour with the shapefile. I want only the portion of contour to be displayed which falls inside the boundary of the shapefile. The data mydata.csv can be found on https://www.dropbox.com/s/khi7nv0160hi68p/mydata.csv The data for shapefile can be

Hello I'm working with a large hydrological data set stored in a netCDF format. The file stores x and y coordinates in the UTM projected coordinate system, yet when I use image to graphically display the z variable, the image is distorted in the sense that it does not plot the map in the correct spatial organization. I'm wondering if I need to define the projection of the netCDF file

Hello All, I have a discontinuous dataset and I used spplot to plot the points. I was wondering whether it is possible to show the legend of the plot as color ramp. I saw that on continuous data. If anyone could tell me that would be great. Another thing I would like to know is how can we create the horizontal legend. By default, the legend is vertical. I was able to shift the legend inside

ComRades: I am geeting the error message Error in ldots[[which(MET)]] : attempt to select less than one element when I try to fit the geostatistical model with the likfit() function of geoR. I have tried with old data for which likfit() successfully maximised the likelihood in previous versions of geoR, and yet the current version fails. I have tried in Windows Vista and Windows XP (I haven't

Hello R Community, I have imported a dataset which contain X Y coordinates and would like to recreate a shape file after some data analysis. What i have done is to import some taxlot data and join them based on some criteria. I want to check to see how well the joining went by reviewing the results in GIS. A couple things. I cant seem to import a shape file correctly using the maptools

I have imported a CSV file: rfishR <- read.csv(file="rfishR.csv",stringsAsFactors = FALSE, strip.white = TRUE, na.strings = c("NA","") ) attach(rfishR) When I call it up in R, it starts with line 2066 rather than 1 and some of the headers (used Headers = TRUE, too) are "masked"? Sample data loc lat lon datum water date obs net

Dear R-helpers, I get GPS readings with bug counts (bugs meaning insects in this case) made along rows in crop fields and use these to make maps of bug distribution. The GPS readings are not quite accurate enough for my purpose, so since I know what row each reading is made in, I adjust the latitudinal coordinate using: grd<-lm(lat~lon+Row,data)

hi, i have a longitude vector (x) a latitude vector (y) and a matrix of bathymetry (z) with the dimensions (x,y). I have already succeeded in plotting it with the image.plot (package 'field') and the contour functions. But now, I want to make a grid in order to extract easily the bathymetry corresponding to a couple of longitude, latitude coordinates. Do you know a function or a

Hi All, I''m having a problem with Rails truncating latitude and longitude values. In my database (mySQL) table, my lat/lng columns are defined as Decimal(25,20), but I get different precision values depending on what I do. If I post, via XML, a value of 37.436490000000006 to my controller, what gets inserted into the db table is 37.43648999999999900000. Then, when I fetch the record in a

I see several mapping libraries but cannot find how to do a simple task. Starting with a vector of zip codes (in the United States and in Canada) I would like to automatically lookup the corresponding coordinates. Can you suggest a package and function that I should turn to? I would also like to automatically pull in the altitude, and some climate data. -- Farrel Buchinsky [[alternative HTML

Hi Team, Need your guidance in building SPDF objects from Latitude, Longitude Information available. I am using package "plotGoogleMaps" which is really awesome but it requires SPDF objects to build the maps. I have a data frame which have Latitude and Longitude information and wanted to convert it in SPDF for making maps. Please help me to find suitable package or function to convert

Hi All, I have a data frame with three columns i.e., pc, lat, lon. The pc column is populated with list of postcodes, and I want to execute R command that can get me the lat and lon for the every item in the pc column and populate the respective lat and lon columns. Is there any package that could be used? Any help will be really appreciated. Many Thanks and Kind Regards -- Muhammad

I have a dataset giving traffic between pairs of ports, and the lat/lon of each port, roughly as follows: sample data as follows > df2VisitTrips[1:4,] Origin Destination NumberOfTrips OriginLatitude OriginLongitude DestinationLatitude DestinationLongitude 1 P1 P16 1 39.45965 -80.15633 40.76111 -79.54583 2 P1 P3

Hello, This is probably very simple but I am trying to select a region from the OIv2 SST dataset that crosses 0 longitude. The longitude variable is arranged from 0 to 360 and I get an error message when trying to go from 330 to 30 for example. The data are from a NOAA netCDF file. *temp = open.ncdf("c:/documents and settings/theoni/desktop/sstoi_v2.nc") print(temp) y1 = get.var.ncdf(

Hi all, I would like to know a way of calculating the distance between 2 points when I already have the longitude and latitude of the points. For example : Point 1 : 52? 9'54.00"N 4?25'8.40"E Point 2 : 52?27'46.80"N 4?33'18.00"E Distance between point 1 and point in km .... Is there any functions already available for this ? Any help will be much

Good day to you all, I have lightning data containing date, time, latitude and longitude. I hope that distribution of latitude and longitude will give number of lightning occurrence in a region. I have used factor function to sum up the number of events on latitude and longitude axis and saved as x and y. But when I tried to plot the two, I had and error message ( Error in image.default(x, y, z) :

Dear List, I am trying to?plot?bathymetry contours around the Hawaiian Islands using the package rgdal and PBSmapping.? I have run into a memory limit when trying to combine two fairly small objects using cbind().? I have increased the memory to 4GB, but am being told I can't allocate a vector of size 240 Kb.? I am running R 2.11.1 on a Dell Optiplex 760 with Windows XP.? I have pasted

Does anyone have any experience retrieving latitutde and longitude for an address from the Google Maps API? I'd like to have an R script that submits a street address, city, state, and zip code and returns the coordinates. So far, I've been submitting the coordinates from another program, then loading the coordinates in R and merging them back into the data frame I want to use.

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

I have a simple problem: I need to load a ERSI shapefile of US states and check whether or not a set of points are within the boundary of these states. I have the shapefile, I have the coordinates but I'm having a great deal of difficulty bringing the two together. The problem is the various GIS packages for R do not play well with each other. sp, shapefiles, maptools, etc all use different