Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1702 1498 (N22W60 L=290)   472    166        1 (1655/1482-1749/1508) -18.4 W55-W64 N20-N23
      2  622 1329 (N30E23 L=207)   271     82        1 ( 581/1315- 663/1344) -18.4 E26-E20 N29-N31
      3 1219 1206 (N14W15 L=245)   603    204        0 (1118/1156-1304/1235) -36.0 W08-W19 N09-N17
      4  178 1192 (N25E63 L=168)  1787    323        4 ( 116/1152- 238/1219)   0.0 E71-E56 N22-N26
      5  450 1057 (N15E37 L=193)  2059    337        5 ( 398/ 873- 524/1218) -74.7 E41-E33 N03-N25
      6  545 1167 (N21E29 L=201)   334     56        0 ( 516/1140- 574/1198) -45.0 E31-E27 N19-N23
      7 1163  963 (N01W07 L=238)   441    126        0 (1096/ 924-1221/1006) -45.0 W03-W10 S02-N04
      8  233  897 (N07E59 L=172)   540    108        2 ( 207/ 872- 264/ 930)  31.0 E62-E55 N05-N08
      9  265  831 (N03E57 L=174)   700    114        3 ( 231/ 789- 293/ 870)  51.3 E61-E53 N00-N05
     10 1026  538 (S24E07 L=223)   372    130        1 (1007/ 473-1041/ 586)  69.4 E08-E05 S28-S20
     11 1425  258 (S58W25 L=256)   556    105        0 (1384/ 246-1467/ 273)  -5.2 W19-W30 S59-S55
  

    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