Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1550 1411 (N13W44 L=341)   497     80        1 (1505/1393-1587/1434) -26.6 W40-W46 N11-N14
      2 1532 1082 (S04W31 L=328)   319     85        1 (1491/1049-1561/1100) -59.0 W28-W33 S06-S01
      3 1451 1045 (S03W25 L=322)   226     65        1 (1428/1004-1480/1067)  21.8 W22-W27 S05-S01
      4  638  995 (N15E23 L=274)   326    108        1 ( 589/ 963- 691/1020)   0.0 E27-E19 N13-N15
      5  835  997 (N10E11 L=286)   527    184        0 ( 747/ 978- 933/1016) -31.0 E16-E05 N06-N12
      6  214  800 (N11E63 L=234)  2425    468        5 ( 108/ 750- 329/ 870) -19.2 E75-E51 N07-N16
      7 1333  760 (S16W11 L=308)   555    148        3 (1293/ 730-1389/ 814) -66.8 W08-W13 S19-S12
      8 1197  563 (S24E02 L=295)   496    195        0 (1159/ 485-1222/ 654)  52.1 E06-W01 S28-S18
      9  439  501 (S10E54 L=244)   348     86        1 ( 424/ 467- 455/ 539)  59.0 E56-E51 S12-S07
     10  502  490 (S12E49 L=249)   598    126        2 ( 473/ 452- 545/ 523) -53.1 E52-E45 S13-S09
     11  423  438 (S14E59 L=239)   642    189        2 ( 382/ 415- 456/ 466) -15.9 E63-E55 S15-S12
  

    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