ระบบสารสนเทศทางภูมิศาสตร์

บางคนสามารถช่วยอธิบายและอ้างอิงบางอย่างเพื่อทำความเข้าใจกับลักษณะของพัลส์คืนภายในระบบ LiDAR

13 remote-sensing  lidar  terminology 

หากฉันต้องการฉายข้อมูลละติจูด & ลองจิจูดสำหรับประเทศสหรัฐอเมริกาที่อยู่ติดกัน (สหรัฐอเมริกาไม่รวมอลาสก้าและฮาวาย) ฉันจะใช้การคาดการณ์แบบใด ฉันชอบระยะทางที่แม่นยำยิ่งขึ้นตามด้วยรูปร่าง

13 coordinate-system  united-states 

ฉันมี Apache ที่ใช้ Geoserver ในอินสแตนซ์เสมือนจริงของฉันและฉันสงสัยว่าไดเรกทอรีข้อมูล Geoserver สามารถชี้ไปที่ที่ฝากข้อมูล Amazon S3 ได้หรือไม่ เป็นไปได้ไหม สิ่งนี้มีผลต่อประสิทธิภาพอย่างไร หาก S3 เป็นการตั้งค่าที่ไม่ดีสำหรับไดเรกทอรีข้อมูลของ Geoserver มีคำแนะนำสำหรับการเพิ่มข้อมูลจำนวนมากในไดเรกทอรีข้อมูล Geoserver หรือไม่? ไดเรกทอรีข้อมูลอาจสูงถึง 100 จาก GB หรืออาจมากกว่านั้น ฉันควรรับเซิร์ฟเวอร์ที่มีฮาร์ดไดรฟ์ขนาดใหญ่หรือเชื่อมต่อกับโซลูชันที่ปรับขนาดได้มากขึ้นเช่นที่ฝากข้อมูล Amazon S3 หรือไม่

13 geoserver  amazon-s3 

ปิด คำถามนี้จะต้องมีมากขึ้นมุ่งเน้น ไม่ยอมรับคำตอบในขณะนี้ ต้องการปรับปรุงคำถามนี้หรือไม่ อัปเดตคำถามเพื่อให้มุ่งเน้นที่ปัญหาเดียวโดยแก้ไขโพสต์นี้ ปิดให้บริการใน2 ปีที่ผ่านมา ฉันมีรูปร่างถนนขนาดใหญ่ (~ 70MB) และต้องการแปลงเป็นแรสเตอร์ด้วยความหนาแน่นของถนนในแต่ละเซลล์ เป็นการดีที่ฉันต้องการทำสิ่งนี้ใน R พร้อมกับเครื่องมือบรรทัดคำสั่งของ GDAL หากจำเป็น วิธีการเริ่มต้นของฉันคือการคำนวณความยาวของส่วนของเส้นตรงในแต่ละเซลล์ตามเธรดนี้โดยตรง สิ่งนี้ให้ผลลัพธ์ที่ต้องการ แต่ค่อนข้างช้าแม้สำหรับรูปร่างที่เล็กกว่าของฉันมาก นี่คือตัวอย่างที่ง่ายมากซึ่งค่าเซลล์ที่ถูกต้องชัดเจน: require(sp) require(raster) require(rgeos) require(RColorBrewer) # Create some sample lines l1 <- Lines(Line(cbind(c(0,1),c(.25,0.25))), ID="a") l2 <- Lines(Line(cbind(c(0.25,0.25),c(0,1))), ID="b") sl <- SpatialLines(list(l1,l2)) # Function to calculate lengths of lines in given raster …

13 raster  r  road  kernel-density  spatstat 

ฉันเจอปัญหาแปลก ๆ ในขณะที่เรียกใช้สคริปต์ประมวลผล PyQGIS แบบสแตนด์อโลน ลำดับการนำเข้าในสคริปต์มีผลต่อการดำเนินการตามปกติ คุณสามารถทำให้เกิดปัญหาได้อีกครั้งโดยเปิดคอนโซล Python และป้อนสคริปต์ต่อไปนี้ (ฉันใช้ GNU / Linux, QGIS 2.6.1, กำลังประมวลผลปลั๊กอิน v.2.2.0-2 และ Python 2.7.3): # Prepare the environment import sys from qgis.core import QgsApplication from PyQt4.QtGui import QApplication app = QApplication([]) QgsApplication.setPrefixPath("/usr", True) QgsApplication.initQgis() # Prepare processing framework sys.path.append('/home/YOUR_USER/.qgis2/python/plugins') from processing.core.Processing import Processing Processing.initialize() …

13 qgis  python  pyqgis  import  qgis-processing 

ฉันมีข้อมูลแบบแรสเตอร์ในตาราง postgres ที่ฉันต้องการใช้เป็น python ในรูปแบบ numpy ฉันใช้ psycopg2 เพื่อเชื่อมต่อกับฐานข้อมูล ฉันสามารถดาวน์โหลดข้อมูล แต่กลับมาเป็นสตริง (อาจเป็นเลขฐานสองต่อเนื่อง) ไม่มีใครรู้วิธีการใช้สายนี้และแปลงเป็นอาร์เรย์ numpy? ฉันสำรวจตัวเลือกอื่น ๆ เพื่อดาวน์โหลดแรสเตอร์เช่นใช้ st_astiff และเข้ารหัสเพื่อดาวน์โหลดไฟล์ hex และใช้ xxd แต่มันไม่ทำงาน ฉันยังคงได้รับข้อผิดพลาด 'rt_raster_to_gdal: ไม่สามารถโหลดไดรเวอร์ GDAL เอาต์พุต' และฉันไม่ได้รับอนุญาตให้ตั้งค่าตัวแปรสภาพแวดล้อมเพื่อให้สามารถเปิดไดรเวอร์ได้ TL, DR: ต้องการนำเข้าข้อมูลแรสเตอร์เข้าสู่อาร์เรย์ numpy (โดยใช้ python)

13 postgis  gdal  numpy  psycopg2  postgis-raster 

ฉันเพิ่งติดตั้ง QGIS 2.6 และติดตั้งปลั๊กอิน OpenLayers แล้ว มันแสดงรายการปลั๊กอิน OpenLayers ในแท็บที่ติดตั้ง แต่ฉันไม่สามารถเพิ่มเลเยอร์ให้กับโครงการของฉันได้เนื่องจากไม่ได้ให้ตัวเลือกในการเพิ่มเลเยอร์ OpenLayers เหมือนเวอร์ชันก่อนหน้า ไม่มีใครรู้ว่าเกิดอะไรขึ้นที่นี่?

13 qgis-openlayers-plugin 

ฉันกำลังมองหาวิธีการประมวลผลภาพระยะไกลและแยกพื้นที่มงกุฎของต้นไม้แต่ละต้นออกจากภาพ ฉันมีทั้งภาพที่มีความยาวคลื่นเป็นภาพและมีข้อมูล Lidar จากพื้นที่ สถานที่ที่สงสัยคือพื้นที่ทะเลทรายดังนั้นต้นไม้จึงไม่หนาแน่นเท่ากับพื้นที่ป่า ความละเอียดของภาพถ่ายทางอากาศคือ 0.5 ฟุตคูณ 0.5 ฟุต ความละเอียด LIDAR อยู่ที่ประมาณ 1 x 1 ฟุต ทั้งข้อมูลภาพและ lidar มาจากชุดข้อมูล Pima County, Arizona ตัวอย่างประเภทภาพถ่ายทางอากาศที่ฉันมีอยู่ที่ส่วนท้ายของโพสต์นี้ คำถามนี้การตรวจสอบต้นไม้แบบเดี่ยวใน ArcMap? ดูเหมือนจะเป็นปัญหาเดียวกัน แต่ดูเหมือนจะไม่มีคำตอบที่ดี ฉันสามารถรับการจำแนกประเภทพืชพรรณที่สมเหตุสมผล (และข้อมูลเกี่ยวกับการครอบคลุมเปอร์เซ็นต์โดยรวม) ในพื้นที่โดยใช้การจัดหมวดหมู่ Iso Cluster ใน Arcmap แต่จะให้ข้อมูลเล็กน้อยเกี่ยวกับต้นไม้แต่ละต้น ที่ใกล้เคียงที่สุดกับสิ่งที่ฉันต้องการคือผลลัพธ์ของการส่งผ่านการจัดหมวดหมู่ isocluster ผ่านคุณสมบัติ Raster to Polygon ใน Arcmap ปัญหาคือว่าวิธีนี้ผสานต้นไม้เข้าด้วยกันเป็นรูปหลายเหลี่ยมเดียว แก้ไข: ฉันควรจะมีรายละเอียดเพิ่มเติมเกี่ยวกับสิ่งที่ฉันมี ชุดข้อมูลดิบที่ฉันมีคือ: ข้อมูลเลเซอร์เต็มรูปแบบและรูปแบบแรสเตอร์ที่สร้างขึ้นจากมัน ภาพที่มองเห็น …

13 remote-sensing  lidar  classification  land-cover  image-segmentation 

(แก้ไข 2019) ST_VoronoiPolygonsมีให้ตั้งแต่ PostGIS v2.3 ! ด้วย PostGIS 2.1+ เราสามารถใช้ST_DelaunayTriangles ()เพื่อสร้างDelaunay triangulationซึ่งเป็นกราฟคู่ของแผนภาพ Voronoiและในทางทฤษฎีพวกเขามีการแปลงที่แน่นอนและย้อนกลับได้ มีสคริปต์มาตรฐาน SQL ที่ปลอดภัย พร้อมอัลกอริธึมที่เหมาะสมสำหรับการแปลง PostGIS2 Delaunay-to-Voronoi นี้หรือไม่? อ้างอิงอื่น ๆ : 1 , 2

13 postgis  sql  voronoi-thiessen 

วิธีที่ดีที่สุดในการรับชื่อ OBJECTID ของคลาสคุณลักษณะคืออะไร ฉันรู้ว่าคุณสามารถใช้ "OID @" ในเคอร์เซอร์เพื่อรับฟิลด์ แต่คุณสามารถรับชื่อจากที่นั่นได้หรือไม่? ฉันกำลังเขียนสคริปต์และต้องคำนึงถึงชื่อ OBJECTID ที่ต่างกันเช่น OBJECTID_2, OBJECTID_12, FID เป็นต้นซึ่งจะใช้ในการสืบค้น ใช้คำอธิบายเป็นวิธีเดียวที่จะได้ชื่อมาหรือไม่ ตัวอย่างเช่น: oid_field = [field.name for field in arcpy.ListFields("FC") if field.type == "OID"] ฉันแค่สงสัยว่ามีวิธีอื่น ขอบคุณ

13 arcpy  objectid 

ฉันอยู่ในโครงการที่ต้องการรวมข้อมูลสำหรับทุกรัฐในสหรัฐอเมริกา หมายความว่าฉันมีชุดข้อมูลถนนสองชุดและฉันจะกำหนดคุณลักษณะของถนนในชุดข้อมูลหนึ่งให้กับถนนสายเดียวกันในชุดข้อมูลอื่น (คล้ายกับการรวมเชิงพื้นที่) แต่ละรัฐเป็น shapefile อิสระและฉันต้องรวมพวกเขาทั้งหมดในตอนท้ายเพื่อสร้างผลิตภัณฑ์ขั้นสุดท้าย ฉันสงสัยว่าระบบพิกัดที่ฉันควรใช้ในกรณีนี้คืออะไร? หรือไม่เป็นไรเพราะมันจะไม่ส่งผลต่อผลลัพธ์ของกระบวนการทางภูมิศาสตร์? สำหรับการทำงานที่ผ่านมาของฉันฉันทำแผนที่พื้นที่ขนาดเล็กเท่านั้นดังนั้นฉันจึงเลือกเครื่องบินของรัฐ NAD 83 ของพื้นที่เฉพาะเป็นระบบที่คาดการณ์ไว้

13 coordinate-system  united-states 

ฉันมีรูปร่างสองไฟล์ที่แทนกริดที่มีค่าที่เป็นสัญลักษณ์ด้วยสีที่สำเร็จการศึกษาดังแสดงด้านล่าง: คุณสามารถเห็นได้ว่ารูปร่างสองไฟล์นั้นรวมช่วงของข้อมูลที่คล้ายกัน แต่กำหนดสีที่แตกต่างกันเล็กน้อยให้กับช่วงในข้อมูลเหล่านี้ ฉันต้องการเป็นสัญลักษณ์ของช่วงในแต่ละไฟล์ทั้งสองโดยใช้สีเดียวกันสำหรับช่วงเดียวกันทำให้การเปรียบเทียบง่ายขึ้นระหว่างไฟล์และอนุญาตให้ใช้คำอธิบายภาพเดี่ยว อย่างไรก็ตามเมื่อฉันพยายามที่จะจำแนกช่วงที่มีช่วงเวลาด้วยตนเอง ArcGIS บังคับให้สิ่งที่ควรเป็นช่วงบนสุดที่จะรวมถึงค่าต่ำสุดจากชุดข้อมูล คุณสามารถดูได้ในช่วง "-81.64 - 10.00" ของกรอบข้อมูลด้านซ้ายด้านล่าง สิ่งนี้บังคับให้เลเยอร์ทั้งหมดถูกทำสัญลักษณ์ด้วยสีนี้ โดยพื้นฐานแล้ว ArcGIS ต้องการให้ใช้ค่าต่ำสุดในชุดข้อมูลเป็นค่าในช่วงใดช่วงหนึ่ง มีวิธีหลีกเลี่ยงคุณลักษณะนี้หรือไม่?

13 arcgis-desktop  arcgis-10.1  cartography  symbology  legend 

ฉันใหม่สำหรับ R และใช้แพ็คเกจแรสเตอร์ ฉันมีปัญหาในการแยกรูปหลายเหลี่ยมจากไฟล์แรสเตอร์ที่มีอยู่ ถ้าฉันใช้ extract(raster, poly_shape) ฟังก์ชั่นบนแรสเตอร์มันมักจะสร้างรายการที่มีข้อมูล สิ่งที่ฉันต้องการคือการแยกไฟล์แรสเตอร์อื่นที่ฉันสามารถโหลดด้วย ArcGIS อีกครั้ง หลังจากอ่านอีกเล็กน้อยฉันคิดว่าฟังก์ชั่นการครอบตัดคือสิ่งที่ฉันต้องการจริงๆ แต่เมื่อลองใช้ฟังก์ชั่นนี้ crop(raster, poly_shape) ฉันได้รับข้อผิดพลาดนี้: Error in .local(x, y, ...) : extents do not overlap In addition: Warning message: In intersect(extent(x), extent(y)) : Objects do not overlap ไฟล์ raster และ poly_shape เหมือนกันทั้งสองฟังก์ชั่น คุณช่วยบอกฉันทีว่าอะไรที่นี่ผิด มันถูกต้องไหมที่ฟังก์ชั่นการครอบตัดสร้างอีกแบบแรสเตอร์และไม่ใช่รายการ? แก้ไข : ฟังก์ชั่นขอบเขต () ไม่ทำงานสำหรับฉัน …

13 raster  shapefile  r  line 

ฉันกำลังทำงานในโครงการ GIS ขนาดใหญ่และการไม่สามารถจัดเรียงรายการในแผงควบคุมเลเยอร์ทำให้ฉันลำบาก ฉันสามารถมีได้ถึง 250 ชั้นในแผงเลเยอร์ ฉันมีรูปแบบการตั้งชื่อพิเศษสำหรับไฟล์แรสเตอร์ TIFF ทั้งหมดที่ฉันนำเข้าและอ้างอิงทางภูมิศาสตร์ใน QGIS แต่เมื่อฉันนำเข้าฉันไม่สามารถนำเข้าทีละตัวตามลำดับตัวอักษรเนื่องจากเพื่อนร่วมงานของฉันกำลังทำงานกับไฟล์บางไฟล์และฉันกำลังทำงานกับคนอื่น ดังนั้นรายการในพาเนลเลเยอร์จึงอยู่ทั่วสถานที่ โปรแกรมไม่ฉลาดพอที่จะเรียงลำดับตามตัวอักษร แต่จะเรียงลำดับพวกเขาตามลำดับที่ฉันนำเข้า ที่นี่ลองดูภาพหน้าจอนี้ ไม่มี QGIS ในตัวที่จะทำให้ฉันสามารถจัดเรียงรายการเหล่านี้ได้หรือไม่? ถ้าไม่มีมีปลั๊กอินที่ฉันสามารถติดตั้งได้หรือไม่ แก้ไข: ฉันพบสิ่งนี้และโค้ดสองสามบรรทัดแรกดูเหมือนจะหลอกลวง from PyQt4.QtCore import * from PyQt4.QtGui import * from qgis.utils import iface mw = iface.mainWindow() lgd = mw.findChild(QTreeWidget, "theMapLegend") # get ref to object by type/objectName lgd.sortItems(0, Qt.AscendingOrder) # …

13 qgis 

ฉันมีรูปหลายเหลี่ยมสองรูปแบบ: รูปหลายเหลี่ยม 1 และรูปหลายเหลี่ยม 2 ฉันต้องการแสดงปริมาณที่โพลีกอน 1 มีสองตัวชี้วัดพื้นที่และความยาวของเส้นรอบวงที่โพลีกอน 1 มีเส้นรอบวง / ขรุขระ / ผิดปกติมากกว่าโพลีกอน 2 รูปหลายเหลี่ยมแต่ละอันมีความยาวเส้นรอบวงเท่ากัน แต่แต่ละรูปครอบคลุมพื้นที่ที่แตกต่างกันมาก ในการหาปริมาณความไม่สม่ำเสมอ / ความขรุขระ / ความผิดปกติของรูปหลายเหลี่ยมแต่ละรูปควรคำนวณดังนี้ area/perimeter หรือ perimeter/area ฉันคิดว่าperimeter/areaแต่แล้วฉันก็พบโพสต์บล็อกนี้ซึ่งใช้area/perimeter: http://www.r-bloggers.com/measuring-the-gerrymander-with-spatstat/

13 spatial-statistics