search for: sptransform

...y tratando de cortar mi mapa de (poner limites en UTM) en un lugar definido como mi area de estudio (en este caso el sur de chile). Pero creo no estar usando bien la función CRS ponendo bien los limites requeridos. > study_area <- readRDS("CHL_adm0.rds") > study_area_UTM <- spTransform(study_area, CRS("+proj=utm +zone=19 +datum=WGS84")) > study_area_UTM <- spTransform(study_area_UTM, CRS( + paste("+x_0=-2000000.0 +y_0=-500000.0 +ellps=GRS80 +units=us-ft +no_defs"))) Error in spTransform(study_area_UTM, CRS(paste("+x_0=-2000000.0 +y_0=-500000.0 +ellp...

Dear all, I am trying to project my LongLat-maps to a plane. The ultimate purpose is to do a search of points in vicinity of other points using overlay-commands (sp) with radius in km. I am applying spTransform (package rgdal) and it gives my some curious results. An example. Let's take a point lying somewhere in Germany, zone=32U x <- 8.968735 y <- 49.454735 After conversion I sould get something like Easting: 426858 (km) Northing: 5427937 (km) sp1 <- SpatialPoints(matrix(c(x,y), ncol=...

Salutations, Long story short: I need to create an SPDF to superimpose on a Google map. I have .shp file which has a long list of locations listed in (what I guess must be) UTM. I need to transform these into lon, lat. My use of spTransform, however, keeps giving me the wrong results. For what it's worth, it only seems to be messing up the latitude. My final line of code looks like this: pointsTransformed <- spTransform(pointsutm, CRS("+proj=longlat +datum=WGS84 +zone=36S")) The line before looks like: pointsutm &l...

...rror. 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")) coordinatesUTM<-spTransform(SP,CRS("+proj=utm")) Error in spTransform(SP, CRS("+proj=utm")) : error in pj_transform: latitude or longitude exceeded limits Can someone explain me what I did wrong? I am using coordinates of Australia so the range of my longitudes is around [110,150] and the one of my lat...

...w, ncol, ncell) resolution : 0.2603037, 0.2916659 (x, y) extent : -180, 180, -85.4581, 85.4581 (xmin, xmax, ymin, ymax) coord. ref. : +proj=longlat +datum=WGS84 names : Retrieved.soil.moisture.value zvar : Soil_Moisture a <- spTransform(d, CRS ("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")) but I got this error: Error in function (classes, fdef, mtable) : unable to find an inherited method for function ?spTransform? for signature ?"RasterLayer", "CRS"?...

Hi, is there a function in R to convert data read with read.shape and which is originally in UTM coordinates into longitude / latitude coordinates? I found the convUL() function from the PBSmapping package but I have no idea how I could apply that to the read.shape object. Many thanks, Werner __________________________________________________ Do sragenden Schutz gegen Massenmails.

I imported the attached tiff file and converted the coordinate system to long lat and graphed it: californiatiff<- readGDAL("california1.tif") proj4string(californiatiff) > rasterprojection <- spTransform(californiatiff, CRS("+proj=longlat") however, when using the plot command for rasterprojection, I get a blob. I can see the outline of the state of california and nevada, but rather than being able to see the geographic features, it is just a monocolor blob. My suspicion is that since...

Hi, I want to use map("state") and have the ecoregion shape (please see link) file projected onto this. ftp://ftp.epa.gov/wed/ecoregions/cec_na/NA_CEC_Eco_Level3.zip ftp://ftp.epa.gov/wed/ecoregions/cec_na/NA_CEC_Eco_Level3.zip Could someone please show me how; I have never messed with this sort of thing. Thanks. -- View this message in context:

..., 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 help me on this? If there is a different way of do...

I am trying to set up a Grass project and need to set up the region so that I can view the map. I can look at a map and find the lat/lon, but the map projection is in UTM NAD38 WGS84 and I need to set the eastings and northings. Is there a package that will help me calculate this in R. thanks -- Stephen Sefick Let's not spend our time and resources thinking about things that are so little

> If you use the spatial objects provided by the > sp-package (http://cran.r-project.org/web/packages/sp/vignettes/sp.pdf) > you transform your data to other projections using the spTransform package. Thus you will need the rgdal package in complement (it actually includes spTransform). This function is extremely convenient: you can manage coordinate transformations extremely easily for common systems (WGS84, UTM) within the R environment.

...into UTM 18. I was wondering if it was possible in R. Furthermore, I looked into 'rgdal' package, but I can't really find out if I'm doing the right thing. So far, here is what I'm doing: library(rgdal) myTiff = 'L5014028_02819840617_B10.TIF' a = readGDAL(myTiff) aa = spTransform(a, CRS("+proj=utm +zone=18 +datum=WGS84")) However, after few seconds, the memory is full and R crashes. I used other proprietary tools to reproject and it takes only like 10 seconds to reproject. Therefore, I suspect I'm not using the 'rgdal' package correctly. Regards, Phi...

Hola grupo, estoy siguiendo una gu?a de la librer?a rLandsat que me la he descargado de: devtools::install_github("socialcopsdev/rLandsat") Y tras hacer los siguiente (obviamente tengo me he registrado previamente en la api correspondiente): product_id = c("LC08_L1TP_145049_20180301_20180308_01_T1", "LC08_L1TP_145049_20170330_20170414_01_T1",

...e data in GCS-WGS-1984 and then reprojected it to Berhmann to ensure equal area distribution of the points. In R, I am using: spPoint <- SpatialPoints(coords=coordinates(Data),proj4string=CRS("+proj=longlat +datum=WGS84")) and then reprojecting it to Berhmann with: spPointReprj=spTransform(Data,CRS("+proj=cea +lon_0=0 +lat_ts=30 +x_0=0 +y_0=0 +ellps=WGS84")) If I put the two outputs of the reprojections in the same map, they are off by few meters in the tropics by up to 40km in the poles. I decided to investigate how the reprojections are done and my calculations ar...

...w that the cause of this error is because somehow, the package is not being loaded in the NAMESPACE, but it should be there because I am using LazyData. That is the code for the PlotBgMap function: PlotBgMap <-function( traj, ... ) { hySplitProj <- CRS(proj4string(traj)) canada <- spTransform(canada.map, hySplitProj) plot(canada.map, border="white", col="lightgrey", ... ) } I'd like to know if I am missing something in the Checking process, or if I should do anything else to fix this error. Thanks. -- View this message in context: http://r.789695.n4.nabbl...

...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', but circles, etc) which I can'...

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

...r?a ayudar? > > > xmn=min(table_2[,1]); xmx=max(table_2[,1]) > ymn=min(table_2[,2]); ymx=max(table_2[,2]) > > b <- data.frame(cbind(table_2$x, table_2$y)) > coordinates(b) = ~X1 + X2 > > proj4string(b) <-CRS("+proj=utm +zone=14 +datum=WGS84") > b1 <- spTransform(b,CRS("+proj=longlat")) > data3 <- table_2 > data3$long <- b1$X1 > data3$lat <- b1$X2 > r1 <- raster(nrows=1386, ncols=1649, > xmn=xmn, xmx=xmx, > ymn=ymn, ymx=ymx ) > ras1 <- rasterize(table_2[,1:2], r1, field = table_2[,...

Hi everyone I have a matrix (let's say that it contains the values of elevation) and want to plot its values on a map using a function such as image.plot or filled.contour. The problem is that my grid is not rectilinear, it is bended. Here is an example lon<-matrix(0,20,25) lat<-matrix(0,20,25) elev<-matrix(0,20,25) for (i in 1:20) { for (j in 1:25) {

Hi everybody: Does anyone know how I can change geodesic coordenates into UTM ones? Thanks in advance for your cooperation, Jimmy M. _________________________________________________________________ Bing brings you maps, menus, and reviews organized in one place. a=TEXT_MFESRP_Local_MapsMenu_Resturants_1x1 [[alternative HTML version deleted]]