QGIS와 Python을 이용해서 raster file 일정 크기의 타일로 분할하기

이 글은 래스터 파일을 일정 크기로 분할하는 것입니다. QGIS를 직접 사용하거나 심지어 Python(GDAL)을 사용하여 이 작업을 수행할 수도 있다.
더 자세히 설명하자면, 256개의 열과 256개의 행의 크기로 분할하려고 할때 유용합니다.

 

이를 위해서 래스터 파일을 마우스 오른쪽 버튼을 클릭해서 Export >> 다른 이름으로 저장 (A)...를 선택한다.

 

 

래스터 레이어를 다른 이름으로 저장... 대화상자에서 다음 그림처럼 실행하면 된다.

 

 

반응형

 

 

최종결과는 tile_map 폴더에 256 X 256 크기로 분할된 파일이 저장된다.

 

이를 Python GDAL을 이용한 코드는 다음과 같다.

import os, gdal
 
in_path = 'C:/python_work/landcover_unet_mydata/mydata/images/'
input_filename = 'sentinel_band.tif'
 
out_path = 'C:/python_work/landcover_unet_mydata/mydata/tile_map/'
output_filename = 'tile_'
 
tile_size_x = 256
tile_size_y = 256
 
ds = gdal.Open(in_path + input_filename)
band = ds.GetRasterBand(1)
xsize = band.XSize
ysize = band.YSize
 
for i in range(0, xsize, tile_size_x):
    for j in range(0, ysize, tile_size_y):
        com_string = "gdal_translate -of GTIFF -srcwin " + str(i)+ ", " + str(j) + ", " + str(tile_size_x) + ", " + str(tile_size_y) + " " + str(in_path) + str(input_filename) + " " + str(out_path) + str(output_filename) + str(i) + "_" + str(j) + ".tif"
        os.system(com_string)