Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  319 1164 (S10E49 L=263)  1108    243        3 ( 238/1097- 383/1257)  50.7 E55-E45 S14-S03
      2  668 1033 (S12E21 L=290)   251     53        0 ( 641/1015- 704/1060) -11.3 E23-E18 S13-S11
      3  460  952 (S20E35 L=277)  2143    495        3 ( 217/ 823- 754/1054)  24.2 E55-E15 S30-S12
      4  137  956 (S23E65 L=247)  1957    244        1 ( 113/ 864- 173/1033)  78.7 E67-E61 S28-S18
      5 1829  939 (N06W56 L=  8)   852    125        5 (1800/ 896-1857/ 984) -51.3 W52-W59 N03-N08
      6 1028  747 (S23W06 L=318)   565    265        0 ( 962/ 647-1094/ 839)  77.0 W03-W08 S29-S16
      7  496  638 (S39E31 L=281)   304     71        0 ( 465/ 618- 534/ 655)  18.4 E33-E27 S41-S37
      8  305  580 (S46E51 L=261)   437     28        0 ( 292/ 554- 322/ 608) -90.0 E52-E48 S48-S44
      9 1615  578 (S19W50 L=  1)   320    105        0 (1561/ 566-1649/ 599)  26.6 W45-W52 S20-S17
     10 1625  517 (S22W53 L=  5)   496    112        1 (1577/ 496-1672/ 550)  45.0 W50-W56 S24-S19
     11 1663  366 (S29W66 L= 18)   641    121        0 (1644/ 341-1680/ 393) -18.4 W62-W69 S29-S27
     12 1584  341 (S33W61 L= 12)  2171    283        3 (1522/ 306-1646/ 379)   4.4 W53-W66 S35-S30
  

    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