dpshelio

import datetime
import calendar
from itertools import product
import json

from sunpy.net import vso

def add_months(sourcedate, months):
    month = sourcedate.month - 1 + months
    year = int(sourcedate.year + month / 12 )

dpshelio

git clone git@github.com:helio-vo/helio.git helio-all
cd helio-all
for i in $(find ./ -not -name '.*' -type d -maxdepth 1); do cp -r ../helio-all \.${i}; cd \.${i}; git filter-branch --prune-empty --subdirectory-filter  ${i} -- --all; cd ../helio-all; done
cd ..
for i in $(find ./ -type d -maxdepth 1'); do cd ${i//.\/}; git remote set-url origin git@github.com:helio-vo/${i//.\/}; git push --all origin; cd .. ; done

dpshelio

dpshelio

dpshelio

dpshelio

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

@author: stuart
"""

import numpy as np
from datetime import timedelta

dpshelio

import astropy.units as u
import matplotlib.pyplot as plt
import numpy as np
from skimage.measure import find_contours

import sunpy.map

xsub = 400 * u.arcsec + 75 * ([-1, 1] * u.arcsec)
ysub = 300 * u.arcsec + 50 * ([-1, 1] * u.arcsec)
map01 = sunpy.map.Map('./hmi.m_45s.2010.08.10_12_01_30_TAI.magnetogram.fits').submap(xsub, ysub)

dpshelio

dpshelio

device, true_color = 24
device, decomposed = 0
device, retain = 2

dpshelio

import numpy as np
import matplotlib.pyplot as plt

import astropy.units as u
from astropy.wcs import WCS
from astropy.io import fits
from wcsaxes import WCSAxes
import sunpy.map

from sunpy.net import vso