Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1606 1652 (N35W59 L=296)   553     97        0 (1573/1636-1638/1667) -26.6 W55-W61 N33-N36
      2 1230 1370 (N24W18 L=255)   441    133        1 (1177/1319-1294/1418) -53.1 W15-W21 N20-N28
      3  753 1186 (N19E15 L=222)  1809    440        4 ( 553/1093- 933/1283) -31.0 E29-E04 N11-N26
      4 1546 1244 (N12W37 L=274)   306     84        1 (1521/1209-1566/1283)  45.0 W34-W39 N09-N13
      5  565 1043 (N12E29 L=208)   535    112        1 ( 511/ 997- 602/1072) -45.0 E32-E26 N08-N13
      6  241  970 (N09E57 L=180)   393     75        1 ( 221/ 934- 257/ 996) -56.3 E59-E55 N07-N10
      7 1789  803 (S20W52 L=289)  2854    348        3 (1686/ 665-1854/ 939)  39.5 W42-W59 S27-S12
      8 1169  821 (S08W06 L=244)   283     31        0 (1146/ 792-1194/ 847)  45.0 W04-W08 S09-S06
      9 1242  749 (S14W10 L=247)  1026    342        1 (1090/ 650-1346/ 847)  31.0 W00-W17 S18-S08
     10 1060  783 (S09E01 L=236)   242    108        0 (1036/ 738-1073/ 836)  80.5 E02-W00 S11-S05
     11  849  802 (S06E13 L=224)   284    101        0 ( 824/ 760- 898/ 835)  36.9 E15-E09 S07-S03
     12 1223  813 (S09W10 L=247)   241     80        0 (1191/ 796-1277/ 827)   0.0 W07-W13 S10-S08
     13  884  727 (S11E12 L=225)   235     73        3 ( 858/ 704- 916/ 754) -53.1 E13-E10 S12-S08
     14 1583  659 (S25W34 L=271)   413     98        0 (1546/ 605-1601/ 701)  63.4 W31-W35 S28-S22
     15  285  595 (S15E62 L=175)  1059    188        1 ( 265/ 517- 315/ 659)  73.3 E64-E58 S20-S10
     16 1487  620 (S26W26 L=263)   417     82        1 (1450/ 573-1518/ 651)  38.7 W23-W28 S28-S24
     17 1594  586 (S31W35 L=272)   325     63        0 (1559/ 566-1618/ 609)  21.8 W32-W37 S31-S29
  

    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