Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1735 1549 (N21W71 L=137)   573     35        1 (1714/1522-1752/1574)  14.0 W67-W76 N19-N21
      2 1679 1373 (N13W53 L=118)   307     24        1 (1659/1351-1697/1392)  14.0 W50-W54 N11-N13
      3  184 1324 (N37E72 L=353)   993     68        2 ( 156/1289- 215/1358) -11.9 E85-E64 N34-N38
      4 1470 1262 (N12W33 L= 98)   267     60        1 (1433/1245-1499/1284) -26.6 W30-W34 N10-N12
      5 1812 1012 (S11W54 L=120)   486     68        0 (1797/ 977-1833/1051) -90.0 W52-W56 S13-S08
      6 1906 1000 (S15W67 L=133)   923     85        1 (1879/ 968-1928/1049)  31.0 W62-W71 S16-S12
      7 1609  972 (S08W36 L=101)   541    137        3 (1549/ 928-1646/1019) -54.5 W32-W38 S11-S04
      8 1033  882 (S02E02 L= 64)   296    130        1 ( 972/ 851-1084/ 920) -51.3 E05-W00 S04-N01
      9  150  812 (N05E74 L=352)   673     66        2 ( 135/ 774- 163/ 846)  68.2 E77-E71 N02-N07
     10  211  699 (S02E69 L=357)   513     36        0 ( 193/ 675- 228/ 721) -71.6 E71-E66 S03-N00
     11 1075  613 (S19E05 L= 61)   410    134        0 (1057/ 549-1101/ 676)  83.7 E05-E02 S23-S14
     12 1726  637 (S33W45 L=111)   732    105        2 (1680/ 620-1777/ 657) -10.3 W39-W51 S35-S31
     13  269  587 (S09E67 L=359)   719     66        1 ( 236/ 562- 300/ 616) -36.9 E70-E62 S11-S07
     14  795  529 (S19E24 L= 42)   412    169        0 ( 762/ 459- 818/ 605)  48.4 E28-E20 S23-S14
     15 1345  548 (S29W12 L= 78)   459     91        0 (1318/ 504-1373/ 590)  45.0 W09-W14 S31-S26
     16  366  504 (S16E60 L=  6)   412     91        1 ( 352/ 461- 389/ 543)  29.1 E64-E55 S18-S13
     17  588  358 (S28E48 L= 18)   518    100        2 ( 536/ 344- 639/ 378)  -5.7 E53-E43 S29-S26
     18  645  255 (S38E52 L= 13)   575     56        1 ( 619/ 235- 667/ 283)  18.4 E56-E47 S39-S35
  

    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