Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  728 1744 (N41E43 L=320)   462     85        1 ( 699/1705- 751/1785) -54.5 E46-E40 N37-N44
      2 1215 1447 (N29W05 L=  8)   330    157        0 (1182/1385-1263/1509) -45.0 W01-W09 N25-N32
      3  499 1425 (N15E43 L=320)   455     97        1 ( 472/1376- 524/1468) -35.5 E46-E39 N12-N17
      4 1525 1255 (N22W29 L= 32)   613    132        3 (1478/1216-1559/1318)  60.3 W25-W31 N19-N26
      5  725  917 (S12E17 L=347)   227    131        0 ( 713/ 861- 747/ 979) -63.4 E18-E14 S14-S08
      6  633  874 (S16E22 L=341)   726    241        2 ( 564/ 807- 688/ 935)  90.0 E27-E18 S20-S12
      7  219  775 (S30E61 L=302)  1741    254        4 ( 158/ 697- 288/ 857)  21.2 E83-E50 S35-S23
      8 1127  614 (S24W14 L= 17)  2809    860        9 ( 810/ 492-1470/ 833)  41.6 E05-W31 S36-S04
      9  576  514 (S41E25 L=338)  1041    203        4 ( 501/ 474- 662/ 555)  23.6 E33-E17 S45-S37
     10 1719  519 (S19W66 L= 69)  1069    143        1 (1658/ 489-1789/ 541)  17.4 W59-W77 S21-S16
     11 1457  410 (S31W47 L= 50)  1013    170        1 (1387/ 379-1532/ 451)  26.6 W40-W52 S34-S27
     12 1247  304 (S43W37 L= 40)   657     67        2 (1225/ 269-1281/ 344) -29.1 W32-W41 S46-S40
  

    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