일상다반사 - 개발

gina-alaska에서 제공하는 dans-gdal-scripts 사용 방법 1. 기본사용 방법    a. GitHub에서 dans-gdal-scripts를 검색하면 기본적으로 우분투를 사용해서 설치하는 방법을 제공한다.    b. 우분투 사용자라면 쉽게 설치를 할 수 있다. 하지만 centos 사용자에 대한 설치 방법은 제공하고 있지 않기 때문에 아래와 같이 하면 해결할 수 있다. 문제 : centos에서 dans-gdal-scripts 설치 안되는 문제 해결방법 : dans-gdal-scripts를 centos에서 사용하려면 기본적으로 gdal 과 boost lib를 설치해야 한다. (gdal의 사용법은 본 블로그에서 검색하면 알 수 있다.) boost lib는 yum으로 설치할 수 있으나 버전의 호환성을 가지기 때문에 나는 소스를 다운로드 받아 설치하였다.  http://www.boost.org/doc/libs/1_46_1/more/getting_started/unix-variants.html 경로를 통하면 boost 소스파일을 다운로드 받을 수 있고 또한 설치 방법도 소개하고 있기 때문에 따라하면 된다. boost 설치방법 1. 다운로드 받은 파일을 압축해제한다. 2. 압축해제한 폴더로 이동한다. 3. ./bootstrap.sh 또는  ./bootstrap.sh --prefix= path / to / installation / prefix 4.  ./bjam install 명령어를 실행하면 정상적으로 설치된다. gdal 과 boost가 정상적으로 설치되면 gans-gdal-scripts를 압축해제한 폴더에서 아래와 같은 명령어를 입력한다. 1. ./autogen.sh 2. ./configure # --prefix=/your/favorite/prefix 3. make 4. make install

> gdal_merge.py -separate -o naples_true_composite.tif LC81890322013210LGN01_B4.TIF LC81890322013210LGN01_B3.TIF LC81890322013210LGN01_B2.TIF 0...10...20...30...40...50...60...70...80...90...100 - done.

Usage: gdallocationinfo [--help-general] [-xml] [-lifonly] [-valonly]                         [-b band]* [-overview overview_level]                         [-l_srs srs_def] [-geoloc] [-wgs84]                         srcfile x y gdallocationinfo LC81150352015290LGN00_B2.TIF -geoloc 439578.358 3990215.231

설치 : sudo pip install pyshp 테스트         sf = shapefile.Reader(path+fname)         shapes = sf.shapes()         print "file len : "+ str(len(shapes))         bbox = shapes[0].bbox         print "bbox value : "         print ['%.3f' % coord for coord in bbox]         print "points len :"+str(len(shapes[0].points))         for index in range(0,100) :             shape = shapes[0].points[index]             #print "point value : 0"             print ['%.3f' % coord for coord in shape]

centos 6.7 에서 python에 gdal lib를 사용할 수 있도록 할 때 : "import gdal" centos 6.7 버전에는 python 버전에 2.6이 기본으로 깔려 있기 때문에 python2.7로 업데이트를 하면 os기본을 2.7로 바꿔줘야 한다. 그렇지 않으면 gdal 설치시(./configure --with-python) 명령어가 정상적으로 실행 되지 않는다. 따라서 gdal 설치시 ./configure 명령을 했을 때 의존성 문제를 항상 체크하고 진행해야 정상적으로 설치할 수 있다.

아래를 실행 했을 때 > gdalwarp -s_srs EPSG:4326 -t_srs EPSG:3857 dest.tif result.tif libproj.so lib가 없다는 에러가 발생했다면, 아래와 같이 해결할 수 있다. ## On Ubuntu sudo apt-get install g++ python-dev libbz2-dev zlib1g-dev ### On CentOS sudo yum install gcc-c++ python-devel bzip2-devel zlib-devel cd /opt/source wget http://download.osgeo.org/proj/proj-4.8.0.tar.gz tar xvfz proj-4.8.0.tar.gz cd proj-4.8.0 mkdir build ./configure --prefix=/opt/source/proj-4.8.0/build # compile make -j$threads # install into build dir make install ### add to /opt/source/scripts/export_paths.sh (see Preliminaries ) ### then run: 'source export_paths.sh' to update these environment paths export LD_LIBRARY_PATH="/opt/source/proj-4.8.0/build/lib:$LD_LIBRARY_PATH" export PATH="/opt/source/proj-4.8.0/build/bin:$PATH # run proj proj # test coordinate transformation # should see coordinates: 70d57'18.075"W 42d15'11.203"N 0.000 echo "338730 4679730...

ubuntu에서 사용할 수있다. 아래 사이트를 참고하면 landsat pansharp를 사용 할 수있다. https://github.com/gina-alaska/dans-gdal-scripts example) gdal_landsat_pansharp -rgb <XXX_B3.TIF> -rgb <XXX_B2.TIF> -rgb <XXX_B1.TIF>       -lum <XXX_B2.TIF> 0.25 -lum <XXX_B3.TIF>.TIF 0.23 -lum <XXX_B4.TIF> 0.52       -pan <XXX_B8.TIF> -ndv 0 -o out.TIF

Tiff 파일중 자신의 관심 영역만을 잘라내는 기능이다. 2가지 방법이 있다. 첫 번째 방법은 이미지 픽셀 값으로 이미지 잘라내기 기능이다. > gdal_translate  -srcwin xoff yoff xsize ysize 원본.tif  복사본.tif          explain) Selects a subwindow from the source image for copying based on pixel/line location. 두 번째 방법으로 좌표점을 알 수 있는 tif 파일에서 절대 좌표를 기준으로 잘라내기 기능이다. > gdal_translate  -projwin ulx uly lrx lry inraster.tif outraster.tif           explain) Selects a subwindow from the source image for copying (like  -srcwin ) but with the corners given in georeferenced coordinates (by default expressed in the SRS of the dataset. Can be changed with -projwin_srs).

문제발생 1.Centos에서 gdal을설치하면 yum으로 지원하지 않아서 rpm으로설치해야한다. 2.이때 발생하는 에러중 하나로 gdal_translate 명령어중 HDF5 파일을 읽어들이는 vrt 옵션이 있는데 여기서 os에 hdf lib 가 설치되어 있지 않다면 파일을 읽어오지 못하는 에러가 발생한다 해결방안 hdf5 Source 설치파일을 다운받는다. 설치명령으로 설치한다 *hdf5 Release로 받으면 설치명령을 실행 할 수 없다. 다시 gdal install 과정을 진행한다 다만 아래내용으로 진행한다 >./configure --with-python --with-hdf5=(hdf5 lib path)

----------우분투에서 gdal을 사용하기위한 두 가지 방법-------------- 1. 커멘드라인으로 실행하기 $sudo apt-get install gdal-bin 2. 우분투 python-gdal을 사용할 때 $sudo apt-get install python-gdal 우분투에서 python-gdal을설치할 때 의존성문제가 발생한다면 아래와같은 명령어로 실행 한다 . ->sudo apt-get -f install 그 다음 다시 설치한다. ----------python에서 gdal 사용하기----------------- Ctrl+Alt+T (Open terminal) cd ~/folder/for/builds/from/source sudo apt-get install build-essential python-all-dev wget http://download.osgeo.org/gdal/1.11.0/gdal-1.11.0.tar.gz tar xvfz gdal-1.11.0.tar.gz cd gdal-1.11.0 ./configure --with-python make sudo make install

1. Longitude vrt 만들기  -->gdal_translate -of VRT HDF5:/home/lee/TEST/1B.GPM.GMI.TB2014.20151012-S000141-E000639.V03C.RT-H5://S1/Longitude /home/lee/TEST/long.vrt 2. Latitude vrt 만들기  -->gdal_translate -of VRT HDF5:/home/lee/TEST/1B.GPM.GMI.TB2014.20151012-S000141-E000639.V03C.RT-H5://S1/Latitude /home/lee/TEST/lat.vrt 3. Img vrt 만들기 gdal_translate -of VRT HDF5:/home/lee/TEST/1B.GPM.GMI.TB2014.20151012-S000141-E000639.V03C.RT-H5://S1/Tb /home/lee/TEST/img.vrt 4.  img.vrt 파일을 열어서    <metadata domain="GEOLOCATION">        <mdi key="X_DATASET">lon.vrt</mdi>        <mdi key="X_BAND">1</mdi>        <mdi key="Y_DATASET">lat.vrt</mdi>        <mdi key="Y_BAND">1</mdi>        <mdi key="PIXEL_OFFSET">0</mdi>        <mdi key="LINE_OFFSET">0</mdi>     ...

1. 원하는 버전의 gdal을 다운로드 받는다.    -->https://trac.osgeo.org/gdal/wiki/DownloadSource 2. make 하기   -->./configure --with-python(이렇게 하면 python에서 lib를 사용할 수 있다.)   -->make   -->sudo make install   -->sudo ldconfig     *   ImportError: libgdal.so.1: cannot open shared object file: No such file or directory   --> /usr/local/lib/libgdal.so.1 을 /usr/lib로 옮겨준다.   --> sudo ldconfig 실행 3. HDF lib는 따로설치해야한다(주의)