[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