Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1707 1500 (N22W60 L=290)   457    154        1 (1663/1485-1752/1508) -11.3 W55-W65 N20-N23
      2 1234 1205 (N14W16 L=246)   595    220        1 (1128/1146-1324/1251) -38.7 W09-W20 N09-N18
      3  184 1194 (N25E62 L=168)  1295    299        3 ( 124/1154- 242/1216)   0.0 E69-E55 N22-N26
      4  455 1062 (N15E37 L=193)  1734    376        5 ( 405/ 875- 530/1213) -77.2 E40-E33 N03-N25
      5  552 1169 (N21E29 L=201)   301     63        0 ( 522/1140- 581/1200) -45.0 E31-E26 N19-N23
      6 1171  965 (N01W08 L=238)   333     89        0 (1120/ 933-1214/1000) -38.7 W05-W10 S01-N03
      7  237  900 (N07E59 L=172)   529    128        3 ( 210/ 873- 277/ 932)  14.0 E62-E54 N05-N09
      8  274  837 (N03E56 L=174)   474     81        2 ( 251/ 806- 297/ 865)  56.3 E58-E53 N01-N04
      9 1032  544 (S24E07 L=223)   288    116        0 (1017/ 483-1044/ 587)  76.0 E08-E05 S28-S20
     10 1429  258 (S58W25 L=255)   562     82        0 (1391/ 246-1471/ 273)  -5.7 W20-W30 S59-S56
  

    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