Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  781 1869 (N47E42 L=198)   746    140        0 ( 765/1838- 801/1902) -54.5 E46-E38 N43-N50
      2 1039 1536 (N25E08 L=232)   698    157        0 ( 996/1455-1089/1597)  79.7 E09-E06 N18-N29
      3 1694 1489 (N36W40 L=280)  1382    197        4 (1611/1452-1759/1525)   7.6 W32-W47 N34-N38
      4 1366 1436 (N24W14 L=254)   777    176        0 (1318/1341-1415/1517) -48.4 W10-W18 N19-N28
      5  371 1411 (N07E50 L=190)   460     83        0 ( 335/1363- 395/1453)  71.6 E52-E48 N04-N10
      6 1677 1217 (N17W38 L=278)   844    182        1 (1640/1136-1699/1301) -84.8 W35-W40 N12-N23
      7  407 1100 (S11E41 L=199)  1429    366        1 ( 331/ 924- 468/1261)  74.7 E45-E37 S22-N00
      8 1191 1208 (N07W07 L=247)   400    102        1 (1150/1164-1233/1245) -26.6 W03-W09 N04-N08
      9  167  797 (S31E61 L=179)  1019    323        0 ( 128/ 728- 210/ 859)  35.0 E66-E56 S34-S27
     10 1344  797 (S14W23 L=263)  2875    747        2 (1016/ 729-1651/ 850)  18.0 W03-W43 S19-S06
     11  617  626 (S37E21 L=219)  1934    553        6 ( 424/ 537- 867/ 696)  26.6 E38-E04 S45-S28
     12 1708  490 (S23W62 L=302)  1899    362        2 (1580/ 457-1795/ 535)  22.8 W50-W69 S26-S17
  

    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