QGIS를 이용하여 shp -> svg 로 변환을 할 수 있지만 매핑된 정보(어트리뷰터)는 제외되어 SVG로 변환된다.

이때문에 jVectorMap에는 id와 name이 모두 빠져있다.
그럼 누군가가 이걸 수동으로 입력을 해야하는데 얼마나 그지 같은가...

구글링 도중 좋은 사이트를 발견하였다.


< 사이트의 초기화면 >


< shp 파일을 읽어오기 위해 왼쪽 하단의 버튼을 눌러 shp와 dbf를 차례로 불러온다. >

10 MB가 베스트라 되어 있지만 50 MB 까지 테스트 해보니 무리 없이 불러와진다.


< Load 버튼을 누르면 지도가 화면에 그려진다 >

 

 화면 중앙에 지도가 그려지고 좌측에 아래와 같은 레이어 창이 나타난다.


SHP와 DBF를 읽어들인 후 속성을 지도에 표시할 것인가에 대한 창이다. 
예를 들어 SIGUNGU_NM을 선택한 후 Okay 를 할 경우 아래와 같이 지도가 표시된다. 
 


얼마나 지저분한가.. 그러므로 하지말고.. 취소를 눌러서 작업을 하도록 하자.. 

 지도를 불러오게 될 경우 좌측 상단에 레이어창이 아래와 같이 나타난다. 


포토샵과 비슷한 레이어 창이다.
각 눈 표시를 누르면 각 레이어에 해당하는 오브젝트들이 사라지고 나타난다.

SVG로 변환할때 유의할 점은 Page와 Graticule를 해제하여 화면에 나타나지 않도록 한다.

[##_http://www.yeory.com/script/powerEditor/pages/1C%7Ccfile9.uf@2248514D51FB4394250C24.png%7Cwidth=%22504%22%20height=%2241%22%20alt=%22%22%20filename=%22%EC%9D%B4%EB%AF%B8%EC%A7%80%209.png%22%20filemime=%22image/jpeg%22%7C_##]화면에 표시된 지도를 SVG로 변환하기 위해 우측 상단에 export 버튼을 클릭

 

대화상자가 나타나는데 include data ~~~ 를 체크해야만 속성도 SVG에 포함되어 변환된다.

파일이름을 변경하고 Generate 버튼을 클릭하면 서버에서 SVG가 생성이되고 다운로드가 가능하다.


다운로드 받은 파일을 에이트플러스 또는 다른 편집툴로 열어보면 아래와 같은 형식이다.


꼭 편집기로 박스로 표시된 태그 전체를 지워주어야 한다.
외곽선과 그리드 선을 표시해주는 것으로써 지우지 않고 jVectorMap으로 변환할 경우 지도에 그대로 표시된다.

위 박스 태그를 모두 지우고 난 후 아래 사이트에서 지도로 변환하면 된다.

http://svgto.jvectormap.com

 


Result


Source

	function convert(latitude, longitude) {
		// lat : 위도, long : 경도
		
		// QGIS > SVG 플러그인 속성에 적은 가로-세로
		var width = 823;
		var height = 823;

		// QGIS > SVG 플러그인 속성에서 SET 된 각 한계선
		var eastLong = 132.054062;		// 우측 한계선
		var westLong = 124.410858;		// 좌측 한계선
		var northLat = 39.686848;		// 상측 한계선
		var southLat = 32.043644;		// 하측 한계선
		
		// 경도 -> x
		var longDiff = eastLong - westLong;
		var lon = (longitude - westLong) * (width / longDiff);
		
		// 위도 -> y
		var latDiff = northLat - southLat;
		var lat = (northLat - latitude) * (height / latDiff);
		
		return [ lon, lat];
	}

	$(document).ready(function() {

		width = $('#jvm_worldMap').width();
		mapCanvasHeight = $('#jvm_worldMap').height();

		// 부산 35.1333, 129.0500
		var markers = [ {
			coords : convert(37.60, 127.00),
			name : 'Seoul'
		}, {
			coords : convert(35.1333, 129.0500),
			name : '부산'
		}, {
			coords : convert(35.80, 128.58),
			name : '대구'
		}, {
			coords : convert(35.1755, 128.10),
			name : '진주'
		} ];

		$('#jvm_worldMap').vectorMap({
			map : 'korea_mill_en',
			markerStyle : {
				initial : {
					fill : '#F8E23B',
					stroke : '#383f47'
				}
			},
			markers : markers
		});
	});


HTML


jVectorMap은 친절히 맵의 Projection을 지원한다. 

국제 표준에 해당하는 Shape 파일을 jVectorMap Converter를 이용할 경우 map projection이 동일하게 변환되지만
한국 지도 행정동 2단계 지도를 변환해 보니 지도 군데 군데 구멍이 뚫리게 되어 이를 사용할 수 없게 되었다. 

관련 글 : http://www.yeory.com/270

 
이 때문에 완벽한 지도를 보여주어야 하기 때문에 QGIS를 이용하여 SVG로 변환하여 지도를 표시하였다.
이 경우 지도는 모두 나오는데 지도 타입(Projection)이 위/경도를 사용할 수 없기 때문에 위/경도로 좌표를 표시할 수 없다.

필요한 것은.... 
  • 지도의 가로 세로 -> 지도 자체의 가로, 세로 사이즈
    (div 등 jVectorMap을 이용하여 그려주는 태그의 사이즈가 아니다.)

  • 지도의 상하좌우 한계점의 위/경도 총 4개
    이 한계점 정보가 정확하지 않으면 X, Y 변환 값이 모두 틀어진다.

이 모든 정보들은 SHP -> SVG로 추출할 때 사용하는 QGIS SVG 플러그인에서 설정이 가능한 정보들이다.

작업 순서는 아래 그림처럼...

< 지도의 가로-세로 사이즈 설정 >

< 가로-세로 설정에 의해 틀어진 지도를 중앙에 위치하도록 설정 >

< 중앙에 맞춰진 지도의 각 상하좌우 위/경도를 설정해 주어야 한다>

마지막 이미지에서 보여지는 정보들이 스크립트에 설정해주어야 하는 정보이다.... 

maxY -> northLat
minY -> southLat
maxX -> eastLong
minY -> westLong

필요한 것이 나왔다면 아래 함수를 이용해 X, Y 변환이 가능하다.
	function convert(latitude, longitude) {
		// lat : 위도, long : 경도
		
		// QGIS > SVG 플러그인 속성에 적은 가로-세로
		var width = 792;
		var height = 611;

		// QGIS > SVG 플러그인 속성에서 SET 된 각 한계선
		var eastLong = 131.754044;		// 우측 한계선
		var westLong = 124.268181;		// 좌측 한계선
		var northLat = 38.749657;		// 상측 한계선
		var southLat = 32.974579;		// 하측 한계선
		
		// 경도 -> x
		var longDiff = eastLong - westLong;
		var lon = (longitude - westLong) * (width / longDiff);
		
		// 위도 -> y
		var latDiff = northLat - southLat;
		var lat = (northLat - latitude) * (height / latDiff);
		
		return [ lon, lat];
	}
경도 1도에 해당하는 픽셀을 구하고... (지도 넓이 / (좌-우)) 표시하고자 하는 경도를 좌측 값으로부터 뺀 후 곱하면 된다. 


모서리에 있는 점은 지도 가로-세로를 표시한 것이고, 위/경도는 이 안쪽을 기준으로 변환이 되어야만 한다.



하루 삽질했으나 원하는 결과가 나왔으니 만족...



업무 상 추세선을 그려야한다.

기존에 사용하던 유료 컴퍼넌트는 추세선을 그려주지 않는다..
(대체 이딴걸 돈주고 사용하는 이유가 뭘까..?)

 유료 컴퍼넌트를 대체하기 위해 찾은 2개의 차트 중에... 

Flot 차트의 플러그인으로 추세선을 그려주는 플러그인이 존재한다. 

jQuery.flot.trendline.js

한번 작업해보고.. 다시 포스팅...



추세선을 그렸다...

기존 플러그인을 조금 수정했는데 기존에는 From - To 형식으로 추세선을 그린다.

보다 정확히 그리고자 데이터 개수와 추세선 값을 1:1로 대응시켜 추세선을 그리게끔 변경하였다.


다운로드 : jquery.flot.trendline.js 

Full Source Download : trendLineChart

아래는 추세선을 표시한 차트 샘플이다.




'프론트엔드 > Flot Chart' 카테고리의 다른 글

Flot 차트에 추세선(trend line) ??  (0) 2013.07.29
ExtJS에서 컴퍼넌트를 만들어 쓴다고 했을때 크게 3가지의 방법이 존재한다. 

그 중에 가장 많이 쓰는 2개의 차이점을 간략하게 적어 놓는다.

 
Ext.widget  
  이를 사용하게 될 경우 기존에 존재하는 컴퍼넌트의 인스턴스 쿼리를 반환하지 않는다. 
  즉 내가 만드는 컴퍼넌트가 다른 방법(create, define)으로 만들어져 있더라도 가져오지도 않고,
  새로 만들지만 다른 곳에서 찾아서 쓸 수 있는 쿼리를 리턴하지 않는 다는 것..
  
Ext.create
  이를 사용하게 될 경우 항상 해당 xtype의 인스턴스로 새로이 생성한다.
  생성과 함께 인스턴스를 찾아 쓸 수 있도록 쿼리가 만들어진다.
  만약 같은 컴퍼넌트를 각각 10곳에서 create 했다면 이론적으로 10개의 인스턴스가 존재하게 되는 것.


무조건 widget만 쓰는것도 안좋고, 그렇다고 create만 죽어라 하는 것도 좋지 않다.

하나의 폼을 10개의 화면에서 사용을 한다면..? 
10개의 화면에서 모두 create를 할것인가..? 아니면 widget을 할 것인가?? 

create와 widget은 임시, 1회용 또는 동적으로 폼의 형태를 변형할 경우에 적합하다고 생각해본다.

아니면 말고 
jVectorMap은 자체적으로 convert를 사용하여 shp 파일의 vector 이미지를 svg로 변환하여 맵 파일을 생성한다.

그러나 이 컨버터를 이용할 경우 svg 변환 상에 오류가 있는것인지 대한민국 지도 상에 구멍이 슝슝 뚫린다.

아래처럼... 


대구와 대전 일부가 없다.

이를 해결하기 위해 제작자에게 메일도 보냈지만 묵묵부답이라 자체적으로 해결 할 수 밖에... 

1. shp 파일을 편집할 수 있는 (읽어 들일 수 있는)  프로그램을 찾는다.
2. 이 프로그램을 이용하여 SVG로 변환할 수 있는지 찾는다.
3. SVG로 변환된 정보를 svgto.jvectormap.com 에서 지도로 변환해 본다.

간단한거 같지만 이게 은근 쉽지가 않았다..
개발자로써 지도 프로그램을 쓰는 것도 그렇고 사용법을 일일이 찾아야 하므로... 

첫번째 사용했던 프로그램 FME Desktop 2013 이란 상용 프로그램이다.

14일간의 트라이얼을 기대하고 설치 했는데... 트라이얼 코드를 웹사이트에 접속해서 등록해야만 코드가 발급된다. ;;
여튼, 이 프로그램으로 좋은 것은.. SHP 파일을 정말 자세하게 까(?) 볼 수 있었다.

  1. FME Data Inspector

    shp 파일은 벡터 정보와 데이터 셋을 같이 가지고 있다. data inspector는 이 데이터 셋을 집중적으로 깔 수 있는 툴이다.

    아래 이미지는 data inspector를 이용하여 kor_adm2.shp 파일을 열었을때의 화면..



  2. FME Universal Viewer
    이 툴은 shp 파일의 벡터 부분을 집중적으로 깔 수 있다. 라인과, 포인트, 텍스트, 채우기 색깔등을 볼 수 있는 것.
    그리고 각 조각을 선택하면 그에 대한 정보를 보여준다. (사실 inspector와 거의 같은 정보..)
    이미지는 생략하고..

    이 툴을 이용해서 shp -> svg로 변환이 가능하다.
    shp가 그린 벡터 조각이 가진 어트리뷰터를 모드 svg path 태그안의 어트리뷰트로 변환이 된다.

    그런데!! 이걸 svgto.jvectormap.com 에서 지도로 변환해보면 지도가 안나온다.

    왜인지 모르겠는데 안나온다...

     
  3. FME Workbench
    이 툴은 전문가 적인 작업이 가능하다.

    내가 불러들여 작업할 파일을 설정하는데 이를 Reader라 하고
    어떤식으로 변환할 것인지를 설정해야하는데 이를 Writer라고 부른다.

    아래는 이에 대한 이미지..


    이렇게 리더와 Writer를 설정하고 OK 누르면 아래 이미지와 같은게 나온다..


    Writer 쪽 ... 을 눌러보면 User Attribute 탭 에서 내가 추출하고자 하는 Reader의 Attribute를 지정할 수 있다.
    Reader와 같은 이름을 지정해야만 SVG로 변환 할 시 나타나며 이름이 틀릴경우 나타나지 않는다.



    대충 설정을 하고 Shift + F5 를 누르게 되면 Reader -> Writer로 변환을 하게 되는데 jVectorMap으로 안나온다.

     
FME Desktop 2013 결론은 jVectorMap에서 사용하는 SVG를 추출할 수 없다.

내가 못하는 거 일 수 있지만 안된다.  그래서 다음 프로그램으로 패스...


두번째 사용한 프로그램은 무료로 사용할 수 있는 Quantum GIS 이다.
Quantum GIS : http://www.qgis.org/
한글 언어를 지원하는 프로그램으로 별 어려움 없이 사용 할 수 있었다.

설치 후 프로그램을 실행시켜 shp를 읽어들이고자 할때는 아래 아이콘을 이용한다. 


벡터 레이어를 추가하여 데이터 셋을 읽어들이면 되는데 사용하고자 하는 shp 파일을 로드하면 된다.

이제 SVG로 변환을 하여야 하는데 기본적인 기능에서는 제공되지 않는다.

그래서 플러그인을 설치하여야 한다. 

메뉴에서 플러그인 -> Python 플러그인을 눌러 플러그인 설치 관리자를 실행 시키고 필터에 Simple SVG 로 검색한다.
하나 밖에 없으므로 눌러서 플러그인 설치를 하면 SVG로 변환할 준비 끝.

설치가 완료되면 플러그인 툴바가 프로그램 하단에 생성된다.
설치된 SVG 플러그인을 누르면 창이 뜨는데 창 가운데 Resize Map or Set Extent를 눌러서 크기를 조절해야한다.




위 이미지에서 처럼 Width와 Height의 화살표를 눌러 값을 바꾸고 우측의 Set을 누르면 지도 전체 크기가 변경된다.
그리고 꼭 !!  빨간색 상자로 표시된 버튼을 눌러 벡터 이미지도 지도 크기에 맞게 사이즈 조절을 해야만 한다.
(창을 닫지 않고도 눌러 지므로 창 닫지는 말자..)

빨간색 상자 버튼을 누르고 난 후 아래 Set Extent 의 Set 버튼을 눌러 표시된 지도의 X,Y를 설정하자.
(누르면 자동으로 맞춰진다.)

OK를 누르고 첫 창으로 돌아와 Browse를 눌러 저장하고자 할 경로와 파일명.svg 를 적어주고 OK를 누르면 변환이 완료된다.

이 변환된 파일을 에디트플러스나 다른 편집 프로그램으로 열어 아래 사이트에서 실제 지도로 변환되는지 확인해보자.

here > http://svgto.jvectormap.com/ 


아래는 변환된 결과. 


구멍 하나 없이 제대로 나온다. 

그런데!! 이 Quantum GIS를 이용하여 SVG로 변환할 경우 Shape가 가지고 있는 attribute를 같이 변환하지 않는다.
Simple SVG는 레이어 단위로 변환을 하기 때문에 id와 name이 따라가려면 레이어로 다 쪼개야하는데... 못하겟다;;

여튼 저 화면에서 Map name과 지역별 id, name을 설정하고 save를 누르면 자바 스크립트 파일로 다운로드가 가능하고
이를 jVectorMap에서 사용할수 있다. 


이제 어떤 파일이든지 Quantum GIS 만 있으면 jVectorMap에서 사용할 수 있다!!
그러나 이름은 손으로 다 입력해야한다. 


SVG Path만 있다면 jVectorMap으로 가능하다. 

아래는 SVG 샘플...  








위 샘플을 아래 사이트에서 확인할 수 있다.

svg to jvectormap : http://svgto.jvectormap.com

 		Ext.Ajax.request({
	        url: '../com/ProxyJsonSimpleOneList.jsp',
	        contentType: 'json',
	        params: {
	        	name: 'value',
	        	today: new Date(),
	        	code: Ext.getCmp( 'brcCode' ).value
	        },
	        success: function(response){
	        	console.log(response);
	        },
	        failure: function(response, opts) {
	        	commonExtAjaxFailure(response, opts);
	        }
	    });





dddd

1. GADM
    site : http://www.gadm.org/country
   가장 많이 보편적으로 이용하는 사이트. 각 국가별로 다운로드 할 수 있으며 다양한 포맷을 지원한다.
   그리고 행정동 레벨을 5 단계로 나누어 놓았고, 레벨별 지원하는 국가도 확인할 수 있다.

2. BIZ-GIS
   site : http://www.biz-gis.com/GISDB/
   국내용으로 다양한 지도의 데이터를 다운로드 할 수 있다. 꼭 좌표계를 확인해서 좌표 종류를 확인하자.

3. 통계청 제공 지도
    site : http://sgis.kostat.go.kr/statbd/statbd_03.vw
    통계청에서 행정구역 경계를 표시하기 위해 제작된 것으로 shp 다운로드를 제공한다.
    이 글 작성일 기준으로 2011년까지 제공.


준비물... 

1. python 2.7 이상 
3. GDAL
2. OSGeo4W


python이야 설치 파일 받아서 설치하면 되므로 간단하다. 정말..

GDAL 또한 아래 사이트에서 자기에 맞는 버전 찾아서 적당한 위치에 압축해제 하면 된다.
donwload : http://download.osgeo.org/gdal/win32/1.5/ 

해제 후 시스템 변수 중 path에 압축해제 한 폴더를 추가 해 주어야 하고.. 
압축 해제 폴더가 gdalwin32-1.7 이라면...  path에 C:\gdalwin32-1.7\bin 를 추가..
 

JAVA_HOME 을 설정할때와 마찬가지로..

Name : GDAL_DATA
Value : C:\gdalwin32-1.7\data

를 넣어주자.. 

이제 필요 한것은  Shapely 아래 사이트에서 다운로드 한 후 압축해제하고... 
downlaod : https://pypi.python.org/pypi/Shapely/1.2.16 

커맨드 창으로 python setup.py install 해주면 된다.

마지막으로.. 필요 한 것이 ogr 패키지 인데... 이건 OSGeo4W를 설치하면 쉽게 해결 할 수 있다.
download : http://trac.osgeo.org/osgeo4w/ (관리자 권한으로 설치 할 것)

installer를 다운받아 설치를 시작하면 몇 가지 플러그인을 설치할 것이냐고 묻는데 아무거나 해도 관계없다.. 

다 설치하게 되면 준비 끝.

shp 파일이 준비되었으면 cmd 창으로 converter.py를 실행시킬것이지만... 에러가 난다.. 분명.. 

시작을 누른후 OSGeo4W 라고 입력하면 프로그램을 실행 시킬 수 있다.
이 쉘을 실행 시킨 후 converter.py 를 실행시키면 아무 문제 없이 실행 시킬 수 있다.

그리고 피씨가 정말 후지지 않는 이상 컨버팅 작업 중 멈추지 않을 것이다..

아래처럼 대충 이렇게... 

python ../converter.py ^

--width 900 ^

--country_name_index 6 ^

--country_code_index 5 ^

--where "ISO = 'KOR'" ^

--name korea ^

--buffer_distance -0.4 ^

--simplify_tolerance 50 ^

KOR_adm2.shp ^

korea-map.js















 

+ Recent posts