Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  636 1721 (N61E09 L=126)   435     90        0 ( 605/1689- 674/1755)  21.8 E12-E02 N59-N63
      2  563 1512 (N46E19 L=115)   324     39        1 ( 543/1484- 588/1530)  33.7 E21-E16 N44-N47
      3 1004 1496 (N32W12 L=147)   342     93        0 ( 958/1474-1043/1519) -51.3 W09-W14 N29-N34
      4 1232 1435 (N22W24 L=159)   229     87        1 (1191/1416-1273/1454)  -9.5 W20-W26 N20-N22
      5  471 1151 (N25E32 L=103)   439    155        1 ( 455/1078- 494/1213)  60.3 E33-E29 N20-N28
      6 1823  977 (S20W52 L=187)   626    134        1 (1791/ 914-1873/1018)  16.7 W48-W59 S22-S18
      7 1763  755 (S32W44 L=179)   699    180        0 (1697/ 710-1830/ 822)  11.3 W37-W52 S33-S29
      8  955  681 (S13E13 L=122)   241     94        2 ( 935/ 645- 977/ 732)  59.0 E14-E10 S14-S09
      9  595  599 (S09E37 L= 98)   675    160        1 ( 571/ 520- 635/ 656)  41.2 E41-E33 S13-S05
     10 1540  552 (S38W22 L=156)   442    128        0 (1490/ 523-1602/ 586)   0.0 W17-W27 S38-S36
     11  636  347 (S26E47 L= 87)  1913    258        5 ( 502/ 308- 745/ 423) -12.8 E63-E40 S29-S20
  

    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