Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1791 1244 (N10W54 L= 89)   390     84        0 (1765/1216-1825/1267) -31.0 W51-W57 N08-N11
      2 1391 1133 (N03W23 L= 57)   248     38        0 (1364/1107-1421/1157) -45.0 W20-W24 N01-N04
      3  271  873 (S10E53 L=341)  1071    211        0 ( 255/ 766- 288/ 976)  81.3 E54-E51 S16-S03
      4 1892  900 (S10W66 L=100)   542    139        0 (1852/ 881-1923/ 920)  11.3 W60-W70 S10-S08
      5  586  849 (S13E28 L=  6)   224     86        0 ( 548/ 825- 627/ 864)  21.8 E30-E25 S14-S11
      6  133  819 (S12E71 L=324)   880    108        0 ( 110/ 787- 153/ 849) -45.0 E75-E67 S14-S10
      7  576  666 (S24E31 L=  3)   350     40        0 ( 550/ 639- 602/ 690)   0.0 E33-E29 S25-S22
      8  364  568 (S29E51 L=343)  1448    250        0 ( 311/ 503- 406/ 641)  34.7 E58-E45 S33-S24
      9 1669  520 (S35W50 L= 85)   407     43        0 (1637/ 499-1695/ 545)  45.0 W48-W52 S36-S32
  

    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