function?fx=PV_MI（AABBx）;

[rowAAcolumnAAhighAA]=size（AA）;
[rowBBcolumnBBhighBB]=size（BB）;
AA=AA+1;
BB=BB+1;%全部加1，避免為0時，作為下標時出錯

Numab=zeros（256256）;
Numa=zeros（1256）;
Numb=zeros（1256）;
count=0;
for?vi=1:highBB
????????????for?vj=1:columnBB
????????????????for?vk=1:rowBB
????????????????????vkm=vk*2.590723;
????????????????????vjm=vj*2.590723;
????????????????????vim=vi*8.000000;
????????????????[xxmyymzzm]=transformation（vkmvjmvimx）;
????????????????xx=xxm/1.250000;
????????????????yy=yym/1.250000;
????????????????zz=zzm/4.000000;
???????????????if?xx>1&xx1&yy1&vk???????????????????count=count+1;
???????????????end
????????????????
?????????????
????????????????grayvaluex=BB（vkvjvi）;
????????????????if?xx<1|xx>rowAA|yy<1|yy>columnAA|zz<1|zz>highAA
????????%if?ptA（vivj）<1&ptB（vivj）<1??????%出界點處理
????????????%Numab（FBB（11））=Numab（FBB（11））+1;
????????%elseif?ptA（vivj）<1&ptB（vivj）>columnBB
????????????%Numab（FBB（1columnBB））=Numab（FBB（1columnBB））+1;
????????????%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
????????????%出界處理一
???????????????%?Numab（grayvaluex1）=Numab（grayvaluex1）+1;
????????????????%Numa（grayvaluex）=Numa（grayvaluex）+1;
????????????????%Numb（1）=Numb（1）+1;
????????????%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
????????????%出界處理二
????????????%1xx?xia?chu?jie
????????????????if?xx<1&yy>=1&yy<=columnAA&zz>=1&zz<=highAA
???????????????????grayvaluey=AA（1round（yy）round（zz））;
???????????????????Numab（grayvaluexgrayvaluey）=Numab（grayvaluexgrayvaluey）+1;?
???????????????????Numa（grayvaluex）=Numa（grayvaluex）+1;
???????????????????Numb（grayvaluey）=Numb（grayvaluey）+1;
????????????????end
????????????%2xx?shang?chu?jie
????????????????if?xx>rowAA&yy>=1&yy<=columnAA&zz>=1&zz<=highAA
???????????????????grayvaluey=AA（rowAAround（yy）round（zz））;
???????????????????Numab（grayvaluexgrayvaluey）=Numab（grayvaluexgrayvaluey）+1;?
???????????????????Numa（grayvaluex）=Numa（grayvaluex）+1;
???????????????????Numb（grayvaluey）=Numb（grayvaluey）+1;
????????????????end
???????????%3yy?xia?chu?jie
????????????????if?yy<1&xx>=1&xx<=rowAA&zz>=1&zz<=highAA
???????????????????grayvaluey=AA（round（xx）1round（zz））;
???????????????????Numab（grayvaluexgrayvaluey）=Numab（grayvaluexgrayvaluey）+1;?
???????????????????Numa（grayvaluex）=Numa（grayvaluex）+1;
???????????????????Numb（grayvaluey）=Numb（grayvaluey）+1;
????????????????end
????????????%4yy?shang?chu?jie
????????????????if?yy>columnAA&xx>=1&xx<=rowAA&zz>=1&zz<=highAA
???????????????????grayvaluey=AA（round（xx）columnAAround（zz））;
???????????????????Numab（grayvaluexgrayvaluey）=Numab（grayvaluexgrayvaluey）+1;?
???????????????????Numa（grayvaluex）=Numa（grayvaluex）+1;
???????????????????Numb（grayvaluey）=Numb（grayvaluey）+1;
????????????????end
????????????%5zz?xia?chu?jie
????????????????if?xx>=1&xx<=rowAA&yy>=1&yy<=columnAA&zz<1
???????????????????grayvaluey=AA（round（xx）round（yy