본문 바로가기

GDAL : ogr2ogr : GeoJSON to PostgreSQL/PostGIS GeoJSON 파일을 PostgreSQL / PostGIS에 하나의 테이블로 입력할 때 사용할 수 있는 툴은 ogr2ogr 이라는 툴입니다. ogr2ogr의 기능은 꽤 다양한데요. 자세한 내용은 아래 웹 페이지를 참고하세요. ogr2ogr description web page : http://www.gdal.org/ogr2ogr.html 기능 중 하나가 GeoJSON과 같은 포맷을 공간정보 데이터베이스에 입력하는 것입니다. 예로 aaa.geojson 이라는 파일이 존재할 때 PostgreSQL / PostGIS에 입력하는 커맨드는 다음과 같습니다. ogr2ogr -s_srs EPSG:3857 -t_srs EPSG:3857 -f PostgreSQL PG:"dbname=mydb user=postgres" a.. 더보기
GDAL build with PostgreSQL GDAL을 빌드해서 설치할 때.. 향후 PostgreSQL과 연계되어 작동시키고 싶을 때는 아래와 같이 몇 가지 과정이 필요합니다. 1. PostgreSQL의 dev 패키지를 설치해줘야 합니다. PostgreSQL의 9.4 버전을 기준으로 아래와 같은 방식으로 패키지를 설치합니다. apt-get install -y postgresql-server-dev-9.4 2. gdal configure를 수행할 때 --with-pg 옵션을 통해 pg_config 파일의 위치를 잡아줍니다. ./configure --with-python --with-pg=/usr/bin/pg_config 위와 같이 PostgreSQL support yes가 나오면 성공! ./configure --with-python --with-pg.. 더보기
GDAL : ogrinfo example : so, where option GDAL : ogrinfo example : al, geom option - http://ngee.tistory.com/644 GDAL : ogrinfo example : so, where option - http://ngee.tistory.com/645 앞선 포스팅에 이어서~~ 계혹 ogrinfo 옵션을 몇 개만 더 정리하려고 합니다. 4. ogrinfo -al -so datasource_name 이번에는 al 옵션과 so 옵션을 같이 사용했습니다. al 옵션은 전체 레이어에 대한 정보(geometry 포함)를 출력하는 옵션이구요. so 옵션은 요약하는 옵션입니다. 이 2개를 같이 사용하는 경우 아래와 같은 결과를 볼 수 있습니다. command : ogrinfo -al -so DoDreamWay05_.. 더보기
GDAL : ogrinfo example : al, geom option command : ogrinfo (enter) Usage: ogrinfo [--help-general] [-ro] [-q] [-where restricted_where] [-spat xmin ymin xmax ymax] [-geomfield field] [-fid fid] [-sql statement] [-dialect sql_dialect] [-al] [-so] [-fields={YES/NO}] [-geom={YES/NO/SUMMARY}][--formats] datasource_name [layer [layer ...]] 이번 포스팅에서는 ogrinfo를 통해서 커맨드라인 명령어를 통해 shape 파일의 속성정보를 확인해보는 예제를 해보려고 합니다.shape를 하나 골라야죠. 저는 서울시 열린데이터 광.. 더보기
CRAN R : rgdal package installation in rstudio 그니깐요. rstudio에서 rgdal package를 설치하려고 할 때 컴퓨터에 gdal이랑 proj가 일단 설치가 되어 있어야 합니다. 엉엉 ㅠ gdal 설치하는 법이랑 proj 설치하는 방법에 대해서는 요 포스팅에서는 다루지 않고요. gdal과 proj를 기본 설정으로 소스코드에서 빌드해서 설치한 후 rstudio에서 다시 rgdal package를 설치하려고 하면 또 에러를 출력하는데요. 혹시 아래와 같은 경우라면 참고하시면 도움이 될 것 같아요. 1. gdal을 못찾네요. 라는 식의 에러 vi /etc/ld.so.conf 에 gdal lib 파일의 위치하는 폴더를 추가해주세요. gdal을 빌드해서 설치하는 경우 prefix 옵션을 주지 않으면 /usr/local/lib에 gdal lib 파일들이.. 더보기
GDAL : gdalwarp example : tr (resolution) command : gdalwarp (enter) Usage: gdalwarp [--help-general] [--formats] [-s_srs srs_def] [-t_srs srs_def] [-to "NAME=VALUE"] [-order n | -tps | -rpc | -geoloc] [-et err_threshold] [-refine_gcps tolerance [minimum_gcps]] [-te xmin ymin xmax ymax] [-tr xres yres] [-tap] [-ts width height] [-wo "NAME=VALUE"] [-ot Byte/Int16/...] [-wt Byte/Int16] [-srcnodata "value [value...]"] [-dstnodata "value [v.. 더보기
GDAL : gdal_translate example : projwin command : gdal_translate (enter) Usage: gdal_translate [--help-general] [--long-usage] [-ot {Byte/Int16/UInt16/UInt32/Int32/Float32/Float64/ CInt16/CInt32/CFloat32/CFloat64}] [-strict] [-of format] [-b band] [-mask band] [-expand {gray|rgb|rgba}] [-outsize xsize[%] ysize[%]] [-unscale] [-scale[_bn] [src_min src_max [dst_min dst_max]]]* [-exponent[_bn] exp_val]* [-srcwin xoff yoff xsize ysize] [-.. 더보기
GDAL : gdal_translate example : output size command : gdal_translate (enter) Usage: gdal_translate [--help-general] [--long-usage] [-ot {Byte/Int16/UInt16/UInt32/Int32/Float32/Float64/ CInt16/CInt32/CFloat32/CFloat64}] [-strict] [-of format] [-b band] [-mask band] [-expand {gray|rgb|rgba}] [-outsize xsize[%] ysize[%]] [-unscale] [-scale[_bn] [src_min src_max [dst_min dst_max]]]* [-exponent[_bn] exp_val]* [-srcwin xoff yoff xsize ysize] [-.. 더보기
GDAL : 래스터 데이터 좌표계 변경 (gdalwarp 사용) 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.2.. 더보기
GDAL : relief image creation using gdaldem GDAL의 커맨드라인 어플리케이션 중 gdaldem을 활용하면 relief image를 만들어 낼 수 있습니다. 먼저 gdaldem에 대한 자세한 설명은 아래 웹 페이지에 나와있습니다. http://www.gdal.org/gdaldem.html 간단하게 아래 파일을 이용해서 relief image를 만들어보는 과정을 정리해보려고 합니다. 먼저 이미지는 쉽게 다운로드가 가능한 SRTM을 사용했습니다. relief image를 만들기 위해서는 color 정보가 들어가 있는 txt 파일이 필요한데요. 아래는 color.txt 파일 예시입니다. 120 255 255 1 81 255 255 1 80 202 227 153 31 124 187 1 30 153 198 232 0 1 114 198 이제 gdaldem을 .. 더보기
GDAL : gdal_grid example GDAL에는 여러 CLI가 있는데요. 이번 포스팅에선 gdal_grid를 활용해서 interpolation하는 간단한 예제를 만들어보려고 합니다. gdal_grid에 대한 전반적인 설명은 아래 url에 나와있습니다. gdal_grid description : http://www.gdal.org/gdal_grid.html 여러 옵션도 존재하는데요. 이번 포스팅에서는 최대한 간단하게 interpolation하는 것을 목표로합니다. interpolation의 전반적인 설명은 아래 url에서 확인하시죠. ArcGIS Help 10.1 - Understanding interpolation analysis : http://resources.arcgis.com/en/help/main/10.1/index.html#//.. 더보기
GDAL : gdalwarp (using mask) gdalwarp할 때 mask를 씌울 수 있는데요. 간단하게 살펴보면 다음과 같습니다. gdalwarp -cutline /path/shp/file -srcnodata value1 -dstnodata -9999 before.tif after.tif 옵션을 잘 살펴보시면 됩니다. -cutline에 마스크로 사용할 shp 파일을 위치까지 잘 작성해주시면 됩니다. -srcnodata value1은 before.tif의 nodata 값을 의미하는 것이고 -dstnodata -9999는 after.tif에 nodata 값을 -9999로 설정하겠다는 것입니다. 다음 src_data dst_data를 작성해주시면 됩니다. 여기서 src_data는 before.tif로 설정한 것이고 dst_data는 after.tif로.. 더보기
GDAL : netCDF to GeoTiff GDAL을 활용해서 netCDF를 GeoTiff 형태로 변환하기 위해서는 다음과 같은 과정이 필요합니다. 1. GDAL 빌드 체크 일단 GDAL을 소스코드를 통해 빌드하신 분들은 netCDF 라이브러리를 설치하고 GDAL을 빌드했는지 부터 체크하셔요. gdalinfo 명령어로 netCDF 파일을 인자로 줬을 때 아래와 같이 에러가 나오면 GDAL 빌드시에 netCDF 라이브러리가 설치되어 있지 않았던 것이죠. aaa@hhh:~# gdalinfo r50.nc ERROR 4: `r50.nc' not recognised as a supported file format. gdalinfo failed - unable to open 'r50.nc'. 저는 dpkg -l | grep netcdf로 확인해보니깐 다음과.. 더보기
GDAL : esr -> geotiff .... using gdal_translate esr -> geotiff gdal_translate -of GTiff esrFile.esr geoTiffFile.tif option : -ot type (ex. -ot Float32) 더보기
GDAL : 리눅스에서 GDAL 소스코드 빌드 이번 포스팅에서는 GDAL 라이브러리를 설치할 때 소스코드로 빌드해서 설치하는 방법에 대해서 간략하게 작성하려고 합니다. 1. GDAL 소스코드 다운 받기 현재 GDAL의 최신 버전은 1.11-0 입니다. - Release 정보를 확인하시려면 다음 페이지에서 확인해주세요. (http://trac.osgeo.org/gdal/) 소스코드를 다운 받으시려면 GDAL 홈페이지에 가셔서 Dolwnload 링크를 클릭하시면 됩니다. GDAL Homepage : http://www.gdal.org/ 다운로드 하시면 아마도 gdal-version.tar.gz 으로 파일 이름이 되어 있죠. version은 여러분들이 다운로드 한 gdal 버전일 것이구요. 그럼 압축을 풀어주시면 됩니다. tar -xzvf gdal-ver.. 더보기
리눅스에서 이미지 정보 알아내기 리눅스에서, 이미지 파일에 대한 정보를 커맨드라인으로 알아내기 위한 몇 가지 간단한 방법을 소개하려구요. 이미지 이름은 white.png를 사용합니다. 1. 진짜 간단하게 file white.png 이라고 커맨드라인을 입력하시면 다음과 같이 이미지에 대한 간략한 정보가 프린트되요. ....$ file white.png white.png: PNG image data, 256 x 256, 8-bit/color RGB, non-interlaced 2. gdalinfo 이용 gdal은 라이브러리인데요. 다양한 이미지 관련 포멧의 입출력을 지원해주는 라이브러리입니다. 사실 이미지라고 하기 보다는 공간정보를 담고 있는 데이터에 대한 라이브러리가 더 알맞은 설명이겠네요. 아무튼 gdal을 설치한 경우 gdalinfo.. 더보기