ersity
> C-1-2-225, Katsura Campus, Kyoto University,
> Nishikyo-ku, 〒615-8530, Kyoto, Japan
> --
> View this message in context:
> http://r-sig-geo.2731867.n2.nabble.com/can-we-use-georob-package-to-make-a-3D-nested-variogram-tp7587522p7587524.html
> Sent from the R-sig-geo mailing li
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Edzer Pebesma, Co-Editor-in-Chief Computers & Geosciences
Institute for Geoinformatics (ifgi), University of Münster
Heisenbergstraße 2, 48149 Münster,
SE)
> plot(migmap,add=TRUE)
> plot(log10(g2GM),zlim=c(1,5),col=colores, breaks=seq(1,5,by=0.5),
> legend=FALSE, add=TRUE)
> #' the map gets distorted depending upon de geometry of the window
>
>
> Thanks,
>
> Agus
>
> _
xecution halted
What does this mean? Why does gUnionCascaded(sp1) not do this?
--
Edzer Pebesma, Co-Editor-in-Chief Computers & Geosciences
Institute for Geoinformatics (ifgi), University of Münster
Heisenbergstraße 2, 48149 Münster, Germany. Phone: +49 251
83 33081 http://ifgi.uni-muenster.de GPG key I
On 01/05/2015 01:53 PM, Roger Bivand wrote:
> On Mon, 5 Jan 2015, Edzer Pebesma wrote:
>
>> When running the following script:
>>
>> library(sp)
>> gt = GridTopology(cellcentre.offset = c(0,0), cellsize = c(1,1),
>> cells.dim = c(3,3))
>>
ort I did this:
>>>
>>> for (i in SB) {
>>> jpeg(paste(SB[[1]],".jpg",sep=""))
>>> bubble(SB[[i]], "mean.Evade.", col=c("#00ff0088", "#00ff0088"),
>>> main = "Evasión" )
>
correct.diagonal = 1.01)
>
> # Perform co-kriging predictions
> ck <- predict.gstat(g1, grid)
>
> Intrinsic Correlation found. Good.
> [using ordinary cokriging]
>
> ERROR: "chfactor.c", line 131: singular matrix in function LDLfactor() Error
>
; Yours sincerely / Med venlig hilsen
>
>
> Frede Aakmann Tøgersen
> Specialist, M.Sc., Ph.D.
> Plant Performance & Modeling
>
> Technology & Service Solutions
> T +45 9730 5135
> M +45 2547 6050
> fr...@vestas.com
> http://www.vestas.com
>
> Comp
compute model values and pass these as model "Tab" to
vgm(), (see ?vgm), or you can write your own cubic model R function and
pass that on to gstat::krige0 (which has more limited capabilities than
gstat::krige).
The meaning of your Structure 1 and 2 listed above requires working
k
TML version deleted]]
>>>
>>> ___
>>> R-sig-Geo mailing list
>>> R-sig-Geo@r-project.org
>>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>
>> ___
>
oints")) coords = coordinates(coords)
>
> When coords is a dataframe, it returns a matrix with the rownames lost.
> Subsequent matching with data is not longer possible. Should'nt the test
> be reversed: if (is(coords, "SpatialPoints")) coords =
> coordinates(
he message. Please let me know how to hide it.
You can't, because it uses Rprintf(). I've changed this to use the R
warning mechanism, so in future versions you can.
--
Edzer Pebesma, Co-Editor-in-Chief Computers & Geosciences
Institute for Geoinformatics (ifgi), University of Mün
om.
>
> _______
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Edzer Pebesma
Institute for Geoinformatics, University of Münster,
Heisenbergstraße 2, 48149 Münster, Germany; +49 251 83
;
>
>
> --
> View this message in context:
> http://r-sig-geo.2731867.n2.nabble.com/block-and-output-grid-size-in-block-kriging-tp7587770p7587773.html
> Sent from the R-sig-geo mailing list archive at Nabble.com.
>
> _______
>
>
>
>
>
>
>
>
> diff= m.bk - est.bk
>
>
>
> --
> View this message in context:
> http://r-sig-geo.2731867.n2.nabble.com/block-and-output-grid-size-in-block-kriging-tp7587770p7587780.html
> Sent from the R-sig-geo mailing list archive at N
-output-grid-size-in-block-kriging-tp7587770p7587782.html
> Sent from the R-sig-geo mailing list archive at Nabble.com.
>
> _______
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Edzer Pebesma
Institute for
59235024 -37.75000 83.58333
322 64257224923 -36.91667 83.58333
422 64257 224923 -36.75000 83.58333
522 64257224923 -36.58333 83.58333
--
Edzer Pebesma
Institute for Geoinformatics (ifgi), Unive
S
>
>
> Thanks in advance for any answer,
>
> Mathieu
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
, or are
you interested in an integral along a line, i.e. is the problem purely
spatial? Is your grid three-dimensional? Is it a grid in latitude/longitude?
--
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster,
Heisenbergstraße 2, 48149 Münster, Germany; +49 251 83 33081
Jo
d Qu.:3322247 3rd Qu.:19.00 WBAR :10095 3rd Qu.:2003 3rd Qu.:1
> Max. :3367519 Max. :19.00 MAIZ : 5069 Max. :2005 Max. :1
>(Other):11061
> nuts_code
> FR25:192621
>
> ___
anée
> *
> Centre d’études et d’expertise sur les risques, l’environnement, la
> mobilité et l’aménagement
> www.cerema.fr <http://www.cerema.fr>
>
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://sta
Google Earth: 7.1.2.2041
> Build-Date: 10/7/2013
> Build-Time: 12:17:00 nachmittags
> Renderer: OpenGL
> OS: Linux (3.13.0.0)
> Graphics Driver: Intel Open Source Technology Center
> -
>
> *Any ideas how to solve this?*
>
> Best,
>
> Guido
>
&g
omething else going on that should not be going
> on. This is where Frede's suggestions might help. You will need to
> carefully track the construction of your locs object and see if it is
> somehow different in the two situations. I don't know of any "designed
> circ
gt;
>> pr1 <- projectRaster(r, crs=newproj)
> Error: 'rawTransform' is not an exported object from 'namespace:rgdal'
>
> Thank you!
>
> T. Hengl
>
> _______
> R-sig-Geo mailing list
> R-sig-Geo@r-proj
r <- spplot(r,
>> sp.layout = list( sp.lines, mp, first = FALSE)
>> )
>>
>> # show the map with the projected raster
>> plot_pr <- spplot(pr,
>> sp.layout = list( sp.lines, mp, first = FALSE)
>> )
>>
>> print(plot_r)
>> print(pl
e_with_hole <- wrld_simpl[5,]
if (avoidGEOS)
gpclibPermit()
shape_with_hole@polygons <- lapply(shape_with_hole@polygons,
checkPolygonsHoles,
avoidGEOS = avoidGEOS)
data(wrld_simpl)
all.equal(x, wrld_simpl)
returns a difference, but TRUE if avoidGEOS = TRUE.
--
Edzer Pebesm
read into R
> tMaxGrid <- readGDAL(fname = unzippedFiles[1])
>
> # delete the file
> file.remove(tmpFile)
> file.remove(unzippedFiles)
>
> # This doesn't work correctly
> states$tmax20140416 <- over(states, tMaxGrid)
this would assign a data.frame to what you want
ction(x) slot(x, name = "ID"))
text(coordinates(i), ids.i, cex = .8, col = 'blue')
# compute & plot area-weighted average of p1[[1]]
ret = aggregatePoly(p1, p2)
spplot(ret["value"])
>
> On Mon, Apr 13, 2015 at 10:47 PM, Edzer Pebesma <
> edzer.pebe...@
;
> Any hints would be greatly appreciated.
>
> Cheers,
> Juta
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r
erform a gaussian conditional simulation on a variable
> that is highly skewed.
> Could someone suggest me how to transform my variable in a gaussian
> distribution?
>
> Thanks,
> Carlo
>
>
--
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster,
H
d then look into copulas.
>
> Carlo
>
>
>
>
> 2015-04-28 22:32 GMT+02:00 Edzer Pebesma <mailto:edzer.pebe...@uni-muenster.de>>:
>
> one way is the rank transform:
>
> x = rgamma(100, 0.1)
> hist(x)
> y = qnorm((order(x)/(length(x) + 0.5)))
&g
you for your help in advance!
>
> Best regards,
>
> Istvan
>
> [[alternative HTML version deleted]]
>
> ___________
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Edze
n
stplot(Produc.st[,'1970/1972' , "unemp"], col.regions = SAGA_pal[[4]],
at = 3:11)
or else update to the spacetime version that fixed this with
devtools::install_github("edzer/spacetime")
>
> Thanks
>
> Edoardo
>
> [[alternative HTML version deleted]]
>
&g
message in context:
> http://r-sig-geo.2731867.n2.nabble.com/stplot-legend-classes-of-categorical-vairable-tp7587860p7588200.html
> Sent from the R-sig-geo mailing list archive at Nabble.com.
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo
ey-Godin
>> Ph.D. Student, Department of Biology, Dalhousie University
>> Industrial Graduate Fellow, WWF-Canada
>> Email: [hidden
>> email] | Web:
>> wormlab.biology.dal.ca<http://wormlab.biology.dal.ca>
>>
>> --
>> View this message in co
14597L,
> 11803L, 18666L, 8734L, 12069L, 14687L, 18756L, 13277L, 13795L,
> 20088L, 13112L, 19159L, 9675L, 7719L, 20285L, 4072L, 20252L,
> 1536L, 1740L, 20418L), class = "data.frame”)
>
> --
>
> Aurelie Cosandey-Godin, PhD
> Dalhousie University, Department of B
. .. .. .. ..@ hole : logi FALSE
> .. .. .. .. ..@ ringDir: int 1
> .. .. .. .. ..@ coords : num [1:5, 1:2] -17.05 -17.05 -8.67 -8.67 -17.05 ...
> .. ..@ plotOrder: int 1
> .. ..@ labpt: num [1:2] -14.3 10.5
> .. ..@ ID : chr "ESH"
> .. ..@ area
e now, so any hints regarding a
> fix for this would be much appreciated.
>
> Best
> Tim
>
--
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster,
Heisenbergstraße 2, 48149 Münster, Germany; +49 251 83 33081
Journal of Statistical Software: http://www.jst
e
> subsetting. In (probably miscoding) the analysis of gRelate of 200,000
> blocks and 2,000 buffer circles, I eat up 6gigs with ease.
>
> Best wishes, again,
> mike
>
> [[alternative HTML version deleted]]
>
> ___
, Edzer Pebesma wrote:
> Mike, you may want to read two recent threads:
>
> https://stat.ethz.ch/pipermail/r-sig-geo/2015-April/022624.html
> https://stat.ethz.ch/pipermail/r-sig-geo/2015-April/022663.html
>
> When you try to aggregate, instead of using over(), have you used
> s
> of samples in ms, duration between samples; from which
> an STI can be constructed. Not much more data than where the eyes were when.
> It would seem that there would be a lot of simple sensor
> data of this sort so I wonder if Track can relax its requirement of STIDF to
> al
;
> So I'm mind-boggled. Why is this happening? Is it just that JS is
> trendy? Is it that there is a surfeit of JS programmers? Is it
> speculative "because I can" development? Should R-spatial data look at
> this as evolution and consider the future of sp classes?
&g
nections code to track distance relative to
> the central fixation point. I assume that @connections$distance is
> the distance from one sampling point/time and the subsequent. And
> then check whether accepted trials were actually valid as well.
> And that will then
les = list(draw = TRUE)
> )
> #######
>
> Any ideas..?
>
> Luc
>
> [[alternative HTML version deleted]]
>
>
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.c
providing an excellent package.
>
> Best,
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Edzer Pebesma
Institute for Geoinfor
(not entirely woken up it seems: I meant "stfdf" where I wrote "stplot")
On 06/04/2015 12:44 PM, Edzer Pebesma wrote:
> try converting variable report into a factor, by
>
> stplot$report = factor(stfdf$report)
>
> and give them meaningful levels if you want t
t; getClassDef('watch_circ_1')
> Class "watch_circ_1" [in ".GlobalEnv"]
>
> Slots:
>
> Name: sptime endTime
> Class: Spatial xts POSIXct
>
> Extends: "ST"
>> ?what am I doing wrong here?
>
> This is prob
with the concept of multiple (user defined) endTimes, or
> it would merely say "NA endTime value not allowed", singular as to
> value, full stop. And with contains in mind and Pebesma (2008)
> "Customizing spatial data
> classes and methods" I start to see a way, ala Michae
;
>
>
> --
> View this message in context:
> http://r-sig-geo.2731867.n2.nabble.com/I-can-not-install-the-package-RHydro-on-my-computer-windows-tp7588448.html
> Sent from the R-sig-geo mailing list archive at Nabble.com.
>
> ___
&
data.frame(v = 1:16, w=5:20,
>> x=factor(1:16)), match.ID = FALSE)
>> p2 = as(gt2, "SpatialPolygons")
>>
>> # all-factor attributes: compute area-dominant factor level:
>> ret = aggregate(p1["x"], p2, areaWeighted = TRUE)
>> # E
HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster,
Heisenbergstraße 2, 48149 Münster, Germany; +49 251 83 33081
> 1280 Main Street West
>
> Hamilton, Ontario, Canada L8S 4M2
>
> phone: 905-525-9140, x. 20384
>
> vic...@mcmaster.ca
>
> www.humanities.mcmaster.ca/~vickup
>
>
>
>
>
>
> [[alternative HTML version deleted]]
>
> _
asses try to represent from the
package vignette / JSS paper.
>
> What am I doing wrong?
>
> Thanks,
> Dan
>
> Dr Dan Bebber
> Senior Research Fellow
> Biosciences
> University of Exeter
>
>
>
>
--
Edzer Pebesma
Institute for Geoinformatics (ifgi), Unive
ty of Agriculture Faisalabad
>> Department of Math & Stat.
>>
>> [[alternative HTML version deleted]]
>>
>> ___
>> R-sig-Geo mailing list
>> R-sig-Geo@r-project.org
>> https://stat.ethz.ch/mailman
> I then tried to use gUnaryUnion with the variable STI (which I expected
> to return the unchanged polygons), which apparently returned exactly the
> same outline as the one obtained with BRFW. Presumable, I am overlooking
> something basic here, but I fail to see what. Can someone help?
pected "Stimmbezirke"
>> BW_EU_BRFW <- gUnionCascaded(BW_EU, "BRFW") ## gUnaryUnion does not
>> make a difference
>> plot(BW_EU_BRFW) ## plots only a single outline
>>
>> I then tried to use gUnaryUnion with the variable STI (which I
>&g
t; Ajoke
>
Please tell the list more what your data are about, what your problem
is, what you have tried so far, and explain why that did not work.
--
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster,
Heisenbergstraße 2, 48149 Münster, Germany; +49 251 83 33081
Journal
ey are of varying extents and resolutions for each time.
> Kind regards,
> Ajoke
>
>
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.eth
tried:
>> x1 = aggregate(pioggieST[,"2008"], "month", mean)
> There were 50 or more warnings (use warnings() to see the first 50)
>> warnings()
> Messaggi di avvertimento:
> 1: In mean.default(X[[i]], ...) :
> argument is not numeric or logical: returning NA
Qu.:2101
> Median :2002-07-02 Median :4201
> Mean :2002-07-02 Mean :4201
> 3rd Qu.:2008-04-01 3rd Qu.:6301
> Max. :2013-12-31 Max. :8401
> [[Data attributes:]]
> values
> Min. : 0.000
> 1st Qu.: 0.000
> Median : 0.000
> Mean :
http://blogs.esri.com/esri/esri-insider/2015/07/20/building-a-bridge-to-the-r-community/#sthash.L6pRPmqc.hpvt
--
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster,
Heisenbergstraße 2, 48149 Münster, Germany; +49 251 83 33081
Journal of Statistical Software: http
com/
>
> The full introductory article including interactive examples is
> published here:
>
> http://environmentalinformatics-marburg.github.io/web-presentations/20150723_mapView.html
>
> I hope this may prove useful for some of you.
> Also, if anyone is keen to get involved in
On 07/24/2015 04:52 PM, MAURICIO ZAMBRANO BIGIARINI wrote:
> On 24 July 2015 at 09:26, Edzer Pebesma wrote:
>> Hi Tim, thanks for starting the discussion.
>
> Thanks Tim for this nice implementation and to all of you for this
> very interesting discussion.
>
>
")
>
>
> It is possible to limit number of digits?
>
Pierluigi, I believe that the files written are binary files, so there
are no number of digits, really. Of course you can apply round() to your
data, prior to exporting, to round your data. Does it concern
attributes, or coordinates?
-
ng but now the colour is not graduated any more but.
> How can I have a graduated palette again?
> Thanks
>
this gave you a warning, which you didn't report.
Use cuts=5 to get 5 classes; use at=c(1,2.5,3) to get colour changes at
these values. And do read ?spplot.
--
Edzer P
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster,
Heisenbergstraße 2, 48149 Münster, Germany; +49 25
Gate : https://www.researchgate.net/profile/Tim_Peterson7
> Google Scholar:
> http://scholar.google.com.au/citations?user=kkYJLF4J&hl=en&oi=ao
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
--
elbourne, 3010 Australia
> T: +61 3 8344 9950, M: +61 0438 385 937
>
> Dept. profile :
> http://www.ie.unimelb.edu.au/people/staff.php?person_ID=141135
> Research Gate : https://www.researchgate.net/profile/Tim_Peterson7
> Google Scholar:
> http://scholar.google.com.au/citat
On 08/28/2015 05:15 AM, Tim Peterson wrote:
>
> On 27/08/15 16:55, Edzer Pebesma wrote:
>>
>> On 08/27/2015 02:45 AM, Tim Peterson wrote:
>>>
>>> Hi all,
>>>
>>> I'm considering editing the GStat kriging file 'predict.c'
On 08/29/2015 10:55 AM, Tim Appelhans wrote:
>
> I am not sure why SpatialPolygons* objects get rendered via levelplot.
because levelplot has the understanding of plotting something in 2D,
using colours according to the continuous variation of a third variable.
--
Edzer Pebesma
Institu
> is input to the Variogram call. Now, while I use *fit.lmc()*, I get
> an error message *"gstatVariogram var1 not present"*. I have done
> some preliminary findings, and there are no pertinent messages in
> this message board that have Raster as an inp
alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster,
Heisenbergstra
, if I think about it it confuses me, so I best not think about that
> first thing in the morning!
xyplot has the concept of one symbol per observation, with a colour, a
size, a symbol type. That is how I thought of plots of
SpatialPointsDataFrame at the time I wrote it.
--
Edzer Pebesm
earchgate.net/profile/Tim_Peterson7
> Google Scholar:
> http://scholar.google.com.au/citations?user=kkYJLF4J&hl=en&oi=ao
>
> [[alternative HTML version deleted]]
>
> _______
> R-sig-Geo mailing list
> R-sig-Geo@
# Re-attach the attribute table to the SPDF
> states@data <- states.df
now, you created an invalid object:
> length(geometry(states))
[1] 56
> dim(states@data)
[1] 55 12
this confirms once more what has been said on this list so often, that
instead of using constructor function
t; m = merge(states, counts) # uses sp::merge.Spatial
> m = m[-which(m$STUSPS == "PR"),]
> length(m)
[1] 55
> Thanks again!
>
> On Fri, Sep 4, 2015 at 9:51 AM Edzer Pebesma
> mailto:edzer.pebe...@uni-muenster.de>>
> wrote:
>
>
>
> On 09/0
Gaussian model.
>
> Best Regards
> Sherien
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Edzer Pebesma
Institute
niversity of Wollongong, Australia
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Edzer Pebesma
Institute for Ge
e that you run update.packages(), or
else make sure you update/reinstall raster AFTER you update sp.
--
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster,
Heisenbergstraße 2, 48149 Münster, Germany; +49 251 83 33081
Journal of Statistical Software: http://www.jst
problem?
>
> Best Regards
> Sherien
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Edzer Pebesma
Institute
t;>> "rgeos_intersection") :
>>>>> TopologyException: no outgoing dirEdge found at 0.5 0
>>>>
>>>> The error message is coming from GEOS - you are welcome to investigate
>>>> further. If you use gUnaryUnion() on the objects first, there is no
>>>> error:
>>>
STATS=1
>> (1 row)
>>
>>
>> gisdb=# select postgis_full_version();
>>
>> postgis_full_vers
>> ion
>> ------
>>
>>
)
> MyLay=SpatialPolygons(list(Pols1,Pols2))
> gIsValid(MyLay)
[1] TRUE
Maybe gUnaryUnion-ing MultiPolygons before intersecting when byid =
FALSE is not such a bad idea after all?
On 25/09/15 21:55, Edzer Pebesma wrote:
> Thanks Vijay; FYI, following OP's example, I see in PostGIS:
&
iguous in terms of the DE-9IM relation with the polygon:
they are both inside the polygon as well as on the polygon boundary.
Still need to sort out why rgeos reports them as valid.
>
> Please pardon the noise.
>
> On Fri, Sep 25, 2015 at 3:55 PM, Edzer Pebesma
> wrote:
>
ssee only. If you are not the intended recipient, you are hereby
> notified that any disclosure, copying, distribution or action in relation to
> the content of this information is strictly prohibited. If you have received
> this e-mail by mistake, please delete the message and any attach
st(pos))
>
> plot(spos)
>
> # Now I want to select (or better get the attribute of) every polygon in "sr"
> that is inside or touches the polygon "spos"?
>
sr does not have attributes.
> I hope this is possible to do with "sp", in orde
;t close in the way I want it. Polygons
> crossing lon +-180 close by crossing the entire globe and not crossing
> the line lon +-180...Is there a way to control this?
I recall there has been discussion about the issue on this list before.
Only 15K messages to go through, according to
http://r-
_____
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
--
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster
Weseler Straße 253, 48151 Münster, Germany. Phone: +
i/Universal_Transverse_Mercator_coordinate_system
might point you to the appropriate zone.
With best regards,
>
>
> Thank you very much for your help!
>
> Kind regards,
>
> Marguerit David
>
> PhD in economics,
>
> University Paris Dauphine, Franc
and somehow to supply
> these to corSpher.
>
> Very grateful for any suggestions or comments.
>
> Dan Bebber
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
--
E
If in the
body of bubble a function is called that you don't understand, read its
documentation.
>
> Any direction (or any documentation where I could find something to go
> ahead) would be welcome.
>
>
>
>
> ___
&g
ach and thence construct a new spdf based upon amore neighborly order
> than a-z.
> Thanks in advance,Chris
>
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mail
can see, it coerces the grid to polygons first, which may
be very inefficient.
Did it work?
--
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster
Weseler Straße 253, 48151 Münster, Germany. Phone: +49 251
8333081, Fax: +49 251 8339763
ents
> of this email and any reply you make may be disclosed by NERC unless
> it is exempt from release under the Act. Any material supplied to
> NERC may be stored in an electronic records management system.
>
>
> ___
> R-sig-Geo mailing
>
> --
> View this message in context:
> http://r-sig-geo.2731867.n2.nabble.com/extracting-row-of-data-in-a-SpatialGridDataFrame-object-tp7579385.html
> Sent from the R-sig-geo mailing list archive at Nabble.com.
>
> ___
> R-si
----Original Message-
> From: Edzer Pebesma [mailto:edzer.pebe...@uni-muenster.de]
> Sent: 29 May 2012 13:23
> To: Rawlins, Barry G.
> Subject: Re: [R-sig-Geo] Plotting with stplot, mode="xy"
>
> svn co https://svn.52north.org/svn/geostatistics/main/spacetime/
&g
sistance would be gratefully received!
>
> Many thanks,
>
> Iain
>
> --
> Iain Dillingham
> School of Informatics
> City University London
> EC1V 0HB
> http://www.soi.city.ac.uk/~abhm372/
> ___
> R-sig-Geo mailing list
minute intervals)
> from a hillslope responding to incident rainfall. I would be happy for you to
> use these as example data if that is useful to you.
>
> Best wishes, Barry
>
> -Original Message-
> From: Edzer Pebesma [mailto:edzer.pebe...@uni-muenster.de]
> Sen
e Professor
> Department of Computer and Mathematical Sciences
> University of Houston - Downtown
> mailto: hodge...@uhd.edu
>
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-pr
edu>
>
>
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
--
Edzer Pebesma
Institute for Geoinformatics (ifg
501 - 600 of 888 matches
Mail list logo
