Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1732 1529 (N12W67 L=131)  1555    150        4 (1682/1467-1782/1582) -69.4 W62-W70 N07-N15
      2 1179 1464 (N25W22 L= 86)  2087    629        5 ( 928/1413-1487/1514) -32.0 W08-W40 N14-N35
      3  284 1474 (N52E56 L=  7)  2089    181        5 ( 196/1439- 362/1509)   0.0 E82-E43 N48-N54
      4  865 1394 (N30W01 L= 65)   253     42        0 ( 839/1374- 893/1412) -71.6 E01-W02 N28-N31
      5  489 1181 (N27E30 L= 34)   221     62        0 ( 452/1171- 525/1190) -21.8 E32-E26 N26-N28
      6  413  999 (N18E39 L= 25)  1506    386        2 ( 381/ 856- 479/1182)  81.5 E42-E34 N09-N29
      7 1581 1098 (S06W35 L= 98)   242     38        0 (1560/1071-1600/1117) -90.0 W33-W36 S07-S04
      8 1835  947 (S23W53 L=117)   492     64        0 (1809/ 913-1858/ 983)  18.4 W50-W56 S24-S20
      9 1260  673 (S21W04 L= 68)   595    249        1 (1143/ 658-1374/ 704) -21.0 E02-W11 S23-S17
     10  448  490 (S12E54 L= 10)   541     76        0 ( 421/ 457- 487/ 517)   6.3 E57-E48 S13-S10
     11  659  355 (S25E45 L= 18)   666    141        0 ( 608/ 312- 684/ 411)  22.6 E52-E40 S27-S21
     12  930  359 (S32E27 L= 37)   435     91        2 ( 905/ 320- 958/ 382)  26.6 E30-E24 S34-S29
  

    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