search for: cellsize

...i suspected global scope), not cx->frame->var_object. It worked. So I think swfdec_as_context_eval_set_property needs to be rechecked. blocky.swf is flash 6, according to swfmill. button 155 { on (press) { if (_root.playgame) { borderXstart = Math.floor((_xmouse - 13) / CellSize) * CellSize + 13; borderYstart = Math.floor((_ymouse - 13) / CellSize) * CellSize + 13; switch (random(3)) { case 0: klik.start(0, 1); break; case 1: klik1.start(0, 1); break; case 2: klik2.sta...

...ear all, A simple understanding question but I cannot find the answer anywhere. I use coordiates (longitude, latitude) in my study. I created a SpatialPoints class from those data and then I created a grid using the 'ascgen' command in the package 'adehabitatMR' like: ascgen(xy,cellsize=40). It worked out fine and my analysis is doing great. But I don't know what unit the 'cellsize' is. Is it the area (squaremeter), or the length of the border of each rastercell (meter)? Or is it some dimension that I cannot grasp from a geografical sense? I can't find the info...

Hello, Is there anyone who could give me an example of how to import an ascii grid (i.e. ArcGIS exported raster) into R. I want to use it with gstat but don't know the appropriate import routine. Thanks very much for your help. Regards, femke Femke Reitsma Graduate Student (ABD) Geography Department 2181 LeFrak Hall University of Maryland College Park, MD 20742 Phone: 301-405-4121

Hello, Could someone please tell if there is a way to append header info to the top of an exported dataframe (exported with write.table). I want to append the following, which are the defininitions for an asciigrid: ncols 532 nrows 999 xllcorner 510465 yllcorner 4766375 cellsize 30 NODATA_value -9999 Thanks very much, femke [[alternative HTML version deleted]]

R-geographers... I'm trying to solve a problem to implement a line-by-line tiled processing using RGDAL (read 1 line of an image, process the one line, write one line of the image to a binary file). Everything except for the final step I'm able to do using a combination of RGDAL and r-base commands. Below is the basic structure, the input file "elev" can be any image file

Dear All, I'v tried to read in data in numerous ways including:- Spain10km<-data.frame(readAsciiGrid("F://RMap//sp10kpointid1.aux")) Error in readAsciiGrid("F://RMap//sp10kpointid1.aux") : object 'cellsize' not found In addition: Warning message: In readLines(t, n = 6) : incomplete final line found on 'F://RMap//sp10kpointid1.aux' Spain10km<-data.frame(readAsciiGrid("F:/GIS.LandcoverEuropeForRisk/Sept10kmmaps/Sp10KPointID.aux")) Error in readAsciiGrid("F:/GIS.Landcove...

...ked but I wanted finer control of the bandwidth. Therefore, I was advised to switch to the package splancs. I've been trying to copy the example provided in the help file but cannot produce the same results for my data. I've pasted my code below and I am not sure how to define the grd1 (cellsize and cells.dim). Here is my code: coords <- S at polygons[[1]]@Polygons[[1]]@coords coord <- as.points(coords) cases1.xy <- as.points(mergedis$jit.x.x, mergedis$pre.fsa.shp.Y2) grd1 <- GridTopology(cellcentre.offset=c(50281627.2, 580082), cellsize=c(0.2, 0.2), cells.dim=c(75,100))...

Dear Happy R-users & experts, I am in need of advice, While working with spatial data (x & y coordinates of seed locations) I have come accross the problem that I need to convert my point data into a matrix or grid system. I then need to count how often a point falls into a certain position in the matrix or grid. I have searched all day online, asked collegeas but nothing works. Sadly

...al("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 <- GridTopology(cellcentre.offset = cc,cellsize = cs, cells.dim = cd) getClass("SpatialGrid") p4s <- CRS(proj4string(rsa2)) x2_SG <- SpatialGrid(rsa2_grd, proj4string = p4s) x2_SP <- SpatialPoints(cbind(x2$X, x2$Y)) v <- variogram(log1p(tsport_ace) ~ 1, x2, cutoff=100, width=9) te<- fit.variogram(v,vgm(0.0437, "Exp&...

Hello R-World, i am trying to calculate data for a DEM (Digital Elavation Model) which i also want to plot in R. i have the coordinates for the lower left corner which look something like this: x<-42,2 y<-50,5 besides i have the cellsize of the grid, which is: z<-1,1 what i do is to calculate the corrdinates of the cells to the right and top, what i can do by specifying the number of rows and columns. these are j<-1:805 (columns) i<-1:616 (rows) well, to calculate them seperatly is easy: r<-x+(i-1)*z h<-y+(j-1)*z...

...at what the problems are in my R script? Firstly, this is an example of data in asc file: NCOLS 241 NROWS 160 XLLCORNER 140.00012207031 YLLCORNER -39.000000000000 CELLSIZE 0.50000000000000E-01 NODATA_VALUE -99.00000 166.30 160.87 155.23 149.33 143.83 138.52 133.29 128.34 123.76 119.21 115.06 110.90 107.22 103.69 100.40 97.29 94.58 92.15 90.00 87.91 86.2...

...1.90e-05 1.90e-05 1.90e-05 ... ..@ coords.nrs : num(0) ..@ grid :Formal class 'GridTopology' [package "sp"] with 3 slots .. .. ..@ cellcentre.offset: Named num [1:2] 345 4278 .. .. .. ..- attr(*, "names")= chr [1:2] "x" "y" .. .. ..@ cellsize : Named num [1:2] 0.1 0.1 .. .. .. ..- attr(*, "names")= chr [1:2] "x" "y" .. .. ..@ cells.dim : Named int [1:2] 2684 1600 .. .. .. ..- attr(*, "names")= chr [1:2] "x" "y" ..@ grid.index : int [1:499105] 2468192 2468...

Dear exports,I have created a hypsometric curve (area-elevation curve) for my watershed by using simple command hypsometric(X,main="Hypsometric Curve", xlab="Relative Area above Elevation, (a/A)", ylab="Relative Elevation, (h/H)", col="blue")It plots the hypsometric curve in "RGraphics window", My question is how can I export values which is used

...og(Z)~X+Y, subground, cutoff=1800, width=80, alpha=c(135, 45)) v.fit = fit.variogram(v, vgm(psill = 1, model="Gau", range=1800, nugget= 0, anis=c(135, 0.5))) plot(v, v.fit, plot.nu=F, pch="+") # create the new grid new.grid <- spsample(subground, type="regular", cellsize=c(1,1)) gridded(new.grid) <- TRUE fullgrid(new.grid) <- TRUE new.grid a grid #using Universal Kriging subground.uk = krige(log(Z)~X+Y, subground, new.grid, v.fit, nmax=40) #ERROR

Is there any R software that create an image from Yahoo maps together with points of known UTM coordinates (or lat/long marked? Note that my region of interest is not covered in sufficient detail by Google maps. It actually does not have to be Yahoo maps as long as it has sufficient coverage of my region but that's the one I have found with coverage of my region. The scale I am interested in

...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.722e-312 Max. :4.774e-312 As you can see, the import has been fully done, yet the attri...

...already created with spplot(). I would like shaded gridboxes, the boxes should be filled with colour. Can anyone help me with this problem? Thank you in advance! This is part of my code: minlon<-4.5 maxlon<-12 minlat<-46.5 maxlat<-52.5 gt=GridTopology(cellcentre.offset=c(5.75,47.5),cellsize=c(2.5,2),cells.dim=c(3,3)) grd=SpatialGrid(gt,proj4string = CRS(as.character(NA))) gridparameters(grd) gd<-as.data.frame(grd,"SpatialGrid") h<-season_djf h<-as.data.frame(h) djf.att=SpatialPointsDataFrame(gd,h) gridded(djf.att)=TRUE spplot(djf.att, col.regions = colorRampPalette...

...A_PER_HA_ is not present in the predicted data ################################################################ The first six lines of the output file contain the header information for an ESRI ASCII grid and look like this ncols 1564 nrows 1062 xllcorner -77876.53335875 yllcorner -3304144.3940076 cellsize 4.000718 NODATA_value -9999 Which is exactly the same as the input ascii grids used during the AsciiGridImpute command. I was wondering if anyone on the R mailing list had ever come across this error and if so how did you correct it. I am using R version 2.7.1 on Ubuntu 7.10 Many thanks for all...

...roj4R package. SpatialDataFrame extends this class; holds a data frame, and the information where in the data frame the coordinates (2D or 3D) are stored. SpatialDataFrameGrid extends SpatialDataFrame for the case where the data are on a regular 2D/3D grid; it contains the offset of the grid, the cellsize in (x,y,z) and the nr of row/cols/layers. One simple way of creating these classes (idea taken from Barry) is: > data(meuse.grid) # data frame with gridded data > coordinates(meuse.grid) = c("x", "y") # promote to SpatialDataFrame > gridded(meuse.grid) = TRUE # prom...

Hi, I am stuck in a looping problem. It might be an easy problem for experienced R users but I have been unable to do it. Any kind of help or advice will be great appreciated. I am creating a spatial grid and have a list of spatial points in a folder. I can read the spatial points separately and get a count of the points in the spatial grid in R. The output should be a table showing the counts of