Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  764 1649 (N40E30 L=124)   644     79        1 ( 733/1613- 800/1680) -15.9 E33-E26 N36-N41
      2 1396 1578 (N42W24 L=178)   843    125        0 (1344/1536-1454/1614) -26.6 W19-W29 N39-N44
      3  349 1535 (N26E64 L= 90)  1559    155        3 ( 327/1466- 377/1605) -26.6 E76-E57 N21-N30
      4 1618 1495 (N38W44 L=198)   500    118        0 (1577/1440-1645/1530) -71.6 W40-W46 N34-N40
      5  797 1402 (N23E19 L=135)   627    203        1 ( 771/1315- 836/1498)  85.6 E20-E18 N16-N29
      6 1631 1246 (N21W40 L=194)   289     97        1 (1612/1195-1650/1291) -80.5 W39-W41 N17-N23
      7  446  869 (S14E37 L=117)  2429    437        4 ( 320/ 697- 562/1010)  46.3 E50-E29 S26-S04
      8 1061  798 (S12W04 L=158)   192     76        1 (1041/ 766-1079/ 829)  76.0 W03-W05 S14-S10
      9 1805  731 (S10W63 L=217)   720    103        3 (1787/ 669-1823/ 781) -74.1 W61-W66 S13-S06
     10 1436  719 (S14W30 L=184)   516    143        1 (1377/ 697-1504/ 737)   0.0 W26-W35 S15-S12
     11 1710  672 (S14W54 L=208)   563    146        2 (1666/ 641-1768/ 697)  18.4 W49-W59 S16-S12
     12  297  650 (S30E56 L= 98)   438     72        1 ( 263/ 632- 326/ 670)  12.5 E61-E52 S32-S28
     13 1392  536 (S26W32 L=186)   396    115        0 (1350/ 502-1435/ 567)  51.3 W29-W34 S28-S23
     14 1736  504 (S25W69 L=223)   949     64        2 (1711/ 473-1767/ 534) -29.7 W63-W75 S27-S23
     15  557  419 (S45E34 L=120)   687    127        1 ( 520/ 383- 617/ 468)  29.7 E40-E26 S48-S39
     16 1401  408 (S35W38 L=192)  1049    141        1 (1328/ 385-1468/ 435)  14.0 W31-W43 S37-S33
  

    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