Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1215 1847 (N50W48 L=116)   586     60        1 (1189/1818-1246/1867)  11.3 W42-W53 N47-N51
      2 1460 1483 (N19W39 L=107)  1553    385        1 (1292/1431-1600/1583)   0.0 W26-W55 N16-N22
      3  944 1517 (N33W06 L= 75)  1644    384        2 ( 798/1442-1111/1569) -32.9 E03-W14 N24-N38
      4  347 1386 (N38E44 L= 25)  1719    180        3 ( 293/1299- 395/1453) -47.5 E52-E37 N31-N43
      5  359 1078 (N18E42 L= 27)  1390    251        3 ( 310/ 982- 411/1179) -66.8 E46-E38 N11-N25
      6  113 1069 (N22E73 L=356)  2856    219        5 (  86/ 966- 150/1152) -59.7 E81-E66 N15-N27
      7 1266  838 (S13W10 L= 79)   415    138        0 (1240/ 770-1288/ 906)  54.5 W07-W12 S16-S09
      8 1734  787 (S27W44 L=113)   477     71        1 (1705/ 758-1761/ 820) -45.0 W41-W46 S29-S24
      9  974  684 (S16E10 L= 59)   522    149        0 ( 924/ 637-1044/ 731) -63.4 E11-E06 S20-S12
     10 1636  636 (S35W35 L=104)   287     54        0 (1608/ 622-1664/ 652) -21.8 W32-W37 S36-S32
     11  440  451 (S20E56 L= 13)  1553    208        6 ( 411/ 356- 484/ 532)  45.0 E65-E49 S26-S14
     12  622  406 (S26E43 L= 26)   427     91        1 ( 601/ 367- 641/ 457)  45.0 E47-E40 S28-S22
  

    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