Skip to content

Instantly share code, notes, and snippets.

View wafels's full-sized avatar

Jack Ireland wafels

View GitHub Profile
@wafels
wafels / gist:6ca1c463c80a9b27d9c92b6e070411ff
Created February 23, 2022 20:51
A small example to show use of spherical screen and composite maps
#
# Use the helioprojective spherical screen to maintain the off-disk emission.
# Use composite map to overlay the two images.
#
import matplotlib.pyplot as plt
from sunpy.map import Map
from sunpy.coordinates import Helioprojective
eit = Map( "/Users/ireland/Data/jp2/TRACE/2008_11_12__08_36_22_272__SOHO_EIT_EIT_284.jp2")
trace = Map("/Users/ireland/Data/jp2/TRACE/2008_11_15__00_11_06_000__TRACE_TRACE_TRACE_195.jp2")
@wafels
wafels / gist:cf79fba38715f9c8e4e0fb7737b6ddb6
Created September 19, 2019 10:07
Experiments in helioviewer JP2 image reprojection
import warnings
#warnings.filterwarnings('ignore')
import copy
import numpy as np
import matplotlib.pyplot as plt
import astropy.units as u
from astropy.coordinates import SkyCoord
@wafels
wafels / turbo_colormap.py
Created August 21, 2019 13:33 — forked from mikhailov-work/turbo_colormap.py
Turbo Colormap Look-up Table
# Copyright 2019 Google LLC.
# SPDX-License-Identifier: Apache-2.0
# Author: Anton Mikhailov
turbo_colormap_data = [[0.18995,0.07176,0.23217],[0.19483,0.08339,0.26149],[0.19956,0.09498,0.29024],[0.20415,0.10652,0.31844],[0.20860,0.11802,0.34607],[0.21291,0.12947,0.37314],[0.21708,0.14087,0.39964],[0.22111,0.15223,0.42558],[0.22500,0.16354,0.45096],[0.22875,0.17481,0.47578],[0.23236,0.18603,0.50004],[0.23582,0.19720,0.52373],[0.23915,0.20833,0.54686],[0.24234,0.21941,0.56942],[0.24539,0.23044,0.59142],[0.24830,0.24143,0.61286],[0.25107,0.25237,0.63374],[0.25369,0.26327,0.65406],[0.25618,0.27412,0.67381],[0.25853,0.28492,0.69300],[0.26074,0.29568,0.71162],[0.26280,0.30639,0.72968],[0.26473,0.31706,0.74718],[0.26652,0.32768,0.76412],[0.26816,0.33825,0.78050],[0.26967,0.34878,0.79631],[0.27103,0.35926,0.81156],[0.27226,0.36970,0.82624],[0.27334,0.38008,0.84037],[0.27429,0.39043,0.85393],[0.27509,0.40072,0.86692],[0.27576,0.41097,0.87936],[0.27628,0.42118,0.89123],[0.27667,0.43134,0.90254],[0.27691,0.44145,0.913
@wafels
wafels / gist:0e9bbee17110bc1fd98e583004a21032
Last active August 21, 2019 12:20
Reproject JP2 files
import warnings
#warnings.filterwarnings('ignore')
import copy
import numpy as np
import matplotlib.pyplot as plt
import astropy.units as u
from astropy.coordinates import SkyCoord
@wafels
wafels / gist:c04bbea5d3beff3558497cb90f75c647
Created January 31, 2018 16:41
Using EIT HEC information
#
# Example EIT HEC to HGS conversion
#
import datetime
import numpy as np
import astropy.units as u
from astropy.coordinates import CartesianRepresentation, HeliocentricTrueEcliptic
import numpy as np
import sunpy.coordinates
from astropy.coordinates import SkyCoord
import astropy.units as u
import matplotlib.pyplot as plt
import sunpy.map
from sunpy.data.sample import AIA_171_IMAGE
# Number of points in great circle
num = 100
layout title author
post
SunPy Update - March 2016
Stuart Mumford

Hello all,

Once again I am a week late, but here is a quick recap.