Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  613 1760 (N67E33 L=196)   496     44        0 ( 587/1742- 643/1778)  -5.2 E37-E26 N65-N69
      2  988 1619 (N45W11 L=240)   383     56        1 ( 961/1591-1016/1642)   0.0 W08-W14 N42-N46
      3 1752 1551 (N21W76 L=304)   779     68        0 (1736/1514-1767/1584)  16.7 W71-W83 N19-N22
      4 1619 1415 (N18W50 L=278)  1709    407        0 (1453/1354-1806/1463) -26.6 W38-W66 N09-N24
      5 1336 1320 (N18W25 L=254)   363    159        1 (1315/1247-1355/1377)  60.3 W22-W27 N14-N22
      6 1040 1326 (N24W07 L=235)   228    124        0 ( 985/1313-1100/1340) -15.9 W03-W10 N22-N26
      7  513 1195 (N25E32 L=197)  1431    331        1 ( 414/1077- 558/1302) -68.2 E40-E28 N17-N32
      8 1623 1154 (N03W41 L=269)   609    151        2 (1600/1087-1640/1224)  76.0 W38-W42 S01-N07
      9 1320 1162 (N09W20 L=249)   375    157        0 (1264/1130-1392/1183)  -7.1 W16-W24 N07-N10
     10  318  936 (N11E50 L=178)   344    106        0 ( 295/ 902- 359/ 960)   9.5 E52-E46 N08-N12
     11  658  615 (S13E31 L=198)   236     84        0 ( 634/ 599- 699/ 645) -45.0 E31-E28 S13-S10
     12  911  402 (S31E20 L=208)  1193    333        2 ( 803/ 311-1023/ 514)  25.6 E32-E09 S37-S24
     13  693  398 (S28E36 L=192)   421    161        1 ( 652/ 334- 740/ 457)  31.6 E43-E30 S32-S24
  

    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