Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  662 1565 (N48E19 L=225)   262     49        1 ( 639/1544- 687/1580) -21.8 E21-E16 N46-N49
      2  815 1553 (N44E05 L=238)   297     37        2 ( 786/1531- 841/1572) -11.3 E07-E02 N42-N45
      3 1333 1371 (N22W27 L=270)  1453    526        1 (1103/1309-1539/1435) -29.1 W12-W39 N13-N29
      4 1290 1221 (N14W20 L=263)   442    180        1 (1207/1193-1370/1245) -29.1 W14-W24 N10-N16
      5  340 1154 (N24E47 L=197)  1459    219        4 ( 308/1063- 391/1240)  84.8 E50-E42 N18-N29
      6 1437 1096 (N03W26 L=270)   234     58        0 (1421/1063-1449/1131)  63.4 W24-W27 N01-N05
      7 1850 1097 (S06W61 L=305)   489     43        0 (1826/1077-1871/1120)  11.3 W58-W64 S07-S05
      8  185  920 (N10E66 L=178)   449     65        0 ( 171/ 888- 204/ 947)  63.4 E67-E62 N08-N12
      9  223  865 (N07E62 L=182)   621     70        1 ( 205/ 835- 245/ 897)  26.6 E64-E58 N05-N08
     10 1098  531 (S25E04 L=239)   310     97        0 (1055/ 495-1142/ 558)  15.9 E07-W00 S27-S23
     11  555  402 (S26E48 L=196)   850    197        1 ( 519/ 321- 621/ 458)  26.6 E58-E38 S32-S22
     12  705  330 (S34E39 L=204)  1046    196        3 ( 637/ 280- 760/ 399)  25.2 E49-E32 S37-S28
     13 1405  289 (S54W18 L=262)   413     63        0 (1377/ 272-1440/ 303) -14.0 W14-W22 S56-S51
  

    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