search for: wgs84

...https://echange-fichiers.inra.fr/get?k=f9DDllPKdUKs5ZNQwfq from the metadata ,the projection is cylindrical and the resolution is 25 km(it is based on authalic sphere based on International 1924 ellipsoid).As I want to compare with other data, I have to make them identical. - my other data are in WGS84 with 0.25*0.25 degree resolution and extent -180, 180, -90, 90. So I want to re-project the data I have here from EASE grid, cylindrical,25 km ,1383 pixel 586 lines to WGS84,0.25*0.25 degree,1440 pixel 720 lines: f=open.ncdf("C:\\Users\\aalyaari\\Desktop\\SM_RE01_MIR_CLF31D...

...nrow=1) RM179 <- matrix(c(-81.890929, 33.317914), nrow=1) RM148 <- matrix(c(-81.7547337, 33.1514072), nrow=1) RM119 <- matrix(c(-81.501919, 32.94038), nrow=1) RM61 <- matrix(c(-81.262388, 32.524739), nrow=1) RM215.sp <- SpatialPoints(RM215, proj4string=CRS("+proj=longlat +datum=WGS84")) d060101 <- as.POSIXct("2006-01-01", tz="EST") study_seq <- seq(from=d060101, length.out=761, by="days") up.215 <- sunriset(RM215.sp, study_seq, direction="sunrise", POSIXct.out=TRUE) down.215 <- sunriset(RM215.sp, study_seq, direction=&quo...

...ults. 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=2), proj4string = CRS("+proj=longlat +datum=WGS84")) sp1Transformed <- spTransform(sp1, CRS("+proj=merc +zone=32u +datum=WGS84")) coordinates(sp1Transformed) coords.x1 coords.x2 [1,] 998395.0133 6319888.068 The result is an obvious nonsense. Well, after some deliberation I swapped the original coordinate columns: x...

...egion 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 into UTM data. To do so, I use package PBSmapping: attr(my.dataframe, "projection") <-"LL" attr(my.dataframe, "zone") <- 30 my.dataframe<-as.PolySet(my.dataframe) my.dataframeUL<-convUL(my.da...

I am interested in making a random sample from a uniform distribution of points over the surface of the earth, using the WGS84 ellipsoid as a model for the earth. I know how to do this for a sphere, but would like to do better. I can supply random numbers, want latitude longitude pairs out. Can anyone point me at a solution? Thanks very much. -- Russell Senior ``I have nine fingers; you have ten.'' s...

...are 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) reso...

...: RasterLayer 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&qu...

...already loaded 2 different line shapefiles with projection then i will see if i can replicate this problem. So, first file: buff class : SpatialLinesDataFrame features : 1 extent : 560525.5, 561302.4, 4698908, 4701110 (xmin, xmax, ymin, ymax) crs : +proj=utm +zone=16 +datum=WGS84 +units=m +no_defs variables : 2 names : ID, buff_dist value : 2, 135 second file: pln class : SpatialLinesDataFrame features : 1 extent : 560615.6, 560705.4, 4698905, 4699180 (xmin, xmax, ymin, ymax) crs : +proj=utm +zone=16 +datum=WGS84 +units=m +no_de...

...names, grid refs etc allStations <- catalogue() allStations This works OK But when I try to run the example codes there are various problems The code is lifted directly from the document btw it?s not mine This bit to convert gives an error message: # Option a: from OS grid reference to WGS84 > osg_parse(gridRef = "NC581062", CoordSystem = "WGS84") # Option b: from OS grid reference to BNG osg_parse(gridRef = "NC581062", CoordSystem = "BNG") ## osg_parse seems to be built into the rnrfa package so I don?t know why it doesn?t work? The...

...oinformation: ---- 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 object in...

...similar error, so I'm stuck with it... Any thoughts on what I could be doing wrong will be very appreciated! See below the code: FOR ONE ANIMAL 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 l...

..."*.hdf",all.files=T,recursive=T,full.names=T) input outname<-paste(substr(input[1],40,45),".tif",sep='') outname reprojectHDF(hdfName=input,filename=outname,MRTpath="~/MRT/bin", resample_type="NEAREST_NEIGHBOR", proj_type="GEO",datum='WGS84') #### Pero me resulta el siguiente mensaje de error: #### > reprojectHDF(hdfName=input,filename=outname,MRTpath="~/MRT/bin", + resample_type="NEAREST_NEIGHBOR", + proj_type="GEO",datum='WGS84') ***********************************************************...

...: 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,filenam...

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"))

...n 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 +ellps=GRS80 +units=us-ft +no_defs"))) : error in ev...

...> xy.sp SpatialPoints: loni lati [1,] -8.85 38.16 [2,] -9.19 37.99 [3,] -9.11 37.97 [4,] -9.06 38.15 [5,] -9.03 37.87 [6,] -9.14 37.81 [7,] -9.09 37.70 [8,] -8.95 37.45 [9,] -9.17 37.37 [10,] -9.00 37.34 Coordinate Reference System (CRS) arguments: +proj=merc +datum=WGS84 > transform(xy.sp, CRS("+proj=utm +zone=29")) SpatialPoints: loni lati [1,] 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,] 1...

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 or so large that all they really do for us is puff us up and make us feel like god...

Hi R-user, I was trying to make a raster map with WGS84 projection in R, but I could not make it. I found one data set in Google that data is almost the same format as of mine. I wanted to make a raster map of temperature with 1 degree spatial resolution for the global scale. I could make it in GIS software but I do have many variables (to be many rast...

...tion, and origin are 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)...

...ding 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, nc...