R/set.fitskey.R

In leeskelvin/astro: Astronomy Functions, Tools and Routines for R

set.fitskey = function(key, value, hdr){
    
    # create key/value vectors
    k = rep(NA, max(length(key),length(value)))
    v = rep(NA, max(length(key),length(value)))
    k[1:length(key)] = key
    v[1:length(value)] = value
    
    # update pre-existing keys
    if(any(k %in% hdr[,"key"])){
        
        # determine order of vectors
        inhdr = which(hdr[,"key"] %in% k)
        inkey = which(k %in% hdr[,"key"])
        matched = inhdr[match(k[inkey],hdr[inhdr,"key"])]
        
        # update header
        hdr[matched,"value"] = v[inkey]
        
    }
    
    # create new keys
    if(any(!k %in% hdr[,"key"])){
        
        # determine new keys
        outhdr = which(!k %in% hdr[,"key"])
        
        # append to header
        hdr = rbind(hdr, cbind(key=k[outhdr], value=v[outhdr], comment=""))
        
    }
    
    # return results
    return(hdr)
    
}