search for: epsg

Dear R People: I'm sure that this is a very simple problem, but I have been wresting with it for some time. I have the following file that has the following one line: CRS("+init=epsg:28992") Fair enough. I scan it into R and get the following: > u [1] "CRS(\"+init=epsg:28992\")" > gsub(pattern='\"',replacement='"',x=u) [1] "CRS(\"+init=epsg:28992\")" I need to get rid of the extra quotation marks a...

...d = "C:\\GDAL_1.9.2_MAPSERVER_6.2.0\\SDKShell_original.bat") Setting environment for using the GDAL and MapServer tools. Hiding the OCI plugin library. D:\R_working_directory\Koordinatenproblem>> source(.trPaths[5], echo=TRUE, max.deparse.length=150) shell( cmd="ogr2ogr -s_srs EPSG:31254 -t_srs EPSG:25832 D:\\R_working_directory\\Koordinatenproblem\\Fliessgewaesser_25832.shp D:\\R_working_directory\\Koordinatenproblem\\Fliessgewaesser.shp", invisible = FALSE, intern = TRUE ) shell(cmd="ogr2ogr -s_srs EPSG:31254 -t_srs EPSG:25832 D:\\R_working_directory\\Koordinate...

...ot work for cokiriging it gives the same error message. require(gstat) require(automap) data_c <- read.csv("cu_s_data.csv", header=T) grid<- read.csv("grid.csv", header=T) coordinates(data_c) <- ~X+Y+Z coordinates(grid) <- ~x+y+z proj4string(grid)=CRS("+init=epsg:28992") proj4string(data_c)=CRS("+init=epsg:28992") g <- gstat(id = "Cu", formula =Cu~1,data = data_c) g <- gstat(g, id = "S", formula = S~1,data =data_c) g <- gstat(g, id = c("Cu", "S"), model = vgm(cov(data_c$Cu,data_c$S) , "Sph...

...mage 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 string and all that. Can anyone...

Dear All, I need to convert some Northing-Easting coordinates from the Gauss-Krüger system into latitude-longitude. Any suggestions on how to do it? Regards, Marco Marco Giannitrapani Statistical Consultant Tel: +44151373 5945 Email: Internet: http://www.shell.com [[alternative HTML version deleted]]

...left coordinates. from command line GDAL: sardinha$ gdalinfo dummy.tif Driver: GTiff/GeoTIFF Files: dummy.tif Size is 660, 420 Coordinate System is: GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.2572235630016, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4326"]] Origin = (-55.500000000000000,-14.500000000000000) Pixel Size = (0.008333333767951...

...d<-readShapePoly('D:/arcGIS/england boundary/england.shp') class(england) #define the projection proj4string(england)<-CRS('+proj=tmerc +lat_0=49 +lon_0=-2 +k=0.999601271625 +x_0=400000 +y_0=-100000 +ellps=airy +units=m +no_defs') # transform the map into the WGS84 projection (epsg:4326): england.wgs<-spTransform(england, CRS('+init=epsg:4326')) plot(england.wgs) #plot data over the map: quilt.plot(data$Longitude, data$Latitude, data$perc_per_year, add=TRUE) My problem is that I would like to be able to change how this data looks (not just 'grid squares',...

...site, I geocoded (using my geonames package) the locations of the CRAN sites returned by getCRANmirrors. I've mapped them here: http://www.maths.lancs.ac.uk/~rowlings/R/Cranography/ The geocoding may have got it wrong, so if CRAN site maintainers want to send me a real lat-long (in WGS 84 aka EPSG:4326) I'll update my map. Maybe eventually lat-long can go in the CRAN csv file distributed with R. Note this is just for fun. No warranties. Maybe I should use a little 'R' as a marker. Maybe I should get a life. Barry

...was trying to run the simple example that I found at the following web page: http://spatial-analyst.net/wiki/index.php?title=Export_maps_to_GE The lines that I can run are: data(meuse.grid) coordinates(meuse.grid) <- ~x+y gridded(meuse.grid) <- TRUE proj4string(meuse.grid) = CRS("+init=epsg:28992") # raster to polygon conversion; writeGDAL(meuse.grid["soil"], "meuse_soil.sdat", "SAGA", mvFlag=-99999) rsaga.geoprocessor(lib="shapes_grid", module=6, param=list(GRID="meuse_soil.sgrd", POLYGONS="meuse_soil....

...machine...) I recently made the switch to Ubuntu 10.04 (lucid), and found that unionSpatialPolygons() gives me the following error when trying to combine polygons within a shapefile: # Read in shapefile and merge polygons foo <- readShapePoly("foo.shp", proj4string = CRS("+init=epsg:2272")) IDs <- as.character(foo$ID) SP_out <- unionSpatialPolygons(foo, IDs) Error: isTRUE(gpclibPermitStatus()) is not TRUE I've seen some posts about a package "rgeos" which may have something to do with this error, and I downloaded the source from: http://r-forge.r-pro...

...******************************************************************** library(sp) Insitu = read.table("test.csv", skip=3,sep=",") Header = read.table("test.csv", skip=1,sep=",",nrows=1) names(Insitu)= as.character(as.matrix(Header)) CRS_4326 = CRS("+init=epsg:4326") Insitu_sp <- SpatialPointsDataFrame(cbind(Insitu[,"Boylam"],Insitu[,"Enlem"]), data=Insitu, proj4string=CRS_4326) plot(Insitu_sp[,"SPM"],col="black",pch=19,cex=0.725*sqrt(Insitu["SPM"]))

Hi I'm trying to analyse some data and need to set the geographic 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

...orum,   I am working in R 2.9.1 and I am trying to sample locations from a network file. Reading in nor plotting is a problem, however when I am trying to sample from the file I get the following message:   nwlim<-readShapeLines("C:/Limburg_nwshape", proj4string=CRS("+init=epsg:31300")) plot(nwlim)   randacc<-spsample(nwlim,n=1000,"random") Error in function (classes, fdef, mtable)  :   unable to find an inherited method for function "coordinates", for signature "numeric"     I found a similar error message in the arch...

...f:data-model" entity="pres:alice at atlanta.example.com"> <dm:device id="target123-1"> <gp:geopriv> <gp:location-info> <gml:location> <gml:Point srsName="urn:ogc:def:crs:EPSG::4326"> <gml:pos>32.86726 -97.16054</gml:pos> </gml:Point> </gml:location> </gp:location-info> <gp:usage-rules> <gbp:retransmission-allowed>false <...

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

Thank you Sarah and Mike for your explanations. My final objective is to produce maps (png image or any kind of extension I can import in LaTeX) where rainfall data are interpolated, using the Inverse Distance method or Kriging. My input file (pointfile.csv in the reported example) reports the station code, lat and long of the meteorological station and the rainfall value (which might be the

...as.numeric(as.vector(data_y_mt[i,])) it_data_mt <- idw(data_it_mt, sta_mt, grid, p = 1) grid <- data.frame(grid, valores = it_data_mt[ , 1]) coordinates(grid) <- ~x+y gridded(grid) <- T grid <- raster(grid, "valores") projection(grid) <- crs("+init=epsg:9377") path <- 'C:/Users/David Gomez/Desktop/TESIS/DATOS/RESULTS/' numeration <- as.character(i) name <- "MEAN_TEMP" exten <- ".tif" writeRaster(grid, paste(path, name, numeration, exten), overwrite = T)}* No entiendo muy bien a qué parte hace...

Dear all can you help me? i have this problem: i have a dataset in a text file in a matrix of 3 columns: x, y, z where x and y are the coordinates and Z are the mesurements. How can i do a variogram with R?

...+a947f6d+99unstable) ... Setting up evolution-data-server (3.26.3-4) ... Setting up qgis-providers (1:2.18.15+git20171224+a947f6d+99unstable) ... Synchronizing CRS database with GDAL/PROJ definitions. Load srs db from: /usr/share/qgis/resources/srs.db Loaded 503/504 from /usr/share/gdal/2.2/gcs.csv EPSG 6966: not imported EPSG 7082: not imported Loaded 4708/4786 from /usr/share/gdal/2.2/pcs.csv Loaded 192/193 from /usr/share/gdal/2.2/vertcs.csv Loaded 156/157 from /usr/share/gdal/2.2/compdcs.csv Loaded 163/164 from /usr/share/gdal/2.2/geoccs.csv Loading epsg.wkt Loading esri_extra.wkt 5722 WKTs lo...

I am running Debian testing and updated my system this morning which included a new kernel: 4.14.0.2-amd64 After this update, my system froze twice -- both times the exact same behavior. In both situations, the system froze while I was resizing a window. The mouse cursor froze, clock stopped, num lock stopped working. I didn't check if the system responded to pings from a remote machine,