Cadair

# -*- coding: utf-8 -*-
"""
:Created on: Tue Aug  7 16:44:50 2012

:author: Stuart Mumford

This is a wrapper for a general rotate method, it will rotate and rescale the
image around an arbitary point and then if requested shift the point to the centre.

"""

Cadair

def fieldlines(self,v1,v2,seeds,dS=1):
    """ Simple Euler fieldline integrator
    v1,v2 - y and x vectors
    seeds - array of coordinate (array indexes) pairs
    dS [optional] - step size
    """
    from scipy import ndimage
    
    field = []
    #Get extent of domain

Cadair

read_sdo, "aia.lev1.193A_2012-06-10T13:00:07.84Z.image_lev1.fits",header,data

rotated_cubic = rot(data,90,cubic=-0.5)

save, data, rotated_cubic

Cadair

# -*- coding: utf-8 -*-
"""
:Created on: Fri Sep 14 15:03:51 2012

:author: Stuart Mumford

Test Sunpy.map.rotate() against idl rot()
"""

import numpy as np

Cadair

# -*- coding: utf-8 -*-
"""
:Created on: Fri Sep 14 15:03:51 2012

:author: Stuart Mumford

Test scipy's spline rotation against sunpy's bicubic rotation.
"""

import numpy as np

Cadair

# -*- coding: utf-8 -*-
"""
Created on Tue Nov 13 15:05:36 2012

@author: stuart

Problem from SunPy mailing list. Save a MASSIVE array to a fits file
"""
import numpy as np
import numpy.lib as nplib

Cadair

Report
======
6627 statements analysed.

Messages by category
--------------------

+-----------+-------+---------+-----------+
|type       |number |previous |difference |
+===========+=======+=========+===========+

Cadair

************* Module sunpy.sun.position
I0011:  7,0: Locally disabling E1101
I0011:  7,0: Locally disabling E1121
R0914: 26,0:position: Too many local variables (22/15)
************* Module sunpy.sun.constants
I0011: 93,0: Locally disabling E0611
C0301:  4,0: Line too long (88/80)
C0301:  5,0: Line too long (83/80)
C0301: 12,0: Line too long (86/80)
C0301: 18,0: Line too long (91/80)

Cadair

    def add_map(self, input_, zorder=None, alpha=1, levels=False):
        """Adds a map to the CompositeMap
        
        Parameters
        ----------
        input_ : {sunpy.map, string}
            Map instance or filepath to map to be added
        zorder : int
            The index to use when determining where the map should lie along
            the z-axis; maps with higher z-orders appear above maps with lower

Cadair

# -*- coding: utf-8 -*-
"""
Created on Sat Jun 29 21:31:58 2013

@author: stuart
"""

import numpy as np
from datetime import timedelta