Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  769 1706 (N30E34 L= 66)  1007    233        2 ( 671/1642- 845/1775) -18.4 E46-E25 N26-N34
      2  422 1671 (N23E62 L= 37)  1058    154        3 ( 342/1639- 478/1699)  21.0 E71-E55 N20-N25
      3 1146 1647 (N34E06 L= 93)   723    106        1 (1110/1598-1177/1686) -26.6 E09-E03 N30-N37
      4 1123 1281 (N10W01 L=100)  1069    426        1 (1085/1108-1159/1481) -66.8 E05-W05 N00-N21
      5 1333 1394 (N21W11 L=111)   228     97        0 (1295/1374-1384/1420)   0.0 W08-W14 N20-N22
      6 1794 1277 (N26W49 L=149)   714     88        2 (1767/1231-1813/1324) -59.0 W46-W51 N23-N28
      7 1334  808 (S13W23 L=123)   517    131        3 (1291/ 773-1375/ 852) -18.4 W20-W26 S14-S10
      8 1463  778 (S11W32 L=131)   601    129        2 (1405/ 740-1523/ 825)  60.3 W28-W35 S14-S07
      9  985  715 (S25W04 L=104)   306     87        0 ( 945/ 678-1023/ 751) -18.4 W01-W07 S27-S23
     10 1497  410 (S31W52 L=152)   388     83        0 (1456/ 395-1537/ 429)   6.3 W47-W56 S31-S30
     11  845  197 (S68W16 L=116)   827    132        2 ( 786/ 160- 901/ 235)  -2.4 W05-W29 S70-S65
  

    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