[R-sig-Geo] Combing attribute table with polygon??

[Robert J. Hijmans]

Zia, something like this might work as long as MUKEY has no NA values:

poly.data = merge(poly.data, data.table, by="MUKEY" ,all.x = TRUE)

Robert

On Wed, Jan 27, 2010 at 1:58 PM, Zia Ahmed <zua3 at cornell.edu> wrote:
> I have created one data table and one polygon in R. Both have a similar type
> of column.  I want to combine or join them with this column or field
> ("MUKEY").  Then  I will convert this polygon  as  a ESRI  shape file.
> But  when  I  used  spCbind,   I got error. I think i miss something here.
> Is there any way join attribute table with polygon? Thanks
> Zia
>
> poly.data1 <- spCbind(poly.data, data.table)
> Error in spCbind(poly.data, data.table) : different numbers of rows
>
>
> names(poly.data) # polygoan
> [1] "AREASYMBOL" "SPATIALVER" "MUSYM"      "MUKEY"    > names(data.table) #
> data
> [1] "MUKEY" "Sand"  "Silt"  "Clay"  "pH"    "CEC"   "EC"    "SAR"   "CaCO3"
> [10] "SOM"   "BD"    "AWC"   "KSAT"  "Kf"    "Kw"    "LL"    "LEP"
>> length(poly.data$MUKEY)
> [1] 76969
>> length(data.table$MUKEY)
> [1] 184
>>
>
>
> poly.data1 <- spCbind(poly.data, data.table)
>
> combinedShp1 <- "D:/test/combined1.shp" writeOGR(poly.data1,
> dsn=combinedShp1, layer="combined1", driver="ESRI Shapefile")
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>