Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1708 1450 (N33W52 L= 74)   671    110        4 (1678/1418-1744/1476)  -9.5 W48-W55 N30-N34
      2  272 1353 (N17E59 L=323)   530    151        1 ( 241/1325- 307/1380) -18.4 E63-E54 N15-N18
      3  271 1124 (N03E55 L=328)   320     79        0 ( 250/1098- 296/1156) -45.0 E57-E51 N01-N05
      4  431  797 (S16E39 L=343)  1869    389        4 ( 391/ 637- 507/ 927)  67.2 E42-E33 S26-S07
      5 1737  756 (S11W53 L= 75)   585    144        0 (1725/ 681-1755/ 822)  90.0 W51-W54 S15-S06
      6 1816  744 (S12W62 L= 84)   804    183        1 (1778/ 712-1842/ 787)  53.1 W58-W66 S13-S09
      7 1652  701 (S15W46 L= 68)   456    197        0 (1585/ 657-1702/ 732) -15.3 W40-W51 S17-S13
      8 1337  538 (S27W25 L= 47)   283     99        0 (1292/ 524-1387/ 553)   0.0 W21-W28 S28-S26
  

    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