similar to: while( ) to get a vector until the sum of elements are <= X ?

I'm trying to create a graph where different cells of a grid (a shapefile) will be painted with a color share scale, where the most easy way is to use gray(). Can I somehow get a vector (gradient) of colors, a vector of colors with other methods but gray()? I'm doing this until now quad_N_sp <-

Hi, I'm not being able to capture a position of a 'string' in a character string. In this example: 'There are 20 species in this grid' I would like to capture the string (number) after 'are' and before 'species'. Consider they do not change. I wouldn't like to use substr because stop position may change. Thanks Paulo

I have a xml exported by Manifold GIS but I'm not being able to import it into R using XLM package. The file have this structure: <?xml version="1.0" encoding="UTF-8" ?> - <layout> <name>Layout 2</name> <pagesByX>1</pagesByX> <pagesByY>1</pagesByY> - <elements> <legend

How can I print a data.frame to a PDF with pdf()...dev.off() ____________ Paulo E. Cardoso [[alternative HTML version deleted]]

I'd like to know how to call a program to run or open a specific file. something like this: system('"C:\\Program Files (x86)\\IrfanView\\i_view32.exe"','-A:\\ teste.jpg') is not working. any help will be appreciated ____________ Paulo E. Cardoso [[alternative HTML version deleted]]

Hi, I'm trying to customize a window with 2 graphs. I'm able to do the first one with something like this general example par(mfrow=c(1,2),cex.axis=0.85,cex.lab=0.80,mai=c(1.3,1,0.5,0),las=3) bplot<-barplot(bar.values,names.arg=cf.names,width=0.5,ylab="% Area held") abline(h=0.3,lty=3,col="red") abline(h=0.1,lty=3,col="blue")

Can we by any change perform in R something similar geoprocessing or, more advanced, Spatial SQL over vectorial data? Paulo [[alternative HTML version deleted]]

Hi Team, Need your guidance in building SPDF objects from Latitude, Longitude Information available. I am using package "plotGoogleMaps" which is really awesome but it requires SPDF objects to build the maps. I have a data frame which have Latitude and Longitude information and wanted to convert it in SPDF for making maps. Please help me to find suitable package or function to convert

Dear All, Please see the short script at the end of the email, which I assembled looking for bits and pieces on the web. It essentially does what I need: it plots several countries as a color-coded map. I just would like to fine-tune a bit the final image, in particular (1) Select my own colors for "painting" the countries (i.e. associate manually a color to every level) (2) Be

Hi, The following script (which I did not develop) is used to calculate and plot a skewed normal curve. The script currently requires the user to input six parameters, rather than reading these directly from a file. I've been spinning wheels here, trying to figure out how to modify the script to automate it. I have four data sets, each in excess of 300 records that I need to process. My

Hi, I would like to do coloring of map regions based on the region values "weight". The approach I am taking is first to break regions into equal intervals, classIntervals(spdf$weight,4)$brks #4 intervals in this case and coloring all regions within the interval with the same color col = brewer.pal(4,"RdYlGn")) The max "weight" is as well the boundary of the

I need to get the exact sentence, with the quotes: OPTIONS COORDSYS("Surface1" AS COMPONENT); but this: options(useFancyQuotes = F) paste("OPTIONS COORDSYS(", dQuote("Surface 1"), "AS COMPONENT);" ) give me this: [1] "OPTIONS COORDSYS( \"Surface 1\" AS COMPONENT);" And it's not readable as SpatialSQL code by

Estimados amigos Estoy dibujando las funciones acf y pacf de una variable de una serie "zoo": > ls.str(pat="T0.5") T0.5 : 'zoo' series from 2017-11-08 23:00:00 to 2017-11-15 06:59:00   Data: num [1:9120, 1:3] 55 49.8 51 50.1 36.5 ...   Index:  POSIXct[1:9120], format: "2017-11-08 23:00:00" "2017-11-08 23:01:00" "2017-11-08

Dear List, I am analyzing the home range area of fish and seem to have lost the individuals ID names during my manipulations, and can't find out how to rename them. I calculated the MCP of the fish using mcp() in Adehabitat. MCP's were converted to spatial points data frame and exported to qGIS for manipulations. At this point the ID names were lost. I brought the manipulated

Thank you Sarah and Mike for your explanations. My final objective is to produce maps (png image or any kind of extension I can import in LaTeX) where rainfall data are interpolated, using the Inverse Distance method or Kriging. My input file (pointfile.csv in the reported example) reports the station code, lat and long of the meteorological station and the rainfall value (which might be the

Hi, I am working with a SpatialPolygonsDataFrame of many islands. There are a lot of polygons (islands) composing my SpatialPolygonsDataFrame. I want to extract the elevation of each island. I need to separate the different polygons (like dissolve function in arcgis), to have the elevation of each island. Do you have any idea how can I do that ? I already read a lot of forum, and read the

Dear Roger, I am trying to use the write.polylistShape() function of maptools for the first time and realize that it handles list of polygons of class 'polylist'. However, it seems that no as.polylist() function exist in the package. The question behind that is: in your opinion, which would be the best way to convert a list of matrix of polygon nodes coordinates into an object of

Hola, # Hacemos el KED. Ver función "krige()": KED.rad <- krige( formula=pluvPcp~layer, # covariable -> radar locations=lluvia.rad.pluv.spdf, newdata=radarGrid, # podría ser cualquier objeto Spatial model=v.fit, # modelo de semivariograma. maxdist=Inf

Hello, I'm interested on implementing the MIPS backend as a student in a compiler course on my university, starting in 2 weeks from now!! I've seen on the mailing list history no one has done it so far (although a lot said they would)! am i right? Thanks! -- Bruno Cardoso Lopes http://www.brunocardoso.org "The knack of flying is learning how to throw yourself at the ground and

Hi Manuel, Thanks for the help. Still persists for me too. Instead of waiting indefinitely, now I get this error: Unhandled Exception ("AphrontDeadlockQueryException") #1205: Lock wait timeout exceeded; try restarting transaction On Thu, Dec 11, 2014 at 11:26 AM, suyog sarda <sardask01 at gmail.com> wrote: > The problem still persist :( > > On 12/11/14, Manuel Klimek