Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  781 1620 (N38E25 L=308)  5751    796       10 ( 470/1516-1131/1771)  -7.0 E79-W03 N32-N46
      2 1236 1432 (N29W11 L=345)   478    120        0 (1194/1379-1270/1485) -54.5 W08-W13 N25-N32
      3  279 1425 (N20E63 L=270)   684     67        1 ( 263/1389- 304/1460) -53.1 E66-E59 N17-N21
      4 1046 1382 (N24E02 L=332)   329    129        0 (1029/1323-1064/1445) -82.9 E02-W00 N20-N28
      5  940 1365 (N22E08 L=325)   558    145        1 ( 913/1299- 978/1439) -58.0 E10-E05 N18-N26
      6 1653 1371 (N28W44 L= 18)   410     83        1 (1613/1343-1697/1395) -14.0 W40-W48 N26-N29
      7  207 1203 (N06E63 L=271)  2717    258        1 ( 181/1076- 236/1320) -70.3 E66-E59 S02-N12
      8 1826 1239 (N20W60 L= 34)   757     97        2 (1794/1188-1844/1277) -90.0 W56-W63 N17-N22
      9 1794 1163 (N15W55 L= 29)   709    121        0 (1748/1141-1849/1196)  -9.5 W50-W62 N13-N17
     10 1801  956 (N03W56 L= 29)   692    112        0 (1763/ 907-1851/ 984)  24.0 W52-W61 N00-N04
     11 1416  933 (S01W25 L=359)   226     26        0 (1397/ 910-1438/ 956) -33.7 W23-W26 S02-N00
     12  469  865 (S12E35 L=298)   237     84        0 ( 445/ 837- 508/ 884)  21.8 E37-E32 S14-S11
     13  385  809 (S17E42 L=291)   471    126        0 ( 343/ 774- 438/ 851) -29.7 E45-E38 S18-S14
     14 1117  674 (S19W09 L=342)   927    207        0 (1020/ 637-1216/ 706)  22.6 W02-W14 S22-S16
     15  829  541 (S31E10 L=324)   862    173        2 ( 742/ 511- 912/ 570)  18.4 E15-E03 S33-S28
     16  554  394 (S45E36 L=297)   464     57        1 ( 522/ 376- 589/ 414)  -9.5 E39-E32 S46-S43
  

    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