python - Fitting a closed curve to a set of points -

-

i have set of points pts form loop , looks this:

enter image description here

this similar 31243002, instead of putting points in between pairs of points, fit smooth curve through points (coordinates given @ end of question), tried similar scipy documentation on interpolation:

values = pts tck = interpolate.splrep(values[:,0], values[:,1], s=1) xnew = np.arange(2,7,0.01) ynew = interpolate.splev(xnew, tck, der=0)

but error:

valueerror: error on input data

is there way find such fit?

coordinates of points:

pts = array([[ 6.55525 ,  3.05472 ],    [ 6.17284 ,  2.802609],    [ 5.53946 ,  2.649209],    [ 4.93053 ,  2.444444],    [ 4.32544 ,  2.318749],    [ 3.90982 ,  2.2875  ],    [ 3.51294 ,  2.221875],    [ 3.09107 ,  2.29375 ],    [ 2.64013 ,  2.4375  ],    [ 2.275444,  2.653124],    [ 2.137945,  3.26562 ],    [ 2.15982 ,  3.84375 ],    [ 2.20982 ,  4.31562 ],    [ 2.334704,  4.87873 ],    [ 2.314264,  5.5047  ],    [ 2.311709,  5.9135  ],    [ 2.29638 ,  6.42961 ],    [ 2.619374,  6.75021 ],    [ 3.32448 ,  6.66353 ],    [ 3.31582 ,  5.68866 ],    [ 3.35159 ,  5.17255 ],    [ 3.48482 ,  4.73125 ],    [ 3.70669 ,  4.51875 ],    [ 4.23639 ,  4.58968 ],    [ 4.39592 ,  4.94615 ],    [ 4.33527 ,  5.33862 ],    [ 3.95968 ,  5.61967 ],    [ 3.56366 ,  5.73976 ],    [ 3.78818 ,  6.55292 ],    [ 4.27712 ,  6.8283  ],    [ 4.89532 ,  6.78615 ],    [ 5.35334 ,  6.72433 ],    [ 5.71583 ,  6.54449 ],    [ 6.13452 ,  6.46019 ],    [ 6.54478 ,  6.26068 ],    [ 6.7873  ,  5.74615 ],    [ 6.64086 ,  5.25269 ],    [ 6.45649 ,  4.86206 ],    [ 6.41586 ,  4.46519 ],    [ 5.44711 ,  4.26519 ],    [ 5.04087 ,  4.10581 ],    [ 4.70013 ,  3.67405 ],    [ 4.83482 ,  3.4375  ],    [ 5.34086 ,  3.43394 ],    [ 5.76392 ,  3.55156 ],    [ 6.37056 ,  3.8778  ],    [ 6.53116 ,  3.47228 ]])

actually, not far solution in question.

using scipy.interpolate.splprep parametric b-spline interpolation simplest approach. natively supports closed curves, if provide per=1 parameter,

import numpy np scipy.interpolate import splprep, splev import matplotlib.pyplot plt  # define pts question  tck, u = splprep(pts.t, u=none, s=0.0, per=1)  u_new = np.linspace(u.min(), u.max(), 1000) x_new, y_new = splev(u_new, tck, der=0)  plt.plot(pts[:,0], pts[:,1], 'ro') plt.plot(x_new, y_new, 'b--') plt.show()

enter image description here

fundamentally, approach not different 1 in @joe kington's answer. although, bit more robust, because equivalent of i vector chosen, default, based on distances between points , not index (see splprep documentation u parameter).


Comments

Post a Comment

Popular posts from this blog

c# - Better 64-bit byte array hash -

-
i need hash algorithm produces 64-bit hash code ( long ) fewer collisions string.gethashcode() , fast (no expensive calls cryptographic functions). here's implementation of fnv still shows 3% of collisions after testing 2 million random strings. need number way lower. void main() { const string chars = "abcdefghijklmnopqrstuvwxyzabcdefghijklmnopqrstuvwxyz!#@$%^&*()_+}{\":?><,./;'[]0123456789\\"; const int n = 2000000; var random = new random(); var hashes = new hashset<long>(); int collisions = 0; for(int = 0; < n; i++) { var len = random.next(chars.length); var str = new char[len]; (int j = 0; j < len; j++) { str[j] = chars[random.next(chars.length)]; } var s = new string(str); if(!hashes.add(get64bithash( s ))) collisions++; } console.writeline("collision percentage after " + n + " random strings: " + ((doubl
Read more

webrtc - Which ICE candidate am I using and why? -

-
[questions in bold below] i have setup kurento media server 5.1.3 in datacenter behind firewall running os ubuntu 14.04. has 2 network cards: 222.222.222.222 (eth0 - private ip) 111.111.111.111 (eth1 - public ip) attached below sdp (setremotedescription) when browser connected kurento media server type: answer, sdp: v=0 o=- 5487318114793304426 0 in ip4 0.0.0.0 s=kurento media server c=in ip4 0.0.0.0 t=0 0 a=group:bundle audio video m=audio 59068 rtp/savpf 111 0 c=in ip4 111.111.111.111 a=rtpmap:111 opus/48000/2 a=rtpmap:0 pcmu/8000 a=sendrecv a=rtcp:59068 in ip4 111.111.111.111 a=rtcp-mux a=ssrc:669011897 cname:user39019747@host-6e83e4c2 a=extmap:3 http://www.webrtc.org/experiments/rtp-hdrext/abs-send-time a=mid:audio b=as:20 a=ice-ufrag:ymdk a=ice-pwd:lylifk5ueqzpwm91ddj37e a=fingerprint:sha-256 ff:0f:81:8c:41:4e:b4:b6:c6:d8:36:f3:d6:5f:09:fd:5f:af:13:b3:9d:fc:12:66:ac:f3:56:d6:5b:0a:73:5d a=candidate:1 1 udp 2013266431 111.111.111.111 55239 typ host a=candidate:2 1 udp
Read more

php - Zend Framework / Skeleton-Application / Composer install issue -

-
i trying create zend framework application using skeleton in netbeans. when run composer install, following error message: "c:\php\php.exe" "c:\composer\composer.phar" "--ansi" "--no-interaction" "update" "--dev" using deprecated option "dev". dev packages installed default now. loading composer repositories package information updating dependencies (including require-dev) requirements not resolved installable set of packages. problem 1 - package requires php >=5.5 php version (5.4.42) not satisfy requirement. problem 2 - installation request zendframework/zendframework 2.5.1 -> satisfiable zendframework/zendframework[2.5.1]. - zendframework/zendframework 2.5.1 requires php >=5.5 -> php version (5.4.42) not satisfy requirement. done. apparently, zend framework only needs php 5.3+ contradicts above. need use php 5.4. my composer.json following: { "name": "
Read more