Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1218 1807 (N52W13 L= 89)   362     49        1 (1193/1786-1247/1826) -21.8 W11-W16 N49-N53
      2 1344 1757 (N48W25 L=100)   445     65        1 (1330/1722-1365/1788)  49.4 W22-W28 N44-N51
      3 1483 1655 (N40W34 L=110)   396     51        1 (1465/1628-1506/1688)  56.3 W31-W36 N37-N43
      4  357 1527 (N26E54 L= 21)   526     78        1 ( 329/1501- 396/1553)   9.5 E57-E49 N24-N27
      5 1569 1478 (N27W37 L=112)   355     66        1 (1539/1448-1595/1500) -26.6 W34-W38 N25-N28
      6  287 1305 (N11E53 L= 22)  1225    298        2 ( 252/1173- 336/1418)  86.2 E56-E51 N03-N18
      7  465 1356 (N14E39 L= 37)   284     82        1 ( 440/1312- 481/1382) -68.2 E41-E37 N11-N16
      8 1854 1091 (N05W60 L=135)   836    108        3 (1827/1039-1883/1150) -63.4 W56-W63 N02-N08
      9 1009  862 (S14E00 L= 75)   266    173        2 ( 975/ 806-1073/ 895)  56.3 E02-W03 S17-S11
     10  215  845 (S16E60 L= 16)   796    126        2 ( 158/ 819- 274/ 877)   7.1 E68-E52 S17-S14
     11 1725  806 (S13W50 L=125)   289     37        0 (1708/ 778-1742/ 828)  90.0 W48-W51 S14-S11
     12 1195  746 (S20W12 L= 88)   620    150        2 (1138/ 694-1251/ 815)  74.1 W08-W15 S23-S15
     13  401  773 (S21E42 L= 33)   300     51        1 ( 381/ 750- 421/ 806) -76.0 E44-E40 S22-S18
     14  191  759 (S22E66 L= 10)   530     81        1 ( 170/ 732- 219/ 792)  63.4 E68-E62 S23-S19
     15  933  731 (S22E05 L= 71)   402    158        0 ( 864/ 704-1004/ 759) -12.5 E09-W00 S23-S20
     16 1665  719 (S19W47 L=122)   501     89        0 (1625/ 692-1704/ 754)  26.6 W43-W49 S20-S16
     17 1329  518 (S34W25 L=101)   489     61        0 (1289/ 492-1367/ 549) -23.2 W21-W28 S35-S32
     18 1017  161 (S68W09 L= 84)   553     66        0 ( 989/ 136-1046/ 187) -26.6 W03-W15 S71-S65
  

    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