Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1165 1687 (N39W32 L=143)   399     51        1 (1135/1663-1189/1713)   0.0 W29-W34 N37-N40
      2 1278 1645 (N33W37 L=148)   599    135        0 (1228/1595-1318/1675) -71.6 W33-W39 N28-N35
      3  817 1626 (N46W05 L=115)   306     50        1 ( 782/1610- 850/1641) -18.4 W01-W07 N44-N47
      4 1331 1567 (N27W36 L=147)   330     61        1 (1301/1545-1352/1598)  11.3 W33-W38 N24-N28
      5  868 1352 (N27E01 L=110)   192     78        1 ( 826/1333- 907/1370) -45.0 E03-W01 N24-N28
      6  430 1288 (N34E34 L= 77)   563    148        2 ( 385/1225- 484/1345) -52.1 E37-E29 N29-N38
      7  779 1275 (N25E08 L=103)   290    105        0 ( 755/1226- 806/1325) -81.9 E09-E07 N21-N28
      8  417 1155 (N26E36 L= 75)   322     99        4 ( 380/1130- 462/1178)   8.1 E39-E32 N25-N27
      9  237 1062 (N25E54 L= 56)   726    103        2 ( 204/1024- 284/1099)  12.5 E58-E49 N22-N27
     10  138  999 (N22E69 L= 42)  1019     79        4 ( 106/ 967- 173/1042) -29.1 E75-E63 N20-N25
     11  190  798 (N10E66 L= 45)   517    107        0 ( 169/ 750- 210/ 847) -90.0 E67-E63 N07-N13
     12  941  768 (S08E11 L=100)   632    171        0 ( 901/ 686- 989/ 840)  37.9 E15-E06 S11-S04
     13 1736  641 (S39W42 L=152)   496     64        1 (1704/ 612-1765/ 669)   8.1 W38-W45 S40-S36
     14 1104  509 (S27E09 L=101)   268     64        0 (1070/ 483-1132/ 530)  21.8 E12-E07 S28-S26
     15  982  460 (S27E18 L= 92)   415     95        2 ( 940/ 438-1037/ 484) -45.0 E20-E15 S29-S24
     16  735  394 (S25E37 L= 73)   505    127        1 ( 673/ 379- 800/ 409) -18.4 E41-E32 S26-S23
     17 1135  191 (S54E25 L= 86)   456     36        1 (1113/ 169-1157/ 213)  14.0 E29-E21 S55-S51
     18 1040  126 (S57E43 L= 68)   606     50        1 (1014/ 108-1061/ 149)  13.0 E49-E36 S59-S54
  

    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