Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1549 1662 (N26W57 L= 81)  3654    481        5 (1370/1547-1767/1724) -31.8 W42-W73 N13-N31
      2 1558 1531 (N18W48 L= 73)   695    201        0 (1470/1482-1643/1564) -45.0 W43-W53 N13-N22
      3  189 1256 (N35E61 L=324)  1093    120        3 ( 136/1222- 252/1290) -13.4 E72-E51 N31-N37
      4 1832 1153 (S09W57 L= 82)  2042    303        5 (1760/1091-1888/1223)  -3.8 W49-W64 S12-S05
      5  269 1155 (N27E49 L=335)   419    115        0 ( 247/1111- 299/1211) -54.5 E51-E46 N23-N30
      6 1383 1115 (N00W22 L= 47)   423    149        0 (1325/1073-1445/1139)  -7.1 W17-W25 S02-N01
      7 1089  920 (S04W01 L= 26)  1732    468        5 ( 957/ 745-1210/1054) -73.3 E03-W04 S16-N05
      8 1406  981 (S08W21 L= 45)   415     88        1 (1359/ 946-1443/1009) -11.3 W18-W23 S09-S05
      9  517  958 (N10E31 L=353)   552    172        1 ( 435/ 938- 596/ 989) -31.0 E36-E26 N07-N13
     10 1369  909 (S11W17 L= 42)   232    103        3 (1341/ 869-1402/ 947) -90.0 W15-W18 S14-S09
     11 1572  728 (S27W28 L= 52)  1468    271        4 (1486/ 682-1688/ 805) -23.6 W21-W37 S30-S20
     12  727  736 (S06E24 L=  1)   247     46        0 ( 705/ 712- 745/ 757) -71.6 E24-E22 S07-S04
     13  856  724 (S10E16 L=  8)   376     94        2 ( 825/ 686- 882/ 754) -90.0 E17-E14 S12-S08
     14 1364  590 (S30W10 L= 35)  1826    403        4 (1303/ 431-1405/ 751)  82.6 W04-W13 S41-S18
     15  694  641 (S11E28 L=356)   497    112        1 ( 643/ 610- 748/ 680) -45.0 E30-E25 S13-S07
     16  631  449 (S21E40 L=344)  3344    612        9 ( 432/ 307- 802/ 624)  15.5 E71-E24 S26-S13
     17  617  595 (S12E35 L=350)   478     99        2 ( 579/ 558- 666/ 623)   8.1 E37-E30 S13-S10
     18 1269  372 (S43E02 L= 23)   969    154        4 (1208/ 317-1325/ 437) -68.2 E04-W02 S47-S36
     19 1202  297 (S47E10 L= 15)  1274    135        3 (1152/ 257-1266/ 342)   0.0 E15-E04 S49-S42
  

    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