MapView 경위도 tm 카텍 변환 알고리즘

다음 맵은 수익이 발생하는 맵에서는 쓸 수 없다고 이용약관에 나와있고 

네이버 맵은 그런 내용이 없는데 ...

그래서 네이버 맵을 사용할려니 이놈은 tm좌표이다 

tm좌표 이건 일본 식민지 시절때 받은거 아직 쓰고 잇단다 

어디서는 관동대지진이후 tm 좌표가 바뀌었다고 하는데 그런 자세한 것까진 모르겠고 

네이버 맵에서 네이버 지역 검색후 나오는 mapx mapy 의 좌표값을 샘플 예제에 넣고 돌리면 

경위도가 나오는 샘플 예제 소스다 

GeoTrans.zip

http://javaexpert.tistory.com/142  이쪽이 원본 소스 

http://www.androidpub.com/1318647 이쪽이 원본 소스를 좀 더 설명한 페이지

두분다 엄청 고마우신 분 ㅠㅠ

objective c 사용 소스

GeoTransPoint *oTM = [[GeoTransPoint alloc]init];

        [oTM setX:301255];

        [oTM setY:540215];

        

        GeoTrans *oGeo = [[GeoTrans alloc]init];

        GeoTransPoint *point = [oGeo convert:KATEC DST:GEO GEOIN:oTM];

       lat = [point getY] * 1E6;

    log = [point getX] * 1E6;

        NSLog(@"lat %f lng %f  %f %f",lat , log,[point getX],[point getY]);

        

        lat = [point getY];

        log = [point getX];

               

       

        [oTM release];

        [oGeo release];

//

//  GeoTrans.m

//  FEvent

//

//  Created by jinuk son on 13. 3. 14..

//  Copyright (c) 2013년 DANAL. All rights reserved.

//

#import "GeoTrans.h"

#import <math.h>

@implementation GeoTrans

-(id)init

{

    self = [super init];

    if(self)

    {

        m_arScaleFactor[GEO] = 1;

        m_arLonCenter[GEO] = 0.0;

        m_arLatCenter[GEO] = 0.0;

        m_arFalseNorthing[GEO] = 0.0;

        m_arFalseEasting[GEO] = 0.0;

        m_arMajor[GEO] = 6378137.0;

        m_arMinor[GEO] = 6356752.3142;

        

        m_arScaleFactor[KATEC] = 0.9999;//0.9999;

        //m_arLonCenter[KATEC] = 2.22529479629277; // 127.5

        m_arLonCenter[KATEC] = 2.23402144255274; // 128

        m_arLatCenter[KATEC] = 0.663225115757845;

        m_arFalseNorthing[KATEC] = 600000.0;

        m_arFalseEasting[KATEC] = 400000.0;

        m_arMajor[KATEC] = 6377397.155;

        m_arMinor[KATEC] = 6356078.9633422494;

        

        m_arScaleFactor[TM] = 1.0;

        //this.m_arLonCenter[TM] = 2.21656815003280; // 127

        m_arLonCenter[TM] = 2.21661859489671; // 127.+10.485 minute

        m_arLatCenter[TM] = 0.663225115757845;

        m_arFalseNorthing[TM] = 500000.0;

        m_arFalseEasting[TM] = 200000.0;

        m_arMajor[TM] = 6377397.155;

        m_arMinor[TM] = 6356078.9633422494;

        

        datum_params[0] = -146.43;

        datum_params[1] = 507.89;

        datum_params[2] = 681.46;

        

        double tmp = m_arMinor[GEO] / m_arMajor[GEO];

        m_Es[GEO] = 1.0 - tmp * tmp;

        m_Esp[GEO] = m_Es[GEO] / (1.0 - m_Es[GEO]);

        

        if (m_Es[GEO] < 0.00001) {

            m_Ind[GEO] = 1.0;

        } else {

            m_Ind[GEO] = 0.0;

        }

        

        tmp = m_arMinor[KATEC] / m_arMajor[KATEC];

        m_Es[KATEC] = 1.0 - tmp * tmp;

        m_Esp[KATEC] = m_Es[KATEC] / (1.0 - m_Es[KATEC]);

        

        if (m_Es[KATEC] < 0.00001) {

            m_Ind[KATEC] = 1.0;

        } else {

            m_Ind[KATEC] = 0.0;

        }

        

        tmp = m_arMinor[TM] / m_arMajor[TM];

        m_Es[TM] = 1.0 - tmp * tmp;

        m_Esp[TM] = m_Es[TM] / (1.0 - m_Es[TM]);

        

        if (m_Es[TM] < 0.00001) {

            m_Ind[TM] = 1.0;

        } else {

            m_Ind[TM] = 0.0;

        }

        

        [self mlfn:[self e0fn:m_Es[TM]] E1:[self e1fn:m_Es[TM]] E2:[self e2fn:m_Es[TM]] E3:[self e3fn:m_Es[TM]] E4:m_arLatCenter[TM]];

        

        

        src_m[GEO] = m_arMajor[GEO] *  [self mlfn:[self e0fn:m_Es[GEO]] E1:[self e1fn:m_Es[GEO]] E2:[self e2fn:m_Es[GEO]] E3:[self e3fn:m_Es[GEO]] E4:m_arLatCenter[GEO]];

        dst_m[GEO] = m_arMajor[GEO] *  [self mlfn:[self e0fn:m_Es[GEO]] E1:[self e1fn:m_Es[GEO]] E2:[self e2fn:m_Es[GEO]] E3:[self e3fn:m_Es[GEO]] E4:m_arLatCenter[GEO]];

        src_m[KATEC] = m_arMajor[KATEC] * [self mlfn:[self e0fn:m_Es[KATEC]] E1:[self e1fn:m_Es[KATEC]] E2:[self e2fn:m_Es[KATEC]] E3:[self e3fn:m_Es[KATEC]] E4:m_arLatCenter[KATEC]];

        dst_m[KATEC] = m_arMajor[KATEC] * [self mlfn:[self e0fn:m_Es[KATEC]] E1:[self e1fn:m_Es[KATEC]] E2:[self e2fn:m_Es[KATEC]] E3:[self e3fn:m_Es[KATEC]] E4:m_arLatCenter[KATEC]];

        src_m[TM] = m_arMajor[TM] * [self mlfn:[self e0fn:m_Es[TM]] E1:[self e1fn:m_Es[TM]] E2:[self e2fn:m_Es[TM]] E3:[self e3fn:m_Es[TM]] E4:m_arLatCenter[TM]];

        dst_m[TM] = m_arMajor[TM] * [self mlfn:[self e0fn:m_Es[TM]] E1:[self e1fn:m_Es[TM]] E2:[self e2fn:m_Es[TM]] E3:[self e3fn:m_Es[TM]] E4:m_arLatCenter[TM]];

    }

    return self;

    

}

-(void)initValue {

    

 }

-(double)D2R:(double) degree {

    return degree* M_PI / 180.0;

}

-(double)R2D:(double) radian {

    return radian * 180.0 / M_PI;

}

-(double)e0fn:(double) x {

    return 1.0 - 0.25 * x * (1.0 + x / 16.0 * (3.0 + 1.25 * x));

}

-(double)e1fn:(double) x {

    return 0.375 * x * (1.0 + 0.25 * x * (1.0 + 0.46875 * x));

}

-(double)e2fn:(double) x {

    return 0.05859375 * x * x * (1.0 + 0.75 * x);

}

-(double)e3fn:(double) x {

    return x * x * x * (35.0 / 3072.0);

}

-(double)mlfn:(double) e0  E1:(double) e1 E2:(double) e2 E3:(double)e3 E4:( double) phi

{

    return e0 * phi - e1 * sin(2.0 * phi) + e2 * sin(4.0 * phi) - e3 * sin(6.0 * phi);

}

-(double)asinz:(double) value {

    if (abs(value) > 1.0) value = (value > 0 ? 1: -1);

    return asin(value);

}

-(GeoTransPoint *) convert:(int) srctype DST:(int)dsttype GEOIN:(GeoTransPoint*) in_pt{

    GeoTransPoint *tmpPt = [[GeoTransPoint alloc]init];

    GeoTransPoint *out_pt = [[[GeoTransPoint alloc]init] autorelease];

    

    

    if (srctype == GEO) {

        tmpPt.x = [self D2R:in_pt.x];

        tmpPt.y = [self D2R:in_pt.y];

    } else {

        [self tm2geo:srctype GEOIN:in_pt GEOOUT:tmpPt];

    }

    

    if (dsttype == GEO) {

        out_pt.x = [self R2D:tmpPt.x];

        out_pt.y = [self R2D:tmpPt.y];

    } else {

        [self geo2tm:dsttype GEOIN:tmpPt GEOOUT:out_pt];

        //out_pt.x = Math.round(out_pt.x);

        //out_pt.y = Math.round(out_pt.y);

    }

    [tmpPt release];

    return out_pt;

}

-(void) geo2tm:(int) dsttype  GEOIN:(GeoTransPoint *)in_pt GEOOUT:(GeoTransPoint *)out_pt {

    double x, y;

    

    [self transform:GEO DST:dsttype GEOPOINT:in_pt];

    double delta_lon = in_pt.x - m_arLonCenter[dsttype];

    double sin_phi = sin(in_pt.y);

    double cos_phi = cos(in_pt.y);

    

    if (m_Ind[dsttype] != 0) {

        double b = cos_phi * sin(delta_lon);

        

        if ((abs(abs(b) - 1.0)) < EPSLN) {

            //Log.d("무한대 에러");

            //System.out.println("무한대 에러");

        }

    } else {

        double b = 0;

        x = 0.5 * m_arMajor[dsttype] * m_arScaleFactor[dsttype] * log((1.0 + b) / (1.0 - b));

        double con = acos(cos_phi * cos(delta_lon) / sqrt(1.0 - b * b));

        

        if (in_pt.y < 0) {

            con = con * -1;

            y = m_arMajor[dsttype] * m_arScaleFactor[dsttype] * (con - m_arLatCenter[dsttype]);

        }

    }

    

    double al = cos_phi * delta_lon;

    double als = al * al;

    double c = m_Esp[dsttype] * cos_phi * cos_phi;

    double tq = tan(in_pt.y);

    double t = tq * tq;

    double con = 1.0 - m_Es[dsttype] * sin_phi * sin_phi;

    double n = m_arMajor[dsttype] / sqrt(con);

    double ml = m_arMajor[dsttype] * [self mlfn:[self e0fn:m_Es[dsttype]] E1:[self e1fn:m_Es[dsttype]] E2:[self e2fn:m_Es[dsttype]] E3:[self e3fn:m_Es[dsttype]] E4: in_pt.y];

    

    

    

    

    

    out_pt.x = m_arScaleFactor[dsttype] * n * al * (1.0 + als / 6.0 * (1.0 - t + c + als / 20.0 * (5.0 - 18.0 * t + t * t + 72.0 * c - 58.0 * m_Esp[dsttype]))) + m_arFalseEasting[dsttype];

    out_pt.y = m_arScaleFactor[dsttype] * (ml - dst_m[dsttype] + n * tq * (als * (0.5 + als / 24.0 * (5.0 - t + 9.0 * c + 4.0 * c * c + als / 30.0 * (61.0 - 58.0 * t + t * t + 600.0 * c - 330.0 * m_Esp[dsttype]))))) + m_arFalseNorthing[dsttype];

}

-(void) tm2geo:(int)srctype GEOIN:(GeoTransPoint *)in_pt GEOOUT:(GeoTransPoint *)out_pt{

    

    

    GeoTransPoint *tmpPt = [[ GeoTransPoint alloc]init];

    

    [tmpPt GeoTransPoint2:[in_pt getX] Y:[in_pt getY]];

    

    int max_iter = 6;

    

    

    if (m_Ind[srctype] != 0) {

        double f = exp(in_pt.x / (m_arMajor[srctype] * m_arScaleFactor[srctype]));

        double g = 0.5 * (f - 1.0 / f);

        double temp = m_arLatCenter[srctype] + tmpPt.y / (m_arMajor[srctype] * m_arScaleFactor[srctype]);

        double h = cos(temp);

        double con = sqrt((1.0 - h * h) / (1.0 + g * g));

        out_pt.y = [self asinz:con];

        

        if (temp < 0) out_pt.y *= -1;

        

        if ((g == 0) && (h == 0)) {

            out_pt.x = m_arLonCenter[srctype];

        } else {

            out_pt.x = atan(g / h) + m_arLonCenter[srctype];

        }

    }

    

    tmpPt.x -= m_arFalseEasting[srctype];

    tmpPt.y -= m_arFalseNorthing[srctype];

    

    double con = (src_m[srctype] + tmpPt.y / m_arScaleFactor[srctype]) / m_arMajor[srctype];

    double phi = con;

    

    int i = 0;

    

    while (true) {

        double delta_Phi = ((con + [self e1fn:m_Es[srctype]] * sin(2.0 * phi) - [self e2fn:m_Es[srctype]] * sin(4.0 * phi) + [self e3fn:m_Es[srctype]] * sin(6.0 * phi)) / [self e0fn:m_Es[srctype]]) - phi;

        phi = phi + delta_Phi;

        

        if (abs(delta_Phi) <= EPSLN) break;

        

        if (i >= max_iter) {

            //Log.d("무한대 에러");

            //System.out.println("무한대 에러");

            break;

        }

        

        i++;

    }

    

    if (abs(phi) < (M_PI / 2)) {

        double sin_phi = sin(phi);

        double cos_phi = cos(phi);

        double tan_phi = tan(phi);

        double c = m_Esp[srctype] * cos_phi * cos_phi;

        double cs = c * c;

        double t = tan_phi * tan_phi;

        double ts = t * t;

        double cont = 1.0 - m_Es[srctype] * sin_phi * sin_phi;

        double n = m_arMajor[srctype] / sqrt(cont);

        double r = n * (1.0 - m_Es[srctype]) / cont;

        double d = tmpPt.x / (n * m_arScaleFactor[srctype]);

        double ds = d * d;

        out_pt.y = phi - (n * tan_phi * ds / r) * (0.5 - ds / 24.0 * (5.0 + 3.0 * t + 10.0 * c - 4.0 * cs - 9.0 * m_Esp[srctype] - ds / 30.0 * (61.0 + 90.0 * t + 298.0 * c + 45.0 * ts - 252.0 * m_Esp[srctype] - 3.0 * cs)));

        out_pt.x = m_arLonCenter[srctype] + (d * (1.0 - ds / 6.0 * (1.0 + 2.0 * t + c - ds / 20.0 * (5.0 - 2.0 * c + 28.0 * t - 3.0 * cs + 8.0 * m_Esp[srctype] + 24.0 * ts))) / cos_phi);

    } else {

        out_pt.y = M_PI * 0.5 * sin(tmpPt.y);

        out_pt.x = m_arLonCenter[srctype];

    }

    [self transform:srctype DST:GEO GEOPOINT:out_pt];

}

-(double)getDistancebyGeo:(GeoTransPoint *) pt1  GEO2:(GeoTransPoint *) pt2 {

    double lat1 = [self D2R:pt1.y];

    double lon1 = [self D2R:pt1.x];

    double lat2 = [self D2R:pt2.y];

    double lon2 = [self D2R:pt2.x];

    

    double longitude = lon2 - lon1;

    double latitude = lat2 - lat1;

    

    double a = pow(sin(latitude / 2.0), 2) + cos(lat1) * cos(lat2) * pow(sin(longitude / 2.0), 2);

    return 6376.5 * 2.0 * atan2(sqrt(a), sqrt(1.0 - a));

}

-(double)getDistancebyKatec:(GeoTransPoint *) pt1  GEO2:(GeoTransPoint *) pt2 {

   

    pt1 =  [self convert:KATEC DST:GEO GEOIN:pt1];

    pt2 =  [self convert:KATEC DST:GEO GEOIN:pt2];

    

    return [self getDistancebyGeo:pt1 GEO2:pt2];

}

-(double)getDistancebyTm:(GeoTransPoint *) pt1  GEO2:(GeoTransPoint *) pt2 {

      

    

    pt1 =  [self convert:TM DST:GEO GEOIN:pt1];

    pt2 =  [self convert:TM DST:GEO GEOIN:pt2];

    

    return [self getDistancebyGeo:pt1 GEO2:pt2];

    

    

}

-(long) getTimebySec:(double) distance {

    return round(3600 * distance / 4);

}

-(long) getTimebyMin:(double) distance {

    return (long)(ceil(    [self getTimebySec:distance] / 60));

}

/*

 Author:       Richard Greenwood rich@greenwoodmap.com

 License:      LGPL as per: http://www.gnu.org/copyleft/lesser.html

 */

/**

 * convert between geodetic coordinates (longitude, latitude, height)

 * and gecentric coordinates (X, Y, Z)

 * ported from Proj 4.9.9 geocent.c

 */

// following constants from geocent.c

double HALF_PI = 0.5 * M_PI;

double COS_67P5  = 0.38268343236508977;  /* cosine of 67.5 degrees */

double AD_C      = 1.0026000 ;

/* Toms region 1 constant */

-(void) transform:(int)srctype DST:(int)dsttype GEOPOINT:(GeoTransPoint *) point {

    if (srctype == dsttype)

        return;

    

    if (srctype != 0 || dsttype != 0) {

        // Convert to geocentric coordinates.

        [self geodetic_to_geocentric:srctype GEOPOINT:point];

        

        // Convert between datums

        if (srctype != 0) {

            [self geocentric_to_wgs84:point];

        }

        

        if (dsttype != 0) {

            [self geocentric_from_wgs84:point];

        }

        

        // Convert back to geodetic coordinates

        [self geocentric_to_geodetic:dsttype GEOPOINT:point];

    }

}

-(Boolean) geodetic_to_geocentric:(int) type  GEOPOINT:(GeoTransPoint *) p {

    

    /*

     * The function Convert_Geodetic_To_Geocentric converts geodetic coordinates

     * (latitude, longitude, and height) to geocentric coordinates (X, Y, Z),

     * according to the current ellipsoid parameters.

     *

     *    Latitude  : Geodetic latitude in radians                     (input)

     *    Longitude : Geodetic longitude in radians                    (input)

     *    Height    : Geodetic height, in meters                       (input)

     *    X         : Calculated Geocentric X coordinate, in meters    (output)

     *    Y         : Calculated Geocentric Y coordinate, in meters    (output)

     *    Z         : Calculated Geocentric Z coordinate, in meters    (output)

     *

     */

    

    double Longitude = p.x;

    double Latitude = p.y;

    double Height = p.z;

    double X;  // output

    double Y;

    double Z;

    

    double Rn;            /*  Earth radius at location  */

    double Sin_Lat;       /*  Math.sin(Latitude)  */

    double Sin2_Lat;      /*  Square of Math.sin(Latitude)  */

    double Cos_Lat;       /*  Math.cos(Latitude)  */

    

    /*

     ** Don't blow up if Latitude is just a little out of the value

     ** range as it may just be a rounding issue.  Also removed longitude

     ** test, it should be wrapped by Math.cos() and Math.sin().  NFW for PROJ.4, Sep/2001.

     */

    if (Latitude < -HALF_PI && Latitude > -1.001 * HALF_PI )

        Latitude = -HALF_PI ;

    else if (Latitude > HALF_PI && Latitude < 1.001 * HALF_PI )

        Latitude = HALF_PI;

    else if ((Latitude < -HALF_PI) || (Latitude > HALF_PI)) { /* Latitude out of range */

        return true;

    }

    

    /* no errors */

    if (Longitude > M_PI)

        Longitude -= (2*M_PI);

    Sin_Lat = sin(Latitude);

    Cos_Lat = cos(Latitude);

    Sin2_Lat = Sin_Lat * Sin_Lat;

    Rn = m_arMajor[type] / (sqrt(1.0e0 - m_Es[type] * Sin2_Lat));

    X = (Rn + Height) * Cos_Lat * cos(Longitude);

    Y = (Rn + Height) * Cos_Lat * sin(Longitude);

    Z = ((Rn * (1 - m_Es[type])) + Height) * Sin_Lat;

    

    p.x = X;

    p.y = Y;

    p.z = Z;

    return false;

} // cs_geodetic_to_geocentric()

/** Convert_Geocentric_To_Geodetic

 * The method used here is derived from 'An Improved Algorithm for

 * Geocentric to Geodetic Coordinate Conversion', by Ralph Toms, Feb 1996

 */

-(void) geocentric_to_geodetic :(int) type GEOPOINT:(GeoTransPoint *) p {

    

    double X = p.x;

    double Y = p.y;

    double Z = p.z;

    double Longitude;

    double Latitude = 0.;

    double Height;

    

    double W;        /* distance from Z axis */

    double W2;       /* square of distance from Z axis */

    double T0;       /* initial estimate of vertical component */

    double T1;       /* corrected estimate of vertical component */

    double S0;       /* initial estimate of horizontal component */

    double S1;       /* corrected estimate of horizontal component */

    double Sin_B0;   /* Math.sin(B0), B0 is estimate of Bowring aux doubleiable */

    double Sin3_B0;  /* cube of Math.sin(B0) */

    double Cos_B0;   /* Math.cos(B0) */

    double Sin_p1;   /* Math.sin(phi1), phi1 is estimated latitude */

    double Cos_p1;   /* Math.cos(phi1) */

    double Rn;       /* Earth radius at location */

    double Sum;      /* numerator of Math.cos(phi1) */

    Boolean At_Pole;  /* indicates location is in polar region */

    

    At_Pole = false;

    if (X != 0.0) {

        Longitude = atan2(Y,X);

    }

    else {

        if (Y > 0) {

            Longitude = HALF_PI;

        }

        else if (Y < 0) {

            Longitude = -HALF_PI;

        }

        else {

            At_Pole = true;

            Longitude = 0.0;

            if (Z > 0.0) {  /* north pole */

                Latitude = HALF_PI;

            }

            else if (Z < 0.0) {  /* south pole */

                Latitude = -HALF_PI;

            }

            else {  /* center of earth */

                Latitude = HALF_PI;

                Height = -m_arMinor[type];

                return;

            }

        }

    }

    W2 = X*X + Y*Y;

    W = sqrt(W2);

    T0 = Z * AD_C;

    S0 = sqrt(T0 * T0 + W2);

    Sin_B0 = T0 / S0;

    Cos_B0 = W / S0;

    Sin3_B0 = Sin_B0 * Sin_B0 * Sin_B0;

    T1 = Z + m_arMinor[type] * m_Esp[type] * Sin3_B0;

    Sum = W - m_arMajor[type] * m_Es[type] * Cos_B0 * Cos_B0 * Cos_B0;

    S1 = sqrt(T1*T1 + Sum * Sum);

    Sin_p1 = T1 / S1;

    Cos_p1 = Sum / S1;

    Rn = m_arMajor[type] / sqrt(1.0 - m_Es[type] * Sin_p1 * Sin_p1);

    if (Cos_p1 >= COS_67P5) {

        Height = W / Cos_p1 - Rn;

    }

    else if (Cos_p1 <= -COS_67P5) {

        Height = W / -Cos_p1 - Rn;

    }

    else {

        Height = Z / Sin_p1 + Rn * (m_Es[type] - 1.0);

    }

    if (At_Pole == false) {

        Latitude = atan(Sin_p1 / Cos_p1);

    }

    

    p.x = Longitude;

    p.y =Latitude;

    p.z = Height;

    return;

} // geocentric_to_geodetic()

/****************************************************************/

// geocentic_to_wgs84(defn, p )

//  defn = coordinate system definition,

//  p = point to transform in geocentric coordinates (x,y,z)

-(void) geocentric_to_wgs84:(GeoTransPoint *) p {

    

    //if( defn.datum_type == PJD_3PARAM )

    {

        // if( x[io] == HUGE_VAL )

        //    continue;

        p.x += datum_params[0];

        p.y += datum_params[1];

        p.z += datum_params[2];

    }

} // geocentric_to_wgs84

/****************************************************************/

// geocentic_from_wgs84()

//  coordinate system definition,

//  point to transform in geocentric coordinates (x,y,z)

-(void) geocentric_from_wgs84:(GeoTransPoint *) p {

    

    //if( defn.datum_type == PJD_3PARAM )

    {

        //if( x[io] == HUGE_VAL )

        //    continue;

        p.x -= datum_params[0];

        p.y -= datum_params[1];

        p.z -= datum_params[2];

        

    }

} //geocentric_from_wgs84()

@end

//

//  GeoTrans.h

//  FEvent

//

//  Created by jinuk son on 13. 3. 14..

//  Copyright (c) 2013년 DANAL. All rights reserved.

//

#import <Foundation/Foundation.h>

#import "GeoTransPoint.h"

#define GEO 0

#define KATEC 1

#define TM 2

#define GRS80 3

#define EPSLN 0.0000000001

@interface GeoTrans : NSObject

{

    

    

    double m_Ind [3];

double m_Es [3];

double m_Esp [3];

double src_m [3];

double dst_m [3];

    

double m_arMajor [3];

double m_arMinor [3];

    

double m_arScaleFactor [3];

double m_arLonCenter [3];

double m_arLatCenter [3];

double m_arFalseNorthing [3];

double m_arFalseEasting [3];

double datum_params [3];

}

-(void)initValue;

-(GeoTransPoint *) convert:(int) srctype DST:(int)dsttype GEOIN:(GeoTransPoint*) in_pt;

@end

//

//  GeoTransPoint.m

//  FEvent

//

//  Created by jinuk son on 13. 3. 14..

//  Copyright (c) 2013년 DANAL. All rights reserved.

//

#import "GeoTransPoint.h"

@implementation GeoTransPoint

@synthesize x,y,z;

-(void) GeoTransPoint2:(double) _x Y:(double) _y

{

    

    x = _x;

    y = _y;

    z = 0;

}

-(void) GeoTransPoint3:(double) _x Y:(double) _y Z:(double)_z {

    

    x = _x;

    y = _y;

    z = 0;

}

-(double) getX {

    return x;

}

-(double) getY {

    return y;

}

@end

//

//  GeoTransPoint.h

//  FEvent

//

//  Created by jinuk son on 13. 3. 14..

//  Copyright (c) 2013년 DANAL. All rights reserved.

//

#import <Foundation/Foundation.h>

@interface GeoTransPoint : NSObject

{

    

}

@property(nonatomic,readwrite)double x;

@property(nonatomic,readwrite)double y;

@property(nonatomic,readwrite)double z;

-(void) GeoTransPoint2:(double) _x Y:(double) _y;

-(void) GeoTransPoint3:(double) _x Y:(double) _y Z:(double)_z ;

-(double) getX;

-(double) getY ;

@end