Hi People,

I am trying to play around with FITS file format. I want to create an image from this file which will include data from center to 20 arc minutes. each pixel is 3x3 arc sec.

import yt

ds = yt.load("signal_wfilt_sn.fits")

ds.print_stats()

ds.field_list

ds.domain_width

ds.derived_field_list

dd = ds.all_data()

sp = ds.sphere(ds.domain_center, (400, 'code_length')) # create sphere of 20 arcminutes with 3x3 pixel size

sp[("fits", "image_0")] # check array of data within the sphere

#plot 1D histogram of data over pixel

pplot = yt.ProfilePlot(sp, "image_0", ["pixel"], weight_field=None, n_bins=1024,

plot_spec=dict(color='red', linestyle="--"), y_log={'pixel':False})

pplot.show()

# convert data to numpy array

import numpy as np

image = np.array(sp['image_0'])

# create image?

imgplot = plt.imshow(image)

imgplot

plt.show()

image

The code is failing at this point with error:

---------------------------------------------------------------------------
TypeError                                 Traceback (most recent call last)
<ipython-input-21-6a83e68fbbf2> in <module>()
----> 1 imgplot = plt.imshow(image)
      2 
      3 imgplot
      4 plt.show()

/home/trina/anaconda2/lib/python2.7/site-packages/matplotlib/pyplot.pyc in imshow(X, cmap, norm, aspect, interpolation, alpha, vmin, vmax, origin, extent, shape, filternorm, filterrad, imlim, resample, url, hold, data, **kwargs)
   3155                         filternorm=filternorm, filterrad=filterrad,
   3156                         imlim=imlim, resample=resample, url=url, data=data,
-> 3157                         **kwargs)
   3158     finally:
   3159         ax._hold = washold

/home/trina/anaconda2/lib/python2.7/site-packages/matplotlib/__init__.pyc in inner(ax, *args, **kwargs)
   1896                     warnings.warn(msg % (label_namer, func.__name__),
   1897                                   RuntimeWarning, stacklevel=2)
-> 1898             return func(ax, *args, **kwargs)
   1899         pre_doc = inner.__doc__
   1900         if pre_doc is None:

/home/trina/anaconda2/lib/python2.7/site-packages/matplotlib/axes/_axes.pyc in imshow(self, X, cmap, norm, aspect, interpolation, alpha, vmin, vmax, origin, extent, shape, filternorm, filterrad, imlim, resample, url, **kwargs)
   5122                               resample=resample, **kwargs)
   5123 
-> 5124         im.set_data(X)
   5125         im.set_alpha(alpha)
   5126         if im.get_clip_path() is None:

/home/trina/anaconda2/lib/python2.7/site-packages/matplotlib/image.pyc in set_data(self, A)
    598         if (self._A.ndim not in (2, 3) or
    599                 (self._A.ndim == 3 and self._A.shape[-1] not in (3, 4))):
--> 600             raise TypeError("Invalid dimensions for image data")
    601 
    602         self._imcache = None

TypeError: Invalid dimensions for image data

I am completely new to image processing, so I am interested to know how you create image for this type of data in yt. Your help is very appreciated. I have also attached the jupyter notebook and the actual data I am using for the processing.

Thanks a lot

Tazkera Haque