Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  970 1616 (N33E06 L=293)   608    103        1 ( 942/1560- 996/1667)  90.0 E08-E04 N29-N36
      2  559 1625 (N33E38 L=262)   348     65        0 ( 525/1609- 596/1641)   0.0 E40-E34 N31-N34
      3 1402 1484 (N26W24 L=323)   702    193        0 (1321/1445-1498/1526) -14.0 W18-W30 N23-N28
      4  826 1433 (N20E14 L=285)   257     73        0 ( 807/1400- 841/1472)  90.0 E15-E12 N17-N22
      5 1578 1312 (N16W35 L=334)   838    290        0 (1451/1232-1692/1369) -22.4 W26-W43 N11-N18
      6  227 1334 (N14E60 L=240)   545    112        0 ( 174/1302- 259/1368)  20.6 E66-E56 N11-N15
      7 1877 1213 (N12W62 L=  2)  1914    237        0 (1850/1131-1935/1326) -71.6 W59-W70 N07-N19
      8  167 1244 (N08E65 L=234)   598     56        0 ( 133/1224- 202/1266) -11.3 E70-E60 N07-N09
      9 1588 1052 (N00W35 L=335)   567    160        0 (1566/ 974-1623/1128) -63.4 W33-W38 S04-N04
     10 1176  900 (S11W10 L=309)   303    130        0 (1142/ 845-1221/ 957)  66.8 W07-W12 S14-S07
     11 1842  824 (S11W60 L=360)  1777    158        6 (1788/ 758-1905/ 880)  26.6 W55-W67 S15-S07
     12 1726  843 (S11W48 L=348)   263     59        2 (1710/ 814-1746/ 871)  56.3 W46-W49 S12-S09
     13  983  672 (S26E01 L=298)  1155    335        0 ( 912/ 559-1121/ 750)  47.5 E05-W07 S32-S20
     14 1651  591 (S27W49 L=348)  1095    159        2 (1567/ 563-1721/ 624)  13.0 W41-W54 S28-S24
     15  258  572 (S33E65 L=234)   645     93        0 ( 224/ 540- 299/ 592)   4.1 E72-E58 S34-S31
     16 1030  537 (S34W02 L=302)   380     79        0 ( 982/ 520-1076/ 557)  -9.5 E01-W05 S35-S32
     17 1477  373 (S43W44 L=343)  2507    393        3 (1297/ 348-1658/ 415)   7.7 W26-W63 S45-S38
     18  664  341 (S50E31 L=268)  1771    294        3 ( 506/ 307- 800/ 386)   3.4 E51-E16 S52-S46
  

    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