[R-sig-Geo] Export subsets of Spatial Points Data Frame to images

[Edzer Pebesma]

On 01/05/2015 04:43 PM, Ariel Fuentesdi wrote:
> Frede,
> 
> This is the reproducible example:
> 
> date <- seq(as.Date("2015/1/1"), as.Date("2015/1/5"), "days")
> stops <- c("aa", "bb", "cc", "dd", "ee")
> F_Up <-  c(1:10)[-seq(from = 4, to = 8, by = 1)]
> x <- c(342328.0, 342374.5, 342413.2, 342348.6, 342394.6)
> y <-  c(6308284, 6308123, 6307915, 6307903, 6307915)
> df <- data.frame(date, stops, F_Up, x, y)
> 
> 
> coordinates(df) <- ~ x + y
> proj4string(obj = df) <- CRS("+init=epsg:32719")
> 
> df2 <- split(df, f = df$date)
> 
> for (i in df2) {

in this loop, i is not an index, but is one of the df2 list elements.

>   jpeg(paste(df2[[1]],".jpg",sep=""))

this does not create valid file names.

>   bubble(df[i], "F_Up", col=c("#00ff0088", "#00ff0088"),
>          main = "df2" )

this does have no effect, because in a loop, bubble needs to be print()ed.

>   dev.off()
> }

The following does something, and may help you further:

for (i in seq(along = df2)) {
  jpeg(paste(names(df2)[i], ".jpg", sep=""))
  print(bubble(df2[[i]], "F_Up", col=c("#00ff0088", "#00ff0088"),
         main = "df2" ))
  dev.off()
}



> 
> 
> 2015-01-05 3:22 GMT-04:00 Frede Aakmann Tøgersen <frtog at vestas.com>:
> 
>> Hi Ariel
>>
>> Please provide a reproducible example. Now I have to guess what the SB
>> object is because you didn't provide any details about that.
>>
>> I guess that the group_by() and summarise() are from the dplyr package.
>>
>> I guess that when you have imported the csv file you have an object SB
>> containing columns Fecha, CodPar, x, and y.
>>
>> If so, your calls to the two dplyr functions are not correct. As first
>> argument they need a data object, which is supposed to be of a tbl type
>> like e.g. data frame.
>>
>> Why didn't you get an error when issuing e.g. group_by(Fecha, CodPar, x,
>> y) instead of group_by(SB, Fecha, CodPar, x, y)?
>>
>> Yours sincerely / Med venlig hilsen
>>
>>
>> Frede Aakmann Tøgersen
>> Specialist, M.Sc., Ph.D.
>> Plant Performance & Modeling
>>
>> Technology & Service Solutions
>> T +45 9730 5135
>> M +45 2547 6050
>> frtog at vestas.com
>> http://www.vestas.com
>>
>> Company reg. name: Vestas Wind Systems A/S
>> This e-mail is subject to our e-mail disclaimer statement.
>> Please refer to www.vestas.com/legal/notice
>> If you have received this e-mail in error please contact the sender.
>>
>>> -----Original Message-----
>>> From: R-sig-Geo [mailto:r-sig-geo-bounces at r-project.org] On Behalf Of
>> Ariel
>>> Fuentesdi
>>> Sent: 5. januar 2015 04:48
>>> To: r-sig-geo at r-project.org
>>> Subject: [R-sig-Geo] Export subsets of Spatial Points Data Frame to
>> images
>>>
>>> I have a csv file, which I converted it to Spatial Points Data Frame.
>> And I
>>> want to split the spatial data and export each map to a jpg file.
>>>
>>>
>>> I did the following:
>>>
>>>     SB <- read.csv("SB_040914.csv", header=TRUE, sep = ";")
>>>     SB <- SB %>% group_by(Fecha, CodPar, x, y) %>%
>>> summarise(mean(Sube),mean(Evade), mean(Baja))
>>>     SB <- data.frame(SB)
>>>     SB <- SB[complete.cases(SB),]
>>>
>>>     coordinates(SB) <- ~ x + y
>>>     proj4string(obj = SB) <- CRS("+init=epsg:32719")
>>>
>>>     SB <- split(SB, f = SB$Fecha)
>>>
>>> The Split function, shows me the next warnings:
>>>
>>> **Warning messages:
>>> 1: In min(x) : no non-missing arguments to min; returning Inf
>>> 2: In max(x) : no non-missing arguments to max; returning -Inf
>>> 3: In min(x) : no non-missing arguments to min; returning Inf
>>> 4: In max(x) : no non-missing arguments to max; returning -Inf**
>>>
>>> And for trying to export I did this:
>>>
>>>     for (i in SB) {
>>>       jpeg(paste(SB[[1]],".jpg",sep=""))
>>>       bubble(SB[[i]], "mean.Evade.", col=c("#00ff0088", "#00ff0088"),
>>>              main = "Evasión" )
>>>       dev.off()
>>>     }
>>>
>>> **Which shows me the next error:
>>> Error in as.character.default(<S4 object of class
>>> "SpatialPointsDataFrame">) :
>>>   no method for coercing this S4 class to a vector
>>> 6 as.character.default(<S4 object of class
>>> structure("SpatialPointsDataFrame", package = "sp")>)
>>> 5 as.character(<S4 object of class structure("SpatialPointsDataFrame",
>>> package = "sp")>)
>>> 4 paste(SB[[1]], ".jpg", sep = "")
>>> 3 gsub("%%", "", s)
>>> 2 checkIntFormat(filename)
>>> 1 jpeg(paste(SB[[1]], ".jpg", sep = ""))**
>>>
>>>       [[alternative HTML version deleted]]
>>>
>>> _______________________________________________
>>> R-sig-Geo mailing list
>>> R-sig-Geo at r-project.org
>>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>
> 
> 	[[alternative HTML version deleted]]
> 
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
> 

-- 
Edzer Pebesma,  Co-Editor-in-Chief Computers & Geosciences
Institute for Geoinformatics (ifgi), University of Münster
Heisenbergstraße 2, 48149 Münster, Germany. Phone: +49 251
83 33081 http://ifgi.uni-muenster.de GPG key ID 0xAC227795

-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 490 bytes
Desc: OpenPGP digital signature
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20150105/6561400a/attachment.bin>