Re: [osgeo-kr:1222] Qgis 사용 관련 질문 드립니다.

1,162 views
Skip to first unread message
Message has been deleted

Sanghee Shin

unread,
Jun 24, 2013, 3:02:16 AM6/24/13
to osge...@googlegroups.com, Sanghee Shin
안녕하세요?

1번 문제는 저장하실 때 encoding이 system이나 euc-kr로 되어 있는지 확인해 보시기 바랍니다.

2번 문제는 해당하는 레이어를 더블클릭하신 뒤 Style에서 Graduated로 분류하신 뒤 각각의 심볼을 더블클릭하셔서 사이즈를 조절하는 방법과, rule-based를 이용해서 한 방에 처리하는 방법이 있을 것 같습니다. 속성필드에 테이블의 필드가 나타나지 않는 것은 분류를 하기 위한 방법으로 graduated와 같은 방법을 쓰신 게 아닌가 싶은데요.. 이런 경우 필드가 숫자로 이뤄져 있어야 합니다. 그렇지 않은 경우에 그와 같은 상황을 만나실 수 있구요.. (말로 설명하려니 잘 전달이 안되네요.. 여하간..)

3번 문제는 heatmap 돌리면 아마 래스터로 나올 텐데요.. 그 레이어를 역시 더블클릭하셔서 스타일의 Render Type에서 Single band pseudocolor 같은 것으로 바꾸신 뒤 색 램프를 먹이시면 바뀌지 않을까 싶습니다.

참고하시기를..

신상희
---
Shin, Sanghee
Gaia3D, Inc. - Helping Earth Sustainable
http://www.gaia3d.com
Tel : +82-(0)2-3397-3475
Fax : +82-(0)2-3397-3478

2013. 6. 24., 오후 2:26, hkha...@gmail.com 작성:

>
> 주로 Arcgis를 사용하다가 최근에 QuantumGis를 새로 사용하기 시작했습니다. 아주 기본적인 기능만 사용하는 편인데,
> 몇가지 문제가 있어서 질문 드립니다.
>
>
> 1. 전국 행정경계화일에 특정 지점들을 지오코딩해 매핑해 중첩시킨 뒤, 서울 지역에 해당하는 부분만 속성 선택해서 별도 shp 화일로
> 저장했습니다.
> 그랬더니 '다른 이름으로 저장하기'한 속성테이블의 한글 필드들이 모두 Null로 표시가 되어 있네요.
> 다른 화일로도 비슷한 작업을 했더니 일부 지역만 지도가 추출돼 저장이 되지만 속성은 날아가버리는 일이 반복되네요. 해법이
> 있을까요?
>
> 2. 서울에서 석면이 함유된 공공건물 위치를 매핑한 뒤에 석면 함유량에 따라 점의 사이즈를 분류해 표시하려고 하는데,
> Arcgis와 달리 점 사이즈를 classification하는 버튼이 안보이네요. 게다가 색상 분류를 하려고 했더니
> classify하는 기준필드 메뉴에 속성테이블의 필드가 하나도 잡히지 않네요.
> 왜 그럴까요?
>
> 3 Heatmap을 그린 뒤에 그려진 등고선의 색상 스케일을 내가 원하는 방식으로 바꿀 수 있는 방법이 있나요?
>
> 참고로 한글 깨지는 문제 때문에 Qgis 1.8 대신 러시아에서 만든 NexusGis를 사용하고 있습니다. 제가 언급한 문제
> 1.같은 경우 러시아판의 버그는 아닌지 모르겠습니다.
> 혹시 저같은 문제 겪으신 분은 없는지 모르겠군요. 감사합니다.
>
> -
>
>
> --
> Google 그룹스 'OSGeo(Open Source GeoSpatial) 한국어 지부 - OSGeo Korean Chapter' 그룹에 가입했으므로 본 메일이 전송되었습니다.
> 이 그룹에서 탈퇴하고 더 이상 이메일을 받지 않으려면 osgeo-kr+u...@googlegroups.com에 이메일을 보내세요.
> 이 그룹에 게시하려면 osge...@googlegroups.com(으)로 이메일을 보내세요.
> http://groups.google.com/group/osgeo-kr에서 그룹을 방문하세요.
> 더 많은 옵션을 보려면 https://groups.google.com/groups/opt_out을(를) 방문하세요.
>
>

mapmatters

unread,
Jun 24, 2013, 8:20:50 AM6/24/13
to osge...@googlegroups.com
안녕하세요, 제가 아는 부분에 대해서만 답변을 드릴까합니다

우선은 1.8버전보다는 곧 나올 2.0 또는 nightly버전인 1.9 버전 사용하시는걸 추천드립니다

개인적으로 1.9를 쓰기 시작한 이후로 Arc를 거의 쓰지 않을 정도로 사용환경이 쾌적해진것 같습니다..

신상희대표님이 올려주신 포스트를 보시면 2.0 nightly 버전 설치하실 수 있습니다



2. size classification

classify 버튼에 있지 않고 고급옵션에 있습니다

필드에 있는 값대로 버블의 크기를 맞춰주게 됩니다. 값이 3000이면 사이즈가 3000이 됩니다..

따라서 마커의 크기는 1 또는 2가 아니라 0.0001 이런식으로 조절해줘야합니다

이미지 파일 한장과, 전에 다른 수업때 만든자료를 첨부해드립니다 참고해주세요



3. heatmap

플러그인 중 sextante를 설치하시고 kernel density 를 사용하시는 것을 권해드립니다

heatmap과 kernel density는 원리가 다르기도 하지만, kernel이 범용적이고 QGIS에서 돌아가는데 훨씬 빠른것 같더군요

주의하실점은 레이어이름, 저장장소 등 모든 것이 영어로 되어있어야한다는 점입니다.

(QGIS를 사용하다보면 자연스레 모든 파일은 영어로 쓰게 되긴 합니다만 ^^)

sextante - kernel density 사용방법은 조만간 자료를 만들 예정인데 공유할수있도록 하겠습니다




2013년 6월 24일 월요일 오후 2시 26분 31초 UTC+9, hkha...@gmail.com 님의 말:
Lecture_tourmap2.pdf
스냅1.PNG
Message has been deleted

hkha...@gmail.com

unread,
Jun 29, 2013, 5:06:20 AM6/29/13
to osge...@googlegroups.com, Sanghee Shin


2013년 6월 24일 월요일 오후 4시 2분 16초 UTC+9, Shin, Sanghee 님의 말:
Message has been deleted

hkha...@gmail.com

unread,
Jun 29, 2013, 5:33:29 AM6/29/13
to osge...@googlegroups.com
  주로 Arcgis를 사용하다가 최근에 QuantumGIS를 새로 사용하기 시작했습니다. 아주 기본적인 기능만 사용하는 편인데, 몇가지 문제가 있어서 질문 드립니다.
  

1. 전국 행정경계화일에 특정 지점들을 지오코딩해 매핑해 중첩시킨 뒤, 서울 지역에 해당하는 부분만 속성 선택해서 별도 shp 화일로 저장했습니다.
     그랬더니 '다른 이름으로 저장하기'한  속성테이블의 한글 필드들이 모두 Null로 표시가 되어 있네요.
     다른 화일로도 비슷한 작업을 했더니 일부 지역만 지도가 추출돼 저장이 되지만 속성은 날아가버리는 일이 반복되네요. 해법이 있을까요?

2.  서울에서 석면이 함유된 공공건물 위치를 매핑한 뒤에 석면 함유량에 따라 점의 사이즈를 분류해 표시하려고 하는데,
     Arcgis와 달리 점 사이즈를 classification하는 버튼이 안보이네요. 게다가 색상 분류를 하려고 했더니 classify하는 기준필드 메뉴에 속성테이블의 필드가 하나도 잡히지 않네요.
     왜 그럴까요?

3   Heatmap을 그린 뒤에 그려진 등고선의 색상 스케일을 내가 원하는 방식으로 바꿀 수 있는 방법이 있나요?

    참고로 한글 깨지는 문제 때문에 Qgis 1.8 대신 러시아에서 만든 NextGIS를 사용하고 있습니다. 제가 언급한 문제 1.같은 경우 러시아판의 버그는 아닌지 모르겠습니다.

hkha...@gmail.com

unread,
Jun 29, 2013, 5:35:40 AM6/29/13
to osge...@googlegroups.com

안녕하세요? 바로 위에 있는 지난 월요일에 올렸던 QuantumGIS 기본 사용법 질문에 대해 신상희 대표님과 mapmatters님 설명대로 해보니 대체로 해답을 찾았습니다. 고맙습니다.

1) 추출해 다시 저장한 지도의 속성이 null로 표시되는 문제는 Quantum GIS 2.0 Nightly version을 설치해 사용해보니 깨끗이 해결됐습니다.
2) 점 사이즈를 조정하는 것도 고급 옵션을 활용하니 간단히 해결되네요.속성 필드에 데이터 필드가 나타나지 않았던 것은 graduated 모드에서 해당 필드 data type이 string으로 읽혀져서 그랬습니다.
    Quantum GIS가 CSV 화일의 데이터를 무조건 string으로 읽곤 한다는 사실을 깜박했군요. csvt 화일로 data type을 지정해주니 해결됐습니다.
3) heatmap은 스타일의 Render Type에서 Single band pseudocolor 로 바꾸면 되는데, 색상 조정은 자유롭지는 않은 것 같더군요.
   sextetante로도 돌려보려고 하는데 나중에 kernel density map 만드는 방법 자료로 정리해주시면 참고해서 잘 활용하겠습니다. 
  

  추가로 궁금한 점은 특정 필드의 data type을 string에서 double이나 real로 바꿔주고 싶을 때 항상 csvt 화일로 설정해줘야 하나요?
  Quantum GIS 사용시 가장 번거로운 점이 excel 화일을 직접 못 불러오고, csv의 data type도 일일이 지정해줘야 한다는 건데 이를 우회할 수 있는 방법은 없나 궁금합니다.
  가령 mmqgis plugin의 text to float 기능을 사용해보니 data type 변환은 되는데  다른 필드의 한글 속성은 모두 null로 표시가 돼 버리네요.
 
  GIS 활용하는 전문가 분들이 주로 찾는 곳 같은데 기초적인 질문을 올려서 죄송한데, 거듭 고맙습니다.

  

Sanghee Shin

unread,
Jul 1, 2013, 9:58:59 AM7/1/13
to osge...@googlegroups.com, Sanghee Shin
안녕하세요?

CSV 파일을 읽어 오실 때 어떻게 읽어 오시는지요?

제가 Layer -> Add Delimited Text Layer 를 통해서 CSV를 읽어 보니 첨부와 그림처럼 Integer, Real, Texㅅ type으로 자동으로 알아서 읽어 오는데요?

그리고 Layer -> Add Vector Layer를 통해서 XLS 포맷은 읽어 오실 수 있을 겁니다. 다만, 이 경우에는 속성값만 읽어올 수 있어서 위치값을 등록하기는 어려울 것입니다. 

그리고 이미 String Type으로 되어 있는 자료는 Field calculator의 Conversion 함수를 이용하시면 Integer나 Real Type 등으로 혹은 그 반대로 변환하실 수 있습니다. 

참고하시기 바랍니다. 

신상희 



---
Shin, Sanghee
Gaia3D, Inc. - Helping Earth Sustainable 
http://www.gaia3d.com 
Tel : +82-(0)2-3397-3475
Fax : +82-(0)2-3397-3478

2013. 6. 29., 오후 6:35, hkha...@gmail.com 작성:

mapmatters

unread,
Jul 1, 2013, 11:53:21 PM7/1/13
to osge...@googlegroups.com
kernel density 자료 올려드립니다. 진작 만들었어야할 것을 이 기회에 만들게되었네요

이외에 구글그룹스 방문하시는 분들께 참고가 될까해서 pdf랑 예제에 사용한 shp파일도 함께 업로드합니다

그럼 수고하세요~


GU 박용재



2013년 6월 29일 토요일 오후 6시 35분 40초 UTC+9, hkha...@gmail.com 님의 말:
Lecture_kerneldensity.pdf
shpfiles.zip

hkha...@gmail.com

unread,
Jul 2, 2013, 4:48:46 AM7/2/13
to osge...@googlegroups.com, Sanghee Shin

고맙습니다.   신 대표님,  며칠 계속 질문만 올렸는데 그러고 보니 제 인사가 늦었네요. 죄송합니다. 저는 YTN에서 방송 기자로 일하는 함형건이라고 합니다.

1년여 전부터 독학해서 GIS를 방송보도에 접목해보려고 노력하고 있는데, Arcgis로는 다양한 공간 분석이 잘 됐었는데 QuantumGIS로는 error 나는 게 너무 많네요. 

신 대표님 말씀하신데로 다시 작업해봤지만 field data type 읽어내는 문제는 제 컴퓨터에서는 그대로네요.

서울시내 공공건축물의 석면분포정보로 mapping 했는데 관련 화일을 첨부했습니다.

다시 csvt 화일 없이 처음부터 불어들이기를 해봐도 원 화일에서는 필드가 number로 정의돼 있는데도, 건물면적 필드가 Quantum GIS에서는 text로 읽혀짐을 확인할 수 있었습니다.


한가지만 더 질문 드리면 이 석면함유 건축물 위치 point map을 서울시 행정경계에 위치를 기반으로 조이닝하려고 하는데요.
행정동 마다 해당 구역 안에 있는 석면건물 속성들은 자동적으로 행정동마다 편입되서 속성표에 별도의 필드로 합산되서 나타나게 하고 싶습니다.
Arcgis에서는 이런 작업이 쉽게 수행됐는데 QuantumGIS에서는 계속 에러가 나네요.
벡터 메뉴->데이터관리도구->위치에 따라 속성 결합으로 들어가서
대화창에서 교차하는 모든 피처 속성정보 이용의 '총계'를 체크한 뒤 실행했습니다.
'대상 벡터 레이어'와 '벡터레이어 조인'에는 석면 point map와 행정경계지도를 각각 입력하고 거꾸로 바꿔서도 입력해봤는데 결과적으로 똑같은 에러메시지가 나옵니다. 

다음과 같은 에러메시지가 뜨는데 무슨 얘기인지는 잘 모르겠습니다. 다른 분들은 QuantumGIS에서 spatial joining이 잘 되시나요?

-----------------------------------------------------------------------------------------
Python 프로그램을 실행하는 동안 오류가 발생했습니다:
Traceback (most recent call last):
  File "C:\OSGeo4W\apps\qgis-dev\python\plugins\fTools\tools\doSpatialJoin.py", line 99, in accept
    res = self.compute(inName, joinName, outPath, summary, sumList, keep, self.progressBar)
  File "C:\OSGeo4W\apps\qgis-dev\python\plugins\fTools\tools\doSpatialJoin.py", line 145, in compute
    fieldList1.extend(fieldList)
TypeError: QgsFields.extend(QgsFields): argument 1 has unexpected type 'list'

Python 버전:
2.7.4 (default, Apr  6 2013, 19:54:46) [MSC v.1500 32 bit (Intel)]


QGIS 버전:
1.9.0-Master Master, 619558d
----------------------------------------------------------------------------------------------

함 형 건 
YTN 보도국 앵커팀 차장. 밤 10시, 11시 '뉴스나이트' 앵커.





서울시공공건축물석면정보_지오코딩2.csv

Sanghee Shin

unread,
Jul 2, 2013, 5:52:01 AM7/2/13
to hkha...@gmail.com, Sanghee Shin, osge...@googlegroups.com
안녕하세요?

주신 자료 원본 엑셀 파일을 받아 볼 수 있을까요? 제 컴퓨터(맥이든 윈도우든)에서는 엑셀로 읽어보거나 혹은 에디터로 읽으나 혹은 QGIS로 읽으나 한글로된 필드명과 레코드가 깨져 보입니다.

한글 깨진 것과 관계없이는 Layer -> Add delimited text layer로 불러 들이니 첨부 파일처럼 수로 되어 있는 것들은 수로 읽어 오는 것 같습니다.

아, 그리고 이 파일의 좌표계는 어떻게 되는지요?

QGIS에서 Spatial Join 잘 되는데요.. 저도 한 번 테스트해 보겠습니다. 요즘 2.0으로 주요한 업그레이드 중이어서 과거 API와 잘 안맞는 경우들이 있습니다.

일단 지금은 오늘 번개 모임에 나가야 해서.. 총총....33

신상희
---
캡처2.PNG
캡처1.PNG

오동하

unread,
Jul 2, 2013, 9:52:02 PM7/2/13
to osge...@googlegroups.com
원본파일을 csv로 만들때 뭔가 문제가 있는 것 같습니다.



2013년 7월 2일 오후 6:52, Sanghee Shin <shs...@gaia3d.com>님의 말:
Shin, Sanghee
Gaia3D, Inc. - Helping Earth Sustainable
http://www.gaia3d.com
Tel : +82-(0)2-3397-3475
Fax : +82-(0)2-3397-3478

2013. 7. 2., 오후 5:48, hkha...@gmail.com 작성:
> <서울시공공건축물석면정보_지오코딩2.csv>



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






--
---------------------------------------------------------------------
오 동하                      Dongha Oh
부산발전연구원       Busan Development Institute
                                   Research Fellow
TEL  82-51-860-8764, 010-2486-4626
--------------------------------------------------------------------

hkha...@gmail.com

unread,
Jul 3, 2013, 1:59:41 AM7/3/13
to osge...@googlegroups.com, Sanghee Shin

신 대표님, 안녕하세요? 함 형건입니다. 이상하게 제 컴퓨터에서는 잘 열리는데,  어제 첨부한 화일은 업로드되면서 실제로 깨져 올라갔네요.왼쪽부터 10번,11번,12번째 필드가 건물 면적과 석면 면적,
그리고 건물 면적대 석면 면적 이렇게 3가지입니다. 모두 number여야 하는데, text, double,text 이렇게 뜨는거죠.  
찾아보니 geocoding한 x,y 좌표 field가 붙어있는 액셀 화일은 갖고 있는게 없고요. 어제 csv 화일이 버전이 2개로 있길래 2개 다 첨부해보고, 이것을 액셀화일로 변환해서 다시 첨부해봤습니다.
x,y 좌표가 없고 data type도 약간 다른 원래 석면 액셀 데이터 화일도 같이 올립니다. 다 깨져 나타나면 할 수 없고요.^^ 
  
제가 좌표계 관련해서도 궁금한 점이 좀 많은데요.
좌표계는 katec으로 통일해서 쓰고 있는데, 같이 쓰는 행정경계지도를 biz gis 사이트에서 카텍좌표로 다운로드 받은 화일을 쓰고 있어서 그렇습니다.
주말에 이민파 선생님도 TMS for Korea 관련 답변에서확인을 해주셨는데, 제가 쓰는 행정경계지도나 석면 지도나 katec이 맞다고 하시거든요.
그런데 이상한 것은 QuantumGIS로 불러들여 속성의 일반정보 탭을 열어보면 EPSG :4326 WGS84라고 항상 뜨네요. 
그래서 보통은 일반정보 탭에서 좌표계 지정을 매번 katec으로 다시 지정해서 작업하든지, 화면 우하단 눌러서 on the fly로 katec으로 해놓고 작업합니다. (호기심에 WGS84 on the fly로 해도 아무 지도가 나타나지 않는데 그건 왜 그런지도 잘 모르겠습니다.)  
좌표 지정하는 방법이 3가지 방법이 있는 것 같던데요. 일반정보탭에서 좌표 지정하는 것과 마우스 우클릭 레이어 CRS 설정 메뉴로 들어가서 katec 지정하는 것과 위의 벡터 메뉴의 데이터 관리도구의 현재 투영정의 메뉴로 좌표 지정하는 것들은 또 어떤 차이가 있는지도 궁금합니다. 현재 투영정의로는 이미 사용자 정의를 정확히 해놓았는데도 katec을 지정하려고 하니 작동이 안되네요.

마지막으로 또 하나는 궁금한 점은 지금까지 네이버 지도는 좌표가 katec이고 Korea 2000/ unified CS EPSG:5179라고 알고 있었는데, 좌표값이 약간 다르더군요.
주말에 이민파 선생님이 가르쳐주신 좌표값으로 katec을 사용자 정의하니 네이버 지도와 잘 맞았습니다.
다른 게시판에 보면 네이버 지도나 katec이 5179라고 써놓은 것을 본 기억이 있는데 이게 잘못된 정보인가요?

아무튼 바쁘실텐데 여러 모로 고맙습니다.~
 

함 형 건 
YTN 보도국 앵커팀 차장. 밤 10시, 11시 '뉴스나이트' 앵커.
개인 블로그 Mapping the News : http://blog.naver.com/barrious12/110168815479

서울시공공건축물석면정보_지오코딩2.csv
서울시공공건축물석면정보_지오코딩.csv
서울시공공건축물석면정보_지오코딩.xls
서울시공공건축물석면정보.xls
서울시동경계_201105.shp

오동하

unread,
Jul 3, 2013, 4:28:09 AM7/3/13
to osge...@googlegroups.com
좌표에 있어서는 초기 값이 WGS84로 QGIS설치시 세팅되어 있습니다. 메뉴중에서 SETTING-OPTION에 들어가시면 CRS를 바꾸시면 항상 KATEC으로 세팅할 수 있을 것 같습니다. 


2013년 7월 3일 오후 2:59, <hkha...@gmail.com>님의 말:

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

오동하

unread,
Jul 3, 2013, 4:32:38 AM7/3/13
to osge...@googlegroups.com
SHP파일은 동일이름의 확장자가 다른 여러 파일을 같이 올려야 합니다..아마 4-5개 동일이름의 확장자가 다른 파일이 있습니다.



2013년 7월 3일 오후 5:28, 오동하 <dong...@gmail.com>님의 말:
Message has been deleted

hkha...@gmail.com

unread,
Jul 5, 2013, 4:29:25 AM7/5/13
to osge...@googlegroups.com

예, 고맙습니다. 좌표는 초기갑에서 setting-option에서 바꿔주니 katec으로 표시돼 나오는군요.
 
이전 석면 화일 관련 data type 읽기 질문 드리면서 참고로 첨부했던 서울시 행정경계화일은 다른 관련 화일도 전부 다시 첨부해 봤습니다.
그런데 prj. 화일은 원래 없네요.  

spatial join이나 data type 읽는 문제는 여전히 제 컴퓨터에서는  error가 나는 것 같은데, 아마도 제 화일에 어떤 문제가 있는지도 모르겠습니다.
그냥 Quantum GIS의 기본 기능이 Arcgis처럼 문제없니 돌아가는지 살펴보면서 하나하나 배워가려고 하는 과정인데 시간을 두고 무엇이 문제인지 생각해보도록 하겠습니다.      
행정동_201105.sbn
행정동_201105.sbx
행정동_201105.shp
행정동_201105.shx
행정동_201105.dbf

hkha...@gmail.com

unread,
Jul 5, 2013, 4:29:57 AM7/5/13
to osge...@googlegroups.com

예, 고맙습니다. 좌표는 초기갑에서 setting-option에서 바꿔주니 katec으로 표시돼 나오는군요.
 
이전 석면 화일 관련 data type 읽기 질문 드리면서 참고로 첨부했던 서울시 행정경계화일은 다른 관련 화일도 전부 다시 첨부해 봤습니다.
그런데 prj. 화일은 원래 없네요.  

spatial join이나 data type 읽는 문제는 여전히 제 컴퓨터에서는  error가 나는 것 같은데, 아마도 제 화일에 어떤 문제가 있는지도 모르겠습니다.
그냥 Quantum GIS의 기본 기능이 Arcgis처럼 문제없니 돌아가는지 살펴보면서 하나하나 배워가려고 하는 과정인데 시간을 두고 무엇이 문제인지 생각해보도록 하겠습니다.      
행정동_201105.sbn
행정동_201105.sbx
행정동_201105.shp
행정동_201105.shx
행정동_201105.dbf

Sanghee Shin

unread,
Jul 5, 2013, 12:00:41 PM7/5/13
to hkha...@gmail.com, Sanghee Shin, osge...@googlegroups.com
안녕하세요?

1. 우선 주신 파일들의 좌표계가 katec 좌표계를 쓰신다고 하셔서 QGIS의 Settings -> Custom CRS을 이용하여 Katec 좌표계를 추가하였습니다. 이런 작업을 수행한 이유는 EPSG(European Petroleum Survey Group)라고 전 세계의 여러 좌표계를 모아서 DB화 하는 곳에 해당 좌표계가 표준 좌표계로 등록되어 있지 않기 때문입니다. Custom CRS를 등록하는 방법은 첨부된 이미지를 참조하십시오. 

asbestos.zip
Custom_CRS.PNG
Overlay2.PNG
FieldType.PNG
건물면적단위오류.PNG
석면면적단위오류.PNG
overlay.PNG

오동하

unread,
Jul 5, 2013, 9:25:49 PM7/5/13
to osge...@googlegroups.com
답변이 깔끔하네요. 혹시 spatial join의 에러가 좌표계 때문에 생기는 것은 아닌지 모르겠습니다.



2013년 7월 6일 오전 1:00, Sanghee Shin <shs...@gaia3d.com>님의 말:
안녕하세요?

1. 우선 주신 파일들의 좌표계가 katec 좌표계를 쓰신다고 하셔서 QGIS의 Settings -> Custom CRS을 이용하여 Katec 좌표계를 추가하였습니다. 이런 작업을 수행한 이유는 EPSG(European Petroleum Survey Group)라고 전 세계의 여러 좌표계를 모아서 DB화 하는 곳에 해당 좌표계가 표준 좌표계로 등록되어 있지 않기 때문입니다. Custom CRS를 등록하는 방법은 첨부된 이미지를 참조하십시오. 


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




2. 일단 주신 csv 파일에서 한글들이 모두 깨져 보여서, 주신 서울시공공건축물석면정보_지오코딩.xls 파일을 일단 csv로 다시 변환하였습니다. 

3. QGIS의 Layer -> Add Delimited Text Layer를 이용하여 변환된 csv 파일을 읽어 들였습니다. 그런데, 말씀하신대로 건물면적, 석면면적 등이 숫자로 들어오지 않고 모두 String 타입으로 들어왔습니다. 



4. 이상하게 생각되어 주신 엑셀 자료를 다시 살펴 보았습니다. 그랬더니 엑셀에 숫자로만 표시되어 있어야 하는 레코드 값에 숫자가 아닌 것들이 다수 포함되어 있더군요. 관련 그림 첨부하였습니다. 





5. 여하간 해당 자료를 불러와서 좌표계를 Custom CRS로 지정한 Katec로 맞춰 줬습니다. 그리고 말씀하신 대로 www.biz-gis.com 의 자료실에서 첨부한 그림과 같이 행정동 경계를 다운로드 받은 뒤 좌표계를 EPSG:5174 로 설정해 줬습니다. 그랬더니 아래와 같이 모든 자료가 일단은 올라오는 것으로 나타났습니다. 속성도 모두 한글과 함께 잘 물고 오구요. 




6. QGIS의 플러그인인 TMS for Korea를 이용해 네이버 하이브리드맵도 깔아보니 아래와 같이 잘 나옵니다. 
QGIS 최근 버젼은 좌표계가 다를 경우 알아서 자동으로 좌표계를 맞춰 줍니다. 




7. 말씀하신 Spatial Join을 해 보니 역시 에러를 뿜어내더군요. 이게 예전 버전에서는 잘 되었고, 올초에 또 업그레이드 중에 버그가 발견되어 수정되었다고 나왔는데 다시 버그가 생긴 것 같기도 하고, 아니면 저나 함 기자님 컴퓨터의 Python이나 QGIS 설정 문제일 수도 있습니다. (여러 버젼을 깔아서 쓰다가 꼬인..) 이 부분은 원인을 좀 찾아보고 리포팅하거나 해야될 것 같습니다. 

8. 좌표계를 항상 설정하기 귀찮으시면 원본 자료를 원하는 좌표계로 저장 한 번 하시면 됩니다. Layer -> Saver Layer를 통해서 원하는 좌표계 지정해 주시고,(이때 주의점은 원래 읽어들일 때 정확한 좌표계를 줘야 한다는 점입니다.) 저장해 주시면 다음부터는 알아서 좌표계 물고 올 것입니다. 

네이버 지도가 KATEC 계열의 좌표계라는 의미이지 과거 KOTI에서 지정한 KATEC 좌표계와 같다는 뜻은 아닙니다. 우리나라 좌표계에 대한 정보는 http://www.osgeo.kr/17 에서 확인하실 수 있습니다. 

덧붙여 제가 변환한 석면자료 SHP 자료를 보내 드립니다. 참고로 SHP 포맷은  파일이름.shp(도형정보), 파일이름.dbf(속성정보), 파일이름.shx(인덱스정보)의 최소 3개 파일이 함께 전달되어야 합니다. 최근에는 파일이름.prj(좌표계정보)를 함께 전달하고는 합니다. 

참고하시기 바라며, spatial join 관련 문제는 조금 더 살펴 보도록 하겠습니다. 어차피 저도 다다음주 월요일에 이화여대에서 QGIS 강의를 해야 하는데, 이 기능과 관련된 내용이 있어서 빨리 살펴보야겠네요. 

신상희 

 
---
Shin, Sanghee
Gaia3D, Inc. - Helping Earth Sustainable 
http://www.gaia3d.com 

2013. 7. 3., 오후 2:59, hkha...@gmail.com 작성:

<서울시공공건축물석면정보_지오코딩2.csv><서울시공공건축물석면정보_지오코딩.csv><서울시공공건축물석면정보_지오코딩.xls><서울시공공건축물석면정보.xls><서울시동경계_201105.shp>


Sanghee Shin

unread,
Jul 7, 2013, 7:45:04 PM7/7/13
to osge...@googlegroups.com, Sanghee Shin
안녕하세요?

이 문제는 최근 업그레이드 중에 다시 발견된 문제라고 하며, 7월 3일(?)인가 이후 버전에서는 해결이 되었다고 하네요. 그런데, 제가 테스트를 해 보니 어떤 경우는 되고 어떤 경우는 QGIS 자체를 죽이는 것 같아 다시 한 번 살펴보고 있습니다.

신상희
---
Shin, Sanghee
Gaia3D, Inc. - Helping Earth Sustainable
http://www.gaia3d.com
Tel : +82-(0)2-3397-3475
Fax : +82-(0)2-3397-3478

2013. 7. 6., 오전 10:25, 오동하 <dong...@gmail.com> 작성:
>> 개인 블로그 Mapping the News : http://blog.naver.com/**barrious12/110168815479<http://blog.naver.com/barrious12/110168815479>
>>
>> <서울시공공건축물석면정보_지오코딩2.csv><서울시공공건축물석면정보_지오코딩.csv><서울시공공건축물석면정보_지오코딩.xls>
>> <서울시공공건축물석면정보.xls><서울시동경계_201105.shp>
>>
>>
>>
>>
>
>
> --
> ---------------------------------------------------------------------
> 오 동하 Dongha Oh
> 부산발전연구원 Busan Development Institute
> Research Fellow
> TEL 82-51-860-8764, 010-2486-4626
> --------------------------------------------------------------------
>
Message has been deleted

hkha...@gmail.com

unread,
Jul 9, 2013, 12:19:35 AM7/9/13
to osge...@googlegroups.com, Sanghee Shin
 안녕하세요?  친절한 답변 감사드립니다.  말씀하신데로  2.0에서 공간 조이닝 할 때 error가 발생한다고  Minoru Akagi라는 분도 QuatumGIS Homepage에 bug 신고했더군요.
신 대표님 말씀데로 말끔히 고쳐지진 않을 듯 한데 어쨌든 나중에 nightly version 2.0 다시 설치해서 확인해 봐야겠네요. 
Arcgis가 만기가 되서 앞으로는 QuantumGIS로만 작업을 해볼까 했는데, spatial joining에서 걸리네요. spatial joining같은 핵심 기능은 언젠가는 제대로 고쳐졌으면 하고 기대해 봅니다.

 data type 읽는 문제는 제가 처음 연습삼아 올렸던 석면 화일에 dirty data가 많아서 그랬던 것 같습니다. 오동하 박사님이 지적하신데로 원본화일을 csv로 변환하면서 문제가 있었네요.
 일부 필드의 내용에 comma가 섞여 있어서 csv로 바꿀 때 셀 내용이 잘못 섞여서 변환된 것이 있었습니다. 건물 면적과 석면면적 필드에 일부 연도, 날짜 데이터 등이 섞여 들어간 셀이 있던데 이것들을 지우고 불러들이니 제대로 double로 다 인식하네요.

 시험삼아 여러 다른 데이터도 QuantumGIS에 불러들여 지도 그리기를 해보다가 느낀 건데, Arcgis처럼 스타일 탭에서 graduated mode의 데이터 구간별 분류를 할 때 histogram이 뜨면 좋을 텐데 하는 생각을
했습니다. 그래프를 보면서 수동으로 구간을 조정하는 기능이 있으면 참 유용할 텐데 말이죠. 1.8 version에서는 stats라는 plugin에서 vector file의 속성을 histogram으로 보여주는 기능이 있어서 이를 참고로
구간을 수동으로 다시 입력하곤 했는데, 2.0에서는 histogram 그리는 plugin이 아직 없더군요.

 현재 언론계에 gis user가 거의 없어서 대부분 혼자서 해결해야 할 부분이 많고 국내 지리관련 공공데이터가 적다보니 적극적으로 활용은 하지 못하고 있는데,  앞으로는 QuamtumGIS를 방송에 간간이 이용해 보려고 합니다.  아낌없이 조언을 주시는 신상희 대표님, 오동하 박사님, 박용재 연구원님 모두 거듭 고맙습니다.
(박용재 연구원님이 올려주신 kernel density algorothm 사용법도 아주 유용한 자료이네요. 고맙습니다. 다만 제가 조작을 잘못한 건지 경우에 따라서 생각데로 안그려질 때도 있던데 이 부분도 나중에 정리해서 질문드려보도록 하겠습니다.)


 -함 형 건- 
YTN 보도국 앵커팀 차장. 밤 10시, 11시 '뉴스나이트' 앵커.
개인 블로그 Mapping the News : http://blog.naver.com/barrious12/110168815479

Sanghee Shin

unread,
Jul 9, 2013, 11:00:32 AM7/9/13
to hkha...@gmail.com, Sanghee Shin, osge...@googlegroups.com
안녕하세요?

제가 테스트를 해 보니 최근 버전에서는 Spatial Join이 잘 되는 것 같습니다.

안되고 죽는 경우가 있어 왜 그런가 살펴보니 제가 만든 asbestos.shp 의 점 중 몇몇이 일종의 outlier로 존재해서 그런 것 같습니다. 어떤 경우는 도저히 말이 안되는 좌표로, 아마도 투영법 한계 밖으로 추정되는 곳에 표시가 되고 있더군요. 그래서 한반도를 벗어난 - 즉, 좌표가 정상적이지 않은 - 점들을 제거하고 가장 최신 버젼으로 Spatial Join을 돌리니 금방 돌아가고 값도 정상적으로 나오는 것 같습니다.

말씀하신 stats 플러그인 같은 경우 아직 버전 2.0의 API 변경에 따른 업그레이드가 이뤄지지 않은 것 같습니다. 다만, 함 기자님께서 말씀하신 내용은 새겨 들을 만한 것으로서 QGIS 커뮤니티에 의견 개진을 하도록 하겠습니다.

그리고 혹시나 QGIS를 이용하셔서 방송을 하시게 되시면 화면 캡쳐라도 하셔서 알려 주시면 감사하겠습니다. 한국의 뉴스 채널에서 QGIS를 사용했다고 QGIS 커뮤니티에 알려 주면 애들 아마도 난리날 것 같아서요. ㅎ

신상희
---
Shin, Sanghee
Gaia3D, Inc. - Helping Earth Sustainable
http://www.gaia3d.com
Tel : +82-(0)2-3397-3475
Fax : +82-(0)2-3397-3478

2013. 7. 8., 오후 5:37, hkha...@gmail.com 작성:

>
> 안녕하세요? 친절한 답변 감사드립니다. 말씀하신데로 2.0에서 공간 조이닝 할 때 error가 발생한다고 Minoru
> Akagi라는 분도 QuatumGIS Homepage에 bug 신고했더군요.
> 신 대표님 말씀데로 말끔히 고쳐지진 않을 듯 한데 어쨌든 나중에 nightly version 2.0 다시 설치해서 확인해 봐야겠네요.
> Arcgis가 만기가 되서 앞으로는 QuantumGIS로만 작업을 해볼까 했는데, spatial joining에서 걸리네요.
> spatial joining같은 핵심 기능은 언젠가는 제대로 고쳐졌으면 하고 기대해 봅니다.
>
> data type 읽는 문제는 지적하신대로 제가 처음 연습삼아 올렸던 석면 화일에 dirty data가 많아서 그랬던 것 같습니다.
> 건물 면적과 석면면적 필드에 일부 연도, 날짜 데이터 등이 섞여 들어간 셀이 있던데
> 이것들을 지우고 불러들이니 제대로 double로 다 읽히네요.
>
> 시험삼아 여러 다른 데이터도 QuantumGIS에 불러들여 지도 그리기를 해보다가 느낀 건데, Arcgis처럼 스타일 탭에서
> graduated mode의 데이터 구간별 분류를 할 때 histogram이 뜨면 좋을 텐데 하는 생각을
> 했습니다. 그래프를 보면서 수동으로 구간을 조정하는 기능이 있으면 참 유용할 텐데 말이죠. 1.8 version에서는 stats라는
> plugin에서 vector file의 속성을 histogram으로 보여주는 기능이 있어서 이를 참고로
> 구간을 수동으로 다시 입력하곤 했는데, 2.0에서는 histogram 그리는 plugin이 아직 없더군요.
>
> 현재 언론계에 gis user가 거의 없어서 대부분 혼자서 해결해야 할 부분이 많고 국내 지리관련 공공데이터가 적다보니 적극적으로
> 활용은 하지 못하고 있는데, 앞으로는 QuamtumGIS를 방송에 간간이 이용해 보려고 합니다. 아낌없이 조언을 주시는 신상희
> 대표님, 오동하 박사님, 박용재 연구원님 모두 거듭 고맙습니다.
> (박용재 연구원님이 올려주신 kernel density algorothm 사용법도 아주 유용한 자료이네요. 고맙습니다. 다만 제가
> 조작을 잘못한 건지 경우에 따라서 생각데로 안그려질 때도 있던데 이 부분도 나중에 정리해서 질문드려보도록 하겠습니다.)
>
>
> -함 형 건-
>
Reply all
Reply to author
Forward
0 new messages