search for: no_def

...equeridos. > 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_0=-2000000.0 +y_0=-50...

...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_defs variables :...

...lp 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))/cs) rsa2_grd <- Gri...

...gt; 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 the work space. The content in...

...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) Gael 3 (-683613, 4459290) Gael 4 (-76...

...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 context: http://r.789695.n4.na...

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

...ws 10 > columns 11 > bands 1 > origin.x 150701.4 > origin.y 7744897 > res.x 250 > res.y 250 > ysign -1 > oblique.x 0 > oblique.y 0 > driver GTiff > projection +proj=utm +zone=23 +south +datum=WGS84 +units=m +no_defs > file > /MOD13Q1.A2001049.h13v11.005.2007002215512.250m_16_days_EVI.tif > apparent band summary: > *GDType* Bmin Bmax Bmean Bsd hasNoDataValue NoDataValue > 1 *Int16* -32768 32767 0 0 FALSE 0 > Metadata: > AREA_OR_POINT=Point > TIFFTAG_SOFTW...

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 plot(SpatialPoints(xy, proj4...

...+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 dirrection. Thanks, Monica _________________________________________________________________ [[alternative HTML version deleted]]

...9 -2.1e+09 -2.1e+09 -2.1e+09 -2.1e+09 -2.1e+09 NA NA NA max values : NA 2.1e+09 2.1e+09 2.1e+09 2.1e+09 2.1e+09 2.1e+09 NA NA NA Step 3 I then provided the projection. projection(predictors)<- "+proj=utm +zone=30 +ellps=WGS84 +datum=WGS84 +units=m +no_defs" Step 4 I've tried numerous ways to get rid of /relace the NA values. #na.action<- #predictors<-predictors[na.rm=FALSE] #99999->predictors[predictors==NA, ] #predictors<-predictors[predictors, na.action=na.omit ] #na.exclude->predictors=NA #na.omit(predictors) multiple co...

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

...ect 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 make my raw input (a csv file)...

...uot;1:4") #bylayer=TRUE, did not work. > VMC1 class : RasterLayer 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 +no_defs +ellps=WGS84 +towgs84=0,0,0 data source : C:\Users\Ahmed\Documents\file.tif names : file values : 2, 80 (min, max) Any help! Ahmed Attia, Ph.D. Agronomist & Soil Scientist [[alternative HTML version deleted]]

...id not work. > >> VMC1 > class : RasterLayer > 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 +no_defs +ellps=WGS84 > +towgs84=0,0,0 > data source : C:\Users\Ahmed\Documents\file.tif > names : file > values : 2, 80 (min, max) > > Any help! > > > Ahmed Attia, Ph.D. > Agronomist & Soil Scientist > > [[alternative HTML version deleted]] > &g...

...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, NA,...

...-0.7966338 -2.144316 #read in the shapefile england<-readShapePoly('D:/arcGIS/england boundary/england.shp') class(england) #define the projection proj4string(england)<-CRS('+proj=tmerc +lat_0=49 +lon_0=-2 +k=0.999601271625 +x_0=400000 +y_0=-100000 +ellps=airy +units=m +no_defs') # transform the map into the WGS84 projection (epsg:4326): england.wgs<-spTransform(england, CRS('+init=epsg:4326')) plot(england.wgs) #plot data over the map: quilt.plot(data$Longitude, data$Latitude, data$perc_per_year, add=TRUE) My problem is that I would like to be able to c...

...rid) Object 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...

...set", "dummy.tif") > GDALinfo("dummy.tif") rows 420 columns 660 bands 1 ll.x -55.5 ll.y -14.5 res.x 0.008333334 res.y 0.008333334 oblique.x 0 oblique.y 0 driver GTiff projection +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs file dummy.tif > nrow(a) [1] 420 > ncol(a) [1] 660 > dummy = getRasterData(a) > nrow(dummy) [1] 660 > ncol(dummy) [1] 420 > test = matrix(0, nrow=420, ncol=660) > putRasterData(a,test) Error in putRasterData(a, test) : GDAL Error 5: Access window out of range in Ras...

...ayer 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=T) My coords are 44 points where s...