Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  559 1700 (N38E44 L=353)  1946    184        4 ( 465/1656- 644/1734)  -9.9 E57-E34 N34-N40
      2 1106 1669 (N38W02 L= 39)   358     92        0 (1059/1654-1153/1683)   9.5 E01-W05 N37-N39
      3  854 1395 (N17E13 L= 24)   677    175        0 ( 801/1320- 895/1472)  68.2 E15-E10 N12-N22
      4 1795 1411 (N26W58 L= 96)  1304     95        2 (1748/1373-1834/1447) -18.4 W53-W64 N24-N28
      5  298 1375 (N14E55 L=343)   968    224        1 ( 231/1317- 377/1435) -21.8 E65-E45 N10-N18
      6 1122 1252 (N10W05 L= 42)   212     64        1 (1088/1236-1158/1273) -14.0 W02-W06 N08-N10
      7  828 1164 (N03E13 L= 25)   283    117        0 ( 786/1122- 867/1200) -38.7 E15-E10 N00-N04
      8  347 1104 (S02E45 L=352)   418     89        0 ( 325/1066- 375/1152)  53.1 E47-E43 S04-N00
      9  143 1025 (S07E67 L=330)   674     55        0 ( 128/ 989- 157/1066)  76.0 E69-E65 S08-S04
     10 1821  969 (S01W56 L= 94)  1236    128        2 (1792/ 897-1846/1041) -74.1 W53-W60 S05-N03
     11 1085  994 (S06W04 L= 41)   202     12        0 (1064/ 975-1107/1013)  26.6 W02-W05 S07-S04
     12  757  864 (S15E16 L= 22)   187    160        0 ( 666/ 849- 815/ 877)   6.3 E21-E12 S16-S14
     13 1642  841 (S11W42 L= 80)   245     34        0 (1620/ 822-1659/ 863)  33.7 W40-W43 S12-S09
     14 1719  787 (S13W50 L= 87)   470    110        1 (1673/ 757-1763/ 809)  -6.3 W45-W54 S15-S11
     15 1648  709 (S19W46 L= 83)   570    100        2 (1600/ 668-1688/ 736)  45.0 W42-W48 S21-S16
     16 1686  662 (S21W51 L= 88)   857    203        1 (1610/ 583-1753/ 725)  51.3 W45-W55 S26-S16
     17 1455  690 (S22W31 L= 68)   287     99        0 (1411/ 661-1498/ 721) -15.9 W27-W34 S23-S19
     18 1354  402 (S42W33 L= 70)   260     33        0 (1337/ 381-1376/ 429) -31.0 W30-W35 S43-S39
     19 1466  405 (S40W43 L= 80)   547    113        1 (1419/ 373-1532/ 424)  10.3 W37-W48 S42-S38
  

    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