Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  567 1727 (N32E56 L=343)  1222    125        2 ( 532/1665- 593/1775) -26.6 E65-E49 N27-N34
      2  389 1616 (N21E64 L=335)  1581    189        5 ( 313/1566- 471/1671)  36.0 E71-E58 N16-N25
      3  171 1219 (S07E68 L=331)  2355    158        9 ( 115/1175- 212/1293)  48.8 E78-E61 S10-S02
      4 1141 1274 (N15W02 L= 41)   221    107        0 (1094/1257-1195/1287)  45.0 W00-W04 N12-N16
      5 1336 1186 (N14W15 L= 54)   405    116        1 (1311/1143-1361/1249) -71.6 W13-W16 N11-N17
      6 1009 1088 (N01E02 L= 37)   995    248        3 ( 936/ 988-1084/1187)  72.9 E04-W00 S05-N08
      7 1693 1136 (N19W42 L= 81)   348     73        0 (1660/1105-1718/1166) -21.8 W38-W44 N17-N19
      8  243 1072 (S14E55 L=344)   428     56        0 ( 218/1043- 261/1102)  45.0 E57-E53 S16-S12
      9 1270  882 (S05W17 L= 56)   544    161        0 (1212/ 837-1338/ 945)  58.0 W14-W20 S08-N00
     10 1760  766 (S01W56 L= 95)  1211    128        2 (1717/ 719-1800/ 829)  71.6 W51-W59 S04-N02
     11  211  788 (S33E60 L=339)  1052    120        2 ( 164/ 760- 260/ 816)   3.8 E67-E52 S34-S30
  

    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