火星02坐标转换为WGS84坐标

import java.io.BufferedReader;
import java.io.File;
import java.io.FileNotFoundException;
import java.io.FileReader;
import java.io.IOException;

public class Gpstest {

    final static int TABLESIZE=660*450;
    //#define ID(i, j)  ((i)  +  660 * (j))

    static double[] TableX=new double[TABLESIZE];
    static double[] TableY=new double[TABLESIZE];
    static boolean bInitTable = false;

    static void InitTable() throws IOException
    {
        File filename = new File("C://Mars2Wgsnew.txt");
        long lX1, lY1;
        String tempString;
        String[] tem;
        bInitTable = false;
        BufferedReader reader = new BufferedReader(new FileReader(filename));
        for(int i=0;;i++){
            tempString = reader.readLine();
            if(tempString==null){
                 break;
             }
    //          System.out.println(tempString);
            tempString=tempString.replaceAll(" ","");
            tem=tempString.split(",");

            TableX[i] = ((double) Long.valueOf(tem[0])) / 100000.0;
            TableY[i] = ((double) Long.valueOf(tem[1])) / 100000.0;
        }
        reader.close();
        bInitTable = true;
    }

    // Result:
//             0 - ok
//            -1 - Table is not initialized
//            -2 - Latitude or Longitude is invalid
    static int Mars2Wgs( double x_mars, double y_mars )
    {
        int i, j, k;
        double x1, y1, x2, y2, x3, y3, x4, y4, xtry, ytry, dx, dy,x_wgs,y_wgs;
        double t, u;

        if(!bInitTable) return -1;

        xtry = x_mars;
        ytry = y_mars;

        for( k=0; k<10; ++k )
        {
            // 只对中国国境内数据转换
            if( xtry < 72 || xtry > 137.9 || ytry < 10 || ytry > 54.9)
            {
                return -2;
            }
            i = (int) ((xtry - 72.0) * 10.0);
            j = (int) ((ytry - 10.0) * 10.0);
            x1 = TableX[i+660*j];
            y1 = TableY[i+660*j];
            x2 = TableX[(i+1)+660*j];
            y2 = TableY[(i+1)+660*j];
            x3 = TableX[(i+1)+660*(j+1)];
            y3 = TableY[(i+1)+660*(j+1)];
            x4 = TableX[i+660*(j+1)];
            y4 = TableY[i+660*(j+1)];
            t = (xtry - 72.0 - 0.1 * i) * 10.0;
            u = (ytry - 10.0 - 0.1 * j) * 10.0;
            dx = (1.0-t)*(1.0-u)*x1 + t*(1.0-u)*x2 + t*u*x3 + (1.0-t)*u*x4 - xtry;
            dy = (1.0-t)*(1.0-u)*y1 + t*(1.0-u)*y2 + t*u*y3 + (1.0-t)*u*y4 - ytry;
            xtry = (xtry + x_mars - dx)/2.0;
            ytry = (ytry + y_mars - dy)/2.0;
        }
        x_wgs = xtry;
        y_wgs = ytry;
        System.out.println(x_wgs+";"+y_wgs);
        return 0;
    }

    public static void main(String[] args) throws IOException
    {
        double x_mars, y_mars, x_wgs, y_wgs;
        File fname = new File("C://2.txt");
        String tempString;
        String[] tem;
        BufferedReader reader = new BufferedReader(new FileReader(fname));
        int result;
        for(int i=0;;i++){
            tempString = reader.readLine();
            if(tempString==null){
                 break;
             }
    //          System.out.println(tempString);
            tempString=tempString.replaceAll(" ","");
            tem=tempString.split(",");

            x_mars = ((double) Double.valueOf(tem[0]));
            y_mars = ((double) Double.valueOf(tem[1]));
            InitTable();
            result = Mars2Wgs(x_mars, y_mars);
        }
    }
}

本例子程序改编自其它网友，其功能是将火星02坐标转换为WGS84坐标

需要一个数据文件Mars2Wgsnew.txt。转换后存在几米到十几米的误差

http://files.cnblogs.com/casicyuan/Mars2Wgsnew.rar