0

RGB私は画像のセットを持っていて、画像にRGBフィットする画像を作成aplpyし、画像にいくつかの輪郭を重ねましたが、画像を切り取り、輪郭のピークが見える小さなフィット画像を作成したいと思います。

import aplpy
import atpy
from pyavm import AVM
import asciitable
import matplotlib
import matplotlib.pyplot as plt
from matplotlib.colors import LogNorm,BoundaryNorm
import montage_wrapper
from astropy.io import fits
import pyfits
from astropy import wcs
fitsfile = 'rgb.fits'
fitsfile_2d = 'rgb_2d.fits'
pngfile = 'rgb_arcsinh_contour.png'

figfile = 'rgb.png'
w = wcs.WCS(naxis=2)

# Set up an "Airy's zenithal" projection
# Vector properties may be set with Python lists, or Numpy arrays
w.wcs.crpix = [5.70000000E+03, 3.05000000E+03]
w.wcs.cdelt = np.array([-6.611111263E-05, 6.611111263E-05])
w.wcs.crval = [23.166667, -7.666667]
w.wcs.ctype = ["RA---TAN", "DEC--TAN"]
w.wcs.cunit =["deg","deg"]

# Print out all of the contents of the WCS object
w.wcs.print_contents()

# Some pixel coordinates of interest.
pixcrd = np.array([[0,0],[24,38],[45,98]], np.float_)

# Convert pixel coordinates to world coordinates
world = w.wcs_pix2world(pixcrd, 1)
print world

# Convert the same coordinates back to pixel coordinates.
pixcrd2 = w.wcs_world2pix(world, 1)
print pixcrd2

# These should be the same as the original pixel coordinates, modulo
# some floating-point error.
assert np.max(np.abs(pixcrd - pixcrd2)) < 1e-6

# Now, write out the WCS object as a FITS header
header = w.to_header()
hdu = pyfits.open(fitsfile)
# header is an astropy.io.fits.Header object.  We can use it to create a new
# PrimaryHDU and write it to a file.
hdu = fits.PrimaryHDU(header=header)


# make rgb image
aplpy.make_rgb_image(fitsfile, pngfile,
                     vmin_r=-0.005, vmax_r=0.2,
                     vmin_g=-0.02, vmax_g=0.1,
                     vmin_b=-0.02,vmax_b=0.04,
                     embed_avm_tags=False)

# make a figure
img = aplpy.FITSFigure(fitsfile_2d)
img.show_rgb(pngfile)
img.set_nan_color('white')
standard_setup(img)

指定されたサイズの画像の指定された座標からサブプロットを作成するにはどうすればよいですか?

4

1 に答える 1

1

APLpy のドキュメントによると、サブプロットを作成できます。

既定では、FITSFigure は Figure 全体を占める単一のサブプロットをもつ Figure を作成します。ただし、APLpy を使用して、既存の matplotlib Figure インスタンスにサブプロットを配置できます。これを行うには、次のように figure= 引数を指定して FITSFigure を呼び出す必要があります。

import aplpy
import matplotlib.pyplot as mpl

fig = mpl.figure()
f = aplpy.FITSFigure('some_image.fits', figure=fig)

そしてあなたの数字を中心に合わせます：

Figure は、ズームとパンによって対話的に探索できます。プログラムで特定の領域を中心に合わせるには、半径を指定して次の方法を使用します。

fig.recenter(33.23, 55.33, radius=0.3) # degrees

または別の幅と高さ:

fig.recenter(33.23, 55.33, width=0.3, height=0.2) # degrees

私もAPLpyを使用しているのでテストしましたが、うまく機能します。

HTH、

ドイツ人。

于 2014-12-02T02:33:09.200 に答える