Sunday, July 1, 2012

Colored 3D Map

In my previous post, I showed how to make a 3D view of an area using the persp function. However, I felt this was not a complete representation, especially for digital elevation. While looking for some reference for my presentation on the use of R for data analysis and visualization, I came across a thread discussing " 3d plot with different levels done in different colors (solution)".  In this thread, a co-blogger Tamas K Papp's shared his levelpersp function. Using this function I tried to plot 3D view of Marinduque, Philippines using terrain and topo color.


3D persp using terrain color

3D persp using topo color



5 comments:

  1. It would be great if you could link to the thread you found, and provide more explanation about what you did to generate these plots.

    ReplyDelete
    Replies
    1. Hi Brandon,

      Please find below:
      https://stat.ethz.ch/pipermail/r-help/2003-July/036151.html

      Delete
  2. sorry but what is the use of a plot without its code and data, is R a drawing competition

    ReplyDelete
  3. This comment has been removed by the author.

    ReplyDelete
  4. Hi Arnold! I really liked your plot and trying to do the same for an elevation map of Morocco.
    I am using the following code to get a persp plot ok:

    alt <- getData("alt", country="MAR")
    adm <- getData("GADM", country="MAR", level=1)
    #mar<-(adm[adm$NAME_1==”Marinduque”,])
    maralt<-crop(alt,adm)
    persp(maralt, exp=0.2,phi=35, xlab="Longitude", ylab="Latitude", zlab="Elevation")


    But I can not figure out how to translate this to use the levelpersp function. How do I assign the x, y & z values?

    Maybe you could share how you did that?

    ReplyDelete