Skip to content

Instantly share code, notes, and snippets.

@wtbarnes

wtbarnes/eui-ar-movie.ipynb

Created January 11, 2024 20:40
Show Gist options
  • Save wtbarnes/4620d83e2804896cd0a219f602a90dbe to your computer and use it in GitHub Desktop.
Save wtbarnes/4620d83e2804896cd0a219f602a90dbe to your computer and use it in GitHub Desktop.
Download ZIP
Making movies of ARs as observed by EUI on Solar Orbiter
Display the source blob
Display the rendered blob
Raw
eui-ar-movie.ipynb
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Display the source blob
Display the rendered blob
Raw
eui-ar-tracking.ipynb
{
"cells": [
{
"cell_type": "markdown",
"id": "d533762b-f88d-47af-b4ee-1bdffc46d105",
"metadata": {},
"source": [
"# EUI AR Tracking"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "1c4a7a64-d188-4eb7-b595-95801757e5f3",
"metadata": {},
"outputs": [],
"source": [
"from sunpy.net import Fido, attrs as a\n",
"import sunpy.map\n",
"import astropy.units as u\n",
"import sunpy_soar\n",
"import matplotlib.pyplot as plt\n",
"from astropy.coordinates import SkyCoord\n",
"from sunpy.coordinates import propagate_with_solar_surface\n",
"from astropy.visualization import ImageNormalize\n",
"\n",
"from IPython.display import HTML"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "84821bbe-84f0-4e74-9f20-bfab05cd9c5a",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"4.1.2\n",
"1.6.2\n"
]
}
],
"source": [
"import sunpy\n",
"print(sunpy.__version__)\n",
"print(sunpy_soar.__version__)"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "afe8629b-586b-49f8-b44d-38b24a323ef3",
"metadata": {},
"outputs": [],
"source": [
"time_range = a.Time('2022-03-02 12:00:22','2022-03-06 11:40:50')"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "a4301311-05f8-4b3a-a18d-9f699b469cd6",
"metadata": {},
"outputs": [],
"source": [
"q = Fido.search(time_range, a.Level(2), a.Provider.soar, a.Instrument.eui, a.soar.Product('EUI-FSI174-IMAGE'))"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "77b952c6-477d-48d0-9c68-5748a9cb29e5",
"metadata": {
"collapsed": true,
"jupyter": {
"outputs_hidden": true
},
"tags": []
},
"outputs": [
{
"data": {
"text/html": [
"Results from 1 Provider:</br></br>567 Results from the SOARClient:</br><div><i>QueryResponseTable length=567</i>\n",
"<table id=\"table5178981344\" class=\"table-striped table-bordered table-condensed\">\n",
"<thead><tr><th>Instrument</th><th>Data product</th><th>Level</th><th>Start time</th><th>End time</th><th>Data item ID</th><th>Filename</th><th>Filesize</th></tr></thead>\n",
"<thead><tr><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th>Mbyte</th></tr></thead>\n",
"<thead><tr><th>str3</th><th>str16</th><th>str2</th><th>str23</th><th>str23</th><th>str43</th><th>str52</th><th>float64</th></tr></thead>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-06 11:40:45.269</td><td>2022-03-06 11:40:55.269</td><td>solo_L2_eui-fsi174-image_20220306T114045269</td><td>solo_L2_eui-fsi174-image_20220306T114045269_V01.fits</td><td>3.623</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-06 11:30:45.247</td><td>2022-03-06 11:30:55.247</td><td>solo_L2_eui-fsi174-image_20220306T113045247</td><td>solo_L2_eui-fsi174-image_20220306T113045247_V01.fits</td><td>3.603</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-06 11:20:45.245</td><td>2022-03-06 11:20:55.245</td><td>solo_L2_eui-fsi174-image_20220306T112045245</td><td>solo_L2_eui-fsi174-image_20220306T112045245_V01.fits</td><td>3.709</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-06 11:10:45.260</td><td>2022-03-06 11:10:55.260</td><td>solo_L2_eui-fsi174-image_20220306T111045260</td><td>solo_L2_eui-fsi174-image_20220306T111045260_V01.fits</td><td>3.649</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-06 11:00:45.243</td><td>2022-03-06 11:00:55.243</td><td>solo_L2_eui-fsi174-image_20220306T110045243</td><td>solo_L2_eui-fsi174-image_20220306T110045243_V01.fits</td><td>3.603</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-06 10:50:45.242</td><td>2022-03-06 10:50:55.242</td><td>solo_L2_eui-fsi174-image_20220306T105045242</td><td>solo_L2_eui-fsi174-image_20220306T105045242_V01.fits</td><td>3.637</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-06 10:40:45.258</td><td>2022-03-06 10:40:55.258</td><td>solo_L2_eui-fsi174-image_20220306T104045258</td><td>solo_L2_eui-fsi174-image_20220306T104045258_V01.fits</td><td>3.635</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-06 10:30:45.240</td><td>2022-03-06 10:30:55.240</td><td>solo_L2_eui-fsi174-image_20220306T103045240</td><td>solo_L2_eui-fsi174-image_20220306T103045240_V01.fits</td><td>3.663</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-06 10:20:45.239</td><td>2022-03-06 10:20:55.239</td><td>solo_L2_eui-fsi174-image_20220306T102045239</td><td>solo_L2_eui-fsi174-image_20220306T102045239_V01.fits</td><td>3.686</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-06 10:10:45.238</td><td>2022-03-06 10:10:55.238</td><td>solo_L2_eui-fsi174-image_20220306T101045238</td><td>solo_L2_eui-fsi174-image_20220306T101045238_V01.fits</td><td>3.658</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-06 10:00:45.237</td><td>2022-03-06 10:00:55.237</td><td>solo_L2_eui-fsi174-image_20220306T100045237</td><td>solo_L2_eui-fsi174-image_20220306T100045237_V01.fits</td><td>3.623</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-06 09:50:45.236</td><td>2022-03-06 09:50:55.236</td><td>solo_L2_eui-fsi174-image_20220306T095045236</td><td>solo_L2_eui-fsi174-image_20220306T095045236_V01.fits</td><td>3.675</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-06 09:40:45.235</td><td>2022-03-06 09:40:55.235</td><td>solo_L2_eui-fsi174-image_20220306T094045235</td><td>solo_L2_eui-fsi174-image_20220306T094045235_V01.fits</td><td>3.689</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-06 09:30:45.234</td><td>2022-03-06 09:30:55.234</td><td>solo_L2_eui-fsi174-image_20220306T093045234</td><td>solo_L2_eui-fsi174-image_20220306T093045234_V01.fits</td><td>3.658</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-06 09:20:45.232</td><td>2022-03-06 09:20:55.232</td><td>solo_L2_eui-fsi174-image_20220306T092045232</td><td>solo_L2_eui-fsi174-image_20220306T092045232_V01.fits</td><td>3.698</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-06 09:10:45.232</td><td>2022-03-06 09:10:55.232</td><td>solo_L2_eui-fsi174-image_20220306T091045232</td><td>solo_L2_eui-fsi174-image_20220306T091045232_V01.fits</td><td>3.646</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-06 09:00:45.231</td><td>2022-03-06 09:00:55.231</td><td>solo_L2_eui-fsi174-image_20220306T090045231</td><td>solo_L2_eui-fsi174-image_20220306T090045231_V01.fits</td><td>3.678</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-06 08:50:45.230</td><td>2022-03-06 08:50:55.230</td><td>solo_L2_eui-fsi174-image_20220306T085045230</td><td>solo_L2_eui-fsi174-image_20220306T085045230_V01.fits</td><td>3.637</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-06 08:40:45.228</td><td>2022-03-06 08:40:55.228</td><td>solo_L2_eui-fsi174-image_20220306T084045228</td><td>solo_L2_eui-fsi174-image_20220306T084045228_V01.fits</td><td>3.678</td></tr>\n",
"<tr><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-02 17:00:17.301</td><td>2022-03-02 17:00:27.301</td><td>solo_L2_eui-fsi174-image_20220302T170017301</td><td>solo_L2_eui-fsi174-image_20220302T170017301_V01.fits</td><td>3.519</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-02 16:45:17.300</td><td>2022-03-02 16:45:27.300</td><td>solo_L2_eui-fsi174-image_20220302T164517300</td><td>solo_L2_eui-fsi174-image_20220302T164517300_V01.fits</td><td>3.433</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-02 16:30:17.297</td><td>2022-03-02 16:30:27.297</td><td>solo_L2_eui-fsi174-image_20220302T163017297</td><td>solo_L2_eui-fsi174-image_20220302T163017297_V01.fits</td><td>3.398</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-02 16:15:17.296</td><td>2022-03-02 16:15:27.296</td><td>solo_L2_eui-fsi174-image_20220302T161517296</td><td>solo_L2_eui-fsi174-image_20220302T161517296_V01.fits</td><td>3.473</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-02 16:00:17.294</td><td>2022-03-02 16:00:27.294</td><td>solo_L2_eui-fsi174-image_20220302T160017294</td><td>solo_L2_eui-fsi174-image_20220302T160017294_V01.fits</td><td>3.398</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-02 15:45:17.293</td><td>2022-03-02 15:45:27.293</td><td>solo_L2_eui-fsi174-image_20220302T154517293</td><td>solo_L2_eui-fsi174-image_20220302T154517293_V01.fits</td><td>3.462</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-02 15:30:17.291</td><td>2022-03-02 15:30:27.291</td><td>solo_L2_eui-fsi174-image_20220302T153017291</td><td>solo_L2_eui-fsi174-image_20220302T153017291_V01.fits</td><td>3.462</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-02 15:15:17.290</td><td>2022-03-02 15:15:27.290</td><td>solo_L2_eui-fsi174-image_20220302T151517290</td><td>solo_L2_eui-fsi174-image_20220302T151517290_V01.fits</td><td>3.421</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-02 15:00:17.288</td><td>2022-03-02 15:00:27.288</td><td>solo_L2_eui-fsi174-image_20220302T150017288</td><td>solo_L2_eui-fsi174-image_20220302T150017288_V01.fits</td><td>3.378</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-02 14:45:17.286</td><td>2022-03-02 14:45:27.286</td><td>solo_L2_eui-fsi174-image_20220302T144517286</td><td>solo_L2_eui-fsi174-image_20220302T144517286_V01.fits</td><td>3.427</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-02 14:30:17.284</td><td>2022-03-02 14:30:27.284</td><td>solo_L2_eui-fsi174-image_20220302T143017284</td><td>solo_L2_eui-fsi174-image_20220302T143017284_V01.fits</td><td>3.384</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-02 14:15:17.283</td><td>2022-03-02 14:15:27.283</td><td>solo_L2_eui-fsi174-image_20220302T141517283</td><td>solo_L2_eui-fsi174-image_20220302T141517283_V01.fits</td><td>3.421</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-02 14:00:17.306</td><td>2022-03-02 14:00:27.306</td><td>solo_L2_eui-fsi174-image_20220302T140017306</td><td>solo_L2_eui-fsi174-image_20220302T140017306_V01.fits</td><td>3.514</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-02 13:45:17.279</td><td>2022-03-02 13:45:27.279</td><td>solo_L2_eui-fsi174-image_20220302T134517279</td><td>solo_L2_eui-fsi174-image_20220302T134517279_V01.fits</td><td>3.384</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-02 13:30:17.281</td><td>2022-03-02 13:30:27.281</td><td>solo_L2_eui-fsi174-image_20220302T133017281</td><td>solo_L2_eui-fsi174-image_20220302T133017281_V01.fits</td><td>3.375</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-02 13:15:17.725</td><td>2022-03-02 13:15:27.725</td><td>solo_L2_eui-fsi174-image_20220302T131517725</td><td>solo_L2_eui-fsi174-image_20220302T131517725_V01.fits</td><td>3.347</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-02 13:00:17.275</td><td>2022-03-02 13:00:27.275</td><td>solo_L2_eui-fsi174-image_20220302T130017275</td><td>solo_L2_eui-fsi174-image_20220302T130017275_V01.fits</td><td>3.381</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-02 12:45:17.546</td><td>2022-03-02 12:45:27.546</td><td>solo_L2_eui-fsi174-image_20220302T124517546</td><td>solo_L2_eui-fsi174-image_20220302T124517546_V01.fits</td><td>3.396</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-02 12:30:17.273</td><td>2022-03-02 12:30:27.273</td><td>solo_L2_eui-fsi174-image_20220302T123017273</td><td>solo_L2_eui-fsi174-image_20220302T123017273_V01.fits</td><td>3.312</td></tr>\n",
"<tr><td>EUI</td><td>EUI-FSI174-IMAGE</td><td>L2</td><td>2022-03-02 12:15:17.346</td><td>2022-03-02 12:15:27.346</td><td>solo_L2_eui-fsi174-image_20220302T121517346</td><td>solo_L2_eui-fsi174-image_20220302T121517346_V01.fits</td><td>3.303</td></tr>\n",
"</table></div></br>"
],
"text/plain": [
"<sunpy.net.fido_factory.UnifiedResponse object at 0x134b0f4f0>\n",
"Results from 1 Provider:\n",
"\n",
"567 Results from the SOARClient:\n",
"\n",
"Instrument Data product Level Start time End time Filesize\n",
" Mbyte \n",
"---------- ---------------- ----- ----------------------- ----------------------- --------\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-06 11:40:45.269 2022-03-06 11:40:55.269 3.623\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-06 11:30:45.247 2022-03-06 11:30:55.247 3.603\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-06 11:20:45.245 2022-03-06 11:20:55.245 3.709\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-06 11:10:45.260 2022-03-06 11:10:55.260 3.649\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-06 11:00:45.243 2022-03-06 11:00:55.243 3.603\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-06 10:50:45.242 2022-03-06 10:50:55.242 3.637\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-06 10:40:45.258 2022-03-06 10:40:55.258 3.635\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-06 10:30:45.240 2022-03-06 10:30:55.240 3.663\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-06 10:20:45.239 2022-03-06 10:20:55.239 3.686\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-06 10:10:45.238 2022-03-06 10:10:55.238 3.658\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-06 10:00:45.237 2022-03-06 10:00:55.237 3.623\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-06 09:50:45.236 2022-03-06 09:50:55.236 3.675\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-06 09:40:45.235 2022-03-06 09:40:55.235 3.689\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-06 09:30:45.234 2022-03-06 09:30:55.234 3.658\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-06 09:20:45.232 2022-03-06 09:20:55.232 3.698\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-06 09:10:45.232 2022-03-06 09:10:55.232 3.646\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-06 09:00:45.231 2022-03-06 09:00:55.231 3.678\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-06 08:50:45.230 2022-03-06 08:50:55.230 3.637\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-06 08:40:45.228 2022-03-06 08:40:55.228 3.678\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-06 08:30:45.227 2022-03-06 08:30:55.227 3.612\n",
" ... ... ... ... ... ...\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-02 17:00:17.301 2022-03-02 17:00:27.301 3.519\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-02 16:45:17.300 2022-03-02 16:45:27.300 3.433\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-02 16:30:17.297 2022-03-02 16:30:27.297 3.398\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-02 16:15:17.296 2022-03-02 16:15:27.296 3.473\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-02 16:00:17.294 2022-03-02 16:00:27.294 3.398\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-02 15:45:17.293 2022-03-02 15:45:27.293 3.462\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-02 15:30:17.291 2022-03-02 15:30:27.291 3.462\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-02 15:15:17.290 2022-03-02 15:15:27.290 3.421\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-02 15:00:17.288 2022-03-02 15:00:27.288 3.378\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-02 14:45:17.286 2022-03-02 14:45:27.286 3.427\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-02 14:30:17.284 2022-03-02 14:30:27.284 3.384\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-02 14:15:17.283 2022-03-02 14:15:27.283 3.421\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-02 14:00:17.306 2022-03-02 14:00:27.306 3.514\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-02 13:45:17.279 2022-03-02 13:45:27.279 3.384\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-02 13:30:17.281 2022-03-02 13:30:27.281 3.375\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-02 13:15:17.725 2022-03-02 13:15:27.725 3.347\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-02 13:00:17.275 2022-03-02 13:00:27.275 3.381\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-02 12:45:17.546 2022-03-02 12:45:27.546 3.396\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-02 12:30:17.273 2022-03-02 12:30:27.273 3.312\n",
" EUI EUI-FSI174-IMAGE L2 2022-03-02 12:15:17.346 2022-03-02 12:15:27.346 3.303\n",
"Length = 567 rows\n"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"q"
]
},
{
"cell_type": "markdown",
"id": "874e7203-0b04-44fe-afa0-8a1207f8255a",
"metadata": {},
"source": [
"Of the >500 available in this time range, I'm only going to download every 100th one (for a total of 6 files) but this method should work on all of them as well."
]
},
{
"cell_type": "code",
"execution_count": 8,
"id": "8ed46ba8-a7cf-4d2f-aeb4-319965fa65dc",
"metadata": {
"tags": []
},
"outputs": [
{
"data": {
"application/vnd.jupyter.widget-view+json": {
"model_id": "13d719d8cbfd479e8f40ed3eae36d987",
"version_major": 2,
"version_minor": 0
},
"text/plain": [
"Files Downloaded: 0%| | 0/6 [00:00<?, ?file/s]"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"files = Fido.fetch(q[0][::-1][::100])"
]
},
{
"cell_type": "markdown",
"id": "2f76dc93-1ff9-41e8-8603-1681a7c9d5f8",
"metadata": {},
"source": [
"Load the maps into sunpy Map objects"
]
},
{
"cell_type": "code",
"execution_count": 9,
"id": "32fc30da-ac97-4219-aa6e-70a6cde9689d",
"metadata": {},
"outputs": [],
"source": [
"maps = sunpy.map.Map(files)"
]
},
{
"cell_type": "markdown",
"id": "619f62a0-4cc3-4d0f-bc1e-8cd94d73081d",
"metadata": {},
"source": [
"I identified this manually the lower left corner of the region of interest using the first image (as indicated by the index 0). You can redo this selection using an image at any timestep. Everything should work the same.\n",
"\n",
"I also estimated a width and a height for the region and used that to compute a center of the box."
]
},
{
"cell_type": "code",
"execution_count": 15,
"id": "da00fa55-00a9-4169-94f8-994cc12ffe09",
"metadata": {},
"outputs": [],
"source": [
"blc = SkyCoord(Tx=-1000*u.arcsec, Ty=-1000*u.arcsec,frame=maps[0].coordinate_frame)\n",
"width = 10 * u.arcsec\n",
"height = 2000 * u.arcsec\n",
"center = SkyCoord(Tx=blc.Tx+width/2, Ty=blc.Ty+height/2, frame=blc.frame)"
]
},
{
"cell_type": "markdown",
"id": "45e3d05e-9d02-413a-8dd2-72de77ea68ea",
"metadata": {},
"source": [
"The following plot shows my identified center and the box denoted by the width and the height. Note that without the rotational transformation applied when we use `propagate_with_solar_surface`, the point will not rotate across the disk (blue point).\n",
"\n",
"The red box denotes the cutout region of interest where the center of the box is the transformed center with rotation included."
]
},
{
"cell_type": "code",
"execution_count": 16,
"id": "43bb0551-5139-40ec-ae13-c836147e7914",
"metadata": {},
"outputs": [
{
"data": {
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment