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
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
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
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
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'
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
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
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...
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
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
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
_
[[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
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,
= 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
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
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
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
dOGR(".", "theme1")
my_SPDF$ab <- my_SPDF$a - my_SPDF$b
writeOGR(my_SPDF, ".", "mynewfile", driver="ESRI Shapefile")
___
R
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
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
___
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
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
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
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
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
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
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
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
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
29 matches
Mail list logo