Following on that I wish to download several single tiles. To do it I declare the tile extension I wish to download and adapting Robert's hint I do the following library(dismo) #I need to load XML as well library(XML) library(maptools) library(rgdal) #declare the extent e = extent(36.1195,36.1205,5.8725,5.8735) #got the map, exp less than 1.0 zooms the result r <- gmap(e, type="satellite", filename="mypath/mynumbered.gif", exp=1.0) #check the result plot(r)
The thing is that now I want to download several tiles getting the extent from a table that looks like that xmin xmax ymin ymax num 1 36.09940 36.10340 5.671173 5.675173 10 2 36.19274 36.19674 5.665194 5.669194 15 3 36.17748 36.18148 5.820774 5.824774 19 4 36.19274 36.19674 5.665194 5.669194 43 5 36.09940 36.10340 5.671173 5.675173 44 6 36.17748 36.18148 5.820774 5.824774 46 so I want pass values from the table to 'extent' and get and numbered gif for each. Any ideas? Cheers Agustin sessionInfo R version 2.11.1 Patched (2010-08-26 r52824) Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit) locale: [1] en_US/en_US/en_US/C/en_US/en_US attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] maps_2.1-4 gpclib_1.5-1 maptools_0.7-38 lattice_0.18-8 foreign_0.8-40 XML_3.2-0 dismo_0.5-11 rJava_0.8-8 raster_1.7-29 [10] rgdal_0.6-28 sp_0.9-76 loaded via a namespace (and not attached): [1] grid_2.11.1 tools_2.11.1 On Feb 3, 2011, at 9:12 PM, Robert Hijmans wrote: > >> Dear all, >> I'm using the package RgoogleMaps to import google earth maps into R. The >> object that is retrieved with function GetMap contains a number of slots >> which I think can be used to translate the google picture into a geoTiff >> to for example import into ArgGis. >> If anyone would have an example of how to accomplish that or what function >> to use, I'd be very greatful. >> Best wishes, >> Erik Meesters >> > > Here is how you can get a single tile > > library(dismo) > library(maptools) > library(rgdal) > > # get a google map > g <- gmap("Australia") > plot(g) > > # overlay boundaries, after projecting to Mercator > data(wrld_simpl) > w = wrld_simpl[wrld_simpl@data[,"NAME"] != "Antarctica", ] > m = spTransform(w, projection(g, asText=F)) > plot(m, add=T) > > # export to png > filename <- 'oz.png' > png(filename, width=ncol(g), height=nrow(g)) > par(mar=c(0,0,0,0)) > image(g, col=g@legend@colortable) > dev.off() > > # create 'world file' for georferencing png > r <- raster(filename) > extent(r) <- extent(g) > hdr(r, 'worldfile', 'pgw') > > # now open the file in ArcMap > > -- > View this message in context: > http://r-sig-geo.2731867.n2.nabble.com/RgoogleMaps-tp5989261p5990287.html > Sent from the R-sig-geo mailing list archive at Nabble.com. > > _______________________________________________ > R-sig-Geo mailing list > R-sig-Geo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo