Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1301 1689 (N37W35 L=239)   292     31        0 (1280/1667-1328/1706)  21.8 W32-W37 N35-N38
      2  625 1603 (N44E20 L=184)   605    127        0 ( 568/1558- 683/1631)  16.7 E25-E15 N41-N46
      3 1692 1492 (N18W57 L=261)   365     77        0 (1654/1465-1722/1517) -38.7 W53-W59 N15-N19
      4  852 1448 (N29E04 L=200)  1596    410        4 ( 712/1349-1043/1510)   9.9 E14-W09 N24-N31
      5  451 1370 (N30E34 L=170)  1768    435        5 ( 263/1329- 670/1422) -16.1 E55-E17 N24-N36
      6 1923 1104 (S09W70 L=274)  1202     98        2 (1886/1056-1947/1158)  63.4 W64-W74 S11-S05
      7 1869 1037 (S12W61 L=265)  1089    172        4 (1824/ 968-1914/1085) -63.4 W55-W67 S16-S08
      8 1824  911 (S19W55 L=259)  1152    250        0 (1782/ 810-1868/ 998)  74.7 W50-W60 S24-S13
      9 1521  889 (S15W29 L=233)   295     77        0 (1478/ 868-1554/ 908)   9.5 W25-W31 S16-S14
     10  993  856 (S09E04 L=200)   564    165        0 ( 955/ 787-1060/ 896)  29.7 E07-W00 S12-S06
     11 1668  811 (S22W40 L=244)   360     89        0 (1654/ 777-1683/ 863) -80.5 W38-W41 S24-S18
     12 1744  751 (S27W47 L=251)   668     87        0 (1695/ 727-1787/ 786)   6.3 W42-W51 S28-S25
     13  559  553 (S20E39 L=165)  1471    448        2 ( 426/ 443- 679/ 685)  22.2 E53-E26 S25-S13
     14  282  630 (S11E60 L=144)   554     59        0 ( 258/ 607- 308/ 658) -56.3 E63-E57 S12-S09
     15  334  458 (S23E66 L=138)  3272    168        3 ( 286/ 360- 385/ 527)  61.4 E75-E59 S29-S18
     16  882  473 (S31E20 L=185)   241     53        0 ( 864/ 446- 900/ 505)  53.1 E21-E18 S32-S28
  

    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