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

ฉันพยายามกำหนดพิกัดที่ฉันได้รับ: N 39,21.440 W 120,48.076 แต่ฉันไม่คุ้นเคยกับรูปแบบ ฉันคุ้นเคยกับ DMS lat / long และรูปแบบมาตรฐาน (?) (เช่น 39.123234, -121.348348) ใครสามารถบอกฉันได้ว่าฉันสามารถแปลงเหล่านี้เป็น DMS หรือมาตรฐานได้อย่างไร

11 latitude-longitude  mapping 

ฉันเคยเห็นคำถามนี้ถามมาก่อน แต่คำตอบที่ฉันพยายามไม่ทำงาน ฉันกำลังค้นหาวิธี OSM ทั้งหมดที่อยู่ในกล่องขอบเขต ข้อมูล OSM ถูกนำเข้าด้วย Mercator ทรงกลมเริ่มต้น ฉันกำลังสอบถามโดยใช้ LAT / LON ดังนั้นการแปลง SELECT * FROM planet_osm_ways WHERE ST_Transform( ST_MakeEnvelope(37.808264, 37.804339, -122.275244, -122.271189, 4326),3857 ); เมื่อฉันเรียกใช้สิ่งนี้ฉันได้รับข้อผิดพลาดต่อไปนี้: ข้อผิดพลาด: อาร์กิวเมนต์ของ WHERE ต้องเป็นชนิดบูลีนไม่พิมพ์เรขาคณิตบรรทัดที่ 3: WHERE ST_Transform (ST_MakeEnvelope (37.808264, 37.804339, -12 ...

11 postgis  extents  osm2pgsql 

ฉันมี dataframe ของ geopandas ซึ่งเป็นพวงของ linestrings ที่มีข้อมูลบางอย่างที่เกี่ยวข้องกับจุดยอด / จุดแต่ละจุด: Point_x = (Lat, Lon, Time, ID, Data1, Data2, Data3) คะแนนจะถูกแปลงเป็น linestrings ตาม ID และเรียงตามเวลา ฉันต้องการทำลาย linestrings ตรงจุดที่ตรงกับเงื่อนไข ตอนนี้คือเมื่อระยะทางระหว่างคะแนนมากกว่าค่าบางอย่าง ในอนาคตอาจเป็นที่ที่ฟังก์ชันของเขตข้อมูลเป็นค่าบางอย่าง ตัวอย่างเช่นแบ่ง linestring เมื่อ Speed ​​ข้าม 5 kph ปัญหาในปัจจุบันคือบางแทร็กเกิดจากจุดที่มี ID ซ้ำกันดังนั้นการวิ่งย้อนกลับไปมาในระยะทางไกลและฉันต้องการให้ขีดแบ่งเส้นแบ่งเหล่านี้ ความคิดใด ๆ เกี่ยวกับวิธีที่ถูกต้องในการจัดโครงสร้างนี้หรือไลบรารี / วิธีการที่อาจใช้ dataframe มีแทร็กมากกว่า 150k ที่มีหลายจุดต่อแทร็คดังนั้นประสิทธิภาพจะดี นี่คือตัวอย่างของแทร็ก DF: ID …

11 python  shapely  geopandas 

ฉันกำลังมองหาการปรับขนาดภาพ TIFF ที่ไม่ได้ลงชื่อ 16 บิตเป็นภาพ 8 บิตโดยใช้ GDAL (OSGeo4W) อย่างไรก็ตามฉันต้องการ จำกัด ค่าพิกเซลเพื่อให้อยู่ในช่วงของภาพ 8 บิต สิ่งที่ฉันหมายถึงคือฉันต้องการให้แน่ใจว่าสัดส่วนภาพในภาพ 16 บิตนั้นแสดงในรูป 8 บิต (0-255) การใช้ gdal_translate เพื่อแปลงจาก 16 บิตเป็น 8 บิตจะตัดค่าพิกเซลและไม่สร้างภาพ RGB ฉันไม่แน่ใจในสิ่งที่ฟังก์ชั่นและการตั้งค่าที่จะโทรจากภายใน gdal_translate หรือถ้ามีตัวเลือกที่ดีกว่า GDALINFO C:\>gdalinfo C:\Projects\Stormwater\ForPCI\images\1537TORO_6cm_04761_RGBI.tif Driver: GTiff/GeoTIFF Files: C:\Projects\Stormwater\ForPCI\images\1537TORO_6cm_04761_RGBI.tif Size is 17310, 11310 Coordinate System is `' Metadata: TIFFTAG_DATETIME=2015:05:03 12:27:13 TIFFTAG_IMAGEDESCRIPTION=UltraCam-Lvl03 …

11 gdal  imagery  gdal-translate 

การใช้ QGIS 2.16.2 ฉันสามารถเลือกบันทึกจากตารางแอ็ตทริบิวต์ shapefile ได้สำเร็จโดยใช้นิพจน์ต่อไปนี้: left("start_time", 10) = '2015-08-01' โปรดทราบว่า start_time เป็นเขตข้อมูลสตริง อย่างไรก็ตามเมื่อฉันคัดลอกและวางนิพจน์เดียวกันลงในคุณสมบัติชั้นของ Shapefile> ทั่วไป> ตัวกรองคุณสมบัติ> ตัวสร้างแบบสอบถามฉันได้รับข้อผิดพลาดดังต่อไปนี้: ฉันคาดหวังว่าการแสดงออกเดียวกันจะทำงานได้ทุกที่ใน QGIS ผมทำอะไรผิดหรือเปล่า?

11 qgis  query  filter  expression 

มีเครื่องมือที่สนับสนุนการแยกส่วนเล็ก ๆ จากฟีด GTFS ที่ใหญ่กว่าหรือไม่? ไลบรารี Transitfeedของ Google มีเครื่องมือผสาน แต่ไม่มีเครื่องมือแยก / คลิปเท่าที่ฉันเห็น ห้องสมุดอื่น ๆ ที่ฉันค้นพบสามารถอ่านและเขียนพล็อตข้อมูลได้ แต่ไม่รองรับการปรับเปลี่ยน ข้อมูลความเป็นมา: Transitfeed Feedvalidator แสดงข้อผิดพลาดมากมายสำหรับไฟล์ GTFS ของฉัน ฉันต้องการแยกภูมิภาคที่สนใจเพื่อประเมินคุณภาพท้องถิ่น คำแนะนำใน Twitter ชี้ให้ฉันเห็น FME ซึ่งมีทั้งผู้อ่าน GTFS และนักเขียน ขั้นตอนการทำงานที่เรียบง่ายของฉันประกอบด้วยตัวอ่าน GTFS (ด้านล่างซ้าย), ตัวอ่าน Shapefile (มุมซ้ายบน, สำหรับพื้นที่รูปหลายเหลี่ยมที่น่าสนใจ), clipper และตัวเขียน GTFS สำหรับเอาท์พุท (เขียนเฉพาะจุดหยุดและไฟล์รูปร่าง จะต้องการไฟล์ทั้งหมด): ปัญหาที่เกิดขึ้นคือการส่งออกรูปร่างเนื้อหา (จุดสีส้ม) ไม่ได้มีจุดอินพุต (สีเขียว) ทั้งหมดซึ่งอยู่ในพื้นที่สีเหลืองที่น่าสนใจ: (แหล่งข้อมูล: https://transitfeeds.com/p/helsinki-regional-transport/735/latest/download …

11 fme  general-transit-feed-specification 

ฉันจะสร้างลักษณะของแผนที่ภูมิประเทศเก่าที่ใช้เส้นบาง ๆ เพื่อแสดงแหล่งน้ำได้อย่างไร ฉันดูในตัวเลือกการเติมและสิ่งที่ใกล้เคียงที่สุดที่ฉันสามารถหาได้คือการเติมแบบต่อเนื่อง แต่มันใช้การไล่ระดับสีแทนการใช้เส้น ภาพนี้แสดงเอฟเฟกต์ที่ฉันพยายามสร้าง: แหล่งรูปภาพWikipedia Commons ฉันคิดว่าชื่อของรูปแบบนี้เรียกว่า "ซับน้ำ" และใช้เพราะแผ่นพิมพ์หินไม่สามารถพิมพ์พื้นที่สีทึบ มีวิธีสร้างสไตล์นี้ใน QGIS หรือไม่?

11 qgis  polygon  style 

เป็นไปได้ไหมที่จะลากจุดยอดหลายจุดของเส้นคู่ขนานบางเส้นใน QGIS พร้อมกัน? ฉันมีเลเยอร์คู่กับเส้นขนาน ในบางภูมิภาคฉันต้องย้ายจุดยอดของเส้นขนานทั้งหมดเพื่อให้ระยะห่างระหว่างกัน ด้วยเครื่องมือโหนดฉันสามารถเลือกจุดยอดของหนึ่งในบรรทัดของฉันฉันต้องการลาก แต่ไม่หลายบรรทัดในเวลาเดียวกัน: ใน Multiline-Layers เป็นไปได้ที่จะลากทุกบรรทัดในเวลาเดียวกัน แต่นั่นไม่ใช่ตัวเลือกเพราะฉันต้องรักษาคุณสมบัติทั้งหมดของบรรทัดเดียวของฉัน มีตัวเลือกอื่น ๆ เพื่อลากเส้นของฉันทั้งหมดในเวลาเดียวกันหรือไม่? แก้ไข 1: เพิ่งพบCadNodeTool ( https://github.com/wonder-sk/CadNodeTool ) ด้วยปลั๊กอินนี้มันเป็นไปได้อย่างน้อยที่สุดแล้วที่จะเลือกจุดยอดของหลายบรรทัดในเวลาเดียวกัน บางทีปลั๊กอินนี้สามารถขยายเพื่อย้ายจุดยอดของเส้นคู่ขนาน

11 qgis  pyqgis  editing  vertices 

ฉันมี และNE lat lng SW lat lngเป้าหมายของฉันง่ายมากฉันใช้ RStudio และฉันต้องการสร้างกล่องกระโดดจากสองด้านบนและด้านล่างจากนั้นสร้างรูปหลายเหลี่ยมและฉันต้องการเพิ่มจุดสุ่มลงในรูปหลายเหลี่ยมสี่เหลี่ยมจัตุรัส ฉันมี NELat / Lng และ SW-Lat / Lng ในไฟล์ csv ที่ฉันสามารถนำเข้าสู่ data.frame ฉันแค่ต้องการคำแนะนำเล็กน้อยในการเริ่มต้นฉันยังใหม่กับภาษา R มันแตกต่างจากที่ฉันคุ้นเคย จนถึงตอนนี้ฉันมีสิ่งนี้: coords = cbind(78.46801, 19.53407) coordsmax = cbind(78.83157, 19.74557 ) sp = SpatialPoints(coordsmax) sp2 = SpatialPoints(coords) r1 = rbind(coords, coordsmax[1, ]) # join P1 = Polygon(r1) …

11 r  rgdal 

ฉันดาวน์โหลด Sentinel-2 ไทล์จากAmazon S3แต่ดูเหมือนว่าพวกเขาจะไม่เรียงกัน 100% ภาพที่ถูกต้องคือแบนด์ที่ 08 Sentinel-2 ถูกเพิ่มลงใน ArcGis Desktop (ฉันไม่ได้ทำอะไรมากไปกว่าการเพิ่ม 08.jp2 ลงในแผนที่), ทิ้งภาพ esri basemap ค่าชดเชยใน Sentinel อยู่ที่ประมาณ 35-40 เมตร ฉันจะแก้ไขให้ถูกต้องได้อย่างไร นี่คือกระเบื้องที่ดาวน์โหลด ระดับ -1C, Sentinel-2: วิธีการแก้: จากคำตอบของ HDunn ฉันทำสิ่งนี้: การใช้เลเยอร์เวกเตอร์ถนนฉันเพิ่มห้าจุดคู่และคำนวณความแตกต่างในพิกัด x และ y ความแตกต่างเฉลี่ยที่ฉันใช้เป็นค่า xvalue และ y ในเครื่องมือ Shift ผลลัพธ์:

11 arcgis-desktop  remote-sensing  sentinel-2 

ฉันกำลังพยายามส่งออกไฟล์ภาพ PNG จาก ArcGIS Pro 1.2.0 โดยตั้งค่าพื้นหลังเป็นแบบโปร่งใส ใน ArcMap 10.3 ง่ายมากมีปุ่มสำหรับตั้งค่าสีโปร่งใสเป็นสีขาวในระหว่างการส่งออก เมื่อฉันไปที่ "แชร์" ภาพใน ArcGIS Pro ฉันไม่เห็นตัวเลือกเดียวกันในการตั้งค่าพื้นหลังหรือสีโปร่งใส ปุ่มนี้ถูกย้ายไปที่อื่นหรือไม่? หรือมันไม่ใช่ตัวเลือกอีกต่อไป? เป้าหมายคือเพื่อค้นหาการตั้งค่าที่คล้ายกับภาพหน้าจอด้านบน สำหรับโครงการนี้ฉันมี basemap เป็นไฟล์ PDF แบบสแตนด์อโลน ฉันจำเป็นต้องวางกล่องไฮไลต์ไว้ด้านบนของไฟล์นั้น - นี่ต้องให้กล่องไฮไลต์มีความใสอยู่รอบตัวไม่ใช่กล่องสีขาว (เหมือนกับที่เกิดขึ้นเมื่อคุณคลิป)

11 arcgis-pro  transparency  png 

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

11 time 

ด้วยรหัสต่อไปนี้ฉันสามารถตั้งค่าความโปร่งใส 100% เป็นค่า raster 0: map=None for lyr in QgsMapLayerRegistry.instance().mapLayers().values(): if lyr.name() == "some_layer": map= lyr tr=None tr=QgsRasterTransparency() tr.initializeTransparentPixelList(0) map.renderer().setRasterTransparency(tr) map.triggerRepaint() ฉันจะกำหนดสองค่า (0 และ 8) ให้โปร่งใส 100% ได้อย่างไร เมื่อฉันทำซ้ำคำสั่งของบล็อกสุดท้ายด้วยค่าแรสเตอร์ 8 ค่าแรก (0) จะถูกลบออกจากเลเยอร์โปร่งใส

11 qgis  python  pyqgis  transparency 

ฉันมีโครงการ QGIS (2.14.3-Essen) ที่ฉันใช้ทั้งที่ทำงานและที่บ้าน ฉันเก็บทั้งหมด แต่สองชั้นกับโครงการ เลเยอร์ที่เหลือทั้งสองนั้นใหญ่มาก - DEM และฮิลล์เชดรวมมากกว่า 20 GB - และพวกมันถูกเก็บไว้ในเส้นทางที่แตกต่างกันในที่ทำงานเทียบกับที่บ้าน ฉันต้องการเพิกเฉยต่อHandle Bad Layersบทสนทนาและรักษาเส้นทางที่แตกสำหรับเลเยอร์เหล่านั้น ฉันมักจะไม่สนใจพวกเขาไม่ดี; ฉันต้องการเปลี่ยนแปลงส่วนอื่น ๆ ของโครงการ มีวิธีทำเช่นนั้นเพื่อเพิกเฉยหรือข้ามบทสนทนาที่จับเลเยอร์ไม่ดีดังนั้นฉันจะไม่สูญเสียการอ้างอิงไปยังเลเยอร์หากฉันไม่อัปเดต

11 qgis  layers 

ฉันถูกถามในหลักสูตรฝึกอบรมเมื่อเร็ว ๆ นี้ว่า QGIS สามารถคำนวณหมายเลขหน้าถัดไป / ก่อนหน้าและด้านบน / ด้านล่างโดยอัตโนมัติสำหรับหนังสือแผนที่ที่สร้างโดยใช้ตัวสร้างแผนที่ ฉันจัดการกับการแสดงออกของฉลากที่สมเหตุสมผลสำหรับกริดปกติถ้าคุณทราบความกว้างและความสูงของกริด แต่จากนั้นเราก็เริ่มนึกถึงตัวอย่างที่เป็นจริงที่เราไม่ต้องการวาดหน้าเว็บที่ไม่มีเขตที่เราสนใจเช่นหนึ่งในเขตของฉัน: ดังนั้นบ่ายนี้ฉันได้เล่นสคริปต์ไพ ธ อนเพื่อหา 4 เพื่อนบ้านที่ฉันสนใจสำหรับแต่ละตารางเซลล์และเพิ่มค่าเหล่านั้นในตารางของฉัน (นี่ขึ้นอยู่กับการสอนของ Ujaval Gandhi ) for f in feature_dict.values(): print 'Working on %s' % f[_NAME_FIELD] geom = f.geometry() # Find all features that intersect the bounding box of the current feature. # We use spatial …

11 qgis  pyqgis  print-composer  atlas