Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1419 1497 (N22W39 L= 15)   548    113        1 (1391/1445-1453/1539)  20.6 W35-W43 N19-N23
      2  891 1485 (N37W05 L=341)   367    101        2 ( 847/1459- 937/1520)   7.1 W01-W09 N35-N37
      3 1176 1476 (N28W23 L=359)   260     93        0 (1133/1460-1223/1491)  -8.1 W19-W26 N26-N28
      4  701 1447 (N40E09 L=327)   393     80        0 ( 671/1418- 750/1467) -31.0 E11-E05 N38-N41
      5 1674 1306 (N04W48 L= 24)   852    180        3 (1602/1260-1749/1372)   0.0 W41-W56 N02-N05
      6  165 1265 (N40E69 L=267)   950     73        1 ( 124/1238- 200/1287)  -8.1 E82-E61 N38-N41
      7  618 1183 (N26E21 L=315)  1431    332        3 ( 461/1146- 761/1212) -24.2 E32-E12 N20-N30
      8  230  921 (N18E58 L=278)  3008    526        4 ( 125/ 781- 457/1016)   2.9 E76-E36 N10-N24
      9 1645  919 (S16W35 L= 10)   960    204        3 (1604/ 850-1687/1000)  36.9 W30-W39 S19-S11
     10 1408  896 (S10W18 L=354)   845    210        0 (1338/ 814-1451/ 986)  31.0 W12-W22 S13-S06
     11  186  718 (N07E70 L=265)   950    100        2 ( 163/ 661- 209/ 760)  80.5 E73-E65 N03-N09
     12 1158  665 (S17E02 L=334)   305    103        1 (1127/ 622-1189/ 713)  38.7 E04-W01 S18-S14
     13  935  651 (S12E15 L=321)   320     70        2 ( 905/ 621- 967/ 677)  26.6 E16-E12 S12-S10
     14  971  612 (S15E14 L=322)   246     16        0 ( 948/ 593- 995/ 631) -18.4 E15-E12 S15-S13
     15  387  411 (S15E65 L=270)  1416    158        5 ( 339/ 353- 425/ 454)  81.9 E71-E59 S19-S11
     16  604  365 (S22E48 L=288)   894    124        1 ( 580/ 299- 637/ 427)  36.9 E54-E42 S26-S17
     17 1183  319 (S41E12 L=323)   557    104        1 (1162/ 268-1203/ 378)  56.3 E15-E08 S45-S36
     18 1404  263 (S56W03 L=338)  1483    172        2 (1320/ 236-1480/ 295) -23.6 E05-W11 S61-S50
  

    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