Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1369 1514 (N31W33 L=321)  1181    386        1 (1201/1438-1523/1561) -27.6 W20-W41 N23-N35
      2  238 1379 (N35E66 L=222)  1083    116        2 ( 185/1357- 297/1397)  -2.4 E80-E56 N33-N35
      3  296 1314 (N30E55 L=233)   347     45        1 ( 270/1289- 314/1340) -45.0 E58-E52 N28-N32
      4  786 1316 (N26E12 L=276)   408    166        0 ( 704/1301- 863/1329)  -9.5 E18-E06 N25-N27
      5 1422 1235 (N13W28 L=316)   284     80        2 (1395/1191-1441/1266) -78.7 W26-W29 N10-N15
      6  611 1131 (N17E25 L=263)   207     57        1 ( 585/1108- 629/1157)  33.7 E27-E24 N15-N18
      7 1884  864 (S18W66 L=354)  4927    400        5 (1817/ 667-1946/1004)  69.1 W60-W78 S30-S09
      8 1599  892 (S10W36 L=324)   740    194        4 (1548/ 823-1651/ 943)  66.8 W31-W39 S14-S06
      9 1641  841 (S14W39 L=327)   340     96        0 (1626/ 800-1662/ 879)  59.0 W37-W40 S16-S11
     10  294  815 (S01E54 L=234)   572    144        2 ( 269/ 769- 317/ 876) -80.5 E56-E51 S03-N03
     11  251  675 (S09E63 L=225)  5856    441       15 ( 169/ 523- 322/ 796)  62.2 E79-E55 S20-S01
     12 1620  719 (S22W37 L=325)   279     77        2 (1583/ 701-1652/ 739) -26.6 W33-W39 S22-S19
     13 1034  576 (S22E05 L=283)  1432    384        5 ( 976/ 471-1149/ 709) -90.0 E08-W02 S30-S13
     14  462  467 (S24E50 L=238)   599    169        3 ( 426/ 419- 491/ 523) -90.0 E52-E47 S28-S20
     15  810  447 (S29E23 L=265)   531    162        1 ( 775/ 393- 839/ 502)  41.2 E27-E19 S32-S25
     16 1253  359 (S42W08 L=296)   884    136        1 (1192/ 328-1312/ 400)  10.3 W02-W13 S44-S39
  

    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