search for: lat_2

...t_point.Lat, lon1 = start_point.Long, lat2 = end_point.Lat, lon2 = end_point.Long) distance_slice = full_distance / numberOfIntervals for(i in 1 : numberOfIntervals) { distance = i * distance_slice ang_dist = distance / 6371.0 lon_1 = Radians(start_point.Long); lat_1 = Radians(start_point.Lat); lat_2 = Radians(end_point.Lat); lon_diff = Radians(end_point.Long - start_point.Long); # First calculate the bearing bearing = atan2( sin(lon_diff) * cos(lat_2), (cos(lat_1) * sin(lat_2)) - (sin(lat_1) * cos(lat_2) * cos(lon_diff)) ); # Then use the bearing and the start point to find the destination...

...on=Lambert conformal conic. # lcc.file <- system.file("extdata/ozone_lcc.ncf", package="M3") lcc.file <- "./ozone_lcc.nc" # unfortunate problem with raster::raster lcc.proj4 <- get.proj.info.M3(lcc.file) lcc.proj4 # debugging # [1] "+proj=lcc +lat_1=33 +lat_2=45 +lat_0=40 +lon_0=-97 +a=6370000 +b=6370000" # Note +lat_0=40 +lat_1=33 +lat_2=45 for maps::map at projection (below) lcc.crs <- sp::CRS(lcc.proj4) lcc.pdf <- "./ozone_lcc.pdf" # for output ## Read in variable with daily ozone. # o3.raster <- raster::raster(x=lcc.file, va...

Hi all, I'm faced with a problem applying the sp package: The projection argument in readShapePoly(Shapefile,proj4string="CRS class argument") e.g.: CRS("+proj=aea +lat_1=46 +lat_2=73 +lat_0=60 +lon_0=84 +x_0=0 +y_0=0 +ellps=clrk66 +units=m +no_defs") doesn't have any impact on the plotted object. I also tested the simple example: xy = cbind(x = 2 * runif(100) - 1, y = 2 * runif(100) - 1) plot(SpatialPoints(xy, proj4string = CRS("+proj=longlat")),xlim=c(...

...aphic coordinate system before I can do the analysis. I've been trying to use the project command in rgdal but keep getting an error message saying: Error in project(locationsMatrix, PROJECTION.OUT) : latitude or longitude exceeded limits ( PROJECTION.OUT <- "+proj=aea +lat_1=-18 +lat_2=-36 +lat_0=0 +lon_0=132 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs") I know that I'm using the right proj4 format (got it from the espy file in the Proj4 data directory) and I know that the datum is correct as I often use it in GIS (GDA94; Australia Albers). Bef...

Hi R users I need the map of France?s ? communes ? (towns) to build a map Is there a way to get it? More generally: How to do to get the map of a country and its different geographical levels? Best regards -- View this message in context: http://old.nabble.com/map-of-a-country-and-its-different-geographical-levels-tp26225645p26225645.html Sent from the R help mailing list archive at

Hi All, I have a netcdf gridded file with LCC projection. I can easily use image.plot to visualize it. However, as the axises are in X,Y, not Lat and Lon, I could not add state or country maps onto it (or lat lon information). I do have a grid2d file that describes the lat and lon for each (X,Y) grid, but the lat and lon are not regularly spaced, so I could not use image.plot. Does anyone know