Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  547 1699 (N24E55 L=221)  2575    338        3 ( 434/1600- 668/1781)  66.4 E59-E49 N15-N32
      2 1130 1510 (N27E07 L=270)  1822    418        5 (1019/1327-1190/1645) -41.6 E18-W01 N16-N34
      3  698 1452 (N12E31 L=246)   321     85        1 ( 675/1413- 729/1488)  90.0 E32-E29 N09-N14
      4 1867 1143 (N26W58 L=335)  3099    396        3 (1795/1013-1946/1302) -47.0 W48-W73 N19-N34
      5  346 1146 (S13E46 L=230)   379    100        0 ( 332/1093- 363/1199) -50.2 E48-E43 S15-S09
      6  170 1024 (S23E64 L=213)  2214    313        3 ( 128/ 938- 240/1122)  47.5 E70-E55 S28-S16
      7  479  977 (S20E32 L=244)   446    127        0 ( 457/ 932- 518/1033) -51.3 E34-E28 S21-S16
      8 1653  883 (N05W42 L=319)   904    204        2 (1584/ 814-1698/ 978)  68.2 W38-W44 N00-N10
      9  503  507 (S50E26 L=251)  1690    267        2 ( 400/ 434- 611/ 592)  -9.5 E37-E13 S53-S45
     10 1726  485 (S14W69 L=346)  1208    144        1 (1659/ 461-1783/ 527)  35.0 W62-W76 S17-S10
  

    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