Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  399 1611 (N46E60 L= 52)   673     55        1 ( 383/1575- 418/1643) -36.9 E66-E55 N42-N48
      2  920 1547 (N39E05 L=108)   481    133        0 ( 855/1528- 976/1563) -12.5 E09-W00 N37-N40
      3  542 1352 (N27E33 L= 80)   470    172        1 ( 511/1286- 565/1412) -83.7 E35-E31 N22-N31
      4  396 1309 (N25E45 L= 68)   732    141        2 ( 380/1244- 416/1376) -60.9 E48-E41 N20-N29
      5  301 1252 (N21E53 L= 60)   441     59        1 ( 285/1220- 326/1283)  90.0 E54-E50 N19-N23
      6 1733 1234 (N12W52 L=165)   369     63        1 (1710/1200-1750/1258) -76.0 W50-W53 N09-N13
      7 1746  967 (S05W50 L=162)   773    190        2 (1675/ 903-1796/1017) -28.6 W43-W54 S08-S01
      8 1709  917 (S07W46 L=159)   339    121        3 (1685/ 860-1721/ 961) -90.0 W44-W47 S10-S04
      9 1766  808 (S15W52 L=165)  1877    298        7 (1673/ 727-1837/ 928)  38.2 W44-W60 S19-S07
     10 1650  867 (S10W41 L=154)   281    116        0 (1593/ 849-1695/ 896)   6.3 W36-W45 S10-S08
     11 1253  800 (S10W13 L=126)   204     46        1 (1229/ 776-1277/ 825)  33.7 W11-W14 S11-S08
     12 1639  781 (S15W41 L=153)   542    146        2 (1583/ 733-1689/ 814)   0.0 W36-W45 S18-S12
     13 1690  652 (S24W47 L=160)  1276    239        6 (1650/ 567-1737/ 735)  68.2 W43-W51 S30-S18
     14  812  539 (S25E17 L= 95)   774    136        0 ( 784/ 473- 845/ 605)  60.9 E20-E14 S29-S20
     15  330  505 (S27E62 L= 50)   571     48        0 ( 301/ 493- 361/ 520)  -7.1 E67-E57 S27-S25
     16  694  426 (S33E29 L= 84)   932    160        2 ( 626/ 396- 775/ 456)   8.7 E35-E22 S35-S30
  

    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