Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1304 1716 (N44W17 L= 61)   614    116        2 (1261/1682-1350/1742)  26.6 W12-W20 N40-N46
      2 1535 1632 (N39W35 L= 79)  2043    343        5 (1396/1559-1651/1687)  -5.2 W23-W46 N34-N42
      3 1114 1606 (N33W02 L= 46)   753    155        1 (1038/1547-1161/1645) -26.6 E03-W05 N28-N35
      4 1668 1540 (N34W45 L= 89)   426     91        2 (1642/1511-1687/1566)  78.7 W42-W47 N31-N36
      5 1763 1495 (N32W53 L= 98)  1259    223        3 (1713/1459-1820/1524)   5.2 W48-W59 N29-N34
      6  282 1353 (N12E55 L=350)   871    219        1 ( 228/1287- 331/1394) -28.6 E60-E48 N07-N14
      7 1040 1283 (N11E01 L= 44)   230     68        1 (1014/1261-1077/1306) -33.7 E02-W01 N09-N12
      8 1410 1100 (N03W23 L= 67)   760    249        1 (1336/ 984-1463/1229) -53.1 W17-W27 S03-N09
      9 1864  870 (S06W62 L=107)   547    183        3 (1826/ 848-1911/ 889)  16.7 W58-W69 S07-S03
     10 1748  631 (S21W56 L=101)   532    160        1 (1692/ 605-1801/ 661)  24.0 W51-W61 S23-S18
     11  532  571 (S36E33 L= 11)   616    115        1 ( 479/ 545- 589/ 602) -20.6 E37-E28 S37-S33
  

    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