Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  709 1610 (N35E23 L=135)   242     37        0 ( 682/1590- 735/1626)  33.7 E24-E20 N33-N35
      2 1038 1589 (N33W01 L=159)   525    112        1 ( 977/1573-1104/1605)  -7.1 E02-W06 N31-N33
      3 1433 1153 (N05W25 L=183)   664     99        2 (1385/1101-1466/1205) -71.6 W22-W27 N01-N07
      4  234 1009 (S02E55 L=103)   782    135        3 ( 207/ 968- 260/1063)  76.0 E58-E52 S04-N00
      5 1364  915 (S10W21 L=179)   904    326        1 (1266/ 728-1450/1024)  58.6 W15-W26 S21-S03
      6  304  949 (S06E49 L=109)   358     80        0 ( 285/ 898- 326/ 985)  59.0 E50-E46 S09-S04
      7  214  823 (S13E59 L= 99)  4292    736       12 ( 148/ 726- 294/ 945) -52.4 E67-E51 S19-S05
      8 1928  858 (S11W71 L=229)   831    152        3 (1909/ 829-1954/ 886)  38.7 W68-W76 S13-S09
      9 1862  797 (S16W63 L=221)  2206    444       12 (1810/ 722-1913/ 881) -84.3 W56-W70 S20-S10
     10 1588  823 (S15W37 L=195)   300    112        2 (1562/ 776-1617/ 874)  63.4 W35-W38 S17-S11
     11  345  748 (S18E48 L=110)  1415    429        7 ( 268/ 638- 414/ 856)  47.3 E56-E40 S25-S12
     12 1681  773 (S18W45 L=203)   574    151        3 (1620/ 732-1739/ 815)  24.0 W40-W49 S20-S15
     13 1763  764 (S18W53 L=210)   640    127        0 (1733/ 709-1805/ 808)  50.2 W50-W56 S21-S15
     14  164  734 (S18E67 L= 90)  1326    204        3 ( 138/ 687- 194/ 793) -66.8 E71-E63 S21-S14
     15  416  630 (S26E44 L=114)   634    252        1 ( 354/ 581- 510/ 678)  18.4 E50-E35 S29-S23
     16  576  520 (S34E34 L=124)   813    226        0 ( 556/ 433- 609/ 611)  63.4 E36-E29 S40-S28
     17  285  555 (S30E59 L= 99)   911    168        2 ( 225/ 497- 323/ 608) -58.0 E63-E56 S34-S26
     18 1553  552 (S32W39 L=197)   294     67        0 (1516/ 534-1588/ 569)  21.8 W36-W41 S33-S31
  

    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