การระบุระบบพิกัดของ Shapefile เมื่อไม่ทราบ?

ฉันมี Shapefile แต่ระบบพิกัดไม่เป็นที่รู้จักและไม่มีไฟล์ * .prj

ฉันจะระบุได้อย่างไรตอนนี้

มีวิธี "กำลังดุร้าย" อยู่เสมอ:

1. ใช้เลเยอร์ด้วยระบบพิกัดที่รู้จักซึ่งควรจะซ้อนทับกับเลเยอร์ที่คุณไม่รู้จัก

2. ทีนี้ลองเดาดูว่าเลเยอร์เลเยอร์ที่ไม่รู้จักนั้นเป็นอย่างไร (UTM, Plate Carree ฯลฯ ) ฉายเลเยอร์ระบบพิกัดที่คุณรู้จักลงในแต่ละการฉายภาพจนกว่าคุณจะพบเลเยอร์ที่ตรงกับเลเยอร์ที่ไม่รู้จักมากที่สุด

ฉันส่งเสริมความคิดเห็นโดย Mapperz และ Brandon Copeland โดยเพิ่มคำตอบที่ใช้เทคนิคของพวกเขา

1. หากคุณมี ArcGIS Desktop เพิ่มข้อมูลอ้างอิงบางส่วนที่ครอบคลุมพื้นที่เดียวกันกับ ArcMap ข้อมูลอ้างอิงนี้ต้องมีข้อกำหนดระบบพิกัด (การฉายภาพ) ที่ถูกต้อง
2. เพิ่มข้อมูลในระบบพิกัดที่ไม่รู้จัก
3. ทำการวิจัยว่าระบบพิกัดใดที่เป็นไปได้ที่จะใช้ในด้านนี้ หนึ่งในสถานที่ที่จะตรวจสอบเป็นRegistry EPSG
4. ตั้งค่าระบบพิกัดของ data frame ให้เป็นไปได้ หาก data-with-a-รู้จัก coordsys สอดคล้องกับ data-without-a-รู้จัก -cordsysys คุณอาจระบุระบบพิกัด

เทคนิคนี้ใช้ประโยชน์จากข้อเท็จจริงที่ว่า ArcMap สามารถแปลงข้อมูลเป็นระบบพิกัดใหม่ในหน่วยความจำ data-without-a-รู้จัก -cordsysys ไม่สามารถแปลงดังนั้นมันเพิ่งแสดง โดยการฉายข้อมูลที่รู้จักในหน่วยความจำเราจะเห็นได้อย่างรวดเร็วว่าพิกัดใดที่มีในระบบพิกัดต่าง ๆ ที่เป็นไปได้

เนื่องจากคำถามนี้ไม่เคยได้รับเก่าฉันสร้างเว็บไซต์ที่ไม่วิธี Brute Force หากคุณลาก zp ที่ซิป + shx ลงบนแผนที่มันจะทำการแมปในระบบพิกัดทุกระบบที่มีอยู่ใน PostGIS สมมติว่าคุณรู้ว่า "ถูกต้อง" ดูเหมือนว่าคุณสามารถซูมไปที่พื้นที่นั้นและคลิกที่รูปหลายเหลี่ยมเพื่อรับไฟล์. prj จาก epsg.io

มีสองลิงค์ที่ยอดเยี่ยมจาก Esri ที่มีรายละเอียดดังนี้:

• ระบุระบบพิกัดที่ไม่รู้จักที่คาดการณ์ไว้โดยใช้ ArcMap
• ระบุการอ้างอิงเชิงพื้นที่การฉายภาพหรือระบบพิกัดของข้อมูล

ในงานสุดท้ายของฉันฉันมีชั้นของข้อมูลธรณีวิทยา (รูปหลายเหลี่ยม) ชื่อ "FSU_Geol.shp" เจ้านายของฉันมอบมันให้ฉันและขอให้ฉันค้นหาสิ่งต่าง ๆ ก่อนอื่นเขาได้ส่งไฟล์รูปร่างนี้โดยลูกค้าและไม่มีไฟล์. prj ดังนั้นเขาจึงต้องการให้ฉันคิดออก เขายังต้องการให้ฉันเข้าใจว่าธรณีวิทยานั้นเป็นอย่างไร ฉันจะไปบนและบน แต่ช่วยลดการล่าสัตว์ .... ฉัน googled มันและฉันจบลงที่นี่

คำย่อในอนุสัญญาการตั้งชื่อหมายถึง "ธรณีวิทยาพื้นผิวของอดีตสหภาพโซเวียต" และการค้นหาของ Google ทำให้ฉัน "ตรง" กับแหล่งที่มา (USGS) ทุกสิ่งที่ฉันจำเป็นต้องรู้เกี่ยวกับ shapefile นี้คือที่ลิงค์ด้านบนที่ฉันกด ฉันไม่ได้บอกว่า Google สามารถหาอะไรได้ทุกอย่างแต่ฉันยังสดใหม่จากมหาวิทยาลัยและแค่ "ถ่ายภาพในที่มืด" และดูความคิดเห็นที่ฉันได้รับ!

ในอินสแตนซ์อื่นบางคนในสำนักงานได้ดาวน์โหลดรูปร่างไฟล์จำนวนหนึ่งโดยใช้ฟังก์ชั่นแบบแบทช์ ฉันจำชื่อของซอฟต์แวร์แพ็กเกจไม่ได้ในตอนนี้ แต่ไฟล์. prj หายไป ฉันเข้าไปที่เมตาดาต้าใน ArcCatalog และฉันก็พบคำตอบที่นั่น อย่าอ้างฉันอีก แต่ฉันคิดว่ามี URL ไปยังแหล่งข้อมูลในนั้น

นับตั้งแต่วันนั้นฉัน "เริ่มต้น" ด้วยการค้นหาของ Google หากไม่มีข้อมูลในเมตาดาต้า!

ถามผู้ผลิต

หากคุณรู้จัก CRS ทั่วไปสำหรับพื้นที่ทางภูมิศาสตร์ของคุณคุณสามารถลองใช้งานได้บ้าง แต่การถามดีกว่า

ผมรักเครื่องมือบนเว็บนี้: http://projfinder.com/ ตรวจสอบไฟล์ของคุณและค้นหาพิกัด ซูมไปที่ตำแหน่งเดียวกันโดยประมาณบนโลกแล้วปล่อยให้เดา

ไม่ใช่เครื่องมือ (ฉันไม่รู้ว่ามีอยู่จริงที่จะให้คุณทำอย่างนั้น) แต่ลองดูคำตอบของ @ mkennedyต่อฉันจะแปลงประเด็นนี้เป็น WKID 4326 ได้อย่างไร . เธออธิบายว่าเธอมาถึงการอ้างอิงเชิงพื้นที่ที่ถูกต้องได้อย่างไร SpatialReference.orgและความอดทนจะเป็นเพื่อนของคุณ

นอกจากนี้ Esri ยังให้คำแนะนำเกี่ยวกับวิธีการคาดเดาระบบพิกัด (แม้ว่าฉันจะชอบวิธีการของ mkennedy ถ้าคุณรู้ข้อมูลเพิ่มเติมเล็กน้อย)

คำตอบที่ดีที่สุดที่ฉันพบสำหรับคำถามนี้คือไม่ใช่ด้านเทคนิค: ค้นหาว่าข้อมูลของคุณมาจากไหน หน่วยงานและองค์กรมีแนวโน้มที่จะสอดคล้องกับการใช้ประมาณการของพวกเขา รู้ว่ามันมาจากรัฐ DOT ของคุณ? ดูข้อมูลที่เหลือของพวกเขาและดูว่ามันบอกอะไรคุณ ไม่รู้ว่ามันมาจากไหน? การเดาที่ได้รับการศึกษานั้นมีแนวโน้มที่จะส่งคุณไปสู่เส้นทางที่ถูกต้อง

อย่างน้อยก็ทำให้การแก้ปัญหาด้วยกำลังดุเดือดมากขึ้นอีกนิด!

ล้าสมัยอย่างจริงจัง ณ จุดนี้ แต่ Werner Flacke และ Birgit Klaus โพสต์Find Projectionบน ArcScripts ในปี 2550 ฉันไม่คิดว่าซอร์สโค้ดจะอยู่ที่นั่นโชคไม่ดี มันใช้ VBA เพื่อให้สามารถใช้งานได้เฉพาะใน ArcGIS Desktop v9.2 และอาจเป็น 9.3 มันรวมไฟล์สองรูปร่างที่มีพื้นที่ที่น่าสนใจจากชุดข้อมูลพารามิเตอร์จีพีดีจีโอเดติกซึ่งสามารถใช้เพื่อ จำกัด ความเป็นไปได้

เครื่องคิดเลขทางภูมิศาสตร์Blue Marbleและ Geographic Transformer มีเครื่องมือประสานงานการกู้คืนระบบ

ลองใช้โปรแกรม ogrinfo ที่ให้มาเป็นส่วนหนึ่งของ GDAL

ดูวิธีเข้าถึงข้อมูลเมตาของ Shapefile โดยใช้ OGR

ดังนั้น ogrinfo อาจไม่ให้ข้อมูลการฉายโดยไม่มีไฟล์. prj แต่ก็ยังเป็นเครื่องมือที่มีประโยชน์ที่จะช่วยคุณตรวจสอบรายการของการคาดการณ์ที่อาจเกิดขึ้น

ตัวอย่างเช่น:

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

การตอบสนองนี้ชี้ให้ฉันเห็นว่า Shapefile กำลังใช้ระบบอ้างอิงโดยอ้างอิงจากเมตรและไม่ใช่เป็นองศา

สมมติว่าคุณรู้ว่าข้อมูลมาจากไหนตอนนี้คุณมีรายการประมาณการที่สั้นกว่า

แนวทางปฏิบัติที่ดีที่สุดอื่น ๆ อาจเป็น:

• ถามคนที่ให้ข้อมูลกับคุณตั้งแต่แรกแม้ว่าพวกเขาจะไม่รู้ว่าพวกเขาควรมีรายการของการคาดการณ์ที่เป็นไปได้ที่สามารถดูได้ในตัวอย่างแรก
• ดูเพื่อดูว่ามีเมทาดาทาที่เกี่ยวข้องหรือไม่

มันอาจจะดีกว่าถ้าเริ่มจากตำแหน่งที่รู้จักกันดีกว่าลองและกำลังดุร้ายตั้งแต่เริ่มต้น

เพื่อเพิ่มคำตอบอื่น ๆ ในโพสต์นี้ฉันจะเพิ่มต่อไปนี้:

วิธีการ: ระบุระบบพิกัดที่ไม่รู้จักที่คาดการณ์ไว้โดยใช้ ArcMap

ระบุการตรวจสอบเฉพาะใน NAD1927

หากพิกัดอยู่ในรูปแบบทศนิยมเช่นระหว่างลองจิจูด -180 ถึง +180 และละติจูด -90 และ +90 ให้ระบุระบบพิกัดทางภูมิศาสตร์ (ตัวเลข) ที่ใช้สำหรับข้อมูล ก่อนหน้าเวอร์ชัน 9.2 ArcMap จะกำหนด GCS_Assumed_Geographic_1 ให้กับข้อมูลตามค่าเริ่มต้น สิ่งนี้จะวางข้อมูลบน NAD_1927 datum

นอกจากนี้

หากข้อมูลอยู่ในสหรัฐอเมริกาและแสดงขอบเขตซึ่งพิกัดทางด้านซ้ายของทศนิยมคือ 6, 7 หรือ 8 หลักข้อมูลนั้นอาจถูกคาดการณ์ไปยังระบบ State Plane หรือระบบพิกัด UTM

คำถามที่พบบ่อย: ข้อมูลเบื้องต้นเกี่ยวกับการฉายภาพ: สิ่งที่ GIS ต้องการทราบ

หน้านี้ให้การตรวจสอบหลายครั้งเกี่ยวกับการพิจารณาการคาดการณ์ซึ่งมีความสำคัญในการช่วยเหลือในการพยายามระบุระบบพิกัดที่ไม่รู้จัก

แนวคิดต่อไปนี้เป็นพื้นฐานในการทำความเข้าใจการใช้แผนที่ใน ArcGIS โปรดทราบว่าหัวข้อของการคาดการณ์นั้นกว้างมากและบทความนี้ไม่สามารถสัมผัสได้ในหัวข้อสำคัญสองสามข้อ

1. ระบบพิกัดที่รู้จักกันในชื่อการคาดการณ์แผนที่เป็นข้อมูลที่กำหนดโดยพลการสำหรับข้อมูลเชิงพื้นที่ วัตถุประสงค์ของพวกเขาคือการให้พื้นฐานทั่วไปสำหรับการสื่อสารเกี่ยวกับสถานที่หรือพื้นที่บนพื้นผิวโลก ปัญหาที่สำคัญที่สุดในการจัดการกับการคาดการณ์แผนที่คือการรู้ว่าการฉายคืออะไรและมีข้อมูลระบบพิกัดที่ถูกต้องที่เกี่ยวข้องกับชุดข้อมูล

2. เมื่อการคาดคะเนแผนที่ครั้งแรกถูกคิดขึ้นมามันไม่ถูกต้องว่าโลกนั้นแบนราบ ต่อมามีการแก้ไขข้อสมมติฐานและโลกก็ถือว่าเป็นทรงกลมที่สมบูรณ์แบบ ในศตวรรษที่ 18 ผู้คนเริ่มตระหนักว่าโลกไม่ได้กลมอย่างสมบูรณ์ นี่คือจุดเริ่มต้นของแนวคิดของการทำแผนที่รูปทรงกลม

3. เพื่อให้ได้ตำแหน่งที่แม่นยำยิ่งขึ้นบนพื้นผิวโลกผู้ทำแผนที่ได้ศึกษารูปร่างของโลก (geodesy) และสร้างแนวคิดของทรงกลม จากนั้นจึงทำการวางแผนระบบพิกัดทางภูมิศาสตร์ (GCS) ซึ่งรวมถึงข้อมูลหน่วยวัดและเส้นลมปราณหลัก ตัวเลขเชื่อมโยงทรงกลมกับส่วนหนึ่งส่วนใดของพื้นผิวโลก datums ล่าสุดได้รับการออกแบบเพื่อให้พอดีกับพื้นผิวโลกทั้งโลก

4. datums ที่ใช้บ่อยที่สุดในอเมริกาเหนือคือ:

• NAD 1927 (อเมริกาเหนือ Datum 1927) โดยใช้ Clarke 1866 Spheroid
• NAD 1983 (อเมริกาเหนือ Datum 1983) โดยใช้ GRS 1980 Spheroid
• WGS 1984 (World Geodetic Survey 1984) โดยใช้ WGS 1984 Spheroid

spheroids รุ่นใหม่ถูกพัฒนาขึ้นจากการตรวจวัดด้วยดาวเทียมและมีความแม่นยำมากกว่าการพัฒนาโดยคล๊าร์คในปี 1866 คำว่า 'ระบบพิกัดทางภูมิศาสตร์' และ 'ข้อมูล' ถูกใช้แทนกันได้ แต่ดังที่ได้กล่าวไว้ข้างต้น และเส้นลมปราณสำคัญ

5. พิกัดของข้อมูลจะเปลี่ยนไปขึ้นอยู่กับ datum และ spheroid ที่พิกัดเหล่านั้นใช้แม้ว่าจะใช้การฉายแผนที่และพารามิเตอร์แบบเดียวกัน

ตัวอย่างเช่นพิกัดทางภูมิศาสตร์ด้านล่างมีจุดเดียวที่ตั้งอยู่ในเมือง Bellingham รัฐวอชิงตันโดยใช้ข้อมูล 3 ชุดที่แตกต่างกัน:

รหัส: DATUM X-Coordinate Y-Coordinate NAD_1927 -122.466903686523 48.7440490722656 NAD_1983 -122.46818353793 48.7438798543649 WGS_1984 -122.46818353793 48.74387934299

6. หลักการของการจัดการข้อมูลที่ดีคือการได้รับพารามิเตอร์การฉายจากแหล่งข้อมูลที่ให้ข้อมูล อย่าคาดเดาเกี่ยวกับการคาดการณ์ข้อมูลเพราะมีฐานข้อมูล GIS ที่ไม่ถูกต้องจะเป็นผลลัพธ์ พารามิเตอร์ที่จำเป็นมีดังต่อไปนี้:

•การฉาย
•หน่วยการวัด
•โซน (สำหรับ UTM)
•โซน FIPS (สำหรับระนาบสถานะ)
•ฐานข้อมูล

อาจจำเป็นต้องใช้พารามิเตอร์อื่น ๆ ทั้งนี้ขึ้นอยู่กับการฉาย ตัวอย่างเช่นการคาดการณ์ของ Albers และ Lambert ต้องการพารามิเตอร์ต่อไปนี้:

•มาตรฐานขนานครั้งที่ 1 ในหน่วยองศานาทีและวินาที (DMS)
•ขนานมาตรฐานครั้งที่ 2 (DMS)
•เซ็นทรัลเที่ยง (DMS)
•ละติจูดของแหล่งกำเนิดการคาดการณ์ (DMS)
•ทิศตะวันออกเท็จและหน่วยวัด
•ทิศเหนือเท็จและหน่วยวัด
• X-shift และหน่วยวัด
• Y-shift และหน่วยวัด

ฉันมีปัญหาเดียวกันในขณะที่ฉันกำลังทำงานในโครงการของฉัน (DEM) หนึ่งใน DEM เป็นพิกัดที่ไม่รู้จักและไม่สามารถคาดการณ์ดังนั้นฉันได้ 'อ้างอิงทางภูมิศาสตร์' ใช้เครื่องมืออ้างอิงทางภูมิศาสตร์ในArcMap หลังจากการอ้างอิงทางภูมิศาสตร์คุณจะได้รับ prj.file ของคุณ

ดูเว็บไซต์นี้:

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

ฉันเป็นผู้พัฒนาเครื่องมือฟรีนี้ บางทีเครื่องมือนี้อาจช่วยให้คุณค้นหาและกำหนดการฉายภาพที่ถูกต้องของไฟล์รูปร่างของคุณ มันขึ้นอยู่กับความคิดของhttp://projfinder.com/โดย Aaron Racicot

คุณสามารถลองใช้โซลูชันพิเศษนี้สำหรับ:

ออสเตรีย (เยอรมัน): https://www.data.gv.at/anwendungen/checkaustrianproject/

ออสเตรเลีย: https://maegger.github.io/map_australia.html

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

คุณสามารถแยกแยะระหว่างพิกัดทางภูมิศาสตร์และการทำแผนที่พิกัดละติจูด / ลองจิจูดครั้งแรกอยู่ในระดับที่ภาคเหนือ / ตะวันออกที่สองอยู่ในหน่วยเมตรหากไฟล์มีรูปร่างคุณสามารถทิ้งพิกัดของรูปทรงเรขาคณิตด้วย shapelib-1.3.0 และดูค่าจากนั้น ..

คุณสามารถทำได้ด้วยฟิโอน่า

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

คุณจะได้คำตอบว่า

    {'init': 'epsg:4326'}