Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  852 1774 (N58W19 L=340)  2988    421        5 ( 696/1653-1015/1889)   5.6 E02-W49 N53-N61
      2 1611 1680 (N23W69 L= 30)   501     71        0 (1585/1654-1638/1704) -76.0 W66-W72 N21-N25
      3  367 1529 (N56E45 L=276)   571     80        1 ( 323/1506- 399/1555)   0.0 E52-E40 N53-N57
      4  950 1429 (N31W07 L=328)   347    119        1 ( 899/1395- 977/1479)  15.9 W03-W10 N28-N33
      5  257 1260 (N39E53 L=268)   544     68        2 ( 218/1237- 293/1275)  -5.7 E58-E48 N36-N39
      6 1609  877 (S19W31 L=352)   542    126        0 (1592/ 819-1627/ 944)  50.2 W28-W33 S21-S15
      7 1448  751 (S21W17 L=338)   324    138        1 (1413/ 693-1481/ 808)  35.5 W13-W20 S23-S18
      8 1509  493 (S41W17 L=338)  1324    222        3 (1474/ 436-1579/ 578)  60.9 W13-W23 S45-S34
      9 1092  484 (S27E12 L=309)   373    124        0 (1044/ 453-1149/ 522) -60.3 E13-E09 S30-S23
     10  919  452 (S24E23 L=297)   374    153        1 ( 855/ 428- 986/ 487)  -6.3 E27-E18 S25-S22
     11  726  318 (S28E43 L=277)  1381    345        2 ( 592/ 262- 870/ 388) -42.9 E50-E35 S35-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