similar to: Using RGDAL to "copy" header information...

Dear list members, I'm getting the following error from R CMD check under Windows only with R 2.10.0 devel: ------------- snip ----------- * installing *source* package 'Rcmdr' ... ** libs making DLL ... gcc -I"c:/R/R-2.10.0dev/include" -O3 -Wall -std=gnu99 -c ismdi.c -o ismdi.o gcc -shared -s -o Rcmdr.dll tmp.def ismdi.o -lRgraphapp -Lc:/R/R-2.10.0dev/bin -lR

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]]

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

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

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'

This hacked up version of figlet contributes a supprisingly large proportion of the Coverity issues found under xen/ (and therefore attributed against Xen) Figlet can be found in all distros, so make use of it. We keep xen.flf (being the Xen figlet font) and replace the hacked up octal transform with a short python script. The Xen Makefile has been tweaked in such a way that it still prints the

Don't know if this has been addressed, but why isn't there a file_include style directive for extensions.conf? I find that my extensions.conf grows a lot, and it would be a lot nicer to have a tree of files rather than one big file to try and navigate. Also, I've got a couple different 'systems' running concurrently on one asterisk box (ie, completely different groups of

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

Hi ALL, Can you please help 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",

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

Hello friends of R, My name is Toni, i'm 25 and I'm working on the Meteorological Investigation team from Balearic Islands. I had contact to you because I have a problem: I done a file for every day since 1912 about precipitation. That file has the following structure: > str(Ast) Loading required package: sp Formal class 'SpatialPixelsDataFrame' [package "sp"]

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

Hi All, I am a PhD student in forestry science and I am working with LiDAR data set (huge data set). I am a brand-new in R and geostatistic (SORRY, my background it?s in forestry) but I wish improve my skill for improve myself. I wish to develop a methodology to processing a large data-set of points (typical in LiDAR) but there is a problem with memory. I had created a subsample data-base but

I've been using * CVS code from May of this year and was able to connect to iConnectHere and receive calls with * being behind NAT. Now that I've upgraded to 1.0 RC1, this no longer works. I've tried setting nat=yes in places, externip, et al with no success .. even though the code I was running from back then worked without that. Any suggestions? BTW, I've gotten DTMF from

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

caveat: I haven't either read actionscript language spec nor really got through swfdec context, frame and such. This is fully based on blocky.swf code inspection. thanks to flare, i can extract the below beautiful code. You don't need to understand the whole thing but there is one thing to notice: border[XY]start are initialized in on(press) and get used in on(release). Currently swfdec

Hello dear 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

In the process of investigating a Cisco ATA 186 that was locked by Vonage, I found that you can still unlock the device yourself. But there's a catch. The device's design has a great plus: a DIP32 *socketed* SST28SF040A flash chip. I found an 8 digit unlock code at 0x03FA71-0x03FA78. I do not know if that is a standard location. If you have the equipment, you're in luck. But

Hi, I have a matrix of data which i can vizualize as an image - for example. I would like to save this image as a geotiff file or at a tiff file with a world file which holds the projection of my data (ultimately the data represent a map of some sort). I know i can save the data as an ESRI grid, but i am not interested in that. I wonder if anybody knows about any code which will help me do