Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1652 1428 (N39W32 L=165)  1634    267        4 (1544/1339-1750/1483)  -5.4 W21-W43 N35-N41
      2 1490 1359 (N28W19 L=152)   692    129        1 (1432/1308-1532/1404) -23.2 W15-W23 N25-N29
      3 1348 1339 (N22W10 L=143)   292    101        1 (1319/1291-1368/1390)  82.9 W09-W11 N18-N26
      4  402 1337 (S02E48 L= 84)  1344    169        6 ( 338/1283- 461/1380) -14.0 E55-E43 S05-N00
      5 1238 1247 (N14W06 L=138)   304    121        0 (1178/1227-1285/1274)  50.2 W02-W08 N10-N16
      6 1812 1186 (N28W51 L=183)   533    122        2 (1767/1136-1846/1219) -12.5 W45-W55 N26-N30
      7 1754 1081 (N20W46 L=178)   526     69        4 (1728/1045-1777/1113)  90.0 W43-W48 N17-N22
      8  141 1064 (S23E71 L= 61)   623     87        1 ( 108/1029- 173/1096)  17.1 E78-E65 S25-S20
      9 1523  871 (N01W33 L=166)   452    173        0 (1505/ 798-1541/ 963) -71.6 W31-W35 S03-N06
     10 1226  876 (S07W15 L=148)  1791    543        6 ( 971/ 782-1363/ 943)  18.4 W01-W25 S15-S01
     11 1675  900 (N07W43 L=176)   264     75        1 (1637/ 881-1713/ 910)  18.4 W40-W46 N05-N08
     12 1238  785 (S12W19 L=151)   339    192        1 (1190/ 738-1272/ 836) -38.7 W15-W21 S14-S09
     13  536  727 (S35E24 L=108)   556    162        1 ( 496/ 659- 583/ 790)  65.6 E26-E21 S40-S29
     14 1620  622 (S09W50 L=183)   291     98        0 (1576/ 602-1672/ 639)  33.7 W47-W53 S11-S07
     15 1729  526 (S11W68 L=200)   950    132        2 (1684/ 478-1764/ 590)  90.0 W64-W71 S13-S06
     16 1486  327 (S30W62 L=194)  1730    187        5 (1426/ 281-1564/ 384)   0.0 W51-W77 S31-S27
  

    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