Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  929 1391 (N11E14 L=243)  2010    357        3 ( 829/1211-1083/1610) -34.7 E27-E01 N04-N22
      2 1141 1537 (N26E06 L=251)   282    118        0 (1123/1478-1157/1590) -50.2 E09-E03 N22-N28
      3  332 1448 (N03E58 L=200)  1294    319        1 ( 222/1405- 434/1482)  13.4 E69-E48 N00-N05
      4  403 1411 (N02E49 L=208)   282     46        0 ( 391/1382- 414/1442) -56.3 E50-E47 N00-N03
      5 1486 1335 (N23W20 L=277)   525     97        2 (1462/1277-1503/1381) -63.4 W17-W21 N19-N25
      6  291 1195 (S13E52 L=205)   597    141        3 ( 246/1170- 344/1217)   6.3 E56-E47 S14-S10
      7  201 1108 (S19E59 L=199)  1942    308        9 ( 152/1030- 255/1194) -52.1 E64-E52 S22-S13
      8 1482 1060 (N06W25 L=283)  1188    231        1 (1441/ 938-1522/1179) -47.7 W20-W30 N00-N12
      9  628  866 (S25E20 L=237)   756    251        1 ( 513/ 835- 740/ 920)  33.7 E28-E13 S29-S19
     10  257  788 (S38E50 L=207)  1617    444        1 ( 166/ 729- 386/ 822)  22.2 E60-E38 S42-S33
     11 1636  413 (S22W62 L=320)  1663    167        1 (1569/ 372-1688/ 483)  60.9 W58-W65 S26-S17
  

    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