R - arsa ile sıkışmış() -

Question

Bir Avustralya İstatistik Bürosu'ndan elde edilen Avustralya'da uzak bölgeleri gösteren şekil dosyası var bir yuva değere dayalı şekil dosyası poligon Boyama: Aynı URL'de

http://www.abs.gov.au/AUSSTATS/abs@.nsf/DetailsPage/1270.0.55.005July%202011?OpenDocument

PDF "ASGS Uzaklık Yapısı Sürümü 2011 PDF Haritalar" - Bu PDF belgesindeki ilk haritayı yeniden oluşturmaya çalışıyorum.

Ben shape okumak ve data yuvaya renk bilgisini eklemiş: arsa haritasıdır yapmamı için

ra <- readShapeSpatial("RA_2011_AUST", delete_null_obj = TRUE) 
ra@data$COLOUR <- "#FFFFFF" 
ra@data$COLOUR[(as.numeric(as.character(ra@data$RA_CODE11)) %% 10) == 0] <- "#006837" 
ra@data$COLOUR[(as.numeric(as.character(ra@data$RA_CODE11)) %% 10) == 1] <- "#31A354" 
ra@data$COLOUR[(as.numeric(as.character(ra@data$RA_CODE11)) %% 10) == 2] <- "#78C679" 
ra@data$COLOUR[(as.numeric(as.character(ra@data$RA_CODE11)) %% 10) == 3] <- "#C2E699" 
ra@data$COLOUR[(as.numeric(as.character(ra@data$RA_CODE11)) %% 10) == 4] <- "#FFFFCC" 

tek şey bıraktı! Sıkıştığım yer burasıdır ...

ra@polygons, her biri ID veri karnesine sahip olan ve her biri ra@data veri çerçevesine bir dizin olan 35 çokgenlerin listesidir. Tek yapmam gereken, ra@data$COLOUR[ID]'daki rengi bulmak için plot()'a anlatmak. Pek iyi değil. 35 çokgenin ("Çokgenler") her birinin kendi çokgen listesi vardır ("Çokgen" sınıfı); Toplamda 6902 çokgen var!

plot() anlayışım, çokgenlerin çizileceği gibi, aynı renkteki bir vektörü aktarmak zorunda olduğumdur. Bu nedenle, ilgili poligonun renk değerini tutan her elemanın 6902 uzunluğundaki bir vektörünü oluşturmam gerektiğine inanıyorum. Şimdiye kadar nasıl yapıyorum?

Çokgenler sırayla çizildiyse bu yeterince kolay olurdu, ama değiller. 35 Çokgenlerin her biri, bir tamsayı vektörü olan bir plotOrder yuvasına sahiptir, bu nedenle, bu vektörlerin her birindeki değerler tarafından renk vektörünün sıralanması gerekecektir.

Bu noktada her şey biraz karmaşık görünüyor. Buradan tamamen uzak mıyım?

Tavsiyeniz için teşekkürler!

Answer 1

Zaten tüm işi yaptın!

plot(ra, col=ra$COLOUR) 

Ve işte bu kadar:

plot(ra, col=ra@data$COLOUR) 

Hatta, @Spacedman olarak öne!

Ve poligon sınırları kurtulmak istiyorsanız:

plot(ra, col=ra$COLOUR, border=NA) 

Düzenleme:

göre: Bu davranışı açıklamak için bir girişim ?SpatialPointsDataFrame:

SpatialPolygonsData Varsayılan kimlik eşlemeli veri çerçevesi, veri çerçeve satır adlarını Çokgen Kimlik yuvalarına karşı denetler. Daha sonra birbirleriyle hemfikir olmalılar ve benzersiz olmalılar (hiçbir Poligon nesnesi kimlikleri paylaşamaz); Çokgen kimlikleriyle eşleşmesi gerekiyorsa veri çerçeve satırları yeniden sıralanacaktır.çokgenler dolayısıyla onların numarası olarak sipariş ve bunun anlamı

yuva @data içinde dataframe sıraları olarak sıralanmıştır.

... 
    polys <- slot(x, "polygons") 
    pO <- slot(x, "plotOrder") 
    if (!is.null(density)) { 
     if (missing(col)) 
      col <- par("fg") 
     if (length(col) != n) 
      col <- rep(col, n, n) 
     if (length(density) != n) 
      density <- rep(density, n, n) 
     if (length(angle) != n) 
      angle <- rep(angle, n, n) 
     for (j in pO) .polygonRingHoles(polys[[j]], border = border[j], 
      xpd = xpd, density = density[j], angle = angle[j], 
      col = col[j], pbg = pbg, lty = lty, ...) 
    } 
... 

col girilen vektör polygons yuvası ile aynı düzeni sahiptir ve bu nedenle numarası:

Şimdi bir noktada bu hatlar bulunmaktadır (getAnywhere(plot.SpatialPolygons) kullanılarak) işlevi plot.SpatialPolygons bakarsak. Hepsi aynı şekilde dizinlemek için plotOrder kullanılır.