Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  650 1731 (N40E48 L=167)   513     45        1 ( 626/1707- 678/1755)  -7.1 E51-E43 N38-N42
      2 1503 1674 (N53W34 L=249)   385     54        0 (1475/1658-1535/1689)  18.4 W30-W39 N50-N54
      3 1358 1600 (N43W16 L=231)   279     62        0 (1328/1580-1386/1617) -11.3 W13-W18 N42-N44
      4 1648 1401 (N33W42 L=256)   363     88        0 (1605/1386-1698/1419)   0.0 W37-W47 N32-N34
      5  450  665 (S31E36 L=178)   378    140        1 ( 397/ 628- 505/ 694)  28.6 E42-E31 S34-S28
      6 1529  459 (S28W49 L=263)   585    114        0 (1475/ 441-1585/ 477)   5.2 W43-W54 S29-S25
      7 1167  362 (S41W24 L=239)   304     63        2 (1131/ 345-1199/ 381)   9.5 W21-W27 S41-S39
      8  745  213 (S66E13 L=201)   740     60        0 ( 718/ 187- 779/ 241) -45.0 E17-E07 S69-S62
  

    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