Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1114 1840 (N52E27 L=154)  1811    167        4 (1031/1802-1181/1875)  26.6 E33-E18 N48-N56
      2  814 1790 (N37E43 L=138)   427     55        2 ( 790/1764- 839/1810)  -9.5 E46-E39 N35-N39
      3  444 1629 (N17E60 L=121)   939    128        2 ( 414/1592- 482/1682) -49.4 E65-E54 N14-N21
      4  492 1512 (N11E48 L=133)   432     95        1 ( 461/1476- 534/1538)  68.2 E49-E45 N08-N13
      5 1712 1416 (N40W40 L=221)   743    132        1 (1648/1400-1774/1432)  22.4 W32-W49 N36-N43
      6 1142 1403 (N21E03 L=178)   330    106        1 (1093/1368-1186/1426)   0.0 E06-W00 N19-N21
      7 1162 1126 (N05W05 L=186)   211     69        1 (1130/1109-1195/1154)  14.0 W03-W07 N04-N06
      8 1733 1075 (N18W44 L=225)   547     45        1 (1706/1041-1754/1107) -63.4 W41-W46 N16-N20
      9  374  991 (S22E42 L=140)   340    123        0 ( 354/ 948- 389/1060) -81.9 E42-E40 S24-S17
     10  887  944 (S12E06 L=175)   471    167        1 ( 858/ 856- 913/1022) -84.3 E07-E04 S17-S07
     11  285  980 (S24E50 L=131)   494    125        2 ( 244/ 943- 344/1011)  21.8 E54-E44 S26-S21
     12  496  924 (S23E30 L=151)   311     57        0 ( 465/ 901- 516/ 949)  33.7 E32-E28 S24-S21
     13 1159  828 (S12W13 L=194)   242     68        0 (1130/ 810-1192/ 847)  33.7 W10-W14 S12-S10
     14  418  737 (S36E35 L=146)   554    121        1 ( 375/ 693- 457/ 784)  48.8 E39-E31 S40-S32
     15 1475  428 (S25W50 L=232)   646    163        0 (1419/ 374-1531/ 483)  83.7 W49-W51 S29-S20
     16  680  412 (S52E06 L=175)   283     26        1 ( 658/ 394- 701/ 433)   0.0 E08-E03 S54-S50
  

    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