gy of
> the observed area. Does any of this make sense are i am way off?
>
> Any ideas
> and/or experiences with such a problem is greatly appreciated,
>
> Regards, m
> [[alternative HTML version deleted]]
>
>
--
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
530.754.7341
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
me objects; this in turn uses stack.data.frame
>
> Bests,
>
>
> On 11/08/2010 05:00 AM, Dylan Beaudette wrote:
> > Interesting demonstration. However, it seems like sp objects and
> > spplot are still more efficient at storing / plotting large grids. The
> > example listed be
Interesting demonstration. However, it seems like sp objects and
spplot are still more efficient at storing / plotting large grids. The
example listed below required over 2 minutes to complete on a 2.4 Ghz
Intel processor. Also, storing the coordinates + data in long format
could easily fill availa
t may be
> involved to do so. Plus, I don't know anything about it. I was hoping to
> use a native R solution. I do have GDAL and rgdal installed.
>
> Scott
Understandable. There is a fairly simple Windows installer available. But if
it is an R-based solution you are lo
with datasets as
large as 14.5 billion points (379GB)."
Some more information from the GRASS Wiki:
http://grass.osgeo.org/wiki/LIDAR
Cheers,
Dylan
--
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
530.754.7341
__
= rbind(ret, coordinates(pp[[j]]))
> > }
> > ret
> > }
> > q = allcoordinates(xx)
> >
> > # check:
> > plot(xx)
> > lines(q, col='blue')
>
> And that's basically what fortify does, except it covers a few more cases.
>
>
list.of.clusters, function(x)
> > {plot(clust.geo.test[which(my.cluster == x),], add = T)})
> >
> > would get you to the point I am now.
> >
> > Hans Ekbrand
> >
> > -BEGIN PGP SIGNATURE-
> > Version: GnuPG v1.4.9 (GNU/Linux)
> >
> > iEYEARECAAYFAkvRsLEACgkQfCyHKnBQYU4JrQCfWMlOgECamin5ZugGoHcUpqx3
> > fuQAoNDoG6iMVlLJCOYfoP5MRbNYHGlI
> > =jKTU
> > -END PGP SIGNATURE-
> >
> > ___
> > R-sig-Geo mailing list
> > R-sig-Geo@stat.math.ethz.ch
> > https://stat.ethz.ch/mailman/listinfo/r-sig-geo
--
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
530.754.7341
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Hi Thomas,
check out the generic lattice argument "par.strip.text". It accepts a
list with common formatting options.
Cheers,
Dylan
On Thu, Apr 22, 2010 at 8:32 AM, Thomas Adams wrote:
> All:
>
> I am using spplot to automatically generate a 21-panel eps graphic within a
> script where the titl
threshold=1000
... to reduce the number of vertices in a map of CA counties, so that I could
plot it with spplot().
Cheers,
Dylan
--
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
530.754.7341
___
nks,
Dylan
--
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
530.754.7341
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
s$name=="GTiff"]
> [1] TRUE
> gdaldrivers$copy[gdaldrivers$name=="GTiff"]
> [1] TRUE
>
> Am I perhaps missing an rgdal installation flag that is required to
> enable GRASS rasters?
>
> --j
Could be. How did you install GRASS/GDAL/rgdal?
Dylan
>
>
Roger can likely point out the problem after some example code is
posted. In the mean time I would suggest checking to make sure your
GDAL is compiled _without_ GRASS support, and that the GDAL GRASS
plugin is correctly installed.
This is how I check on my machine:
gdal_translate --formats | grep
On Monday 22 February 2010, Edzer Pebesma wrote:
> Dylan Beaudette wrote:
> > On Monday 22 February 2010, Edzer Pebesma wrote:
> >> Carlos, please try:
> >>
> >> library(gstat)
> >> data(meuse)
> >> meuse$rn = 1:nrow(meuse)
> >> coordi
late record
> > number with nmax=1 and compute a table of the resulting "prediction"
> > grid'.
> >
> > Could you provide me some examples for each one of this methods for an
> > R user,please?.
> >
> > Regards,
> >
> > Carlos Gutierrez C
time series? Or would I be better to stick with the Stineman algorithm
> for the time-series interpolation?
>
> Any feedback on this (including any warnings that I might be abusing
> things with my time series and IDW approaches) would be very much
> appreciated. I am definitely in
On Monday 08 February 2010, Karl Ove Hufthammer wrote:
> On Sun, 7 Feb 2010 17:30:23 -0800 Dylan Beaudette
>
> wrote:
> > How about an equal-area map projection?
> >
> > http://en.wikipedia.org/wiki/Map_projection#Equal-area
>
> Good idea. This will easily give y
How about an equal-area map projection?
http://en.wikipedia.org/wiki/Map_projection#Equal-area
Cheers,
Dylan
On Sun, Feb 7, 2010 at 1:37 PM, Matthew Vavrek
wrote:
> Hello,
> I have a question that I have yet to find a definitive answer for on this
> mailing list (or elsewhere as of yet) about a
max
> > r1 2.560776 6.403847
> > r2 49.493613 51.505419
> >
> > Thansk,
> > B. Coulmont
This type of work is most efficiently done in a standard GIS. I would suggest
installing GRASS GIS or something like that. This is a standard vector
operation, and should be v
TRUE)
>
>
>finalSPDF<-SpatialPixelsDataFrame(finalSPntsDF, data.frame(interp$z),
> tolerance = sqrt(.Machine$double.eps),
> proj4string = CRS(as.character(NA)))
>
> ...
> On a related note - I have tried using the 64 bit version of Revolution
> but u
On Sat, Oct 17, 2009 at 8:13 AM, Barry Rowlingson
wrote:
> On Sat, Oct 17, 2009 at 3:46 PM, Ezra Boyd wrote:
>
>> Will I be able to complete the tasks described above using the the Open
>> Source GIS without a bunch of time consuming complications? Or, should I
>> fork over some cash to extend m
t you are looking for. The resulting file is a
raster, decimated to some set resolution. Some screen shots:
http://hamish.bowman.googlepages.com/grassfiles#xyz
Cheers,
Dylan
--
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
t providing pdf not working?
> >
> > Thanks a lot
> >
> > Matthieu Stigler
> >
> > ___
> > R-sig-Geo mailing list
> > R-sig-Geo@stat.math.ethz.ch
> > https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
> __
> Wesley
> >
> >
> > Wesley Roberts MSc.
> > Researcher: Earth Observation
> > Natural Resources & the Environment (NRE)
> > CSIR
> > Tel: +27 (0)21 888-2490
> > Fax: +27 (0)21 888-2693
> > "To know the road ahead, ask tho
.frame(x=sample(1:10, 5), y=sample(1:10, 5))
dist(d)
1234
2 3.605551
3 6.403124 7.071068
4 9.219544 7.211103 5.830952
5 8.062258 5.099020 7.211103 3.162278
Cheers,
Dylan
--
Dylan Beaudette
Soil Resource
r an approach towards automatic trend detection: we may be able to skip
parametric specification of the trend (i.e. poly(2,3,4,5, ...)) by using some
kind of smoother: restricted cubic splines comes to mind (see rcs() in the
Design package). The script could then check the various elements of an RC
; ***
> Dr. Jeffrey W. Hollister
> US EPA
> Atlantic Ecology Division
> 27 Tarzwell Drive
> Narragansett, RI 02882
> (401) 782-9655
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@stat.math.ethz.ch
> https:
stat.math.ethz.ch
> >https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
> ***
> Marcelino de la Cruz Rot
> Depto. Biologia Vegetal
> EUIT Agricola
> Universidad Politecnica de Madrid
> 28040-Madrid
> SPAIN
> ***
On Friday 05 June 2009, Roger Bivand wrote:
> On Fri, 5 Jun 2009, Dylan Beaudette wrote:
> > On Fri, Jun 5, 2009 at 1:55 AM, Roger Bivand wrote:
> >> On Thu, 4 Jun 2009, Dylan Beaudette wrote:
> >>> On Monday 01 June 2009, Roger Bivand wrote:
> >>>&
On Fri, Jun 5, 2009 at 1:55 AM, Roger Bivand wrote:
> On Thu, 4 Jun 2009, Dylan Beaudette wrote:
>
>> On Monday 01 June 2009, Roger Bivand wrote:
>>>
>>> On Mon, 1 Jun 2009, Dylan Beaudette wrote:
>>>>
>>>> On Friday 29 May 2009, Torleif Marku
On Monday 01 June 2009, Roger Bivand wrote:
> On Mon, 1 Jun 2009, Dylan Beaudette wrote:
> > On Friday 29 May 2009, Torleif Markussen Lunde wrote:
> >> Hi
> >>
> >> First try updating glibc, and install spgrass6 again. If this does not
> >> work y
the same error... Very strange, as I did not encounter
this problem on a very similar machine.
Could this have anything to do with configure arguments such as -march ?
Cheers,
Dylan
> On Friday 29 May 2009 11:26:56 pm Dylan Beaudette wrote:
> > Hi,
> >
> > Just upgraded to
local/lib/R/bin/exec/R: double free or corruption
(fasttop): 0x096ac2b0 ***
At this point R stops responding, and I have to issue a `kill -9` to make
regain control of the shell.
Any ideas?
Dylan
--
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.
On Thu, May 28, 2009 at 1:21 AM, Anne Ghisla Insubriae
wrote:
> -BEGIN PGP SIGNED MESSAGE-
> Hash: SHA1
>
> Dylan Beaudette ha scritto:
> [snip]
>> Hi. It is alway nice to see people working GRASS SoC projects, thank you!
>>
>> I would just add-- I thi
On Sat, May 23, 2009 at 10:59 AM, Anne Ghisla Insubriae
wrote:
> -BEGIN PGP SIGNED MESSAGE-
> Hash: SHA1
>
> Ebrahim Jahanshiri a écrit:
>> Dear Anne,
>
> Ebrahim,
>
>> Thanks for the enquiry. I think we should have done this long before. I
>> havent done any geostatistics with GRASS but I
See the clara() function from the cluster package. It scales fairly
well to larger-sizes data sets.
Cheers,
Dylan
On Mon, May 11, 2009 at 5:35 AM, Wesley Roberts wrote:
> Hi Dan,
>
> Thanks for the advice. I want to classify my data into three classes; canopy,
> non-canopy and ground based on s
Hi,
What version of GRASS are you using? Have you asked on the GRASS list
regarding the slow r.clump [-ing] ?
Cheers,
Dylan
On Sun, May 10, 2009 at 2:30 PM, milton ruser wrote:
> Dear all,
>
> sorry for this double-posting, I know it is not a good idea,
> but may be this list is more appropria
Hi,
You can use the v.generalize command in GRASS to reduce the complexity
of vectors like these.
Cheers,
Dylan
On Sun, Apr 26, 2009 at 3:53 PM, hadley wickham wrote:
> Hi Enrico,
>
> I have some code to do map generalisation (reducing map resolution
> without visible loss in detail) at
> http:
j in 1:length(columns)) {
> >>>> sset<-subset(wind, wind$vt==times[i], select=c(columns[j]))
> >>>> dOut[i,columns[j]]<-gks(lat,lon,sset,columns[j])
> >>>> }
> >>>> }
> >>>> dOut<-cbind(dOut, m
this possible?
>
> brian Cooper
>
> [[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
>
> ___
> R-sig-Geo mailing list
&
On Wed, Feb 25, 2009 at 12:53 AM, Roger Bivand wrote:
> On Tue, 24 Feb 2009, Heuvelink, Gerard wrote:
>
>> Dear list,
>>
>> The stratified random sampling problem that I submitted a few days ago has
>> already been solved, with the help of several of you, notably Edzer Pebesma.
>> Edzer came up wi
this help:
http://casoilresource.lawr.ucdavis.edu/drupal/node/644
Cheers,
Dylan
--
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
530.754.7341
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
ystem
2) the data all contain enough SRS information to convert to a common
coordinate system
I'd ask on the GDAL list.
Dylan
--
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
530.754.7341
___
On Wed, Jan 14, 2009 at 1:14 AM, Roger Bivand wrote:
> On Tue, 13 Jan 2009, Dylan Beaudette wrote:
>
>> On Tuesday 13 January 2009, Edzer Pebesma wrote:
>>>
>>> This sounds very familiar; I ran into this about 6 years ago or so.
>>
>> Coul
u.Max.NA's
> >> -32770284465167718 11600 32770 6483714
> >>
> >> The results from readAsciiGrid are pretty close (Max = 37,080) to the
> >> original grid, but readGDAL results in much different values. I would
> >> expect both functions to create an object with exactly the same values
> >
> > as
> >
> >> the Ascii grid. Can anyone explain why this occurs, and if there is a
> >> different way to get my grid into R (I eventually need an 'im' to work
> >> with)? I am new to spatial analysis using R, so if there is something
> >> obvious I am missing I apologize. Thanks in advance for any help!
> >>
> >> Cheers,
> >> -Brian
> >>
> >>
> >> ===
> >> Brian R. Miranda
> >> USDA Forest Service
> >> Northern Research Station
> >> brmira...@fs.fed.us
> >>
> >> Institute for Applied Ecosystem Studies
> >> 5985 Highway K
> >> Rhinelander, WI 54501
> >> phone: 715-362-1186
> >> fax: 715-362-1166
> >>
> >>
> >>[[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
--
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
530.754.7341
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
On Tue, Dec 30, 2008 at 3:41 AM, francesco wrote:
> Dear List,
>
> I have been banging my head on a problem for three days now, so I believe the
> time has come to look for help.
>
> I have a polygon shapefile with an irregular grid and I would like to turn
> the polygons in the shapefile into t
graphy
> School for Marine Science and Technology
> University of Massachusetts
> 200 Mill Road, Suite 325
> Fairhaven, MA 02719
> jloeh...@umassd.edu
> T 508-910-6393
> F 508-910-6396
>
> _______
> R-sig-Geo mailing list
> R-
t;
> ___
> R-sig-Geo mailing list
> R-sig-Geo@stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
--
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
her solution, based on the definition of the geometric mean
w/respect to logarithms.
geom-mean = e^{ 1/n * sum{ln x} }
# only works for positive coordinates !!
x <- runif(n=100, min=1, max=100) ; y <- runif(n=100, min=1, max=100)
# plot points
plot(x,y)
# plot the arithmetic mean point
point
; > know it. Thanks.
> >
> >
> > Jim Bouldin, PhD
> > Associate in Agricultural Experiment Station
> > Department of Plant Sciences, UC Davis
> > Davis CA, 95616
> > 530-554-1740
> >
> > ___
> > R-si
he world of spatial applications of R and am
> > wondering about this trivial question:
> >
> > What is the easiest way to subset a SpatialGridDataFrame, based on 4
> > corner coordinates?
> >
> >
> >
> > The reason I ask is that I want to do an o
On Fri, Oct 24, 2008 at 8:24 AM, Edzer Pebesma
<[EMAIL PROTECTED]> wrote:
> Tomislav Hengl wrote:
>>
>> In GIS terms, getting a raster from one grid topology to (any) other grid
>> topology (also subgrids)
>> is referred to 'resampling' (I do not think that this is yet implemented
>> in any R packa
oesn't yet have the packages.
>
> Thanks for your help.
>
> Regards,
> Malin
Also-- some tips on iteratively sampling several 'geometries' contained in a
single file.
http://casoilresource.lawr.ucdavis.edu/drupal/node/644
Cheers,
Dylan
--
Dylan Beaudette
S
ps)
>
> So far, my searches on r-sig-geo have come up negative. Has anyone
> dealt with data like this before? Can someone recommend the
> appropriate packages/functions? PostGIS sounds like it may be another
> option if R doesn't yet have the packages.
>
> Thanks fo
data...
>
> Well, yes, in that case animation becomes compelling. Aren't there any
> GIS clients that can do animation with time series of vector maps? I
> only know one for gridded data.
If you are just making maps, and lots of them GRASS might be useful for
automating the proces
7;10s70s0e150e.sgrd'))
> This are bathymetry data from the GEBCO website.
>
> Thanks
> Anne
>
GDAL
--
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
530.754.7341
_
On Wednesday 08 October 2008, Roger Bivand wrote:
> On Wed, 8 Oct 2008, Dylan Beaudette wrote:
> > On Tuesday 07 October 2008, Dylan Beaudette wrote:
> >> On Tuesday 07 October 2008, Roger Bivand wrote:
> >>> On Mon, 6 Oct 2008, Dylan Beaudette wrote:
> >>&
On Tuesday 07 October 2008, Dylan Beaudette wrote:
> On Tuesday 07 October 2008, Roger Bivand wrote:
> > On Mon, 6 Oct 2008, Dylan Beaudette wrote:
> > > Hi,
> > >
> > > I have noticed that saving data to files that include a DBF, result in
> > > bogus
; >>
> >> LIOV
> >
> > You can try the functions in package geoRglm.
> > You can set binomial counts with spatial coordinates and fit a
> > geostatistical logistic model in geoRglm.
> > HTH
> > R.
--
Dylan Beaudette
Soil Resource Laboratory
ht
On Tue, Oct 7, 2008 at 6:26 PM, Hamish <[EMAIL PROTECTED]> wrote:
> Dylan:
>> It looks like the limiting factor in this equation is the
>> code used in v.out.ogr.
>
> maybe a silly question, but is a 3rd party format even needed here?
>
> $ ogrinfo --formats | grep -i grass
> -> "GRASS" (readonly)
On Tuesday 07 October 2008, Roger Bivand wrote:
> On Mon, 6 Oct 2008, Dylan Beaudette wrote:
> > Hi,
> >
> > I have noticed that saving data to files that include a DBF, result in
> > bogus data where there were NA. Using the write.dbf() function from
> > the forei
Hi,
I have noticed that saving data to files that include a DBF, result in
bogus data where there were NA. Using the write.dbf() function from
the foreign package seems to work a little better, but I still get odd
results in numeric columns. Writing to GRASS with the methods in the
spgrass6 packag
their actual values. Have a look at this
thread:
https://stat.ethz.ch/pipermail/r-sig-geo/2007-October/002694.html
Cheers,
Dylan
--
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
530.754.7341
On Friday 03 October 2008, Dylan Beaudette wrote:
> On Thursday 25 September 2008, Edzer Pebesma wrote:
> > demo(uisim)
>
> This seems work work with the demo code, but not with another example using
> my own data:
>
> # setup
> library(gstat)
> library(spgrass6)
&g
prediction grid:
data(meuse.grid)
gridded(meuse.grid) = ~x+y
# define variable as dummy data
v = vgm(.25, "Sph", 900)
g = gstat(NULL, "var1", x~1, beta = .5, nmax = 20, model = v, dummy = TRUE)
# simulation of a single variable
out = predict(g, meuse.grid, nsim = 20, indicat
ify the
> periodicity/regularity of the patches in said landscapes.
>
>
>
> Can anyone point me in the right direction?
>
Forgot to mention the RandomFields package.
require(RandomFields)
example(soil)
--
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis
On Wednesday 24 September 2008, Satterthwaite, Paul wrote:
> Hello,
>
>
>
> I am a very new user of R, but have used S-Plus for about a year or so.
> I am curious about methods of creating simulated patch landscapes in R,
> and how to parameterize patch cover values and modify the
> periodicity/reg
RS as 'proj=lonlat datum=...'
3. invoke the projection method, with the target projection = "+proj=lcc
+a=637000 +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=-7.5"
4. the result should be in the target projection.
Cheers,
Dylan
--
Dylan Beaudette
Soil Resource Laboratory
http://caso
dows installer for GRASS, or you
can just install the latest QGIS: this comes with GRASS!
Cheers,
Dylan
--
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
530.754.7341
___
R-sig-Geo mail
On Mon, Jul 7, 2008 at 2:48 AM, Paulo Cardoso <[EMAIL PROTECTED]> wrote:
> Can R perform matrix algebra between ascii converted Esri grids? I'm
> convinced that its trivial.
>
>
>
> The question is a little more complex (from my point of view and
> inexperience): I'm interested in perform cost surf
If you are interested in a (supervised) imagery classification routine
that takes spatial arrangement into consideration, check out the
i.smap command in GRASS GIS.
Cheers,
Dylan
On Thu, Jun 12, 2008 at 4:57 AM, Laura Poggio <[EMAIL PROTECTED]> wrote:
> thank you. It seems now solved!
>
> Laura
On Tue, Apr 22, 2008 at 6:49 AM, Tomislav Hengl <[EMAIL PROTECTED]> wrote:
>
> Dear list,
>
> I know that much has already been said about the memory limit problems. If
> there is any progress
> about this problem, we would be interested to hear.
>
> In our project, we are importing 24 maps/ba
On Wednesday 12 March 2008, Roger Bivand wrote:
> On Wed, 12 Mar 2008, Dylan Beaudette wrote:
> > On Thursday 06 March 2008, Thomas Adams wrote:
> >> Dylan,
> >>
> >> I think a solution using GRASS can be found on pages 110-111 of "Open
> >> S
using two rasters with spplot() and see what happens.
Cheers,
Dylan
> Dylan Beaudette wrote:
> > On Tuesday 04 March 2008, Edzer Pebesma wrote:
> >> I find it hard to imagine how you want to plot two raster files on top
> >> of each other. Do you want some form of t
e using some kind of union
> process (I don't know how to do this in R) to show the quantitative raster
> data in the categorical classes. Then you can delete the larger context
> raster which does not contain any information and reduce the window
> (spatial extent) of the overlai
lly covers
about 30% of the region of interest. The other raster is just contextual
data, and thus would be useful to plot "behind" the first raster.
Ideas?
Dylan
> Dylan Beaudette wrote:
> > Hi,
> >
> > Is it possible to plot two raster images using spplot()
rial photo, used only for context,
while the second one is one with interesting z-values. The second raster is
masked and thus does not cover the entire region.
Thanks,
Dylan
--
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of
On Tuesday 26 February 2008, Roger Bivand wrote:
> On Tue, 26 Feb 2008, Dylan Beaudette wrote:
> > Hi,
> >
> > Is there a simple approach to selecting random "blocks" using the sp
> > classes and methods?
> >
> > Say I have a 10x10 grid of poly
(I hope!).
Cheers,
Dylan
--
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
530.754.7341
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
ocessor(lib="geostatistics_kriging", module=3,
> param=list(GRIDS="DEM.sgrd;SLOPE.sgrd;PLANC.sgrd;TWI.sgrd;SINS.sgrd;SMU1.sg
>rd;SMU3.sgrd;SMU4.sgrd;SM U9.sgrd", GRID="SOLUM_rk.sgrd",
> SHAPES="baranja.shp", FIELD=0, MODEL=1, NUGGET=0, SILL=200, RANGE=500,
> INTERPOL=0))
> # Does not work yet. Possibly a bug in the saga_cmd.exe?
>
> ***
>*
>
> The complete script and datasets are available at:
>
> http://spatial-analyst.net/GRK/examplesRSAGA.zip (400 KB)
>
> So the only real problem is the import/export from R to SAGA, which I guess
> could be solved very easily if the next version of rgdal would support SAGA
> format.
>
>
This is all very interesting, but doesn't the GRASS-R combination already do
these things- and quite well ? As far as I can tell GRASS can handle the
massive grid operations, and R+gstat can do the statistical modeling, etc.
But maybe I should check out SAGA again- it wouldn't compile last time...
Thanks for the post,
--
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
530.754.7341
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
lylist but probably I need to have "consistant"
> > polygon for each element of the list ?
> > -is there any distance method supported by a class "polygon" which can
> > compute the distance to another geometry?
> >
> > thanks
>
Not sure if it
On Friday 25 January 2008, Roger Bivand wrote:
> On Thu, 24 Jan 2008, Dylan Beaudette wrote:
> > On Thursday 24 January 2008, Roger Bivand wrote:
> >> On Thu, 24 Jan 2008, Dylan Beaudette wrote:
> >>> Hi Roger. Thanks for contributing some answers to this.
> >&g
On Thursday 24 January 2008, Roger Bivand wrote:
> On Thu, 24 Jan 2008, Dylan Beaudette wrote:
> > Hi Roger. Thanks for contributing some answers to this.
> >
> > I was recently working with a colleague on developing some sample
> > exercises for new students. Since join
# some example codes
veg_codes <- data.frame(code=1:4, meaning=c('code 1','code 2','code 3','code
4'))
# join the original data table with the veg codes table
combined <- merge([EMAIL PROTECTED], y=veg_c
ailing list
> > R-sig-Geo@stat.math.ethz.ch
> > https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
--
Dy
On Friday 07 December 2007 11:13:17 am Markus Loecher wrote:
> Dear geo list members,
> does there exist an R package or might one of you have an algorithm that
> partitions a given area into 2dimensional "bins" (could be rectangles,
> squares, polygons) adaptively, i.e. modifying the coarseness of
On Friday 07 December 2007 05:20:59 am Anne Ghisla wrote:
> Dear all,
>
> during kernel LSCV calculation for home range estimation (kernelUD, package
> adehabitat) I encountered a problem related to topology.
> One homerange is made of 5 polygons, but one of them has a null area: it is
> in fact a
try and post an examples of this.. in the mean time check out Roger's
sp website- it should be in the manual page for the sp package. There are
numerous mapping examples in there.
For more complex maps I tend to favor GMT. Examples:
http://casoilresource.lawr.ucdavis.edu/drupal/node/130
&g
On Saturday 24 November 2007 02:04:17 pm Tom Sgouros wrote:
Hi Tom, sorry to hear that you haven't discovered what you were looking for
yet. Here are some comments / suggestions.
> Hello all:
>
> I was referred to this list when I asked a question on the r-help list
> about mapping. Unfortunate
7;ogr2ogr' useful in this context. Here are some examples:
# convert shapefile whith defined projection to UTM Zone 10:
ogr2ogr -t_srs "+proj=utm +zone=10 +datum=NAD83" utm.shp original.shp
Cheers,
Dylan
--
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
530.754.7341
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
On Friday 26 October 2007, Roger Bivand wrote:
> On Fri, 26 Oct 2007, Dylan Beaudette wrote:
> > On Friday 26 October 2007, Roger Bivand wrote:
> >> On Fri, 26 Oct 2007, Dylan Beaudette wrote:
> >>> Hi,
> >>>
> >>> noticed that when I impor
Hi,
noticed that when I import a shapefil with readOGR() the coordinates are
truncated to integers. Is there some way to prevent this type of behaviour?
cheers,
Dylan
--
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
#x27;m not sure how much sense
> they make without my witty narrative that goes along with it :) Did
> anyone record my talk for posterity?
>
> Link here:
>
> http://www.foss4g2007.org/presentations/view.php?abstract_id=135
>
> Barry
>
Excellent presentation, on an exce
92 matches
Mail list logo
