Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  237 1178 (N26E52 L=247)   371     62        1 ( 212/1155- 255/1208)  36.9 E54-E49 N24-N27
      2  156 1178 (N27E61 L=238)   704    119        1 ( 122/1153- 184/1207)   8.1 E64-E57 N25-N28
      3 1834 1009 (S15W55 L=353)   525    129        0 (1803/ 966-1863/1061)  26.6 W50-W58 S17-S12
      4  148  910 (N11E66 L=232)   512    121        2 ( 120/ 882- 172/ 943) -53.1 E69-E62 N09-N13
      5  412  824 (N02E42 L=256)   483     66        1 ( 393/ 785- 433/ 863)  63.4 E44-E40 N00-N04
      6 1674  762 (S27W38 L=336)  1365    235        3 (1592/ 684-1765/ 843)  21.8 W30-W46 S30-S23
      7  344  718 (S03E51 L=247)   398     74        1 ( 317/ 688- 371/ 753)  31.0 E53-E47 S04-S01
      8 1605  397 (S53W35 L=334)   659    135        0 (1584/ 370-1635/ 439) -45.0 W31-W40 S56-S48
      9  668  280 (S37E46 L=252)   873    137        1 ( 618/ 262- 726/ 307) -14.0 E50-E41 S38-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