Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1164 1739 (N45W08 L= 98)   321     47        0 (1134/1720-1184/1761)  36.9 W05-W09 N43-N46
      2 1338 1717 (N44W23 L=113)   382     48        2 (1310/1697-1363/1741)  11.3 W20-W25 N42-N46
      3  611 1608 (N32E33 L= 56)   342     89        2 ( 581/1568- 641/1650) -40.6 E37-E30 N29-N35
      4 1396 1485 (N26W24 L=113)   641    217        0 (1331/1412-1472/1567) -45.0 W19-W28 N21-N31
      5 1734 1422 (N25W52 L=141)   381     67        0 (1705/1393-1761/1449) -36.9 W49-W54 N23-N26
      6 1808 1398 (N24W60 L=150)   942    106        1 (1774/1351-1848/1444) -39.8 W56-W64 N21-N26
      7  248 1276 (N10E57 L= 33)   585     98        0 ( 234/1219- 265/1322) -80.5 E59-E54 N06-N12
      8 1839 1092 (N05W58 L=147)  3021    449        5 (1777/ 944-1893/1268) -86.8 W51-W64 S04-N15
      9  214 1237 (N08E60 L= 30)   359     21        0 ( 197/1215- 229/1262)  71.6 E61-E58 N06-N09
     10 1771 1224 (N12W51 L=141)   498     88        4 (1732/1189-1801/1252) -15.9 W47-W54 N10-N14
     11 1719 1145 (N07W46 L=136)   439     99        3 (1694/1088-1742/1183)  71.6 W43-W48 N03-N09
     12  791  877 (S14E14 L= 76)   399    165        0 ( 747/ 818- 854/ 914)  40.6 E16-E09 S17-S11
     13  149  840 (S16E69 L= 20)  1025    132        1 ( 101/ 801- 205/ 868)   5.0 E83-E60 S17-S14
     14  998  778 (S19E01 L= 89)   432     64        2 ( 976/ 740-1026/ 822)  78.7 E02-W01 S21-S16
     15  952  733 (S22E03 L= 86)   198     69        1 ( 913/ 713- 983/ 750) -11.3 E06-E01 S23-S20
     16  835  665 (S27E11 L= 78)   234     95        0 ( 806/ 637- 871/ 710) -51.3 E13-E08 S28-S23
     17 1144  534 (S34W11 L=100)   341     71        1 (1112/ 514-1179/ 558) -21.8 W08-W13 S35-S32
  

    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