cedtest.asv

超声图像实时消噪实验程序

    a = dicomread( 'F:\photos\contextvision\dicom\1.0011\1.2.756.9999.999.99.9.0011.0001\1.2.840.10008.5.1.4.1.1.7.1.1.2.0011.001.dcm' );
    figure(5);
    imshow(a,[]);
%     a = im2double( a );
%     ad = a( 51:125, 76:150 );
    ad = double(a);
    %ad = log(ad+1);
    bd = ced(ad, 0.4, 1, 1, 0.15, 4, 1);
    %bd = exp(bd);
    figure(6);
    imshow(bd,[]);
    bdmax = max(max(bd));
    bdmin = min(min(bd));
%     bd = uint16(bd);
%     fid = fopen('d:\chuliqian.raw','wb');
%     count = fwrite(fid,bd,'uint16');
%     fclose(fid);

% write to bmp file
    bd = uint8(bd*255/(bdmax-bdmin));
    imwrite(bd,'e:\us\0.4, 1, 1, 0.1541.bmp');

% h = fspecial('average');
% h = fspecial('gaussian',4,0.3);
% re = filter2(h,ad,'same');
% figure(7);
% imshow(re,[]);
% title('average Image');