Hi? subset(spGRID, raster==1) or subset(spGRID, data$raster==1) Can I do this or something similar?? i.e I'd like to avoid to convert raster to regular data.frame and use subset on dataf.rame due to long time of conversion tahnx Jarek