[R-sig-Geo] spacetime::stplot: problem overlaying a SpatialLines object over a STIDF points map
Edzer Pebesma
edzer.pebesma at uni-muenster.de
Thu Aug 25 13:56:05 CEST 2011
Thanks; for your case, with current spacetime,
stplot(ciotat.t, panel=function(x,y,...) {
panel.xyplot(x,y,...)
sp.lines(contour)
}
)
might do what you need; I'll look into having sp.layout work the same
way as it does in spplot, for stplot'ing STIDF objects.
Best wishes,
On 08/25/2011 01:28 PM, Mathieu Rajerison wrote:
> Hi Edzer,
>
>
> Thanks for your answer.
>
> stplot(ciotat.t) gives a good result:
> http://www.imagup.com/data/1128936911.html
>
> It represents evolution of urbanization along time, each point being the
> date and location of the first building for the city "La Ciotat".
>
> The object is irregular because each spatial location appears one time. I'd
> like to overlay the city contour of "La Ciotat" to better vizualize it.
> http://data.imagup.com/10/1128936911.png
>
> I've sent on your email my STIDF object ciotat.t and my SpatialLines object
> contour I want to be overlaid, if you have some time to look after it.
>
> Thanks a lot,
>
>
> Mathieu
>
> 2011/8/25 Edzer Pebesma <edzer.pebesma at uni-muenster.de>
>
>> Hi Mathieu,
>>
>> the idea behind stplot is indeed that it should be as easy as spplot,
>> hence its name. This idea is however much easier than practice, as the
>> number of possible (and useful) plots for S/T data might be quite large.
>>
>> I would be helpful if you could provide (or generate) a simple data set
>> of your type (STIDF), explain what it represents, and tell how you'd
>> want it to be plotted. Does stplot(ciotat.t) give what you want, i.e. is
>> passing the sp.layout the only problem?
>>
>> Best regards,
>> --
>> Edzer
>>
>>
>> On 08/25/2011 10:12 AM, Mathieu Rajerison wrote:
>>> Hi,
>>>
>>>
>>> I have a STIDF points object named ciotat.t that I've tried to plot with
>> a
>>> city contour named contour.
>>>
>>> I tried this but it doesn't work:
>>>> stplot(ciotat.t, sp.layout= list("sp.lines", contour))
>>>
>>> Also, this doesn't work:
>>>> stplot(ciotat.t, panel=function(...) {
>>> sp.lines(contour)
>>> panel.pointsplot(...)
>>> })
>>>
>>> I thought it was as easy as with spplot.The documentation tells that
>>> arguments are passed on to spplot in case of stplot.STFDF, and to xyplot
>> in
>>> case of
>>> stplot.STIDF. So, the method may be a little different in my case
>> (STIDF)?
>>>
>>> Any help would be appreciated,
>>>
>>> Mathieu
>>>
>>> [[alternative HTML version deleted]]
>>>
>>> _______________________________________________
>>> R-sig-Geo mailing list
>>> R-sig-Geo at 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: +49 251
>> 8333081, Fax: +49 251 8339763 http://ifgi.uni-muenster.de
>> http://www.52north.org/geostatistics e.pebesma at wwu.de
>>
>> _______________________________________________
>> R-sig-Geo mailing list
>> R-sig-Geo at r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>
>
> [[alternative HTML version deleted]]
>
>
>
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at 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: +49 251
8333081, Fax: +49 251 8339763 http://ifgi.uni-muenster.de
http://www.52north.org/geostatistics e.pebesma at wwu.de
More information about the R-sig-Geo
mailing list