def getDSSImage(ra,dec,radius=1.0,xsize=800,**kwargs):
    """
    Download Digitized Sky Survey images

    https://archive.stsci.edu/cgi-bin/dss_form    
    https://archive.stsci.edu/cgi-bin/dss_search

    Image is in celestial orientation (RA increases to the right)
    https://archive.stsci.edu/dss/script_usage.html

    ra (r) - right ascension
    dec (d) - declination
    equinox (e) - equinox (B1950 or J2000; default: J2000)
    height (h) - height of image (arcminutes; default: 15.0)
    width (w) - width of image (arcminutes; default: 15.0)
    format (f) - image format (FITS or GIF; default: FITS)
    compression (c) - compression (UNIX, GZIP, or NONE; default: NONE; compression 
         applies to FITS only)
    version (v)  - Which version of the survey to use:
         1  - First Generation survey (garden variety)
         2  - Second generation survey (incomplete)
         3  - Check the 2nd generation; if no image is available,
              then go to the 1st generation.
         4  - The Quick V survey (whence came the Guide Stars Catalog;
              used mostly for Phase II proposal submission)
    save (s)  - Save the file to disk instead of trying to display.
         (ON (or anything) or not defined; default: not defined.)

    For the skyview service, see:
    https://skyview.gsfc.nasa.gov/current/docs/batchpage.html
    """
    import subprocess
    import tempfile
    service = 'skyview'

    if service == 'stsci':
        url="https://archive.stsci.edu/cgi-bin/dss_search?"
        scale = 2.0 * radius * 60.
        params=dict(ra='%.3f'%ra,dec='%.3f'%dec,width=scale,height=scale,
                    format='gif',version=1)
        #v='poss2ukstu_red'
    elif service == 'skyview':
        url="https://skyview.gsfc.nasa.gov/cgi-bin/images?"
        params=dict(survey='DSS',position='%.3f,%.3f'%(ra,dec),scaling='Linear',
                    Return='GIF',size=2*radius,projection='Car',pixels=xsize)
    else:
        raise Exception("Unrecognized service.")

    query='&'.join("%s=%s"%(k,v) for k,v in params.items())
    
    tmp = tempfile.NamedTemporaryFile(suffix='.gif')
    cmd='wget --progress=dot:mega -O %s "%s"'%(tmp.name,url+query)
    subprocess.call(cmd,shell=True)
    im = plt.imread(tmp.name)
    tmp.close()
    if service == 'stsci' and xsize: 
        im = scipy.misc.imresize(im,size=(xsize,xsize))
    return im