Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  641 1797 (N38E59 L=142)  1119    118        2 ( 595/1760- 687/1843)  -8.7 E77-E51 N35-N41
      2  630 1673 (N29E47 L=155)  1017    202        1 ( 539/1634- 711/1717)  45.0 E52-E42 N24-N33
      3 1450 1540 (N40W18 L=219)  1215    245        4 (1361/1449-1498/1611) -31.6 W09-W22 N34-N44
      4  342 1539 (N14E65 L=136)   642     97        1 ( 317/1497- 370/1586) -71.6 E67-E62 N11-N17
      5  912 1517 (N25E19 L=182)   335    175        0 ( 855/1451- 968/1572) -24.4 E24-E13 N21-N27
      6 1364 1410 (N29W13 L=214)  1399    501        0 (1218/1207-1451/1544) -36.0 W00-W22 N17-N35
      7  958 1338 (N15E11 L=190)   777    259        5 ( 846/1300-1064/1380)   0.0 E18-E03 N12-N17
      8 1708 1205 (N24W43 L=244)   849    147        1 (1683/1131-1735/1266) -54.5 W40-W45 N20-N28
      9  234 1142 (S12E59 L=143)  2224    286        4 ( 151/1042- 313/1246)  57.3 E70-E51 S18-S04
     10 1903 1012 (N18W69 L=270)   814    112        0 (1891/ 959-1917/1074) -60.3 W66-W72 N14-N21
     11  250 1018 (S19E55 L=147)   400     31        1 ( 230/ 991- 265/1043) -71.6 E57-E53 S20-S17
     12  378  953 (S20E41 L=160)   713    195        3 ( 315/ 885- 423/1036)  57.5 E47-E38 S25-S14
     13 1293  788 (S10W21 L=223)   300     93        1 (1267/ 753-1320/ 825) -36.9 W18-W23 S11-S07
     14 1182  687 (S18W17 L=219)   417    166        1 (1110/ 668-1245/ 718)  26.6 W12-W20 S20-S15
     15  523  617 (S39E27 L=175)  2067    397        6 ( 374/ 560- 684/ 702)  -4.1 E41-E12 S41-S35
     16  296  659 (S41E52 L=149)  1140    116        3 ( 237/ 639- 359/ 696)  11.3 E63-E43 S44-S37
     17  533  464 (S50E26 L=175)   337     66        1 ( 505/ 447- 569/ 486)   8.1 E29-E22 S51-S48
     18 1236  347 (S39W37 L=238)   476     94        0 (1215/ 319-1270/ 390) -20.6 W32-W40 S40-S36
  

    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