Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1530 1757 (N33W67 L=  3)  1311    149        3 (1456/1735-1600/1776) -30.3 W59-W76 N29-N36
      2  524 1683 (N59E28 L=269)   790     88        2 ( 487/1640- 557/1721)  12.5 E32-E23 N55-N62
      3 1586 1691 (N27W64 L=  1)   455     85        1 (1544/1670-1622/1708) -23.2 W60-W69 N25-N29
      4 1179 1496 (N27W22 L=318)   501    225        1 (1098/1447-1272/1536)   0.0 W15-W28 N25-N28
      5 1004 1184 (N13W02 L=299)   417    164        1 ( 972/1125-1021/1263)  69.4 W00-W04 N09-N17
      6 1806 1103 (S11W54 L=351)   373     54        1 (1782/1077-1829/1129)  33.7 W50-W56 S12-S09
      7  388  565 (S10E54 L=243)   513    100        1 ( 365/ 517- 414/ 607)  26.6 E57-E49 S12-S07
      8  895  375 (S33E26 L=270)  1500    464        3 ( 730/ 275-1058/ 530)  12.7 E41-E10 S37-S26
      9  849  445 (S27E26 L=271)   317    107        1 ( 820/ 417- 879/ 481)  31.0 E28-E22 S28-S25
  

    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