Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1171 1797 (N49W27 L=336)   323     63        1 (1148/1775-1193/1820)  21.8 W24-W29 N47-N51
      2 1698 1596 (N27W60 L=  9)   476     82        0 (1672/1582-1722/1614)   0.0 W56-W62 N25-N27
      3  194 1309 (N27E59 L=250)  3575    466        9 ( 137/1217- 245/1407) -85.2 E65-E53 N20-N32
      4  344 1283 (N23E43 L=266)   663    219        0 ( 258/1242- 408/1355)  17.1 E51-E38 N21-N26
      5 1758 1160 (S01W51 L=360)   716    203        5 (1718/1112-1816/1189)  15.3 W46-W57 S03-N01
      6 1221 1101 (N02W12 L=321)   193     42        0 (1205/1070-1234/1127)  71.6 W11-W13 N00-N03
      7  199  853 (S01E61 L=248)   469     63        1 ( 178/ 823- 221/ 887) -90.0 E63-E58 S02-N01
      8  325  734 (S09E51 L=258)  1364    279        1 ( 246/ 664- 398/ 818) -42.0 E57-E45 S13-S03
      9  923  661 (S21E11 L=298)   991    239        2 ( 838/ 569- 978/ 765)  39.8 E17-E05 S25-S15
     10 1205  681 (S23W07 L=316)   310     80        1 (1168/ 650-1242/ 728)  38.7 W04-W10 S24-S20
     11 1885  695 (S31W63 L= 12)   844     60        0 (1863/ 656-1908/ 727)  59.0 W61-W65 S33-S28
     12  248  666 (S13E61 L=248)  1286    168        2 ( 205/ 634- 297/ 711) -24.0 E67-E56 S14-S09
     13 1457  632 (S30W24 L=333)   403     78        1 (1421/ 597-1495/ 667)  26.6 W21-W27 S31-S27
     14  305  601 (S17E58 L=251)  1320    199        4 ( 255/ 556- 341/ 652) -39.8 E62-E53 S20-S14
     15 1685  602 (S35W44 L=353)   530    104        0 (1667/ 562-1704/ 644) -71.6 W42-W46 S37-S31
     16  358  524 (S23E56 L=253)   798    176        2 ( 332/ 458- 378/ 591)  53.1 E58-E52 S26-S18
  

    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