Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  958 1695 (N36E24 L=135)  2567    310        4 ( 851/1579-1038/1789) -23.5 E36-E13 N30-N42
      2  529 1713 (N27E58 L=101)   449     66        0 ( 511/1687- 553/1746) -21.8 E63-E53 N25-N29
      3 1593 1604 (N50W30 L=190)  1352    180        3 (1522/1563-1672/1655)  -3.0 W22-W41 N47-N52
      4  402 1329 (N01E47 L=112)  1573    197        3 ( 365/1242- 451/1412) -52.1 E51-E41 S04-N05
      5 1263 1341 (N21W07 L=166)   413    168        0 (1201/1282-1307/1392) -24.0 W01-W10 N18-N22
      6  520  977 (S17E30 L=130)   247     43        2 ( 496/ 952- 546/1001)  14.0 E31-E27 S18-S15
      7 1020  842 (S14W04 L=163)  1052    363        0 ( 873/ 768-1178/ 889)  42.9 E02-W12 S21-S07
      8 1474  803 (S05W32 L=191)   263     96        0 (1439/ 761-1494/ 845) -31.0 W28-W34 S06-S03
      9 1141  612 (S24W18 L=177)  1265    291        1 (1016/ 589-1284/ 653)  26.6 W09-W26 S27-S19
     10  449  578 (S44E33 L=127)   536     81        2 ( 408/ 555- 497/ 612)  21.8 E37-E27 S47-S41
  

    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