similar to: Reading in .aux (ESRI raster files) into R

I started to use the rgdal package in order to work with DEM files. Since these files are based on non-rectangular quads, there are lots of NA-values in the corresponding datasets. Does anybody know how to remove the NA-values from an object of class GDALDataset in order to get a rectangular dataset? Also, the function saveDataset doesn't seem to work for these kind of data (or at least for

Hi, I'm trying to calculate a sum of differences between two ascii grids. I have the file names of the two grid files in a data.frame along with other data. I'm then trying to add a new column to the data.frame with the result of that calculation. eg. library(maptools) difcount <- function(surfA, surfB) { A <- readAsciiGrid(surfA, colname="a") B <-

I'm not sure I'm barking up the right tree here, but would I need to make use of groups to plot two separate datasets within ONE panel in xyplot? The desired end result is a single xy plot of two separate (but similar in values and ranges). Full code follows, xyplot code at bottom #########Determine Frequencies ##########coastal_slope #needs the maptools package to read ESRI grid

Hi all, I have a set of 54 files that I need to convert from ASCII grid format to .shp files to .bnd files for BayesX. I have the following R code to operate on those files: library(maptools) library(Grid2Polygons) library(BayesX) library(BayesXsrc) library(R2BayesX) readfunct <- function(x) { u <- readAsciiGrid(x) } modfilesmore <- paste0("MaxFloodDepth_", 1:54,

>From the following: basin.map <- readAsciiGrid("c:/temp/area.asc", colname="area") I have a SpatialGridDataFrame which has the x and y cordinate of a cell, and the drainage area of that cell. There are many cells with a low drainage area (in my case, 33000 with an area of 37.16) and one cell with the highest drainage area (again, in my case, a drainage area of of

Hi All, I have a greyscale image that I am reading in through RGDAL and placing in a pixmap object. As an example use the logo.jpg file that comes with the RGDAL package: #Read the file logo <- system.file("pictures/logo.jpg", package="rgdal")[1] x <- new("GDALReadOnlyDataset", logo) #Create the pixmap object xGrey <- getPixmapGDAL(x) Now I would

Dear All, I am a research student in environment. I have only little programming knowledge. I am currently doing the last project about rainfall impact on ground water quality in an area. It happens that I have to use R to read rainfall data (3 dimension) from ASC file (*.asc), and then write them into one NCDF file (*.nc). I have been working very hard on study R, but I

Does Mapstraction support ESRI (arcwebservices.com)? Thanks, - Michael

Hello, I have an ESRI ArGis geo database directory which I would like to convert to geojson or a shape file (or anything else that I can read into R). Unfortunately that does not work out of the box with rgdal, because it does not come with the fileGDB or openfileGDB driver. I could successfully install gdal and the fileGDB driver/extension, but it seems that i can use gdal only to convert the

We just upgraded our samba server from 3.0.14 to 3.0.21a and have discovered a subtle but very painful problem. Users of ESRI ArcMap (on fully patched Windows XP clients) are unable to read complex layered file data. Arcmap works fine with 3.0.14, it works fine with the data on a windows share, and it works fine when I set up a server to run the version of samba shipping with Solaris 9 (3.0.2! ).

Hi, I am working on a library of map functions for R and I have started with code to read shapefiles (the format is not great but it is ubiquitous). I am not a windows person, but where I am working I don't have access to any UNIX machines easily. So here's the problem I wrote C code to read the information out of the shapefile header and I compiled it using the script given in Ripley and

Well Dr. Ripley helped me solve the first problem, next question I need to do some debugging now that the function loads properly into R. I checked the S documentation and to use the C I/O functions you need to include newredef.h, this file doesn't exist in /src/include. How does this work in R? under S3.4 on a dec station i could use fprintf(stdout,..); fflush(stdout); thanks, Nicholas PS

Hello- Does anyone know if R can read in ESRI file geodatabase feature classes, directly from the file geodatabase? I currently have to export the data out of Arc, to a csv file, then read the csv file. I'd like to bypass this step and read directly from the file geodatabase. Thanks for your help. -al

I just uploaded two packages to CRAN. shapefiles_0.1.tar.gz - functions to read and write ESRI shapefiles (including dbfs) emme2_0.1.tar.gz - functions to read binary data from an EMME/2 databank data (EMME/2 is a transportation modeling program) Please let me know if you find any bugs or have some suggestions. Thanks. Regards, Benjamin Stabler Transportation Planning Analysis Unit Oregon

Did not hear from anyone on this, so maybe we are the only ESRI users who got bit, but this turned out to be an issue with oplocks. Disabling oplocks on the server eliminated the problem. On 1/27/06, Elizabeth Schwartz <betsy.schwartz@gmail.com> wrote: > > We just upgraded our samba server from 3.0.14 to 3.0.21a and have > discovered a subtle but very painful problem. > Users of

thanks to all the help so far I have code now that works on windows (albeit very clunky) that can read the pertinent info from the shape header and return a list with all the necessary parameters to read info out of the shape file. Now come the big design questions 1) I have considered building a map object, basically a list that would hold the info from the shape file and a data frame as part of

Hi all, who knows how to get an ESRI Shapefile for the NUTS-2 Regions of the enlarged European Union? Particularly I want to draw maps of Germany, Poland, Czech Republik, Hungary and Austria. I've found Shapefiles for the US, Russia and other countries elsewhere in the web, but for Europe it seems really difficult. With many thanks for any hint Albrecht

Sorry, I did not make my question clear. Since I have a point theme with many points, some of them may clump together. the problems here are: 1. how to find clumps in a point theme? 2. the convex-hull extension I found only deal with all the points in a theme at each time? how to make each convex hull around each point clump automatically? Thanks. Xiaohua On 5/12/06, Bob Booth <bbooth

Hello, I'm attempting to evaluate the accuracy of the probability predictions for my model. As previously discussed here, the AUC is not a good measure as I'm not concerned with classification accuracy but probability accurcy. It was suggested to me that the loess function would be a good measure to look at. I can see some libraries (Design) will plot the loess function as a curve

Dear R user, I am running R under Windows 2000. I am looking for a routine for importing - shape files (ESRI) into R - dbase files (FOXPRO) into R and I am looking for time-space models for description and prediction of Bernoulli-, Binomial- and Poissonvaraibles. Thank's a lot for a reply. Sincerely yours, Ekkehardt Altpeter Swiss Federal Office of Public