Tuesday, July 6, 2010

Jaxa Arctic Ice data and UAH AMSU temperatures

2 comments:

  1. Nick

    Welcome as an Arctic Sea Ice observer. I am finding it incredibly interesting and a real data visualization challenge.

    Here's a snippet of R code to help you display the month names rather than day of year in your trend chart.


    #### To plot day of year axis #############################
    ### Day of year axis setup
    ## Set up basic day of year vectors (mon_names, 1st day of mon, mo name plot pos)
    mon_names <- c("Jan", "Feb", "Mar", "April", "May", "June", "July", "Aug", "Sept", "Oct","Nov","Dec")
    mon_doy <- c(1,32,60,91,121,151,182,213,244,274,305,335)
    mon_pos <- c(16, 46, 75, 106,135, 165, 200, 228, 255, 289, 320, 355)

    ## Example DOY plot using month labels
    plot(DOY, sie, type="l", col = "red", axes=F)
    axis(side = 1, at=mon_doy, labels=F)
    axis(side=1, at= mon_pos, labels=mon_names, tick=F, line=F)
    axis(side=2, at=NULL, labels=T)

    Here's how my SIE chart looks using the month names:

    http://chartsgraphs.files.wordpress.com/2010/06/jaxa_sie_07_10_diff_7_06.png

    Kelly

    ReplyDelete
  2. Thanks, Kelly
    I should have checked your site before embarking on this. I've added it to the blogroll to make sure I do so in the future. I'll experiment with your script.

    ReplyDelete