Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1123 1886 (N57E34 L=  1)   966     67        2 (1095/1857-1150/1919) -18.4 E42-E27 N53-N61
      2 1490 1706 (N56W14 L= 49)   593     79        1 (1467/1660-1514/1736) -66.8 W11-W17 N51-N60
      3 1698 1573 (N53W42 L= 76)   730     85        1 (1674/1535-1719/1608) -90.0 W37-W46 N49-N56
      4 1643 1249 (N26W33 L= 68)   426    130        0 (1587/1210-1668/1298)  58.0 W29-W35 N21-N29
      5 1733 1079 (N19W44 L= 78)   825    260        5 (1645/1012-1809/1159) -14.9 W35-W52 N15-N21
      6 1895  964 (N18W65 L=100)  1099    120        3 (1869/ 914-1923/1015)  60.3 W61-W70 N14-N21
      7  255  928 (S29E53 L=342)   733    198        2 ( 226/ 817- 282/ 999) -70.0 E56-E49 S35-S24
      8 1726  952 (N11W46 L= 81)   414     83        4 (1681/ 931-1764/ 975)  15.9 W42-W49 N10-N13
      9 1841  884 (N11W60 L= 95)  1102    134        3 (1795/ 829-1870/ 958) -50.2 W56-W64 N07-N16
     10  365  853 (S31E40 L=354)   318    114        1 ( 317/ 827- 422/ 882)  26.6 E45-E35 S33-S28
     11 1211  402 (S35W34 L= 69)   610    100        1 (1169/ 385-1266/ 430)  15.9 W30-W38 S36-S32
     12  510  393 (S61E26 L=  9)   418     55        2 ( 481/ 365- 533/ 415) -23.2 E30-E22 S63-S59
  

    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