Skip to content

Instantly share code, notes, and snippets.

from pathlib import Path
import shutil
import sys
import dask.array as da
import glymur
from numcodecs import Blosc
from loguru import logger
import stack_to_chunk
"""
Code for generating and working with neuroglancer links.
"""
import ast
import json
from pathlib import Path
import neuroglancer
from neuroglancer import CoordinateSpace
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.

The future of time-series data in sunpy

David Stansby

Introduction

In late 2022 I got a small development grant from NumFocus to scope the future of time-series data in sunpy. The sucessful application can be read on the sunpy wiki - this contains context that I won't repeat here.

The current document will be the key outcome of the small development grant, with a record of what I did, the recommendations I made, and any decisions we came to as a community.

In it's current form, please feel free to leave comments.

#!/usr/bin/env bash
#
# Bash script to sample Intel RAPL inteface and save numbers to a file.
# This script will run until manually interrupted.
#
# The CPU_POWER_OUTPUT_FILE environment variable must be set before
# running this script.
# Make sure CPU_POWER_OUTPUT_FILE variable is defined
set -u
HPLinpack benchmark input file
Innovative Computing Laboratory, University of Tennessee
HPL.out output file name (if any)
6 device out (6=stdout,7=stderr,file)
1 # of problems sizes (N)
4000 Ns
1 # of NBs
256 NBs
0 PMAP process mapping (0=Row-,1=Column-major)
1 # of process grids (P x Q)
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
import sunpy.map
from astropy.wcs import WCS
from reproject import reproject_interp
import matplotlib.pyplot as plt
import astropy.units as u
# Load WISPR map
m = sunpy.map.Map('~/Downloads/psp_L1_wispr_20200125T000229_V1_2302.fits')
# Create a new WCS in a helioprojective (HPLN/HPLT) coordinate system with a
from astropy.coordinates import SkyCoord
import astropy.units as u
from matplotlib.colors import LogNorm
import matplotlib.pyplot as plt
import sunpy.data.sample
import sunpy.map
from skimage.transform import warp_polar