Displaying 7 results from an estimated 7 matches for "grs80".
2016 Jul 20
4
un solo un favor
...o 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("+x_...
2007 Sep 06
1
write geotiff with projection - RGDAL package
...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 dirrec...
2010 Oct 12
2
Memory limit problem
...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
? .. .. ..$ : chr...
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
...ne 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]]
2011 Oct 18
1
problem with project command in rgdal
...s. 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
+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 co...
2011 Nov 04
2
How to write a shapefile with projection
...p <- SpatialPointsDataFrame(try[,1:2], try,
proj4string=CRS(crs))
?
summary(crest.sp)
Object of class SpatialPointsDataFrame
Coordinates:
??????? min?????? ??????? max
x? 610235.1? ?????? 610354.1
y 3374862.4 ?????? 3375751.4
Is projected: TRUE
proj4string :
[+proj=utm +zone=15 +ellps=GRS80 +datum=NAD83 +units=m
+no_defs +towgs84=0,0,0]
Number of points: 890
Data attributes:
????? ???????? ?x????????? y???????????????? z1??????
?Min.?? :? 610235?? Min. : 3374862?? Min.: 6.966?
?1st
Qu.:610269?? 1st Qu.:3375085?? 1st Qu.:7.570?
?Median
:610298?? Median :3375307?? Median :7.901...