NAR lon/lat datum: documented EPSG:4326 vs measured EPSG:4269

Đang mở
  • Có gì không ổn?
#4 0 bình luận 0 reaction 0 người được giao Xem trên GitHub

Chưa có ai nhận issue này.

Đánh giá

Độ khó
4/5
Thời gian dự kiến
3-5 ngày
Mức phù hợp với người mới
48/100
Loại issue
Lỗi
Độ rõ ràng
Khá rõ ràng
Mức độ hoạt động
Ít trao đổi
Công nghệ
r
Lĩnh vực
data

Hướng nghiên cứu

Bắt đầu với nar_lonlat_crs() và các entry point bị ảnh hưởng nar_point(), nar_lon(), nar_lat() và collect_nar(crs = ...). Kiểm tra StatCan NAR User Guide và bằng chứng round-trip đã được báo cáo cho EPSG:4326/EPSG:4269, sau đó xác nhận datum đã công bố với StatCan. Được xem là hoàn tất khi lựa chọn CRS đã được giải quyết và nar_lonlat_crs() cùng các hàm gọi bị ảnh hưởng phản ánh nhất quán câu trả lời đó.

Do mô hình lập chỉ mục viết ra từ nội dung của issue.

Mô tả

Issue

The StatCan NAR User Guide documents BG_LATITUDE/BG_LONGITUDE (and BF_REPPOINT_LATITUDE/_LONGITUDE) as EPSG:4326, but the data behaves as EPSG:4269 (NAD83). nar_lonlat_crs() is currently set to EPSG:4269 on that basis.

Evidence

Using NAR locations with exactly one address (so BG_X/BG_Y and the lon/lat describe the same point), re-deriving the projected coordinate from the lon/lat and comparing against BG_X/BG_Y over ~300k records:

assumed lon/lat CRS median residual
EPSG:4269 0.057 m (just NAR's 6-decimal rounding)
EPSG:4326 0.057 m (identical here — see below)
OGC:CRS84 1.08 m

StatCan evidently derived the projected and geographic coordinates without a datum shift between them.

Why it currently doesn't bite

On the dev machine PROJ lacks the NAD83↔WGS84 hgridshift grids, so EPSG:4326 and EPSG:4269 transform identically. On a system with those grids installed the two labels diverge by ~1 m. EPSG:4269 keeps the round-trip exact either way, which is why it was chosen.

Scope: affects nar_point(), nar_lon(), nar_lat() and collect_nar(crs = ...). geom is built from BG_X/BG_Y and is unaffected.

Follow-up

Confirm with StatCan which datum the published lon/lat columns are actually on, and whether the EPSG:4326 label in the User Guide is nominal. Revisit nar_lonlat_crs() once answered.

Reference: NAR User Guide, https://www150.statcan.gc.ca/n1/pub/46-26-0002/462600022026001-eng.htm

Ngôn ngữ chính
R
Star
0
Fork
0
Không có pull request nào được merge trong 30 ngày

Chuẩn bị môi trường

Dự án này không cung cấp dev container, Dockerfile hay hướng dẫn đóng góp, nên bạn cần tự thiết lập môi trường: hãy bắt đầu từ README và xem hướng dẫn đóng góp lần đầu của chúng tôi để biết các bước chung.

Bắt đầu từ đâu

  1. Đọc hết issue, rồi đọc hướng dẫn đóng góp của dự án.
  2. Bình luận trên issue rằng bạn sẽ nhận — tránh hai người làm cùng một việc.
  3. Fork repository và làm thay đổi trên một nhánh.
  4. Mở pull request có tham chiếu số hiệu của issue.

Issue khác của mountainMath/cangeocode

Tất cả issue của mountainMath/cangeocode

Issue tương tự

Thêm issue về R

Nhận issue mới trong hộp thư của bạn

Bản tóm tắt ngắn những issue GitHub phù hợp với người mới.