Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1011 1760 (N38E27 L=128)   661    134        0 ( 961/1718-1059/1796) -11.3 E31-E21 N36-N40
      2  689 1729 (N27E46 L=109)   390    100        1 ( 671/1682- 710/1775) -33.7 E51-E42 N23-N29
      3  778 1708 (N27E38 L=116)   493    172        1 ( 695/1664- 826/1761)  73.3 E41-E36 N22-N32
      4  555 1680 (N20E53 L=102)   628     93        2 ( 521/1648- 602/1715)   0.0 E58-E48 N18-N21
      5 1141 1691 (N37E15 L=140)   344     70        0 (1118/1665-1178/1711)  68.2 E16-E12 N34-N39
      6 1128 1603 (N31E11 L=143)   234     25        0 (1113/1578-1143/1628) -76.0 E12-E10 N28-N32
      7  490 1407 (N02E43 L=112)  1954    387        6 ( 417/1268- 611/1534)  86.6 E46-E37 S05-N12
      8 1292 1449 (N26W03 L=158)   339     98        1 (1261/1399-1324/1491) -38.7 W00-W06 N23-N28
      9  801 1327 (N05E21 L=134)   347    126        1 ( 776/1266- 829/1382)  90.0 E21-E19 N01-N08
     10 1705  992 (N12W42 L=197)   269    116        0 (1695/ 934-1714/1051) -74.1 W40-W43 N08-N15
     11 1525  373 (S25W59 L=214)  1001    149        2 (1464/ 339-1598/ 410)  53.1 W55-W63 S29-S21
     12 1235  150 (S49W67 L=222)  1119    121        1 (1183/ 117-1287/ 186)   4.1 W55-W84 S50-S47
  

    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