Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  745 1299 (N23E17 L=196)   213     72        2 ( 719/1270- 771/1332) -53.1 E19-E15 N21-N25
      2  735  935 (N01E19 L=194)   214     37        0 ( 714/ 919- 758/ 959) -26.6 E19-E16 N00-N02
      3 1251  896 (S04W13 L=226)   650    206        3 (1153/ 863-1331/ 941) -24.4 W07-W18 S06-S01
      4 1562  835 (S11W34 L=247)   376    111        0 (1520/ 819-1613/ 862) -23.2 W31-W38 S12-S08
      5  215  772 (S09E64 L=148)  1032     92        3 ( 186/ 732- 248/ 806) -45.0 E68-E60 S11-S06
      6 1675  705 (S21W45 L=258)   526    115        0 (1640/ 653-1702/ 750)  56.3 W42-W47 S23-S17
      7  659  613 (S20E27 L=185)   690    181        1 ( 610/ 547- 725/ 675)  26.6 E32-E22 S23-S15
      8  577  509 (S27E36 L=176)   541     75        1 ( 558/ 468- 608/ 546)  71.6 E39-E33 S29-S23
      9  569  404 (S35E42 L=171)   358     51        1 ( 552/ 377- 596/ 423)  26.6 E44-E38 S36-S33
  

    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