Displaying 20 results from an estimated 51 matches for "ellps".
Did you mean:
ellis
2016 Jul 20
4
un solo un favor
...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(&qu...
2011 Jan 05
1
Prediction error for Ordinary Kriging
...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))/...
2007 Sep 06
1
write geotiff with projection - RGDAL package
...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...
2010 Oct 12
2
Memory limit problem
...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
? .. .. ..$...
2018 May 24
2
Problem with adding a raster and a brick
...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...
2018 Mar 28
2
Fwd: netCDF to GeoTIFF by layer in r
...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",ove...
2011 Jan 24
0
writeRaster with raster package
...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
obje...
2012 Oct 30
0
lapply and kernelUD (adehabitatHR package): Home Range kernel estimation for a list of individuals
...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) G...
2008 Aug 19
4
converting coordinates from utm to longitude / latitude
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.
2018 Mar 28
0
netCDF to GeoTIFF by layer in r
...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="...
2017 Aug 19
4
My very first loop!! I failed. May I have some start-up aid?
...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)
kern...
2013 Jan 18
2
How to re-project ease( Equal Area Scalable Earth) grid with a ~25 km cylindrical projection to WGS84 0.25 degree?
...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 con...
2012 Jun 05
2
Converiting longitude/latitude to utm
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"))
2010 Apr 01
3
Using GIS data in R
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
2017 May 12
1
spTransform
...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]]
2005 Jun 09
1
Using transform on spproj package.
...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
2006 Feb 01
1
Gauss-Krüger coordinates system
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]]
2012 May 01
0
Plotting shapefiles on existing maps
...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(&q...
2017 Aug 19
0
My very first loop!! I failed. May I have some start-up aid?
...( 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 = "...
2007 Jan 17
1
sp: proj4string has no impact
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
plo...