Cut out poles from synoptic magnetograms

[Paolo Pagano]

sorry in advance for our silly this may sound.

I do have a synoptic magnetogram map extracted from a fits file in this way

hmi_map = sunpy.map.Map('/Users/paolopagano/Desktop/Solar/SolO_Eruption_042021/hmi.synoptic_mr_polfil_720s.2246.Mr_polfil.fits')

in this way I get a 2D array of 

> print('shape: ', hmi_map.data.shape)

shape:  (1440, 3600)

from which I would like to cut off the northernmost and southernmost bands, so that it will become, let's say

 (1440, 3400)

if I remove the 100 northernmost pixels and the 100 southernmost pixels.

anything I have tried hasn't succeeded.

Any suggestion?

Paolo 

[Albert Y. Shih]

Hi, Paolo,

    Sorry for the slow response!  Can you elaborate on what you have tried? You’ll likely want to use the `submap()` method of your `hmi_map` object:

https://docs.sunpy.org/en/stable/generated/api/sunpy.map.GenericMap.html#sunpy.map.GenericMap.submap

Also, I’ll point out that the dimension with 3600 pixels is the longitude direction, so that will not be the dimension that gets smaller if you cut off the north and south poles.

Albert

--
You received this message because you are subscribed to the Google Groups "SunPy" group.
To unsubscribe from this group and stop receiving emails from it, send an email to sunpy+un...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/sunpy/a599c103-4eec-4d96-99ea-9731230734bdn%40googlegroups.com.