Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1369 1760 (N37W49 L= 73)   353    109        1 (1347/1729-1400/1801)  63.4 W47-W50 N34-N40
      2 1689 1529 (N15W59 L= 84)  1038    343        2 (1573/1505-1783/1568) -22.8 W51-W70 N11-N19
      3 1495 1427 (N14W39 L= 64)   401    140        1 (1418/1376-1555/1473) -56.3 W35-W41 N09-N18
      4 1796 1372 (N03W62 L= 87)   914    234        2 (1753/1326-1844/1420)   3.6 W55-W71 N01-N05
      5 1455 1097 (S03W26 L= 51)   405    100        0 (1423/1046-1488/1156) -81.9 W25-W27 S06-N01
      6 1903 1056 (S17W62 L= 87)   744    119        2 (1891/ 996-1921/1109)  50.2 W60-W65 S20-S14
      7  683  860 (N02E23 L=  1)   536    152        1 ( 608/ 823- 749/ 896)  -5.7 E28-E18 N00-N03
      8  819  734 (S08E19 L=  6)   253     35        1 ( 797/ 707- 841/ 754)  33.7 E20-E16 S09-S06
      9  416  482 (S14E55 L=329)  1342    255        3 ( 313/ 425- 502/ 540) -36.9 E64-E49 S19-S09
  

    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