Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1365 1677 (N47W35 L=340)   344     60        0 (1332/1657-1400/1696)   7.1 W31-W39 N45-N48
      2 1482 1627 (N42W44 L=349)   371     46        1 (1448/1609-1513/1639)   0.0 W40-W47 N40-N42
      3  654 1415 (N29E26 L=280)   230     31        1 ( 633/1390- 672/1437) -71.6 E27-E24 N27-N30
      4 1475 1225 (N15W30 L=336)   426    140        0 (1409/1184-1528/1260) -35.5 W26-W34 N12-N17
      5 1124 1141 (N11W07 L=312)   439     73        2 (1098/1105-1154/1185) -78.7 W04-W08 N08-N13
      6 1800  808 (S13W58 L=  3)   809    127        2 (1771/ 753-1834/ 855) -41.6 W53-W63 S17-S09
      7 1739  781 (S14W51 L=357)   426    105        0 (1721/ 722-1750/ 827)  81.9 W49-W52 S18-S11
      8  388  620 (S21E48 L=257)   385    124        0 ( 361/ 563- 401/ 667)  81.9 E50-E47 S25-S18
      9  640  586 (S23E28 L=278)   353     87        0 ( 612/ 549- 666/ 622)  45.0 E30-E25 S25-S20
     10 1526  481 (S34W38 L=344)   440     97        0 (1489/ 463-1570/ 512) -16.7 W33-W43 S35-S31
  

    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