Thanks for the comments, but it's still not working. The version below
uses vectors instead of a matrix for the attributes, but it's still making
two plots all of the same color. Time1 should be one color, Time2 should
be the other. As should be clear from my attribute names, I am trying to
use the lattice-like paneling capability to represent values at different
times of the same variable.
my.text <-
c("POLYGON ((0 2.69, 4.483057 2.86, 4.483057 2.76, 0 2.59, 0 2.69))",
"POLYGON ((0 2.59, 4.483057 2.76, 4.483057 2.26, 0 2.09, 0 2.59))",
"POLYGON ((0 2.09, 4.483057 2.26, 4.483057 1.26, 0 1.09, 0 2.09))",
"POLYGON ((0 1.09, 4.483057 1.26, 4.483057 -0.74, 0 -0.91, 0 1.09))")
my.poly <- st_as_sfc(my.text)
my.poly2 <- st_sf(Time1 = rep(1,4), Time2 = rep(2,4), geom = my.poly)
plot(my.poly2, breaks = c(0.5, 1.5, 2.5))
--Scott Waichler
------------------------------
*From:* Zivan Karaman <zivan.karaman at gmail.com>
*Sent:* Friday, September 27, 2024 9:37 AM
*To:* Waichler, Scott R <Scott.Waichler at pnnl.gov>
*Cc:* r-sig-geo at r-project.org <r-sig-geo at r-project.org>
*Subject:* Re: [R-sig-Geo] unexpected attribute plotting with sf polygons
Check twice before you click! This email originated from outside PNNL.
You have defined two attributes that are constant across your geometries,
so they are shown in the same color.
The code `mat <- matrix(c(rep(1, 2), rep(2, 2)), byrow=FALSE, ncol=1)`
will define one attribute with two distinct values that will be shown in
two colors.
Hope this helps!
On Fri, Sep 27, 2024 at 4:37?PM Waichler, Scott R via R-sig-Geo <
r-sig-geo at r-project.org> wrote:
Hi, When I plot the "sf" and "data.frame" object my.poly2 below, all of
the polygons are the same color, despite my setting two different values
for the Attribute. What have I done wrong?
my.text <-
c("POLYGON ((0 2.69, 4.483057 2.86, 4.483057 2.76, 0 2.59, 0 2.69))",
"POLYGON ((0 2.59, 4.483057 2.76, 4.483057 2.26, 0 2.09, 0 2.59))",
"POLYGON ((0 2.09, 4.483057 2.26, 4.483057 1.26, 0 1.09, 0 2.09))",
"POLYGON ((0 1.09, 4.483057 1.26, 4.483057 -0.74, 0 -0.91, 0 1.09))")
my.poly <- st_as_sfc(my.text)
mat <- matrix(c(rep(1, 4), rep(2, 4)), byrow=F, ncol=2)
my.poly2 <- st_sf(Attribute = mat, geom = my.poly)
plot(my.poly2, breaks = c(0.5, 1.5, 2.5))
Thank you,
Scott Waichler, PhD
Pacific Northwest National Lab
[[alternative HTML version deleted]]