Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1032 1906 (N54E44 L= 85)   819     63        1 (1004/1876-1055/1933) -45.0 E51-E39 N49-N58
      2  495 1673 (N19E59 L= 70)  1121    172        2 ( 458/1607- 544/1733) -68.2 E63-E55 N14-N24
      3  544 1596 (N15E49 L= 81)   772    228        0 ( 500/1505- 573/1649) -28.3 E55-E42 N09-N18
      4  700 1567 (N17E36 L= 93)   254     81        2 ( 668/1549- 729/1590)  63.4 E38-E33 N15-N19
      5  799 1523 (N17E28 L=101)   536    121        3 ( 752/1472- 851/1556)   0.0 E31-E23 N15-N19
      6  746 1490 (N14E30 L=100)   276     35        2 ( 727/1465- 769/1517) -56.3 E31-E28 N12-N15
      7  956  675 (S27W06 L=136)   250    106        0 ( 911/ 657-1009/ 691)  38.7 W04-W09 S28-S24
      8 1708  651 (S05W56 L=186)   394     32        1 (1681/ 631-1729/ 667)  26.6 W53-W58 S05-S03
      9 1549  505 (S18W51 L=181)   312     49        1 (1520/ 490-1570/ 525)  56.3 W49-W53 S18-S15
     10 1107  450 (S35W25 L=155)   233     40        0 (1082/ 437-1138/ 465)  33.7 W23-W27 S36-S34
     11  669  428 (S52E05 L=124)   596     76        0 ( 635/ 402- 711/ 460)   7.1 E09-E01 S54-S50
     12 1496  422 (S24W53 L=183)   410     71        0 (1459/ 392-1523/ 447)  78.7 W51-W54 S26-S21
  

    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