Re: [R-sig-Geo] RSAGA multi.focal.function

2010-06-04 Thread Alexander Brenning
ow—globalsoilmap.net West Virginia University Division of Plant and Soil Sciences National Soil Survey Center Geospatial Research Unit ph: +1 304 293 9839 [[alternative HTML version deleted]] -- Alexander Brenning brenn...@uwaterloo.ca - T +1-519-888-4567 ext 35783 Department of Geography and Env

Re: [R-sig-Geo] Rpy with ArcMap or Spatial Data

2010-02-17 Thread Alexander Brenning
h.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Alexander Brenning brenn...@uwaterloo.ca - T +1-519-888-4567 ext 35783 Department of Geography and Environmental Management University of Waterloo 200 University Ave. W - Waterloo, ON - Canada N2L 3G1 http://www.fes.uwaterloo.ca/geograph

Re: [R-sig-Geo] RPyGeo Query -- probably easy

2009-12-18 Thread Alexander Brenning
arten Maarten I'm not sure about the single quotes but possibly you could add an extra special character break (the \ ) to indicate to python that the double quotes surrounding your field name are not the end of a string. You could try the following code: rpygeo.geoprocessor("s

Re: [R-sig-Geo] Including SAGA grid as GDAL-supported format

2009-10-16 Thread Alexander Brenning
a[[1]]) spplot(sdat.sp) Of course, it would be much easier to have this in a single line: meuse.grid <- readGDAL("meuse_soil.sgrd") or writeGDAL(meuse.grid["soil"], "meuse_soil.sgrd", "SAGA") PS: A new version of SAGA has just been released few days

Re: [R-sig-Geo] RSAGA, Bayes' Theorem and summary stats

2009-09-30 Thread Alexander Brenning
Hi Steve, I recommend checking the data type of all involved rasters in SAGA GIS (GUI) itself, maybe one of the grids involved is integer so everything gets rounded to zero? After loading a grid into SAGA GIS, select it in the 'Data' tab of the 'Workspace' window, and examine the 'Description'

[R-sig-Geo] RSAGA under Linux (was: Slope and Aspect calculations in R)

2009-07-10 Thread Alexander Brenning
nfo/r-sig-geo -- Thomas E Adams National Weather Service Ohio River Forecast Center 1901 South State Route 134 Wilmington, OH 45177 EMAIL: thomas.ad...@noaa.gov VOICE: 937-383-0528 FAX:937-383-0033 ___ R-sig-Geo mailing list R-sig-Geo@stat.math.eth

Re: [R-sig-Geo] matching points to a raster

2009-06-05 Thread Alexander Brenning
y help. Many thanks, Juliane -- Alexander Brenning brenn...@uwaterloo.ca - T +1-519-888-4567 ext 35783 Department of Geography and Environmental Management University of Waterloo 200 University Ave. W - Waterloo, ON - Canada N2L 3G1 http://www.fes.uwaterloo

Re: [R-sig-Geo] RPyGeo Query -- probably easy

2009-04-13 Thread Alexander Brenning
ok at ?rpygeo.geoprocessor - I am a bit in a hurry now because I am on a trip to Chile and ready for a delicious dinner ;-) So try using c( ) (after having a closer look at the rpygeo.geoprocessor help / command syntax), or use separate rpygeo.geoprocessor calls for the two gp. commands...

Re: [R-sig-Geo] R interface to SEXTANTE?

2009-04-06 Thread Alexander Brenning
t, Nick ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Alexander Brenning brenn...@uwaterloo.ca - T +1-519-888-4567 ext 35783 Department of Geography and Environmental Ma

Re: [R-sig-Geo] raster package and "?focal" FUN

2009-04-06 Thread Alexander Brenning
ction or package for to do this but the R world is very huge... es: mobile window 10 m pixelmax <- value NO pixelmax <- 0 thanks Gianni [[alternative HTML version deleted]] ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat

Re: [R-sig-Geo] RPyGeo Query -- probably easy

2009-04-01 Thread Alexander Brenning
ticking point is hrwk05 <> 0. In Python my query would be double quote, field name, close double quote, comparison, number e.g, "hrwk05" <> 0 but I'm not sure how to put this into RPyGeo If I run like so: rpygeo.geoprocessor("makefeaturelayer('X:\\abc.shp

Re: [R-sig-Geo] extract values to points

2009-03-05 Thread Alexander Brenning
_ [[alternative HTML version deleted]] ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Alexander Brenning brenn...@uwaterloo.ca - T +1-519-888-4567 ext 35783 Department of Geography and Environmental

Re: [R-sig-Geo] RPyGeo: merge geoprocessing problems with "field_mappings" option

2009-03-03 Thread Alexander Brenning
with ArcGIS Geoprocessing tools. Cheers Alex gianni lavaredo wrote: Hi All I am testing RPyGeo (old version) developed to Prof. Alexander Brenning. I wish to merge three points shapefile each with this database: ID, X, Y, V. The Scripting syntax in Python is "Merge_management (inputs,

Re: [R-sig-Geo] testing DEM geoprocessing in RSAGA, SAGA GIS and ARCGIS: different of results in value range

2009-03-02 Thread Alexander Brenning
= 4) rsaga.aspect("DEM","ASPECT") rsaga.sgrd.to.esri("ASPECT", prec= 4) but there are some difference when I convert the ASPECT.sgrd in ASCII file. ex ASPECT.sgrd: min value= 0.020869 max value = 360.1 (the same Value range I have in SAGA GIS or in ARCGIS) ASPE

Re: [R-sig-Geo] HELP with RPyGeo and RSAGA slope function

2009-03-02 Thread Alexander Brenning
ut.ascii.file)) %in% c(".asc", ".txt"))) stop("'out.ascii.file' must have extension '.asc' or '.txt'.\n") args = list(in.raster, out.ascii.file) rpygeo.geoprocessor(fun = "RasterToASCII_conversion", arg

Re: [R-sig-Geo] RSAGA .asc morphology loop question

2009-03-01 Thread Alexander Brenning
y processes one file and stops. My directory now includes the following files: [1] "SLOPEtmp1_220E-15N2005-03-02" "tmp1_220E-15N2005-03-02" "tmp10_220E-15N2005-05-13" "tmp100_220E-15N2007-05-13" Thanks again for suggestions on how to process the whole d

Re: [R-sig-Geo] RSAGA .asc morphology loop question

2009-03-01 Thread Alexander Brenning
g in the grid files and would appreciate advice. Best regards, Tim Sippel University of Auckland [[alternative HTML version deleted]] ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Alexa

Re: [R-sig-Geo] difficult to subtract two column of my shapefile

2009-02-24 Thread Alexander Brenning
dOGR(".", "theme1") my_SPDF$ab <- my_SPDF$a - my_SPDF$b writeOGR(my_SPDF, ".", "mynewfile", driver="ESRI Shapefile") ___ R

Re: [R-sig-Geo] R+SAGA: Table calculator for shapes

2009-02-23 Thread Alexander Brenning
Shapes: canopy_height.shp Formula: a-b Field Name: a-b [[alternative HTML version deleted]] ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Alexander Brenning brenn...@uwaterloo.ca - T +1-519-888-4567 ext 35783 Dep

Re: [R-sig-Geo] problem with merge grids in R+SAGA

2008-11-13 Thread Alexander Brenning
R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Alexander Brenning [EMAIL PROTECTED] - T +1-519-888-4567 ext 35783 Department of Geography and Environmental Management University of Waterloo 200 University Ave. W - Waterloo, ON - Cana

Re: [R-sig-Geo] Extracting data/value

2008-10-31 Thread Alexander Brenning
___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Alexander Brenning [EMAIL PROTECTED] - T +1-519-888-4567 ext 35783 Department of Geography and Environmental Management Unive

Re: [R-sig-Geo] Add Grid Values to Points

2008-10-31 Thread Alexander Brenning
rpolation [2] Inverse Distance Interpolation [3] Bicubic Spline Interpolation [4] B-Spline Interpolation 'add_grid.shp' is not a correct numeric value for option 'INTERPOL'. I can`t explain why R is trying to understand the expression "add_grid

Re: [R-sig-Geo] Tiled processing...

2008-10-30 Thread Alexander Brenning
ing -- am I really just creating a large in-memory matrix (ds2) and then writing the entire output image out to disk, or is this somehow writing line-by-line? If the former is true, how do I modify this to write within the loop, rather than all-at-once? Thanks! --j -- Alexander Brenning

Re: [R-sig-Geo] R: suggestion to write a more elegant code

2008-10-16 Thread Alexander Brenning
ula = ~^power[i]) Errore: unexpected symbol in: "for (i in 1:length(power)){ rsaga.grid.calculus(in.grids = c("DCM_1.sgrd"),out.grid =paste("DCM_1_power"power" } Errore: unexpected '}' in "}" -Messaggio originale- Da: Alexander Brenning [mailto:[EMAIL PROTEC

Re: [R-sig-Geo] suggestion to write a more elegant code

2008-10-16 Thread Alexander Brenning
1.sgrd"),out.grid = "DCM_1_power15.sgrd", formula = ~^3.6) rsaga.grid.calculus(in.grids = c("DSM_1.sgrd"),out.grid = "DCM_1_power15.sgrd", formula = ~^3.7) rsaga.grid.calculus(in.grids = c("DSM_1.sgrd"),out.grid = "DCM_1_power15.sgrd", formula = ~^3.8) rsaga.grid.calculus(in.gri

Re: [R-sig-Geo] R+SAGA formula to subtract two raster

2008-10-16 Thread Alexander Brenning
AIL PROTECTED],2])) # filter the missing values using neighbours: rsaga.geoprocessor(lib="grid_tools", module=7, param=list(INPUT="DSM_1.sgrd", RESULT="DSM_1_f.sgrd", THRESHOLD=0.1)) #DCM = DSM - DEM Rsaga.geoprocessor(???) Thanks for help Ale

Re: [R-sig-Geo] opening a grd file (Surfer)

2008-10-09 Thread Alexander Brenning
GEBCO website. Thanks Anne ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Alexander Brenning [EMAIL PROTECTED] - T +1-519-888-4567 ext 35783 Department of Geography and Environmental Management University o

Re: [R-sig-Geo] CODE for spatial logistic regression

2008-10-09 Thread Alexander Brenning
regression Hello I need some guide about spatial logistic regression, Is it available a code in R? {} LIOV ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Alexander Brenning [EMAIL PROT

Re: [R-sig-Geo] from ArcGIS GRID to .ASC

2008-09-19 Thread Alexander Brenning
s01jan, "mes01jan.asc", drivername = "AAIGrid", type="Byte") Any idea? miltinho, brazil [[alternative HTML version deleted]] ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/list