Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  679 1779 (N65E03 L= 39)  1245    229        2 ( 573/1735- 757/1818)  -6.6 E19-W07 N60-N67
      2 1127 1772 (N46W36 L= 78)   329     49        1 (1100/1751-1157/1787)   8.1 W32-W39 N44-N47
      3 1555 1672 (N26W60 L=102)  1306    231        5 (1509/1597-1616/1768)  53.1 W55-W69 N20-N32
      4  518 1719 (N65E29 L= 13)   645     65        2 ( 496/1689- 543/1745)  90.0 E34-E24 N61-N68
      5  657 1659 (N54E09 L= 33)   985    134        4 ( 607/1602- 712/1710)  28.6 E14-E03 N50-N57
      6 1411 1646 (N29W46 L= 88)   554    104        3 (1361/1623-1457/1673)  -5.7 W41-W51 N27-N30
      7  798 1586 (N43W01 L= 43)   358    131        1 ( 760/1550- 833/1629)  23.2 E02-W05 N41-N45
      8  996 1241 (N17W04 L= 45)   277    121        0 ( 975/1189-1013/1293)  71.6 W01-W05 N14-N20
      9  456 1217 (N28E32 L= 10)   426    143        2 ( 404/1156- 491/1256) -49.4 E36-E29 N24-N31
     10 1065 1173 (N11W06 L= 48)   640    180        2 (1028/1085-1099/1255)  53.1 W01-W08 N07-N15
     11  733 1191 (N20E13 L= 29)   406    173        0 ( 648/1162- 801/1211) -29.1 E18-E09 N17-N23
     12  159 1059 (N25E65 L=337)   974    131        5 ( 113/1016- 200/1097) -15.9 E73-E59 N21-N26
     13 1812  886 (S25W51 L= 92)   848    126        2 (1786/ 824-1835/ 953) -77.5 W48-W53 S29-S20
     14 1653  792 (S26W34 L= 76)   341     66        1 (1615/ 777-1685/ 813) -26.6 W31-W37 S27-S24
     15  425  485 (S14E55 L=346)   577     81        3 ( 389/ 456- 462/ 513) -18.4 E58-E52 S15-S11
     16 1675  446 (S53W40 L= 82)   720     96        1 (1638/ 427-1728/ 462) -17.5 W33-W52 S56-S49
     17 1558  362 (S56W26 L= 67)   436     69        2 (1528/ 341-1584/ 386)  31.0 W21-W29 S57-S53
     18 1245  169 (S63E18 L= 24)   667     79        1 (1197/ 152-1279/ 185) -45.0 E21-E13 S66-S58
  

    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