In [1]:
Copied!
import mrcfile
import numpy as np
import matplotlib.pyplot as plt
from src.pickpoints import *
from src.bezierfit import *
import mrcfile
import numpy as np
import matplotlib.pyplot as plt
from src.pickpoints import *
from src.bezierfit import *
Set the paramters¶
In [2]:
Copied!
filename_2Daverage = './db/templates_selected_bezierfit.mrc'
section = 16
initial_data_points_num = 500
filename_2Daverage = './db/templates_selected_bezierfit.mrc'
section = 16
initial_data_points_num = 500
Generate the data points¶
In [3]:
Copied!
with mrcfile.open(filename_2Daverage) as f:
image = f.data[section]
max_filtered_image = max_filter(image)
# 生成数据点
data_points = generate_data_points(image, initial_data_points_num)
with mrcfile.open(filename_2Daverage) as f:
image = f.data[section]
max_filtered_image = max_filter(image)
# 生成数据点
data_points = generate_data_points(image, initial_data_points_num)
Show the original image, max-filteed image, and the data points
In [4]:
Copied!
plt.figure(figsize=(10, 5))
plt.subplot(1, 3, 1)
plt.imshow(image, cmap='gray', origin='lower')
plt.title('Original Image')
plt.axis('off')
plt.subplot(1, 3, 2)
plt.imshow(max_filtered_image, cmap='gray', origin='lower')
plt.title('Max Filtered Image')
plt.axis('off')
plt.subplot(1, 3, 3)
plt.imshow(image, cmap='gray', origin='lower')
plt.title('Data Points')
plt.axis('off')
plt.plot(data_points[:, 0], data_points[:, 1], 'b.')
plt.tight_layout()
plt.show()
plt.figure(figsize=(10, 5))
plt.subplot(1, 3, 1)
plt.imshow(image, cmap='gray', origin='lower')
plt.title('Original Image')
plt.axis('off')
plt.subplot(1, 3, 2)
plt.imshow(max_filtered_image, cmap='gray', origin='lower')
plt.title('Max Filtered Image')
plt.axis('off')
plt.subplot(1, 3, 3)
plt.imshow(image, cmap='gray', origin='lower')
plt.title('Data Points')
plt.axis('off')
plt.plot(data_points[:, 0], data_points[:, 1], 'b.')
plt.tight_layout()
plt.show()
Coarse fit: generate intial control points¶
In [5]:
Copied!
bezier_degree = 3
coarse_iteration = 700
bezier_degree = 3
coarse_iteration = 700
In [6]:
Copied!
coarsefit = Coarsefit(image=image, num_points=initial_data_points_num, degree=bezier_degree, iteration=coarse_iteration)
initial_control_points = coarsefit()
print('Initial Control Points:', initial_control_points)
coarsefit = Coarsefit(image=image, num_points=initial_data_points_num, degree=bezier_degree, iteration=coarse_iteration)
initial_control_points = coarsefit()
print('Initial Control Points:', initial_control_points)
gen nevals avg std min max 0 100 1.93298e+07 7.45036e+06 3.93846e+06 3.83652e+07 1 77 1.10273e+07 4.44083e+06 3.38866e+06 2.60016e+07 2 67 6.60018e+06 2.38678e+06 2.6819e+06 1.65696e+07 3 73 4.40417e+06 1.62709e+06 1.53115e+06 1.0306e+07 4 67 2.97223e+06 940933 1.31142e+06 7.97569e+06 5 71 2.13214e+06 530194 1.17533e+06 3.58231e+06 6 80 1.63262e+06 480155 855839 3.43829e+06 7 77 1.22557e+06 270201 557494 2.1998e+06 8 71 955983 213393 447932 1.68143e+06 9 63 770947 184040 411039 1.37217e+06 10 66 608170 164867 363731 1.63209e+06 11 76 471268 84517.5 335800 875004 12 75 402217 65320.7 282429 764520 13 72 364355 64332.5 249205 593900 14 78 307629 40947 249374 473903 15 77 274175 39973.4 241838 568967 16 69 257194 21318.6 225375 377353 17 77 249405 18557.8 217120 369865 18 80 237679 13909.3 207754 290000 19 62 227083 12381.8 203763 275653 20 68 216277 8398.24 203234 256801 21 69 208990 5482.95 194853 228478 22 72 203920 3707.32 191126 216292 23 78 200288 5395.62 189291 217650 24 66 194851 4505.88 184760 210028 25 72 190041 3696.42 175608 201486 26 69 185691 4200.33 174370 199368 27 74 180568 4039.82 170599 187977 28 74 176310 3329.89 169470 186728 29 71 173450 2645.35 167552 183546 30 73 170948 2351.78 166573 176966 31 77 168841 1920.22 163660 173963 32 65 167178 1643.46 162596 174764 33 70 165686 1702.06 159819 171237 34 76 163864 2085.02 157636 169755 35 83 161373 1939.5 153670 166597 36 68 159108 1755.73 153288 162774 37 71 156907 2050.2 150813 162282 38 74 154345 2146.01 149616 159125 39 78 151961 1356.37 147496 155552 40 76 150348 1543.68 144023 154020 41 66 148530 2311.19 140963 154154 42 77 145529 2794.44 138214 151384 43 73 142033 2174.68 137179 151633 44 74 139879 1739.68 136464 146888 45 81 137951 1158.93 135144 141423 46 72 136780 1065.08 134223 140680 47 72 135583 1058.09 129142 139174 48 82 134635 1438.91 128823 138762 49 60 132632 2278.95 127047 136828 50 77 129798 1858.74 125466 135096 51 68 127739 1603 123044 131698 52 62 126072 1387.45 122668 129465 53 71 124506 1497.42 119741 128509 54 64 122802 1916.03 117433 126666 55 61 120207 2182.29 113964 126650 56 76 117660 2370.47 112090 123520 57 73 114754 1940.27 111012 122515 58 69 112642 1190.88 109217 116409 59 78 111145 1074.21 108904 114498 60 74 109843 938.884 107156 111928 61 67 108823 1039.61 106055 111548 62 57 107595 911.344 105650 109873 63 73 106565 717.252 104323 108820 64 62 105981 891.732 103899 109912 65 73 105068 834.759 102613 107161 66 78 104049 796.307 102298 105983 67 77 103221 622.306 101504 104484 68 73 102443 794.3 99378.6 104484 69 74 101749 1155.43 97897.9 108187 70 71 100614 1188.58 97278.7 104646 71 72 99234.5 1425.43 95211.1 103105 72 69 97197 1454.65 94069.6 102040 73 68 95665 1119.02 92599.9 97950.3 74 77 94407.9 878.002 92131.7 96861.2 75 73 93487 678.177 92131.7 95604.5 76 70 92758.3 622.941 90809 94536.3 77 73 92008.2 597.618 90026 93689.8 78 63 91293 704.382 89147.4 93374.1 79 73 90589.3 623.21 89147.4 92264.6 80 80 89953 640.25 88761.3 92712.9 81 73 89301.6 440.816 87904.1 90659.8 82 78 88887.8 527.33 87591.1 90590.2 83 75 88306 560.258 86823 90620.2 84 66 87709.6 553.403 85709.8 89652.5 85 72 87120.5 631.092 85142.4 88409 86 73 86526.1 698.785 85072.6 88680.2 87 66 85839.1 570.544 84695.9 88142.9 88 72 85307.7 410.745 84674.7 86729.9 89 84 84980.7 429.339 83863.3 87034.8 90 76 84583.6 378.761 83430.7 85646.3 91 67 84151.6 377.758 83321.5 85721.9 92 82 83797.3 338.27 82929 84955.8 93 72 83517.6 298.837 82894.5 84508.1 94 73 83238.4 355.472 82144.2 84289.1 95 74 82941.2 368.185 81876.5 84481.8 96 84 82668.4 472.755 81423 84542.5 97 84 82205.7 480.11 81158.3 83573.7 98 72 81654.9 371.628 80876.3 82751.4 99 76 81314.4 284.989 80792.9 82222.2 100 75 81087.1 271.145 80401.1 81941.6 101 80 80856.4 260.979 79910.8 81897.6 102 73 80605.8 296.911 79880.3 81622.7 103 72 80275.7 312.174 79656.9 81516.6 104 69 79985 209.016 79599.7 81042.3 105 65 79843.6 185.629 79491.8 80803.8 106 67 79683.2 138.869 79174.3 80066.3 107 72 79548.6 209.701 78931.3 80455.4 108 65 79355.1 247.657 78600.4 80328.4 109 74 79134.9 235.616 78358.7 79874.4 110 80 78898.1 258.626 78353 79952 111 71 78693.8 267.205 78200.7 79940.6 112 71 78477.1 224.25 77806.3 79994.1 113 82 78300.2 302.094 77376.7 79363.9 114 69 78048.6 363.613 77057.4 79283 115 77 77776.3 334.823 77057.4 79055.7 116 74 77475.6 242.687 77038.7 78366.7 117 66 77288.8 190.993 76788.5 77973.1 118 74 77139.2 150.006 76809.8 77556.3 119 65 76996.5 152.493 76636.2 77696.2 120 70 76875.8 184.349 76515.9 77612.4 121 64 76701.6 147.436 76258.6 77256.3 122 79 76586.6 161.884 76006.8 77074.5 123 72 76487.9 194.356 76090 77170.2 124 72 76313.4 192.52 75795.2 76950.3 125 82 76157.7 169.537 75695.7 76903.8 126 62 76000.1 199.078 75550.9 76849.5 127 66 75810.6 202.931 75052.6 76422.2 128 80 75654.7 215.375 75096.2 76619.6 129 68 75478.1 173.679 74893.9 76119.7 130 79 75362.7 219.425 74801 76429.2 131 77 75205 219.06 74771.6 76286 132 76 75058.3 279.466 74624.6 76656.3 133 65 74857.9 166.29 74469.8 75460.9 134 81 74764.9 188.31 74182.4 75473.3 135 69 74651.2 285.532 74182.4 76137.3 136 73 74489 234.828 74120 75674.3 137 66 74262.4 159.169 73788.5 74847.5 138 64 74150.8 189.427 73788.5 75066 139 73 74011.4 171.374 73612.3 74491.6 140 73 73858.1 217.064 73367.4 75066.3 141 72 73716 190.081 73236.2 74548.1 142 69 73519.8 224.442 73090.9 74305.6 143 75 73349.5 204.245 72897.9 74222 144 68 73205.4 137.624 72897.9 73960.5 145 70 73127.3 134.154 72875.6 73798.4 146 67 73052 187.043 72675.9 73944.7 147 73 72917.9 173.014 72442.5 73767.7 148 75 72813.3 188.49 72367.4 73512.8 149 72 72676.1 173.453 72367.4 73293 150 82 72570.6 155.449 72316.2 73323.1 151 70 72451.3 92.5147 72274.3 72735.4 152 72 72383.3 99.189 72157.9 72922.9 153 70 72338.9 181.753 72083.2 73519 154 80 72253.5 129.241 71949.7 72677.2 155 55 72147.1 114.626 71946.4 72750.4 156 71 72081.8 106.28 71776 72359.4 157 75 71990 132.03 71728.4 72708.2 158 66 71914.5 134.683 71622.2 72623.5 159 70 71816.7 122.371 71622.2 72335.8 160 73 71755.3 146.342 71556.6 72323.4 161 78 71688.9 170.932 71531.7 72894.5 162 69 71623.7 98.9562 71477.6 72140.1 163 80 71605.1 140.025 71379.7 72096.9 164 58 71521.2 120.2 71271.6 71948.5 165 71 71455.5 144.596 71271.6 72456.3 166 78 71388.8 116.098 71219.2 71984.1 167 76 71332.2 92.2036 71188.1 71687.5 168 74 71287.6 122.409 71172.3 71961.9 169 72 71249.3 131.708 71110.3 72038.7 170 73 71205.3 103.642 71097.8 71647.4 171 68 71154.8 66.6753 71065.3 71480.5 172 71 71125.4 57.9064 71047 71392.2 173 63 71106 84.4234 71004.1 71674.8 174 77 71099.1 98.8177 70988.2 71691.4 175 63 71066.6 100.781 70982.3 71501.9 176 73 71038.3 75.1214 70927.1 71364 177 76 71001.8 64.062 70898.1 71302.5 178 69 70982.6 73.3536 70886.5 71364.4 179 73 70983.7 148.255 70880.6 72075.5 180 68 70959.2 184.416 70866.2 72628.8 181 75 70932.3 103.697 70840.1 71521.9 182 74 70905.6 95.5782 70816.1 71557.8 183 72 70876.5 73.4116 70776.7 71224.3 184 77 70866.4 107.68 70764.6 71630.4 185 74 70835.2 104.662 70723.5 71490.2 186 75 70802.2 78.905 70695.5 71182.1 187 78 70784.4 72.8487 70680.5 71125.9 188 71 70753.6 81.3207 70640.3 71210.4 189 68 70732.1 73.0191 70639.7 71092.3 190 76 70706.6 82.2918 70593.7 71063 191 72 70672.3 56.0005 70583.1 70859.4 192 72 70663.5 75.6656 70562.3 70996.7 193 79 70649.9 124.451 70541.2 71344.2 194 76 70613.4 74.9445 70530.6 70993.8 195 64 70578.1 43.1494 70519.3 70784.5 196 74 70580.3 108.2 70491.9 71175.5 197 71 70566.2 90.0578 70481 70971.4 198 72 70549.4 91.5063 70481 70972.6 199 65 70529.9 95.0373 70479.2 71083.6 200 71 70507.1 68.998 70462.6 71053.5 201 63 70502.8 62.0261 70458.3 70804.5 202 64 70481.5 59.8537 70448.6 71015.1 203 68 70492.3 78.1776 70431.7 70953.2 204 67 70506.1 135.181 70421 71212.4 205 74 70495.1 141.103 70405.4 71257.6 206 67 70463 93.4677 70372.2 70978.4 207 71 70478.7 148.901 70358.9 71294.5 208 63 70435.9 88.4804 70358.9 70876.3 209 72 70453.6 176.644 70344.1 71378.4 210 81 70402.3 61.9158 70335.7 70722.2 211 78 70406.8 96.9284 70326.8 71038 212 64 70371.5 51.5666 70326.8 70621.3 213 82 70372.5 102.967 70319.7 71094.3 214 85 70358.2 79.2217 70314.8 70757.1 215 80 70357.8 90.3797 70290.3 71017.2 216 75 70345.4 69.2047 70278.1 70835.6 217 67 70328.2 71.1097 70278.1 70756.2 218 56 70321.4 87.5809 70275.6 70783.5 219 59 70297.2 40.1293 70265.8 70480.1 220 66 70316 116.568 70261.3 70912.1 221 73 70306.1 122.749 70261.3 71417.2 222 75 70311.6 117.059 70259.6 71204.7 223 69 70292.4 91.1066 70256.3 71007.4 224 76 70307.2 94.4793 70233.6 70676.5 225 80 70285.7 72.6069 70233.6 70657.6 226 59 70285.8 88.9897 70228.3 70844.1 227 71 70259.3 41.1383 70228.3 70492.9 228 76 70266.4 85.3976 70225.7 70772.1 229 69 70258.3 70.1542 70220.5 70667.6 230 65 70252.5 72.5794 70217.5 70760.1 231 79 70262.8 111.637 70214.4 71016.7 232 67 70245.5 84.9601 70214.4 70809.4 233 74 70242.8 75.7363 70206.8 70740.1 234 73 70237.4 71.9866 70206.6 70626.5 235 75 70220.7 28.0785 70201.7 70381.2 236 78 70248.4 100.9 70201 70749.8 237 75 70232.4 62.4606 70198.7 70542.8 238 76 70239.7 111.486 70193.7 70947.8 239 62 70229.2 87.7483 70193.7 70750.6 240 81 70227.7 67.7814 70191.6 70620 241 68 70215.2 63.3902 70190.6 70606.7 242 77 70232.2 125.661 70186.3 71068.9 243 80 70209.3 60.1309 70184.5 70614.1 244 67 70233 116.073 70183.3 70815 245 78 70216.4 88.1889 70180 70717.9 246 68 70217.7 87.4873 70180 70664.7 247 64 70207.7 80.2926 70179.7 70738.9 248 73 70211.3 76.279 70177.4 70691.8 249 69 70209.9 80.6683 70173.1 70578.8 250 69 70201.7 64.395 70172.4 70573 251 80 70218.7 106.731 70172.4 70922.1 252 75 70200.4 59.7435 70172.4 70451.7 253 57 70193.5 48.3375 70172.4 70404.7 254 63 70210.8 112.838 70170.6 71026 255 72 70198.4 65.7852 70170.6 70554 256 66 70208.9 102.264 70169 70825.1 257 72 70197.2 64.227 70166.2 70593 258 70 70204.4 81.8356 70161.8 70588.5 259 74 70200.5 77.8624 70160.2 70626.8 260 77 70185.5 48.1584 70159.8 70432.8 261 68 70177.4 53.3886 70159.8 70608.4 262 66 70185 53.543 70158.4 70417.3 263 73 70193.6 77.5774 70157.8 70507.5 264 71 70186.5 73.4263 70156.4 70576.2 265 76 70183.9 77.1089 70154.8 70711 266 63 70190 85.8046 70154.8 70680.5 267 70 70179 66.7107 70153 70505.2 268 72 70180.2 77.6124 70152.5 70726.4 269 71 70176 63.3907 70152.5 70578 270 66 70186.9 76.3367 70152.5 70535.8 271 75 70178.8 65.465 70152 70576.9 272 73 70182.2 77.5626 70152 70584.8 273 79 70180.7 100.419 70151.8 70876.8 274 79 70186 91.7746 70151.7 70814.3 275 82 70179.5 79.6825 70151.2 70774.6 276 71 70177.4 78.8808 70151.2 70735.7 277 64 70189.2 104.864 70150.9 70815.5 278 70 70170.2 53.2298 70150.8 70462.6 279 75 70184.1 78.2798 70150.4 70551.3 280 74 70182 87.0134 70150.3 70749.7 281 64 70165.3 39.8781 70150.1 70372.3 282 71 70189.8 100.093 70150 70873.1 283 70 70172.8 68.9089 70149.7 70603.7 284 77 70179.9 75.0207 70149.2 70588.6 285 63 70179 70.9523 70149.1 70546.6 286 69 70160.4 42.3938 70148.4 70505.7 287 72 70191.6 114.214 70148.3 70843.9 288 64 70185.9 94.4247 70148.2 70713.1 289 76 70171.1 49.4649 70148 70404.1 290 77 70174.2 81.209 70147.8 70740.3 291 67 70162.5 52.1755 70147.6 70504.2 292 77 70181.8 90.3273 70147.6 70703 293 66 70175.8 71.7222 70147.5 70623.3 294 76 70172.7 83.8332 70147.5 70809.9 295 66 70163.7 62.6232 70147.4 70631.7 296 65 70170.8 70.717 70147.4 70603.7 297 71 70182.8 109.816 70147.3 70914.2 298 69 70182 74.2484 70147.3 70557.2 299 62 70177.8 102.13 70147.3 70833.6 300 71 70169.9 63.4763 70147 70611 301 78 70175.2 81.5195 70146.9 70664.4 302 67 70167.7 55.7338 70146.8 70400.9 303 70 70171.1 73.1817 70146.7 70501.6 304 80 70165.4 55.1772 70146.7 70469.6 305 73 70182.1 96.2128 70146.6 70857 306 79 70180.3 97.5704 70146.6 70770.3 307 78 70177.5 81.3719 70146.1 70727 308 69 70179.5 88.1878 70145.8 70655.1 309 69 70173.1 78.1677 70145.8 70643.7 310 77 70173.4 87.4665 70145.8 70692.8 311 74 70180.8 102.723 70145.6 70892.7 312 74 70160.5 55.3707 70145.6 70589.8 313 68 70169.9 55.1585 70145.5 70435 314 70 70173.2 83.9731 70145.4 70641.2 315 75 70193.4 115.977 70145.3 70722.4 316 76 70164.8 70.7144 70145.3 70625.5 317 72 70158.5 39.2913 70145.2 70375.3 318 79 70184.6 99.7582 70144.5 70700.2 319 68 70172.4 65.3169 70144.5 70519.6 320 75 70165.2 44.2132 70144.5 70374.8 321 66 70167.5 78.8841 70144.5 70726.1 322 70 70174.6 109.226 70144.4 71086.7 323 69 70168.8 59.9534 70144.4 70530.5 324 70 70175.6 123.278 70144.4 71219.4 325 76 70182.1 101.507 70144.4 70711.4 326 79 70171.4 59.0811 70144.3 70444.6 327 80 70178.9 90.5959 70144.3 70601.9 328 72 70172.7 81.3126 70143.6 70621.8 329 69 70174.5 77.6646 70143.8 70580 330 66 70168.7 75.1296 70143.8 70783 331 79 70164.7 55.3305 70143.5 70399.9 332 72 70164.3 59.4877 70143.4 70514 333 79 70176.8 104.497 70143.4 70782.9 334 61 70161.3 53.7778 70143.4 70495.6 335 74 70164.5 59.2257 70143.3 70535 336 69 70165.8 65.9298 70142.9 70569 337 70 70172.1 90.4292 70142.4 70745.4 338 57 70179.7 105.077 70142.3 70729 339 65 70169 92.1286 70142.3 70874.2 340 67 70183.3 81.7284 70142.3 70519.9 341 66 70169.1 84.7801 70142.3 70779 342 71 70161.3 61.8475 70141.3 70618.4 343 76 70180.5 138.229 70141.3 71353.7 344 69 70179.6 115.269 70141.3 71097.1 345 73 70160.4 49.8964 70141.2 70454.5 346 63 70167.5 90.8644 70141.1 70853.2 347 71 70165.8 68.8576 70141.1 70569.9 348 75 70185.1 107.069 70141 70760.1 349 73 70170.7 88.4331 70141 70667.8 350 73 70164.4 65.4151 70140.8 70553.3 351 81 70170.4 76.1254 70140.8 70690.2 352 65 70160.2 51.3828 70140.8 70463.3 353 73 70159.8 53.5856 70140.8 70447.1 354 71 70176.2 133.248 70140.8 71228.9 355 73 70173 101.976 70140.7 70716.3 356 78 70174.4 108.879 70140.7 70899.7 357 61 70160.8 60.047 70140.7 70633.5 358 71 70174.9 79.4858 70140.7 70613.6 359 70 70184.6 104.751 70140.7 70863.2 360 64 70163.3 73.5177 70140.7 70555 361 78 70172.9 91.9752 70140.7 70816.3 362 61 70158.7 51.6952 70140.7 70414.9 363 77 70179.1 99.4498 70140.7 70658.9 364 68 70163.2 69.8087 70140.7 70558.5 365 74 70164.6 72.6055 70140.7 70532.9 366 72 70172.5 104.073 70140.5 70988.8 367 79 70182.7 117.922 70140.4 70732.1 368 68 70171.8 102.916 70140.4 70820.1 369 68 70168.3 78.9463 70140.4 70723.3 370 70 70159 47.8688 70140.4 70360.6 371 65 70165.2 83.3526 70140.3 70636.1 372 73 70177.6 108.901 70140.3 70923.3 373 70 70180.9 121.652 70140.3 70963.3 374 74 70160.5 88.0695 70140.3 70875.2 375 72 70159.9 57.4204 70140.3 70528.5 376 69 70171.9 89.6398 70140.3 70822.2 377 73 70174.7 102.043 70140.3 70852.2 378 82 70154.2 39.1806 70140.3 70371 379 66 70172.7 96.397 70140.3 70742.4 380 70 70158.6 44.0683 70140.3 70369.2 381 71 70150.9 35.4356 70140.3 70335.2 382 73 70170.9 84.7181 70140.3 70592.6 383 70 70165.5 67.9163 70140.3 70543.7 384 67 70162.8 62.2117 70140.3 70528.6 385 66 70155 55.4458 70140.3 70532.6 386 74 70153.7 41.1611 70140.3 70430.2 387 71 70192 193.471 70140.2 71884.7 388 71 70163.2 86.4693 70140.2 70862.4 389 76 70177.3 104.086 70139.8 70909.2 390 63 70158.9 66.4424 70139.8 70708.6 391 74 70165.2 58.4125 70139.8 70463.7 392 71 70166.6 72.6105 70139.7 70578.1 393 77 70165.2 70.5379 70139.7 70605.2 394 71 70154.5 65.5804 70139.7 70620.5 395 67 70157.2 58.8287 70139.7 70526.3 396 78 70170.5 92.2738 70139.7 70720.5 397 76 70160.5 55.6704 70139.7 70474.6 398 70 70159.7 59.5535 70139.7 70476 399 67 70176 97.318 70139.7 70742.1 400 74 70171.8 100.138 70139.7 70750.5 401 76 70170.7 95.6087 70139.7 70921.9 402 75 70172.1 91.1559 70139.7 70663.7 403 73 70163.6 73.5009 70139.7 70536.5 404 75 70167.6 71.5575 70139.7 70610.2 405 77 70162.1 70.0827 70139.4 70660.3 406 64 70167.5 74.8699 70139.4 70616.1 407 64 70162.5 68.2297 70139.4 70690.2 408 73 70167.6 73.8299 70139.3 70525.2 409 82 70163.9 77.6862 70139.3 70675 410 79 70165.7 69.2492 70139.3 70529.6 411 72 70161.5 81.8093 70139.3 70828.1 412 67 70167.3 83.9043 70139.3 70677.2 413 64 70171 120.126 70139.3 71196.5 414 69 70162.6 62.9534 70139.3 70439.9 415 72 70168.4 75.7922 70139.3 70587.2 416 74 70158.1 57.8689 70139.3 70498.6 417 76 70162.2 61.3149 70139.3 70548.2 418 74 70158.4 54.1552 70139.2 70420 419 76 70161.6 60.9598 70139.2 70539.2 420 69 70178.8 101.035 70139.2 70717.7 421 72 70159.2 58.652 70139.2 70486.6 422 82 70169.6 75.3349 70139.2 70506.5 423 69 70156.9 47.8082 70139.1 70405.4 424 72 70169.2 82.9178 70139.1 70671 425 79 70164.6 70.9574 70139.1 70571.8 426 75 70170.1 84.3139 70139.1 70566 427 68 70164.4 76.8647 70139.1 70711.1 428 66 70166.9 68.5605 70139.1 70461.6 429 77 70171.8 100.954 70139.1 70995.8 430 70 70170.1 114.175 70139.1 71104.4 431 79 70156.6 56.1303 70139.1 70507.2 432 80 70168.3 71.1203 70139.1 70530.9 433 72 70161.7 51.8344 70139.1 70379.1 434 74 70162.7 63.2119 70139.1 70529.2 435 72 70170.2 79.2474 70139.1 70550.9 436 57 70170 117.051 70139.1 70894.2 437 77 70178.6 106.206 70139.1 70692.5 438 68 70166.3 62.5192 70139.1 70381.6 439 71 70169.4 95.6304 70139.1 70748.6 440 75 70160.8 66.1806 70139 70491.7 441 73 70167.7 76.0811 70139 70615.8 442 59 70152.4 39.3667 70139 70325.2 443 66 70163.6 71.5049 70139 70592.2 444 63 70161.6 55.1092 70139 70461.5 445 69 70152.8 51.5603 70139 70479.5 446 76 70152.9 50.242 70139 70423.7 447 60 70164.3 66.2144 70139 70533.9 448 66 70162.5 67.7063 70138.9 70506.5 449 69 70157.1 54.0561 70138.9 70506.4 450 71 70167.8 82.8506 70138.9 70615.1 451 66 70164.1 75.0728 70138.9 70605.9 452 77 70173.8 80.84 70138.9 70478.3 453 69 70163.9 78.3589 70138.8 70768.7 454 71 70161.1 57.8966 70138.8 70444.1 455 65 70163.9 74.8285 70138.8 70563.1 456 68 70171.1 99.3456 70138.7 70895.7 457 70 70164.7 59.543 70138.6 70435.6 458 66 70168.4 71.1733 70138.6 70485.5 459 85 70177.3 107.878 70138.5 70766.7 460 80 70175.1 82.5823 70138.5 70527.8 461 79 70167.9 73.2038 70138.5 70529.9 462 65 70165 76.2297 70138.5 70595.5 463 67 70161.8 69.0007 70138.5 70612.1 464 78 70169 74.4983 70138.5 70612.4 465 78 70174 85.96 70138.4 70622.3 466 76 70168.1 73.5738 70138.4 70585.2 467 78 70161.5 66.1374 70138.3 70514.1 468 69 70174 101.059 70138.3 70812.8 469 63 70164.7 77.2192 70138.3 70660.4 470 68 70165.6 73.6991 70138.3 70543.8 471 70 70156 54.896 70138.3 70527.5 472 75 70154.9 53.6431 70138.3 70502.5 473 76 70159.3 52.6088 70138.3 70402.6 474 69 70162.7 70.681 70138.3 70543.8 475 70 70160 54.3368 70138.2 70403.1 476 70 70167.9 100.083 70138.2 70840 477 78 70167.9 77.7815 70138.2 70560.9 478 69 70175.3 113.13 70138.2 70839.7 479 72 70154.8 49.6274 70138.2 70452.2 480 64 70170.9 124.539 70138.2 71173.1 481 58 70178.5 103.23 70138.2 70680.4 482 73 70183.3 111.95 70138.1 70851.6 483 67 70158.1 62.0517 70138.1 70556 484 67 70163.4 78.74 70138.1 70761.8 485 66 70167.1 65.9573 70138.1 70406 486 68 70158.1 54.1677 70138.1 70404.1 487 69 70159.6 58.593 70138.1 70473.9 488 62 70168.6 84.0944 70138.1 70581.4 489 66 70153.8 47.2283 70138.1 70445.5 490 71 70162.7 61.4336 70138.1 70496.2 491 73 70154.1 51.0064 70138 70470.5 492 85 70179 104.541 70138.1 70682.7 493 77 70163.2 75.2092 70138 70577.2 494 65 70165.2 64.8869 70138 70472 495 81 70174.4 92.7472 70138 70793.5 496 76 70160 57.9795 70138 70432.5 497 82 70159.3 62.8622 70138 70493.6 498 64 70161.4 56.9193 70138 70410.8 499 74 70171 81.091 70138 70624.7 500 71 70147.5 30.5645 70138 70337.3 501 69 70181.6 91.4833 70138 70510.8 502 67 70160.9 69.7878 70137.9 70671.8 503 73 70160.4 61.6774 70137.9 70494.3 504 77 70167.7 68.7897 70137.9 70548.5 505 80 70161.1 55.353 70137.9 70525.6 506 73 70162.3 58.8107 70137.9 70480.7 507 73 70153.4 56.174 70137.9 70527.7 508 75 70160.7 65.2258 70137.9 70497.5 509 67 70164.9 79.6173 70137.9 70653.9 510 72 70155.3 45.985 70137.9 70380 511 60 70159.9 64.4995 70137.9 70431.5 512 80 70168.8 89.3191 70137.9 70726.3 513 69 70153.1 52.9749 70137.9 70549.1 514 76 70158 55.3327 70137.9 70421.8 515 75 70174.8 81.9963 70137.9 70504.2 516 72 70165.4 70.5034 70137.8 70644.9 517 70 70162.7 105.336 70137.8 71070.8 518 75 70174.2 123.294 70137.8 71005.2 519 77 70175 90.6284 70137.8 70647.2 520 71 70173.4 105.694 70137.8 70924 521 66 70157.8 62.1638 70137.8 70500.7 522 66 70173.2 79.6981 70137.8 70521 523 70 70154.5 55.7756 70137.8 70499.8 524 76 70173.7 148.339 70137.8 71499.2 525 83 70165.5 71.1112 70137.7 70602 526 78 70169.9 88.1529 70137.7 70624.8 527 75 70170.1 111.268 70137.7 70993.8 528 77 70163.9 66.0686 70137.7 70527.8 529 67 70161 79.8608 70137.6 70794.6 530 80 70162.2 80.2539 70137.6 70700.7 531 53 70155.6 62.3289 70137.6 70522.3 532 77 70158.8 54.6457 70137.6 70415 533 77 70177.4 97.731 70137.5 70842.5 534 69 70165.5 68.1788 70137.5 70523.1 535 72 70172.5 91.5547 70137.5 70632.3 536 77 70179.1 99.1261 70137.5 70807.8 537 81 70186.3 118.671 70137.5 70866.8 538 80 70177.5 89.9349 70137.5 70609.4 539 74 70159.3 62.6763 70137.5 70544 540 82 70174.3 82.0334 70137.4 70584.6 541 74 70162.9 68.0124 70137.4 70606.9 542 72 70172.3 104.3 70137.4 70732.6 543 70 70162.8 78.7445 70137.4 70622.2 544 76 70154.7 38.3628 70137.4 70322.2 545 78 70165.5 68.1223 70137.4 70509.9 546 70 70159.8 60.0244 70137.4 70536 547 65 70164.8 89.169 70137.4 70774.7 548 82 70158.2 51.5465 70137.4 70372.3 549 74 70158.3 72.9944 70137.4 70716.2 550 69 70157.2 61.2997 70137.4 70582.4 551 76 70166.8 81.4973 70137.4 70612.1 552 73 70190.2 123.798 70137.4 70859.4 553 68 70149.4 41.7843 70137.4 70431.7 554 83 70172.5 83.294 70137.4 70546.3 555 78 70183.4 142.003 70137.4 71252.5 556 68 70161.1 77.5149 70137.3 70743.2 557 74 70170.3 92.5275 70137.3 70757.6 558 74 70159.6 59.6096 70137.3 70494.6 559 76 70157.2 49.9222 70137.3 70396.8 560 76 70159.3 81.1824 70137.3 70747 561 70 70149.4 45.9094 70137.3 70455.6 562 71 70172.9 84.6141 70137.3 70615.7 563 66 70181.4 163.188 70137.3 71586.9 564 67 70160.3 72.9573 70137.2 70619.9 565 77 70178.9 102.906 70137.2 70704.7 566 76 70160.4 66.0117 70137.2 70534.3 567 74 70166.1 85.5509 70137.2 70746.7 568 69 70172.3 85.5886 70137.2 70706.7 569 75 70154.2 56.3759 70137.2 70474.2 570 79 70176.1 101.626 70137.2 70618.7 571 75 70170.9 81.4996 70137.2 70621.2 572 66 70157.6 61.2941 70137.2 70416.1 573 61 70153.7 53.1772 70137.2 70439.9 574 73 70166.1 69.5667 70137.1 70553.4 575 74 70165.4 94.298 70137.1 70805 576 70 70172.7 121.495 70137.1 71169.1 577 77 70157.4 60.8174 70137.1 70449.5 578 72 70179.5 90.7561 70137.1 70581.2 579 70 70155.6 64.2027 70137.1 70660.7 580 63 70190.3 126.599 70137.1 70823.8 581 76 70167 64.8967 70137.1 70494.5 582 64 70161.8 80.2849 70137.1 70656.6 583 78 70164.7 72.4289 70137.1 70613.9 584 74 70166 96.7491 70137.1 70977 585 71 70167.6 83.6529 70137.1 70645.8 586 77 70161.4 57.9884 70137.1 70470.4 587 70 70164.3 77.8791 70137.1 70536.2 588 70 70164.1 62.8292 70137.1 70404.5 589 67 70152.8 66.3735 70137.1 70729.1 590 77 70162.9 62.3333 70137.1 70469.2 591 73 70170.5 79.7733 70137.1 70573.3 592 68 70160.4 65.0277 70137.1 70497.6 593 69 70159.6 65.2639 70137.1 70581.2 594 67 70163.6 67.189 70137.1 70498.7 595 68 70160 58.9842 70137.1 70435.6 596 70 70153.8 45.5843 70137.1 70385.4 597 80 70177.3 95.0123 70137.1 70702.6 598 80 70155.6 43.7745 70137.1 70344 599 63 70172.2 98.1005 70137.1 70604.6 600 74 70171.2 76.9045 70137.1 70533.7 601 67 70153.9 49.3591 70137.1 70476.2 602 73 70164 72.3833 70137.1 70594 603 70 70171.5 97.4033 70137 70764.2 604 67 70168.5 89.2462 70137 70702.4 605 76 70171.4 92.7616 70137 70689.8 606 75 70168.5 72.6218 70137 70496 607 71 70164.9 79.5708 70137 70587.4 608 62 70173.7 128.041 70137 70953.9 609 71 70170 94.2697 70137 70715.2 610 66 70154.3 80.026 70137 70783.7 611 78 70156.2 46.9764 70137 70366.7 612 75 70169.5 92.3331 70137 70780 613 66 70157.3 81.7952 70137 70768.1 614 73 70171.1 78.9997 70137 70621.7 615 66 70166.6 97.2714 70137 70945.8 616 71 70173.4 96.8556 70137 70678.3 617 67 70172.2 78.8178 70137 70531 618 74 70170.8 93.8532 70137 70782.3 619 73 70161.2 69.662 70137 70560.3 620 70 70150.9 41.2179 70137 70416.9 621 76 70163.9 71.8269 70137 70484.9 622 77 70173.1 83.3165 70137 70543.3 623 78 70161.5 78.9345 70137 70774.4 624 72 70157.4 60.976 70137 70585 625 76 70172.2 83.2434 70137 70587.8 626 75 70169.4 87.5514 70137 70805 627 70 70157 57.9351 70137 70426.6 628 76 70153.9 47.9487 70137 70441.3 629 64 70174.7 99.6321 70137 70755.3 630 66 70166.2 119.641 70137 71217.9 631 59 70160.2 63.9956 70137 70456.3 632 84 70167.9 86.7756 70137 70674.3 633 62 70163.7 82.3981 70137 70739 634 61 70165.5 84.7712 70137 70648.7 635 77 70163 71.856 70136.9 70635.9 636 71 70164.7 86.3632 70136.9 70686.5 637 75 70171.8 107.958 70136.9 70786.4 638 71 70165.2 95.3845 70136.9 70941.2 639 82 70162 77.3401 70136.9 70602.6 640 73 70169.5 95.7018 70136.9 70869.6 641 67 70163.9 72.1967 70136.9 70617.3 642 67 70159.8 76.9098 70136.9 70742.3 643 80 70174.3 90.3211 70136.9 70646 644 75 70157.4 69.4626 70136.9 70554.7 645 71 70169.5 90.1247 70136.9 70736 646 70 70169.7 104.93 70136.9 70861.5 647 76 70172.7 92.9982 70136.9 70698.5 648 74 70153.5 42.3523 70136.9 70332.5 649 67 70157.6 63.0184 70136.9 70583.2 650 79 70170.2 98.5144 70136.9 70795 651 72 70170.7 80.6287 70136.9 70594.6 652 67 70160.6 71.7557 70136.9 70668.3 653 69 70156.8 52.5452 70136.9 70435.3 654 70 70156.8 55.6613 70136.9 70465.7 655 79 70183.1 113.046 70136.9 70701.7 656 81 70165.7 64.4007 70136.9 70406.5 657 78 70178.6 116.684 70136.9 71048.3 658 73 70157.4 66.3146 70136.9 70568.1 659 71 70187.3 122.195 70136.8 71051.1 660 71 70163.3 85.1603 70136.8 70747 661 72 70168.7 108.262 70136.8 70902.5 662 71 70159.6 65.2315 70136.8 70543.7 663 74 70165.4 77.5646 70136.8 70536.8 664 76 70163.2 78.2632 70136.8 70551.5 665 73 70191.4 127.288 70136.8 70797.3 666 70 70174.6 97.9403 70136.8 70688 667 74 70159.9 70.4017 70136.7 70698.4 668 78 70160.4 62.5995 70136.7 70463.2 669 72 70168.5 89.5142 70136.7 70644.8 670 73 70170.1 75.6184 70136.7 70495.6 671 85 70178.3 111.323 70136.7 70918.5 672 75 70165 76.8132 70136.7 70612.7 673 71 70155.9 55.7982 70136.7 70448.6 674 69 70163.6 97.2091 70136.7 70970.7 675 63 70166 90.2744 70136.7 70646 676 72 70157.4 51.5447 70136.7 70408 677 71 70152.9 43.4944 70136.7 70356.1 678 72 70159.3 61.1652 70136.7 70475.2 679 76 70182.9 109.385 70136.7 70658.4 680 74 70160.1 56.9964 70136.7 70440.9 681 76 70185.6 100.239 70136.7 70533.9 682 72 70164.3 69.3518 70136.6 70532.7 683 65 70156.3 76.6309 70136.6 70783.9 684 64 70149.1 42.6181 70136.6 70469.4 685 77 70157.6 67.4755 70136.6 70544.7 686 74 70183.5 140.68 70136.6 71071.4 687 75 70170.9 95.6097 70136.6 70811.5 688 80 70165.8 112.034 70136.6 70898.5 689 62 70163.3 83.5683 70136.6 70708.1 690 78 70163.2 72.7427 70136.6 70580 691 66 70157.6 62.9074 70136.6 70579.4 692 70 70155.6 62.7523 70136.6 70598.9 693 69 70158.6 64.298 70136.6 70509.3 694 82 70181 113.585 70136.6 70924 695 77 70167.3 78.6337 70136.5 70589.9 696 65 70182.4 122.777 70136.5 70962.1 697 73 70174.5 94.5763 70136.5 70808.9 698 73 70157.2 71.1066 70136.5 70687.3 699 76 70158.8 56.6713 70136.5 70497.9 700 69 70165.3 71.9741 70136.5 70593.9 Initial Control Points: [[ 31.08594556 97.07734616] [151.35331041 128.65688751] [104.42517432 119.65593613] [225.76573489 138.73866387]]
In [7]:
Copied!
initial_fitted_curve_points, initial_t_values = generate_curve_within_boundaries(initial_control_points, image.shape, 0.01)
plt.figure(figsize=(10, 5))
plt.subplot(1, 3, 1)
plt.plot(data_points[:, 0], data_points[:, 1], 'b.')
plt.title('Data Points')
plt.imshow(image, cmap='gray', alpha=0, origin='lower')
plt.subplot(1, 3, 2)
plt.plot(data_points[:, 0], data_points[:, 1], 'b.')
plt.plot(initial_fitted_curve_points[:, 0], initial_fitted_curve_points[:, 1], 'r-')
plt.plot(initial_control_points[:, 0], initial_control_points[:, 1], 'g.')
plt.title('Initial Fitted Curve')
plt.imshow(image, cmap='gray', alpha=0, origin='lower')
plt.subplot(1, 3, 3)
plt.imshow(image, cmap='gray', origin='lower')
plt.plot(initial_fitted_curve_points[:, 0], initial_fitted_curve_points[:, 1], 'r-')
plt.plot(initial_control_points[:, 0], initial_control_points[:, 1], 'g.')
plt.title('Initial Fitted Curve on Image')
plt.axis('off')
plt.tight_layout()
plt.show()
initial_fitted_curve_points, initial_t_values = generate_curve_within_boundaries(initial_control_points, image.shape, 0.01)
plt.figure(figsize=(10, 5))
plt.subplot(1, 3, 1)
plt.plot(data_points[:, 0], data_points[:, 1], 'b.')
plt.title('Data Points')
plt.imshow(image, cmap='gray', alpha=0, origin='lower')
plt.subplot(1, 3, 2)
plt.plot(data_points[:, 0], data_points[:, 1], 'b.')
plt.plot(initial_fitted_curve_points[:, 0], initial_fitted_curve_points[:, 1], 'r-')
plt.plot(initial_control_points[:, 0], initial_control_points[:, 1], 'g.')
plt.title('Initial Fitted Curve')
plt.imshow(image, cmap='gray', alpha=0, origin='lower')
plt.subplot(1, 3, 3)
plt.imshow(image, cmap='gray', origin='lower')
plt.plot(initial_fitted_curve_points[:, 0], initial_fitted_curve_points[:, 1], 'r-')
plt.plot(initial_control_points[:, 0], initial_control_points[:, 1], 'g.')
plt.title('Initial Fitted Curve on Image')
plt.axis('off')
plt.tight_layout()
plt.show()
In [8]:
Copied!
pixel_size = 2.170
penalty_threshold = 0.05
dither_range = 70
iterations = 700
pixel_size = 2.170
penalty_threshold = 0.05
dither_range = 70
iterations = 700
In [9]:
Copied!
ga_refine = GA_Refine(image, pixel_size=pixel_size, penalty_threshold=penalty_threshold, dithering_range=dither_range, iterations=iterations)
refined_control_points = ga_refine(initial_control_points, image)
refined_control_points = np.array(refined_control_points)
print(f'Control points: {refined_control_points}')
ga_refine = GA_Refine(image, pixel_size=pixel_size, penalty_threshold=penalty_threshold, dithering_range=dither_range, iterations=iterations)
refined_control_points = ga_refine(initial_control_points, image)
refined_control_points = np.array(refined_control_points)
print(f'Control points: {refined_control_points}')
/Users/zhenhuang/opt/anaconda3/envs/cryoem/lib/python3.9/site-packages/deap/creator.py:185: RuntimeWarning: A class named 'Individual' has already been created and it will be overwritten. Consider deleting previous creation of that class or rename it.
warnings.warn("A class named '{0}' has already been created and it "
gen nevals avg std min max 0 31 -2.19361e+11 6.37697e+11 -3.19596e+12 2268.36 1 20 1807.92 553.378 666.08 2268.36 2 26 2253.73 21.2737 2196.9 2285.43 3 26 2271.17 14.0597 2230.47 2303.17 4 22 2286.59 21.6743 2227.5 2327.79 5 26 2311.58 17.6472 2247.04 2330.3 6 22 2326.56 2.99077 2316.32 2330.36 7 28 2328.49 2.05926 2321.55 2331.87 8 28 2330.2 2.04154 2322.19 2333.5 9 26 2330.65 2.82643 2320.42 2333.5 10 21 2332.66 1.96151 2327.78 2339.26 11 22 2334.27 2.51638 2326.51 2340.58 12 23 2336.4 3.57707 2326.09 2340.58 13 24 2338.03 3.76101 2325.24 2341.43 14 22 2339.96 2.04535 2331.77 2341.43 15 22 2340.57 1.30779 2336.45 2341.62 16 25 2340.98 1.22573 2336.78 2342.17 17 24 2341.68 0.236803 2341.17 2342.17 18 25 2341.51 0.918996 2338.22 2342.17 19 22 2341.83 0.807699 2337.95 2342.17 20 20 2341.86 1.07593 2336.86 2342.51 21 19 2341.96 0.834305 2338.39 2342.54 22 20 2342.18 0.65352 2339.58 2342.54 23 23 2342.35 0.682342 2338.68 2342.56 24 27 2342.02 1.54013 2334.68 2342.59 25 22 2341.33 3.57911 2324.25 2342.68 26 23 2341.4 2.44796 2334.03 2342.86 27 26 2342.25 1.30374 2335.86 2342.96 28 24 2341.7 3.08845 2327.03 2343.05 29 24 2342.2 1.93555 2332.4 2343.5 30 23 2342.49 1.33994 2336.91 2343.5 31 25 2342.52 2.12565 2333.02 2343.67 32 22 2342.51 2.44694 2331.66 2343.91 33 26 2342.97 1.39644 2338.92 2344.23 34 23 2343.72 0.570253 2341.75 2344.51 35 23 2344.05 0.467681 2342.75 2344.75 36 21 2342.88 4.17304 2324.21 2344.75 37 23 2344.13 1.15838 2339.8 2344.75 38 23 2344.17 1.3303 2338.75 2344.77 39 24 2344.16 1.91597 2335.47 2344.77 40 24 2344.05 2.61882 2330.3 2344.8 41 23 2344.01 2.20027 2334 2344.86 42 25 2344.12 1.56191 2337.61 2344.86 43 21 2344.38 1.28652 2337.92 2344.92 44 23 2344.34 1.26586 2339.02 2345.24 45 26 2344.48 1.19441 2338.64 2345.32 46 22 2344.58 1.07059 2341.51 2345.32 47 20 2344.27 2.71701 2335.23 2345.32 48 24 2344.37 2.74581 2332.03 2345.37 49 22 2344.98 1.22507 2338.48 2345.39 50 25 2344.21 3.20105 2328.02 2345.46 51 24 2344.62 1.30323 2341.61 2345.46 52 24 2345.11 0.740261 2342.37 2345.49 53 20 2344.78 1.37462 2339.73 2345.71 54 24 2345.21 0.693663 2342.41 2345.71 55 21 2345.27 0.823774 2342.66 2345.71 56 21 2345.25 0.976412 2342.31 2345.78 57 26 2345.14 1.66991 2337.19 2345.78 58 20 2345.54 0.869234 2340.89 2345.94 59 21 2345.09 1.78127 2338.27 2346.16 60 25 2344.43 6.36313 2309.98 2346.16 61 27 2345.62 1.01332 2341.41 2346.16 62 23 2345.83 0.646896 2343.29 2346.26 63 17 2345.62 1.82049 2337.85 2346.26 64 23 2345.78 1.5271 2337.66 2346.37 65 18 2346 0.653186 2343.7 2346.37 66 27 2344.82 2.37379 2337.34 2346.49 67 27 2346.08 1.03714 2340.64 2346.61 68 24 2345.85 1.72215 2337.28 2346.64 69 24 2346.05 1.88719 2336.15 2346.87 70 23 2346.34 0.84763 2343.59 2346.87 71 25 2346.56 0.729116 2343.66 2346.9 72 24 2346.53 0.914584 2342.32 2346.9 73 24 2346.74 0.58324 2343.77 2347.02 74 22 2346.67 0.637278 2344.03 2347.02 75 28 2346.82 0.367388 2345.07 2347.08 76 22 2346.7 0.618555 2344.71 2347.16 77 27 2346.09 2.04481 2338.89 2347.33 78 23 2346.41 1.58224 2340.61 2347.4 79 25 2346.18 2.48222 2336.92 2347.41 80 24 2346.79 0.924359 2344.19 2347.46 81 25 2346.84 1.72536 2339.43 2347.43 82 23 2346.67 1.4363 2342.07 2347.46 83 27 2347.1 1.09479 2341.91 2347.69 84 24 2346.48 1.85259 2341.53 2347.69 85 23 2347.45 0.425707 2345.74 2348.24 86 22 2346.95 2.07363 2340.03 2348.77 87 24 2347.81 0.578239 2346.09 2348.77 88 21 2348.19 0.600388 2345.76 2348.9 89 23 2348.31 0.978953 2344.18 2348.9 90 24 2348.02 1.73059 2341.03 2348.96 91 27 2348.29 2.08712 2337.08 2348.98 92 23 2348.42 1.99271 2337.72 2349.14 93 23 2348.4 1.52927 2342.27 2349.27 94 18 2348.54 1.62455 2340.38 2349.2 95 19 2348.71 0.899674 2345.04 2349.26 96 25 2348.91 0.909627 2345.14 2349.44 97 26 2348.94 0.90263 2345.23 2349.6 98 18 2349.34 0.288851 2348.31 2349.6 99 28 2348.85 1.85424 2339.98 2349.6 100 25 2349.16 1.39092 2342.1 2349.75 101 21 2349.35 0.723362 2346.51 2349.81 102 25 2349.29 1.23939 2344.58 2349.83 103 24 2349.11 2.03808 2339.94 2349.96 104 20 2349.53 0.696406 2347.26 2350.01 105 21 2349.31 1.94314 2340.48 2350.09 106 23 2349.6 1.38158 2342.56 2350.34 107 24 2348.69 5.42149 2322.31 2350.35 108 29 2349.4 1.69555 2344.02 2350.44 109 23 2349 3.16637 2338.03 2350.67 110 20 2349.61 3.29694 2332.42 2350.67 111 23 2349.71 2.15072 2339.53 2350.73 112 27 2350.24 1.22348 2344.25 2350.8 113 23 2350.3 1.01242 2347.23 2350.95 114 23 2350.45 1.14609 2345.44 2350.95 115 24 2350.08 1.92183 2342.32 2350.99 116 25 2349.97 2.58454 2340.64 2351.17 117 22 2350.47 1.51112 2343.29 2351.17 118 23 2350.42 1.93679 2340.15 2351.14 119 18 2350.95 0.402603 2349.05 2351.14 120 28 2350.73 0.977932 2346.99 2351.27 121 22 2350.11 3.0736 2336.08 2351.27 122 19 2350.75 0.826592 2348.33 2351.29 123 27 2350.37 1.70712 2343.86 2351.35 124 26 2350.12 3.79407 2330.1 2351.38 125 21 2350.9 1.42446 2343.65 2351.38 126 23 2351.01 0.883519 2347.71 2351.38 127 26 2350.97 0.913061 2347.12 2351.43 128 21 2350.57 1.93939 2341.92 2351.47 129 25 2350.61 1.32524 2346.13 2351.47 130 19 2351.22 0.550225 2349.14 2351.5 131 22 2351.08 1.72998 2341.71 2351.54 132 21 2350.62 1.80806 2342.93 2351.65 133 25 2350.11 4.81842 2327.29 2351.67 134 22 2351.15 0.938377 2347.13 2351.67 135 24 2351.05 1.65222 2343.79 2351.79 136 19 2351.09 1.59101 2343.92 2351.83 137 24 2351.04 1.77843 2343.32 2352.26 138 22 2351.64 0.522271 2349.49 2352.29 139 27 2350.73 4.5204 2326.87 2352.29 140 27 2351.33 2.75636 2337.12 2352.42 141 26 2351.84 1.22212 2345.6 2352.3 142 22 2351.66 1.14215 2347.49 2352.3 143 23 2351.49 1.92837 2343.13 2352.4 144 23 2352.03 0.743291 2348.49 2352.4 145 22 2351.09 4.70626 2326.08 2352.4 146 21 2352.32 0.336056 2350.49 2352.48 147 24 2351.95 1.83798 2342.48 2352.49 148 26 2351.9 1.3392 2346.14 2352.5 149 21 2352.07 1.54985 2343.85 2352.52 150 19 2351.77 1.84518 2342.98 2352.76 151 27 2351.39 3.85921 2335.2 2352.54 152 21 2352.35 0.474929 2350.26 2352.56 153 21 2352.27 0.879779 2347.7 2352.56 154 20 2352.41 0.458361 2350.21 2352.65 155 22 2351.95 1.45031 2346.91 2352.65 156 18 2352.19 1.60552 2343.77 2352.7 157 28 2351.9 1.70374 2346.38 2352.73 158 26 2352.42 0.941341 2348.29 2352.98 159 18 2352.04 2.34972 2340.14 2352.98 160 26 2352.38 1.73602 2343.34 2353.36 161 23 2352.08 2.85988 2341.72 2353.36 162 24 2351.99 2.95927 2342.18 2353.36 163 24 2352 3.449 2336.36 2353.44 164 26 2352.59 2.03448 2342.8 2353.43 165 23 2352.92 1.47636 2346.18 2353.51 166 24 2353.12 0.792561 2350.17 2353.56 167 24 2352.77 2.29677 2341.53 2353.56 168 23 2353.4 0.324405 2351.8 2353.59 169 24 2353.05 1.37975 2346.75 2353.91 170 28 2352.53 2.33472 2343.46 2353.92 171 27 2353.2 1.29429 2348.35 2353.92 172 19 2353.06 2.02003 2344.43 2353.92 173 23 2353.1 1.72035 2347.23 2354.02 174 23 2353.75 0.587538 2350.84 2354.02 175 19 2352.57 3.69566 2336.5 2354.02 176 18 2353.85 0.387261 2352.26 2354.03 177 23 2353.52 2.00383 2342.7 2354.03 178 25 2352.55 5.84633 2321.56 2354.21 179 18 2353.88 0.533175 2351.67 2354.21 180 25 2353.72 0.737473 2350.74 2354.3 181 25 2353.09 2.19061 2343.12 2354.3 182 21 2353.96 0.8613 2350.6 2354.3 183 18 2353.78 0.979832 2350.64 2354.34 184 23 2353.8 1.67971 2345.34 2354.35 185 28 2353.62 1.84441 2344.07 2354.59 186 28 2353.29 3.38665 2335.4 2354.59 187 25 2354.19 0.769999 2350.23 2354.59 188 28 2353.71 3.37142 2335.58 2354.7 189 18 2354.13 1.29605 2347.7 2354.75 190 27 2354.01 1.91619 2346.27 2354.77 191 27 2354.34 0.898417 2350.55 2354.79 192 20 2353.57 3.34831 2341.01 2354.81 193 22 2353.44 4.28453 2336.85 2354.97 194 19 2354.69 0.448941 2352.42 2354.97 195 25 2354.08 2.29195 2346.77 2355.31 196 20 2354.58 1.58677 2346.21 2355.31 197 21 2354.92 0.629982 2351.77 2355.6 198 22 2355.16 0.754817 2351.27 2355.6 199 23 2354.47 2.35243 2343.38 2355.96 200 24 2354.64 2.009 2346.23 2356.14 201 16 2355.81 0.350557 2354.2 2356.14 202 23 2355.65 0.932627 2352.84 2356.14 203 27 2355.08 2.2461 2345.18 2356.29 204 22 2355.67 0.858154 2353.43 2356.29 205 23 2355.92 0.6859 2353.82 2356.29 206 23 2355.79 1.63804 2347.31 2356.47 207 17 2356.09 0.897211 2351.63 2356.47 208 22 2356.32 0.324151 2354.77 2356.51 209 23 2356.38 0.388148 2354.63 2356.61 210 23 2356.17 1.3004 2349.62 2356.61 211 21 2355.9 2.84562 2340.55 2356.63 212 23 2355.8 3.17191 2338.79 2356.63 213 22 2356.22 1.10633 2351.32 2356.63 214 21 2356.27 0.867002 2353.3 2356.8 215 25 2356.39 0.913014 2352.2 2357.05 216 21 2356.44 0.771046 2353.77 2357.22 217 27 2356.03 2.1043 2345.81 2357.09 218 28 2356.61 0.946014 2352.06 2357.09 219 23 2356.54 0.947815 2354 2357.24 220 22 2356.44 1.26678 2352.75 2357.24 221 25 2356.8 0.956379 2351.97 2357.31 222 26 2356.47 1.80673 2348.86 2357.36 223 20 2356.96 0.958404 2352.72 2357.74 224 27 2357.04 1.56964 2348.55 2357.74 225 24 2357.17 0.854546 2353.46 2357.91 226 24 2356.29 5.25253 2328.4 2357.97 227 19 2356.68 3.33593 2340.97 2358.04 228 20 2357.18 2.21888 2345.73 2358.06 229 23 2357.05 2.67136 2345.09 2358.06 230 25 2357.01 2.96339 2342.24 2358.29 231 24 2357.75 1.69786 2348.92 2358.57 232 19 2358.07 0.737747 2354.91 2358.57 233 23 2358.27 0.588425 2356.03 2358.6 234 23 2358.25 0.725806 2355.88 2358.98 235 18 2358.01 2.05874 2348.56 2358.98 236 22 2358.32 0.870523 2354.63 2359.01 237 27 2358.27 1.48802 2351.29 2359.11 238 23 2358.58 1.47135 2351.08 2359.15 239 24 2357.16 3.96455 2345.86 2359.21 240 25 2358.97 0.666878 2355.52 2359.43 241 26 2358.65 1.63274 2351.16 2359.43 242 22 2359.05 0.535297 2356.34 2359.28 243 22 2358.64 1.93264 2349.8 2359.76 244 22 2358.75 1.58771 2352.46 2359.78 245 21 2358.66 2.38393 2346.99 2359.78 246 17 2359.31 1.16683 2354.19 2359.84 247 26 2358.75 2.33518 2350.04 2359.84 248 23 2359.02 3.81047 2338.2 2359.84 249 23 2359.69 0.428475 2358.04 2360.05 250 27 2357.93 3.46758 2347.87 2360.05 251 25 2359.33 1.65877 2352.6 2360.05 252 22 2359.59 0.803414 2355.79 2360.14 253 23 2359.07 3.90766 2338.17 2360.44 254 20 2358.62 4.25333 2339.02 2360.44 255 26 2359.43 1.29274 2354.92 2360.55 256 25 2359.69 2.01309 2351.32 2361.3 257 27 2360.02 1.87631 2351.88 2361.67 258 22 2360.86 0.679494 2358.56 2361.67 259 20 2360.49 4.04474 2338.69 2361.84 260 24 2360.82 2.34551 2350.2 2361.84 261 22 2361.51 0.731675 2357.83 2361.98 262 24 2361.26 2.06404 2350.53 2361.99 263 22 2361.73 0.600404 2358.68 2361.99 264 27 2360.15 5.26231 2334.64 2362.03 265 19 2361.6 1.41355 2356.01 2362.04 266 25 2361.37 1.75791 2354.2 2362.04 267 21 2361.71 0.801087 2358.08 2362.05 268 24 2361.05 2.96139 2347.44 2362.73 269 27 2361.15 2.9927 2346.17 2362.73 270 20 2362.11 0.57206 2359.7 2362.73 271 27 2361.69 2.56143 2351.63 2362.89 272 23 2360.54 9.34858 2309.7 2362.89 273 24 2362.2 1.49778 2355.53 2363.06 274 22 2362.62 0.97312 2357.68 2363.29 275 23 2362.61 0.998282 2357.88 2363.29 276 25 2362.88 0.815393 2359.68 2363.68 277 21 2362.45 2.57794 2349.79 2363.68 278 24 2362.37 2.52927 2352.71 2363.73 279 24 2363.25 0.880535 2359.73 2363.73 280 21 2362.8 2.29814 2354.31 2363.73 281 21 2363.31 1.15327 2358.45 2364.09 282 24 2363.22 2.14342 2352.24 2364.09 283 28 2362.51 2.47248 2353.96 2364.1 284 25 2363.12 2.83356 2348.8 2364.26 285 30 2362.88 2.85983 2351.74 2364.33 286 21 2363.98 0.630375 2361.64 2364.9 287 26 2364.23 0.750658 2360.71 2365.2 288 27 2363.95 1.7869 2357.1 2365.22 289 26 2364.49 1.49869 2358.08 2365.43 290 27 2364.04 2.87363 2351.65 2365.62 291 25 2364.9 1.24408 2359.02 2365.76 292 23 2365.27 0.633756 2363.18 2365.76 293 17 2365.13 1.67937 2356.51 2365.76 294 25 2365.1 1.36068 2358.41 2365.76 295 25 2364.83 2.92276 2350.12 2365.83 296 19 2364.68 3.87238 2344.22 2365.83 297 24 2364.43 5.09476 2341.65 2366 298 22 2365.51 0.923915 2361.62 2366 299 25 2365.39 2.51882 2351.63 2366.04 300 29 2365.32 2.19857 2355.92 2366.3 301 23 2365.72 0.767842 2363.02 2366.31 302 28 2364.91 4.65057 2340.23 2366.46 303 29 2365.45 2.68121 2354.91 2367.37 304 27 2365.67 2.50324 2356.48 2367.44 305 20 2366.58 1.33849 2361.64 2367.72 306 27 2366.52 2.08566 2358.34 2367.73 307 26 2366.77 1.7302 2360.36 2367.84 308 27 2367.33 0.777778 2363.95 2367.99 309 22 2367.39 1.54601 2359.7 2368.8 310 23 2367.29 2.73863 2353.34 2368.8 311 22 2367.59 2.55826 2355.71 2368.81 312 24 2367.53 3.87914 2347.26 2368.9 313 22 2368.3 1.32131 2363.07 2368.9 314 23 2368.28 1.50034 2360.83 2368.96 315 20 2368.46 1.09434 2363.99 2368.97 316 24 2367.9 3.87969 2347.75 2369.07 317 23 2368.44 1.67129 2360.31 2369.11 318 23 2368.98 0.325 2367.72 2369.28 319 26 2368.72 1.04081 2364.37 2369.35 320 24 2368.73 1.14218 2363.72 2369.35 321 27 2367.43 5.13636 2342.14 2369.56 322 19 2369.16 0.73377 2365.61 2369.69 323 22 2367.95 4.7256 2346.17 2369.69 324 24 2368.43 2.46984 2358.87 2369.72 325 24 2369.17 1.65834 2360.28 2369.72 326 25 2369.41 0.658804 2367.32 2369.72 327 27 2368.9 1.94107 2361.36 2370.15 328 26 2368.73 2.32529 2359.35 2370.15 329 25 2369.42 1.56773 2361.2 2370.15 330 23 2369.36 1.62707 2361.83 2370.15 331 22 2367.53 9.68226 2319.96 2370.29 332 23 2369.7 1.42649 2363.74 2370.48 333 21 2369.57 1.57945 2364.57 2370.48 334 22 2370.01 0.770532 2366.42 2370.48 335 23 2369.75 1.8044 2362.19 2370.48 336 27 2369.53 1.93696 2363.12 2370.48 337 26 2369.9 1.85263 2360.55 2371.24 338 22 2369.12 5.12624 2341.62 2371.24 339 24 2370.23 1.7348 2362.41 2371.36 340 18 2370.26 2.5595 2357.8 2371.36 341 23 2371.06 0.501816 2369.46 2371.54 342 26 2370.69 1.90719 2362.89 2371.79 343 26 2371.26 0.772797 2368.28 2372.09 344 26 2370.19 3.40561 2360.37 2372.45 345 25 2371.17 4.02829 2349.37 2372.61 346 23 2370.77 4.52826 2350.85 2372.8 347 26 2371.04 4.68541 2353.19 2372.8 348 19 2372.56 0.45573 2370.99 2372.88 349 21 2372.23 1.65813 2365.14 2372.89 350 26 2372.57 0.723529 2369.46 2373.09 351 22 2372.44 1.24482 2367.29 2373.2 352 22 2372.72 0.913307 2368.2 2373.31 353 23 2372.06 2.37549 2364.15 2373.62 354 27 2372.53 2.27523 2362.18 2373.76 355 30 2372.78 1.93963 2363.25 2373.97 356 21 2373.4 1.04308 2368.82 2374.22 357 23 2373.56 0.916165 2370.33 2374.22 358 25 2373.39 1.58204 2366.66 2374.24 359 22 2373.29 2.43058 2363.06 2374.23 360 28 2373.7 1.64857 2366.46 2374.75 361 24 2374.06 0.79299 2370.69 2374.75 362 25 2373.81 1.44037 2368.59 2374.84 363 25 2374.19 1.42067 2368.03 2374.84 364 25 2373.59 3.5744 2357.99 2374.99 365 24 2374.27 2.43869 2361.22 2375.04 366 28 2374.06 2.59681 2363.85 2375.1 367 19 2374.73 1.3483 2367.43 2375.13 368 21 2374.55 1.64652 2366.12 2375.68 369 23 2374.79 1.7283 2366.79 2375.83 370 18 2375.01 1.73109 2366.79 2376.27 371 19 2375.08 2.36017 2363.61 2376.27 372 21 2375.86 0.670142 2373.73 2376.5 373 25 2375.47 3.05026 2359.3 2376.71 374 19 2375.51 3.80844 2355.17 2376.9 375 22 2375.35 3.24968 2361.23 2376.91 376 21 2376.14 1.52832 2369.71 2377.11 377 24 2376.38 2.11474 2364.87 2377.24 378 23 2376.15 2.20204 2367.53 2377.42 379 24 2376.89 0.871082 2373.43 2377.61 380 27 2376.79 1.89767 2366.81 2377.68 381 21 2375.66 4.78973 2355.01 2377.75 382 22 2376.93 2.04089 2368.11 2377.75 383 22 2377.65 0.232459 2376.61 2377.82 384 23 2376.37 5.02832 2350.13 2377.82 385 24 2377.06 2.82064 2361.83 2377.82 386 19 2377.31 1.38416 2370.69 2378.03 387 21 2377.38 1.66381 2368.85 2378.44 388 20 2376.78 3.3902 2364.06 2378.75 389 26 2377.14 2.48883 2368.03 2378.75 390 21 2378.2 0.866992 2374.78 2379.01 391 23 2377.76 2.83913 2366.57 2379.23 392 23 2378.49 1.11568 2374.13 2379.73 393 28 2378.39 1.78229 2372.09 2379.86 394 26 2378.98 1.52716 2374.21 2380.2 395 26 2377.66 7.84818 2337.54 2380.65 396 19 2378.63 6.35402 2345.78 2380.65 397 29 2379.94 1.22154 2375.7 2381.37 398 20 2380.06 1.39126 2375.87 2381.53 399 23 2380.02 2.39516 2370.28 2381.62 400 23 2380.96 1.06554 2377.08 2381.92 401 29 2381.06 1.28531 2377.67 2381.96 402 20 2380.93 2.46289 2372.05 2382.22 403 24 2381.75 0.759266 2378.91 2382.44 404 27 2380.95 4.65119 2356.26 2382.43 405 25 2382.04 0.691915 2379.56 2382.6 406 25 2381.49 3.96093 2360.19 2382.67 407 28 2380.68 4.30432 2367.63 2383.16 408 23 2381.88 3.01103 2366.21 2383.16 409 26 2382.13 1.74648 2376.58 2383.44 410 20 2382.69 1.71611 2376.06 2383.78 411 26 2382.91 1.80825 2375.66 2384.49 412 18 2383.61 0.80759 2381.38 2385.31 413 24 2383.47 2.32518 2376.48 2385.38 414 21 2383.83 3.70525 2364.5 2385.89 415 22 2384.74 1.34764 2379.15 2385.89 416 23 2385.22 1.06984 2380.37 2386.22 417 23 2385.1 1.93081 2378.89 2386.76 418 24 2383.69 8.5751 2338.89 2386.76 419 22 2385.55 1.76607 2379.04 2386.76 420 20 2385.72 3.30549 2367.82 2386.92 421 19 2386.63 0.574262 2383.77 2387 422 25 2386.38 1.29744 2381.54 2387.26 423 25 2386.97 0.325994 2385.54 2387.66 424 24 2386.93 1.01729 2383.59 2388.03 425 22 2387.14 1.45773 2381.71 2388.24 426 26 2370.94 90.0697 1877.66 2388.52 427 24 2387.61 1.45393 2381.19 2388.78 428 25 2372 88.4542 1887.53 2389.52 429 20 2354.05 127.307 1867.5 2389.52 430 20 2372.46 89.7863 1880.69 2389.52 431 28 2389.01 1.45371 2381.42 2389.61 432 24 2352.86 133.112 1813.85 2390 433 22 2389.27 1.27361 2383.65 2390.13 434 22 2373.35 88.5541 1888.34 2390.89 435 24 2353.86 137.269 1821.04 2390.9 436 28 2389.93 1.66872 2382.42 2391.1 437 27 2354.5 132.559 1848.34 2391.47 438 23 2357.08 128.511 1865.94 2391.71 439 21 2340.74 152.717 1868.96 2391.71 440 23 2357.22 127.971 1859.06 2391.81 441 26 2357.93 125.25 1876.92 2391.84 442 24 2306.66 193.275 1830.44 2391.89 443 24 2358.73 124.169 1883.61 2392.25 444 23 2375.58 88.8156 1889.13 2392.25 445 20 2325.29 172.169 1871.77 2392.46 446 28 2290.82 206.523 1827.42 2392.46 447 23 2375.15 93.7702 1861.55 2392.47 448 27 2289.58 209.262 1842.23 2392.82 449 23 2309.7 187.458 1879.73 2392.82 450 22 2325.28 171.462 1875.77 2393.23 451 26 2358.27 130.613 1852 2393.23 452 27 2324.98 175.314 1814.08 2393.24 453 27 2274.8 218.831 1835.16 2393.26 454 20 2292.12 205.39 1856.22 2393.67 455 24 2340.82 159.346 1818.92 2393.76 456 23 2292.21 206.391 1852.53 2393.76 457 27 2309.92 190.204 1867.51 2394.31 458 20 2326.9 172.544 1862.12 2394.81 459 23 2324.78 178.325 1827.5 2394.81 460 21 2377.36 92.2337 1872.21 2394.85 461 20 2325.17 178.608 1830.12 2394.85 462 23 2360.23 130.317 1843.01 2395.3 463 25 2309.71 193.515 1842.31 2396.44 464 22 2327.03 175.566 1849.5 2396.44 465 20 2326.79 176.24 1841.81 2396.44 466 19 2377.57 95.2159 1857.85 2396.5 467 25 2328.19 173.095 1864.24 2396.5 468 21 2380.06 89.6814 1888.86 2397.29 469 22 2328.63 175.755 1840.82 2397.34 470 18 2330.23 172.316 1870.46 2397.43 471 25 2293.54 209.788 1801.55 2397.48 472 23 2364.61 124.146 1891.79 2397.72 473 24 2294.55 205.868 1838.47 2397.85 474 25 2348.01 151.048 1880.41 2398.02 475 23 2232.19 239.819 1852.95 2398.46 476 26 2315.49 187.826 1872.12 2398.52 477 20 2347.77 152.997 1859.91 2399.03 478 20 2280.36 217.778 1833.33 2399.03 479 20 2298.96 202.283 1874.56 2399.03 480 20 2315.57 190.403 1854.73 2399.67 481 16 2346.61 158.896 1836.62 2399.67 482 19 2349.89 150.197 1879.73 2399.67 483 20 2282.37 216.453 1850.21 2399.67 484 19 2330.28 179.865 1828.78 2399.67 485 21 2314.82 191.408 1862.46 2399.67 486 20 2281.75 217.632 1841.98 2399.75 487 28 2265.89 224.043 1855.57 2399.78 488 26 2281.82 217.176 1843.8 2399.75 489 27 2348.93 148.672 1894.55 2399.75 490 24 2313.87 190.806 1859.24 2399.9 491 23 2332.2 175.222 1851.72 2400.13 492 21 2245.52 238.245 1829.3 2400.13 493 26 2268.24 222.281 1876.91 2400.13 494 22 2348.89 156.293 1821.2 2400.13 495 23 2349.4 153.328 1873.29 2400.13 496 27 2243.98 244.633 1796.99 2400.13 497 24 2348.42 156.404 1863.81 2400.17 498 30 2331.33 174.454 1869.62 2400.17 499 18 2283.47 215.642 1868.71 2400.17 500 19 2350.13 152.655 1878.71 2400.17 501 16 2350.04 152.307 1876.56 2400.17 502 22 2267.07 225.015 1868.13 2400.17 503 22 2347.47 159.837 1836.59 2400.17 504 23 2300.39 203.01 1863.07 2400.17 505 24 2282.01 217.462 1827.49 2400.17 506 28 2283.48 214.938 1871.06 2400.19 507 25 2351.25 148.576 1892.78 2400.19 508 24 2315.24 192.927 1828.45 2400.19 509 20 2349.41 154.623 1864.37 2400.19 510 22 2332.78 172.719 1877.48 2400.19 511 25 2301.07 201.698 1876.35 2400.19 512 25 2333.8 170.796 1885.27 2400.2 513 23 2284.79 212.831 1878.52 2400.2 514 17 2235.27 238.936 1874.99 2400.2 515 26 2283.05 215.987 1834.82 2400.2 516 25 2298.76 203.133 1869.79 2400.2 517 23 2363.66 137.838 1824.81 2400.2 518 28 2316.93 189.845 1860.73 2400.2 519 27 2382.27 92.2912 1876.99 2400.2 520 20 2317.43 188.222 1869.74 2400.2 521 27 2383.35 88.9428 1896.26 2400.2 522 25 2315.78 191.894 1839.6 2400.2 523 20 2281.35 218.29 1849.64 2400.2 524 26 2250.94 233.083 1868.93 2400.2 525 21 2299.96 203.318 1872.82 2400.2 526 27 2283.36 216.406 1861.46 2400.2 527 24 2267.45 224.197 1848.74 2400.2 528 25 2316.02 190.809 1850.77 2400.2 529 24 2316.02 191.8 1860.65 2400.2 530 25 2332.43 176.016 1822.68 2400.2 531 20 2316.05 190.294 1843.61 2400.24 532 21 2333.23 173.762 1858.07 2400.24 533 21 2333.07 171.214 1879.74 2400.32 534 23 2250.34 234.39 1856.91 2400.41 535 25 2316.28 191.187 1867.35 2400.41 536 22 2333.1 173.617 1871.47 2400.41 537 20 2316.23 190.495 1850.52 2400.41 538 26 2346.8 157.513 1833.44 2400.43 539 27 2298 206.307 1838.43 2400.43 540 19 2382.79 93.8737 1868.69 2400.43 541 17 2296.72 211.643 1826.74 2400.58 542 24 2299.16 206.84 1847.14 2400.98 543 26 2299.32 205.833 1840.67 2400.98 544 18 2351.14 150.366 1885.77 2401.07 545 24 2301.28 201.157 1876.31 2401.07 546 26 2318.1 186.992 1870.5 2401.07 547 17 2351.7 150.371 1888.77 2401.07 548 17 2298.84 208.568 1818.89 2401.07 549 26 2351.58 149.918 1883.46 2401.12 550 24 2332.5 177.185 1827.85 2401.13 551 21 2383.51 95.3786 1861.11 2401.2 552 23 2335.3 170.431 1878.31 2401.2 553 25 2215.04 250.737 1826.78 2401.2 554 24 2367.67 125.599 1881.5 2401.2 555 20 2349.95 155.767 1863.16 2401.2 556 27 2318.48 188.169 1871.64 2401.22 557 26 2317.66 188.705 1854.43 2401.73 558 21 2299.28 208.205 1830.3 2401.73 559 24 2368.46 124.353 1893.54 2401.73 560 25 2351.82 151.59 1869.43 2401.85 561 19 2334.84 173.57 1861.47 2401.95 562 26 2334.15 174.193 1838.18 2402.35 563 21 2268.53 224.854 1855.11 2402.25 564 24 2334.8 170.002 1881.4 2402.25 565 25 2369.08 124.71 1894.21 2402.25 566 26 2285.72 214.265 1862.07 2402.6 567 23 2367.35 130.671 1855 2402.47 568 26 2264.4 232.929 1841.96 2402.26 569 21 2335.89 171.292 1861.69 2402.38 570 28 2315.25 198.005 1853.97 2402.38 571 19 2334.52 175.355 1850.37 2402.69 572 27 2333.21 176.872 1860.51 2402.69 573 26 2283.54 219.962 1857.85 2402.87 574 18 2367.79 128.747 1872.41 2402.87 575 28 2385.74 92.7628 1877.66 2403.04 576 26 2336.65 170.278 1872.6 2403.04 577 23 2302.63 202.291 1863.3 2403.18 578 26 2321.14 185.597 1890.03 2403.18 579 25 2354.08 147.986 1901.56 2403.29 580 17 2353.14 151.426 1877.57 2403.56 581 27 2335.3 176.164 1845.67 2403.56 582 26 2337.22 171.178 1886.7 2404.4 583 20 2302.24 207.071 1829.76 2404.86 584 24 2317.41 194.7 1861.84 2404.86 585 24 2336.94 172.473 1877.78 2404.93 586 21 2352.19 155.293 1858.97 2404.86 587 17 2387.77 90.6382 1891.33 2404.86 588 18 2322.26 186.689 1885.22 2404.89 589 25 2334.91 180.283 1834.92 2404.89 590 26 2270.03 227.688 1842.29 2404.89 591 25 2337.84 174.022 1854.15 2404.95 592 22 2337.49 174.479 1871.63 2404.96 593 24 2335.11 179.622 1853.89 2404.96 594 25 2336.72 176.933 1845.5 2404.97 595 24 2371.68 125.257 1888.59 2404.97 596 26 2304.69 204.169 1855.98 2404.97 597 22 2338.12 173.517 1878.54 2404.97 598 22 2355.88 149.831 1892.32 2404.97 599 16 2306.1 201.124 1886.41 2405.38 600 22 2354.84 150.191 1890.76 2405.38 601 23 2339.89 168.796 1898.87 2405.38 602 18 2355.81 149.808 1890.81 2405.38 603 24 2370.92 130.101 1850.71 2405.43 604 25 2388.96 88.5125 1904.18 2405.43 605 25 2387.6 90.6746 1891.81 2405.45 606 20 2354.43 151.717 1887.09 2405.51 607 27 2320.37 190.841 1850.15 2405.58 608 24 2322.07 189.534 1877.69 2405.84 609 22 2172.08 257.244 1873.84 2405.84 610 26 2305.61 201.55 1879.27 2405.92 611 22 2354.79 152.042 1879.84 2405.92 612 25 2339.26 172.323 1878.9 2405.97 613 22 2323.33 188.222 1887.98 2405.98 614 20 2354.73 151.495 1887.98 2406.13 615 16 2373.12 124.645 1894.98 2406.13 616 26 2304.88 203.888 1864.84 2406.13 617 22 2356.24 150.423 1892.35 2406.16 618 21 2339.8 170.506 1888.31 2406.16 619 24 2354.79 156.861 1849.37 2406.17 620 15 2355.73 153.285 1886.04 2406.17 621 27 2274.59 222.98 1891.12 2406.17 622 18 2354.62 153.312 1859.12 2406.18 623 21 2355.07 154.374 1878.05 2406.18 624 25 2340.73 169.625 1893.75 2406.19 625 15 2370.83 132.86 1861.74 2406.19 626 24 2356.54 150.781 1886.27 2406.19 627 26 2303.1 209.449 1841.82 2406.46 628 19 2323.59 187.514 1892.63 2406.47 629 22 2324.24 186.485 1885.12 2406.47 630 23 2307.91 201.152 1884.48 2406.47 631 23 2355.1 155.647 1836.03 2406.47 632 16 2340.08 172.463 1866.85 2406.47 633 23 2389.91 88.6538 1904.42 2406.47 634 25 2405.91 1.15662 2401.25 2406.47 635 20 2340.17 171.879 1867.94 2406.47 636 27 2270.61 229.396 1823.5 2406.47 637 25 2306.44 203.636 1835.58 2406.47 638 25 2221.61 249.207 1850.95 2406.47 639 20 2355.8 152.924 1877.76 2406.47 640 21 2354.95 153.382 1870.12 2406.73 641 25 2338.63 175.03 1867.78 2407.16 642 19 2357.77 149.213 1895.22 2407.16 643 26 2338.39 175.954 1814.6 2407.39 644 22 2341.06 171.035 1881.49 2407.39 645 21 2273.89 225.787 1876.92 2407.39 646 25 2322.75 192.166 1821.15 2407.39 647 24 2324.9 187.84 1872.93 2407.63 648 19 2274.84 224.571 1874.9 2407.63 649 23 2324.91 187.301 1878.85 2407.63 650 26 2307.39 201.977 1866.65 2407.63 651 27 2390.32 90.5919 1894.4 2407.63 652 24 2357.32 153.222 1886.93 2407.67 653 24 2291.43 215.085 1878.04 2407.67 654 24 2372.62 132.027 1865.08 2407.67 655 21 2374.48 123.511 1902.46 2407.92 656 22 2273.74 226.901 1868.55 2407.92 657 26 2340.65 170.454 1884.48 2407.92 658 25 2357.05 151.574 1879.47 2407.92 659 24 2407.43 2.43342 2394.14 2408.14 660 20 2357.86 151.788 1878.01 2408.14 661 18 2372.02 135.47 1855.56 2408.18 662 25 2256.97 235.325 1871.63 2408.19 663 24 2308.06 203.247 1865.57 2408.19 664 29 2225.94 245.573 1867.12 2408.19 665 26 2259.23 231.871 1873 2408.19 666 21 2391.59 89.3743 1902.09 2408.46 667 25 2375.21 123.387 1903.33 2408.46 668 28 2374.8 124.776 1896.22 2408.61 669 24 2390.73 92.1097 1886.4 2408.61 670 20 2325.34 189.228 1875.74 2408.61 671 24 2340.84 171.893 1877.5 2408.62 672 22 2274.77 226.663 1868.23 2408.62 673 14 2374.99 126.91 1883.31 2408.64 674 25 2275.13 226.19 1865.81 2408.64 675 28 2241.21 241.476 1848.93 2408.64 676 24 2242.92 239.763 1854.53 2408.64 677 28 2375.08 125.923 1886.01 2408.78 678 21 2340.86 174.63 1856.05 2408.78 679 26 2325.7 188.7 1861.91 2408.78 680 23 2307.62 206.658 1828.58 2408.78 681 24 2310.43 200.009 1897.49 2408.78 682 25 2340.84 175.54 1877.27 2408.78 683 21 2341.21 175.03 1858.99 2408.78 684 22 2309.35 201.44 1873.94 2408.78 685 27 2322.7 195.393 1840.23 2408.78 686 26 2276.2 223.623 1883.06 2408.78 687 22 2341.83 171.566 1885.13 2408.78 688 20 2340.83 175.086 1876.73 2408.78 689 23 2324.55 188.959 1870.79 2408.78 690 22 2392.01 88.9493 1904.86 2408.78 691 21 2358.23 153.751 1865.7 2408.78 692 25 2342.45 172.234 1872.68 2408.84 693 25 2291.3 216.621 1870.04 2408.84 694 23 2342.03 172.675 1879.57 2408.88 695 21 2326.56 187.011 1887.74 2408.9 696 23 2308.39 204.859 1871.42 2408.9 697 25 2323.79 190.338 1873.79 2408.9 698 21 2323.31 194.977 1829.1 2408.92 699 20 2358.46 153.271 1864.66 2408.92 700 25 2340.6 175.656 1871.43 2409 Control points: [[ 29.98404996 91.88591948] [226.49665573 154.72143214] [153.818419 105.18835249] [216.16069477 136.30931099]]
In [10]:
Copied!
fitted_curve_points, t_values = generate_curve_within_boundaries(refined_control_points, image.shape, 0.01)
plt.figure(figsize=(10, 5))
plt.subplot(1, 2, 1)
plt.imshow(image, cmap='gray', origin='lower')
plt.plot(initial_fitted_curve_points[:, 0], initial_fitted_curve_points[:, 1], 'r-')
plt.plot(initial_control_points[:, 0], initial_control_points[:, 1], 'g.')
plt.title('Initial Fitted Curve on Image')
plt.axis('off')
plt.subplot(1, 2, 2)
plt.imshow(image, cmap='gray', origin='lower')
plt.plot(fitted_curve_points[:, 0], fitted_curve_points[:, 1], 'r-')
plt.plot(refined_control_points[:, 0], refined_control_points[:, 1], 'g.')
plt.title('Refined Fitted Curve on Image')
plt.axis('off')
plt.tight_layout()
plt.show()
fitted_curve_points, t_values = generate_curve_within_boundaries(refined_control_points, image.shape, 0.01)
plt.figure(figsize=(10, 5))
plt.subplot(1, 2, 1)
plt.imshow(image, cmap='gray', origin='lower')
plt.plot(initial_fitted_curve_points[:, 0], initial_fitted_curve_points[:, 1], 'r-')
plt.plot(initial_control_points[:, 0], initial_control_points[:, 1], 'g.')
plt.title('Initial Fitted Curve on Image')
plt.axis('off')
plt.subplot(1, 2, 2)
plt.imshow(image, cmap='gray', origin='lower')
plt.plot(fitted_curve_points[:, 0], fitted_curve_points[:, 1], 'r-')
plt.plot(refined_control_points[:, 0], refined_control_points[:, 1], 'g.')
plt.title('Refined Fitted Curve on Image')
plt.axis('off')
plt.tight_layout()
plt.show()
Membrane Average¶
In [11]:
Copied!
membrane_averaging = MemAverage(image=image, control_points=refined_control_points, pixel_size=2.170)
membrane_mask, averaged_2d = membrane_averaging.mem_average()
membrane_averaging = MemAverage(image=image, control_points=refined_control_points, pixel_size=2.170)
membrane_mask, averaged_2d = membrane_averaging.mem_average()
In [12]:
Copied!
fitted_curve_points, t_values = generate_curve_within_boundaries(refined_control_points, image.shape, 0.01)
plt.figure(figsize=(10, 5))
plt.subplot(1, 3, 1)
plt.imshow(image, cmap='gray', origin='lower')
plt.plot(fitted_curve_points[:, 0], fitted_curve_points[:, 1], 'r-')
plt.plot(refined_control_points[:, 0], refined_control_points[:, 1], 'g.')
plt.title('Refined Fitted Curve on Image')
plt.axis('off')
plt.subplot(1, 3, 2)
plt.imshow(membrane_mask, cmap='gray', origin='lower')
plt.title('Membrane Mask')
plt.axis('off')
plt.subplot(1, 3, 3)
plt.imshow(averaged_2d*membrane_mask, cmap='gray', origin='lower')
plt.title('Averaged 2D')
plt.axis('off')
plt.tight_layout()
plt.show()
fitted_curve_points, t_values = generate_curve_within_boundaries(refined_control_points, image.shape, 0.01)
plt.figure(figsize=(10, 5))
plt.subplot(1, 3, 1)
plt.imshow(image, cmap='gray', origin='lower')
plt.plot(fitted_curve_points[:, 0], fitted_curve_points[:, 1], 'r-')
plt.plot(refined_control_points[:, 0], refined_control_points[:, 1], 'g.')
plt.title('Refined Fitted Curve on Image')
plt.axis('off')
plt.subplot(1, 3, 2)
plt.imshow(membrane_mask, cmap='gray', origin='lower')
plt.title('Membrane Mask')
plt.axis('off')
plt.subplot(1, 3, 3)
plt.imshow(averaged_2d*membrane_mask, cmap='gray', origin='lower')
plt.title('Averaged 2D')
plt.axis('off')
plt.tight_layout()
plt.show()
