Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1301 1759 (N53W38 L= 54)   782    104        1 (1260/1729-1348/1777)   0.0 W31-W44 N50-N54
      2 1544 1633 (N40W53 L= 70)   486     68        1 (1518/1608-1579/1661) -68.2 W49-W56 N37-N42
      3  546 1489 (N37E36 L=341)   523    137        0 ( 511/1430- 586/1534)  63.4 E37-E33 N32-N40
      4  688 1433 (N32E22 L=354)   947    193        1 ( 628/1379- 766/1502) -36.9 E27-E15 N28-N37
      5 1520 1353 (N21W37 L= 54)   365    146        1 (1501/1300-1548/1415)  74.1 W34-W40 N18-N25
      6  824 1194 (N16E12 L=  5)   800    373        1 ( 691/1092- 965/1263) -25.2 E20-E03 N09-N20
      7 1236 1144 (N11W14 L= 31)   323     85        1 (1200/1112-1272/1169) -53.1 W11-W16 N08-N12
      8 1226 1078 (N07W13 L= 29)   405     68        3 (1192/1048-1259/1107)  36.9 W10-W14 N05-N08
      9 1798  820 (S15W56 L= 73)   616    150        1 (1760/ 770-1827/ 858)  39.8 W52-W60 S17-S12
     10 1580  817 (S12W36 L= 52)   516    113        1 (1535/ 779-1630/ 850)  33.7 W32-W39 S14-S10
     11  319  512 (S27E64 L=312)  1219    153        2 ( 276/ 443- 356/ 567) -90.0 E69-E61 S31-S22
  

    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