python - find average shape from seed points -


i wondering if there way can take array of xy coordinates , find best fit shape those. example have circular map of angles ranging -90 90 in disk-like shape negative angles ring on outside, 0 in middle , positive toward inside. however, these not uniform or circular, , want able fit similar values shape. far have bit of code pluck out values close zero: 

slicedim = 7 x = r_[61] y = r_[55] in range (61,140): j in range (55,131):      if fa[i,j,slicedim] < 0.2:         continue     if ha_array[i,j,slicedim] < 10 , ha_array[i,j,slicedim] > -10:         i_array = np.array([i])         j_array = np.array([j])         x = concatenate((x,i_array))         y = concatenate((y,j_array))

so have array of x , y values correlate zero. , far have been able use these values give best fit circle. , generate image:

fit circle

but can see more of oval shape. there make shape more accurately reflects data?


Comments

Post a Comment

Popular posts from this blog

Linux vanilla kernel on QEMU and networking with eth0 -

i have downloaded , compiled vanilla linux kernel (3.7.1) used busybox ramdisk booted using qemu. qemu command line qemu-system-i386 -kernel bzimage -initrd ramdisk.img -append "root=/dev/ram rw console=ttys0 rdinit=/bin/ash" -nographic -net nic -net user everything goes well. however, can't use networking on vanilla kernel busybox. 'ifup eth0' tells me / # ifup eth0 ip: siocgifflags: no such device i googled internet can' clue... advice nice thank in advance. most there no driver (in example should e1000) loaded or device has name. in /sys/class/net/ should find listing of available net-devices. if there none (besides lo) driver not loaded. in qemu monitor type "info pci" , show pci-address of ethernet card. should this: ... bus 0, device 3, function 0: ethernet controller: pci device 8086:100e ... this device corresponds /sys/devices/pci0000:00/0000:00:03.0/. files "vendor" , "device" must c...
Read more

rdbms - what exactly the undo information lives in oracle? -

some docs says, there "undo tablespaces" implies should create tablespace undo imformation. some docs says, there "undo segment" implies undo information lives in normal tablespaces , use of it's segment. so, undo information lives , organized? earlier releases of oracle database used rollback segments store undo. oracle9i introduced automatic undo management, simplifies undo space management eliminating complexities associated rollback segment management. oracle recommends (oracle 9i , on words) use undo tablespace (automatic undo management) manage undo rather rollback segments. when creating undo tablespace, these automatically created: n undo segments (based on sessions parameter value) named _syssmun$ owned public (usable ops configuration) not manually manageable
Read more

clojure - 'get' replacement that throws exception on not found? -

i'd access values in maps , records, throwing exception when key isn't present. here's i've tried. there better strategy? this doesn't work because throw evaluated every time: (defn get-or-throw-1 [rec key] (get rec key (throw (exception. "no value.")))) maybe there's simple method using macro? well, isn't it; has same problem first definition, if evaluation of throw happens later: (defmacro get-or-throw-2 [rec key] `(get ~rec ~key (throw (exception. "no value.")))) this 1 works letting get return value (in theory) never generated other way: (defn get-or-throw-3 [rec key] (let [not-found-thing :i_would_never_name_some_thing_this_021138465079313 value (get rec key not-found-thing)] (if (= value not-found-thing) (throw (exception. "no value.")) value))) i don't having guess keywords or symbols never occur through other processes. (i use gensym generate special value of...
Read more