Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  653 1750 (N62E06 L= 52)   430     46        1 ( 626/1726- 678/1770) -23.2 E09-E02 N59-N64
      2  456 1665 (N62E35 L= 23)  1678    200        3 ( 383/1605- 550/1710)   0.0 E49-E20 N57-N65
      3 1533 1263 (N05W36 L= 94)   355    144        1 (1479/1244-1586/1288) -23.2 W33-W40 N02-N07
      4 1790 1235 (S04W55 L=114)  1193    141        3 (1732/1189-1851/1274)   4.1 W48-W62 S06-S02
      5 1851 1144 (S11W58 L=117)   604     86        1 (1826/1112-1888/1179) -39.8 W56-W62 S13-S08
      6  324 1105 (N25E44 L= 14)   612     72        1 ( 295/1063- 356/1146)  68.2 E47-E41 N22-N27
      7  848 1119 (N14E08 L= 51)   243     79        1 ( 806/1101- 874/1142)   0.0 E10-E05 N12-N15
      8 1106 1114 (N07W07 L= 65)   208     59        1 (1073/1096-1137/1137) -14.0 W04-W08 N06-N08
      9 1103 1067 (N05W05 L= 64)   212     24        1 (1080/1046-1124/1087)   0.0 W04-W06 N03-N06
     10  887  813 (S04E13 L= 46)  1015    324        0 ( 758/ 723- 993/ 925)  15.5 E22-E04 S05-N00
     11 1684  571 (S41W36 L= 94)   542     74        0 (1656/ 542-1717/ 608)  26.6 W32-W39 S43-S39
     12  550  522 (S13E43 L= 16)   321     66        1 ( 537/ 487- 564/ 563)  63.4 E44-E40 S14-S10
     13  488  482 (S14E50 L=  9)   874    239        2 ( 454/ 369- 542/ 552)  72.9 E53-E46 S22-S09
     14  337  443 (S14E68 L=351)   756    101        0 ( 322/ 393- 355/ 487)  41.2 E72-E62 S17-S10
     15 1404  391 (S45W08 L= 66)   296     56        0 (1381/ 374-1432/ 411)  14.0 W05-W10 S46-S43
     16 1145  297 (S44E15 L= 43)   612    120        1 (1082/ 277-1214/ 322) -16.7 E20-E09 S45-S40
     17  920  224 (S42E36 L= 22)   507     86        0 ( 884/ 198- 963/ 249)  11.3 E41-E31 S43-S40
  

    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