Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1557 1684 (N25W59 L=182)  1933    272        5 (1470/1637-1631/1711) -36.0 W52-W65 N20-N29
      2  738 1328 (N30E09 L=114)   807    356        1 ( 574/1298- 908/1358) -22.2 E20-W02 N25-N35
      3 1708 1255 (S02W49 L=172)  1057    244        5 (1650/1211-1785/1302) -24.4 W43-W55 S05-N01
      4 1014 1186 (N14W04 L=127)   386     71        1 ( 983/1153-1049/1223) -53.1 W02-W05 N12-N16
      5 1647 1164 (S04W41 L=164)   381    150        1 (1612/1098-1694/1204)  18.4 W36-W45 S07-S02
      6 1847 1080 (S16W56 L=179)   576    108        1 (1812/1058-1874/1114)   0.0 W51-W59 S16-S14
      7 1611  946 (S15W33 L=156)   570    186        1 (1596/ 866-1640/1041)  68.2 W30-W35 S20-S10
      8 1332  781 (S16W11 L=134)   389     89        2 (1290/ 749-1378/ 817) -63.4 W08-W13 S19-S13
      9  942  351 (S32E27 L= 97)   449    143        3 ( 918/ 290- 964/ 400)  36.9 E31-E23 S36-S28
     10  499  360 (S20E56 L= 67)   436     84        0 ( 484/ 328- 518/ 388)  33.7 E59-E53 S22-S18
     11 1183  268 (S47E16 L=107)  1243    219        4 (1129/ 216-1218/ 329)  51.3 E21-E11 S51-S41
  

    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