Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1185 1509 (N33W22 L= 60)   613    158        3 (1130/1439-1238/1549)  20.0 W15-W26 N29-N34
      2 1204 1237 (N16W15 L= 53)   266    100        2 (1167/1209-1251/1270) -36.9 W13-W17 N13-N18
      3 1366 1064 (N03W21 L= 60)   232     67        0 (1335/1041-1392/1093)  26.6 W18-W23 N01-N03
      4 1870  981 (S15W60 L= 99)  1396    254        4 (1843/ 886-1899/1062)  74.7 W57-W64 S20-S09
      5  175  992 (N16E64 L=334)   646    113        3 ( 149/ 957- 200/1015)  33.7 E67-E60 N13-N17
      6 1705  971 (S10W43 L= 82)   275     74        1 (1691/ 931-1716/1009)  63.4 W41-W45 S12-S08
      7  232  758 (N01E63 L=335)   768    179        3 ( 189/ 733- 263/ 787)   5.2 E70-E59 N00-N03
      8  309  583 (S10E61 L=337)  1216    185        3 ( 258/ 523- 355/ 655)  60.9 E65-E57 S14-S05
      9 1717  626 (S34W43 L= 81)   461     96        1 (1689/ 600-1750/ 654)  18.4 W40-W46 S35-S31
     10 1592  525 (S38W32 L= 70)   393     65        0 (1556/ 502-1625/ 554)   9.5 W28-W34 S39-S35
     11 1418  390 (S44W17 L= 55)   437     92        0 (1387/ 374-1457/ 410) -23.2 W13-W20 S45-S41
  

    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