« first day (2134 days earlier)      last day (1111 days later) » 

7:24 AM
@AndrasDeak someone just sent me this in response: extremelearning.com.au/…
3
(has a lot of pretty pictures)
 
nice!
Wonder why it's not an actual publication...
 
Why go through all the trouble if you can just make a blog post, and possibly reach more people? Since he's probably not in academia anymore he doesn't have to count the publications...
 
Because if he wants it to be included in things like scipy it helps a lot if there's a peer-reviewed publication with citations
 
8:17 AM
@flawr wow that is a read
havent finished it, but added to the list
 
@AndrasDeak ok that's a good point
but then again an blog post like this might undergo a lot more scrutiny than it would in any peer review:)
 
It might, or it might not
 
in this case it seems it did judging by the comments
 
A dozen comments saying it's shit could be deleted. Ten citing papers saying it's shit can not.
 
sure but I like ranting about the current system
I like the idea behind this one: en.wikiversity.org/wiki/WikiJournal_of_Science
 
8:32 AM
I ran into joss.theoj.org a few days ago
 
Ah yes I've heard good things about this one!
 
There is also SoftwareX, but Elsevier
 
ewwww
 
hehe expected reaction, and fair one
 
9:22 AM
@AndrasDeak what?!? You can get fined for falling off your bike? Was the charge something like unlawfully holding up other traffic?
 
 
2 hours later…
11:16 AM
posted on May 04, 2021 by Steve Eddins

I have published more than 560 blog posts here since 2006, and I estimate that about 98% of them started out as MATLAB scripts.Recently, I've started writing my blog posts as live scripts. Live... read more >>

 
Sam
11:48 AM
@flawr So pretty!
 
@Adriaan no, it was something like "not considering the quality of the pavement and not respecting the weather conditions". The kid wasn't sure of the official wording but it was nonsense.
 
12:39 PM
@AnderBiguri If we want to make a picture of a single slice using ct (one source), how does the resolution of the sensor influence the resolution of the reconstruction?
I guess more = better, but how exactly? I had the impression that the resolution of the reconstruction can somehow be chosen arbitrarily but there obviously must be some trade-off
 
generally speaking, for fan-beam type of CT (i.e. point source) your image resolution is, geometrically speaking detector_pixel_mm * distance_source_object/distance_source_detector
this last ratio is often logically named "magnification"
 
detector_pixel_mm is the spacing of the detector "pixels" right?
 
now, because there is a PFS to the source, the truth is that your real resolution is more determined by the source PFS than anything else. Then you just make sure that your image pixel size is at least 3 times smaller than more or less the FWHM of the gaussian (its poisson, but whatever) of the spot size
@flawr correct. So essentially you make your image pixels project into the detector 1-to-1
anything smaller you are just wasting computational power really, as you are not measuring it
 
I see
and I guess ideally the sensor has the same distance to the source everywhere?
 
going back to the PFS, this is generally something you don't worry (unless you have a very flexible machine) because detectors are generally designed such that the pixels are 3~5 the PFS.
@flawr mm yes, and this is true for many fan-beam medical CT machines, but in cone beam, detectors are flat
this is corrected algorithmically
and with flat-field corrections
 
12:48 PM
ah I just thought it hsould be possible to calibrate that away
you know I thought I'd build my own CT
 
yeah, flat field corrections do account for most of that.
 
its all good fun until you try to acquire materials to build a radioactive source, and the goverment knocks on your door
 
hehe
I guess I have to find a really strong banana
I thought it would be a nice tool to peek into those mechanical puzzles
 
hahaha dang
tha fact that you can do this
 
12:53 PM
but I guess 12x12 pixels is not gonna cut it
 
there is a group at MIT I think that are tuning TIGRE to do "reverse tomography"
they put a photoreactive liquid and with a projector, they project X-ray measurements
so it hardens in the shape on the 3D object
its crazy XD
 
oh that sounds fun!
 
the idea is so wild, yet so ovbious when you are told about it
 
obvious to someone with a phd in ct reconstruction stuff
 
@flawr i guess you are right :D
 
12:56 PM
Isn't this similar to like the planning of radiotherapy stuff?
 
I guess it is in some sense, as you need to modulate the radiating beam to account for the body
its just the radiotherapy "end shape" its generally a blob
 
ok yes, probably not quite as sophisticated
 
Sam
@AnderBiguri btw - I did apply for a phd
 
neat! Any response?
 
Sam
Not until June
(deadline)
 
1:00 PM
Many people apply to these things so dont put all your baskets in one egg
or something like that :D
What is it about?
 
sounds very fancy!
 
neat!
solid PhD!
good luck 😀 let us know how it goes!
 
Sam
Thanks, I'll keep you all posted
 
1:56 PM
@flawr But then why not use a grid? It has been shown that “uniform random sampling” (sampling with a grid placed randomly in the domain) is a lot more efficient than random sampling (efficient == you need fewer samples to get the same amount of information out).
Yes, I have only skimmed over that post, will read it in detail later, but so far it looks like a grid with the samples ordered not in the easy sequence. I’m wondering what the benefit of that is.
 
because if you do a grid then its not really Monte Carlo anymore
I think its a way of better generating "random" samples for Monte Carlo experiments
so at samples->Inf the reandom and the ones that flwars shows are the same, thus MC techniques are still valid, but at lower number of samples, those are still random enough, yet sample better the ailable space
 
But in fact you can produce better results with Mone Carlo simulations if you generate the samples in a grid rather that randomly.
 
that is what I understood, at least
 
As long as the grid is placed randomly (I.e. the origin and rotation is randomly chosen).
 
hum I see
Don't really know enough
 
2:04 PM
I imagine scrambling the order and allowing for infinitely many samples lets you stop the simulation at any time, or continue going as long as you want. With a grid you have to decide up front how many samples you want to draw.
...makes it harder to do p-hacking if you have to decide up front, an obvious problem.
 
yeah that makes sense. With these pseudo-random generators you can run it until some arbitrary, perhaps unknown, criteria, and stop, and the results would be valid, while with a grid you need to run the grid
 
 
4 hours later…
6:31 PM
@AnderBiguri sounds like rapid prototyping with extra steps ;)
 
7:24 PM
@CrisLuengo I guess the property of grids is that they have a "low discrepancy", but unfortunately they are finite (if you consider a finite domain)
 
@flawr but you can make them as dense as you need. Not sure how this is an issue.
 
@CrisLuengo well for many applications you don't know in advance when to stop
I think
 
Right. That is the one thing that I can imagine as a reason for not using a grid.
But the whole concept of sampling till a condition is met sounds bad to me.
I can prove that a coin is biased by sampling throws until I have 5 more heads than tails. The outcome is guaranteed!
 
Then again if you have to rely on monte carlo I think the problem is already bad enough:)
 
I'm sure you prefer Las Vegas algorithms. Here in the world of reality we make do with what we have :P
 
7:39 PM
@AndrasDeak Is gambling really the best investment then? :P
 
8:20 PM
@AndrasDeak its just 3D printing!
 
 
1 hour later…
9:20 PM
So the Cauchy distribution has no mean and no variance, the student t with n=2 has a mean but no variance. Usually the definition of variance relies on the value of a mean, but could we generalize this definition and find a distribution that has no mean but a (generalized) variance?
 
@flawr no
Variance is exactly second moment minus firsst moment (mean) squared. This is not the definition, but a theorem.
So unless the generalized variance is not the variance, then no :P
 
V[X] = E[X^2] - E[X]^2 right?
what if both terms are diverge?
the difference could still converge no?
(assuming X is some kind of limit of a sequence of distributions in some space with some metric yadayadayada)
but yes I guess any other "generalized" notion of variance must be pretty stupid
but it still would be fun
@AndrasDeak have you ever heard about p-adic analyis?
it's quite wacky
 
@flawr that's Cauchy
@flawr no, because you can't subtract two infinities
@flawr and no
 
@AndrasDeak hence the generalized definition:P
@AndrasDeak so in real analysis you base any notion of convergence on the usual absolute value
(you actually construct the real number like that from the rational numbers)
 
I don't like where this is going
if someone starts talking about rationals in analysis, something's wrong :P
you can't fool me this easily, discrete math
 
9:36 PM
but on the rationals you can also define the p-adic absolute value: any ratio can be writen as p^n * a/ b (where a,b are coprime and not divisible by n)
 
yeah, no
 
and then abs(p^n * a/b) = p^-n
(p is a prime obviously)
so now you can construct the p-adic numbers
they are just as valid as real numbers
 
I guarantee I can't
@flawr no they aren't
 
not with that attitude
In number theory, Ostrowski's theorem, due to Alexander Ostrowski (1916), states that every non-trivial absolute value on the rational numbers Q {\displaystyle \mathbb {Q} } is equivalent to either the usual real absolute value or a p-adic absolute value. == Definitions == Raising an absolute value to a power less than 1 always results in another absolute value. Two absolute values | ⋅ | {\displaystyle |\cdot |} and...
^ nice guy
 
"In number theory..."
 
9:38 PM
one could get the impression you don't like numbers
but deep in your heart I know you love number theory
cause even when using your fancy machines you always use integers or rational numbers
 
powers of two are just barely numbers
 
imagine this, the sequence 1,2,4,8,16, etc converges to zero in the 2-adic numbers!
 
one more reason to mistrust this silliness
Lol, almost all survey suggestions have at least one downvote, except meta.stackoverflow.com/a/407375/5067311 meta.stackoverflow.com/a/407379/5067311 meta.stackoverflow.com/a/407303/5067311, two of which come from the same user
 
 
2 hours later…
11:29 PM
@flawr That sequence results in 0 when you overflow. Obviously.
 
protip: use uint8 for quicker convergence
 

« first day (2134 days earlier)      last day (1111 days later) »