Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  648 1789 (N35E53 L= 95)  1022    208        1 ( 607/1720- 689/1845) -63.4 E57-E49 N28-N40
      2 1032 1738 (N42E22 L=125)   681    100        0 (1003/1692-1054/1776) -40.6 E26-E19 N38-N44
      3  434 1549 (N14E54 L= 93)  1962    327        6 ( 357/1473- 509/1602)  55.0 E60-E47 N08-N18
      4  599 1505 (N15E39 L=108)   437     91        1 ( 574/1453- 625/1542) -63.4 E41-E36 N11-N17
      5  756 1449 (N15E27 L=120)   224    126        0 ( 687/1418- 789/1490)  60.3 E30-E25 N11-N18
      6  848 1391 (N14E19 L=128)   274     94        1 ( 820/1345- 884/1424)  90.0 E21-E17 N11-N16
      7  717 1312 (N06E25 L=122)   313     97        0 ( 697/1263- 731/1365) -74.1 E26-E23 N02-N09
      8 1117 1086 (N02W04 L=151)   445    100        1 (1076/1037-1157/1139)  74.1 W01-W05 S01-N06
      9  902  841 (S17E03 L=145)   320     79        1 ( 874/ 801- 930/ 878) -45.0 E05-W00 S18-S14
     10 1723  530 (S13W64 L=212)  1567    214        1 (1654/ 471-1794/ 597)  48.0 W60-W70 S17-S07
     11  662  403 (S53E10 L=138)   368     83        1 ( 611/ 385- 702/ 422)   0.0 E15-E05 S53-S51
     12 1112  375 (S40W27 L=174)  1003    241        1 (1007/ 349-1240/ 403)  12.5 W18-W36 S42-S37
  

    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