Single pixel event image processing in python -


i have had lonng term problem proceccing ccd image of xray source, attached hereccd image after theresholding arbitrary value, thing need subtract multi-pixel events image. need count number of pixels belongs single pixels.

*- multi pixel mean pixels has nonzero value in sorounding pixels.

i have code working using pil.image.open() read in list , analysis pixel pixel! i'm looking standard image prcessing routin more reliable , better results. appriciate if can give me how it.

cheers

you can quite imagemagick installed on linux distros , available osx , windows. there c/c++, python, perl, php, ,net , other bindings available. i'll @ command line here.

so, here start image:

enter image description here

first, let's threshold @ arbitrary 80%:

convert ccd.png -threshold 80% result.png

enter image description here

now let's find white pixels have no other white pixels around them:

convert ccd.png -threshold 80% -morphology hmt peaks:1.9 result.png

enter image description here

the technique "hit-and-miss morphology" , described here in anthony thyssen's excellent imagemagick usage pages. search on page "peak" find exact paragraph.

now let's count them:

convert ccd.png -threshold 80% -morphology hmt peaks:1.9 -format "%[fx:int(mean*w*h)]" info: 975

let's check multi-pixel events removed:

convert ccd.png -threshold 80% \( +clone -morphology hmt peaks:1.9 \) -compose difference -composite result.png

enter image description here


-

Comments

Post a Comment

Popular posts from this blog

serialization - Convert Any type in scala to Array[Byte] and back -

i have following question: i have variable value in program declared value. i want convert value byte array.. how can serialize byte array , back? found examples related other types such double or int, not any. this should need. it's pretty similar how 1 in java. import java.io.{bytearrayinputstream, bytearrayoutputstream, objectinputstream, objectoutputstream} object serialization extends app { def serialise(value: any): array[byte] = { val stream: bytearrayoutputstream = new bytearrayoutputstream() val oos = new objectoutputstream(stream) oos.writeobject(value) oos.close stream.tobytearray } def deserialise(bytes: array[byte]): = { val ois = new objectinputstream(new bytearrayinputstream(bytes)) val value = ois.readobject ois.close value } println(deserialise(serialise("my test"))) println(deserialise(serialise(list(1)))) println(deserialise(serialise(map(1 -> 2)))) println(deserialise(seri
Read more

matplotlib support failed in PyCharm on OSX -

any idea how fix ? how can debug might problem ? not seem : https://youtrack.jetbrains.com/issue/py-15520 i in console: python 2.6.9 (unknown, jul 14 2015, 19:46:31) [gcc 4.2.1 compatible apple llvm 6.0 (clang-600.0.39)] on darwin >>> import matplotlib.pyplot plt matplotlib support failed the backend : >>> import matplotlib >>> print(matplotlib.rcparams['backend']) macosx
Read more

python - Matplotlib: TypeError: 'AxesSubplot' object is not callable -

i new programming please pardon me if rookie mistake. i writing python script have defined functions create subplots of distributions. i have slider widget end user change value of variables(mu , sigma) , see effect on distribution. here's code import numpy np print np.__version__ import scipy print scipy.__version__ scipy.stats import norm, lognorm, uniform import matplotlib print matplotlib.__version__ import matplotlib.pyplot plt matplotlib.widgets import slider, button, radiobuttons, checkbuttons #####calculating mean , standard deviation##### global mu_a1 mu_a1 = 1 global mu_b1 mu_b1 = 10 global mu_c1 mu_c1 = -13 global sigma_a1 sigma_a1 = 0.14 global sigma_b1 sigma_b1 = 1.16 global sigma_c1 sigma_c1 = 2.87 mu_x01 = -11 sigma_x01 = 1.9 #####_____##### #####generating random data##### a1 = 0.75*mu_a1 + (1.25 - 0.75)*sigma_a1*np.random.sample(10000) b1 = 8*mu_b1 + (12 - 8)*sigma_b1*np.random.sample(10000) c1 = -12*mu_c1 + 2*sigma_c1*np.random.sample(10000)
Read more