Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1617 1577 (N21W58 L=182)   614    154        1 (1593/1515-1630/1624)  39.8 W54-W61 N18-N23
      2  158 1212 (N35E62 L= 62)   592    160        0 ( 124/1188- 183/1232)  -9.5 E65-E59 N33-N36
      3 1792 1144 (S07W54 L=178)  1079    285        1 (1712/1007-1830/1232) -78.7 W49-W56 S15-N00
      4 1670  680 (S31W35 L=159)  1737    359        2 (1530/ 556-1793/ 798)  15.6 W21-W46 S35-S26
      5 1132  622 (S19E03 L=121)   243    109        1 (1080/ 607-1192/ 635) -26.6 E06-W00 S20-S17
      6  820  537 (S17E25 L=100)   456    144        2 ( 784/ 490- 860/ 574) -90.0 E26-E20 S20-S14
      7  716  544 (S14E31 L= 93)   230     39        0 ( 693/ 527- 739/ 558)  14.0 E32-E28 S15-S13
      8  935  476 (S23E20 L=105)   532    131        0 ( 879/ 446-1000/ 528) -58.0 E21-E15 S26-S18
      9  685  443 (S20E37 L= 87)   725    194        4 ( 618/ 370- 730/ 519) -85.2 E38-E34 S25-S13
     10  437  437 (S16E57 L= 67)   397     72        0 ( 412/ 411- 458/ 470) -90.0 E58-E55 S18-S13
     11  526  406 (S20E51 L= 74)  1357    278        6 ( 478/ 363- 586/ 455)   4.8 E57-E45 S23-S16
  

    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