有的时候,在生物学分析时,需要的大脑的分辨率可能不需要太高,只需要固定在某一分辨率等级下能够看清楚就行,因此不需要terafly模式,从而减少观察的等待时间。
在高分辨率下的手动绘制的swc文件,可能需要进行down_sample
%para
V3DIOPATH='D:\v3d_2013\v3d_external\matlab_io_basicdatatype';
%输入文件
inputfile='17302_053_stamp_2019_03_14_16_24_222.swc_resampled.swc';
output_suufix='_2.eswc';
%输入分辨率
res_input_x=54600;
res_input_y=34412;
res_input_z=9847;
%输出分辨率
res_output_x=1706;
res_output_y=1075;
res_output_z=307;
%程序开始
if ~exist(V3DIOPATH,'dir')
%disp(['error , empty V3DIOPATH in ' V3DIOPATH])
error(['error , empty V3DIOPATH in ' V3DIOPATH]);
else
addpath(V3DIOPATH)
end
%
inputswc=load_v3d_swc_file(inputfile);
outswc=inputswc;
outswc(:,3)=inputswc(:,3)/res_input_x*res_output_x;
outswc(:,4)=inputswc(:,4)/res_input_y*res_output_y;
outswc(:,5)=inputswc(:,5)/res_input_z*res_output_z;
outputfile=[inputfile output_suufix];
save_v3d_swc_file(outswc,outputfile);
disp(['save on ' outputfile])
实际上只需要坐标对应的除以比例就行了。