Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  835 1706 (N35E25 L=196)   707     97        2 ( 798/1662- 882/1738)  51.3 E26-E21 N31-N37
      2  926 1585 (N27E15 L=206)   417    103        0 ( 893/1531- 966/1621)  74.1 E15-E12 N23-N30
      3  742 1204 (N01E19 L=202)   835    257        1 ( 698/1087- 784/1330) -58.4 E23-E15 S05-N08
      4 1596 1217 (N14W33 L=254)   305     51        2 (1569/1196-1623/1241) -14.0 W31-W35 N12-N15
      5 1892 1086 (N12W61 L=282)  3132    389        5 (1840/ 933-1932/1199) -71.6 W55-W66 N03-N18
      6  547  973 (S14E29 L=192)  1019    325        2 ( 437/ 863- 676/1087)  47.0 E36-E21 S21-S06
      7  362  992 (S15E43 L=178)   745    260        0 ( 266/ 936- 454/1041)  21.8 E51-E36 S18-S11
      8 1195  820 (S16W13 L=234)   376    138        1 (1130/ 763-1225/ 862) -39.8 W08-W15 S18-S13
      9  487  654 (S35E33 L=188)   364     43        0 ( 462/ 632- 513/ 683) -45.0 E35-E31 S36-S32
     10  862  214 (S65W02 L=223)   545     63        1 ( 833/ 193- 895/ 241)  -7.1 E01-W07 S67-S62
  

    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