numberformat.py

# Formats a number tag in the OEC XML format.
# <tag errorplus="0.1" errorminus="0.2">1.0</tag>
from math import *

def toFixed(f,digits):
    formatstring = "%%.%df"%digits
    return formatstring % f

notAvailableString = "<span class=\"NA\"><span class=\"numericvalue\"></span>N/A</span>"

def renderText(tag):
    if tag is None:
        return notAvailableString
    return tag.text

def renderFloat(tag, factor=1.):
    if tag is None:
        return notAvailableString
    if "errorplus" in tag.attrib:
        error_plus = float(tag.attrib["errorplus"])*factor
    else:
        error_plus = None
    if "errorminus" in tag.attrib:
        error_minus = float(tag.attrib["errorminus"])*factor
    else:
        error_minus = None
    if "upperlimit" in tag.attrib:
        upperlimit = float(tag.attrib["upperlimit"])*factor
    else:
        upperlimit = None
    if "lowerlimit" in tag.attrib:
        lowerlimit = float(tag.attrib["lowerlimit"])*factor
    else:
        lowerlimit = None
    if tag.text is not None:
        value = float(tag.text)*factor
    else:
        value = None

    hideErrors = False
    dataString = ""
    exp0String = ""
    exp1String = ""
    exp2String = ""
    errorMinusString = ""
    errorPlusString = ""
    errorString = ""
    
    isExp = True
    exponent = ""
    dexponent = ""
    if value is not None:
        if value!=0.:
            exponent = floor(log10(fabs(value)))
            dexponent = pow(10, exponent)
        else:
            exponent = -2
            dexponent = pow(10, exponent)

    elif upperlimit is not None and upperlimit!=0.0:
        exponent = floor(log10(upperlimit))
        dexponent = pow(10, exponent)
    elif lowerlimit is not None and lowerlimit!=0.0:
        exponent = floor(log10(lowerlimit))
        dexponent = pow(10, exponent)
    else:
        return notAvailableString
            
    if (exponent>-4 and exponent<5) or (value is not None and value<=0.):
        isExp = False
    else:
        diexponent = 1./dexponent
        if value is not None:
            value		*=diexponent
        if error_minus is not None:
            error_minus	*=diexponent
        if error_plus is not None:
            error_plus	*=diexponent
        if upperlimit is not None:
            upperlimit	*=diexponent
        if lowerlimit is not None:
            lowerlimit	*=diexponent
    
    significantdigits = 2
    if value is not None:
        if value!=0.:
            significantdigits = -floor(log10(fabs(value/2.)))
        else:
            significantdigits = 0
        if error_plus is not None and error_minus is not None and error_plus!=0. and error_minus!=0.:
            error_significantdigits = max(-floor(log10(fabs(error_plus/2.))),-floor(log10(fabs(error_minus/2.))))
            significantdigits = max(error_significantdigits,significantdigits)
        else:
            significantdigits +=2
        
    else:
        # upper/lower limits
        significantdigits = 4
    
    significantdigits = min(significantdigits,4.)
    significantdigits = max(significantdigits,0.)

    hasErrors = False
    if (error_plus is not None or error_minus is not None) and (not hideErrors):
        hasErrors = True
    
    
    if value is not None:
        dataString = "<span class=\"numericvalue\">"+toFixed(value,significantdigits)+"</span>"
    else:
        if lowerlimit is not None:
            dataString = "&gt; <span class=\"numericvalue\">"+toFixed(lowerlimit,significantdigits)+"</span>"
        elif upperlimit is not None:
            dataString = "&lt; <span class=\"numericvalue\">"+toFixed(upperlimit,significantdigits)+"</span>"
        
    
    if isExp:
        exp2String = "%d"%int(exponent)
        if hasErrors:
            exp1String = " )&#183;10"
        else:
            exp1String = " &#183;10"
        
    else:
        exp1String = ""
        exp2String = ""
    
    
    if hasErrors:
    	if error_minus is not None:
    		errorMinusString = "&#8722;"+toFixed(error_minus,significantdigits)
    	else:
    		errorMinusString = ""
    	
    	if error_plus is not None:
    		errorPlusString = "&#43;"+toFixed(error_plus,significantdigits)
    	else:
    		errorPlusString = ""
    	
    	if error_minus is not None and error_plus is not None:
            if errorMinusString[7:]==errorPlusString[5:]:
                errorString = "&#177;"+ errorMinusString[7:]
                errorMinusString = ""
                errorPlusString = ""
            else:
                errorString = ""
    	else:
            errorString = ""
    
    else:
        errorString		= ""
        errorMinusString	= ""
        errorPlusString		= ""
    
    
    if isExp and hasErrors:
        exp0String = "( "
    else:
        exp0String = ""
    
    returnString = ""
    returnString += exp0String
    returnString += dataString
    returnString += errorString
    if errorPlusString or errorMinusString:
        returnString += "<span class=\"errorbars\">"
        returnString += "<span class=\"errorplus\">"+errorPlusString+"</span>"
        returnString += "<span class=\"errorminus\">"+errorMinusString+"</span>"
        returnString += "</span>"
    
    if exp2String or exp1String:
        returnString += exp1String+"<sup>"+exp2String+"</sup>"
    
    
    return returnString;

def getFloat(obj,tag,default=None):
    v = obj.find(tag)
    if v is not None:
        try:
            return float(v.text)
        except:
            return default
    return default

def getText(obj,tag,default=None):
    v = obj.find(tag)
    if v is not None:
        return v.text
    return default 

def getCoordinates(obj):
    ra = obj.find("./rightascension")
    dec = obj.find("./declination")
    if ra is None or dec is None:
        return None
    ra = ra.text.split(" ")
    dec = dec.text.split(" ")
    raf = float(ra[0])
    decf = float(dec[0])
    sign = 1.
    if decf<0.:
        sign = -1.
    raf += float(ra[1])/24.
    decf += sign*float(dec[1])/60.
    raf += float(ra[2])/24./60.
    decf += sign*float(dec[2])/60./60.
    return raf/24.*2.*pi,decf/360.*2.*pi