Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1154 1650 (N37W29 L=102)   277     68        1 (1131/1618-1179/1671)  11.3 W26-W31 N35-N38
      2 1684 1607 (N17W65 L=137)   734    208        0 (1633/1588-1744/1627) -22.6 W60-W72 N14-N19
      3 1438 1416 (N15W36 L=108)   701    173        1 (1345/1366-1518/1453) -48.4 W31-W39 N10-N19
      4  624 1111 (N20E21 L= 51)   436    110        2 ( 569/1086- 670/1133) -39.8 E24-E18 N17-N22
      5 1907 1081 (S18W63 L=135)  1412    195        2 (1864/1037-1942/1124) -48.8 W58-W67 S21-S13
      6 1532  826 (S20W25 L= 97)   504    158       -1 (1505/ 792-1569/ 865)  36.9 W22-W27 S22-S17
      7 1638  635 (S36W30 L=102)  2010    259        1 (1545/ 523-1698/ 767)  40.6 W20-W35 S40-S28
      8 1104  539 (S25E09 L= 63)  1072    271        2 (1063/ 414-1146/ 646)  65.0 E12-E03 S33-S18
      9  422  600 (S05E50 L= 23)   318     58        0 ( 408/ 569- 438/ 628)  76.0 E50-E47 S06-S02
     10  779  415 (S24E33 L= 39)  1007    239        0 ( 718/ 344- 810/ 508)  36.9 E40-E27 S27-S18
     11  503  422 (S17E52 L= 20)  1944    221        2 ( 457/ 342- 552/ 506)  65.6 E55-E46 S22-S11
  

    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