Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  920 1468 (N15E18 L= 98)   507    159        1 ( 891/1397- 950/1545) -45.0 E22-E14 N11-N19
      2 1045 1436 (N17E10 L=106)   256     77        0 (1012/1422-1084/1447)  11.3 E12-E07 N15-N18
      3 1140 1347 (N14E02 L=114)   236     90        0 (1121/1299-1158/1389) -59.0 E04-W00 N11-N16
      4 1825 1051 (N18W52 L=168)   401     82        0 (1811/1008-1842/1096) -45.0 W50-W55 N16-N20
      5 1339 1042 (N02W17 L=133)   244     53        1 (1313/1022-1368/1064)   0.0 W15-W18 N01-N03
      6 1590  883 (N01W37 L=153)   421     52        1 (1558/ 855-1618/ 912)  11.3 W34-W39 N00-N02
      7  840  582 (S37W02 L=118)   973    212        1 ( 741/ 554- 949/ 612)  33.7 E05-W07 S41-S32
      8 1058  478 (S37W20 L=136)   429     88        1 (1040/ 439-1081/ 533) -63.4 W17-W23 S39-S33
      9 1485  434 (S25W51 L=167)   669    125        0 (1442/ 385-1531/ 493) -83.7 W49-W52 S29-S20
  

    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