Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1367 1660 (N35W45 L=296)  1156    197        0 (1282/1626-1446/1683) -21.0 W37-W51 N31-N38
      2 1611 1661 (N27W67 L=318)  1110    126        3 (1544/1643-1673/1683) -12.5 W59-W78 N24-N29
      3  340 1467 (N49E51 L=201)   567     59        1 ( 322/1437- 363/1503) -32.0 E54-E46 N46-N51
      4 1591 1356 (N12W44 L=296)   493    126        1 (1535/1308-1647/1389) -54.5 W41-W46 N07-N15
      5  276 1326 (N40E55 L=197)   793    113        1 ( 244/1288- 315/1371)  35.5 E58-E50 N37-N42
      6 1255 1251 (N15W19 L=270)   569    185        1 (1219/1159-1294/1312)  83.7 W16-W22 N09-N18
      7  707  946 (N09E20 L=231)   277    114        0 ( 676/ 889- 732/ 993) -81.9 E21-E19 N05-N12
      8  837  821 (S01E15 L=236)   935    326        0 ( 762/ 686- 951/ 885)  21.8 E22-E07 S07-N02
      9 1702  763 (S25W39 L=291)   310     51        0 (1678/ 736-1722/ 788) -59.0 W37-W41 S27-S22
     10 1258  747 (S14W07 L=259)   218     46        0 (1236/ 722-1288/ 767)   0.0 W05-W09 S15-S12
     11  606  618 (S08E35 L=217)   581    145        0 ( 562/ 576- 650/ 674) -74.1 E36-E32 S11-S04
     12 1602  501 (S41W30 L=281)  2078    364        4 (1477/ 403-1710/ 631)  25.5 W17-W40 S46-S34
     13  286  588 (S06E64 L=188)  1329     91        2 ( 245/ 556- 327/ 624) -29.7 E68-E59 S08-S03
     14  388  552 (S09E55 L=197)   716     93        1 ( 373/ 501- 413/ 599)  60.3 E57-E52 S12-S05
     15  441  501 (S13E52 L=199)   550     20        0 ( 417/ 475- 463/ 528)  56.3 E55-E50 S14-S11
     16 1124  477 (S28E07 L=245)   312     94        0 (1078/ 465-1161/ 502) -11.3 E09-E04 S29-S26
     17  753  434 (S23E32 L=219)   792    125        3 ( 719/ 390- 787/ 484)  82.9 E34-E29 S26-S18
     18  470  410 (S20E56 L=196)   882    138        3 ( 434/ 354- 496/ 477)  26.6 E65-E48 S23-S15
  

    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