Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1513 1573 (N21W50 L=209)  1222    206        2 (1433/1519-1610/1602) -10.0 W41-W58 N17-N23
      2 1139 1542 (N31W24 L=183)   712    224        0 (1056/1516-1268/1564) -23.2 W18-W32 N27-N33
      3  632 1382 (N37E15 L=144)   314     88        0 ( 581/1368- 677/1393) -20.6 E19-E11 N35-N38
      4 1735 1322 (N01W54 L=213)   906    147        2 (1699/1247-1770/1386)  66.8 W52-W56 S03-N04
      5  135  898 (N19E72 L= 86)   683     83        2 ( 115/ 849- 149/ 933) -90.0 E75-E70 N15-N20
      6 1260  750 (S16W06 L=165)   645    262        1 (1163/ 671-1337/ 841)  21.0 E01-W12 S18-S12
      7  813  662 (S09E22 L=137)   921    314        1 ( 734/ 539- 880/ 781) -90.0 E23-E20 S17-N00
      8 1752  655 (S39W42 L=201)   609     54        1 (1720/ 630-1779/ 681)   0.0 W38-W45 S40-S37
      9 1602  582 (S38W26 L=184)   305     46        1 (1578/ 561-1622/ 605)  14.0 W23-W27 S39-S36
     10 1303  544 (S30W03 L=162)   424     97        1 (1251/ 527-1357/ 567) -18.4 W00-W06 S31-S27
     11  700  486 (S16E35 L=123)   303     65        0 ( 679/ 457- 727/ 514) -76.0 E36-E34 S18-S14
     12  961  335 (S33E27 L=132)   794    130        1 ( 920/ 276- 988/ 383)  41.2 E31-E23 S37-S29
  

    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