Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1818 1104 (N07W55 L= 42)   371    121        1 (1785/1070-1854/1161) -29.7 W51-W58 N05-N09
      2 1661 1067 (N03W41 L= 28)   321    109        0 (1617/1042-1698/1100)  39.8 W37-W43 N00-N05
      3 1091 1055 (S02W04 L=351)   234     60        2 (1069/1017-1122/1089)  90.0 W02-W05 S04-N00
      4  635  925 (S13E24 L=324)   212     78        0 ( 604/ 882- 660/ 970) -59.0 E25-E21 S14-S09
      5 1440  884 (S10W27 L= 14)   547    251        1 (1324/ 845-1548/ 917)  23.2 W19-W33 S13-S07
      6  128  813 (S19E69 L=278)  2511    226        1 ( 103/ 721- 148/ 889)  73.3 E72-E66 S24-S14
      7  297  734 (S25E51 L=296)  1527    440        5 ( 251/ 565- 348/ 881)  57.0 E59-E45 S35-S15
      8  862  846 (S16E09 L=338)   337    107        0 ( 805/ 822- 921/ 868)  23.2 E12-E05 S17-S14
      9 1583  746 (S17W39 L= 26)   333    112        0 (1553/ 707-1626/ 774)  45.0 W37-W42 S19-S14
     10  647  659 (S29E24 L=324)   268     89        2 ( 615/ 626- 685/ 699)  45.0 E25-E20 S31-S26
     11 1273  637 (S26W19 L=  6)   316     97        1 (1223/ 599-1321/ 678)  50.2 W16-W21 S29-S23
     12 1455  534 (S31W35 L= 22)  2941    622        7 (1168/ 423-1714/ 664)  28.4 W15-W52 S41-S20
  

    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