2019.3.28 医学图像处理行业留下的点点滴滴

前两个礼拜弄出来的二维成像的处理：
（matlab）源代码3.28
clc; clear all; close all;
fidin=fopen(‘F:\上班的工具汇总\3.14 三维超声\3.14 采集水的图像数据\水经过降采样得到的数据\1.bin.txt’,‘rb’);
%% matlab读的数据转化为16位
A=fread(fidin,‘uint16’);%该语句的核心
A=A(1:232128);
B=reshape(A,[936,248]);
figure;
% imshow(B,[0,65535])
imtool(uint16(B));
imwrite(B,‘CB.jpg’)
d=B;
%% 进行坐标变换 3.25修改
%下面进行坐标变换，从直角坐标系转换到扇形坐标系下
%设新的图像的大小仍为原数据矩阵大小,则每线有m个数据，一共有n线
[m n] = size(d);
data = zeros(m,n);
data = data - 1;
%根据仪器构造，扫描角度为68°，设角度圆心到图像的起始半径为r（像素），整个扇形有效数据的半径为n+r（像素）
deltasita = 68/248/180pi;
Fsita = 68/180pi;
r = 70;
Y = rcos(34/180pi);
X = fix(n/2);
for sita = deltasita:deltasita:Fsita
for i = 1:m
R = r + i;
alpha = (270-3