Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  972 1923 (N65W46 L= 10)   950    104        2 ( 936/1879-1008/1964)   9.2 W34-W66 N61-N67
      2  884 1860 (N62W25 L=349)   829    105        2 ( 847/1816- 924/1892) -39.8 W20-W32 N57-N65
      3 1102 1748 (N45W31 L=355)  1290    215        2 (1007/1701-1205/1786)  -3.0 W21-W40 N42-N46
      4 1316 1735 (N37W45 L=  9)  1129    140        2 (1245/1699-1398/1759) -31.0 W40-W50 N33-N40
      5 1668 1621 (N20W66 L= 30)   725     41        0 (1636/1600-1701/1649) -36.9 W61-W70 N18-N21
      6  804 1594 (N43W01 L=325)   380     99        0 ( 782/1544- 831/1634)  56.3 W00-W04 N39-N45
      7 1754 1496 (N11W65 L= 29)   574     90        0 (1726/1449-1780/1536)  71.6 W62-W66 N07-N13
      8 1578 1110 (S05W34 L=358)  2125    465        2 (1522/ 886-1612/1312)  58.2 W26-W40 S15-N06
      9 1888 1159 (S11W64 L= 28)   502     41        0 (1865/1135-1913/1180) -31.0 W60-W67 S12-S09
     10  137 1109 (N27E67 L=257)  1624    144        2 (  86/1054- 185/1162) -22.8 E80-E59 N23-N31
     11  758 1085 (N13E14 L=310)   702    252        1 ( 679/1032- 857/1160) -57.5 E17-E09 N08-N19
     12  385 1062 (N20E39 L=285)   762    107        1 ( 356/1006- 408/1112) -90.0 E41-E37 N16-N23
     13  304  958 (N16E48 L=276)   582    177        0 ( 245/ 902- 366/ 997)  14.0 E54-E42 N13-N16
     14  342  867 (N10E47 L=277)   338     55        0 ( 314/ 836- 360/ 887) -63.4 E48-E44 N07-N11
     15  253  729 (N03E59 L=265)   920    167        1 ( 200/ 677- 294/ 769) -45.0 E64-E55 N00-N06
     16 1272  610 (S26W05 L=329)   513    190        0 (1213/ 544-1340/ 674)  26.6 W00-W10 S28-S23
     17 1755  572 (S44W46 L= 10)  1078    159        0 (1689/ 526-1823/ 622)   6.7 W37-W55 S45-S40
     18 1382  376 (S46W07 L=331)   326     58        1 (1347/ 357-1409/ 401)   0.0 W03-W09 S47-S43
     19  738  253 (S36E46 L=278)   831    115        0 ( 693/ 216- 780/ 302)  14.9 E53-E38 S38-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