search for: proj4

...to get a 1km by lkm grid raster image using my csv data. If I call latitude=a, longitude=b and preciptation=c. a<-(1,2,3,4,5) b<-(6,7,8,9,10) c<-(10,20, 30,40, 50) Then I found an example in r help which goes like pts = read.table("file.csv",......) library(sp) library(rgdal) proj4string(pts)=CRS("+init=epsg:4326") # set it to lat-long pts = spTransform(pts,CRS("insert your proj4 string here")) gridded(pts) = TRUE r = raster(pts) projection(r) = CRS("insert your proj4 string here") Because I am new to R, I have no idea what to put into the proj4...

An embedded and charset-unspecified text was scrubbed... Name: n?o dispon?vel URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20080916/9ba9430c/attachment.pl>

...ge 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). Before using the project command I have a series of commands to make my raw input (a csv file) into a matrix with only the two columns of co-ordi...

I am having trouble getting the rgal package to work on Mac OS X. I suspect it has something to do with the link between the rgal package and the associated programs GDAL and PROJ4. Here is the message I get when I attempt to load rgal > library("rgdal") Error in dyn.load(file, DLLpath = DLLpath, ...) : unable to load shared object '/Library/Frameworks/R.framework/Versions/2.15/Resources/library/rgdal/libs/x86_64/rgdal.so': dlopen(/Library/Framewor...

...quot;) # http://cran.r-project.org/web/packages/rasterVis/ ## Use an example file with projection=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" # f...

Hello, We (two mac users) have been attempting to install rgdal from "http://www.kyngchaos.com/software:frameworks", given that it is not available as a binary on the CRAN (binaries) of the "Package Installer". I have also tried to solve this problem by looking on the net for an old question, and though I have found it, the answers do not help very much because our porgraming

I run R on a Fedora 23 64bit , I encounter this problem *WARNING: no proj_defs.dat in PROJ.4 shared files Linking to sp version: 1.2-2 * when trying to load library("rgdal", lib.loc="~/R/x86_64-redhat-linux-gnu-library/3.2") and also ogr . can any one help me in fixing this issue please Thanks in advance Irenge B. Jules MSC student @ Liverpool University email :

...ordinates 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 with rgdal or proj4 routines first before I send it to the graphics device. My code is as follows: q1_2001 <- open.ncdf("H:\\SKF_DESKTOP FILES\\My Documents\\EDEN\\EDEN\\Surfaces\\2000_q1.nc", readunlimi=FALSE) #opens ncdf file for reading wat.data2000q1 <- get.var.ncdf(q1_2001, verbose=FALSE )...

...o the difficult thorough access in the study area. I could use a code to create the points throughout the area, but I need to create a code to make the R only comes up with points close to this road. I've been using the packages "raster", "shapefile", "clhs", "proj4" and "rgdal" and I have the shapefiles of the road, the raster of distance from the road, and so forth. Thanks for any kind of help. Sergio -- View this message in context: http://r.789695.n4.nabble.com/Conditioned-Latin-Hypercube-Sampling-within-determined-distance-tp4633974.htm...

Dear list, I have to compare two digital elevation models in raster format (geotiff). I then have to calculate the differences in altitude for each cell and make some statistics (basic as mean, median, std, range but also more "advanced" as RMSE) on that. I do not know very much how to proceed: 1) is it possible to import the geotiff in R? If so with which package? if not which is the

...y value : 560615.564407584, 4698904.71208346 Please observe that both files have identical projections. p1 <- gIntersection(buff, pln) Warning message: In RGEOSBinTopoFunc(spgeom1, spgeom2, byid, id, drop_lower_td, unaryUnion_if_byid_false, : spgeom1 and spgeom2 have different proj4 strings This warning is baffling since both files have same projection: crs(buff) CRS arguments: +proj=utm +zone=16 +datum=WGS84 +units=m +no_defs crs(pln) CRS arguments: +proj=utm +zone=16 +datum=WGS84 +units=m +no_defs Probably i am doing something wrong but for the life of me i cannot see i...

...ent by reviewing the results in GIS. A couple things. I cant seem to import a shape file correctly using the maptools package and the readShapeSpatial. I have tried Building=file("data/input/BuildingShape/Building.shp") Bldg<-readShapeSpatial(fn=data/input/BuildingShape/Building,proj4string=NAD83) #---------------------- Bldg<-readShapeSpatial(data/input/BuildingShape/Building,proj4string=NAD83) #--------------- Building=file("data/input/BuildingShape/Building.shp") Bldg<-readShapeSpatial(Building,proj4string=NAD83) I know i am mis interpreting the documentation...

Dear groups: I have tried to install rgdal in my Ubuntu, but it is not as straight as installing "sp". It demand the pre-installed GDAL and proj.4. I searched the source code in http://trac.osgeo.org/proj/ and http://trac.osgeo.org/gdal/wiki/DownloadSource. but don't know how to compile them. Who will give me a hint? Thanks a lot. Ethan Chang Centers for Disease Control &

I am using R1.2.0 on Windows NT4.0 and am looking for information on drawing geographical maps. [The things I'd like to do at the moment are rather simple, just like plotting the values of some variable (e.g. a SMR) on a given map.] Searching the r-help mailing list archives, it seems this topic comes up from time to time. For example, in November 1999 Ray Brownrigg writes about an "R

Is it possible to suppress compiling for the x64 arch when doing 'Rcmd build <mypkg>? Such that the pre-compiled package for binary distribution only contains a 'lib/i386/mypkg.dll'. For CHECK (and also INSTALL) there is a '--no-multiarch' argument, but I didn't find something similar for build. Otherwise, I suppose, I could - handle the problem in Makevars or

Hi there, I'd like to plot some maps. Is it possible for me to use SAS map data in R? Thank you. Best, Junjie _________________________________________________________________ your life. [[alternative HTML version deleted]]

Dear R-users, I would like to process some spectroscopic data with R, and I was hoping some people might have some example code on how to do this. I would like to be able to do the following things: * Detect outlier-spectra -> This can be done by using scoreplot from the pls package * Determine the range of the spectrum to be used -> For this, I should be able to calculate the regression

...th = /opt/gmp read only = No [localopt] comment = Witham Local Specific Applications path = /local/opt read only = No [proj2] comment = Witham Project Area 2 path = /export/project2 read only = No [proj3] comment = Witham Project Area 3 path = /export/project3 read only = No [proj4] comment = Witham Project Area 4 path = /export/project4 read only = No [proj5] comment = Witham Project Area 5 path = /export/project5 read only = No [pwise] comment = ProjectWise Data path = /export/pwise read only = No browseable = No [indexes] comment = Technical Index...

#################################################################################### We have used maptools to construct state, county, township, census-tract, and zipcode level R maps with an Albers Conical Equal Area projection. We would like to be able to plot the location of weather stations or other point locations on the maps. The data the point locations are in latitude-longitude units

R Help, I posted a question on StackOverflow yesterday regarding an issue I've been having with the RgoogleMaps packages' displaying of axes. ?Here is the text of that submission: http://stackoverflow.com/questions/6258408/rgooglemaps-axes "I can't find any documentation of the following problem I'm having with the axis labels in RGoogleMaps: library(RgoogleMaps) datas