Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1624 1541 (N21W56 L=133)  2131    419        4 (1475/1482-1765/1597) -32.5 W46-W68 N13-N28
      2  326 1424 (N46E48 L= 30)   728    169        1 ( 292/1388- 361/1468) -33.7 E52-E43 N43-N49
      3 1299 1296 (N16W23 L=100)   538    171        1 (1220/1261-1342/1350)   6.3 W17-W26 N14-N18
      4  333 1135 (N27E45 L= 32)   664    113        0 ( 294/1100- 369/1175) -39.8 E48-E42 N24-N29
      5  918  617 (S15E15 L= 62)   895    197        1 ( 868/ 524- 963/ 697)  45.0 E20-E10 S19-S10
      6 1720  603 (S38W41 L=118)   655    188        1 (1690/ 576-1764/ 653) -41.6 W37-W46 S41-S32
      7  326  557 (S08E61 L= 16)  1982    218        3 ( 259/ 506- 401/ 619) -35.0 E68-E54 S12-S03
      8 1710  537 (S43W41 L=118)  1207    301        2 (1622/ 515-1780/ 565) -22.8 W31-W51 S47-S39
      9  660  325 (S29E44 L= 33)   699    236        1 ( 624/ 279- 679/ 410)  48.0 E49-E39 S32-S22
  

    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