Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  940 1785 (N49E28 L=252)   485     72        0 ( 912/1761- 975/1806)  71.6 E30-E25 N46-N52
      2  760 1685 (N37E35 L=245)   875    177        0 ( 721/1621- 819/1754) -26.6 E43-E27 N33-N41
      3 1042 1698 (N44E14 L=266)   535     91        0 ( 999/1663-1077/1728) -20.6 E18-E10 N41-N46
      4  518 1625 (N28E52 L=228)   902    195        1 ( 435/1596- 583/1658)   8.1 E58-E44 N24-N30
      5 1208 1355 (N23W06 L=286)   610    250        0 (1114/1287-1287/1430) -24.8 E01-W12 N19-N26
      6 1662 1267 (N25W41 L=321)   330     80        0 (1627/1239-1694/1291) -18.4 W37-W43 N23-N26
      7  405 1170 (S02E43 L=237)   713    159        1 ( 342/1129- 465/1217) -15.3 E49-E37 S03-N00
      8 1698  644 (S12W55 L=335)   521     66        1 (1671/ 609-1719/ 675)  78.7 W53-W57 S14-S09
      9 1089  521 (S30W14 L=294)   402     84        1 (1051/ 502-1138/ 537)  26.6 W11-W17 S32-S29
  

    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