Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1415 1689 (N32W49 L=  1)  1279    336        1 (1300/1643-1563/1723) -25.3 W37-W58 N25-N36
      2  404 1629 (N58E45 L=268)  1084    134        1 ( 347/1583- 458/1671)  21.0 E54-E37 N55-N61
      3 1100 1304 (N17W11 L=324)   244    109        2 (1043/1271-1137/1325) -45.0 W08-W12 N14-N19
      4 1187 1250 (N12W15 L=327)   213    105        0 (1144/1219-1236/1275)   8.1 W11-W18 N11-N12
      5  762 1085 (N12E14 L=299)   279    117        1 ( 731/1035- 778/1145)  45.0 E16-E11 N10-N15
      6 1641 1055 (S09W38 L=351)   321     25        1 (1618/1029-1663/1075) -45.0 W36-W40 S10-S07
      7 1771  626 (S40W49 L=  2)  1028    174        0 (1711/ 583-1838/ 691)  11.3 W41-W58 S43-S37
      8  401  527 (S12E55 L=258)   419     64        1 ( 378/ 499- 426/ 560) -78.7 E56-E53 S14-S09
      9  745  328 (S32E40 L=273)  1559    387        3 ( 622/ 232- 879/ 459)  15.9 E58-E23 S37-S26
     10 1002  275 (S43E25 L=287)   638    182        0 ( 955/ 213-1037/ 362)  32.0 E33-E17 S48-S37
  

    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