python - To find the pixel coodinates for a given ra-dec -

-

i'm new python, , presently working on labeling galaxy images. have set of 512x512 pix^2 images scaled 0.2 arcsec/pix. each 1 of images have 2 objects labeled - 1 @ centre (256pix,256pix) @ other @ offset (i've projected separation between objects ra-dec values both objects). task encircle objects in each of images. easy central objects, i'm stuck on how same thing done other object. following snippet encircling of central object, need figure out how other one:

i have ra-dec values both objects (in degrees).
have separation between objects (in kpc).

for image 512x512 px @ .2 ''/px:

this iterates through objids in list - each objid have file objid.png in folder called images_fin - loads correct image, , labels objid.

for galaxy in range(0, len(objid)):

![enter image description here][1]im=imread('images_fin/'+objid[galaxy]+'.jpeg') imshow(im) a=gca() print a.text(20, 480, "objid:", color ='w', fontsize='10') print a.text(20, 490, "w1-w2:", color ='w', fontsize='10') print a.text(20, 500, "sep1 (in kpc):", color ='w', fontsize='10') print a.text(20, 510, "sep2 (in kpc):", color ='w', fontsize='10') print a.text(52, 480, objid[galaxy], color='w', fontsize='10') print a.text(52, 490, w[galaxy], color='w', fontsize='10') print a.text(78, 500, sep1[galaxy], color='w', fontsize='10') print a.text(78, 510, sep2[galaxy], color='w', fontsize='10')  print a.hlines(y=20, xmin=480, xmax=495, linewidth=1, color = 'w') print a.vlines(x=480, ymin=18, ymax=22, linewidth=1, color = 'w') print a.vlines(x=495, ymin=18, ymax=22, linewidth=1, color = 'w')  print a.text(485, 18, "3\"", color='w', fontsize='10')  a.xaxis.set_visible(false) a.yaxis.set_visible(false) circle=plt.circle((256,256),12,color='w',fill=false) plt.gcf().gca().add_artist(circle)

i'm pasting example how labeling should when done (this done manually)


Comments

Post a Comment

Popular posts from this blog

python - Healpy: From Data to Healpix map -

-
Image
i have data grid rows represent theta (0, pi) , columns represent phi (0, 2*pi) , f(theta,phi) density of dark matter @ location. wanted calculate power spectrum , have decided use healpy. what can not understand how format data healpy use. if provide code (in python obvious reasons) or point me tutorial, great! have tried hand @ doing following code: #grid dimensions nrows*ncols (subject change) theta = np.linspace(0, np.pi, num=grid.shape[0])[:, none] phi = np.linspace(0, 2*np.pi, num=grid.shape[1]) nside = 512 print "pixel area: %.2f square degrees" % hp.nside2pixarea(nside, degrees=true) pix = hp.ang2pix(nside, theta, phi) healpix_map = np.zeros(hp.nside2npix(nside), dtype=np.double) healpix_map[pix] = grid but, when try execute code power spectrum. specifically, : cl = hp.anafast(healpix_map[pix], lmax=1024) i error: typeerror: bad number of pixels if point me tutorial or edit code great. more specifications: data in 2d np array , can change numrows/...
Read more

c - Bitwise operation with (signed) enum value -

-
i using enumerator values flags: typedef enum { = 0x00, b = 0x01u, // u has no influence, expected c = 0x02u, // u has no influence, expected ... } enum_name; volatile unsigned char* reg = someaddress; *reg |= b; according misra-c:2004 bitwise operations shall not done signed type. unfortunately, compiler iar use signed int (or short or char) underlying type of enums, , option can find relates size, not signedness ("--enum-is-int"). according iar c/c++ development guide arm , pages 169 , 211, can define type of enum s if enable iar language extensions ( -e command-line option, or project > options > c/c++ compiler > language > allow iar extensions in ide). in particular, should define "sentinel" value, make sure compiler chooses correct type. prefers signed types, , uses smallest possible integer type, sentinel should largest positive integer corresponding unsigned integer type can describe. example, typedef enum ...
Read more

xslt - Unnest parent nodes by child node -

-
in xml there items 0-n attributes , item should copied each attribute new item 1 attribute. given xml this: <?xml version="1.0" encoding="utf-8" ?> <items> <item> <name>a</name> <attributes> <attribute> <key>attribute1</key> <value>1</value> </attribute> </attributes> </item> <item> <name>b</name> </item> <item> <name>c</name> <attributes> <attribute> <key>attribute1</key> <value>5</value> </attribute> <attribute> <key>attribute2</key> <value>2</value> </attribute> <attribute> <key>attri...
Read more