Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1201 1838 (N48W50 L=135)   438     67        0 (1180/1810-1235/1857)   5.2 W45-W56 N46-N49
      2 1372 1816 (N40W61 L=146)   518     49        1 (1356/1788-1394/1843)  16.7 W56-W66 N38-N42
      3 1509 1737 (N31W64 L=149)  1707    241        2 (1426/1656-1578/1812)  65.2 W59-W68 N24-N37
      4  409 1516 (N51E37 L= 48)   352     62        0 ( 392/1491- 430/1547) -45.0 E40-E34 N48-N53
      5  698 1428 (N36E10 L= 75)   313    119        0 ( 645/1399- 750/1454) -45.0 E13-E07 N32-N38
      6  357 1423 (N45E42 L= 43)   788    125        4 ( 308/1391- 414/1446)  -8.7 E48-E35 N42-N47
      7  949 1424 (N29W06 L= 91)   209    113        0 ( 896/1397- 986/1438)   0.0 W02-W09 N27-N29
      8  614 1293 (N30E19 L= 66)   208     83        0 ( 596/1260- 647/1316)  26.6 E20-E16 N27-N30
      9  287 1182 (N31E48 L= 37)   418     51        3 ( 269/1154- 315/1208)  33.7 E50-E45 N28-N32
     10  322 1088 (N24E45 L= 40)   722    163        4 ( 268/1053- 393/1149) -35.0 E50-E39 N21-N28
     11  410  952 (N14E40 L= 45)   616    125        1 ( 380/ 887- 453/ 999)  81.9 E42-E35 N10-N17
     12 1322  942 (S08W15 L=100)   301    110        1 (1301/ 886-1351/ 984)  45.0 W12-W17 S10-S05
     13 1291  678 (S22W07 L= 92)   306    108        0 (1280/ 617-1304/ 736)  76.0 W05-W08 S26-S18
     14 1849  694 (S40W58 L=143)   667     41        0 (1819/ 679-1881/ 708) -13.0 W52-W66 S41-S37
     15 1083  543 (S25E09 L= 76)   217     87        0 (1034/ 530-1121/ 554)  -9.5 E12-E06 S26-S23
     16  533  311 (S27E60 L= 25)  1345    162        2 ( 484/ 256- 584/ 357)  13.0 E77-E51 S30-S24
     17 1313  282 (S52E03 L= 82)   633     74        2 (1270/ 261-1352/ 301) -15.9 E07-W00 S55-S49
     18 1209  189 (S58E19 L= 66)   698     50        2 (1185/ 159-1231/ 217)  33.7 E24-E15 S61-S54
  

    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