使用类和构造函数实现输入地面点坐标，在给定内外方位元素施输出像点坐标的C++代码

以下是一个简单的示例代码，其中包含一个名为"Point"的类和一个名为"Camera"的类，用于计算像点坐标。

#include <iostream>
#include <cmath>
using namespace std;

class Point {
public:
    double x, y, z;
    Point(double xx, double yy, double zz) : x(xx), y(yy), z(zz) {}
};

class Camera {
public:
    double f, k1, k2, px, py, omega, phi, kappa;
    Camera(double ff, double kk1, double kk2, double pxx, double pyy, double om, double ph, double kap) : f(ff), k1(kk1), k2(kk2), px(pxx), py(pyy), omega(om), phi(ph), kappa(kap) {}
    Point WorldToCamera(Point p);
    Point CameraToImage(Point p);
};

Point Camera::WorldToCamera(Point p) {
    double x1 = p.x, y1 = p.y, z1 = p.z;
    double x2 = cos(phi)*x1 + sin(phi)*z1, y2 = y1, z2 = -sin(phi)*x1 + cos(phi)*z1;
    double x3 = x2, y3 = cos(omega)*y2 - sin(omega)*z2, z3 = sin(omega)*y2 + cos(omega)*z2;
    double x4 = cos(kappa)*x3 + sin(kappa)*y3, y4 = -sin(kappa)*x3 + cos(kappa)*y3, z4 = z3;
    Point p2(x4, y4, z4);
    return p2;
}

Point Camera::CameraToImage(Point p) {
    double x1 = p.x / p.z, y1 = p.y / p.z;
    double r = sqrt(x1*x1 + y1*y1);
    double x2 = x1*(1 + k1*r*r + k2*r*r*r*r);
    double y2 = y1*(1 + k1*r*r + k2*r*r*r*r);
    double u = f * x2 + px, v = f * y2 + py;
    Point p2(u, v, 0);
    return p2;
}

int main() {
    Point p1(100, 200, 300);
    Camera c1(1000, 0.01, 0.001, 500, 500, 0.1, 0.2, 0.3);
    Point p2 = c1.WorldToCamera(p1);
    Point p3 = c1.CameraToImage(p2);
    cout << "World coordinate: (" << p1.x << ", " << p1.y << ", " << p1.z << ")" << endl;
    cout << "Camera coordinate: (" << p2.x << ", " << p2.y << ", " << p2.z << ")" << endl;
    cout << "Image coordinate: (" << p3.x << ", " << p3.y << ")" << endl;
    return 0;
}

在此示例代码中，Point 类表示三维坐标点，Camera 类表示相机，包括焦距、径向畸变系数、主点坐标、内外方位元素等参数。WorldToCamera 函数将世界坐标转换为相机坐标，CameraToImage 函数将相机坐标转换为像点坐标。在 main 函数中，我们创建一个 Point 对象表示世界坐标，创建一个 Camera 对象表示相机，然后使用 WorldToCamera 函数将世界坐标转换为相机坐标，再使用 CameraToImage 函数将相机坐标转换为像点坐标。最后将结果输出到控制台。

请注意，此代码仅为示例，实际应用中需要根据具体的需求进行修改和优化。