Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1686 1628 (N19W69 L=173)  1898    138        4 (1660/1562-1720/1679)  15.5 W62-W83 N15-N21
      2 1300 1663 (N33W39 L=143)   373     88        0 (1254/1648-1346/1675) -26.6 W35-W41 N31-N34
      3 1310 1562 (N27W34 L=138)   693    199        1 (1243/1485-1381/1596)   4.4 W26-W39 N23-N28
      4 1473 1486 (N18W41 L=145)   276     76        1 (1439/1459-1501/1511)   9.5 W37-W43 N16-N18
      5 1391 1415 (N16W32 L=136)   213     67        1 (1366/1376-1414/1443)  78.7 W30-W33 N13-N18
      6 1121 1346 (N19W14 L=118)   199     89        0 (1087/1311-1146/1384) -80.5 W12-W14 N16-N22
      7  844 1282 (N23E04 L=100)   837    210        1 ( 791/1195- 924/1357) -71.6 E06-W00 N15-N27
      8 1735 1152 (S07W47 L=152)  2128    407        1 (1679/ 966-1775/1308)  63.4 W40-W52 S16-N02
      9  583 1218 (N25E22 L= 82)  2067    553        5 ( 361/1079- 789/1293) -34.2 E39-E11 N12-N32
     10  481 1079 (N19E32 L= 72)   249    117        2 ( 425/1062- 538/1095) -26.6 E35-E27 N17-N21
     11 1543  843 (S19W26 L=130)  1422    331        2 (1412/ 740-1667/ 919)   8.1 W15-W36 S22-S16
     12 1781  825 (S28W46 L=150)  2871    462        6 (1564/ 685-1918/ 918)   1.4 W25-W65 S30-S24
     13  366  771 (N04E47 L= 57)   655    158        1 ( 338/ 699- 401/ 832)  45.0 E51-E43 N00-N07
     14  745  719 (S07E23 L= 81)   322    122        0 ( 717/ 672- 763/ 770)  59.0 E26-E21 S08-S03
     15  366  585 (S06E55 L= 49)  6447      0       17 ( 242/ 440- 459/ 726) -71.6 E66-E45 S16-N02
     16  533  524 (S13E44 L= 60)   267     61        0 ( 504/ 503- 564/ 548)   9.5 E46-E40 S14-S12
     17  644  348 (S27E46 L= 58)   831    164        3 ( 577/ 308- 704/ 405)   7.6 E53-E38 S28-S24
     18 1407  235 (S62W05 L=109)   573     57        0 (1383/ 212-1442/ 260) -52.1 W02-W09 S67-S57
  

    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