안녕하세요.

저는 Hadoop기반에서 R패키지로 경매, 토지실거래가 및 각종 공간데이터를 이용하여 토지가치분석결과를 표출하는 대시보드를 만들고 있습니다.


지금은 공간정보포털에서 '연속지적도(shp)'를 geojson으로 변경하여 hadoop/hive에 올리고(geometry는 binary로 저장),

그것을 다시 hive에서 string 포멧으로 가져와서 spatialPolygon으로 만들고,

이를 OpenstreetMap에 leaflet을 이용하여 대시보드에 표출하려고 합니다.

그러나, 공간정보포털의 shp파일을 "좌표변환"하는데, 혼자서 해결이 되지 않아서 질문을 올립니다.(polygon의 회전에 문제가 있는 듯...)


Qgis에서 이런 저런 CRS정보를 참조하여 설정을 변경하면서 적용을 해 보았습니다.

EPSG:5181로 CRS를 적용하면 다음map이든 어디든 잘 올라가는 듯 햇는데..


R패키지에서 장흥군의 지적도를 Bassel/TM으로 잡고 OpenstreetMap에 올려보니, 그나마 비슷하게 올려지긴 했는데, 한참 부족하고...

여기저기 구글링을 해봐도 완전한 해결방법을 찾을 수가 없네요.


[ spatialPolygon으로 생성하기 위한 기본 설정 ] - http://www.biz-gis.com/index.php?mid=QnA&document_srl=222803

(EPSG:5173) +proj=tmerc +lat_0=38 +lon_0=125.0028902777778 +k=1 +x_0=200000 +y_0=500000 +ellps=bessel +units=m +no_defs +towgs84=-115.80,474.99,674.11,1.16,-2.31,-1.63,6.43


[ OpenstreetMap에 올리기 위해 좌표변환에 적용한 설정 ] - https://www.osgeo.kr/17

(EPSG:3857) +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs


결과를 보면 위치도 조금 다른 것 같아, lon/lat과 축적값을 조금씩 변경하며 테스트 해보면 근접까지는 되는데,

Polygon의 회전은 안되네요...

혹시, 공간정보포털의 공개데이터로 좌표변환 작업에 대한 경험이 있으신 분이 계시다면 조언 부탁드립니다.



* 아래 이미지는 위의 설정값으로 R/leaflet 패키지를 이용하여 법정경계(장흥군)를 올려 본 최종 화면입니다.

장흥군.PNG : 공간정보포털의 좌표체계에 대해 문의 드립니다.

profile

곡근

2019.08.19 07:24
*.144.157.134

자문자답...

직접 답을 찾아 댓글을 달게 될 줄은.........  ㅜㅜ


적절한 방법 중에 두가지가 있네요.


첫째 : 직접 회전함수를 적용하여 데이터 변경

coords.90 <- Rotation(coords, 90*pi/180)

# coords : A 2-columns matrix or data frame containing a set of X and Y coordinates.


둘째 : CRS의 프로젝트 스트링에 회전값을 직접 적용

  # convert to lat-long
  boxLL = bbox(spTransform(boxpts, CRS("+init=epsg:4326")))
  # find the centre
  llc = apply(boxLL, 1, mean)
  # construct the proj4 string
  prj = paste0("+proj=omerc +lat_0=", llc[2], " +lonc=", llc[1], " +alpha=",
               angle
, " +gamma=0.0 +k=1.000000 +x_0=0.000 +y_0=0.000 +ellps=WGS84 +units=m ")
  # return as a CRS:
  CRS(prj)

List of Articles
번호 제목 글쓴이 날짜 조회 수
공지 질문올리 실때 조금만 질문 내용을 조금만 더 자세히 올려주시면 답변받기 수월합니다. + 11 김한국 2008-01-04 239732
6901 포인트(또는 그리드) 內 버퍼가 몇 개 포함되어있는지 카운트해서 새로운 필드를 만들고 싶습니다. file + 2 onnie 2019-07-29 169
» R패키지를 이용하여, 공간포털의 연속지적도(shp)를 openstreetMap에 올리기 위한 좌표변환에 대하여 file + 1 곡근 2019-07-28 268
6899 arcgis10.3에서 평균경사도를 구하고 싶습니다. file + 1 박수영 2019-07-28 293
6898 Biz-gis 자료실에 자연과학 GIS DB 한글이 깨집니다. + 1 공주지질 2019-07-27 102
6897 Loading Table Data 오류 관련 질문 드립니다. (Open Table Attribute) file + 1 다이라트스 2019-07-26 75
6896 폴리곤과 교차되는 포인트 개수 세기 + 2 다릇 2019-07-25 145
6895 인터섹트나 뭐 이런걸 하면 예전엔 가운데 진행률이 떳는데 그게 지금 안뜹니다 어떻게 해결하나요 + 1 노망 2019-07-25 76
6894 도로명주소 전자지도와의 좌표 변환 문제 문의드립니다...ㅠ file + 1 닉닉젠 2019-07-23 203
6893 피처 크기에 따라 속성값 분배하는 법 file + 1 다릇 2019-07-22 86
6892 다른 공간데이터로부터 속성정보를 조인하려고 합니다. file + 1 바나나모양똥 2019-07-22 171
6891 항공사진을 이용한 DEM 구축방법....??? + 2 Micia 2019-07-22 258
6890 gis를 이용하여 풍수해과업에서 하천재해 및 토사재해를 할려고 합니다 수자원빠돌이 2019-07-22 75
6889 래스터 클립 후 Intersect 이용 file + 3 다릇 2019-07-22 150
6888 arcgis ArcToolbox 에서 원하는 필드만 export 하는게 있나요? + 1 우르르르르 2019-07-19 117
6887 GIS를 이용한 고지형복원(3D) 질문드립니다 file 고고학자 2019-07-18 92
6886 포인트피쳐를 중심으로 폴리곤 생성 및 속성 분배하기 file + 2 다릇 2019-07-18 177
6885 완전 초보가 좌표계 관련해서 질문합니다. file + 5 청소하자 2019-07-17 200
6884 버퍼가 작동되지 않아요 ㅜㅜ file + 1 민소리 2019-07-17 108
6883 추출한 버텍스의 좌표값을 위경도값으로 변환하는 방법 문의드립니다. file + 3 마름모꼴 2019-07-17 287
6882 폴리곤 분할 file + 1 뱀파이어디 2019-07-17 155
6881 연속지적도 SHP 파일에서 지번별 POLYGON 좌표값을 추출하는 방법이 있을까요? file + 1 마름모꼴 2019-07-16 319
6880 두 포인트의 간격이 몇 미터 이하일 경우만 추출하는 방법 가르쳐주세요 file + 1 발루아블레 2019-07-15 64
6879 WGS84 좌표계를 이용해서 1km x 1km 그리드 만드는 부분에 대해서 질문드립니다. + 1 걸음마수준123 2019-07-13 112
6878 아크gis 10.3 인데 틴생성이 되질 않습니다. + 2 이명호 2019-07-11 126
6877 지형도랑 지적도를 겹치려고 하는데 file + 3 달다인생 2019-07-11 347
6876 도곽 좌표계에 대한 질문 file 조형 2019-07-11 74
6875 속성정보 조인시 숫자 및 문자 인식 오류 file + 1 gis_beginner 2019-07-09 164
6874 좌표관련 문의드립니다 dofla 2019-07-09 61
6873 MXD 파일 저장속도 및 Repair Geometry 레이어 초기화 질문 드립니다. file + 3 다이라트스 2019-07-09 166
6872 R에서 격자화 한 후 데이트 값 입력하는 방법에 대해서 질문드립니다. 걸음마수준123 2019-07-08 110
6871 쉐입파일 한번에 Clip하는 방법 문의드립니다 + 2 퇴사의꿈 2019-07-07 179
6870 QGIS 포인트 사이 거리 측정하기 file 댕석멍박 2019-07-06 204
6869 QGIS 설치 오류 file + 1 댕석멍박 2019-07-06 125
6868 국토지리정보원에서 받은 DEM 해상도 낮은데 여러분들은 어디서 제공 받나요...? + 2 Micia 2019-07-05 335
6867 QGIS 툴박스 오류 댕석멍박 2019-07-05 108
6866 [ArcGIS Pro] 기준좌표 모를때 좌표정의 방법 file + 5 다릇 2019-07-03 379
6865 6 자리 정수는 무슨 좌표계일까요? + 1 littleprince 2019-06-29 280
6864 버퍼와 레스터 관련 질문입니다 file 송문수 2019-06-28 112
6863 slope시 Z 계수 오류가 납니다. file glauconite 2019-06-27 170
6862 tin 설정 후 경계편집 file + 1 glauconite 2019-06-27 198
6861 좌표 없는 정사영상 좌표 설정 어떻게 하나요...? + 5 Micia 2019-06-27 357
6860 Arcgis에서 두 레이어가 있을 때 객체들이 서로 겹처있는 경우를 찾는 방법이 있나요? file + 1 발루아블레 2019-06-23 167
6859 gis 툴을 찾고싶어요 + 1 klkllk 2019-06-23 124
6858 Error exporting data.. + 2 하루사이 2019-06-21 118
6857 위도가 뒤집혀 있어서 질문드립니다. file 송문수 2019-06-20 133
6856 arcmap 출력시 오른쪽하단에 글씨 없애는방법 알려주세요 + 1 양승찬 2019-06-19 100
6855 고도값을 고려한 거리계산 jayz 2019-06-18 99
6854 안녕하세요 CSV파일과 SHP파일을 이용하여 맵에다가 표시를 할 수 있을지 질문드립니다. + 1 이인입니다 2019-06-18 245
6853 수계추출 및 면적계산 방법 좀 알려주세요~~ file xiu 2019-06-14 138
6852 이미지를 넣으면 진하게 나와요 file JYYYYY 2019-06-14 103