Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  506 1715 (N22E62 L=191)   637     80        2 ( 479/1668- 523/1750) -33.7 E66-E57 N18-N24
      2  993 1699 (N34E25 L=229)   420    124        0 ( 972/1638-1011/1749) -45.0 E29-E21 N29-N37
      3  451 1663 (N17E63 L=191)  2004    158        9 ( 410/1603- 506/1728) -45.0 E72-E54 N12-N22
      4  776 1578 (N19E32 L=222)   431    104        2 ( 739/1527- 806/1618) -63.4 E34-E29 N16-N22
      5  458 1514 (N08E51 L=202)  1243    163        3 ( 384/1474- 536/1540)  31.0 E56-E45 N04-N11
      6 1778 1285 (N34W45 L=299)   508     93        1 (1765/1231-1796/1327)  76.0 W43-W48 N29-N37
      7  603 1132 (S10E27 L=226)   286     92        1 ( 562/1108- 643/1160)  -9.5 E30-E24 S11-S09
      8  400 1095 (S17E41 L=212)   454    110        2 ( 377/1052- 418/1154) -45.0 E44-E39 S19-S14
      9 1842 1022 (N20W56 L=310)  1500    321        2 (1801/ 863-1884/1150) -45.0 W49-W65 N12-N26
     10 1741 1041 (N17W45 L=299)   300     88        0 (1716/1006-1759/1082)  68.2 W43-W46 N14-N19
     11  343  946 (S27E44 L=210)   812    231        0 ( 258/ 908- 451/1019)  30.5 E52-E35 S31-S20
     12  211  950 (S30E59 L=195)   630    104        2 ( 176/ 917- 246/ 986) -18.4 E64-E54 S31-S28
     13 1823  881 (N11W58 L=311)   403     85        0 (1810/ 835-1837/ 923) -78.7 W56-W59 N08-N14
     14  835  660 (S31E00 L=253)   283     75        1 ( 821/ 625- 857/ 700) -45.0 E02-W01 S32-S29
     15  864  500 (S40W07 L=261)   292     58        0 ( 846/ 464- 886/ 527) -78.7 W05-W08 S43-S38
     16 1127  442 (S35W27 L=281)   241     46        1 (1107/ 419-1149/ 468) -76.0 W25-W29 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