Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1177 1596 (N35W29 L=156)  1835    441        5 (1003/1541-1384/1632) -17.8 W15-W44 N28-N38
      2  956 1468 (N34W09 L=136)   491    111        2 ( 938/1417- 976/1521)  66.8 W06-W11 N30-N37
      3  962 1368 (N28W06 L=133)   233     44        1 ( 933/1351- 990/1386) -18.4 W03-W07 N26-N29
      4  326 1013 (N22E47 L= 80)   286     87        0 ( 296/ 984- 365/1042) -39.8 E49-E43 N19-N24
      5 1854  950 (S22W56 L=183)   785    122        0 (1815/ 923-1897/ 977) -15.3 W51-W62 S23-S18
      6 1746  910 (S20W44 L=171)   371    108        0 (1734/ 861-1760/ 965)  81.9 W42-W45 S23-S16
      7 1530  789 (S20W24 L=151)   232     53        1 (1509/ 762-1547/ 819)  56.3 W22-W25 S21-S18
      8 1610  694 (S29W29 L=156)   596    107        1 (1553/ 669-1663/ 730) -12.5 W24-W33 S30-S26
      9 1838  671 (S40W56 L=183)   582    101        0 (1805/ 633-1872/ 709)  11.3 W51-W62 S42-S38
     10 1520  579 (S34W20 L=147)   850    123        3 (1477/ 526-1554/ 638)  35.5 W16-W23 S36-S30
     11  402  384 (S18E68 L= 59)  1027     79        3 ( 369/ 350- 427/ 426)  80.5 E72-E64 S20-S14
  

    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