搬家啦!请去goldengrape.org找我
  • 购买了“数字游牧计划| Digital Nomads Project”的服务,建立了独立Blog,

    “金色葡萄的精华区”。地址是http://goldengrape.org

    服务器地址不在国内,教育网的读者,请继续访问在blogbus上的镜像“金色葡萄的国内精华区”,或者使用能够访问国外网站的工具。比如sogo浏览器。

    我的Feed订阅地址,仍然保持为http://feed.feedsky.com/zhen 我会把feed重定向到新的位置上。对于blogbus本身提供的那个feed,我似乎无法操纵它,请读者更换到feedsky的地址上。至于Feedburn上的地址,由于长期处于废弃状态,也请您更换服务到feedsky。谢谢

    在重定向Feed,以及对blog进行导入导出的过程中,有可能会将你以往读过的文章重新发布一次。如果您发现有大量过往文章出现在您的订阅中,请直接将它们标记为已读。造成的不便,非常抱歉。

    今后的blog文章,绝大部分会同时发布在两个站点中。此处作为镜像,我会逐渐关闭评论功能。您的评论,最好发布在http://goldengrape.org上。

    感谢数字游牧计划| Digital Nomads Project的帮助

    感谢shizhao的介绍和ooof的帮助

    最后,感谢Blogbus长期以来周到的服务。

  • 字母Jaeger视力表,详解

    注意:本文保留所有权利。基于下述实验方法发文章的话,要和我联系,根据文章的实验目的,我可能要求作为第一或者第二至第N作者。如果出于非商业和非科研的目的,您可以自由利用文中的方法和程序。

    Intro:给人查视力是个很烦的事情。一个人要花去好几分钟。不过,视力对眼科的疾病诊断确实非常重要。病人来看病,拒绝查视力,就好像看中医,不让号脉一样。具体解释,请看雪儿医生的《看病为什么要查视力》和《查视力<>验光》。还有,就是查视力是个线性过程,看完低的看高的,查完一个病人,再查另一个病人,与我急躁的高通量处理病人的原则不符。故改进之。

    不失一般性,以查近视力的Jager视力表为例。利用不同空间频率的条纹作为字母的笔画书写字母,背景色采用条纹的平均颜色,如果被试者能够分辨条纹,则能够发现背景色块中的字母,反之则不能。将不同的字母以不同空间频率的条纹来书写,随着视力的不同,能够写下的字符也不尽相同。通过读取被试者写下的字符,就能够测量被试者的视力。

    方法:

    第一步,Matlab程序,没兴趣的请直接略过:

    clear all, close all

    % 清理工作区

    w=600; % 产生一个长宽都是600个像素的图片,这次做大一点

    r_dpi=600; % 设定打印机的分辨率。

    jr=[30,45,60,90,120,180,300]; % Jaeger视力表的定义,距离多少cm处的分辨角为1′

    kjr=tan(pi/180/30)jr/2.54r_dpi; % 一个周期占多少个点,除30是因为分辨角为1′的话,一个周期包含一个白条和一个黑条,所以是2′,

    ks=w./kjr; % 在W个像素下,到底转了多少个周期

    -------------

    http://goldengrape.org/2009/04/char_jaeger_visual_acuity/

     

  • 自杀预防的资料

    2009-04-24

    Tag: 发现


    WHO Suicide prevention (SUPRE)

    http://www.who.int/mental_health/prevention/suicide/suicideprevent/en/index.html

    北京心理危机研究与干预中心:
    其中有WHO SUPRE的翻译资料:
    WHO自杀预防供通科医生参考:
    北京心理危机研究与干预中心热线:
    全国免费接听热线:            
    800-810-1117
    手机和电话卡用户拨打号码:  
    010-8295-1332
    —-
    一个就已经太多


  •  

     

    进宝,第一季第一集

    进宝,第一季第一集

     

  • 给猫测视力和对比敏感度,详解

    注意:本文保留所有权利。基于下述实验方法发文章的话,要和我联系,根据文章的实验目的,我可能要求作为第一或者第二至第N作者。如果出于非商业和非科研的目的,您可以自由利用文中的方法和程序。

    基本原理:猫是好奇的动物,本能会追逐移动的物体。利用计算机产生一个条纹图形,它的对比度和空间频率已知,然后放置在一个背景色块之上移动,背景色块的颜色是前景色的平均值。如果猫能够看到,那么就会有追逐行为发生,反之则没有。为了一开始吸引猫咪的注意力,前景条纹还要在一个背景色明显不同的地方先移动一下。

    第一步:产生前景和背景图片:

    下面是Matlab程序,用来产生前景条纹和背景色块。

    clear all; close all

    % 清理工作区

    w=200;% 产生一个长宽都是200个像素的图片

    k=72;    % 空间频率,这里是在200个像素的长度里面,有k个周期,实际的定义需要根据具体显示器的大小而定,

    n=256^2*1/3 

    % 65536是屏幕的最大亮度,也就是全白

    http://goldengrape.org/2009/04/measure_cat_visual_acuity_and_csf/

  • 花钱买健康

    2009-03-27

    Tag: 创意

    花钱买健康

    2009年3月23日

    保持健康是需要成本的。要享受到健康,有几种花钱的方式。

    第一,选择良好的生活环境。这基本是不可能做到的。北京,根本就不适合人类居住。人太多的地方,人口聚居的因素主要是因为发展机会多,并非由于适合居住。不过反之,论证则是成立的,像什么香格里拉之类的仙境,必然是因为不适合人类才能被神仙占据。一个中小城市才是适宜健康的地方,可惜我们都不会去的。

    其后的顺序是运动,食品,药品,理疗,手术,最不推荐的是保健品。如果把宗教信仰也纳入其中的话,应该也排在保健品之前。

    先从几个常见的逻辑谬误来说。难以治疗的疾病,和名贵稀有的药材之间没有任何关系。

    黄金、钻石稀有,吃下去是要死人的,铀、鈈也是稀有的,那是用来制造核燃料的。。

    运动:运动的目的以锻炼心肺功能为主,要做有氧运动,(像跑百米那种上气不接下气的是无氧运动)。运动之中关节最好少受压力。因此,最好的运动有三种:游泳,游泳,游泳。同时,太极拳并不是很好的健身运动,尤其是你打得好的话,磨损膝关节,还是慢慢的研磨。

    。。。。原文:花钱买健康

    金色葡萄创意



  • 2009年3月22日

    Blogger,是Blog Service Provider( BSP)的用户,也是BSP的资源。一个Blogger脱离了BSP,自己租用主机,建立独立的Blog,成为一个独立Blogger,可以视为是BSP的损失。

    据我猜测,一个blogger,要从bsp脱离,有如下特征可以识别,目前还没有想好权重的评分,以后更新。

    • 使用词汇“Blog”而不使用“博客”的人,最终一定会走向独立

    • 写作时间超过3年以上


    欢迎补充

    金色葡萄发现

  • 理想的GUI

    2009-03-19

    Tag: 创意

    原文链接:理想的GUI

    2009年3月19日

    想必你已经看过MIT的第六感,也看了微软的2019。实在都是非常的震撼。

    However,(我真的喜欢however这个词)…

    先看微软的2019,要实现这样的效果,必须四处都布满微软的显示器和触摸屏,当然这是M$希望看到的,却违反了商业常识。即使显示器发明了这么多年了,至今 也没有公用的,未来也不会有公用的显示设备,第一,没法收费,第二,它会用来播放广告。即使是最便宜的显示设备,“纸”,也没有能够在大街上能够随便拿来一张用。在完全实现公有制之前,不可能出现可以供个人使用的公共媒体。

    未来的随身设备,必然属于且仅仅属于个人。你不得不自己带着输入、输出设备。

    再看,MIT的第六感,这个有实际的演示样品,比微软那个靠谱多了,而且人家自带投影仪,看得人心潮澎湃。如果有这么一套,价钱在可接受的范围内,我一定要搞一套。不过投影的问题是隐私,虽然屏幕够大,看的够爽,不过身边的人也都知道你在看什么。尤其是那个把投影投到对方身上的身份识别程序,其实我觉得那个非常的不礼貌。如果可能的话,我更期望一个眼镜,用来显示。当然难度大一些,因为离得太近显示其实不爽,正确的做法是在眼镜上显示整个物体发光的波相位,也就是全息。

    ......原文链接:理想的GUI

    ideal GUI

    金色葡萄创意


     

  • 创意:wii绑腿

    2009-03-15

    Tag: 创意

    手工DIY。。

    用大号的可乐瓶子,切掉头部和底部,从中间剖开,边角剪成圆角,四角连接绳子,绑住大腿。剖线对侧切开4个小口,用绳子绑住wii,

    打算找个裁缝,做一个布的。更进一步,可以量产成wii运动服。

    ...

    http://goldengrape.org/2009/03/wii_leg_band/

  • 很少在blog中提起正在研究的问题。抛砖一下:

    首先,请大致了解一下image deblurring

    上一篇我贴出的程序,其实我想做这样一件事情,在一张照片上,远景和近景调换一下,图像near,是聚焦于近景,图像far,是聚焦于远景,(感谢进宝小朋友的友情出演,虽然两张照片的拍摄时间相隔几分钟,但是进宝依然镜头感很好,保持几乎不变的姿势),两者之间相差的距离,根据相机的镜头,可以换算成屈光度。于是可以计算出两者之间PSF的差距。

    那么,在不考虑噪声的条件下,如果认为far是理想图像,那么near就是离焦状态,可以基于离焦PSF矫正之。反之亦然。

    接下来,出现了个令人纠结的问题。

    可以根据任何程度的离焦图像,重建清晰的图像。即使2000度的瓶子底级别的近视,理论上,他不戴眼镜,也获得了所有图像的信息,经过适当的处理,他看到的和我一样清楚。真是不甘心啊。问题是,如何解释离焦的MTF与没有离焦的MTF之间的差别呢?更严重的是,MTF的截止频率如何解释。

    再接下来,距离和离焦是可以互相换算的。 

    ......

    http://goldengrape.org/2009/03/from_image_deblurring_to_3d/

  • 这事不能只吓唬我一个人,大家都要分享一下。以下现象在Matlab 7.0 for Mac上实现:

    clear all

    close all

    %清除工作区,关闭所有图片窗口,注意,我已经清除了工作区了!所有的变量默认值都应该是空

    near=imread(’I1.tif’);

    far=imread(’I2.tif’);

    near=rgb2gray(near);

    far=rgb2gray(far);

    imagepic=far;

    %读取两个图片文件,转换成黑白,交给一个变量imagepic,注意变量名是imagepic,不是image

    [lens_sizex,lens_sizey]=size(near);

    D=zeros(lens_sizex,lens_sizey); 

    % 建立一个尺寸与读取图形大小一样的0矩阵

    centerx=floor(lens_sizex/2)+1;

    centery=floor(lens_sizey/2)+1;

    D(centerx,centery)=1;

    % 取0矩阵的中心点=1

    R=bwdist(D);

    % 以欧氏距离返回各点到1点的距离,也就是说(x,y)点的值=半径,相当于直角坐标系转换成极坐标系。

    for DefocusInDiopters=2;

    % 这个for循环在Matlab中是允许的,相当于for i=2 to 2.

        PARAMS3=6;

        c4=(1/(4sqrt(3)))DefocusInDiopters*((PARAMS3/2000)^2);

        wavefront=c4sqrt(3)((2)R.R+(-1)*R.^0);

        % Zernike函数的转换,没啥特殊,把第四项离焦的像差计算出来

        pupilfunc=exp(-i2pi/0.550.*wavefront);

        Hamp=fft2(pupilfunc);

        Hint=(Hamp .* conj(Hamp));

        PSF = real(fftshift(Hint));

        %计算离焦产生的点扩散函数

        wnr1 = deconvwnr(image,PSF);

        %用图像处理工具箱之中的维纳滤波复原,试图使用离焦的点扩散函数复原图像

        %注意!!这里拼写错误,应该使用imagepic,而写成了image,image应该是没有定义的!!!

        figure;imshow(wnr1);

        %显示图像

    end

    显然程序报错:

    ============

    http://goldengrape.org/2009/03/ghost_in_matlab/