Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1276 1429 (N29W08 L=228)  2122    647        4 (1091/1179-1442/1585) -32.3 E07-W23 N15-N35
      2 1150 1209 (N12W04 L=223)   700    234        1 (1054/1176-1256/1249)   4.8 E02-W10 N10-N13
      3 1763 1078 (N16W50 L=270)   661     97        2 (1743/1032-1785/1122)  63.4 W48-W52 N13-N19
      4  252  942 (S20E54 L=165)  2887    403        5 ( 132/ 809- 366/1083)  31.7 E78-E44 S30-S09
      5 1794  994 (N12W55 L=274)   385     67        1 (1770/ 972-1821/1022) -21.8 W51-W57 N10-N13
      6 1777  875 (N05W55 L=274)   659    127        2 (1750/ 806-1808/ 923)  74.1 W53-W57 N00-N07
      7 1098  724 (S18W10 L=230)   334     48        1 (1069/ 700-1126/ 747) -26.6 W08-W12 S19-S16
      8 1624  666 (S11W48 L=267)   984    148        2 (1547/ 635-1679/ 699)  26.6 W42-W52 S13-S08
      9  946  596 (S29W04 L=223)   286    126        0 ( 892/ 575-1008/ 610)  29.7 W00-W07 S31-S26
     10 1600  539 (S19W52 L=271)   532    165        0 (1545/ 497-1663/ 581)  45.0 W48-W56 S22-S15
     11  781  293 (S57E01 L=218)  1360    264        4 ( 644/ 269- 890/ 328)  21.4 E19-W10 S64-S52
  

    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