Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  992 1845 (N59W38 L=290)   595     85        1 ( 961/1809-1022/1877)   7.6 W31-W46 N57-N61
      2 1384 1725 (N36W52 L=305)   608     73        3 (1340/1697-1417/1752)  -6.3 W48-W57 N34-N37
      3 1642 1372 (N09W49 L=301)   338     67        0 (1604/1361-1682/1385) -26.6 W45-W51 N07-N10
      4 1799 1223 (S04W56 L=309)   957    148        0 (1773/1150-1827/1285)  32.0 W51-W61 S06-N00
      5  863 1237 (N23E04 L=248)   288     86        1 ( 821/1209- 900/1278) -63.4 E06-E02 N20-N26
      6  829 1149 (N19E08 L=244)   620    254        2 ( 722/1105- 943/1179) -15.9 E15-E01 N15-N20
      7  865  982 (N08E10 L=243)   237    162        0 ( 835/ 912- 904/1049)  40.6 E13-E06 N05-N11
      8 1400  997 (S04W21 L=273)   282     80        0 (1371/ 957-1416/1026)  56.3 W19-W21 S05-S02
      9  274  878 (N14E54 L=198)   816    155        1 ( 243/ 801- 305/ 943) -77.5 E56-E51 N09-N18
     10  181  862 (N14E66 L=187)   897    134        2 ( 151/ 814- 210/ 906) -80.5 E69-E61 N11-N17
     11 1740  801 (S26W41 L=294)   619     99        2 (1699/ 759-1775/ 840)  23.2 W37-W45 S27-S23
     12  153  769 (N09E74 L=179)   641     49        3 ( 135/ 738- 174/ 799) -53.1 E76-E69 N06-N10
     13  319  502 (S09E66 L=187)   783     87        3 ( 295/ 455- 337/ 548)  50.2 E68-E62 S12-S06
     14  462  423 (S16E56 L=197)  3287    332       12 ( 370/ 324- 554/ 528) -48.0 E65-E49 S23-S09
     15  979  466 (S24E18 L=235)   686    158        1 ( 917/ 414-1040/ 521) -73.3 E19-E15 S29-S19
     16  653  471 (S17E38 L=215)   524     77        1 ( 620/ 429- 685/ 505)  80.5 E39-E35 S19-S13
     17 1258  343 (S42E05 L=248)  2951    612        5 (1041/ 204-1500/ 502)  11.3 E28-W17 S47-S36
     18  617  387 (S21E45 L=207)   602    119        2 ( 578/ 347- 639/ 447)  36.9 E49-E40 S24-S17
  

    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