Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  830 1841 (N38E44 L= 51)   663    126        2 ( 799/1808- 872/1870)  -6.3 E48-E39 N36-N40
      2 1087 1685 (N35E18 L= 78)   969    249        3 (1022/1588-1147/1756) -28.1 E25-E10 N29-N38
      3  710 1471 (N10E31 L= 65)   436     91        2 ( 664/1448- 756/1487)  18.4 E34-E28 N08-N12
      4 1703 1117 (N19W39 L=134)   361     66        1 (1688/1076-1719/1152)  90.0 W36-W40 N16-N21
      5  372  815 (S35E38 L= 57)  1192    364        2 ( 252/ 791- 510/ 843)  19.2 E49-E26 S39-S31
      6  829  817 (S23E05 L= 90)   445    177        0 ( 738/ 791- 916/ 839)  26.6 E10-W00 S26-S20
      7 1032  786 (S19W07 L=103)   303     73        0 ( 991/ 763-1072/ 812)  11.3 W04-W09 S19-S18
      8 1765  665 (S02W60 L=156)   583    126        2 (1721/ 645-1806/ 689)  26.6 W56-W64 S04-N00
  

    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