Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  410 1392 (N30E46 L=194)  1438    313        3 ( 332/1307- 536/1447)  19.4 E54-E34 N24-N34
      2  971 1394 (N28E01 L=239)   372    200        0 ( 877/1360-1041/1416) -15.3 E08-W03 N25-N30
      3  255 1377 (N29E65 L=176)   680     88        2 ( 220/1349- 294/1395)  11.3 E70-E59 N27-N30
      4  282 1318 (N25E58 L=182)   406     85        2 ( 256/1284- 317/1350) -24.0 E61-E52 N23-N27
      5  455 1321 (N25E40 L=200)   686    114        1 ( 397/1296- 522/1340)  10.3 E45-E34 N23-N26
      6  196 1186 (N17E66 L=175)  1543    243        5 ( 176/1082- 220/1292) -77.0 E71-E61 N10-N23
      7 1797  914 (S08W56 L=296)   619    103        1 (1776/ 873-1819/ 970)  51.3 W53-W58 S10-S05
      8 1302  934 (S02W17 L=257)   430    101        1 (1255/ 886-1335/ 964)  26.6 W13-W19 S05-N00
      9 1749  857 (S11W51 L=291)   481     55        3 (1721/ 833-1775/ 878) -21.8 W48-W53 S13-S09
     10  980  594 (S22E05 L=235)   265    145        0 ( 918/ 559-1029/ 637)  26.6 E09-E01 S24-S19
     11  752  570 (S23E21 L=219)   657    137        3 ( 727/ 505- 777/ 632)  71.6 E23-E18 S27-S18
     12  521  551 (S24E40 L=201)   342    112        0 ( 473/ 527- 560/ 584)  18.4 E44-E35 S25-S21
     13  823  249 (S50E27 L=214)   421     82        2 ( 778/ 234- 865/ 267)  11.3 E32-E22 S51-S48
  

    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