-
2 5.4mm 5.27mm 1600px 5400um 0.00527m 1600px
样例输出
-
Case 1: 1639.47px Case 2: 1639.47px
描述
一般来说,我们采用针孔相机模型,也就是认为它用到的是小孔成像原理。
在相机坐标系下,一般来说,我们用到的单位长度,不是“米”这样的国际单位,而是相邻像素的长度。而焦距在相机坐标系中的大小,是在图像处理领域的一个非常重要的物理量。
假设我们已经根据相机参数,得到镜头的物理焦距大小(focal length),和相机胶片的宽度(CCD width),以及照片的横向分辨率(image width),则具体计算公式为:
Focal length in pixels = (image width in pixels) * (focal length on earth) / (CCD width on earth)
比如说对于Canon PowerShot S100, 带入公式得
Focal length in pixels = 1600 pixels * 5.4mm / 5.27mm = 1639.49 pixels
现在,请您写一段通用的程序,来求解焦距在相机坐标系中的大小。
输入
多组测试数据。首先是一个正整数T,表示测试数据的组数。
每组测试数据占一行,分别为
镜头的物理焦距大小(focal length on earth)
相机胶片的宽度(CCD width on earth)
照片的横向分辨率大小(image width in pixels),单位为px。
之间用一个空格分隔。
输出
每组数据输出一行,格式为“Case X: Ypx”。 X为测试数据的编号,从1开始;Y为焦距在相机坐标系中的大小(focallength in pixels),保留小数点后2位有效数字,四舍五入取整。
数据范围
对于小数据:focal length on earth和CCD width on earth单位都是毫米(mm)
对于大数据:长度单位还可能为米(m), 分米(dm), 厘米(cm), 毫米(mm), 微米(um),纳米(nm)
代码如下:
//source her #include<vector> #include<map> #include<iostream> #include<algorithm> #include<string> using namespace std; int main() { int T; double b; double a=0; double danwei=0; string temp,dw1; int dw=0; cin>>T; vector<double> *p=new vector<double>[T]; for(int i=0;i<T;++i) { for(int j=0;j<3;++j) { cin>>temp; //cout<<temp; dw1.clear(); dw=0; a=0; b=1; for(string::iterator k=temp.end()-1;!(*k>='0'&&*k<='9');k--) { dw+=*k; dw1+=*k; } switch(dw) { case 'm'+'m':danwei=1;break; case 'm'+'u':danwei=0.001;break; case 'm':danwei=1000;break; case 'm'+'d':danwei=100;break; case 'm'+'c':danwei=10;break; case 'm'+'n':danwei=0.000001;break; case 'p'+'x':danwei=1;break; default :break; } string::iterator pos=find(temp.begin(),temp.end(),'.'); if(pos==temp.end()) { b=1; for(string::iterator pos1=temp.end()-dw1.size();(--pos1)!=temp.begin();) { a+=(*pos1-'0')*b; b*=10; } a+=(*temp.begin()-'0')*b; } else { for(string::iterator pos1=pos;(--pos1)!=temp.begin();) { a+=(*pos1-'0')*b; b*=10; } a+=(*temp.begin()-'0')*b; b=0.1; for(string::iterator pos2=pos;(++pos2)!=temp.end()-dw1.size();) { a+=((double)(*pos2-'0'))*b; b*=0.1; } } a*=danwei; p[i].push_back(a); } double result1=(p[i][2]*p[i][0])/p[i][1]; int result2=result1*100; if(((int)(result1*1000)%10)>=5) ++result2; result1=(result2/100.0); cout<<"Case "<<i+1<<": "<<result1<<"px"<<endl; } //system("pause"); return 0; }