Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  574 1702 (N30E53 L=346)  1267    213        0 ( 537/1577- 599/1771) -31.8 E63-E42 N21-N34
      2  392 1615 (N21E64 L=335)  1120    178        2 ( 314/1568- 475/1661)  31.6 E71-E57 N16-N24
      3  177 1218 (S07E67 L=332)  1882    156        9 ( 118/1177- 215/1276)  32.5 E78-E61 S09-S02
      4 1344 1188 (N14W16 L= 54)   426    125        1 (1318/1142-1369/1247) -63.4 W13-W16 N11-N17
      5 1015 1081 (N01E02 L= 37)  1047    264        2 ( 940/ 975-1094/1176)  77.9 E04-W01 S06-N08
      6 1700 1133 (N19W42 L= 81)   351     69        0 (1665/1101-1726/1161) -21.8 W38-W44 N17-N19
      7  248 1071 (S14E55 L=344)   416     52        1 ( 223/1044- 272/1101)  53.1 E57-E52 S16-S12
      8 1277  877 (S05W18 L= 56)   618    158        0 (1219/ 834-1339/ 941)  58.0 W14-W20 S08-N00
      9 1763  761 (S01W56 L= 95)  1198    145        2 (1712/ 715-1801/ 826)  80.5 W52-W59 S04-N02
     10  206  781 (S33E61 L=338)   789     84        0 ( 169/ 757- 246/ 808)  -5.2 E66-E55 S34-S31
  

    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