Dear Tomislav and Virgilio,

It is funny,... I wanted to try it, because I was going through your 
reports and manuals :-)

I am not sure, but I think the projection and reference system are ok, 
because when I create the klm it works fine. Only with the png, does not 
function.


summary(grid.bord[1])

Object of class SpatialPixelsDataFrame
Coordinates:
          min       max
long -8.19592 -7.392709
lat  36.94938 37.177010
Is projected: FALSE
proj4string :
[+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0]
Number of points: 2895
Grid attributes:
     cellcentre.offset    cellsize cells.dim
long         -8.193378 0.005083612       158
lat          36.950000 0.001237113       184
Data attributes:
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
  59.73   78.67   98.95   96.99  115.40  130.50


I will have a look at the resources you gave.

Thank you,
BW,
Marta
 





Tomislav Hengl wrote:
> Dear Marta,
>
> I am not sure what is the projection system you use for the Algarve case study? If it is different
> from CRS("+proj=longlat +datum=WGS84"), then you do need to first reproject your maps.
>
> You should really look at the examples from the book by Bivand et al.:
> http://www.asdar-book.org/code.php?chapter=3&figure=1
> http://www.asdar-book.org/code.php?chapter=3&figure=2 
>
> Here you will find small examples of how to export any raster map to GE (it reprojects a map using
> the MapResample method of ILWIS GIS):
> http://geomorphometry.org/R.asp 
>
> I personally do not use much the maptools to generate GE KMLs, but instead write directly a KML file
> from R e.g. to export a time-series of interpolated temperature images to KML in a loop:
>
>   
>> filename = "croclim12.kml"
>> write("<?xml version=\"1.0\" encoding=\"UTF-8\"?>", filename)
>> write("<kml xmlns=\"http://earth.google.com/kml/2.2\">", filename, append = TRUE)
>> write("<Folder>", filename, append = TRUE)
>> write("\t<name>Mean Daily TEMP</name>", filename, append = TRUE)
>> write("\t<open>1</open>", filename, append = TRUE)
>> for (i in slices) {
>>     
>  slice = floor(unclass(as.POSIXct(i))/86400)[[1]]
>  write("\t\t<GroundOverlay>", filename, append = TRUE)
>  write(paste("\t\t<name>", i, "</name>", sep = ""),
>  filename, append = TRUE)
>  write("\t\t<TimeSpan>", filename, append = TRUE)
>  write(paste("\t\t\t<begin>", slice, "</begin>", sep = ""),
>  filename, append = TRUE)
>  write(paste("\t\t\t<end>", slice + 1, "</end>", sep = ""),
>  filename, append = TRUE)
>  write("\t\t</TimeSpan>", filename, append = TRUE)
>  write("\t\t<color>99ffffff</color>", filename, append = TRUE)
>  write("\t\t<Icon>", filename, append = TRUE)
>  write(paste("\t\t\t<href>", getwd(), "/MDTEMP", slice, ".png</href>", sep = ""), filename, append =
> TRUE)
>  write("\t\t\t<viewBoundScale>0.75</viewBoundScale>",
>  filename, append = TRUE)
>  write("\t\t</Icon>", filename, append = TRUE)
>  write("\t\t<altitude>50</altitude>", filename, append = TRUE)
>  write("\t\t<altitudeMode>relativeToGround</altitudeMode>",
>  filename, append = TRUE)
>  write("\t\t<LatLonBox>", filename, append = TRUE)
>  write(paste("\t\t\t<north>", grids.kml$ylim[[2]],
>  "</north>", sep = ""), filename, append = TRUE)
>  write(paste("\t\t\t<south>", grids.kml$ylim[[1]],
>  "</south>", sep = ""), filename, append = TRUE)
>  write(paste("\t\t\t<east>", grids.kml$xlim[[2]],
>  "</east>", sep = ""), filename, append = TRUE)
>  write(paste("\t\t\t<west>", grids.kml$xlim[[1]],
>  "</west>", sep = ""), filename, append = TRUE)
>  write("\t\t</LatLonBox>", filename, append = TRUE)
>  write("\t</GroundOverlay>", filename, append = TRUE)
> }
>   
>> write("<ScreenOverlay>", filename, append = TRUE)
>> write(paste(" <name>values from", MDTEMPxlim[[1]], " to ", MDTEMPxlim[[2]], "</name>", sep = ""),
>>     
> filename, append = TRUE)
>   
>> write(" <Icon>", filename, append = TRUE)
>> write(paste("\t\t\t<href>", getwd(), "/legend.png</href>", sep = ""), filename, append = TRUE)
>> write(" </Icon>", filename, append = TRUE)
>> write(" <overlayXY x=\"0\" y=\"1\" xunits=\"fraction\" yunits=\"fraction\"/>", filename, append =
>>     
> TRUE)
>   
>> write(" <screenXY x=\"0\" y=\"1\" xunits=\"fraction\" yunits=\"fraction\"/>", filename, append =
>>     
> TRUE)
>   
>> write(" <rotationXY x=\"0\" y=\"0\" xunits=\"fraction\" yunits=\"fraction\"/>", filename, append =
>>     
> TRUE)
>   
>> write(" <size x=\"0\" y=\"0\" xunits=\"fraction\" yunits=\"fraction\"/>", filename, append = TRUE)
>> write(" </ScreenOverlay>", filename, append = TRUE)
>> write("</Folder>", filename, append = TRUE)
>> write("</kml>", filename, append = TRUE)
>>     
>
>
> I guess you need to read the KML tutorial
> (http://code.google.com/apis/kml/documentation/kml_tut.html) to find out what exactly you need to
> use.
>  
> hth
>
> Tom Hengl 
> http://spatial-analyst.net 
>
>
> -----Original Message-----
> From: r-sig-geo-bounces@stat.math.ethz.ch [mailto:r-sig-geo-bounces@stat.math.ethz.ch] On Behalf Of
> Marta Rufino
> Sent: Friday, October 24, 2008 11:50 AM
> To: r-sig-geo@stat.math.ethz.ch
> Subject: [R-sig-Geo] bad position of png & klm file in google earth
>
> Dear list,
>
> I would like to start by acknowledging the wonderful work of this 
> support list :-)
> I have  a question, which should be prety dummy,... but I was looking 
> around in the net, in reports, etc. and could not sort it out.
>
> I am trying to export a image plot to google earth. I manage to do it 
> with  points data... ok. but when I do it with the images, they simply 
> go to Fiji (which I am sure it is a nice place... but a bit far from the 
> Algarve, where it should be).
> Even the example given in the help file, the image goes to the same 
> place. What am I doing wrong?
>
>
> require(maptools); require(rgdal)
> # This is the example, in the file
> opt_exask <- options(example.ask=FALSE)
> qk <- SpatialPointsDataFrame(quakes[, c(2:1)], quakes)
> proj4string(qk) <- CRS("+proj=longlat")
> tf <- tempfile()
> SGqk <- GE_SpatialGrid(qk)
> png(file=paste(tf, ".png", sep=""), width=SGqk$width, 
> height=SGqk$height, bg="transparent")
> par(mar=c(0,0,0,0), xaxs="i", yaxs="i")
> plot(qk, xlim=SGqk$xlim, ylim=SGqk$ylim, setParUsrBB=TRUE)
> dev.off()
> kmlOverlay(SGqk, paste(tf, ".kml", sep=""), paste(tf, ".png", sep=""))
>
> # than the files produced are in
> tf()
>
> I go to the file, pick it the two of them and drag it in google earth... 
> and there they are, in Fiji!!! I am sure it is silly, but it is my first 
> steps in googleearth--R :-(
> If I open it through google earth, it is the same...
>
> A second doubt I have is if we can set apriori, in R the markers 
> projected in google earth... how do I do if instead of the pins, I want 
> a ball? Do I have to change one by one?
>
>
> Thank you very much in advance,
> Best wishes,
> Marta
>
> I am in windows, R 2.7.1, maptools, version 0.7-15, 2008-08-14...
> GDAL 1.5.2, released 2008/05/29
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo@stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
>
>
>   

-- 
.......................................................................
Marta M. Rufino (PhD)

.....
Instituto Nacional de Investigação Agrária e das Pescas (INIAP/IPIMAR),
Centro Regional de Investigação Pesqueira do Sul (CRIPSul)
Avenida 5 de Outubro s/n
P-8700-305 Olhão, Portugal
+351 289 700 541

..... 
Institut de Ciències del Mar - CMIMA (CSIC)
Passeig Marítim de la Barceloneta, 37-49      
08003 BARCELONA - Catalunya
Spain 






	[[alternative HTML version deleted]]


