Find physical coordinates of a pixel in a fits file with python -

-

I am stuck in a python script to get the physical sky coordinates of a given pixel from the wcs of astropy I would like to use, but I will do anything from within the dragon

I have tried these two snippets of code. Astropy.io Import fits from 

  to WCS def astropymethod1 (IMG) from ASTROP. WESSES: # from http://astropy.readthedocs.org/en/latest/wcs/w = WCS (IMG) lawn, lat = w.all_pix2world (100., 100., 1) Print Lawn, Let Df Astropaimeth 2 (IMG): # from http://astropy.readthedocs.org/en/latest/wcs/hdu = Fits.open (IMG) w = WCS (HDU [0] .hader) Lone, LET = w.wcs_pix2world (100., 100., 1) Print Lone, LET   
 me The first time I get the error when I try to use WCS for the first time and the result is only pixel value at any time. I've gotten. Warning: FITSFixed Warning: WCS conversion is associated with more axes (2) images (0) [astropy.wcs.wcs]    < P> 
 
 The problem is that you have a multi-extension FITS file. Here's an example session showing how you can gain access to the appropriate WCS:  
 in : fits into astropy.io import from [2]: h = fits . Structure ('SN1415_F625W_1_drz Fits') [3]: f = fits.open ('SN1415_F625W_1_drz.fits') [4]: ​​F Out [4] [& LT; Astropy.io.fits.hdu.image.PrimaryHDU at 0x106735490>, & lt; Astropy.io.fits.hdu.image.ImageHDU at 0x106749750>, & lt; Astropy.io.fits.hdu.image.ImageHDU> 0x106751310, & lt; Astropy.io.fits.hdu.image.ImageHDU 0x106751d10 & gt;; & lt; 0x1067dfdd0 & gt; astropy.io.s..hdu.table.BinTableHDU;] [5]: Astropy Import from WCS [6]: w = wcs.WCS (f [0] .header) WARNING: FITSFixed Warning: In WCS conversion In the image (0) [astropy.wcs.wcs] [7]: w.wcs.naxis [7]: 2 [8]: F. [0] has more axis than the image associated with [0]. In data [9]: w = wcs.WCS (F [1] .header) [10]: w.wcs.naxis out [10]: 2 in [11]: F [1] .Data out [11] : Array ([[0.0198678, -0.04018363, 0.03330525, ..., 0., 0., 0.], [0.0695872, -0.00 9 7, 143, 0.00147662, ..., 0., 0 ., 0.], [-0.09292094, 0.02481506, -0.01057338, ..., 0., 0., 0.] ... ..., [0., 0., 0., ..., 0.02375774 , 0.0389731, 0.03825707], [0., 0., 0., ..., -0.015709-18, -0.01053802, 0.00461219], [0., 0., 0., ..., -0.0638448, -0.0240754 , 0.02679451]], dtype = float32 [12] :. W.wcs_pix2world (100, 100, 1) out [12] [array (6.113076380801787), array (0.616758775753701)]   
 So you probably want to redefine your method: / P> 
  def astropymethod2 (IMG, HDID = 1): # from http://astropy.readthedocs.org/en/latest/wcs/hdu=fit (IMG) w = wcs (hdu [hduid]. Header) lon, lat = w.wcs_pix2world (100., 100., 1) Print lawn, late

Comments

Post a Comment

Popular posts from this blog

ios - Adding an SKSpriteNode to SKScene from a child SKSpriteNode -

-
I have a SKScene where I am adding an SKSpriteNode I have subclassed the SKSpriteNode class to create this node. In the subclass I am defining some SKActions on Sprite. What do I want to do when the end of this phantom ends the SKAction sequence, then I add a new Prayer node to this scene. How is this possible. The following code is mine: The code for the sequence is that I'm running on the Scansprinitode subclass (TEMissileNode): - SKAction * moveDown = [SKAction moveToY: Self.position.y - 20 Duration: 0.2]; SKAction * animation = [SKAction animateWithTextures: textures timePerFrame: time / 7]; SKAction * moveMissileProjectile = [SKAction Moving: Pointoffscreen Term: Time]; SKAction * group = [SKAction group: @ [animation, movement projectile]]; SKAction * Sequence = [SKAction Sequence: @ [Hilldown, Group, [SKAction removalFromParent]]]; [Self run action: sequence]; From the main scene, I call those actions that execute these tasks TEMissileNode * missile = [tmissil...
Read more

Matlab transpose a table vector -

-
An embarrassing simple question seems to be, but how can I move a Matlab table vector? aTableT = aTable '; I tried a standard syntax for the simple interaction of a line vector to a line vector aTable : ATableT = reshape (aTable, 1, Height (aTable)); and aTableT = rot90 (aTable); According to the head, the last time table should work for table, see. However, I get this error code: Error type in using table / permit (line 396) Undefined function 'permute' for the input arguments. Error in rot 90 (line 29) b = magnitude (b, [2 1 3: ndims (a)]); NB: fliplr is not useful either pretty sure I have covered clear angles - any ideas? Thanks! Try changing your table into an array, move it, then back to a table In other words, try doing this: aTableArray = table2array (aTable); ATableT = array 2table (aTableArray. '); I have also read the document for rot90 , and it says that rot90 is definitely working for tables Do, and I find you...
Read more

c# - Textbox not clickable but editable -

-
I have a small form with 10 text boxes, I have set them in the correct tab order. I want them on my tab service I was wondering if there is a way to set the text box so that they can be selected for editing until they are tabbed. That is ... I do not want the end user to edit the text box to click on them, I just want to make them editable through tabu. this should do the trick public partial category poor text box: Text Box {Safe Override Zero WndProc (Ref: Message Message) {if (m.msg == (Int) WM.LBUTTONDOWN) {Return; // mount down events} base.WndProc (ref m); }} Window message enum can be found. How to do this without anherizing: text box : class EatMouseDown: netwindows {safe override zero WndProc (ref message message) {if (m.msg == (int) WM.LBUTTONDOWN) {return; } Base.WandProc (Ref M); }} Safe Override Zero Onload (EventEurge E) {base.OnLoad (e); How to do this without any part: left a clean part, whatever is important, this can be a buggy but it works R...
Read more