search for: ellp

...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 +ellps=GRS80 +units=us-ft +no_defs"))) : error in evaluating the argument 'CRSobj' in selecting a method for function 'spTransform': Error in CRS(paste(&q...

...L, Can you please help me on how to determine the prediction error for ordinary kriging?Below are all the commands i used to generate the OK plot: rsa2 <- readShapeSpatial("residentialsa", CRS("+proj=tmerc +lat_0=39.66666666666666 +lon_0=-8.131906111111112 +k=1 +x_0=0 +y_0=0 +ellps=intl +units=m +no_defs")) x2 <- readShapeSpatial("ptna2", CRS("+proj=tmerc +lat_0=39.66666666666666 +lon_0=-8.131906111111112 +k=1 +x_0=0 +y_0=0 +ellps=intl +units=m +no_defs")) bb <- bbox(rsa2) cs <- c(1, 1) cc <- bb[, 1] + (cs/2) cd <- ceiling(diff(t(bb))...

...here i am doing a mistake when i define the projection in R. My code follows: data.grid <- read.csv(x, header=TRUE) gridded(data.grid) = ~East.m.+North.m. proj4string(data.grid) = CRS("+proj=tmerc +lat_0=0.00000000000 +lon_0=-81.00000000000 +k=0.99960000 +x_0=500000.000000 +y_0=0.0000000 +ellps=GRS80 +units=m") tr <- "e:\\JELA_veg\\test_gtiff\\test.tif" writeGDAL(data.grid["class.pca"], tr) mg3 <- readGDAL(tr) proj4string(mg3) [1] " +proj=utm +zone=17 +ellps=GRS80 +units=m +no_defs" I will really appreciate if anybody can point me in the right...

...cannot allocate vector of size 240 Kb > memory.limit() [1] 4000 > memory.size() [1] 1971.68 > summary(hi.to.utm) Object of class SpatialPoints Coordinates: ??????? min?????? max x? 708745.5? 923406.7 y 2046153.1 2327910.9 Is projected: TRUE proj4string : [+proj=utm +zone=4 +datum=NAD83 +ellps=GRS80 +towgs84=0,0,0] Number of points: 15328 > str(hi.to.utm) Formal class 'SpatialPoints' [package "sp"] with 3 slots ? ..@ coords???? : num [1:15328, 1:2] 708746 710482 712218 713944 715681 ... ? .. ..- attr(*, "dimnames")=List of 2 ? .. .. ..$ : NULL ? .. .. ..$...

...e same, so according to my understand it should work. The objects look like so: > obs.clim class : RasterLayer dimensions : 60, 200, 12000 (nrow, ncol, ncell) resolution : 0.5, 0.5 (x, y) extent : -70, 30, 50, 80 (xmin, xmax, ymin, ymax) coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 data source : /home/mpayne/Documents/Predictability_engine/scratch/Bluefin/HadISST/ obs_climatology.nc names : sst z-value : 1988-02-15 zvar : sst > mdl.anom class : RasterBrick dimensions : 60, 200, 12000, 1 (nrow, ncol, ncell, nlayers) resolution...

...asterLayer band : 1 (of 4 bands) #the bands here are the four soil layers dimensions : 16800, 43200, 725760000 (nrow, ncol, ncell) resolution : 0.00833333, 0.00833333 (x, y) extent : -180, 179.9998, -56, 83.99995 (xmin, xmax, ymin, ymax) coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 data source : C:\FAPSEP_Eucalyptus\FAPSEP\Soil_Weather_data\Agro_IBIS_Eucalipto\VMC21\VMC21.nc names : Volumetric.water.content.at..33.kPa z-value : 4.5 zvar : VMC2 VMC1<-writeRaster(importnetcdf,filename="file.tiff",format="GTiff",ov...

...ing geoinformation: ---- R output > mask class : RasterLayer filename : E:/Daten/FE/HyMAP/Luxembourg_2010/Kehlen_Useldange/mosaik/LUX_LC_noOverlap_mask nrow : 5198 ncol : 2813 ncell : 14621974 min value : 0 max value : 255 projection : +proj=utm +zone=31 +ellps=WGS84 +datum=WGS84 +units=m +no_defs +towgs84=0,0,0 extent : 710938, 722190, 5501612, 5522404 (xmin, xmax, ymin, ymax) resolution : 4, 4 (x, y) --- If I now use writeRaster to write the file in Envi format, the geoinformation is lost for the object written to the file BUT not for the obj...

...NAMED "Gael". IT WORKED PERFECTLY! ># Read a shapefile and convert it into a SpatialPointsDataFrame with its corresponding CRS >Gael_WGS84_WorldM <- readShapePoints("900_Gael_WGS84_WorldM", proj4string=CRS("+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs")) ># Remove all the columns except the name of the animal to use the kernelUD function. My data looks like this: > head(Gael_WGS84_WorldM[-c(2:25)]) coordinates Name 0 (-683614, 4459280) Gael 1 (-769563, 4516660) Gael 2 (-721607, 4431310)...

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.

...1 (of 4 bands) > #the bands here are the four soil layers > dimensions : 16800, 43200, 725760000 (nrow, ncol, ncell) > resolution : 0.00833333, 0.00833333 (x, y) > extent : -180, 179.9998, -56, 83.99995 (xmin, xmax, ymin, ymax) > coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 > data source : > C:\FAPSEP_Eucalyptus\FAPSEP\Soil_Weather_data\Agro_IBIS_Eucalipto\VMC21\VMC21.nc > names : Volumetric.water.content.at..33.kPa > z-value : 4.5 > zvar : VMC2 > > > > VMC1<-writeRaster(importnetcdf,filename=&quo...

...mysample <- myframe[sample(1:nrow(myframe), 3,replace=FALSE),] mysample mysampletime <- as.POSIXct (strptime(as.character(mysample$Timestamp), "%d.%m.%Y %H:%M"), tz="GMT") mysamplexy <- project (cbind (mysample$Longitude, mysample$Latitude), "+proj=utm +zone=32 +ellps=WGS84") colnames(mysamplexy) <- c ("xCord", "yCord") ID <- mysample$ID datltraj <- as.ltraj(mysamplexy, mysampletime, id=ID) Ddat <- BRB.D(datltraj, Tmax=21600, Lmin=36) BRBdat <- BRB(datltraj, D= Ddat,type=c("UD"),Tmax=21600,Lmin=36, hmin=100) ker...

...603037, 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"? -- View this message in co...

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

...nal line of code looks like this: pointsTransformed <- spTransform(pointsutm, CRS("+proj=longlat +datum=WGS84 +zone=36S")) The line before looks like: pointsutm <- SpatialPoints(cbind(pointsmapspatial$Point_Y, pointsmapspatial$Point_X), proj4string = CRS("+proj=utm +zone=36S +ellps=GRS80") Can you see what I did wrong? [[alternative HTML version deleted]]

...1505638 38.62347 [2,] 1505638 38.45141 [3,] 1505638 38.43116 [4,] 1505638 38.61335 [5,] 1505638 38.32995 [6,] 1505638 38.26922 [7,] 1505638 38.15788 [8,] 1505638 37.90485 [9,] 1505638 37.82388 [10,] 1505638 37.79351 Coordinate Reference System (CRS) arguments: +proj=utm +zone=29 +ellps=WGS84 What am I doing wrong ? Thanks EJ

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]]

...prj file accompanying the *.shp), try using the transform() method in the rgdal package to do this. If there is a *.prj file, this should work: library(rgdal) my_lines <- readOGR(<directory>, <shapefile_without_extension>) my_lines1 <- transform(my_lines, CRS("+proj=longlat +ellps=WGS84")) *I believe the above line should be: my_lines1 <- spTransform(my_lines, CRS("+proj=longlat +ellps=WGS84")) plot(my_lines1, add=TRUE) If not, you'll need to work out exactly what the projection is, and set it in PROJ.4 notation by: proj4string(my_lines) <- CRS(&...

...( mysample[ , c( "Longitude" , "Latitude" ) ] ) , "+proj=utm +zone=32 +ellps=WGS84" ) colnames( mysamplexy ) <- c( "xCord", "yCord" ) datltraj <- as.ltraj( mysamplexy, mysample$mysampletime, id=mysample$ID ) Ddat <- BRB.D( datltraj, Tmax=21600, Lmin=36 ) BRBdat <- BRB( datltraj, D = Ddat, type = &quot...

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(-1,1),ylim=c(-1,1)) looks exactly like pl...