Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1420 1716 (N37W46 L=239)   341     65        0 (1399/1691-1452/1736)   9.5 W43-W49 N35-N38
      2  740 1642 (N44E10 L=182)   494    122        2 ( 692/1604- 794/1670)   6.3 E14-E05 N42-N46
      3  998 1492 (N29W06 L=199)  2352    510       11 ( 793/1382-1195/1579)  10.1 E08-W20 N24-N32
      4 1697 1503 (N18W58 L=250)   868    141        3 (1646/1455-1758/1555) -41.2 W54-W62 N14-N22
      5  521 1424 (N33E29 L=164)  1330    271        5 ( 415/1393- 665/1456) -13.4 E38-E17 N29-N36
      6 1916  981 (S16W67 L=260)  2181    283        2 (1882/ 843-1948/1088)  57.3 W63-W73 S23-S09
      7 1838 1012 (S13W57 L=249)   544     74        1 (1822/ 968-1856/1045) -68.2 W55-W59 S15-S10
      8 1677  930 (S15W41 L=233)   308     87        0 (1633/ 904-1708/ 956)   9.5 W37-W43 S15-S13
      9 1185  906 (S09W08 L=200)   359    157        0 (1140/ 839-1245/ 945)  35.5 W04-W11 S12-S07
     10 1821  806 (S25W55 L=248)  1591    241        3 (1765/ 756-1893/ 898) -48.0 W49-W66 S28-S18
     11  708  594 (S20E27 L=165)  1246    485        4 ( 553/ 480- 855/ 732)  22.2 E41-E14 S25-S14
     12  382  509 (S21E57 L=136)   349     48        1 ( 366/ 485- 409/ 529)  15.9 E59-E52 S21-S19
  

    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