gdalwarp를 사용하는 방법은  http://www.gdal.org/gdalwarp.html 에 자세히 나와있습니다.


요 포스팅에서는 gdalwarp를 이용해서 좌표계를 변경하는 방법에 대해서 한글로 간단하게 정리하는데 의미가 있습니다. 


설명을 위해 A.tif 라는 파일을 생성해뒀습니다.  

A.tif 라는 파일의 좌표계는 EPSG 4326입니다. 

gdalinfo 명령어로도 확인이 가능하죠 아래와 같이 말이에요.


gdalinfo  A.tif


Driver: GTiff/GeoTIFF

Files: Original00.tif

Size is 2000, 2000

Coordinate System is:

GEOGCS["WGS 84",

    DATUM["WGS_1984",

        SPHEROID["WGS 84",6378137,298.257223563,

            AUTHORITY["EPSG","7030"]],

        AUTHORITY["EPSG","6326"]],

    PRIMEM["Greenwich",0],

    UNIT["degree",0.0174532925199433],

    AUTHORITY["EPSG","4326"]]


EPSG 4236 좌표계로 설정되어 있는 A.tif 파일을 EPSG 3857로 변경하기 위해서는 gdalwarp 명령어를 사용해야 합니다.


gdalwarp에 help를 보시면


Usage: gdalwarp [--help-general] [--formats]

    [-s_srs srs_def] [-t_srs srs_def] [-to "NAME=VALUE"]

    .........................

    srcfile* dstfile



일케 나와 있는데요. 이중에서 좌표계 변경은 s_srs와 t_srs 옵션을 사용해서 처리합니다. 

s_srs는 입력 영상의 좌표계 정보를 입력하는 것이고, t_srs는 목표 영상의 좌표계 정보를 입력하는 것입니다. 

즉 EPSG 4236은 s_srs에 입력해야 하는 정보이고, t_srs는 EPSG 3857이라는 값을 입력해야 목표하는 일을 처리할 수 있습니다. 


다음과 같이 말이죠.


gdalwarp -t_srs EPSG:3857 A.tif B.tif 


위 명령어에서는 s_srs를 사용하지 않았는데요. s_srs는 이미 A.tif에 입력되어 있기 때문에 사용하지 않아도 괜찮습니다. 

위 명령어를 통해 생성되는 B.tif는 EPSG 3857 좌표계로 설정됩니다. gdalinfo로 보면 다음과 같이 출력됩니다. 


gdalinfo B.tif


Driver: GTiff/GeoTIFF

Files: B.tif

Size is 5230, 5158

Coordinate System is:

PROJCS["WGS 84 / Pseudo-Mercator",

    GEOGCS["WGS 84",

        DATUM["WGS_1984",

            SPHEROID["WGS 84",6378137,298.257223563,

                AUTHORITY["EPSG","7030"]],

            AUTHORITY["EPSG","6326"]],

        PRIMEM["Greenwich",0],

        UNIT["degree",0.0174532925199433],

        AUTHORITY["EPSG","4326"]],

    PROJECTION["Mercator_1SP"],

    PARAMETER["central_meridian",0],

    PARAMETER["scale_factor",1],

    PARAMETER["false_easting",0],

    PARAMETER["false_northing",0],

    UNIT["metre",1,

        AUTHORITY["EPSG","9001"]],

    EXTENSION["PROJ4","+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext  +no_defs"],

    AUTHORITY["EPSG","3857"]]





+ Recent posts