Message-ID: <alpine.LRH.2.00.1105301028120.15211@reclus.nhh.no>
Date: 2011-05-30T08:31:02Z
From: Roger Bivand
Subject: newbie bbox type question
In-Reply-To: <BANLkTik2S+K-nRxbNQS_yfpCMCo_HBUdig@mail.gmail.com>
On Mon, 30 May 2011, Barry Rowlingson wrote:
> On Sun, May 29, 2011 at 8:16 PM, Hodgess, Erin <HodgessE at uhd.edu> wrote:
>
>> I want to have one map (via writeOGR) which has the whole thing, and then have several maps by days.
>>
>> However, I want the whole area as defined by the entire data set to show every day.
>>
>> I looked at bbox but I can't set a bbox.
>
> I think the problem here is that writeOGR doesn't create maps, it
> creates data files.
>
> If you wnt to create a map plot that is based on the bbox of a given
> sp object, you can set up an empty plot from its bbox this way:
>
> > plot(SpatialPoints(t(bbox(m))),cex=0)
Equivalently:
plot(as(m, "Spatial"))
which drops all the unnecessary stuff. I wasn't sure what was needed,
either. It could be an output file with the data reshaped to "wide" from
"long", but this is just another guess.
Roger
>
> two things to note: I use cex=0 so the points dont show up (because
> type="n" doesnt work) and the bbox is transposed because of the way
> its rows and columns are arranged.
>
> Then you can add stuff to the plot area, eg the first polygon:
>
> > plot(m[1,],add=TRUE)
>
> Not sure if this is what you are really after since its not clear what
> your 'map' is and where writeOGR fits into this...
>
> Barry
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Roger Bivand
Economic Geography Section, Department of Economics, Norwegian School of
Economics and Business Administration, Helleveien 30, N-5045 Bergen,
Norway. voice: +47 55 95 93 55; fax +47 55 95 95 43
e-mail: Roger.Bivand at nhh.no