Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  674 1750 (N29E50 L=235)  1572    147        5 ( 629/1703- 731/1815) -18.4 E60-E42 N25-N34
      2  437 1665 (N18E65 L=220)   518     35        1 ( 414/1647- 461/1680)   0.0 E69-E61 N16-N19
      3  382 1575 (N11E63 L=223)  2100    244        8 ( 336/1450- 412/1676) -48.8 E73-E56 N02-N18
      4 1693 1568 (N53W40 L=325)   978     97        2 (1655/1528-1732/1598)  -5.2 W34-W46 N49-N56
      5  549 1540 (N13E45 L=240)   362     91        0 ( 528/1500- 577/1577) -78.7 E46-E43 N10-N15
      6 1583 1280 (N26W28 L=313)   280    108        2 (1544/1246-1622/1322)  -8.1 W24-W31 N24-N27
      7 1732 1219 (N27W42 L=327)   497     84        0 (1709/1168-1753/1251) -45.0 W39-W43 N24-N29
      8  276 1096 (S19E53 L=233)   402     70        1 ( 249/1065- 301/1126) -26.6 E55-E49 S20-S17
      9 1744  916 (N10W49 L=334)   430     55        0 (1729/ 880-1758/ 954) -36.9 W46-W51 N08-N12
     10  414  816 (S33E36 L=250)   532    147        0 ( 354/ 779- 476/ 858)  36.0 E41-E30 S36-S28
     11  821  605 (S34E00 L=285)   384     79        2 ( 787/ 582- 858/ 628)  68.2 E01-W01 S36-S31
     12 1580  497 (S17W54 L=340)  3380    390        4 (1493/ 335-1714/ 617) -73.3 W49-W65 S27-S07
  

    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