Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1270 1674 (N35W40 L=148)  1437    226        1 (1173/1634-1369/1698) -31.6 W33-W46 N29-N38
      2  362 1393 (N45E41 L= 67)  1362    170        3 ( 278/1357- 441/1432) -19.2 E54-E31 N41-N49
      3 1842 1212 (S09W61 L=169)  1697    125        5 (1795/1158-1876/1263) -63.4 W56-W66 S12-S04
      4 1510 1124 (S02W31 L=138)   373     71        0 (1494/1084-1531/1170)  76.0 W29-W31 S04-N00
      5 1798 1026 (S17W50 L=158)  1251    249        0 (1712/ 939-1866/1111)  14.0 W40-W60 S19-S14
      6  146 1055 (N28E68 L= 39)  1859    142        2 (  89/ 998- 194/1105) -17.7 E85-E60 N23-N31
      7 1563  646 (S32W24 L=131)  1598    350        3 (1484/ 523-1601/ 793)  60.6 W15-W26 S39-S21
      8  391  746 (N05E47 L= 61)   369     68        0 ( 358/ 725- 428/ 776) -38.7 E49-E44 N03-N07
      9  228  731 (N07E64 L= 43)   483     23        1 ( 207/ 711- 247/ 753)   9.5 E67-E61 N06-N08
     10  881  583 (S15E21 L= 87)   209     55        0 ( 869/ 551- 892/ 614)  63.4 E21-E19 S17-S13
     11  742  447 (S20E35 L= 73)   515    128        0 ( 687/ 405- 789/ 497)  16.7 E39-E29 S21-S17
     12 1364  409 (S42W03 L=111)   297     76        0 (1327/ 392-1407/ 428)  -8.1 W00-W07 S42-S40
  

    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