대용량 정사영상의 projection 방법?

[오동하]

국토지리정보원에서 정사영상을 다운 받았습니다.

대략 20기가(300여개의 파일) 정도되는 분량인데

문제는 정사영상에 좌표가 없습니다.

따라서 정사영상에 좌표를 주어  QGIS에 올리고 싶습니다. 궁극적으로는 하나의 큰 파일로 만들고 싶습니다.

정사영상은 1:5,000 도엽으로 나누어져 제공되고 있습니다.

이 작업을 손쉽게 할 수 있는 방법이 있는지요?

[Jonghun Jin]

안녕하세요.

DEM tiff 파일에 좌표가 잘못들어가 있어서  gdal_translate로 수정했던 기억이 있습니다.

( http://www.gdal.org/gdal_translate.html )

gdal 유틸리티 중에서 gdal_merge.py를 이용하면 파일 하나로 합치는 것도 (각 도엽에 메타데이터가 제대로 들어가 있다는 가정하에) 쉽게 가능합니다.

몇 년 전 기억이라… 더 좋은 방법이 있을것 같기도 합니다.

진종훈 드림

2016. 3. 16., 오전 10:23, 오동하 <dong...@gmail.com> 작성:

--
이 메일은 Google 그룹스 'OSGeo(Open Source GeoSpatial) 한국어 지부 - OSGeo Korean Chapter' 그룹에 가입한 분들에게 전송되는 메시지입니다.
이 그룹에서 탈퇴하고 더 이상 이메일을 받지 않으려면 osgeo-kr+u...@googlegroups.com에 이메일을 보내세요.
이 그룹에 게시하려면 osge...@googlegroups.com에 이메일을 보내세요.
https://groups.google.com/group/osgeo-kr에서 이 그룹을 방문하세요.
더 많은 옵션을 보려면 https://groups.google.com/d/optout을(를) 방문하세요.

[Sanghee Shin]

안녕하세요? 박사님, 

다른 사람들도 같은 어려움을 겪지 않으려면, 해당 문제점을 그대로 지리원에 알려 주는 게 좋을 듯합니다. 

감사합니다. 

신상희
---
Shin, Sanghee
Gaia3D, Inc. - The GeoSpatial Company
http://www.gaia3d.com 

2016. 3. 16., 오전 10:23, 오동하 <dong...@gmail.com> 작성:

[오동하]

예, 일단 국립지리정보원 질문에 올려놓았습니다. 아마 뭐라 답변이 있을 것 같습니다.

2016년 3월 16일 (수) 오후 12:53, Sanghee Shin <shs...@gaia3d.com>님이 작성:

--

오 동하

부산발전연구원 연구위원

해양환경연구실

mobile: 010-2486-4626

office: 051-860-8764

[오동하]

좌표가 있다면 가능한데 좌표가 아예 없습니다. 

해당 도곽정보를 읽어와서 projection을 해야 할 것 같습니다..

그런데 영상 매수가 많다보니 이 작업을 좀 손쉽게 할 방법을 찾고 있습니다.

2016년 3월 16일 (수) 오전 11:45, Jonghun Jin <jinjo...@gmail.com>님이 작성:

[Seongkyu Lee]

박사님, 

두 가지 정도 방법이 있을 것 같습니다.

1) listgeo와 geotifcp를 이용한 방법

> 좌표가 입력된 정사영상을 한 장 만듭니다

> listgeo xxxx.tif 명령어로 등록된 좌표 정보를 text 파일로 저장

> text 파일에서 'Corner Coordinates' 항목의 Upper/Lower Left, Upper/Lower Right 의 좌표값 변경

> geotifcp -g xxx.txt input.tif output.tif 명령어를 이용하여 좌표 정보를 등록

참고: http://egis3.lacounty.gov/dataportal/wp-content/uploads/2012/05/geotiff_info.pdf

2) gdal_translate를 이용한 방법

참고: http://gis.stackexchange.com/questions/38987/how-can-i-add-all-the-4-corners-as-gcp-in-an-image-that-i-want-to-reproject

프로그램 모두 QGIS에 포함되어 있습니다

파이썬 등으로 위 작업을 자동으로 수행하도록 하면 원하시는 작업이 가능할 것 같습니다.

Yours sincerely,

                                                                                                                 

Seongkyu Lee, Ph.D. in Spatial Information Engineering

Research Fellow,

Climate Informatics and Application Team,

APEC Climate Center

12 Centum 7-ro, Haeundae-gu,

Busan 612-020, Republic of Korea

Publications: http://scholar.google.co.kr/citations?hl=ko&user=XvFALoIAAAAJ

2016년 3월 16일 오후 2:00, 오동하 <dong...@gmail.com>님이 작성:

[오동하]

이 박사님 

알려준 방법 감사합니다. 

어째든 300장 영상의 도엽외곽 좌표는 필요하네요.

그런후 배치파일이나 파이썬 등으로 작업을 해야 할 것 같네요..

안되면 포토샵 같은데서 영상을 붙인 후에 한번에 projection을 하면 될 것 같기도 합니다.

이방법은 많은 지역적 에러를 낼 것 같습니다.

2016년 3월 16일 (수) 오후 3:29, Seongkyu Lee <geosl...@gmail.com>님이 작성:

[Sanghee Shin]

안녕하세요?

이승규 박사님이 첫번째 제안하신 방법이 아마 TFW(Tiff World File) 만드는 방법일 겁니다. 근데, 저렇게 복잡하게 하지 않고 그냥 6줄짜리 같은 이름의 .tfw만 추가로 노트패드 같은 걸로 만들어 주면 됩니다. 도엽명을 파싱해서 자동으로 .tfw를 만들어주는 배치스크립트 같은 걸 짤 수도 있겠지만 이게 약간 문제를 야기할 수도 있습니다. 아시다시피 도곽은 직사각형이 아니라 약간 한쪽이 더 기울어진 사다리꼴입니다. 따라서 도곽의 시작점과 정사각형 래스터 픽셀의 시작 위치가 정확히 맞지 않을 확률이 매우 높습니다. 

가장 좋은 방법은 지리원에 tfw 파일이나 혹은 GeoTiff로 달라고 하시는 방법입니다. 이후에는 gdal_merge.py 같은 걸로 합치시면 됩니다. 

참고하시기 바랍니다. 

신상희
---
Shin, Sanghee
Gaia3D, Inc. - The GeoSpatial Company
http://www.gaia3d.com 

2016. 3. 16., 오후 4:34, 오동하 <dong...@gmail.com> 작성:

[Seongkyu Lee]

안녕하세요?

신상희 대표님 의견에 약간 덧붙이면,

정사영상은 도엽 간 경계가 약간 중첩되어 있었습니다. 

다른 경우이기는 하지만 MODIS 위성영상을 gdal_merge로 작업해보면 잘 안되는 경우도 있고, 특히 파일이 많은 경우에는 먼저 vrt (virtual raster) 파일을 만들고 merge  작업을 하는 게 편할 수 있습니다.

일단 tfw나 geotiff 파일부터 받고 나서ㅎㅎ

내 삼성 디바이스에서 보냈습니다

-------- 원본 메시지 --------
발신: Sanghee Shin <shs...@gaia3d.com>
날짜: 2016/03/16 17:58 (GMT+09:00)
수신: osge...@googlegroups.com
참조: Sanghee Shin <shs...@gaia3d.com>
제목: Re: [osgeo-kr:2599] 대용량 정사영상의 projection 방법?

[서상일]

안녕하세요..

 

모니터링 회원입니다..^^;

 

국토지리정보원에서 정사영상을 다운 받으셨으면..51cm급 정사영상인것으로 사료됩니다..

 

잘 아시다시피 정사영상은 Tiff 포맷으로 되어있으며..Tiff 포맷의 좌표정보는 2가지 방법으로 입력할 수 있습니다..

 

첫번째는 Tif 영상 헤더 부분에 입력하는 방법, 두번째는 World 파일을 생성하는 방법입니다..

 

제가 알기론 지리원의 51cm 정사영상의 Tiff 파일에는 헤더에 좌표를 입력하지 않고 World파일만 생성하는것으로 알고 있습니다..

 

정확한 내용은 지리원에 확인하셔야 겠지만.. 영상만 다운 받으셨으면..당연히 좌표정보는 없을 것입니다.

 

영상의 좌표를 획득하시려면..원에서 World 파일을 받으셔야 합니다.

 

그리고 1/5,000 정사영상 300도엽은 통판재작이 가능하지만..파일포맷을 IMG로 생성하셔야 합니다..

 

작업시간은 고성능 WS 경우 6시간 정도??, 용량은 엄청 큽니다..QGIS 열수 있을지..도 의문입니다..

 

이상 짧은 지식으로 말씀드렸습니다..

--