คำถามติดแท็ก coordinate-system

กรอบอ้างอิงประกอบด้วยชุดของจุดเส้นและ / หรือพื้นผิวและชุดของกฎใช้เพื่อกำหนดตำแหน่งของจุดในอวกาศในสองหรือสามมิติ แท็กนี้ยังหมายถึงการคาดคะเนแบบทำแผนที่และการแปลงระบบพิกัด

1
ระบบพิกัดของฉันคืออะไร
ฉันพยายามที่จะเข้าใจระบบพิกัดข้อมูลตัวเลขประมาณการและอื่น ๆ ได้ดีขึ้นฉันเข้าใจว่าบางระบบรักษาสิ่งต่าง ๆ ไว้ ระยะทางรูปร่างมุมพื้นที่ ฉันสงสัยว่าคุณจะบอกได้อย่างไรว่าระบบบางอย่างกำลังคงอยู่ ส่วนใดของชื่อ CRS ที่ระบุสิ่งนี้ โดยเฉพาะฉันสงสัยเกี่ยวกับระบบ UTM ในอัลเบอร์ตาแคนาดา NAD83 (CSRS) / อัลเบอร์ตา 10-TM (ป่าไม้)

3
การคาดการณ์ของ Sentinel-2
ฉันเพิ่งสังเกตเห็นการเปลี่ยนแปลงในประมาณการทางภูมิศาสตร์ของข้อมูล Sentinel 2 โปรดดูภาพหน้าจอด้านล่างสำหรับการชี้แจง (Granule T33UXP, แบนด์ 8 เท่านั้น) ภาพแรกถูกถ่ายใน 20160314 ภาพที่สอง (เลื่อน) ใน 20160403 ฉันเพิ่มรูปหลายเหลี่ยมเป็นข้อมูลอ้างอิง นี่เป็นข้อมูลที่ยังไม่ได้ประมวลผลซึ่งดาวน์โหลดจากฮับ API S2A_OPER_MSI_L1C_TL_SGS__20160314T171554_A003794_T33UXP_N02.01 S2A_OPER_MSI_L1C_TL_SGS__20160403T152610_A004080_T33UXP_N02.01 anyobdy สังเกตเห็นสิ่งนี้หรือไม่? นี่เป็นปัญหาที่ทราบหรือเกิดขึ้นเฉพาะในกรณีของฉันหรือไม่

1
กำลังล้างข้อมูลระบบพิกัดอ้างอิงที่ใช้ล่าสุดโดยใช้ QGIS
หลังจากใช้ SRC ที่ไม่ปกติเพื่อวัตถุประสงค์ในการสาธิตฉันต้องการล้างรายการ "ระบบพิกัดที่ใช้ล่าสุดของฉัน" เพื่อประโยชน์ของ praticity เป็นไปได้ไหม? ฉันใช้ QGIS 2.10.1

1
ความแตกต่างระหว่าง gdalwarp และ projectRaster
ฉันพยายามฉายภาพแรสเตอร์ ใน R มีprojectRaster()ฟังก์ชั่นสำหรับสิ่งนี้ (ด้านล่างเป็นตัวอย่างที่ทำซ้ำได้อย่างสมบูรณ์): # example Raster require(raster) r <- raster(xmn=-110, xmx=-90, ymn=40, ymx=60, ncols=40, nrows=40) r <- setValues(r, 1:ncell(r)) projection(r) # project to newproj <- "+init=epsg:4714" # using raster package to reproject pr1 <- projectRaster(r, crs = CRS(newproj), method = 'bilinear') ซึ่งใช้งานได้ดี อย่างไรก็ตามมันค่อนข้างช้า เพื่อเพิ่มความเร็วฉันแม้ว่าจะใช้gdalwarpแทน (ด้วย SSD ค่าใช้จ่ายในการอ่านและเขียนจาก / …


1
ไฟล์ Georeferencing LAS ในพิกัดทางภูมิศาสตร์ (lat / lon) หรือใน cartesian (ที่คาดการณ์)?
ก่อนอื่นฉันเป็นมือใหม่ใน GIS ฉันมีคลาวด์จุด lidar กับ GPS และข้อมูล telemetry ที่เกี่ยวข้อง จากนั้นฉันเปลี่ยนพิกัดของทุกจุดให้เป็นระบบพิกัด ECEF และในที่สุดฉันก็เปลี่ยนพิกัดเชิงมุมเป็นระบบพิกัดคาร์ทีเซียนโดยใช้การประมาณ Local Tangent Plane (LTP) (การประมาณนี้ใช้พิกัด GPS ของจุดแรกเป็นจุดกำเนิดของระบบอ้างอิงจากนั้นอ้างอิงจุดที่เหลือตามลำดับแรก) ฉันกำลังเขียนไฟล์ LAS โดยใช้ liblas และฉันไม่รู้ว่าจะตั้งค่าข้อมูลอ้างอิงทางภูมิศาสตร์เป็นไฟล์ LAS ได้อย่างไร ฉันกำลังเขียนพิกัดของทุกจุดในการประมาณ LTP ดังนั้นฉันจึงมี 2 คำถาม: มันถูกต้องหรือไม่ที่จะเขียนพิกัดจุดที่เกี่ยวข้องกับการประมาณ LTP หรือฉันควรจะเขียนพิกัดในหน่วย lat / lon (gps)? หากสิ่งนี้ถูกต้องฉันจะตั้งค่าการอ้างอิงทางภูมิศาสตร์ให้สอดคล้องกับ LTP ได้อย่างไร หาก '1' ไม่ถูกต้องและฉันต้องเขียนพิกัดในหน่วย lat / lon ฉันจะตั้งค่าการอ้างอิงทางภูมิศาสตร์ได้อย่างไร การเขียน VLR …

4
การคาดการณ์ ESRI WKT แตกต่างจากการคาดการณ์ของ OGC WKT อย่างไร
ไม่มีใครทราบรายการความแตกต่างที่แม่นยำระหว่างสตริงรูปแบบการฉาย ESRI WKT และ OGC WKT หรือไม่ ฉันรู้ว่ามีเครื่องมือหลายอย่างที่ช่วยแปลงจาก ESRI WKT เป็น OGC WKT รวมถึงอรรถประโยชน์ของ GDAL และบริการเว็บไซต์ต่างๆ แต่คำถามของฉันไม่ใช่แบบที่ใช้งานได้จริงฉันเพียงต้องการเข้าใจความแตกต่างของการจัดรูปแบบ / ไวยากรณ์ที่บริการเหล่านี้ใช้งาน คำถาม Stackexchange ก่อนหน้านี้ได้พูดคุยเกี่ยวกับความแตกต่างในตัวอย่างเฉพาะเท่านั้นหรือเกี่ยวกับเครื่องมือและบริการที่มีให้เท่านั้น แม้ว่าคุณจะเพิ่งทราบความแตกต่างมันก็คงจะดีถ้าคุณเพิ่งโพสต์มัน จากประสบการณ์ของฉันเองควรมีความแตกต่างเล็กน้อย ความแตกต่างที่ฉันรู้คือ: รายการข้อความส่วนใหญ่ในการกำหนด esri ใช้ขีดล่างที่ ogc ใช้พื้นที่ ข้อความที่กำหนด datum ใน esri wkt เหมือนกับ ogc wkt ยกเว้นว่าเริ่มต้นด้วย "D_" บางครั้งตัวระบุข้อความสำหรับ PROJCS ที่กำหนดไว้ล่วงหน้า, PROJECTION, GEOGCS และ DATUM จะถูกเขียนต่างกัน (เช่น "NAD83" …

1
แปลงพิกัดจาก readShapePoly ใน R เป็น long-lat พิกัด
ฉันกำลังพยายามหาเซนทรอยด์บางส่วนของไฟล์รูปร่างสำหรับเทศบาลเดนมาร์กและหลังจากนั้นหาเวลาขับรถระหว่างพวกเขา ผมใช้RของreadShapePolyฟังก์ชั่นจากmaptoolsรวมกับgCentroidฟังก์ชั่นจากrgeosและผลงานทุกอย่าง อย่างไรก็ตามฉันได้รับ spatialpoints เช่น SpatialPoints: x y 1 571860.7 6225016 Coordinate Reference System (CRS) arguments: NA ซึ่งไม่ใช่สิ่งที่ชัดเจนที่ฉันสามารถใช้ใน Google เพื่อคว้าเวลาเดินทาง ฉันกำลังมองหาวิธีการแปลงตัวเลขเหล่านี้เป็นลองจิจูดละติจูด แต่ไม่มีความคิดวิธี เมื่อฉันอ่านข้อมูลโดยใช้readOGRจากrgdalห้องสมุดฉันได้รับพิกัดเดียวกัน แต่มันบอกฉันเกี่ยวกับสิ่งที่ฉันถือว่าเป็นเส้นโครง (แต่พิกัดเหมือนกัน) Slot "proj4string": CRS arguments: +proj=utm +zone=32 +ellps=intl +units=m +no_defs ตัวอย่างที่ทำซ้ำได้: ฉันใส่ข้อมูลสำหรับตัวอย่างที่นี่: https://github.com/sebastianbarfort/shapefiles สิ่งนี้ควรทำให้เกิดปัญหาซ้ำ: library(maptools) library(rgdal) library(rgeos) map = readShapePoly("~/Downloads/shapefiles-master/kommuner1983.shp") centroid = gCentroid(map) centroid

1
พล็อตแบบอคติของ Mercator ในแผ่นพับ
ฉันได้รับแรงบันดาลใจจากภาพวิกิพีเดียนี้จากหน้าโปรเจคเตอเรเตอร์เพื่อแฮกใบปลิว แนวคิดก็คือให้เคอร์เซอร์ของคุณเคลื่อนที่ไปตามทวีปและจะได้รับการแก้ไขสำหรับการฉายภาพ ดู JSFiddle ต่อไปนี้สำหรับสถานะ sofar; http://jsfiddle.net/cantdutchthis/NrAEL/1/ ฉันเป็นคนใหม่ในโลกฉายและถึงแม้ว่าตัวอย่างจะใช้ได้ดีพอสำหรับละติจูด เมื่อฉันเปรียบเทียบequirectangularและMercator เงื้อมมือบนวิกิพีเดียฉันได้รับความประทับใจว่าอคติของเส้นโครงนั้นอยู่ที่ 'แกน y' เมื่อคุณเปรียบเทียบการฉายสองภาพนี้ ฉันได้รับความประทับใจว่าฉันไม่ได้แก้ไขลองจิจูดเนื่องจากรูปภาพแรกของวิกิพีเดีย ฉันต้องการสนับสนุนสิ่งนี้ แต่ไม่รู้ว่าจะทำอย่างไรดีที่สุด ในปัจจุบันวิธีที่ฉันสามารถแสดงอคติของการฉายภาพของ Mercator คือการเปลี่ยนพิกัดไปที่ผู้ชี้ขาดทุกครั้งที่เมาส์เคลื่อนที่ map.on("mousemove", function(d){ Object.keys(myLayer._layers).map(function(d){ myLayer.removeLayer( myLayer._layers[d] ) }) var plotCoords = []; // green is the coordinate array for the landmass var green = [ [-46.76373672, 82.62796021],... ] for ( i in …

4
ใช้ไลบรารี PROJ.4 เพื่อเปลี่ยนจากระบบพิกัดในท้องถิ่นเป็นระบบพิกัดทั่วโลกโดยใช้จุดควบคุมภาคพื้นดินหรือไม่?
ฉันมีคลาวด์พอยท์ซึ่งมีพิกัดที่เกี่ยวกับระบบพิกัดในพื้นที่ ฉันยังมีจุดควบคุมภาคพื้นดินพร้อมค่า GPS ฉันสามารถแปลงพิกัดท้องถิ่นเหล่านี้เป็นระบบพิกัดทั่วโลกโดยใช้ PROJ.4 หรือห้องสมุดอื่น ๆ ได้หรือไม่? รหัสใด ๆ ใน Python สำหรับปัญหาที่ระบุข้างต้นจะเป็นความช่วยเหลือที่ดี

3
การวัดระยะทางข้ามโซน UTM: ใช้วิธีทางภูมิศาสตร์หรือภาพถ่าย
ฉันมีตารางสำรวจกระจายทั่วทั้ง 3 โซน UTM (36N, 36S, 37S) ฉันต้องการค้นหาระยะทางที่ใกล้ที่สุด (หรือสั้นที่สุด) ของเซนทรอยด์ของกริดเหล่านี้ไปยังถนนและจุดต่าง ๆ ระหว่าง ดูเหมือนว่ามีการประนีประนอมมากเกินไปเมื่อใช้การฉายภาพถ่ายแบบใด ๆ (อ่าน: เกี่ยวกับการรักษาระยะห่างระหว่างจำนวนคะแนนใด ๆ บนแผนที่ ) เราควรลืมเกี่ยวกับการใช้การฉายภาพในกรณีนี้และใช้เทคนิคแบบ goedesic หรือ ellipsoidal (อ่าน: ภูมิศาสตร์) มีใครบ้างที่มีความรู้เกี่ยวกับเทคนิคระนาบที่จะรักษาระยะห่างระหว่างจำนวนคะแนนใด ๆ บนแผนที่? ดูเหมือนว่าฉันจะไม่สามารถใช้การฉายภาพระยะเท่ากันยกเว้นการฉายแบบ gnomonic ถูกต้องหรือไม่

2
จุดที่ใกล้ที่สุดบนเส้น (ประมาณการทรงกลม / Mercator)
ฉันมีเส้น (Axe, Ay - Bx, By)เหนือเส้นโครงแผนที่ปรอท (google maps) และจุดสุ่ม (Cx, Cy)ใกล้กับเส้นนั้นฉันอยากจะรู้จุดที่ใกล้ที่สุด (ภาพโปร่งใสสีน้ำเงิน) บรรทัดนั้นชี้ไปที่ (สีน้ำเงินในภาพ) แก้ไข: เพื่อชี้แจงว่าสิ่งนี้อยู่ในการฉายภาพ Mercator (การฉายทรงกลม)

3
การฉายแผนที่ลูกบาศก์ทรงกลมแบบสี่เหลี่ยมด้านเท่านั้นเหมือนกับการฉายแผนที่พื้นที่ลูกบาศก์ของสไนเดอร์หรือไม่?
การฉายแผนที่ทรงกลมลูกบาศก์สี่เหลี่ยมจตุรัส [1] เหมือนกับการฉายแผนที่พื้นที่ลูกบาศก์ของสไนเดอร์ [2] หรือไม่? ฉันรู้ว่าทั้งคู่เป็นพื้นที่ประมาณการ polyhedral แผนที่เท่ากัน ( http://www.progonos.com/furuti/MapProj/Dither/ProjPoly/projPoly.html ) ตามคิวบ์ แต่ฉันไม่สามารถเข้าถึงสิ่งพิมพ์เพื่อตอบคำถามด้วยตัวเอง . ขอบคุณ [1] EM O'Neill และ RE Laubscher ศึกษาเพิ่มเติมเกี่ยวกับฐานข้อมูลลูกบาศก์โลกทรงกลม quadrilateralized cube, Tech รายงาน 3-76, Computer Science Corp Silver Md System Sciences Div, 1976 [2] สไนเดอร์, JP 1992 การฉายแผนที่พื้นที่เท่ากันสำหรับลูกโลกโพลีฮีด Cartographica 29 (1): 10-21

2
วิธีการกำหนดพารามิเตอร์การฉายภาพเมื่อปรับแต่งการฉายภาพ
ฉันกำลังพยายามปรับแต่ง Albers และการฉายภาพ Hotine Oblique Mercator (HOM) เพื่อลดการบิดเบือนในภูมิภาคที่ฉันกำลังวิเคราะห์ ภูมิภาคนี้มีความกว้างตั้งแต่ 51 ถึง 62 องศาละติจูดครอบคลุมพื้นที่ประมาณขนาดของยูเครน ภูมิภาคนี้เป็นภูมิภาคที่มุ่งเน้น NW - SE ฉันต้องการที่จะให้แน่ใจว่าฉันกำลังใช้วิธีการที่ถูกต้องสำหรับการกำหนดพารามิเตอร์ที่สองฉาย: lat / ยาวศูนย์ฉายและที่ราบกลาง ฉันใช้ ArcMap v10 นี่คือขั้นตอนที่ฉันได้ปฏิบัติตาม: สร้างรูปหลายเหลี่ยมเดียวที่กำหนดขอบเขตการวิเคราะห์ (โดยทั่วไปแล้วสร้างเปลือกนูนรอบขอบเขตของแหล่งต้นน้ำที่ครอบคลุมพื้นที่) รูปหลายเหลี่ยมนี้เป็นพื้นที่ที่ฉันปรับแต่งการฉายภาพสำหรับ ฉายรูปหลายเหลี่ยมเป็น Geographic / NAD 83 ใช้เครื่องมือของ Jeff Jenness สำหรับกราฟิกและรูปทรง ( http://www.jennessent.com/arcgis/shapes_graphics.htm ) เพื่อกำหนดจุดศูนย์กลางมวลของรูปหลายเหลี่ยมบนทรงกลม GRS80 พิกัดที่ได้คือสิ่งที่ฉันใช้สำหรับพารามิเตอร์ "การฉายภาพกลาง" ในการพิจารณา centerline azimuth ฉันแรกคาดว่ารูปหลายเหลี่ยมกับการฉายภาพระยะเท่ากันแบบ azimuthal โดยระบุศูนย์กลางการฉายที่พิกัดที่กำหนดไว้ในขั้นตอนที่ …

1
ตัวอย่าง Tissot ellipse สำหรับการฉาย Equirectangular?
ฉันกำลังพยายามคำนวณการบิดเบี้ยวของการฉายภาพ Equirectangular ผ่าน Tissot Indicatrices ฉันได้ลองทำตามคำแนะนำในโพสต์นี้ (เหนือสิ่งอื่นใด) แต่มันเกินความเข้าใจของฉันในฐานะมือสมัครเล่น ดังนั้นฉันสงสัยว่าใครบางคนจะใจดีกับการคำนวณรูปวงรี Tissot เดียวสำหรับตัวอย่าง equirectangular lat / long (แล้วแต่จำนวนใดที่คุณชื่นชอบและบิดเบี้ยวในการฉายภาพแบบสองทิศทาง) ฉันไม่เข้าใจว่าตัวแปรคืออะไรและมาจากที่ใดดังนั้นการเห็นสมการในการดำเนินการจะมีประโยชน์มาก ฉันพยายามที่จะเข้าใจสมการเหล่านี้เพื่อที่ฉันจะได้เสียบเข้ากับโปรแกรมการแมปที่ฉันกำลังเขียนโปรแกรมอยู่ ฉันถามคำถามทั่วไปจำนวนมากในชุดข้อความนี้แต่ฉันคิดว่าตัวอย่างเฉพาะจะช่วยให้ฉันคิดออกส่วนที่เหลือ ขอบคุณมากเช่นเคย NCashew

โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.