Xác định hệ tọa độ của Shapefile khi chưa biết?

Tôi có một Shapefile nhưng hệ thống tọa độ của nó là Unknown và không có tệp * .prj.

Làm thế nào tôi có thể xác định nó bây giờ?

Luôn có phương pháp "vũ phu":

1. Lấy một lớp với một hệ tọa độ đã biết được cho là phủ lên lớp không xác định của bạn.

2. Bây giờ hãy thực hiện một số dự đoán có giáo dục về việc chiếu lớp chưa biết có thể là gì. (UTM, Tấm Carree, v.v.). Chiếu lớp hệ tọa độ đã biết của bạn vào từng phép chiếu cho đến khi bạn tìm thấy lớp phù hợp với lớp chưa biết càng nhiều càng tốt.

Tôi đang quảng bá ý kiến ​​của Mapperz và Brandon Copeland bằng cách thêm câu trả lời sử dụng kỹ thuật của họ.

1. Nếu bạn có ArcGIS Desktop, hãy thêm một số dữ liệu tham chiếu có cùng diện tích vào ArcMap. Dữ liệu tham chiếu này phải có định nghĩa hệ tọa độ (phép chiếu) hợp lệ.
2. Thêm dữ liệu trong một hệ tọa độ không xác định.
3. Nghiên cứu những gì hệ thống tọa độ có thể được sử dụng trong lĩnh vực này. Một nơi để kiểm tra là Sổ đăng ký EPSG .
4. Đặt hệ tọa độ của khung dữ liệu theo các khả năng. Nếu dữ liệu hợp âm dữ liệu khớp với dữ liệu không có dữ liệu hợp âm, có lẽ bạn đã xác định được hệ tọa độ.

Kỹ thuật này lợi dụng thực tế là ArcMap có thể chuyển đổi dữ liệu sang hệ thống tọa độ mới trong bộ nhớ. Các coordsys dữ liệu không có dữ liệu có thể được chuyển đổi để nó chỉ hiển thị. Bằng cách chiếu dữ liệu đã biết vào bộ nhớ, chúng ta có thể thấy nhanh chóng tọa độ của nó trong các hệ tọa độ khác nhau có thể.

Vì câu hỏi này không bao giờ cũ, tôi đã xây dựng một trang web thực hiện phương pháp Brute Force . Nếu bạn kéo shp đã nén + shx lên bản đồ, nó sẽ ánh xạ nó trong mọi hệ tọa độ có sẵn trong PostGIS. Giả sử bạn biết "chính xác" trông như thế nào, bạn có thể phóng to khu vực đó và nhấp vào đa giác để lấy tệp .prj từ epsg.io.

Có hai liên kết tuyệt vời từ Esri đi sâu vào chi tiết về điều này:

• Xác định một hệ tọa độ dự kiến ​​chưa biết bằng ArcMap
• Xác định tham chiếu không gian, phép chiếu hoặc hệ tọa độ của dữ liệu

Ở công việc cuối cùng của tôi, tôi có một lớp dữ liệu Địa chất (đa giác) có tên là "FSU_Geol.shp". Ông chủ của tôi đã đưa nó cho tôi và yêu cầu tôi tìm hiểu một số điều. Trước hết, anh ta được khách hàng trao chiếc shapefile này và không có tệp .prj, vì vậy anh ta muốn tôi tìm ra nó. Anh ấy cũng muốn tôi tìm ra các loại Địa chất. Tôi có thể tiếp tục, nhưng hãy cắt ngang để đuổi theo .... Tôi đã googled nó, và tôi đã kết thúc TẠI ĐÂY .

Từ viết tắt trong quy ước đặt tên là "Địa chất bề mặt của Liên Xô cũ" và tìm kiếm Google đã dẫn tôi "trực tiếp" đến nguồn (USGS). Tất cả mọi thứ tôi có thể cần biết về shapefile này là ở liên kết hàng đầu tôi nhấn. Tôi không nói rằng Google có thể tìm thấy bất cứ thứ gì & mọi thứ , nhưng tôi mới ra khỏi trường đại học, và chỉ cần "bắn trong bóng tối", và nhìn vào phản hồi tôi nhận được!

Trong một trường hợp khác, một người nào đó trong văn phòng đã tải xuống một số hình dạng bằng cách sử dụng hàm bó. Tôi không nhớ tên của gói phần mềm tại thời điểm này, nhưng các tệp .prj bị thiếu. Tôi chỉ đơn giản là đi vào siêu dữ liệu trong ArcCatalog và tôi thực sự tìm thấy câu trả lời của mình ở đó. Một lần nữa, đừng trích dẫn tôi, nhưng tôi nghĩ rằng có một URL đến nguồn trong đó.

Kể từ ngày đó, tôi "bắt đầu" với một tìm kiếm google nếu không có manh mối nào trong siêu dữ liệu!

Hỏi nhà sản xuất.

Nếu bạn biết CRS phổ biến cho khu vực địa lý của mình, bạn có thể thử một số trong số họ. Nhưng hỏi thì tốt hơn.

Tôi thực sự yêu thích công cụ web này: http://projfinder.com/ . Kiểm tra tập tin của bạn và tìm kiếm một số tọa độ. Thu phóng đến cùng một vị trí trên toàn cầu và để nó đoán.

Không phải là một công cụ (tôi không biết có một công cụ nào sẽ cho phép bạn làm điều đó), nhưng hãy xem trả lời của @ mkennedy về Làm thế nào tôi có thể chuyển đổi điểm này thành WKID 4326? . Cô giải thích làm thế nào cô đến tham chiếu không gian chính xác. SpatialReference.org và sự kiên nhẫn sẽ là bạn của bạn.

Ngoài ra, Esri cung cấp một hướng dẫn về cách đoán hệ tọa độ (mặc dù tôi thích phương pháp của mkennedy hơn nếu bạn biết thêm một chút về dữ liệu).

Câu trả lời tốt nhất tôi tìm thấy cho câu hỏi này là phi kỹ thuật: tìm hiểu dữ liệu của bạn đến từ đâu. Các cơ quan và tổ chức có xu hướng phù hợp với việc sử dụng các dự báo của họ. Biết nó đến từ DOT tiểu bang của bạn? Nhìn vào phần còn lại của dữ liệu của họ và xem những gì nó nói với bạn. Không biết nó đến từ đâu? Một phỏng đoán có giáo dục cũng giống như có thể đưa bạn đi đúng đường.

Ít nhất nó làm cho việc giải quyết vấn đề với lực lượng vũ phu hơn một chút có thể làm được!

Đã lỗi thời nghiêm trọng tại thời điểm này, nhưng Werner Flacke và Birgit Klaus đã đăng Tìm dự đoán trên ArcScripts vào năm 2007. Tôi không nghĩ rằng mã nguồn ở đó, thật không may. Nó dựa trên VBA nên chỉ có thể sử dụng được trong ArcGIS Desktop v9.2 và có thể là 9.3. Nó bao gồm hai shapefile với các lĩnh vực quan tâm từ Bộ dữ liệu tham số trắc địa EPSG có thể được sử dụng để thu hẹp các khả năng.

Blue Marble Geographic Calculator và Geographic Transformer có các công cụ phục hồi hệ thống phối hợp.

Hãy thử chương trình ogrinfo được cung cấp như một phần của GDAL.

Xem Cách truy cập siêu dữ liệu Shapefile bằng OGR?

Vì vậy, ogrinfo có thể không cung cấp cho bạn thông tin chiếu mà không có tệp .prj, nhưng nó vẫn là một công cụ hữu ích để giúp bạn điều tra danh sách các dự đoán tiềm năng.

Ví dụ:

Geometry: Polygon
Feature Count: 269 
Extent: (320000.000000, 505000.000000) - (323000.000000, 511000.000000) 

Phản hồi này gợi ý cho tôi rằng shapefile đang sử dụng một hệ thống tham chiếu dựa trên mét chứ không phải theo độ.

Giả sử bạn biết dữ liệu được lấy từ đâu, bây giờ bạn có một danh sách ngắn hơn các phép chiếu có thể.

Các thực hành tốt nhất khác có thể là:

• Hỏi những người đã cung cấp cho bạn dữ liệu ngay từ đầu, ngay cả khi họ không biết chính xác họ nên có một danh sách các dự đoán có thể có trong trường hợp đầu tiên.
• Nhìn để xem nếu có bất kỳ siêu dữ liệu liên quan.

Có lẽ tốt hơn để bắt đầu từ một số vị trí được biết đến, thay vì cố gắng và vũ phu từ đầu.

Để thêm vào các phản hồi khác trong bài viết này, tôi sẽ thêm vào như sau:

Cách thực hiện: Xác định hệ tọa độ dự kiến ​​chưa biết bằng ArcMap

Cung cấp một kiểm tra cụ thể về NAD1927

Nếu tọa độ ở độ thập phân, chẳng hạn như giữa kinh độ -180 và +180 và vĩ độ -90 và +90, hãy xác định Hệ thống tọa độ địa lý (mốc) được sử dụng cho dữ liệu. Trước phiên bản 9.2, ArcMap chỉ định GCS_Assumed_Geographic_1 cho dữ liệu theo mặc định. Điều này đặt dữ liệu trên mốc dữ liệu NAD_1927

Ngoài ra

Nếu dữ liệu ở Hoa Kỳ và hiển thị phạm vi tọa độ bên trái của số thập phân là 6, 7 hoặc 8 chữ số, dữ liệu có thể được chiếu vào hệ thống tọa độ Máy bay Nhà nước hoặc UTM.

Câu hỏi thường gặp: Khái niệm cơ bản về dự báo: Những gì chuyên gia về GIS cần biết

Trang này cung cấp một số kiểm tra về việc xác định các phép chiếu quan trọng trong việc hỗ trợ cố gắng xác định một hệ tọa độ không xác định.

Các khái niệm sau đây là cơ bản để hiểu việc sử dụng các phép chiếu bản đồ trong ArcGIS. Xin lưu ý rằng mặc dù chủ đề của các phép chiếu là rất rộng, và bài viết này không thể làm gì hơn là chạm vào một vài chủ đề quan trọng.

1. Các hệ tọa độ, còn được gọi là phép chiếu bản đồ, là các chỉ định tùy ý cho dữ liệu không gian. Mục đích của họ là cung cấp một cơ sở chung để liên lạc về một địa điểm hoặc khu vực cụ thể trên bề mặt trái đất. Vấn đề quan trọng nhất trong việc xử lý các phép chiếu bản đồ là biết được phép chiếu là gì và có thông tin hệ tọa độ chính xác liên quan đến tập dữ liệu.

2. Khi các phép chiếu bản đồ đầu tiên được phát minh, người ta cho rằng trái đất không bằng phẳng. Sau đó, giả định đã được sửa đổi, và trái đất được coi là một hình cầu hoàn hảo. Vào thế kỷ 18, mọi người bắt đầu nhận ra rằng trái đất không tròn hoàn hảo. Đây là sự khởi đầu của khái niệm hình cầu bản đồ.

3. Để thể hiện chính xác hơn các vị trí trên bề mặt trái đất, các nhà chế tạo bản đồ đã nghiên cứu hình dạng của trái đất (trắc địa) và tạo ra khái niệm hình cầu. Sau đó, các hệ tọa độ địa lý (GCS) đã được nghĩ ra, bao gồm mốc, đơn vị đo và kinh tuyến gốc. Một mốc thời gian liên kết một hình cầu với một phần cụ thể của bề mặt trái đất. Các mốc gần đây được thiết kế để phù hợp với toàn bộ bề mặt trái đất.

4. Các mốc được sử dụng phổ biến nhất ở Bắc Mỹ là:

• NAD 1927 (Bắc Mỹ Datum 1927) sử dụng hình cầu Clarke 1866
• NAD 1983 (Bắc Mỹ Datum 1983) sử dụng hình cầu GRS 1980
• WGS 1984 (Khảo sát trắc địa thế giới 1984) sử dụng hình cầu WGS 1984

Các nhân vật siêu anh hùng mới hơn được phát triển từ các phép đo vệ tinh và chính xác hơn so với các nhân vật do Clarke phát triển năm 1866. Thuật ngữ 'hệ tọa độ địa lý' và 'mốc' được sử dụng thay thế cho nhau, nhưng như đã lưu ý ở trên, GCS bao gồm thước đo, hình cầu, đơn vị đo và một kinh tuyến gốc.

5. Các tọa độ cho dữ liệu thay đổi tùy thuộc vào mốc và hình cầu mà các tọa độ đó dựa trên, ngay cả khi chúng đang sử dụng cùng một phép chiếu và tham số bản đồ.

Ví dụ: tọa độ địa lý bên dưới dành cho một điểm nằm trong thành phố Bellingham, Washington, sử dụng 3 mốc khác nhau:

Mã số: Tọa độ Y-tọa độ X-tọa độ NAD_1927 -122.466903686523 48.7440490722656 NAD_1983 -122.46818353793 48.7438798543649 WGS_1984 -122.46818353793 48.74387942

6. Một nguyên tắc quản lý dữ liệu tốt là lấy các tham số chiếu từ nguồn dữ liệu cung cấp dữ liệu. Đừng đưa ra phỏng đoán có giáo dục về việc chiếu dữ liệu, bởi vì cơ sở dữ liệu GIS không chính xác sẽ là kết quả. Các thông số cần thiết là như sau:

• Phép chiếu
• Đơn vị đo
• ZONE (cho UTM)
• Khu vực (đối với mặt phẳng trạng thái)
• Datum

Các thông số khác có thể được yêu cầu, tùy thuộc vào phép chiếu. Ví dụ, các phép chiếu Albers và Lambert yêu cầu các tham số sau:

• Song song tiêu chuẩn thứ nhất, tính bằng độ, phút và giây (DMS)
• Song song tiêu chuẩn thứ 2 (DMS)
• Kinh tuyến trung tâm (DMS)
• Vĩ độ của nguồn gốc dự đoán (DMS)
• Hướng đông sai và đơn vị đo lường
• Hướng bắc sai và đơn vị đo
• X-shift và đơn vị đo
• Y-shift và đơn vị đo

Tôi đã có cùng một vấn đề trong khi tôi đang làm việc trong dự án của mình (DEMs). Một trong những DEM là tọa độ không xác định và không thể chiếu, vì vậy tôi đã sử dụng công cụ 'Georeferences' sử dụng công cụ Georeferences trong ArcMap . Sau Georeferences, bạn sẽ nhận được prj.file của mình.

Nhìn vào trang web này:

http://www.egger-gis.at/shapefile-projectionfinder/

Tôi là nhà phát triển của công cụ miễn phí này. Có lẽ công cụ này có thể giúp bạn tìm và xác định hình chiếu chính xác của shapefile của bạn. Nó dựa trên ý tưởng về http://projfinder.com/ của Aaron Racicot.

Bạn cũng có thể thử các giải pháp đặc biệt này cho:

Áo (tiếng Đức): https://www.data.gv.at/anwendungen/checkaustrianprojection/

Úc: https://maegger.github.io/map_australia.html

UTM - Khu: https://maegger.github.io/map_utm.html

Bạn có thể phân biệt giữa tọa độ địa lý và bản đồ, lat / lon đầu tiên tính theo độ Bắc / Đông thứ hai tính bằng mét Nếu tệp có hình dạng, bạn có thể kết hợp tọa độ của hình học với shapelib-1.3.0 và sau đó xem các giá trị. ..

Bạn có thể làm điều đó với Fiona.

    import fiona as f
    a = fiona.open("C:\QGIS_ShapeFile\qgis\shafile_XXX.shp")
    print(a.crs)

Bạn sẽ nhận được câu trả lời là

    {'init': 'epsg:4326'}