Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  676 1733 (N58E09 L=248)  3414    548        5 ( 466/1623- 860/1832)   5.7 E35-W16 N54-N62
      2 1626 1611 (N22W61 L=317)  5818    345       14 (1521/1499-1735/1729)  83.7 W50-W73 N13-N31
      3 1510 1569 (N23W47 L=304)   533    102        1 (1477/1514-1539/1615)  63.4 W45-W49 N19-N25
      4  295 1380 (N40E50 L=207)  1211    143        2 ( 267/1310- 329/1441)  90.0 E54-E46 N35-N44
      5 1697 1327 (N05W51 L=307)   474    114        1 (1651/1292-1740/1360) -39.8 W47-W54 N02-N07
      6  602  793 (S02E30 L=227)   564     95        2 ( 580/ 737- 623/ 842)  71.6 E31-E28 S05-N01
      7  907  761 (S10E12 L=244)   225     70        0 ( 876/ 723- 941/ 788)  18.4 E15-E09 S11-S08
      8 1374  552 (S34W12 L=268)  1922    472        2 (1185/ 445-1526/ 709)  12.5 E02-W25 S37-S26
      9 1518  444 (S46W23 L=280)   887    141        0 (1463/ 407-1577/ 498)  15.3 W17-W28 S48-S42
     10  409  415 (S21E62 L=194)  2119    209        2 ( 343/ 371- 500/ 455)  -2.3 E78-E53 S23-S18
     11 1091  165 (S58E28 L=229)   462     33        1 (1067/ 149-1115/ 187) -76.0 E31-E24 S60-S55
  

    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