Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1388 1726 (N35W54 L=141)  2419    306        2 (1325/1607-1461/1823)  19.0 W44-W75 N27-N39
      2 1726 1044 (S12W44 L=132)   282     37        1 (1709/1020-1743/1070) -56.3 W42-W45 S13-S10
      3 1622  893 (S17W33 L=120)  1017    219        0 (1597/ 790-1641/ 998)  57.5 W29-W36 S22-S11
      4 1897  972 (S23W62 L=150)   666     90        1 (1864/ 946-1921/ 998)   8.1 W57-W67 S24-S20
      5  445  972 (N18E37 L= 50)   248     52        0 ( 417/ 952- 469/ 996) -36.9 E39-E35 N16-N19
      6  231  848 (N14E60 L= 27)  1202    171        0 ( 196/ 776- 262/ 925) -78.7 E63-E56 N09-N19
      7  292  743 (N07E57 L= 31)  1055    196        6 ( 254/ 670- 334/ 839)  59.0 E62-E50 N02-N12
      8  579  620 (S06E38 L= 50)   244     57        0 ( 549/ 601- 604/ 642)  14.0 E39-E35 S06-S04
      9 1456  546 (S34W14 L=102)   323     58        1 (1428/ 519-1475/ 576)  26.6 W11-W16 S35-S31
     10  968  425 (S27E21 L= 66)   416     60        0 ( 933/ 397- 995/ 453)  18.4 E24-E18 S28-S25
     11 1297  258 (S52E07 L= 80)   449     86        0 (1255/ 233-1336/ 276)  -6.3 E12-E03 S53-S50
     12  792  246 (S35E44 L= 43)   481     44        0 ( 765/ 223- 821/ 271)  14.0 E48-E40 S36-S33
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y