Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  837 1865 (N55E48 L=204)   678     60        1 ( 801/1846- 873/1879)   5.2 E54-E41 N53-N56
      2 1178 1729 (N50E03 L=249)   360     57        0 (1160/1700-1208/1750)  59.0 E04-W00 N47-N52
      3  487 1678 (N32E60 L=191)   883    107        2 ( 448/1642- 518/1721) -11.9 E72-E53 N30-N34
      4 1339 1625 (N45W13 L=265)   425     44        1 (1311/1603-1369/1646)   0.0 W10-W16 N42-N46
      5  792 1553 (N29E26 L=226)   416    166        0 ( 715/1524- 858/1573)  29.1 E30-E21 N25-N31
      6 1410 1400 (N29W20 L=272)   261    103        0 (1370/1362-1453/1426) -18.4 W17-W23 N27-N30
      7 1809 1358 (N34W61 L=313)   479     50        0 (1779/1338-1831/1380)  15.3 W55-W66 N31-N35
      8 1614 1241 (N22W37 L=288)   675    226        0 (1524/1175-1677/1317) -24.8 W29-W42 N19-N25
      9  192  829 (S25E64 L=188)   738     84        4 ( 165/ 794- 223/ 859)  15.3 E70-E59 S27-S23
     10  738  469 (S41E11 L=240)   257     72        0 ( 713/ 454- 772/ 496) -11.3 E13-E08 S42-S39
     11 1505  427 (S30W49 L=301)   723    114        0 (1451/ 401-1559/ 454)  26.6 W45-W53 S32-S27
     12 1061  196 (S58W30 L=282)   697    101        0 (1026/ 156-1088/ 241) -25.2 W21-W38 S63-S54
  

    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