Displaying 20 results from an estimated 22 matches for "towgs84".
2010 Oct 12
2
Memory limit problem
...te 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 [1:2] &qu...
2018 May 24
2
Problem with adding a raster and a brick
...cording 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 : 0.5, 0.5 (x...
2018 Mar 28
2
Fwd: netCDF to GeoTIFF by layer in r
...nd : 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",overwrite=TRUE,byl...
2011 Jan 24
0
writeRaster with raster package
...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
object in the work space. The content in the file...
2018 Mar 28
0
netCDF to GeoTIFF by layer in r
...nds)
> #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"...
2018 May 24
0
Problem with adding a raster and a brick
...the 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 : 0.5, 0.5 (x, y)
extent : -70, 30, 50,...
2018 Jun 01
0
Problem with adding a raster and a brick
...d 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)
&g...
2012 Oct 18
1
spTransform longlat to utm
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
2011 Oct 18
1
problem with project command in rgdal
...9;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 commands to...
2013 Mar 18
4
Why stacking rasters return NAs?
...ing the error.
saf
class : RasterStack
dimensions : 720, 1440, 1036800, 601 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : -180, 180, -90, 90 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +ellps=WGS84 +towgs84=0,0,0,0,0,0,0
+no_defs
names : Vol_025_H//00_1_wgs84, Vol_025_H//00_1_wgs84,
Vol_025_H//00_1_wgs84, Vol_025_H//00_1_wgs84, Vol_025_H//00_1_wgs84,
Vol_025_H//00_1_wgs84, , ...
min values : NA, NA,
NA,...
2007 Jul 06
0
import DTM with readRAST6()
...ct of class SpatialGridDataFrame
Coordinates:
min max
coords.x1 1611600 1729700
coords.x2 5058500 5158000
Is projected: TRUE
proj4string :
[+proj=tmerc +lat_0=0.0000000000 +lon_0=9.0000000000 +k_0=0.9996000000
+x_0=1500000.0000000000 +y_0=0.0000000000 +a=6378388 +rf=297 +no_defs
+towgs84=-225.000,-65.000,9.000]
Number of points: 2
Grid attributes:
cellcentre.offset cellsize cells.dim
1 1611650 100 1181
2 5058550 100 995
Data attributes:
DTM_trentino_100m
Min. :1.192e-319
1st Qu.:1.065e-314
Median :1.066e-314
Mean :1.285e-312
3rd Qu.:2...
2012 Feb 23
0
resistanceDistance representation
...Y_transitionclass : TransitionLayer
dimensions : 2181, 1648, 3594288 (nrow, ncol, ncell)
resolution : 100.0049, 99.98945 (x, y)
extent : 1460708, 1625516, 4947383, 5165460 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=tmerc +lat_0=0 +lon_0=9 +k=0.9996 +x_0=1500000
+y_0=0 +ellps=intl +towgs84=-225,-65,9,0,0,0,0 +units=m +no_defs
values : conductance
matrix class: dsCMatrix
I applied two different geocorrections to this same transition:
> MY_correction<-geoCorrection(MY_transition,"c",F,scl=T)
> MY_Rcorrection<-geoCorrection(MY_transition,"r",F,scl...
2012 May 15
0
Problem with legend and RasterLayer
...my problem:
####### Data
> img
class : RasterLayer
dimensions : 5, 5, 25 (nrow, ncol, ncell)
resolution : 0.002245, 0.002245 (x, y)
extent : -48.33875, -48.32753, -20.35756, -20.34634 (xmin, xmax,
ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84
+towgs84=0,0,0
values : in memory
min value : 2116
max value : 36070
##### plot Rasterlayer
> plot(img, col=rev(terrain.colors(25)), axes = TRUE, xlab="x", ylab="y",
useRaster=TRUE)
Erro em as.double(y) :
cannot coerce type 'S4' to vector of type 'double...
2013 Mar 19
1
How to get the rolling standard deviation in rasters?
...> saf
class : RasterStack
dimensions : 720, 1440, 1036800, 12 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : -180, 180, -90, 90 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84
+towgs84=0,0,0
names : SM_RE01_M//7_amenlast, SM_RE01_M//7_amenlast,
SM_RE01_M//7_amenlast, SM_RE01_M//7_amenlast, SM_RE01_M//7_amenlast,
SM_RE01_M//7_amenlast, SM_RE01_M//7_amenlast, SM_RE01_M//7_amenlast,
SM_RE01_M//7_amenlast, SM_RE01_M//7_amenlast, SM_RE01_M//7_amenlast,
--
View t...
2012 Oct 30
0
lapply and kernelUD (adehabitatHR package): Home Range kernel estimation for a list of individuals
...ordinates Name
0 (-683614, 4459280) Gael
1 (-769563, 4516660) Gael
2 (-721607, 4431310) Gael
3 (-683613, 4459290) Gael
4 (-765266, 4502750) Gael
5 (-683602, 4459280) Gael
Coordinate Reference System (CRS) arguments: +proj=merc +lon_0=0 +k=1
+x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs
+towgs84=0,0,0
># Href Fixed Kernel Density Estimator
>K_Gael_estUDm <- kernelUD(Gael_WGS84_WorldM[-c(2:25)], h="href", grid=500,
kern="bivnorm")
FOR THE LIST OF ANIMALS. EVERYTHING SEEMS TO WORK FINE UNTIL I CALL THE
lapply FUNCTION, THEN IT GIVES AN ERROR:
> # Create a...
2006 Oct 18
1
conversion of LL coordenates to UTM problems (ED50-WGS84 format)
Hi R-Users,
I have plotted a region whose polygon coordinates are given in shp format ED50 UTM (zone=30) ) using "readShapePoly" in library(maptools).
Now I need to plot a set of points in that region (my.dataframe, with X and Y geographic coordinates), which have been read using GPS in Longitud-Latitud form (using WGS84 system), so I first need to convert these Longitud-Latitud data
2011 Nov 04
2
How to write a shapefile with projection
...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?
?Mean?? :610300??
Mean? ?:3375307??...
2012 Jun 06
0
randomForest Species Distribution Modelling
...ue).
> predictors
class : RasterBrick
dimensions : 1745, 3909, 6821205, 10 (nrow, ncol, ncell, nlayers)
resolution : 10, 10 (x, y)
extent : 417085.1, 456175.1, 5556329, 5573779 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=utm +zone=30 +ellps=WGS84 +datum=WGS84 +units=m +no_defs
+towgs84=0,0,0
values : in memory
min values : -1 -59 -10 -487 -26 -5 -1199 -14 -16 0
max values : 358 0 19 863 32 11 2551 16 14 34
Step 5
create the 'model'/rf.
model<-factor(G1)~ Slope+Aspect+Curvature+Rugosity+ Plan_Curv+ Prof_Curv+
BS_BPI+ BS_...
2018 Nov 03
1
odd NOTE in R CMD check of data; Rcpp module related
...kage) :
unable to find required package ?raster?
> readRDS("x.rds")
class : RasterLayer
dimensions : 10, 10, 100 (nrow, ncol, ncell)
resolution : 36, 18 (x, y)
extent : -180, 180, -90, 90 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
data source : in memory
names : layer
values : 1, 100 (min, max)
Any idea about what I can do to get rid of this?
Thanks,
Robert Hijmans
[[alternative HTML version deleted]]
2014 Jun 10
2
Como controlar la altura de "colorkey" en levelplot de RatserVis
...: RasterLayer
dimensions : 780, 1296, 1010880 (nrow, ncol, ncell)
resolution : 30, 30 (x, y)
extent : 530156.6, 569036.6, 1092167, 1115567 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=tmerc +lat_0=0 +lon_0=-84 +k=0.9999 +x_0=500000 +y_0=0
+datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0
data source :
C:\Users\Manuel\Documents\Documentos_personales\Proyectos_en_R\Panthera\cobertura.tif
names : cobertura
values : 1, 11 (min, max)
attributes :
ID cobertura
from: 1 Bosque
to : 11 Vegetación anegada
levelplot(cobertura, margin = F, c...