Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1366 1301 (N15W27 L=  7)  2064    520        3 (1258/1163-1544/1489) -78.7 W23-W35 N03-N28
      2  716 1295 (N30E14 L=326)  1198    400        5 ( 574/1227- 884/1380) -36.9 E23-E03 N22-N37
      3 1821 1234 (S02W61 L= 40)   559     50        2 (1801/1203-1842/1258)  36.9 W58-W63 S03-N00
      4 1013 1158 (N15W02 L=342)   288    208        3 ( 951/1088-1084/1214)  18.4 E02-W07 N11-N17
      5 1661 1095 (S05W41 L= 21)   267     69        1 (1637/1076-1694/1123) -53.1 W39-W43 S06-S02
      6  366  991 (N18E44 L=295)  2094    399        3 ( 199/ 952- 520/1060)  -5.5 E62-E31 N15-N21
      7  653  916 (N09E25 L=315)   942    213        2 ( 581/ 875- 744/ 968)   0.0 E29-E17 N06-N10
      8 1180  579 (S23E01 L=338)   261     88        0 (1137/ 561-1219/ 597) -45.0 E03-W01 S24-S20
      9  388  522 (S11E57 L=282)   405     17        1 ( 373/ 499- 407/ 544)  56.3 E59-E54 S12-S09
     10 1218  452 (S33E02 L=337)   464     99        0 (1193/ 416-1237/ 510)  74.1 E03-W00 S35-S28
     11  476  400 (S20E57 L=283)   598    116        0 ( 439/ 348- 515/ 444)  82.9 E58-E55 S24-S16
  

    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