Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  955 1513 (N38W08 L=126)  1099    287        2 ( 921/1365- 992/1631)  54.8 W01-W14 N28-N45
      2 1183 1595 (N37W27 L=145)   380     67        1 (1153/1564-1222/1615) -21.8 W24-W30 N35-N38
      3 1399 1575 (N30W41 L=159)   583    106        1 (1347/1544-1462/1597) -11.3 W37-W47 N28-N31
      4 1555 1483 (N20W48 L=166)   316    107        0 (1502/1464-1606/1499) -39.8 W43-W50 N17-N22
      5 1014 1277 (N22W05 L=123)   469     86        0 ( 969/1244-1056/1302) -51.3 W02-W07 N19-N24
      6  155 1144 (N28E69 L= 49)  3269    177        5 ( 111/1068- 190/1256)  85.2 E78-E63 N23-N35
      7  234  884 (N12E59 L= 59)  1216    266        1 ( 186/ 820- 279/ 965)  48.8 E64-E52 N08-N16
      8 1770  811 (S24W47 L=165)  2078    337        1 (1682/ 668-1828/ 945)  42.9 W37-W53 S31-S17
      9  293  845 (N09E53 L= 65)   374     66        1 ( 275/ 813- 314/ 880) -63.4 E54-E51 N07-N11
     10 1338  810 (S12W14 L=132)   363    117        3 (1302/ 769-1376/ 870)  45.0 W11-W17 S14-S09
     11  663  768 (S01E27 L= 91)   239    114        0 ( 652/ 711- 671/ 824)  63.4 E28-E25 S04-N02
     12 1478  718 (S21W21 L=139)   297     95        2 (1435/ 682-1518/ 758)  18.4 W18-W24 S22-S18
     13 1339  724 (S17W12 L=130)   505    135        2 (1312/ 668-1394/ 757)  23.2 W09-W16 S20-S14
     14 1022  625 (S16E09 L=109)   305    167        0 (1002/ 559-1063/ 682)  83.7 E10-E05 S21-S12
     15 1200  657 (S18W02 L=121)   281     68        1 (1167/ 627-1233/ 676) -63.4 W00-W04 S20-S16
     16 1090  584 (S20E06 L=112)   344    137        1 (1051/ 548-1112/ 635)  71.6 E07-E03 S22-S16
     17 1739  583 (S39W46 L=164)   540     83        0 (1709/ 560-1781/ 605)   0.0 W41-W51 S40-S37
     18  464  444 (S18E54 L= 64)  1411    214        2 ( 369/ 411- 547/ 506)   0.0 E64-E44 S19-S14
     19 1046  327 (S38E17 L=101)  4147    650        7 ( 787/ 212-1267/ 469)   4.1 E42-W01 S42-S31
  

    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