Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  366 1537 (N20E54 L=220)   395     81        0 ( 330/1516- 400/1553)  18.4 E56-E50 N18-N21
      2 1128 1272 (N10W03 L=277)   832    209        5 (1072/1205-1187/1349) -45.0 W00-W07 N06-N14
      3 1826 1300 (N23W54 L=328)   637    163        0 (1785/1270-1858/1341)  40.6 W49-W58 N20-N26
      4 1730 1022 (N04W46 L=320)   472    107        1 (1701/ 982-1767/1066) -26.6 W43-W49 N02-N06
      5 1786 1044 (N06W51 L=325)   323     96        0 (1754/1015-1812/1063)   0.0 W48-W53 N04-N07
      6  125  986 (S14E68 L=206)   513    147        1 ( 102/ 960- 153/1003)   8.1 E71-E63 S15-S13
      7  409  943 (S16E39 L=234)   435    163        0 ( 367/ 891- 443/1000) -56.3 E42-E36 S19-S12
      8 1870  912 (N00W62 L=336)   405     97        1 (1843/ 882-1895/ 950)  45.0 W59-W65 S01-N02
      9  494  784 (S25E33 L=241)   578    216        4 ( 420/ 737- 544/ 847)  37.9 E38-E28 S28-S21
     10 1735  620 (S19W57 L=330)   839    168        3 (1687/ 575-1794/ 647)  36.9 W53-W62 S22-S16
     11 1344  441 (S37W33 L=306)  1789    402        5 (1181/ 413-1550/ 491)  21.0 W20-W46 S40-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